- Open Access
Network analysis of potential risk genes for psoriasis
Hereditas volume 158, Article number: 21 (2021)
Psoriasis is a complex chronic inflammatory skin disease. The aim of this study was to analyze potential risk genes and molecular mechanisms associated with psoriasis.
GSE54456, GSE114286, and GSE121212 were collected from gene expression omnibus (GEO) database. Differentially expressed genes (DEGs) between psoriasis and controls were screened respectively in three datasets and common DEGs were obtained. The biological role of common DEGs were identified by enrichment analysis. Hub genes were identified using protein–protein interaction (PPI) networks and their risk for psoriasis was evaluated through logistic regression analysis. Moreover, differentially methylated positions (DMPs) between psoriasis and controls were obtained in the GSE115797 dataset. Methylation markers were identified after comparison with the common genes.
A total of 118 common DEGs were identified, which were mainly involved in keratinocyte differentiation and IL-17 signaling pathway. Through PPI network, we identified top 10 degrees as hub genes. Among them, high expression of CXCL9 and SPRR1B may be risk factors for psoriasis. In addition, we selected 10 methylation-modified genes with the higher area under receiver operating characteristic curve (AUC) value as methylation markers. Nomogram showed that TGM6 and S100A9 may be associated with an increased risk of psoriasis.
This suggests that immune and inflammatory responses are active in keratinocytes of psoriatic skin. CXCL9, SPRR1B, TGM6 and S100A9 may be potential targets for the diagnosis and treatment of psoriasis.
Psoriasis is an immune-mediated inflammatory skin disease with important physiological and psychosocial consequences . Psoriasis has complex characteristics, affecting about 2% of the general population, and the prevalence of psoriasis varies from country to country [2, 3]. There are several clinical variants of this disease, namely plaque psoriasis, erythrodermic psoriasis and pustular psoriasis, among which plaque psoriasis accounts for more than 90% . Clinically, red patches of silver-white multilayered scales are characteristic of psoriasis, and the patient's epidermis is thickened, clearly demarcating from adjacent non-banded skin . The disease affects more than 10% of the body surface area and usually requires systemic drug therapy.
Inflammatory infiltration consisting of dermal dendritic cells, macrophages, T cells, and neutrophils is present in psoriatic plaques . Help T cells have long been recognized as an important pathogenic factor in psoriasis. Further evidence suggests that there is a strong T cell component to maintain psoriasis through Th1, Th17, and Th22 cells and their derived cytokines . The complex interaction of pro-inflammatory cytokines, chemokines, growth factors and chemical mediators initiated by Th17 cells may be the key factors to induce keratinocyte proliferation, angiogenesis and neutrophil influx, ultimately leading to keratinocyte over-proliferation and the characteristics of psoriatic plaques . Recently, mast cells have been shown to be the major producers of interleukin-22 in psoriasis and atomic dermatitis . Among various molecules associated with psoriasis, tumor necrosis factor (TNF)-α, IL-23, IL-17 and IL-22 are important regulators of psoriasis . Emerging evidence suggests that different phenotypes have different immunogenetic characteristics, which may affect treatment options .
Epidermal gene modification is considered as an essential factor in the pathogenesis of psoriasis . Epigenetic changes also play an important role in the differentiation of CD4( +) T lymphocyte subsets and in the pathogenesis of psoriasis . Among all epigenetic mechanisms, DNA methylation is one of the important factors in keratinocyte differentiation [14, 15]. Common drugs used to treat psoriasis, such as methotrexate, have been reported to interfere with the methyl transfer function of folic acid, thereby restoring normal methylation status . However, researchers still need to make greater efforts to elucidate how abnormal epigenetic modifications affect the pathogenesis of psoriasis. In addition, some environmental factors, including modifiable variables such as physical trauma, drug reversal, psychological stress, obesity, are associated with the development and deterioration of psoriasis .
In this study, psoriasis-related sequencing data in public databases were used to explore the molecular mechanisms and potential risk factors of the disease.
Materials and methods
We collected psoriasis-related datasets from the gene expression omnibus (GEO) database. GSE54456 included gene expression profile from skin samples of 92 psoriatic and 82 normal punch biopsies. The genes annotated in RefSeq was used to quantify gene expression levels. The gene expression was normalized to the number of reads per kilobase per million mapped reads (RPKM). GSE114286 included gene expression profile from 9 normal skins from healthy volunteers and 18 lesional skins from patients with psoriasis. The gene expression was normalized through RPKM. GSE121212 included gene expression profile of skin tissues obtained from a carefully matched and tightly defined cohort of 28 psoriasis patients, and 38 healthy controls. Paired reads were mapped to the human reference genome (b37), number of reads for each gene was counted using HTSeq.
Identification of differentially expressed genes
The differentially expressed genes (DEGs) between psoriasis and normal were obtained through limma R software package. The log2 fold change (FC) ≥ 2 and P-values < 0.05 were considered as a statistically significant difference.
The enrichment analysis was performed using clusterProfiler R software package for the DEGs, including gene ontology (GO) functional analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The cellular component (CC), biological process (BP) and molecular function (MF) terms belonged to GO function. P-value of < 0.05 was considered statistically significant.
Construction of protein–protein interaction (PPI) network
The DEGs were put into the online tool STRING (https://string-db.org) to construct the PPI network. A combined score of ≥ 0.5 was considered significant. The hub genes were chosen based on a higher number of associations with other genes (degree) in the PPI network.
DNA methylation analysis
GSE115797 included DNA methylation profile of 24 psoriatic disease skin tissue and adjacent normal skin samples. The signal intensities along with the detection P-values were calculated for each CpG probes after BMIQ normalization. Normal skinized average beta values were calculated using ChAMP R software. Differentially methylated positions (DMPs) between psoriasis and normal were achieved by limma R software package. Set screening threshold with false discovery rate (FDR) < 0.05.
Statistical analyses were carried out with SPSS Statistics 21.0. The differences between two groups were compared by Student t test. LRM function in RMS R software package was used for logistic regression analysis.
Differentially expressed genes in psoriasis
By comparing the differences between psoriasis and healthy controls, we obtained differentially expressed genes (DEGs). A total of 401, 1857, and 466 DEGs were obtained in GSE54456, GSE114286 and GSE121212, respectively (Fig. 1A). Among these DEGs, we found 118 genes simultaneously present in three datasets and then defined them as common DEGs (Fig. 1B).
Biological functions of differentially expressed genes
In the enrichment analysis results, 195 biological processes (BP), 14 cellular components (CC), and 38 molecular functions (MF) were enriched in 118 common DEGs (Fig. 2A, B, C). It mainly included keratinocyte differentiation, epidermal cell differentiation, regulation of leukocyte chemotaxis and inflammatory response. In addition, 22 KEGG signaling pathways were also enriched (Fig. 2D). Such as “IL-17 signaling pathway”, “NOD-like receptor signaling pathway”, “Toll-like receptor signaling pathway”, and “TNF signaling pathway”.
PPI network of differentially expressed genes
After screening using a combined score, a total of 86 DEGs entered to PPI network (Fig. 3A). Among them, we identified the top 10 genes with the highest degree as hub genes (Fig. 3B). Including IL1B, CXCL10, S100A7, IL17A, CCL20, SPRR1B, CXCL1, PI3, LCN2 and CXCL9. Surprisingly, they were up-regulated in all three datasets (Fig. 3C). The AUC values of hub genes were all greater than 0.7, which may have the ability to judge psoriasis, especially SPRR1B (Fig. 3D). Importantly, we performed logistic regression analysis and presented the risk prediction of hub genes for psoriasis through a nomogram (Fig. 3E). Results suggested that the higher the expression of CXCL9 and SPRR1B, the greater the risk of psoriasis. Hub genes were also mainly enriched in the IL-17 signaling pathway.
Identification of methylation markers
By analyzing the methylation modification differences between psoriasis and controls, we obtained 153,537 differential methylation probes (CpGs) for 3111 differentially methylated positions (DMPs) (Fig. 4A). These CpGs were mainly concentrated in body (Fig. 4B). Among CpGs with significant differences, HyperProbe was much more than HypoProbe (Fig. 4C). Among these DMPs, we found that 17 genes overlapped with common genes, which may be the genes regulated by methylation (Fig. 4D). Among them, CXCL1 was also the hub gene we identified. We screened the top 10 methylation markers with the higher AUC values, including cg00187686 (TCN1), cg06051311 (TRIM15), cg06261066 (TGM6), cg17217296 (AIM2), cg19913971 (TNIP3), cg22946974 (SPRR2F), cg02772121 (TRIM15), cg19391247 (TGM6), cg16139316 (S100A9), and cg17515347 (AIM2) (Table 1). The logistic regression analysis and nomogram suggested that the lower the methylation level of cg06261066 (TGM6), the higher risk for psoriasis, and the opposite was true for cg16139316 (S100A9) (Fig. 4E).
High-throughput sequencing studies have revealed the impact of large numbers of gene expression and epigenetic modification changes on psoriasis through transcriptome and epigenetic modification studies. However, the complementarity of multiple omics data has not been fully utilized for comprehensive systematic analysis to obtain more views on disease regulation. In this study, three datasets of psoriasis were analyzed to obtain more accurate psoriasis-related genes. The molecular dysregulation mechanism and potential risk factors of psoriasis were further explored by enrichment and network analysis. Combined with methylation analysis, we obtained the methylation modified markers, which expanded the regulatory mechanism of psoriasis. In addition to previous molecular mechanism studies, our study provided a more in-depth understanding of psoriasis epidermal gene regulation.
Among the three groups of DEGs, we found 118 common genes, which may be closely related to the dysregulation mechanism of psoriasis. These DEGs were significantly involved in psoriasis-related biological functions and signaling pathways. As one of the most famous immune processes underlying the pathogenesis of psoriasis, interleukin-17 (IL-17) pathway, showed a strong enrichment effected on psoriasis-related genes and epigenetic variation . In addition, there is an imbalance between the differentiation and proliferation of keratinocytes in patients with psoriasis, and IL-17A can promote the proliferation of epidermal keratinocytes [19, 20]. This promotes the development of thickened skin lesions infiltrated with a variety of inflammatory cells . The inhibitory effect of anti-IL-17A on psoriasis plays an important role in the early clinical, histopathological and molecular treatment of psoriasis . In general, the main function of TLRs is to mediate the inflammatory response, which has also been proved to be involved in the development of psoriasis . Keratinocytes also recognize pathogens and endogenous cellular stress signals through NOD-like receptor (NLR), thereby mediating the immune response .
Compared with healthy controls, the hub genes we identified were up-regulated in three datasets of psoriasis. Interleukin-1β interferes with epidermal homeostasis by inducing insulin resistance, thereby participating in the development of psoriasis . IL-17A activates the expression of CXCL1, CCL20 and S100A7 in keratinocytes, and activating the innate immune system . The chemokines CXCL 9 and CXCL 10 released by epidermal keratinocytes have a strong chemotactic effect on the key cell monocytes, neutrophils and T cells in psoriasis [26, 27]. Consistent with our analysis, SPRR1B has also been identified as a potential biomarker of psoriasis by other studies [28, 29]. The score of psoriasis area and severity index was positively correlated with the expression of PI3 . Lipocalin-2 (LCN2) is significantly higher in patients with psoriasis than in healthy controls, and may be used as a clinical indicator of psoriatic pruritus .
In addition, psoriasis patients have different degrees of DNA methylation compared with healthy human skin and are associated with immune activation or immunoinflammatory diseases . The presence of differentially methylated positions in patients with psoriasis may be of great interest. Transglutaminase 6 (TGM6) has been reported to be associated with keratinocyte differentiation . S100A9 can be an effective target for the treatment of psoriasis . The level of S100A9 is related to the severity of psoriasis, which induces the production of inflammatory mediators by activating TLR 2 and 4 .
The limitation of this study is that we only analyzed the expression of mRNA level in limited epidermal specimens, but lacked the detection of protein expression in specimens. Secondly, we only explored the methylation modification of genes, and the specific role in the pathogenesis of psoriasis remains to be further studied. In addition, more experiments are needed to determine whether the key genes have a role in the diagnosis and treatment of psoriasis.
This study identified multiple differentially expressed genes associated with psoriasis. Importantly, CXCL9 and SPRR1B, as well as methylation markers TGM6 and S100A9, have an impact on the risk of psoriasis. These genes were mainly enriched in keratinocyte differentiation and IL-17 signaling pathway. This suggests that these genes have great potential for the diagnosis and treatment of psoriasis.
Availability of data and materials
Data is available at NCBI GEO: GSE54456, GSE114286 and GSE121212.
Area under receiver operating characteristic curve
Differentially expressed genes
Differentially methylated positions
False discovery rate
Gene expression omnibus
Kyoto Encyclopedia of Genes and Genomes
Reads per kilobase per million mapped reads
Tumor necrosis factor
Strober BE, van der Walt JM, Armstrong AW, Bourcier M, Carvalho AVE, Chouela E, et al. Clinical goals and barriers to effective psoriasis care. Dermatol Ther. 2019;9(1):5–18.
Boehncke WH, Schon MP. Psoriasis. Lancet. 2015;386(9997):983–94.
Ding X, Wang T, Shen Y, Wang X, Zhou C, Tian S, et al. Prevalence of psoriasis in China: a population-based study in six cities. Eur J Dermatol. 2012;22(5):663–7.
Tampa M, Sarbu MI, Mitran MI, Mitran CI, Matei C, Georgescu SR. The pathophysiological mechanisms and the quest for biomarkers in psoriasis, a stress-related skin disease. Dis Markers. 2018;2018:5823684.
Huang TH, Lin CF, Alalaiwe A, Yang SC, Fang JY. Apoptotic or antiproliferative activity of natural products against keratinocytes for the treatment of psoriasis. Int J Mol Sci. 2019;20(10):2558.
Rendon A, Schakel K. Psoriasis pathogenesis and treatment. Int J Mol Sci. 2019;20(6):1475.
Golden JB, McCormick TS, Ward NL. IL-17 in psoriasis: implications for therapy and cardiovascular co-morbidities. Cytokine. 2013;62(2):195–201.
Ayala-Fontanez N, Soler DC, McCormick TS. Current knowledge on psoriasis and autoimmune diseases. Psoriasis (Auckl). 2016;6:7–32.
Mashiko S, Bouguermouh S, Rubio M, Baba N, Bissonnette R, Sarfati M. Human mast cells are major IL-22 producers in patients with psoriasis and atopic dermatitis. J Allergy Clin Immunol. 2015;136(2):351-9.e1.
Woo YR, Cho DH, Park HJ. Molecular mechanisms and management of a cutaneous inflammatory disorder: psoriasis. Int J Mol Sci. 2017;18(12):2684.
Mahil SK, Capon F, Barker JN. Update on psoriasis immunopathogenesis and targeted immunotherapy. Semin Immunopathol. 2016;38(1):11–27.
Nedoszytko B, Szczerkowska-Dobosz A, Stawczyk-Macieja M, Owczarczyk-Saczonek A, Reich A, Bartosinska J, et al. Pathogenesis of psoriasis in the “omic” era. Part II. Genetic, genomic and epigenetic changes in psoriasis. Postepy Dermatol Alergol. 2020;37(3):283–98.
Potaczek DP, Harb H, Michel S, Alhamwe BA, Renz H, Tost J. Epigenetics and allergy: from basic mechanisms to clinical applications. Epigenomics. 2017;9(4):539–71.
Chatterjee R, Vinson C. CpG methylation recruits sequence specific transcription factors essential for tissue specific gene expression. Biochim Biophys Acta. 2012;1819(7):763–70.
Rishi V, Bhattacharya P, Chatterjee R, Rozenberg J, Zhao J, Glass K, et al. CpG methylation of half-CRE sequences creates C/EBPalpha binding sites that activate some tissue-specific genes. Proc Natl Acad Sci USA. 2010;107(47):20311–6.
Kim YI, Logan JW, Mason JB, Roubenoff R. DNA hypomethylation in inflammatory arthritis: reversal with methotrexate. J Lab Clin Med. 1996;128(2):165–72.
Arnone M, Takahashi MDF, Carvalho AVE, Bernardo WM, Bressan AL, Ramos AMC, et al. Diagnostic and therapeutic guidelines for plaque psoriasis - Brazilian Society of Dermatology. An Bras Dermatol. 2019;94(2 Suppl 1):76–107.
Zhao Y, Jhamb D, Shu L, Arneson D, Rajpal DK, Yang X. Multi-omics integration reveals molecular networks and regulators of psoriasis. BMC Syst Biol. 2019;13(1):8.
Bowcock AM, Krueger JG. Getting under the skin: the immunogenetics of psoriasis. Nat Rev Immunol. 2005;5(9):699–711.
Furue M, Furue K, Tsuji G, Nakahara T. Interleukin-17A and Keratinocytes in Psoriasis. Int J Mol Sci. 2020;21(4):1275.
Hawkes JE, Yan BY, Chan TC, Krueger JG. Discovery of the IL-23/IL-17 signaling pathway and the treatment of psoriasis. J Immunol. 2018;201(6):1605–13.
Krueger JG, Wharton KA Jr, Schlitt T, Suprun M, Torene RI, Jiang X, et al. IL-17A inhibition by secukinumab induces early clinical, histopathologic, and molecular resolution of psoriasis. J Allergy Clin Immunol. 2019;144(3):750–63.
Yang BY, Cheng YG, Liu Y, Liu Y, Tan JY, Guan W, et al. Datura Metel L. Ameliorates imiquimod-induced psoriasis-like dermatitis and inhibits inflammatory cytokines production through TLR7/8-MyD88-NF-kappaB-NLRP3 inflammasome pathway. Molecules. 2019;24(11):2157.
Tervaniemi MH, Katayama S, Skoog T, Siitonen HA, Vuola J, Nuutila K, et al. NOD-like receptor signaling and inflammasome-related pathways are highlighted in psoriatic epidermis. Sci Rep. 2016;6:22745.
Buerger C, Richter B, Woth K, Salgo R, Malisiewicz B, Diehl S, et al. Interleukin-1beta interferes with epidermal homeostasis through induction of insulin resistance: implications for psoriasis pathogenesis. J Invest Dermatol. 2012;132(9):2206–14.
Albanesi C, De Pita O, Girolomoni G. Resident skin cells in psoriasis: a special look at the pathogenetic functions of keratinocytes. Clin Dermatol. 2007;25(6):581–8.
Nickoloff BJ, Xin H, Nestle FO, Qin JZ. The cytokine and chemokine network in psoriasis. Clin Dermatol. 2007;25(6):568–73.
Aydin B, Arga KY, Karadag AS. Omics-driven biomarkers of psoriasis: recent insights, current challenges, and future prospects. Clin Cosmet Investig Dermatol. 2020;13:611–25.
Luo Y, Luo Y, Chang J, Xiao Z, Zhou B. Identification of candidate biomarkers and pathways associated with psoriasis using bioinformatics analysis. Hereditas. 2020;157(1):30.
Xu M, Deng J, Xu K, Zhu T, Han L, Yan Y, et al. In-depth serum proteomics reveals biomarkers of psoriasis severity and response to traditional Chinese medicine. Theranostics. 2019;9(9):2475–88.
Aizawa N, Ishiuji Y, Tominaga M, Sakata S, Takahashi N, Yanaba K, et al. Relationship between the degrees of itch and serum lipocalin-2 levels in patients with psoriasis. J Immunol Res. 2019;2019:8171373.
Surace AEA, Hedrich CM. The role of epigenetics in autoimmune/inflammatory disease. Front Immunol. 2019;10:1525.
Chang KC, Shen Q, Oh IG, Jelinsky SA, Jenkins SF, Wang W, et al. Liver X receptor is a therapeutic target for photoaging and chronological skin aging. Mol Endocrinol. 2008;22(11):2407–19.
Schonthaler HB, Guinea-Viniegra J, Wculek SK, Ruppen I, Ximenez-Embun P, Guio-Carrion A, et al. S100A8-S100A9 protein complex mediates psoriasis by regulating the expression of complement factor C3. Immunity. 2013;39(6):1171–81.
Choudhary V, Uaratanawong R, Patel RR, Patel H, Bao W, Hartney B, et al. Phosphatidylglycerol inhibits toll-like receptor-mediated inflammation by danger-associated molecular patterns. J Invest Dermatol. 2019;139(4):868–77.
We thank Zhen Chen for advice and support on bioinformatics analysis.
This research was supported by the Military logistics scientific research program in 2020 (CLJ20J033).
Ethics approval and consent to participate
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wang, H., Chen, W., He, J. et al. Network analysis of potential risk genes for psoriasis. Hereditas 158, 21 (2021). https://doi.org/10.1186/s41065-021-00186-w
- Differentially expressed genes
- Differentially methylated positions