Integrated computational approach to the analysis of RNA-seq data reveals new transcriptional regulators of psoriasis

Polo-like kinase 1 (PLK1) is a novel CARD14-binding protein in keratinocytes

Caspase recruitment domain (CARD)-containing protein 14 (CARD14) is an intracellular protein that mediates nuclear factor-kappa B (NF-ĸB) signaling and proinflammatory gene expression in skin keratinocytes. Several hyperactivating CARD14 mutations have been associated with psoriasis and other inflammatory skin diseases. CARD14-induced NF-ĸB signaling is dependent on the formation of a CARD14-BCL10-MALT1 (CBM) signaling complex, but upstream receptors and molecular mechanisms that activate and regulate CARD14 signaling are still largely unclear. Using unbiased affinity purification and mass spectrometry (AP-MS) screening, we discover polo-like kinase 1 (PLK1) as a novel CARD14-binding protein. CARD14-PLK1 binding is independent of the CARD14 CARD domain but involves a consensus phospho-dependent PLK1-binding motif in the CARD14 linker region (LR). Expression of the psoriasis-associated CARD14(E138A) variant in human keratinocytes induces the recruitment of PLK1 to CARD14-containing signalosomes in interphase cells, but does not affect the specific location of PLK1 in mitotic cells. Finally, disruption of the PLK1-binding motif in CARD14(E138A) increases CARD14-induced proinflammatory signaling and gene expression. Together, our data identify PLK1 as a novel CARD14-binding protein and indicate a negative regulatory role for PLK1 in CARD14 signaling.

The Aryl Hydrocarbon Receptor Regulates Epidermal Differentiation through Transient Activation of TFAP2A

The aryl hydrocarbon receptor (AHR) is an evolutionary conserved environmental sensor identified as an indispensable regulator of epithelial homeostasis and barrier organ function. Molecular signaling cascade and target genes upon AHR activation and their contribution to cell and tissue function are however not fully understood. Multiomics analyses using human skin keratinocytes revealed that upon ligand activation, AHR binds open chromatin to induce expression of transcription factors, for example, TFAP2A, as a swift response to environmental stimuli. The terminal differentiation program, including upregulation of barrier genes, FLG and keratins, was mediated by TFAP2A as a secondary response to AHR activation. The role of AHR-TFAP2A axis in controlling keratinocyte terminal differentiation for proper barrier formation was further confirmed using CRISPR/Cas9 in human epidermal equivalents. Overall, the study provides additional insights into the molecular mechanism behind AHR-mediated barrier function and identifies potential targets for the treatment of skin barrier diseases.

Cutaneous lesions in psoriatic arthritis are enriched in chemokine transcriptomic pathways

OBJECTIVES

Skin from people with psoriasis has been extensively studied and is assumed to be identical to skin from those with psoriatic arthritis (PsA). Chemokines and the CC chemokine scavenger receptor ACKR2 are upregulated in uninvolved psoriasis. ACKR2 has been proposed as a regulator of cutaneous inflammation in psoriasis. The aim of this study was to compare the transcriptome of PsA skin to healthy control (HC) skin and evaluate ACKR2 expression in PsA skin.

METHODS

Full-thickness skin biopsies from HC, lesional and uninvolved skin from participants with PsA were sequenced on NovaSeq 6000. Findings were validated using qPCR and RNAscope.

RESULTS

Nine HC and nine paired PsA skin samples were sequenced. PsA uninvolved skin was transcriptionally similar to HC skin, and lesional PsA skin was enriched in epidermal and inflammatory genes. Lesional PsA skin was enriched in chemokine-mediated signalling pathways, but uninvolved skin was not. ACKR2 was upregulated in lesional PsA skin but had unchanged expression in uninvolved compared with HC skin. The expression of ACKR2 was confirmed by qPCR, and RNAscope demonstrated strong expression of ACKR2 in the suprabasal layer of the epidermis in PsA lesions.

CONCLUSION

Chemokines and their receptors are upregulated in lesional PsA skin but relatively unchanged in uninvolved PsA skin. In contrast to previous psoriasis studies, ACKR2 was not upregulated in uninvolved PsA skin. Further understanding of the chemokine system in PsA may help to explain why inflammation spreads from the skin to the joints in some people with psoriasis.

Vitamin D and Systems Biology

The biological actions of the vitamin D receptor (VDR) have been investigated intensively for over 100 years and has led to the identification of significant insights into the repertoire of its biological actions. These were initially established to be centered on the regulation of calcium transport in the colon and deposition in bone. Beyond these well-known calcemic roles, other roles have emerged in the regulation of cell differentiation processes and have an impact on metabolism. The purpose of the current review is to consider where applying systems biology (SB) approaches may begin to generate a more precise understanding of where the VDR is, and is not, biologically impactful. Two SB approaches have been developed and begun to reveal insight into VDR biological functions. In a top-down SB approach genome-wide scale data are statistically analyzed, and from which a role for the VDR emerges in terms of being a hub in a biological network. Such approaches have confirmed significant roles, for example, in myeloid differentiation and the control of inflammation and innate immunity. In a bottom-up SB approach, current biological understanding is built into a kinetic model which is then applied to existing biological data to explain the function and identify unknown behavior. To date, this has not been applied to the VDR, but has to the related ERα and identified previously unknown mechanisms of control. One arena where applying top-down and bottom-up SB approaches may be informative is in the setting of prostate cancer health disparities.

Transcription Factors Runx1 and Runx3 Suppress Keratin Expression in Undifferentiated Keratinocytes

The Runt-related transcription factor (Runx) family has been suggested to play roles in stem cell regulation, tissue development, and oncogenesis in various tissues/organs. In this study, we investigated the possible functions of Runx1 and Runx3 in keratinocyte differentiation. Both Runx1 and Runx3 proteins were detected in primary cultures of mouse keratinocytes. Proteins were localized in the nuclei of undifferentiated keratinocytes but translocated to the cytoplasm of differentiated cells. The siRNA-mediated inhibition of Runx1 and Runx3 expression increased expression of keratin 1 and keratin 10, which are early differentiation markers of keratinocytes. In contrast, overexpression of Runx1 and Runx3 suppressed keratin 1 and keratin 10 expression. Endogenous Runx1 and Runx3 proteins were associated with the promoter sequences of keratin 1 and keratin 10 genes in undifferentiated but not differentiated keratinocytes. In mouse skin, the inhibition of Runx1 and Runx3 expression by keratinocyte-specific gene targeting increased the ratios of keratin 1- and keratin 10-positive cells in the basal layer of the epidermis. On the other hand, inhibition of Runx1 and Runx3 expression did not alter the proliferation capacity of cultured or epidermal keratinocytes. These results suggest that Runx1 and Runx3 likely function to directly inhibit differentiation-induced expression of keratin 1 and keratin 10 genes but are not involved in the regulation of keratinocyte proliferation.

Transcriptomic Profiling of Peripheral Edge of Lesions to Elucidate the Pathogenesis of Psoriasis Vulgaris

Elucidating transcriptome in the peripheral edge of the lesional (PE) skin could provide a better understanding of the molecules or signalings that intensify inflammation in the PE skin. Full-thickness biopsies of PE skin and uninvolved (UN) skin were obtained from psoriasis patients for RNA-seq. Several potential differentially expressed genes (DEGs) in the PE skin compared to those in the UN skin were identified. These DEGs enhanced functions such as angiogenesis, growth of epithelial tissue, chemotaxis and homing of cells, growth of connective tissues, and degranulation of myeloid cells beneath the PE skin. Moreover, the canonical pathways of IL-17A, IL-6, and IL-22 signaling were enriched by the DEGs. Finally, we proposed that inflammation in the PE skin might be driven by the IL-36/TLR9 axis or IL-6/Th17 axis and potentiated by IL-36α, IL-36γ, IL-17C, IL-8, S100A7, S100A8, S100A9, S100A15, SERPINB4, and hBD-2. Along with IL-36α, IL-17C, and IκBζ, ROCK2 could be an equally important factor in the pathogenesis of psoriasis, which may involve self-sustaining circuits between innate and adaptive immune responses via regulation of IL-36α and IL-36γ expression. Our finding provides new insight into signaling pathways in PE skin, which could lead to the discovery of new psoriasis targets.

Determining the immune environment of cutaneous T-cell lymphoma lesions through the assessment of lesional blood drops

Detailed analysis of the cells that infiltrate lesional skin cannot be performed in skin biopsy specimens using immunohistochemistry or cell separation techniques because enzyme treatments applied during the isolation step can destroy small amounts of protein and minor cell populations in the biopsy specimen. Here, we describe a method for isolating T cells from drops of whole blood obtained from lesions during skin biopsy in patients with cutaneous T-cell lymphoma. Lesional blood is assumed to contain lesional resident cells, cells from capillary vessels, and blood overflowing from capillary vessels into the lesion area. The lesional blood showed substantial increases in distinct cell populations, chemokines, and the expression of various genes. The proportion of CD8+CD45RO+ T cells in the lesional blood negatively correlated with the modified severity-weighted assessment tool scores. CD4+CD45RO+ T cells in the lesional blood expressed genes associated with the development of cancer and progression of cutaneous T-cell lymphoma. In addition, CD8+CD45RO+ T cells in lesional blood had unique T-cell receptor repertoires in lesions of each stage. Assessment of lesional blood drops might provide new insight into the pathogenesis of mycosis fungoides and facilitate evaluation of the treatment efficacy for mycosis fungoides as well as other skin inflammatory diseases.

H3K27Ac modification and gene expression in psoriasis

BACKGROUND

Numerous alterations in gene expression have been described in psoriatic lesions compared to uninvolved or healthy skin. However, the mechanisms which induce this altered expression remain unclear. Epigenetic modifications play a key role in regulating genes' expression. Only three studies compared the whole-genome DNA methylation of psoriasis versus healthy skin. The present is the first study of genome-wide comparison of histone modifications between psoriatic to healthy skins.

OBJECTIVE

Our objective was to explore the pattern of H3K27Ac modifications in psoriatic lesions compared to uninvolved psoriatic and healthy skin, in order to identify new genes involved in the pathogenesis of psoriasis.

METHOD

Using ChIP-seq with anti H3K27Ac we compared the acetylation of lysine 27 on histone 3 (H3K27Ac) modification between psoriatic to healthy skins, combined with mRNA array.

RESULTS

We found a differential H3K27Ac pattern between psoriatic compared to uninvolved or healthy skins. We found that many of the overexpressed and H3K27Ac enriched genes in psoriasis, harbor a putative GRHL transcription factor-binding site.

CONCLUSIONS

In the most overexpressed genes in psoriasis, there is an enrichment of H3K27Ac. However, the loss of H3K27 acetylation modification does not correlate with decreased gene expression. GRHL appears to play an important role in the pathogenesis of psoriasis and therefore, might be a new target for psoriasis therapeutics.

Predicted regulatory SNPs reveal potential drug targets and novel companion diagnostics in psoriasis

Psoriasis is an autoimmune disease associated with interleukins, their receptors, key transcription factors and more recently, antimicrobial peptides (AMPs). Cathelicidin LL-37 is an AMP proposed to play a fundamental role in psoriasis etiology. With our proprietary software SNPClinic v.1.0, we analyzed 203 common SNPs (MAF frequency ​> ​1%) in proximal promoters of 22 genes associated with psoriasis. These include nine genes which protein products are classic drug targets for psoriasis (TNF, IL17A, IL17B, IL17C, IL17F, IL17RA, IL12A, IL12B and IL23A). SNPClinic predictions were run with DNAseI-HUP chromatin accessibility data in eight psoriasis/epithelia-relevant cell lines from ENCODE including keratinocytes (NHEK), T_H1 and T_H17 lymphocytes. Results were ranked quantitatively by transcriptional relevance according to our novel Functional Impact Factor (FIF) parameter. We found six rSNPs in five genes (CAMP/cathelicidin, S100A7/psoriasin, IL17C, IL17RA and TNF) and each was confirmed as true rSNP in at least one public eQTL database including GTEx portal and ENCODE (Phase 3). Predicted regulatory SNPs in cathelicidin, IL17C and IL17RA genes may explain hyperproliferation of keratinocytes. Predicted rSNPs in psoriasin, IL17C and cathelicidin may contribute to activation and polarization of lymphocytes. Predicted rSNPs in TNF gene are concordant with the epithelium-mesenchymal transition. In spite that these results must be validated in vitro and in vivo with a functional genomics approach, we propose FOXP2, RUNX2, NR2F1, ELF1 and HESX1 transcription factors (those with the highest FIF on each gene) as novel drug targets for psoriasis. Furthermore, four out of six rSNPs uncovered by SNPClinic v.1.0 software, could also be validated in the clinic as companion diagnostics/pharmacogenetics assays for psoriasis prescribed drugs that block TNF-α (e.g. Etanercept), IL-17 (e.g. Secukinumab) and IL-17 receptor (Brodalumab).

•

We found six putative regulatory SNPs in cathelicidin (LL-37), psoriasin (S100A7), IL17C, IL17RA and TNF genes.

•

These rSNPs could be validated also as companion diagnostics/pharmacogenetics assays for most approved psoriasis drugs.

•

Regulatory SNPs in TNF gene are concordant with the epithelial-mesenchymal transition.

•

Regulatory SNPs in IL17C and IL17RA may partially explain hyperproliferation of keratinocytes.

•

Regulatory SNP rs12049559 in psoriasin (S100A7) may contribute to T-cell polarization.

Using Micro- and Macro-Level Network Metrics Unveils Top Communicative Gene Modules in Psoriasis

(1) Background: Psoriasis is a multifactorial chronic inflammatory disorder of the skin, with significant morbidity, characterized by hyperproliferation of the epidermis. Even though psoriasis’ etiology is not fully understood, it is believed to be multifactorial, with numerous key components. (2) Methods: In order to cast light on the complex molecular interactions in psoriasis vulgaris at both protein–protein interactions and transcriptomics levels, we studied a set of microarray gene expression analyses consisting of 170 paired lesional and non-lesional samples. Afterwards, a network analysis was conducted on the protein–protein interaction network of differentially expressed genes based on micro- and macro-level network metrics at a systemic level standpoint. (3) Results: We found 17 top communicative genes, all of which were experimentally proven to be pivotal in psoriasis, which were identified in two modules, namely the cell cycle and immune system. Intra- and inter-gene interaction subnetworks from the top communicative genes might provide further insight into the corresponding characteristic interactions. (4) Conclusions: Potential gene combinations for therapeutic/diagnostics purposes were identified. Moreover, our proposed workflow could be of interest to a broader range of future biological network analysis studies.

The miR-155/GATA3/IL37 axis modulates the production of proinflammatory cytokines upon TNF-α stimulation to affect psoriasis development

Psoriasis is a recrudescent chronic immune-mediated inflammatory dermatosis; the production and release of proinflammatory cytokines/chemokines such as TNF-α has been regarded as critical issues during psoriasis pathogenesis. Based on online microarray profiles, the expression of the transcription factor GATA3 was downregulated in psoriasis lesion tissues. In the present study, we searched for miRNAs that might be related to TNF-α and GATA3 to investigate an in-depth understanding of psoriasis pathogenesis. Herein, higher TNF-α and GATA3 protein levels were observed in psoriasis lesion tissues and that GATA3 overexpression significantly reverses TNF-α-induced increases within the production of IL-6 and CXCL8 in keratinocytes. TNF-α stimulation increases miR-155 expression dose-independently, and the miR-155 inhibitor significantly reverses TNF-α-induced suppression of GATA3 protein levels and increases IL-6 and CXCL8 production. miR-155 could suppress the expression of GATA3 by targeting its 3'UTR, while GATA3 could activate the transcription of IL37 by targeting its promoter region. miR-155 overexpression reduces IL37 protein and increases CXCL8 production; GATA3 overexpression might significantly attenuate the effects of miR-155 overexpression. In contrast to GATA3, miR-155 expression is significantly upregulated in psoriasis lesion tissue and is negatively correlated with GATA3 and IL37. In summary, the miR-155/GATA3/IL37 axis modulates the production of IL-6 and CXCL8 upon TNF-α stimulation to affect psoriasis development. Thus, miR-155/GATA3/IL37 may be potent targets for psoriasis treatment, which needs further in vivo and clinical investigation.

Review-Current Concepts in Inflammatory Skin Diseases Evolved by Transcriptome Analysis: In-Depth Analysis of Atopic Dermatitis and Psoriasis

During the last decades, high-throughput assessment of gene expression in patient tissues using microarray technology or RNA-Seq took center stage in clinical research. Insights into the diversity and frequency of transcripts in healthy and diseased conditions provide valuable information on the cellular status in the respective tissues. Growing with the technique, the bioinformatic analysis toolkit reveals biologically relevant pathways which assist in understanding basic pathophysiological mechanisms. Conventional classification systems of inflammatory skin diseases rely on descriptive assessments by pathologists. In contrast to this, molecular profiling may uncover previously unknown disease classifying features. Thereby, treatments and prognostics of patients may be improved. Furthermore, disease models in basic research in comparison to the human disease can be directly validated. The aim of this article is not only to provide the reader with information on the opportunities of these techniques, but to outline potential pitfalls and technical limitations as well. Major published findings are briefly discussed to provide a broad overview on the current findings in transcriptomics in inflammatory skin diseases.

Identification of Novel Genes in Osteoarthritic Fibroblast-Like Synoviocytes Using Next-Generation Sequencing and Bioinformatics Approaches

Synovitis in osteoarthritis (OA) the consequence of low grade inflammatory process caused by cartilage breakdown products that stimulated the production of pro-inflammatory mediators by fibroblast-like synoviocytes (FLS). FLS participate in joint homeostasis and low grade inflammation in the joint microenvironment triggers FLS transformation. In the current study, we aimed to identify differentially expressed genes and potential miRNA regulations in human OA FLS through deep sequencing and bioinformatics approaches. The 245 differentially expressed genes in OA FLS were identified, and pathway analysis using various bioinformatics databases indicated their enrichment in functions related to altered extracellular matrix organization, cell adhesion and cellular movement. Moreover, among the 14 dysregulated genes with potential miRNA regulations identified, src kinase associated phosphoprotein 2 (SKAP2), adaptor related protein complex 1 sigma 2 subunit (AP1S2), PHD finger protein 21A (PHF21A), lipoma preferred partner (LPP), and transcription factor AP-2 alpha (TFAP2A) showed similar expression patterns in OA FLS and OA synovial tissue datasets in Gene Expression Omnibus database. Ingenuity Pathway Analysis identified the dysregulated LPP participated in cell migration and cell spreading of OA FLS, which was potentially regulated by miR-141-3p. The current findings suggested new perspectives into understanding the novel molecular signatures of FLS involved in the pathogenesis of OA, which may be potential therapeutic targets.

FRA1 mediates the activation of keratinocytes: Implications for the development of psoriatic plaques

In this study we investigated the role of FRA1, a transcription factor from the AP-1 family, in the regulation of keratinocyte characteristics important for the development of psoriatic plaques. FRA1 is characterized by elevated expression in the skin of psoriasis patients, thus leading us to predict it to be one of the major regulators of keratinocyte phenotype during the development of psoriatic lesions. Pathway analysis of RNAseq data allowed us to identify FRA1-mediated signaling cascades leading to the manifestation of the most prominent skin characteristics of the disease: the development of inflammation, epithelial-mesenchymal transition, activation of metalloproteases, and keratinocyte proliferation and migration. We have confirmed that FRA1-overexpressing keratinocytes produce elevated amounts of proinflammatory cytokines and active matrix metalloproteases, leading to the induction of the autoinflammatory loop and paracrine activation in neighbor cells. Therefore, the elevated expression of FRA1 and its altered transcriptional regulation in the skin of patients with psoriasis is an important driving factor in the development of psoriatic plaques.

RNAI-MEDIATED SILENCING OF MATRIX METALLOPROTEINASE 1 IN EPIDERMAL KERATINOCYTES INFLUENCES THE BIOLOGICAL EFFECTS OF INTERLEUKIN 17A

Matrix metalloproteinases (MMPs) are important for the pathogenesis of psoriasis and other autoimmune disorders. In the extracellular matrix, accumulation of proinflammatory cytokines, such as interleukin 17A (IL-17A), leads to induction of several MMPs, including MMP1. MMPs change the composition and other properties of the extracellular matrix. These changes facilitate tissue remodeling and promote the development of psoriatic plaques. The aim of this study was to explore how MMP1 silencing might influence the biological effects of IL-17A on migration and proliferation of human epidermal keratinocytes and the expression of genes involved in their division and differentiation. The experiments were performed with MMP1-deficient and control epidermal keratinocytes, HaCaT-MMP1 and HaCaT-KTR, respectively. Cell proliferation and migration were assessed by comparative analysis of the growth curves and scratch assay, respectively. To quantify cell migration, representative areas of cell cultures were photographed at the indicated time points and compared to each other. Changes in gene expression were analyzed by real-time PCR. The obtained results demonstrated that MMP1 silencing in the cells treated with IL-17A resulted in downregulation of MMP9 and -12, FOSL1, CCNA2, IVL, KRT14 and -17 as well as upregulation of MMP2, CCND1 and LOR. Moreover, MMP1 silencing led to a decrease in cell proliferation and an impairment of cell migration. Thus, MMP1-deficiency in epidermal keratinocytes can be beneficial for psoriasis patients that experience an accumulation of IL-17 in lesional skin. Knocking MMP1 down could influence migration and proliferation of epidermal keratinocytes in vivo, as well as help to control the expression of MMP1, -2, -9 и -12, CCNA2, CCND1, KRT14 and -17 that are crucial for the pathogenesis of psoriasis.

A catalog of potential putative functional variants in psoriasis genome-wide association regions

Psoriasis is a common inflammatory skin disease, with considerable genetic contribution. Genome-wide association studies have successfully identified a number of genomic regions for the risk of psoriasis. However, it is challenging to pinpoint the functional causal variants and then further decipher the genetic mechanisms underlying each region. In order to prioritize potential functional causal variants within psoriasis susceptibility regions, we integrated the genetic association findings and functional genomic data publicly available, i.e. histone modifications in relevant immune cells. We characterized a pervasive enrichment pattern of psoriasis variants in five core histone marks across immune cells/tissues. We discovered that genetic alleles within psoriasis association regions might influence gene expression levels through significantly affecting the binding affinities of 17 transcription factors. We established a catalog of 654 potential functional causal variants for psoriasis and suggested that they significantly overlapped with causal variants for autoimmune diseases. We identified potential causal variant rs79824801 overlay with the peaks of five histone marks in primary CD4+ T cells. Its alternative allele affected the binding affinity of transcription factor IKZF1. This study highlights the complex genetic architecture and complicated mechanisms for psoriasis. The findings will inform the functional experiment design for psoriasis.

Transcriptomic profiling of human skin biopsies in the clinical trial setting: A protocol for high quality RNA extraction from skin tumours

Transcriptomic profiling of skin disease using next generation sequencing allows for detailed information on aspects of RNA biology including gene expression, non-coding regulatory elements and gene splicing. The application of RNA sequencing to human skin disease and cancer is often hampered by degraded RNA. Here we describe a protocol that allows for consistently intact RNA to be extracted from snap frozen skin biopsy samples, which has been validated in a clinical trial setting.

Human skin tumour punch biopsies (n=28) ranging from 4-6mm in diameter were obtained from 14 patients with an inherited skin tumour syndrome (CYLD cutaneous syndrome) and frozen in liquid nitrogen prior to being stored at -80°C. These samples were then subject to cyrostat sectioning, allowing for histological assessment, and were homogenised using a bead-based lysis platform. RNA extraction was performed using a silica column-based system. RNA concentration was measured using fluorescent quantitation and RNA integrity assessed using microfluidic gel electrophoresis. We also processed normal skin biopsies using the same protocol (n=10). The mean RNA integrity score of the tumour and normal samples was 9.5, and the quantity of RNA obtained from the small amounts of tissue used exceeded requirements for RNA-seq library generation.

We propose that the method of RNA extraction suggested here allows for transcriptomic profiling from small pieces of human tissue without the need for PCR amplification during library preparation. This protocol could be utilised in healthy and diseased skin to improve mechanistic understanding in a range of human skin diseases.

Bioinformatic approaches to interrogating vitamin D receptor signaling

Bioinformatics applies unbiased approaches to develop statistically-robust insight into health and disease. At the global, or "20,000 foot" view bioinformatic analyses of vitamin D receptor (NR1I1/VDR) signaling can measure where the VDR gene or protein exerts a genome-wide significant impact on biology; VDR is significantly implicated in bone biology and immune systems, but not in cancer. With a more VDR-centric, or "2000 foot" view, bioinformatic approaches can interrogate events downstream of VDR activity. Integrative approaches can combine VDR ChIP-Seq in cell systems where significant volumes of publically available data are available. For example, VDR ChIP-Seq studies can be combined with genome-wide association studies to reveal significant associations to immune phenotypes. Similarly, VDR ChIP-Seq can be combined with data from Cancer Genome Atlas (TCGA) to infer the impact of VDR target genes in cancer progression. Therefore, bioinformatic approaches can reveal what aspects of VDR downstream networks are significantly related to disease or phenotype.

Identification of Transcriptional Regulators of Psoriasis from RNA-Seq Experiments

Psoriasis is a common inflammatory skin disease with complex etiology and chronic progression. To provide novel insights into the molecular mechanisms of regulation of the disease we performed RNA sequencing (RNA-Seq) analysis of 14 pairs of skin samples collected from psoriatic patients. Subsequent pathway analysis and an extraction of transcriptional regulators governing psoriasis-associated pathways was executed using a combination of MetaCore Interactome enrichment tool and cisExpress algorithm, and followed by comparison to a set of previously described psoriasis response elements. A comparative approach has allowed us to identify 42 core transcriptional regulators of the disease associated with inflammation (NFkB, IRF9, JUN, FOS, SRF), activity of T-cells in the psoriatic lesions (STAT6, FOXP3, NFATC2, GATA3, TCF7, RUNX1, etc.), hyperproliferation and migration of keratinocytes (JUN, FOS, NFIB, TFAP2A, TFAP2C), and lipid metabolism (TFAP2, RARA, VDR). After merging the ChIP-seq and RNA-seq data, we conclude that the atypical expression of FOXA1 transcriptional factor is an important player in psoriasis, as it inhibits maturation of naive T cells into this Treg subpopulation (CD4+FOXA1+CD47+CD69+PD-L1(hi)FOXP3-), therefore contributing to the development of psoriatic skin lesions.

Analysis of cis-Regulatory Elements in Gene Co-expression Networks in Cancer

Analysis of gene co-expression networks is a powerful "data-driven" tool, invaluable for understanding cancer biology and mechanisms of tumor development. Yet, despite of completion of thousands of studies on cancer gene expression, there were few attempts to normalize and integrate co-expression data from scattered sources in a concise "meta-analysis" framework. Here we describe an integrated approach to cancer expression meta-analysis, which combines generation of "data-driven" co-expression networks with detailed statistical detection of promoter sequence motifs within the co-expression clusters. First, we applied Weighted Gene Co-Expression Network Analysis (WGCNA) workflow and Pearson's correlation to generate a comprehensive set of over 3000 co-expression clusters in 82 normalized microarray datasets from nine cancers of different origin. Next, we designed a genome-wide statistical approach to the detection of specific DNA sequence motifs based on similarities between the promoters of similarly expressed genes. The approach, realized as cisExpress software module, was specifically designed for analysis of very large data sets such as those generated by publicly accessible whole genome and transcriptome projects. cisExpress uses a task farming algorithm to exploit all available computational cores within a shared memory node.We discovered that although co-expression modules are populated with different sets of genes, they share distinct stable patterns of co-regulation based on promoter sequence analysis. The number of motifs per co-expression cluster varies widely in accordance with cancer tissue of origin, with the largest number in colon (68 motifs) and the lowest in ovary (18 motifs). The top scored motifs are typically shared between several tissues; they define sets of target genes responsible for certain functionality of cancerogenesis. Both the co-expression modules and a database of precalculated motifs are publically available and accessible for further studies.