Mrinal K. Sarkar

Analysis of long non-coding RNAs highlights tissue-specific expression patterns and epigenetic profiles in normal and psoriatic skin

BackgroundAlthough analysis pipelines have been developed to use RNA-seq to identify long non-coding RNAs (lncRNAs), inference of their biological and pathological relevance remains a challenge. As a result, most transcriptome studies of autoimmune disease have only assessed protein-coding transcripts.ResultsWe used RNA-seq data from 99 lesional psoriatic, 27 uninvolved psoriatic, and 90 normal skin biopsies, and applied computational approaches to identify and characterize expressed lncRNAs. We detect 2,942 previously annotated and 1,080 novel lncRNAs which are expected to be skin specific. Notably, over 40% of the novel lncRNAs are differentially expressed and the proportions of differentially expressed transcripts among protein-coding mRNAs and previously-annotated lncRNAs are lower in psoriasis lesions versus uninvolved or normal skin. We find that many lncRNAs, in particular those that are differentially expressed, are co-expressed with genes involved in immune related functions, and that novel lncRNAs are enriched for localization in the epidermal differentiation complex. We also identify distinct tissue-specific expression patterns and epigenetic profiles for novel lncRNAs, some of which are shown to be regulated by cytokine treatment in cultured human keratinocytes.ConclusionsTogether, our results implicate many lncRNAs in the immunopathogenesis of psoriasis, and our results provide a resource for lncRNA studies in other autoimmune diseases.

Cellular dissection of psoriasis for transcriptome analyses and the post-GWAS era

BackgroundGenome-scale studies of psoriasis have been used to identify genes of potential relevance to disease mechanisms. For many identified genes, however, the cell type mediating disease activity is uncertain, which has limited our ability to design gene functional studies based on genomic findings.MethodsWe identified differentially expressed genes (DEGs) with altered expression in psoriasis lesions (n = 216 patients), as well as candidate genes near susceptibility loci from psoriasis GWAS studies. These gene sets were characterized based upon their expression across 10 cell types present in psoriasis lesions. Susceptibility-associated variation at intergenic (non-coding) loci was evaluated to identify sites of allele-specific transcription factor binding.ResultsHalf of DEGs showed highest expression in skin cells, although the dominant cell type differed between psoriasis-increased DEGs (keratinocytes, 35%) and psoriasis-decreased DEGs (fibroblasts, 33%). In contrast, psoriasis GWAS candidates tended to have highest expression in immune cells (71%), with a significant fraction showing maximal expression in neutrophils (24%, P < 0.001). By identifying candidate cell types for genes near susceptibility loci, we could identify and prioritize SNPs at which susceptibility variants are predicted to influence transcription factor binding. This led to the identification of potentially causal (non-coding) SNPs for which susceptibility variants influence binding of AP-1, NF-κB, IRF1, STAT3 and STAT4.ConclusionsThese findings underscore the role of innate immunity in psoriasis and highlight neutrophils as a cell type linked with pathogenetic mechanisms. Assignment of candidate cell types to genes emerging from GWAS studies provides a first step towards functional analysis, and we have proposed an approach for generating hypotheses to explain GWAS hits at intergenic loci.

Proteogenomic analysis of psoriasis reveals discordant and concordant changes in mRNA and protein abundance

BackgroundPsoriasis is a chronic disease characterized by the development of scaly red skin lesions and possible co-morbid conditions. The psoriasis lesional skin transcriptome has been extensively investigated, but mRNA levels do not necessarily reflect protein abundance. The purpose of this study was therefore to compare differential expression patterns of mRNA and protein in psoriasis lesions.MethodsLesional (PP) and uninvolved (PN) skin samples from 14 patients were analyzed using high-throughput complementary DNA sequencing (RNA-seq) and liquid chromatography-tandem mass spectrometry (LC-MS/MS).ResultsWe identified 4122 differentially expressed genes (DEGs) along with 748 differentially expressed proteins (DEPs). Global shifts in mRNA were modestly correlated with changes in protein abundance (r = 0.40). We identified similar numbers of increased and decreased DEGs, but 4-fold more increased than decreased DEPs. Ribosomal subunit and translation proteins were elevated within lesions, without a corresponding shift in mRNA expression (RPL3, RPS8, RPL11). We identified 209 differentially expressed genes/proteins (DEGPs) with corresponding trends at the transcriptome and proteome levels. Most DEGPs were similarly altered in at least one other skin disease. Psoriasis-specific and non-specific DEGPs had distinct cytokine-response patterns, with only the former showing disproportionate induction by IL-17A in cultured keratinocytes.ConclusionsOur findings reveal global imbalance between the number of increased and decreased proteins in psoriasis lesions, consistent with heightened translation. This effect could not have been discerned from mRNA profiling data alone. High-confidence DEGPs were identified through transcriptome-proteome integration. By distinguishing between psoriasis-specific and non-specific DEGPs, our analysis uncovered new functional insights that would otherwise have been overlooked.

Psoriasis drug development and GWAS interpretation through in silico analysis of transcription factor binding sites

BackgroundPsoriasis is a cytokine-mediated skin disease that can be treated effectively with immunosuppressive biologic agents. These medications, however, are not equally effective in all patients and are poorly suited for treating mild psoriasis. To develop more targeted therapies, interfering with transcription factor (TF) activity is a promising strategy.MethodsMeta-analysis was used to identify differentially expressed genes (DEGs) in the lesional skin from psoriasis patients (n = 237). We compiled a dictionary of 2935 binding sites representing empirically-determined binding affinities of TFs and unconventional DNA-binding proteins (uDBPs). This dictionary was screened to identify “psoriasis response elements” (PREs) overrepresented in sequences upstream of psoriasis DEGs.ResultsPREs are recognized by IRF1, ISGF3, NF-kappaB and multiple TFs with helix-turn-helix (homeo) or other all-alpha-helical (high-mobility group) DNA-binding domains. We identified a limited set of DEGs that encode proteins interacting with PRE motifs, including TFs (GATA3, EHF, FOXM1, SOX5) and uDBPs (AVEN, RBM8A, GPAM, WISP2). PREs were prominent within enhancer regions near cytokine-encoding DEGs (IL17A, IL19 and IL1B), suggesting that PREs might be incorporated into complex decoy oligonucleotides (cdODNs). To illustrate this idea, we designed a cdODN to concomitantly target psoriasis-activated TFs (i.e., FOXM1, ISGF3, IRF1 and NF-kappaB). Finally, we screened psoriasis-associated SNPs to identify risk alleles that disrupt or engender PRE motifs. This identified possible sites of allele-specific TF/uDBP binding and showed that PREs are disproportionately disrupted by psoriasis risk alleles.ConclusionsWe identified new TF/uDBP candidates and developed an approach that (i) connects transcriptome informatics to cdODN drug development and (ii) enhances our ability to interpret GWAS findings. Disruption of PRE motifs by psoriasis risk alleles may contribute to disease susceptibility.

Differentiation-dependent requirement of Tsix long non-coding RNA in imprinted X-chromosome inactivation

Imprinted X-inactivation is a paradigm of mammalian transgenerational epigenetic regulation resulting in silencing of genes on the paternally-inherited X-chromosome. The pre-programmed fate of the X-chromosomes is thought to be controlled in cis by the parent-of-origin-specific expression of two long non-coding RNAs, Tsix and Xist, in mice. Exclusive expression of Tsix from the maternal–X has implicated it as the instrument through which the maternal germline prevents inactivation of the maternal–X in the offspring. Here, we show that Tsix is dispensable for inhibiting Xist and X-inactivation in the early embryo and in cultured stem cells of extra-embryonic lineages. Tsix is instead required to prevent Xist expression as trophectodermal progenitor cells differentiate. Despite induction of wild-type Xist RNA and accumulation of histone H3-K27me3, many Tsix-mutant X-chromosomes fail to undergo ectopic X-inactivation. We propose a novel model of lncRNA function in imprinted X-inactivation that may also apply to other genomically imprinted loci.

A gene network regulated by the transcription factor VGLL3 as a promoter of sex-biased autoimmune diseases

Autoimmune diseases affect 7.5% of the US population, and they are among the leading causes of death and disability. A notable feature of many autoimmune diseases is their greater prevalence in females than in males, but the underlying mechanisms of this have remained unclear. Through the use of high-resolution global transcriptome analyses, we demonstrated a female-biased molecular signature associated with susceptibility to autoimmune disease and linked this to extensive sex-dependent co-expression networks. This signature was independent of biological age and sex-hormone regulation and was regulated by the transcription factor VGLL3, which also had a strong female-biased expression. On a genome-wide level, VGLL3-regulated genes had a strong association with multiple autoimmune diseases, including lupus, scleroderma and Sjögrens syndrome, and had a prominent transcriptomic overlap with inflammatory processes in cutaneous lupus. These results identified a VGLL3-regulated network as a previously unknown inflammatory pathway that promotes female-biased autoimmunity. They demonstrate the importance of studying immunological processes in females and males separately and suggest new avenues for therapeutic development.Autoimmune diseases affect 7.5% of the U.S. population, and are among the leading causes of death and disability. A striking feature of many autoimmune diseases is their increased prevalence in females, but the underlying mechanisms have remained unclear. Using high-resolution global transcriptome analyses we demonstrate a female-biased molecular signature associated with autoimmune disease susceptibility, and linked to extensive sex-dependent, co-expression networks. This signature was independent of biological age and sex-hormone regulation, and regulated by the transcription factor VGLL3, which also had a strong female biased expression. On a genome-wide level, VGLL3-regulated genes had a strong association with multiple autoimmune diseases including lupus, scleroderma and Sjögren’s syndrome and had a prominent transcriptomic overlap with inflammatory processes in cutaneous lupus. These results identify VGLL3-regulated gene network as a novel inflammatory pathway promoting female-biased autoimmunity, they demonstrate the importance of studying immunological processes in females and males separately, and open up new avenues for therapeutic development.

Psoriasis: a mixed autoimmune and autoinflammatory disease

In recent years marked progress has been made in our understanding of the critical biologic and immunologic pathways involved in psoriasis. Genetic studies have demonstrated that susceptibility to psoriasis involves components of both the adaptive and innate immune system and not surprisingly activation of both of these arms of the immune system is found in psoriatic skin. While adaptive immune responses predominate in chronic plaque psoriasis, innate and autoinflammatory responses dominate in pustular forms of psoriasis, with other clinical subtypes extending on a spectrum between plaque and pustular psoriasis. This makes psoriasis a unique disease where both autoimmune and autoinflammatory responses co-exist, with the balance between the two being critical in shaping its clinical presentation.

Lupus Skin Is Primed for IL-6 Inflammatory Responses through a Keratinocyte-Mediated Autocrine Type I Interferon Loop

Cutaneous lupus erythematosus is a disfiguring and common manifestation in systemic lupus erythematosus, and the etiology of this predisposition for cutaneous inflammation is unknown. Here, we sought to examine the keratinocyte as an important source of IL-6 and define the mechanism for its increased production in cutaneous lupus erythematosus. Evaluation of discoid and subacute cutaneous lupus erythematosus lesions showed significant epidermal up-regulation of IL-6 compared with control via real-time PCR and immunohistochemistry. Keratinocytes from unaffected skin of lupus patients produced significantly more IL-6 compared with healthy control subjects after exposure to toll-like receptor 2, 3, or 4 agonists or exposure to UVB radiation. Pretreatment with type I interferons (IFN-α and IFN-κ) increased IL-6 production by control keratinocytes, and type I IFN blockade decreased IL-6 secretion by lupus keratinocytes. Secretion of keratinocyte-specific IFN-κ was significantly increased after toll-like receptor 2 and UVB treatment in lupus keratinocytes, and neutralization of IFN-κ decreased IL-6 production by lupus keratinocytes. Thus, lupus keratinocytes are primed for IL-6 hyperproduction in a type I IFN-dependent manner. Increased production of IFN-κ by lupus keratinocytes drives this response, indicating that IFN-κ may play a pathogenic role in cutaneous lupus erythematosus and serve as a target for treatment.

Imiquimod has strain-dependent effects in mice and does not uniquely model human psoriasis.

BackgroundImiquimod (IMQ) produces a cutaneous phenotype in mice frequently studied as an acute model of human psoriasis. Whether this phenotype depends on strain or sex has never been systematically investigated on a large scale. Such effects, however, could lead to conflicts among studies, while further impacting study outcomes and efforts to translate research findings.MethodsRNA-seq was used to evaluate the psoriasiform phenotype elicited by 6 days of Aldara (5% IMQ) treatment in both sexes of seven mouse strains (C57BL/6 J (B6), BALB/cJ, CD1, DBA/1 J, FVB/NJ, 129X1/SvJ, and MOLF/EiJ).ResultsIn most strains, IMQ altered gene expression in a manner consistent with human psoriasis, partly due to innate immune activation and decreased homeostatic gene expression. The response of MOLF males was aberrant, however, with decreased expression of differentiation-associated genes (elevated in other strains). Key aspects of the IMQ response differed between the two most commonly studied strains (BALB/c and B6). Compared with BALB/c, the B6 phenotype showed increased expression of genes associated with DNA replication, IL-17A stimulation, and activated CD8+ T cells, but decreased expression of genes associated with interferon signaling and CD4+ T cells. Although IMQ-induced expression shifts mirrored psoriasis, responses in BALB/c, 129/SvJ, DBA, and MOLF mice were more consistent with other human skin conditions (e.g., wounds or infections). IMQ responses in B6 mice were most consistent with human psoriasis and best replicated expression patterns specific to psoriasis lesions.ConclusionsThese findings demonstrate strain-dependent aspects of IMQ dermatitis in mice. We have shown that IMQ does not uniquely model psoriasis but in fact triggers a core set of pathways active in diverse skin diseases. Nonetheless, our findings suggest that B6 mice provide a better background than other strains for modeling psoriasis disease mechanisms.

An Xist-activating antisense RNA required for X-chromosome inactivation

The transcriptional imbalance due to the difference in the number of X chromosomes between male and female mammals is remedied through X-chromosome inactivation, the epigenetic transcriptional silencing of one of the two X chromosomes in females. The X-linked Xist long non-coding RNA functions as an X inactivation master regulator; Xist is selectively upregulated from the prospective inactive X chromosome and is required in cis for X inactivation. Here we discover an Xist antisense long non-coding RNA, XistAR (Xist Activating RNA), which is encoded within exon 1 of the mouse Xist gene and is transcribed only from the inactive X chromosome. Selective truncation of XistAR, while sparing the overlapping Xist RNA, leads to a deficiency in Xist RNA expression in cis during the initiation of X inactivation. Thus, the Xist gene carries within its coding sequence an antisense RNA that drives Xist expression.