Next Article in Journal
Cell Death by Entosis: Triggers, Molecular Mechanisms and Clinical Significance
Next Article in Special Issue
Spatiotemporal Gene Expression Regions along the Anterior–Posterior Axis in Mouse Embryos before and after Palatal Elevation
Previous Article in Journal
The PDE4 Inhibitor Tanimilast Restrains the Tissue-Damaging Properties of Human Neutrophils
Previous Article in Special Issue
Low-Dose SARS-CoV-2 S-Trimer with an Emulsion Adjuvant Induced Th1-Biased Protective Immunity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Suphagan Boonpethkaew
1,
Jitlada Meephansan
1,*,
Onjira Jumlongpim
1,
Pattarin Tangtanatakul
2,
Wipasiri Soonthornchai
3,
Jongkonnee Wongpiyabovorn
4,
Ratchanee Vipanurat
5 and
Mayumi Komine
6
1
Division of Dermatology, Chulabhorn International College of Medicine, Rangsit Campus, Thammasat University, Klong Luang, Pathum Thani 12120, Thailand
2
Department of Transfusion Medicine and Clinical Microbiology, Faculty of Allied Health Sciences, Chulalongkorn University, Bangkok 10330, Thailand
3
School of Science, University of Phayao, Phayao 56000, Thailand
4
Center of Excellence in Immunology and Immune-Mediated Disease, Department of Microbiology, Faculty of Medicine, Chulalongkorn University, Bangkok 10330, Thailand
5
Division of Dermatology, Department of Medicine, Rajavithi Hospital, Ministry of Public Health, Bangkok 10400, Thailand
6
Department of Dermatology, Jichi Medical University, 3311-1 Yakushiji, Shimotsuke, Tochigi 329-0498, Japan
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(9), 4983; https://doi.org/10.3390/ijms23094983
Submission received: 18 April 2022 / Revised: 27 April 2022 / Accepted: 28 April 2022 / Published: 30 April 2022

Abstract

:
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.

1. Introduction

Psoriasis is a chronic skin disease influenced by the interplay of genetic predisposition, environmental triggers, and immune response. The classic feature of chronic plaque or psoriasis vulgaris is a well-demarcated salmon-pink plaque covered with silvery scales [1]. In addition to psoriasis vulgaris, inverse, guttate, pustular, and erythrodermic forms have been reported as the different phenotypes of psoriasis, which may be associated with comorbidities such as arthritis, metabolic syndrome, and cardiovascular diseases [1,2,3]. Various genes involved in antigen presentation (HLA-Cw6), skin barrier formation, innate immunity, and adaptive immunity have been identified as psoriasis-susceptible genes that increase the risk of developing the disease [4,5]. Individuals who carry these susceptible genes will be more likely to develop psoriasis under the influence of environmental factors, which initiate an inflammatory response leading to psoriasis phenotypes [1,5,6,7].
In the initial phase, injured keratinocytes release self-DNA/RNA and cathelicidin (LL-37), which form an LL-37/self-RNA/DNA complex that activates plasmacytoid dendritic cells (pDCs) via Toll-like receptor (TLR)7/9. The activated pDCs produce interferon (IFN)-α and -β to activate myeloid dendritic cells (mDCs). In addition to pDC-derived IFN-α, activated keratinocyte-derived interleukin (IL)-1β, IL-6, and tumor necrosis factor (TNF) can mature and activate myeloid dendritic cells (mDCs) to produce IL-23 [5,8]. Subsequently, the mDC- derived IL-23 drives the differentiation of T helper (Th)17 and T cytotoxic (c)17 cells. Both Th17 and Tc17 cells provide a substantial source of IL-17A and IL-22, which are the key cytokines for antimicrobial peptide (AMP) induction and keratinocyte (KC) activation to produce T cell- and neutrophil-chemokines and promote KC proliferation, ultimately resulting in the typical features of psoriasis. Thus, the IL-23/Th17 axis has been suggested as the main psoriasis driver [5,7,8]. In addition to the three main cells (DCs, T cells, and KCs), neutrophils, innate lymphoid cells (ILCs), Langerhans cells, macrophages, and natural killer cells could contribute to the pathogenesis of psoriasis [5,8]. Currently, there are many effective biologics available for treatment, including TNF inhibitors, IL-12 and IL-23 blockers, and anti-IL-17. However, the adverse effects of these drugs necessitate long-term monitoring, and there is currently no therapeutic option to cure psoriasis due to the complex and not yet fully elucidated pathogenesis [1,9,10].
The first RNA-sequencing (RNA-seq) analysis in psoriatic tissue was performed in 2012, which identified psoriasis-associated differentially expressed genes (DEGs), and has facilitated further studies to gain a better understanding of psoriasis pathogenesis [11]. Although the number of genes associated with psoriasis has been increasing, most of these have been identified from lesional (LS) skin samples [11,12,13,14,15,16,17,18,19,20,21]. Recently, transcriptomic studies in psoriasis have been focusing on more cell-specific or single-cell RNA-seq approaches; these approaches provide insight into psoriasis pathogenesis [22,23,24]. However, the peripheral edge of lesional (PE) skin RNA-seq is still unexplored. In the active edge lesion, psoriatic vascular change precedes epidermal and immune cell changes [25,26]. Furthermore, the interaction between dendritic cells and T cells is thought to initiate early inflammation and plaque propagation beneath the perilesional (PR) skin, a normal-appearing skin adjacent to the PE skin [27,28]. Hence, these observations suggested that the skin may have some molecules or signaling pathways that focus on intensifying inflammation in the transitional zone, which might contribute to plaque propagation. We here present the first transcriptomic profile of the PE skin along with a proposed network. These findings could provide new insight to expand the molecular network contributing to the pathogenesis of psoriasis toward identifying new therapeutic options for patients with psoriasis vulgaris.

2. Results

2.1. DEGs in the PE Skin

A total of 1202 DEGs were identified. Of these, 653 (54%) were upregulated, and 549 (46%) were downregulated in PE skin compared to the uninvolved (UN) skin (Figure 1). A heatmap of the differentially regulated transcripts between the PE and UN skin is shown in Supplementary Figure S1. The top 50 upregulated and downregulated DEGs according to the log2FC value are listed in Table 1. The top 100 DEGs are shown in Supplementary Tables S1 and S2.

2.2. DEG Analysis

2.2.1. Disease Associations and Biological Functions

The 1202 DEGs were analyzed for downstream effects in terms of associated diseases; among them, 135 genes have previously been reported to be associated with psoriasis based on the Ingenuity Knowledge Base, 86 of which showed upregulated in the PE skin. Thus, we hypothesized that the proteins encoded by these 86 DEGs might be involved in inflammation in the PE skin to contribute to psoriasis (Figure 1 and Supplementary Table S3). Furthermore, the 1202 DEGs were analyzed for their functions. These DEGs enhanced functions are to enrich angiogenesis (EDNRB, LAMA2, S100A8, S100A9, SPINK5, TNFRSF1A, and VAV3), growth of epithelial tissue (ADAM17, KRT16, IKBKB, IL36G, ITGB1, RICTOR, and XDH), chemotaxis and homing of cells (A2M, C10orf99, DEFB4A/DEFB4B, GFRA1, PLEC, S100A7, and S100A7A), growth of connective tissue (ADAR, ELN, IL36A, KLK6, PLAT, STAT1, and STAT3) and degranulation of myeloid cells (CD59, CXCL8, DSC1, DSG1, LCN2, and SERPINB3). The list of all DEGs in each function is provided in Supplementary Table S4.
In addition, the upregulated- and downregulated DEGs were analyzed separately as pathway and function enrichments using Metascape. Notably, cornified envelope formation, VEGFA-VEGFR2 signaling, interferon signaling, and neutrophil degranulation were enriched by upregulated DEGs (Figure 2A), whereas epithelial differentiation was enriched by downregulated DEGs (Figure 2B) such as KRT10, KRT2, LCE1C, LCE1D, LCE2B, and LORICRIN. The list of all DEGs in each DEG-enriched pathway and function is provided in Supplementary Tables S5 and S6. We reviewed and selected the potential DEGs based on their biofunctions and psoriasis associations; they are shown in Table 2 and Figure 3. Moreover, we compared and mapped PE vs. UN skin DEGs from previous studies [11,12,13,14,22,23,24,29,30]; the details are provided in Supplementary Table S7.

2.2.2. Canonical Pathways

We analyzed the known canonical pathways enriched by the DEGs of PE skin and the activity prediction of these pathways based on the QIAGEN Knowledge Base. To include both significant and highly predicted activity pathways, we determined the threshold of minimum significance (p-value of overlap > 0.05 or −log (p-value) > 1.3) and absolute z-score (|z-score| ≥ 2). Only seven canonical pathways are shown in Figure 4A,B. It is noteworthy that the role of IL-17A in psoriasis, eukaryotic initiation factor 2 (eIF2) signaling, and the coronavirus pathogenesis pathway had the top three highest ratios (mapped DEGs divided by total molecules in a given pathway = 0.36, 0.30, and 0.25, respectively). Moreover, these three pathways are involved in the stress response and inflammation, especially the coronavirus pathogenesis pathway, the hallmark of which is the cytokine storm [31]. The ratios of all seven pathways are listed in Supplementary Table S8.
To explore the role of cytokine signaling in the PE skin in further detail, we analyzed the canonical pathway of cytokines, with a −log (p-value) > 1.3 indicating significance. A total of 11 cytokine signaling pathways were among the DEG-enriched pathways (Figure 4C,D), including the IL-8, IL-17A, IL-6, C-C motif chemokine receptor 3(CCR3), IL-15, granulocyte-macrophage colony-stimulating factor (GM-CSF), IL-22, and formyl-methionyl-leucyl-phenylalanine (fMLP) signaling pathways.

2.2.3. Upstream Regulators

We further identified the upstream regulators of the DEGs in the PE skin. We used the upstream analysis function according to a statistically significant overlap between the dataset genes and the genes that are regulated by a potential regulator (p-value overlap < 0.05) and |z-score| ≥ 2. The lists of potential regulators with their predicted activity are provided in Supplementary Table S9. Among the potential regulators, seven molecules were identified as DEGs, including CCND1, DDX21, IL36A, KDM5A, NUPR1, RICTOR, and STAT1. Some relationships between the potential regulators and their target DEGs are shown in Figure 5 and Supplementary Figure S2.

2.2.4. Regulator Effect

Next, we established a causal hypothesis between upstream regulators and their downstream effects using the regulator effect function to correlate the significant upstream regulators with their downstream functions. The useful predicted networks were prioritized and ranked according to their consistency scores. The two highest consistency score networks are shown in Figure 6. Notably, IL-36α, IL-17C, nuclear factor-kappa B inhibitor zeta (IκBζ), and Rho-associated kinase 2 (ROCK2) emerged as potentially important regulators in the network.

2.2.5. Validation of Selected DEGs by Multiplex Real-Time PCR

Finally, to validate the expression of DEGs such as SERPINB4, IL36A, PLSCR1, CXCL8, DEFB4A/B, DEFB103A/B, S100A7, S100A7A, S100A8, and S100A9, we performed quantitative multiplex Real-Time PCR. Furthermore, we compared the mRNA expression of upstream regulators IL-17C, IκBζ, and ROCK2 in the PE skin to those in the UN skin. The results are presented in Figure 7.

3. Discussion

To the best of our knowledge, this is the first study to establish the transcriptomic profile at the edge of psoriatic lesions. We revealed that the increased functions in the PE skin were involved in angiogenesis, growth of epithelial tissue, chemotaxis and homing of cells, growth of connective tissue, and the degranulation of myeloid cells. The upregulated DEGs were enriched in functions such as cornified envelope formation, VEGFA-VEGFR2 signaling, interferon signaling, and neutrophil degranulation. Moreover, we showed that the genes involved in epithelial differentiation were downregulated in the PE skin compared with the UN skin, which might limit the normal differentiation capacity of KC and contribute to plaque formation. We discussed some potential DEGs that were related to AMPs, angiogenesis, autoantigens, and chemotaxis that might play a pivotal role in the pathogenesis of PE skin.
AMPs activate the innate immune response to create an inflammatory environment in the skin through antimicrobial, chemotactic, angiogenic, and proinflammatory properties [32,33]. We identified various AMP-coding genes that were upregulated in the PE skin. We confirmed that DEFB4A/B mRNA was significantly upregulated in PE skin and its fold change was the highest among the AMP-coding genes. It was previously reported to be expressed in psoriatic plaques [34,35] and more highly in lesional skin than in UN skin [36]. Its protein, human β-defensin 2 (hBD-2), promotes pDCs to uptake self-DNA and secrete IFN-α in a TLR-9–dependent manner [37]. Moreover, we confirmed that S100A7, S100A7A, S100A8, and S100A9 mRNA were significantly upregulated in PE skin compared to those in UN skin. Interestingly, among the psoriasis-associated S100 protein-coding genes, S100A7A (also known as S100A15) was the most upregulated gene in the PE skin. Its protein koebnerisin is expressed in psoriatic lesions along with psoriasin (encoded by S100A7), S100A8, and S100A9 [38,39] by induction of psoriatic-derived Th17/Th22 cytokines (IL-17, IL-22, and TNF-α) [40,41,42]. Of particular interest, apart from their antimicrobial properties of them [33,43], both koebnerisin and psoriasin have been proposed to be responsible for the inflammation priming of KCs [41], as well as S100A8 and S100A9 also function as proangiogenic factors [39]. However, further study is needed to determine their role in the PE skin, in particular the very-high-fold change of S100A7A.
We also identified the upregulation of SERPINB3 and SERPINB4 in the PE skin and confirmed that SERPINB4 was significantly upregulated in the PE skin. Their proteins, serpin family B members 3 and 4, were reported to be expressed in psoriatic lesions [44] by induction of IL-17 and IL-22 [45]. Pso p27 antigen, a proteolyzed product of SERPINB4 and SERPINB3 [46,47], was detected in the PE skin, which is considered to be equivalent to its detection on the lesion, but was not detected on the UN skin [48]. The pso p27 forms an immune complex with anti-pso p27 to activate the complement cascade at the psoriatic scale [49]. However, the roles of SERPINB3 and SERPINB4 as antigens in psoriasis are controversial and require further investigation, especially SERPINB4, which was highly upregulated in the PE skin.
We further used the canonical pathway, regulator analysis, and regulator effect functions of QIAGEN’s IPA tool to generate the network among the molecules. The proposed pathogenic network is schematically displayed in Figure 8. We found that IL36A and IL36G expression were upregulated in the PE skin and then confirmed that IL36A expression was significantly upregulated in the PE skin. Their proteins IL-36α and IL-36γ are expressed in psoriatic plaque [50,51,52], and their mRNA levels were reported to be increased in KCs after injury [53,54]. Moreover, IL-36α and IL-36γ alone or synergistically with IL-17A could amplify their own mRNA and proteins in KCs [55]. We identified and confirmed the significant upregulation of PLSCR1 (encoding a TLR-9 transporter) in the PE skin. A previous RNA-seq analysis of whole blood from psoriasis patients implicated the pDC-mediated release of IFN-α after induction by IL-36 in a TLR-9-dependent manner, evidenced by the upregulation of PLSCR1 in pDCs [56]. The pDC-derived IFN-α promotes mDC maturation and function, which amplifies and sustains T cells for the psoriatic inflammatory response [7,8]. Therefore, we suggest that the IL-36/TLR-9 axis, which upregulates IFN-α production, might also drive psoriatic inflammation in the PE skin (Figure 8).
IL-36α and IL-36γ also activate mDCs to mature and secrete IL-6 [57]. Dermal DC-derived IL-6 then drives Th17 and Th22 cell differentiation [7,58]. Furthermore, IL-6 is also secreted by KCs after IFN-γ treatment via the suppression of miR-149, which amplifies skin inflammation [59]. Consistently, we found that the DEGs in the PE skin were enriched in the IL-6 canonical pathway, suggesting that the IL-6/Th17 axis might be another axis driving psoriatic inflammation in the PE skin (Figure 8). However, despite the lack of IL-6, KCs expressed various additional proinflammatory cytokines after the transient suppression of psoriasiform skin manifestation, delineating the ineffective treatment of psoriasis plaque by IL-6 blockage [60]. Our results showed that the DEGs were enriched in the IL-17A and IL-22 canonical pathways and that IL-22 was an activated potential regulator, thus suggesting the crucial roles of IL-17A and IL-22 in the PE skin. IL-17A and IL-22 potentiate the psoriatic milieu by inducing AMPs, recruiting other immune cells, and increasing KC proliferation [7,8]. In addition, a synergistic effect of the Th17 cytokines IL-17A and IL-22 induces the expression of IL-36α and IL-36γ mRNA in KCs [50,61]. IL-22, in particular, promotes KC proliferation and suppresses KC terminal differentiation [42]. Previous studies comparing LS skin to healthy skin suggested that IL-22 functions better when the levels of soluble scavenging receptor IL-22 binding protein (encoded by IL22RA2) are decreased. However, we could not identify the alteration of IL22RA2 as DEG of PE vs. UN skin. Furthermore, serum IL-22 level was higher in psoriasis patients compared with healthy individuals and was also positively correlated with disease severity [42,62].
Considering IL-36, IL-36α was identified as an activated potential regulator in the PE skin, which induces the expression of IL36G, IL17C, S100A9, and NFKBIZ in KCs [63]. Notably, the mouse skin showed enhanced chemokine expression, leukocyte infiltration, and acanthosis after the injection of IL-36α [57]. Moreover, imiquimod induced much milder psoriasiform dermatitis in Il36a−/− mice compared with that induced in wild-type mice [64].
IL-17C was identified as an activated potential regulator. We established that its mRNA expression was not significantly upregulated in the PE skin because similar relative expression levels were also detected in the UN skin. These results suggest that IL-17C plays a role during the early stages of psoriasis inflammation or inflammatory priming for plaque formation. A previous report demonstrated that IL-17C transgenic mice developed psoriasiform skin. Moreover, IL-17C and its mRNA significantly increase in UN and LS skin compared to healthy skin [65]. IL-17C is mainly secreted from KCs and not leukocytes [66]. It has also been detected in psoriatic plaques at a higher concentration than IL-17A [65,67]. IL-17C stimulates the release of TNF-α from monocytes [68] and also simulates neutrophil migration, especially in synergy with TNF-α [67]. IL-17C could amplify itself in KCs, especially in synergy with TNF-α [65,67]. In addition, IL-36α and IL-36γ alone or synergistically with IL-17A increased IL-17C levels from KCs [55]. Moreover, IL-17C induced the expression of S100A15, S100A7, S100A8, S100A9, DEFB4A, CXCL8, IL36G, NFKBIZ, and TNIP3 in KCs [66,67]. In particular, IL-17C synergistically with TNF-α induced the expression of S100A7, S100A8, S100A9, LCN2 (coding lipocalin-2), CXCL8, IL36G, and DEFB4 in KCs [65].
IκBζ was identified as an activated potential regulator in the PE skin network. Moreover, we also observed that NFKBIZ mRNA was significantly upregulated in the PE skin. A recent study reported that NFKBIZ was co-expressed with IL36G in the stratum spinosum of psoriasis epidermis [24]. IκBζ is a nuclear transcriptional cofactor that plays a role as the mediator of the expression of IL-36α-induced genes, including S100A7, S100A8, S100A9, LCN2, DEFB4, IL36G, SERPINB4, and CXCL8 [63], IL-17A-induced genes, including DEFB4, S100A7, NFKBIZ, and LCN2 [69,70], and IL-17A/TNF-α-induced gene, including DEFB4, LCN2, CXCL8, IL17C, NFKBIZ, and IL36G [69,70,71,72] in KCs.
We identified significant upregulation of CXCL8 in the PE skin. It encodes IL-8, which attracts T lymphocytes to the psoriasis plaque [73], promotes angiogenesis [74,75], and promotes epidermal proliferation [76]. Interestingly, IL-8 is secreted by accumulated neutrophils (Munro’s microabscess) and is the chemokine for other neutrophils to aggregate into the plaque in an autocrine manner [77]. Moreover, IL-8 is secreted from KCs after scratching [54] or the induction of IL-22 or induction of IL-36α and IL-36γ, especially in synergy with IL-17A [55,57,61,78]. Particularly, KC-derived IL-8 is capable of neutrophil chemotaxis, but this capacity was reduced when miR-146a, an IL-17 counteractant, which is expressed in psoriasis plaque, was overexpressed [79].
We identified that ROCK2, as an activated potential regulator, could be equally important as IL-36α, IL-17C, and IκBζ in the network. ROCK2 is located in the nucleus and cytoplasm of the human T- cell [80]. It plays an important role in keratinocyte differentiation [81] as well as in angiogenesis [82]. Moreover, ROCK2 was found to regulate CD4+ T cells to secrete IL-17 [83]. These processes are important for the formation of psoriasis plaques [7]. In 2017, the activity of ROCK2 was studied in moderately severe plaque psoriasis patients using selective oral ROCK2 inhibitors; 71% of the patients achieved a reduction of 50% in the psoriasis area and severity index (PASI 50) at 12 weeks, peripheral blood IL-17 and IL-23 significantly decreased, and ROCK2 staining decreased with decreasing epidermal thickness [84]. As the predicted regulator in the PE skin, ROCK2 was associated with the upregulation of nine DEGs (CXCL8, DSC1, DSC2, DSG1, DSG3, SPINK5, PI3, S100A7, and S100A8) (Figure 5). We further investigated the ROCK2 mRNA expression in PE skin and compared it to those in the UN skin. Our results showed that ROCK2 mRNA was significantly downregulated. However, activation of the ROCK2 kinase protein depends on post-translational modification steps. Previous reports suggest that activation of ROCK2 kinase needs phosphorylation [85,86]. Interestingly, its structure is capable of self-phosphorylation (auto-phosphorylation) [87]. Thus, we suggest that the downregulation of ROCK2 mRNA in PE skin might be the result of autoregulation overfunction. In addition, ROCK2 kinase has been proposed to cooperate with phosphorylated signal transducer and activator of transcription 3 (STAT3) to regulate Th17 function [80]. Due to the complexity of psoriasis, further investigation of ROCK2 mRNA expression in the PE skin compared to the mature lesion or healthy skin could provide new information toward a better understanding of the mechanism of oral ROCK2 inhibitor action and the role of ROCK2 in the pathogenesis of psoriasis.
Finally, based on our results and these previous findings, we propose that IL-36α, IL-36γ, IκBζ, IL-17C, IL-8, S100A7, S100A8, S100A9, S100A15, SERPINB4, hBD-2, and ROCK2 potentiate inflammation in the PE skin and that IL-36α and IL-36γ may run a self-sustaining circuit between the innate and adaptive immune response, creating a psoriatic milieu in the PE skin (Figure 8).
The limitation of this study is the small sample size. Ideally, more samples of a larger size need to be evaluated to determine the PE skin transcriptomic profile.
In conclusion, this study identified a set of upregulated DEGs, including PLSCR1, IL36A, CXCL8, S100A7, S100A8, S100A9, S100A15, DEFB4A/B, and SERPINB4, and the molecules that may regulate these DEGs including IL-17C, IL-36α, IκBζ, and ROCK2. These DEGs and upstream regulators may potentiate the inflammation in PE skin as the proposed network might be different from that seen in the LS skin. This study provides new insight into the molecular mechanisms occurring in the PE skin and could lead to the identification of specific molecules that could serve as therapeutic targets in psoriasis vulgaris. Further, it would be helpful to compare these molecules with those observed in the center of LS skin, normal-appearing skin in patients with psoriasis, and healthy skin, to further confirm changes in the levels of other DEGs and regulators of interest at both the mRNA and protein levels, and to further examine therapeutic modalities that might target these molecules by comparing the before and after treatment transcript levels in the PE skin. These insights could provide a better understanding of the pathogenesis of psoriasis and lead to the discovery of novel therapeutic strategies.

4. Materials and Methods

4.1. Patients

We enrolled voluntary patients with plaque-type psoriasis (aged ≥18 years) with moderate-to-severe psoriasis [Psoriasis Area and Severity Index (PASI) score ≥ 10] diagnosed by either classical clinical features or a confirmed pathological examination. Patients who had been treated with topical agents (steroids, vitamin D analogs) within 2 weeks or had undergone phototherapy and systemic medication (prednisolone, methotrexate, biologics) within 4 weeks of the date of tissue sampling were excluded. Based on previous RNA-seq profiling studies in human psoriasis skin, three samples were deemed to be adequate for analysis [11,88]; therefore, three eligible patients were included in this study. The PASI scores of the three patients are shown in Supplementary Table S10. All three patients were informed of all study protocols and signed the consent forms. The protocol was designed in accordance with the principles of the Declaration of Helsinki and was approved by the ethics committee of Thammasat University, Thailand.

4.2. Tissue Sampling

Full-thickness skin biopsies were collected from the three patients using a 6 mm punch biopsy under local anesthesia (2% lidocaine with adrenaline). The PE skin was taken at the edge of the plaque (Figure 9), and UN skin was collected from phenotypically normal skin a minimum of 10 cm away from the plaque from the same patient. The biopsies were preserved in RNAlater (ThermoFisher, AMBION, Austin, TX, USA) at −80 °C until gene expression profiling.

4.3. High-Throughput Next-Generation Sequencing (NGS)

Total RNA was extracted from each skin biopsy of the volunteers (PE skin, n = 3, UN skin, n = 2) using TRIzol Reagent (Ambion, Austin, TX, USA). NGS was performed by Macrogen (Seoul, Korea). Briefly, the purity and structural integrity of RNA were analyzed using an Agilent 2100 Bioanalyzer (RNA integrity number > 7). Sequence libraries were constructed using the SMARTer Universal Low Input RNA Kit and TruSeq RNA Sample Prep Kit v2 (pair-end). NGS was performed using the NovaSeq 6000 S4 Reagent Kit on the NovaSeq 6000 System (Illumina Inc., San Diego, CA, USA) with 100-bp paired-end reads. The raw data generated in this study have been deposited in the NCBI Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo, accessed on 8 September 2021) and are accessible through the GEO Series accession number GSE183732.

4.4. Sequencing Data Analysis

Basic data analysis was conducted by Macrogen. Briefly, overall read qualities, total bases, total reads, guanine-cytosine content (%GC content), and basic statistics were calculated. Trimmed reads were mapped to the reference genome using HISAT2, a splice-aware aligner. The transcripts were assembled using StringTie with aligned reads. Expression profiles were repeated as reading counts and normalized based on the transcript length and depth of coverage. The fragments per kilobase of transcript per million mapped reads (FPKM) value was used as a normalization value. Furthermore, we used the read counts to statistically analyze the expression profile with edgeR Bioconductor statistical library version 3.13 [89,90] on R Studio [91]. The DEGs between the PE skin and UN skin were determined according to p < 0.01 (0.026 ≤ FDR ≤ 0.194) and fold change (FC) ratio (|log2FC|) ≥ 1.5.

4.5. DEGs Analysis

All DEGs with expression profiles were uploaded into QIAGEN’s Ingenuity Pathway Analysis (IPA) software (QIAGEN Redwood City, CA, USA; www.qiagen.com/ingenuity, accessed on 1 October 2021) using the “core analysis” of DEGs. Disease and biofunction, canonical pathways, upstream regulators, and regulator effects of the DEGs were analyzed based on the Ingenuity Knowledge Base. Pathway and function enrichment analysis was additionally analyzed using Metascape Online Tool (https://metascape.org, accessed on 1 October 2021) [92].

4.6. Quantitative Multiplex Real-Time PCR (Multiplex qPCR)

For validation experiments, UN and PE skin samples were biopsied from moderate to severe psoriasis vulgaris patients. RNA was extracted and reverse transcribed to complementary DNA using the ImProm-IITM Reverse Transcription Kit (Promega Corp., Madison, WI, USA). TagMan® Gene Expression Assay predesigned qPCR primers of the selected DEGs (IL36A, CXCL8, PLSCR1, SERPINB4, DEFB4A/B, DEFB103A/B, S100A7, S100A7A, S100A8, and S100A9), upstream regulators (IL17C, NFKBIZ, and ROCK2), and of housekeeping gene (GAPDH) for qPCR were purchased from Gene Plus Co., Ltd. (ThermoFisher Scientific, West Sacramento, CA, USA). Their assay IDs are listed in Supplementary Table S11. The multiplex qPCR was performed using the QuantiStudioTM 6 Flex Real-Time PCR system (Applied Biosystem). A statistical difference in the mean of ∆Ct (targeted gene vs GAPDH) between the UN and PE skin groups was calculated by a paired t-test. Graphs created by GraphPad Prism 9.3.1 were shown for relative mRNA expression using the 2−∆Ct method and for fold change using the 2−∆∆Ct method [93].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23094983/s1, References [11,12,13,14,22,23,24,29,30] are cited in the supplementary materials.

Author Contributions

S.B. wrote all sections of the manuscript. J.M. contributed to the conception and design of the study. S.B., O.J. and R.V. collected the data. S.B., P.T. and W.S. contributed to the data and performed the analysis. J.M., J.W. and M.K. contributed to manuscript revisions. All authors have read and agreed to the published version of the manuscript.

Funding

Funding was provided by Thammasat University Research Fund under the TU Research Scholar (contract No. 1/2560) and Chulabhorn International College of Medicine (contract No. T1/2563).

Institutional Review Board Statement

This study was approved by the Human Ethics Committee of Thammasat University No. MTU-EC-OO-2-027/60, conducted according to the principles of the Declaration of Helsinki.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The raw data generated in this study have been deposited in the NCBI Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo (accessed on 1 March 2022)) and are accessible through the GEO Series accession number GSE183732.

Acknowledgments

The authors gratefully acknowledge Poom Pansritum, who draws excellent figures and Pannaphan Makarathut, who assisted in the multiplex qPCR lab.

Conflicts of Interest

The authors have no conflict of interest to declare.

References

  1. Griffiths, C.E.M.; Armstrong, A.W.; Gudjonsson, J.E.; Barker, J.N.W.N. Psoriasis. Lancet 2021, 397, 1301–1315. [Google Scholar] [CrossRef]
  2. Dastoli, S.; Nisticò, S.P.; Morrone, P.; Patruno, C.; Leo, A.; Citraro, R.; Gallelli, L.; Russo, E.; De Sarro, G.; Bennardo, L. Colchicine in Managing Skin Conditions: A Systematic Review. Pharmaceutics 2022, 14, 294. [Google Scholar] [CrossRef] [PubMed]
  3. Amoruso, G.; Nisticò, S.; Iannone, L.; Russo, E.; Rago, G.; Patruno, C.; Bennardo, L. Ixekizumab May Improve Renal Function in Psoriasis. Healthcare 2021, 9, 543. [Google Scholar] [CrossRef] [PubMed]
  4. Harden, J.L.; Krueger, J.G.; Bowcock, A.M. The immunogenetics of Psoriasis: A comprehensive review. J. Autoimmun. 2015, 64, 66–73. [Google Scholar] [CrossRef] [Green Version]
  5. Di Meglio, P.; Villanova, F.; Nestle, F.O. Psoriasis. Cold Spring Harb. Perspect. Med. 2014, 4, a015354. [Google Scholar] [CrossRef] [Green Version]
  6. Zeng, J.; Luo, S.; Huang, Y.; Lu, Q. Critical role of environmental factors in the pathogenesis of psoriasis. J. Dermatol. 2017, 44, 863–872. [Google Scholar] [CrossRef] [Green Version]
  7. Perera, G.K.; Di Meglio, P.; Nestle, F.O. Psoriasis. Annu. Rev. Pathol. Mech. Dis. 2012, 7, 385–422. [Google Scholar] [CrossRef]
  8. Kim, J.; Krueger, J.G. The Immunopathogenesis of Psoriasis. Dermatol. Clin. 2015, 33, 13–23. [Google Scholar] [CrossRef]
  9. Hawkes, J.E.; Chan, T.C.; Krueger, J.G. Psoriasis pathogenesis and the development of novel targeted immune therapies. J. Allergy Clin. Immunol. 2017, 140, 645–653. [Google Scholar] [CrossRef] [Green Version]
  10. Iannone, L.F.; Bennardo, L.; Palleria, C.; Roberti, R.; De Sarro, C.; Naturale, M.D.; Dastoli, S.; Donato, L.; Manti, A.; Valenti, G.; et al. Safety profile of biologic drugs for psoriasis in clinical practice: An Italian prospective pharmacovigilance study. PLoS ONE 2020, 15, e0241575. [Google Scholar] [CrossRef]
  11. Jabbari, A.; Suárez-Fariñas, M.; Dewell, S.; Krueger, J.G. Transcriptional Profiling of Psoriasis Using RNA-seq Reveals Previously Unidentified Differentially Expressed Genes. J. Investig. Dermatol. 2012, 132, 246–249. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Li, B.; Tsoi, L.C.; Swindell, W.R.; Gudjonsson, J.E.; Tejasvi, T.; Johnston, A.; Ding, J.; Stuart, P.E.; Xing, X.; Kochkodan, J.J.; et al. Transcriptome Analysis of Psoriasis in a Large Case–Control Sample: RNA-Seq Provides Insights into Disease Mechanisms. J. Investig. Dermatol. 2014, 134, 1828–1838. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Swindell, W.R.; Xing, X.; Voorhees, J.J.; Elder, J.T.; Johnston, A.; Gudjonsson, J.E. Integrative RNA-seq and microarray data analysis reveals GC content and gene length biases in the psoriasis transcriptome. Physiol. Genom. 2014, 46, 533–546. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Ahn, R.; Yan, D.; Chang, H.-W.; Lee, K.; Bhattarai, S.; Huang, Z.-M.; Nakamura, M.; Singh, R.; Afifi, L.; Taravati, K.; et al. RNA-seq and flow-cytometry of conventional, scalp, and palmoplantar psoriasis reveal shared and distinct molecular pathways. Sci. Rep. 2018, 8, 11368. [Google Scholar] [CrossRef] [Green Version]
  15. Yu, Z.; Gong, Y.; Cui, L.; Hu, Y.; Zhou, Q.; Chen, Z.; Yu, Y.; Chen, Y.; Xu, P.; Zhang, X.; et al. High-throughput transcriptome and pathogenesis analysis of clinical psoriasis. J. Dermatol. Sci. 2020, 98, 109–118. [Google Scholar] [CrossRef]
  16. Kõks, S.; Keermann, M.; Reimann, E.; Prans, E.; Abram, K.; Silm, H.; Kõks, G.; Kingo, K. Psoriasis-Specific RNA Isoforms Identified by RNA-Seq Analysis of 173,446 Transcripts. Front. Med. 2016, 3, 46. [Google Scholar] [CrossRef] [Green Version]
  17. Zolotarenko, A.; Chekalin, E.; Mesentsev, A.; Kiseleva, L.; Gribanova, E.; Mehta, R.; Baranova, A.; Tatarinova, T.; Piruzian, E.S.; Bruskin, S. Integrated computational approach to the analysis of RNA-seq data reveals new transcriptional regulators of psoriasis. Exp. Mol. Med. 2016, 48, e268. [Google Scholar] [CrossRef] [Green Version]
  18. Ahn, R.; Gupta, R.; Lai, K.; Chopra, N.; Arron, S.T.; Liao, W. Network analysis of psoriasis reveals biological pathways and roles for coding and long non-coding RNAs. BMC Genom. 2016, 17, 841. [Google Scholar] [CrossRef] [Green Version]
  19. Suárez-Fariñas, M.; Li, K.; Fuentes-Duculan, J.; Hayden, K.; Brodmerkel, C.; Krueger, J.G. Expanding the Psoriasis Disease Profile: Interrogation of the Skin and Serum of Patients with Moderate-to-Severe Psoriasis. J. Investig. Dermatol. 2012, 132, 2552–2564. [Google Scholar] [CrossRef] [Green Version]
  20. Zolotarenko, A.; Chekalin, E.; Mehta, R.; Baranova, A.; Tatarinova, T.V.; Bruskin, S. Identification of Transcriptional Regulators of Psoriasis from RNA-Seq Experiments. Methods Mol. Biol. 2017, 1613, 355–370. [Google Scholar] [CrossRef]
  21. De Oliveira, P.S.S.; Pereira, M.C.; De Paula, S.K.S.; Lima, E.V.A.; Lima, M.M.D.A.; De Arruda, R.G.; De Oliveira, W.L.M.; Duarte, L.B.P.; Pitta, I.D.R.; Rêgo, M.J.M.B.; et al. Increased IL17A, IFNG, and FOXP3 Transcripts in Moderate-Severe Psoriasis: A Major Influence Exerted by IL17A in Disease Severity. Mediat. Inflamm. 2016, 2016, 4395276. [Google Scholar] [CrossRef] [PubMed]
  22. Pasquali, L.; Srivastava, A.; Meisgen, F.; Mahapatra, K.; Xia, P.; Landén, N.; Pivarcsi, A.; Sonkoly, E. The Keratinocyte Transcriptome in Psoriasis: Pathways Related to Immune Responses, Cell Cycle and Keratinization. Acta Derm. Venereol. 2019, 99, 196–205. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Swindell, W.R.; Sarkar, M.; Liang, Y.; Xing, X.; Baliwag, J.; Elder, J.T.; Johnston, A.; Ward, N.L.; Gudjonsson, J.E. RNA-seq identifies a diminished differentiation gene signature in primary monolayer keratinocytes grown from lesional and uninvolved psoriatic skin. Sci. Rep. 2017, 7, 18045. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Kim, J.; Lee, J.; Kim, H.J.; Kameyama, N.; Nazarian, R.; Der, E.; Cohen, S.; Guttman-Yassky, E.; Putterman, C.; Krueger, J.G. Single-cell transcriptomics applied to emigrating cells from psoriasis elucidate pathogenic versus regulatory immune cell subsets. J. Allergy Clin. Immunol. 2021, 148, 1281–1292. [Google Scholar] [CrossRef] [PubMed]
  25. Goodfield, M.; Hull, S.M.; Holland, D.; Roberts, G.; Wood, E.; Reid, S.; Cunliffe, W. Investigations of the ‘active’ edge of plaque psoriasis: Vascular proliferation precedes changes in epidermal keratin. Br. J. Dermatol. 1994, 131, 808–813. [Google Scholar] [CrossRef] [PubMed]
  26. Hull, S.M.; Goodfield, M.; Wood, E.J.; Cunliffe, W.J. Active and Inactive Edges of Psoriatic Plaques: Identification by Tracing and Investigation by Laser-Doppler Flowmetry and Immunocytochemical Techniques. J. Investig. Dermatol. 1989, 92, 782–785. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Komine, M.; Karakawa, M.; Takekoshi, T.; Sakurai, N.; Minatani, Y.; Mitsui, H.; Tada, Y.; Saeki, H.; Asahina, A.; Tamaki, K. Early Inflammatory Changes in the “Perilesional Skin” of Psoriatic Plaques: Is there Interaction between Dendritic Cells and Keratinocytes? J. Investig. Dermatol. 2007, 127, 1915–1922. [Google Scholar] [CrossRef] [Green Version]
  28. Alshenawy, H.A.; Hasby, E.A. Immunophenotyping of dendritic cells in lesional, perilesional and distant skin of chronic plaque psoriasis. Cell. Immunol. 2011, 269, 115–119. [Google Scholar] [CrossRef]
  29. Nattkemper, L.A.; Tey, H.L.; Valdes-Rodriguez, R.; Lee, H.; Mollanazar, N.K.; Albornoz, C.; Sanders, K.M.; Yosipovitch, G. The Genetics of Chronic Itch: Gene Expression in the Skin of Patients with Atopic Dermatitis and Psoriasis with Severe Itch. J. Investig. Dermatol. 2018, 138, 1311–1317. [Google Scholar] [CrossRef] [Green Version]
  30. Swindell, W.R.; Sarkar, M.; Liang, Y.; Xing, X.; Gudjonsson, J.E. Cross-Disease Transcriptomics: Unique IL-17A Signaling in Psoriasis Lesions and an Autoimmune PBMC Signature. J. Investig. Dermatol. 2016, 136, 1820–1830. [Google Scholar] [CrossRef] [Green Version]
  31. Chen, R.; Lan, Z.; Ye, J.; Pang, L.; Liu, Y.; Wu, W.; Qin, X.; Guo, Y.; Zhang, P. Cytokine Storm: The Primary Determinant for the Pathophysiological Evolution of COVID-19 Deterioration. Front. Immunol. 2021, 12, 589095. [Google Scholar] [CrossRef] [PubMed]
  32. Büchau, A.S.; Gallo, R.L. Innate immunity and antimicrobial defense systems in psoriasis. Clin. Dermatol. 2007, 25, 616–624. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Takahashi, T.; Yamasaki, K. Psoriasis and Antimicrobial Peptides. Int. J. Mol. Sci. 2020, 21, 6791. [Google Scholar] [CrossRef] [PubMed]
  34. de Jongh, G.J.; Zeeuwen, P.L.; Kucharekova, M.; Pfundt, R.; van der Valk, P.G.; Blokx, W.; Dogan, A.; Hiemstra, P.S.; van de Kerkhof, P.C.; Schalkwijk, J. High Expression Levels of Keratinocyte Antimicrobial Proteins in Psoriasis Compared with Atopic Dermatitis. J. Investig. Dermatol. 2005, 125, 1163–1173. [Google Scholar] [CrossRef]
  35. Steinz, K.; Schubert, S.; Harder, J.; Gerdes, S.; Mrowietz, U.; Gläser, R. Bacterial soft tissue infection in psoriasis despite induction of epidermal antimicrobial peptides. Exp. Dermatol. 2014, 23, 862–864. [Google Scholar] [CrossRef]
  36. Bierkarre, H.; Harder, J.; Cuthbert, R.; Emery, P.; Leuschner, I.; Mrowietz, U.; Hedderich, J.; McGonagle, D.; Gläser, R. Differential expression of antimicrobial peptides in psoriasis and psoriatic arthritis as a novel contributory mechanism for skin and joint disease heterogeneity. Scand. J. Rheumatol. 2015, 45, 188–196. [Google Scholar] [CrossRef]
  37. Tewary, P.; De La Rosa, G.; Sharma, N.; Rodriguez, L.G.; Tarasov, S.G.; Howard, O.M.Z.; Shirota, H.; Steinhagen, F.; Klinman, D.M.; Yang, D.; et al. β-Defensin 2 and 3 Promote the Uptake of Self or CpG DNA, Enhance IFN-α Production by Human Plasmacytoid Dendritic Cells, and Promote Inflammation. J. Immunol. 2013, 191, 865–874. [Google Scholar] [CrossRef] [Green Version]
  38. Wolf, R.; Mirmohammadsadegh, A.; Walz, M.; Lysa, B.; Tartler, U.; Remus, R.; Hengge, U.; Michel, G.; Ruzicka, T. Molecular cloning and characterization of alternatively spliced mRNA isoforms from psoriatic skin encoding a novel member of the S100 family. FASEB J. 2003, 17, 1969–1971. [Google Scholar] [CrossRef] [Green Version]
  39. Lee, Y.; Jang, S.; Min, J.-K.; Lee, K.; Sohn, K.-C.; Lim, J.-S.; Im, M.; Lee, H.-E.; Seo, Y.-J.; Kim, C.-D.; et al. S100A8 and S100A9 are messengers in the crosstalk between epidermis and dermis modulating a psoriatic milieu in human skin. Biochem. Biophys. Res. Commun. 2012, 423, 647–653. [Google Scholar] [CrossRef]
  40. Hegyi, Z.; Zwicker, S.; Bureik, D.; Peric, M.; Koglin, S.; Batycka-Baran, A.; Prinz, J.C.; Ruzicka, T.; Schauber, J.; Wolf, R. Vitamin D Analog Calcipotriol Suppresses the Th17 Cytokine–Induced Proinflammatory S100 “Alarmins” Psoriasin (S100A7) and Koebnerisin (S100A15) in Psoriasis. J. Investig. Dermatol. 2012, 132, 1416–1424. [Google Scholar] [CrossRef] [Green Version]
  41. Wolf, R.; Mascia, F.; Dharamsi, A.; Howard, O.M.Z.; Cataisson, C.; Bliskovski, V.; Winston, J.; Feigenbaum, L.; Lichti, U.; Ruzicka, T.; et al. Gene from a Psoriasis Susceptibility Locus Primes the Skin for Inflammation. Sci. Transl. Med. 2010, 2, 61ra90. [Google Scholar] [CrossRef] [PubMed]
  42. Jiang, Q.; Yang, G.; Xiao, F.; Xie, J.; Wang, S.; Lu, L.; Cui, D. Role of Th22 Cells in the Pathogenesis of Autoimmune Diseases. Front. Immunol. 2021, 12, 688066. [Google Scholar] [CrossRef] [PubMed]
  43. Rosen, T.; Nolan, E.M. Metal Sequestration and Antimicrobial Activity of Human Calprotectin Are pH-Dependent. Biochemistry 2020, 59, 2468–2478. [Google Scholar] [CrossRef] [PubMed]
  44. Takeda, A.; Higuchi, D.; Takahashi, T.; Ogo, M.; Baciu, P.; Goetinck, P.F.; Hibino, T. Overexpression of Serpin Squamous Cell Carcinoma Antigens in Psoriatic Skin. J. Investig. Dermatol. 2002, 118, 147–154. [Google Scholar] [CrossRef] [Green Version]
  45. Izuhara, K.; Yamaguchi, Y.; Ohta, S.; Nunomura, S.; Nanri, Y.; Azuma, Y.; Nomura, N.; Noguchi, Y.; Aihara, M. Squamous Cell Carcinoma Antigen 2 (SCCA2, SERPINB4): An Emerging Biomarker for Skin Inflammatory Diseases. Int. J. Mol. Sci. 2018, 19, 1102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Iversen, O.-J.; Lysvand, H.; Hagen, L. The autoantigen Pso p27: A post-translational modification of SCCA molecules. Autoimmunity 2010, 44, 229–234. [Google Scholar] [CrossRef]
  47. Lysvand, H.; Hagen, L.; Klubicka, L.; Slupphaug, G.; Iversen, O.-J. Psoriasis pathogenesis—Pso p27 is generated from SCCA1 with chymase. Biochim. Biophys. Acta (BBA) Mol. Basis Dis. 2014, 1842, 734–738. [Google Scholar] [CrossRef] [Green Version]
  48. Dalaker, T.J.M. Expression of the Psoriasis-associated Antigen, Pso p27, is Inhibited by Cyclosporin A. Acta Derm. Venereol. 1999, 79, 281–284. [Google Scholar] [CrossRef]
  49. Åsbakk, K.; Bergh, K.; Iversen, O. The psoriasis-associated antigen, pso p27, participates in the formation of complement activating immune-complexes in psoriatic scale. APMIS 1990, 98, 143–149. [Google Scholar] [CrossRef]
  50. Johnston, A.; Xing, X.; Guzman, A.M.; Riblett, M.; Loyd, C.M.; Ward, N.L.; Wohn, C.; Prens, E.P.; Wang, F.; Maier, L.E.; et al. IL-1F5, -F6, -F8, and -F9: A Novel IL-1 Family Signaling System That Is Active in Psoriasis and Promotes Keratinocyte Antimicrobial Peptide Expression. J. Immunol. 2011, 186, 2613–2622. [Google Scholar] [CrossRef] [Green Version]
  51. Keermann, M.; Koks, S.; Reimann, E.; Abram, K.; Erm, T.; Silm, H.; Kingo, K. Expression of IL-36 family cytokines and IL-37 but not IL-38 is altered in psoriatic skin. J. Dermatol. Sci. 2015, 80, 150–152. [Google Scholar] [CrossRef] [PubMed]
  52. Boutet, M.; Bart, G.; Penhoat, M.; Amiaud, J.; Brulin, B.; Charrier, C.; Morel, F.; Lecron, J.; Rolli-Derkinderen, M.; Bourreille, A.; et al. Distinct expression of interleukin (IL)-36α, β and γ, their antagonist IL-36Ra and IL-38 in psoriasis, rheumatoid arthritis and Crohn’s disease. Clin. Exp. Immunol. 2016, 184, 159–173. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Jiang, Z.; Liu, Y.; Li, C.; Chang, L.; Wang, W.; Wang, Z.; Gao, X.; Ryffel, B.; Wu, Y.; Lai, Y. IL-36γ Induced by the TLR3-SLUG-VDR Axis Promotes Wound Healing via REG3A. J. Investig. Dermatol. 2017, 137, 2620–2629. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Furue, K.; Ito, T.; Tanaka, Y.; Yumine, A.; Hashimoto-Hachiya, A.; Takemura, M.; Murata, M.; Yamamura, K.; Tsuji, G.; Furue, M. Cyto/chemokine profile of in vitro scratched keratinocyte model: Implications of significant upregulation of CCL20, CXCL8 and IL36G in Koebner phenomenon. J. Dermatol. Sci. 2019, 94, 244–251. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Miura, S.; Garcet, S.; Salud-Gnilo, C.; Gonzalez, J.; Li, X.; Murai-Yamamura, M.; Yamamura, K.; Rambhia, D.; Kunjravia, N.; Guttman-Yassky, E.; et al. IL-36 and IL-17A Cooperatively Induce a Psoriasis-Like Gene Expression Response in Human Keratinocytes. J. Investig. Dermatol. 2021, 141, 2086–2090. [Google Scholar] [CrossRef]
  56. Catapano, M.; Vergnano, M.; Romano, M.; Mahil, S.K.; Choon, S.-E.; Burden, A.D.; Young, H.S.; Carr, I.M.; Lachmann, H.J.; Lombardi, G.; et al. IL-36 Promotes Systemic IFN-I Responses in Severe Forms of Psoriasis. J. Investig. Dermatol. 2020, 140, 816–826.e3. [Google Scholar] [CrossRef] [Green Version]
  57. Foster, A.M.; Baliwag, J.; Chen, C.S.; Guzman, A.M.; Stoll, S.W.; Gudjonsson, J.E.; Ward, N.L.; Johnston, A. IL-36 Promotes Myeloid Cell Infiltration, Activation, and Inflammatory Activity in Skin. J. Immunol. 2014, 192, 6053–6061. [Google Scholar] [CrossRef] [Green Version]
  58. Akdis, M.; Palomares, O.; Van De Veen, W.; Van Splunter, M.; Akdis, C.A. TH17 and TH22 cells: A confusion of antimicrobial response with tissue inflammation versus protection. J. Allergy Clin. Immunol. 2012, 129, 1438–1449. [Google Scholar] [CrossRef]
  59. Srivastava, A.; Luo, L.; Lohcharoenkal, W.; Meisgen, F.; Pasquali, L.; Pivarcsi, A.; Sonkoly, E. Cross-talk between IFN-γ and TWEAK through miR-149 amplifies skin inflammation in psoriasis. J. Allergy Clin. Immunol. 2021, 147, 2225–2235. [Google Scholar] [CrossRef]
  60. Fritz, Y.; Klenotic, P.A.; Swindell, W.R.; Yin, Z.Q.; Groft, S.G.; Zhang, L.; Baliwag, J.; Camhi, M.I.; Diaconu, D.; Young, A.B.; et al. Induction of Alternative Proinflammatory Cytokines Accounts for Sustained Psoriasiform Skin Inflammation in IL-17C+IL-6KO Mice. J. Investig. Dermatol. 2017, 137, 696–705. [Google Scholar] [CrossRef] [Green Version]
  61. Carrier, Y.; Ma, H.-L.; Ramon, H.E.; Napierata, L.; Small, C.; O’Toole, M.; Young, D.A.; Fouser, L.A.; Nickerson-Nutter, C.; Collins, M.; et al. Inter-Regulation of Th17 Cytokines and the IL-36 Cytokines In Vitro and In Vivo: Implications in Psoriasis Pathogenesis. J. Investig. Dermatol. 2011, 131, 2428–2437. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Wawrzycki, B.; Pietrzak, A.; Grywalska, E.; Krasowska, D.; Chodorowska, G.; Roliński, J. Interleukin-22 and Its Correlation with Disease Activity in Plaque Psoriasis. Arch. Immunol. Ther. Exp. 2019, 67, 103–108. [Google Scholar] [CrossRef] [Green Version]
  63. Müller, A.; Hennig, A.; Lorscheid, S.; Grondona, P.; Schulze-Osthoff, K.; Hailfinger, S.; Kramer, D. IκBζ is a key transcriptional regulator of IL-36–driven psoriasis-related gene expression in keratinocytes. Proc. Natl. Acad. Sci. USA 2018, 115, 10088–10093. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Hashiguchi, Y.; Yabe, R.; Chung, S.-H.; Murayama, M.A.; Yoshida, K.; Matsuo, K.; Kubo, S.; Saijo, S.; Nakamura, Y.; Matsue, H.; et al. IL-36α from Skin-Resident Cells Plays an Important Role in the Pathogenesis of Imiquimod-Induced Psoriasiform Dermatitis by Forming a Local Autoamplification Loop. J. Immunol. 2018, 201, 167–182. [Google Scholar] [CrossRef]
  65. Johnston, A.; Fritz, Y.; Dawes, S.M.; Diaconu, D.; Al-Attar, P.M.; Guzman, A.M.; Chen, C.S.; Fu, W.; Gudjonsson, J.E.; McCormick, T.S.; et al. Keratinocyte Overexpression of IL-17C Promotes Psoriasiform Skin Inflammation. J. Immunol. 2013, 190, 2252–2262. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Ramirez-Carrozzi, V.; Sambandam, A.; Luis, E.; Lin, Z.; Jeet, S.; Lesch, J.; Hackney, J.; Kim, J.; Zhou, M.; Lai, J.; et al. IL-17C regulates the innate immune function of epithelial cells in an autocrine manner. Nat. Immunol. 2011, 12, 1159–1166. [Google Scholar] [CrossRef]
  67. Lauffer, F.; Jargosch, M.; Baghin, V.; Krause, L.; Kempf, W.; Absmaier-Kijak, M.; Morelli, M.; Madonna, S.; Marsais, F.; Lepescheux, L.; et al. IL-17C amplifies epithelial inflammation in human psoriasis and atopic eczema. J. Eur. Acad. Dermatol. Venereol. 2020, 34, 800–809. [Google Scholar] [CrossRef]
  68. Li, H.; Chen, J.; Huang, A.; Stinson, J.; Heldens, S.; Foster, J.; Dowd, P.; Gurney, A.L.; Wood, W.I. Cloning and characterization of IL-17B and IL-17C, two new members of the IL-17 cytokine family. Proc. Natl. Acad. Sci. USA 2000, 97, 773–778. [Google Scholar] [CrossRef] [Green Version]
  69. Johansen, C.; Mose, M.; Ommen, P.; Bertelsen, T.; Vinter, H.; Hailfinger, S.; Lorscheid, S.; Schulze-Osthoff, K.; Iversen, L. IκBζ is a key driver in the development of psoriasis. Proc. Natl. Acad. Sci. USA 2015, 112, E5825–E5833. [Google Scholar] [CrossRef] [Green Version]
  70. Muromoto, R.; Hirao, T.; Tawa, K.; Hirashima, K.; Kon, S.; Kitai, Y.; Matsuda, T. IL-17A plays a central role in the expression of psoriasis signature genes through the induction of IκB-ζ in keratinocytes. Int. Immunol. 2016, 28, 443–452. [Google Scholar] [CrossRef] [Green Version]
  71. Ovesen, S.; Schulze-Osthoff, K.; Iversen, L.; Johansen, C. IkBζ is a Key Regulator of Tumour Necrosis Factor-a and Interleukin-17A-mediated Induction of Interleukin-36g in Human Keratinocytes. Acta Derm. Venereol. 2021, 101, adv00386. [Google Scholar] [CrossRef]
  72. Johansen, C.; Bertelsen, T.; Ljungberg, C.; Mose, M.; Iversen, L. Characterization of TNF-α– and IL-17A–Mediated Synergistic Induction of DEFB4 Gene Expression in Human Keratinocytes through IκBζ. J. Investig. Dermatol. 2016, 136, 1608–1616. [Google Scholar] [CrossRef] [Green Version]
  73. Barker, J.N.; Jones, M.L.; Mitra, R.S.; Crockett-Torabe, E.; Fantone, J.C.; Kunkel, S.L.; Warren, J.S.; Dixit, V.M.; Nickoloff, B.J. Modulation of keratinocyte-derived interleukin-8 which is chemotactic for neutrophils and T lymphocytes. Am. J. Pathol. 1991, 139, 869–876. [Google Scholar] [PubMed]
  74. Koch, A.E.; Polverini, P.J.; Kunkel, S.L.; Harlow, L.A.; DiPietro, L.A.; Elner, V.M.; Elner, S.G.; Strieter, R.M. Interleukin-8 as a macrophage-derived mediator of angiogenesis. Science 1992, 258, 1798–1801. [Google Scholar] [CrossRef] [PubMed]
  75. Nickoloff, B.J.; Mitra, R.S.; Varani, J.; Dixit, V.M.; Polverini, P.J. Aberrant production of interleukin-8 and thrombospondin-1 by psoriatic keratinocytes mediates angiogenesis. Am. J. Pathol. 1994, 144, 820–828. [Google Scholar] [PubMed]
  76. Tuschil, A.; Lam, C.; Haslberger, A.; Lindley, I. Interleukin-8 Stimulates Calcium Transients and Promotes Epidermal Cell Proliferation. J. Investig. Dermatol. 1992, 99, 294–298. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Duan, H.; Koga, T.; Kohda, F.; Hara, H.; Urabe, K.; Furue, M. Interleukin-8-positive neutrophils in psoriasis. J. Dermatol. Sci. 2001, 26, 119–124. [Google Scholar] [CrossRef]
  78. Guilloteau, K.; Paris, I.; Pedretti, N.; Boniface, K.; Juchaux, F.; Huguier, V.; Guillet, G.; Bernard, F.-X.; Lecron, J.-C.; Morel, F. Skin Inflammation Induced by the Synergistic Action of IL-17A, IL-22, Oncostatin M, IL-1α, and TNF-α Recapitulates Some Features of Psoriasis. J. Immunol. 2010, 184, 5263–5270. [Google Scholar] [CrossRef] [Green Version]
  79. Srivastava, A.; Nikamo, P.; Lohcharoenkal, W.; Li, D.; Meisgen, F.; Landén, N.X.; Ståhle, M.; Pivarcsi, A.; Sonkoly, E. MicroRNA-146a suppresses IL-17–mediated skin inflammation and is genetically associated with psoriasis. J. Allergy Clin. Immunol. 2017, 139, 550–561. [Google Scholar] [CrossRef] [Green Version]
  80. Chen, W.; Nyuydzefe, M.S.; Weiss, J.M.; Zhang, J.; Waksal, S.D.; Zanin-Zhorov, A. ROCK2, but not ROCK1 interacts with phosphorylated STAT3 and co-occupies TH17/TFH gene promoters in TH17-activated human T cells. Sci. Rep. 2018, 8, 1–10. [Google Scholar] [CrossRef]
  81. Lock, F.E.; Hotchin, N.A. Distinct Roles for ROCK1 and ROCK2 in the Regulation of Keratinocyte Differentiation. PLoS ONE 2009, 4, e8190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Bryan, B.A.; Dennstedt, E.; Mitchell, D.C.; Walshe, T.E.; Noma, K.; Loureiro, R.; Saint-Geniez, M.; Campaigniac, J.; Liao, J.K.; Patricia, D.A. RhoA/ROCK signaling is essential for multiple aspects of VEGF-mediated angiogenesis. FASEB J. 2010, 24, 3186–3195. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Zanin-Zhorov, A.; Weiss, J.M.; Nyuydzefe, M.S.; Chen, W.; Scher, J.U.; Mo, R.; Depoil, D.; Rao, N.; Liu, B.; Wei, J.; et al. Selective oral ROCK2 inhibitor down-regulates IL-21 and IL-17 secretion in human T cells via STAT3-dependent mechanism. Proc. Natl. Acad. Sci. USA 2014, 111, 16814–16819. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Zanin-Zhorov, A.; Weiss, J.M.; Trzeciak, A.; Chen, W.; Zhang, J.; Nyuydzefe, M.S.; Arencibia, C.; Polimera, S.; Schueller, O.; Fuentes-Duculan, J.; et al. Cutting Edge: Selective Oral ROCK2 Inhibitor Reduces Clinical Scores in Patients with Psoriasis Vulgaris and Normalizes Skin Pathology via Concurrent Regulation of IL-17 and IL-10. J. Immunol. 2017, 198, 3809–3814. [Google Scholar] [CrossRef] [Green Version]
  85. Chuang, H.-H.; Yang, C.-H.; Tsay, Y.-G.; Hsu, C.-Y.; Tseng, L.-M.; Chang, Z.-F.; Lee, H.-H. ROCKII Ser1366 phosphorylation reflects the activation status. Biochem. J. 2012, 443, 145–151. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Lowery, D.M.; Clauser, K.; Hjerrild, M.; Lim, D.; Alexander, J.; Kishi, K.; Ong, S.-E.; Gammeltoft, S.; A Carr, S.; Yaffe, M.B. Proteomic screen defines the Polo-box domain interactome and identifies Rock2 as a Plk1 substrate. EMBO J. 2007, 26, 2262–2273. [Google Scholar] [CrossRef] [PubMed]
  87. Couzens, A.L.; Saridakis, V.; Scheid, M.P. The hydrophobic motif of ROCK2 requires association with the N-terminal extension for kinase activity. Biochem. J. 2009, 419, 141–148. [Google Scholar] [CrossRef]
  88. Vacharanukrauh, P.; Meephansan, J.; Tangtanatakul, P.; Soonthornchai, W.; Wongpiyabovorn, J.; Serirat, O.; Komine, M. High-Throughput RNA Sequencing Reveals the Effect of NB-UVB Phototherapy on Major Inflammatory Molecules of Lesional Psoriasis. Psoriasis Targets Ther. 2021, 11, 133–149. [Google Scholar] [CrossRef]
  89. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. EdgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [Green Version]
  90. McCarthy, D.J.; Chen, Y.; Smyth, G.K. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012, 40, 4288–4297. [Google Scholar] [CrossRef] [Green Version]
  91. Loraine, A.E.; Blakley, I.C.; Jagadeesan, S.; Harper, J.; Miller, G.; Firon, N. Analysis and Visualization of RNA-Seq Expression Data Using RStudio, Bioconductor, and Integrated Genome Browser. In Plant Functional Genomics: Methods and Protocols; Alonso, J.M., Stepanova, A.N., Eds.; Springer: New York, NY, USA, 2015; pp. 481–501. [Google Scholar]
  92. Zhou, Y.; Zhou, B.; Pache, L.; Chang, M.; Khodabakhshi, A.H.; Tanaseichuk, O.; Benner, C.; Chanda, S.K. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. Commun. 2019, 10, 1523. [Google Scholar] [CrossRef] [PubMed]
  93. Schmittgen, T.D.; Livak, K.J. Analyzing real-time PCR data by the comparative CT method. Nat. Protoc. 2008, 3, 1101–1108. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The differentially expressed genes (DEGs). A total of l202 genes was identified as the DEGs in the peripheral edge of the lesional (PE) skin vs. uninvolved (UN) skin. Of these, 653 (54%) were upregulated and 549 (46%) were downregulated. Among these DEGs, 135 DEGs are overlapped with previously reported psoriasis-associated genes or proteins. Red presents upregulation in the PE skin and green represents downregulation in the PE skin. Proteins of the 86 upregulated DEGs might be potential molecules involved in psoriasis pathogenesis in the PE skin. The DEGs are listed based on their molecule type.
Figure 1. The differentially expressed genes (DEGs). A total of l202 genes was identified as the DEGs in the peripheral edge of the lesional (PE) skin vs. uninvolved (UN) skin. Of these, 653 (54%) were upregulated and 549 (46%) were downregulated. Among these DEGs, 135 DEGs are overlapped with previously reported psoriasis-associated genes or proteins. Red presents upregulation in the PE skin and green represents downregulation in the PE skin. Proteins of the 86 upregulated DEGs might be potential molecules involved in psoriasis pathogenesis in the PE skin. The DEGs are listed based on their molecule type.
Ijms 23 04983 g001
Figure 2. DEG-enriched pathway and function analysis (A) The top 20 upregulated DEG-enriched pathways and functions. (B) The top 20 downregulated DEG-enriched pathways and functions. A −log10 (10) value is the p-value in log base 10 to create the height of each bar chart. The associated pathways and functions are annotated with *. CORUM, comprehensive resource of mammalian protein complexes; DEG, differentially expressed gene; GO, gene ontology; ko, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway; M, canonical pathways; R-HSA, reactome gene sets; WP, WikiPathways.
Figure 2. DEG-enriched pathway and function analysis (A) The top 20 upregulated DEG-enriched pathways and functions. (B) The top 20 downregulated DEG-enriched pathways and functions. A −log10 (10) value is the p-value in log base 10 to create the height of each bar chart. The associated pathways and functions are annotated with *. CORUM, comprehensive resource of mammalian protein complexes; DEG, differentially expressed gene; GO, gene ontology; ko, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway; M, canonical pathways; R-HSA, reactome gene sets; WP, WikiPathways.
Ijms 23 04983 g002
Figure 3. Potential upregulated differentially expressed genes (DEGs). The circles represent the fold change of expressed genes; the larger the circle, the more upregulated the expression of DEGs. A2M, alpha 2 macroglobulin; A2ML1, alpha 2 macroglobulin-like 1; ADAM17, ADAM metallopeptidase domain 17; C10orf99, chromosome 10 open reading frame 99; CXCL8, C-X-C motif chemokine ligand 8; DEFB4A/DEFB4B, defensin beta 4A; DSC2, desmocollin2; DSG3, desmoglein 3; EDNRB, endothelin receptor type B; FLG, filaggrin; GFRA1, GDNF family receptor alpha 1; GJB2, gap junction protein beta 2; IL36A, interleukin 36 alpha; ITGB1, integrin subunit beta 1; KLK7, kallikrein related peptidase 7; KRT16P1, keratin 16 pseudogene 1; LAMA2, laminin subunit alpha 2; LCE3A, late cornified envelope 3A; LCN2, lipocalin 2; mDC, myeloid dendritic cells; NAMPT, nicotinamide phosphoribosyltransferase; pDC, plasmacytoid dendritic cell; PLAT, plasminogen activator, tissue type; PLEC, plectin; PLSCR1, phospholipid scramblase 1; PMN, polymorphonuclear leukocytes; S100A7A, S100 calcium binding protein A7A; SERPINB4, serpin family B member 4; Th 0, T helper 0 cell; TNFRSF1A, TNF receptor superfamily member 1A; VAV3, vav guanine nucleotide exchange factor 3; XDH, xanthine dehydrogenase.
Figure 3. Potential upregulated differentially expressed genes (DEGs). The circles represent the fold change of expressed genes; the larger the circle, the more upregulated the expression of DEGs. A2M, alpha 2 macroglobulin; A2ML1, alpha 2 macroglobulin-like 1; ADAM17, ADAM metallopeptidase domain 17; C10orf99, chromosome 10 open reading frame 99; CXCL8, C-X-C motif chemokine ligand 8; DEFB4A/DEFB4B, defensin beta 4A; DSC2, desmocollin2; DSG3, desmoglein 3; EDNRB, endothelin receptor type B; FLG, filaggrin; GFRA1, GDNF family receptor alpha 1; GJB2, gap junction protein beta 2; IL36A, interleukin 36 alpha; ITGB1, integrin subunit beta 1; KLK7, kallikrein related peptidase 7; KRT16P1, keratin 16 pseudogene 1; LAMA2, laminin subunit alpha 2; LCE3A, late cornified envelope 3A; LCN2, lipocalin 2; mDC, myeloid dendritic cells; NAMPT, nicotinamide phosphoribosyltransferase; pDC, plasmacytoid dendritic cell; PLAT, plasminogen activator, tissue type; PLEC, plectin; PLSCR1, phospholipid scramblase 1; PMN, polymorphonuclear leukocytes; S100A7A, S100 calcium binding protein A7A; SERPINB4, serpin family B member 4; Th 0, T helper 0 cell; TNFRSF1A, TNF receptor superfamily member 1A; VAV3, vav guanine nucleotide exchange factor 3; XDH, xanthine dehydrogenase.
Ijms 23 04983 g003
Figure 4. DEG-enriched canonical pathways. (A,C) The height of the bar represents the level of significance. The predicted activation pathways are shown by orange bars and the predicted inhibition pathways are shown by blue bars. The more intense the color, the stronger the activity prediction. The orange line connecting square points in each bar indicates the ratio of the dataset DEGs that are members of each known pathway divided by the total molecules’ pathway. (B,D) Percentage of DEGs overlapping in each pathway. Percentage of the non-overlapped pathway’s molecules is shown by the white segment of the bar. The total number of molecules in each canonical pathway are shown on the top of the bar. CCR3, C-C Motif Chemokine Receptor 3; CNTF, ciliary neurotrophic factor; CXCL8, C-X-C Motif Chemokine Ligand 8; DEG, differentially expressed gene; eIF2, eukaryotic initiation factor 2; fMLP, formyl-methionyl-leucyl-phenylalanine; GM-CSF, granulocyte-macrophage colony-stimulating factor; HER-2, human epidermal growth factor receptor; IL, interleukin; JAK, Janus Kinase; LXR/RXR, liver X receptor/retinoid X receptor; NER, Nucleotide excision repair.
Figure 4. DEG-enriched canonical pathways. (A,C) The height of the bar represents the level of significance. The predicted activation pathways are shown by orange bars and the predicted inhibition pathways are shown by blue bars. The more intense the color, the stronger the activity prediction. The orange line connecting square points in each bar indicates the ratio of the dataset DEGs that are members of each known pathway divided by the total molecules’ pathway. (B,D) Percentage of DEGs overlapping in each pathway. Percentage of the non-overlapped pathway’s molecules is shown by the white segment of the bar. The total number of molecules in each canonical pathway are shown on the top of the bar. CCR3, C-C Motif Chemokine Receptor 3; CNTF, ciliary neurotrophic factor; CXCL8, C-X-C Motif Chemokine Ligand 8; DEG, differentially expressed gene; eIF2, eukaryotic initiation factor 2; fMLP, formyl-methionyl-leucyl-phenylalanine; GM-CSF, granulocyte-macrophage colony-stimulating factor; HER-2, human epidermal growth factor receptor; IL, interleukin; JAK, Janus Kinase; LXR/RXR, liver X receptor/retinoid X receptor; NER, Nucleotide excision repair.
Ijms 23 04983 g004
Figure 5. Potential regulators and their relationship to the target dataset’s DEGs. In each radiant, the molecule at the center is the predicted potential regulator and the outer circular molecules are the dataset’s DEGs. DEGs, differentially expressed genes; PE skin, peripheral edge of the lesional skin.
Figure 5. Potential regulators and their relationship to the target dataset’s DEGs. In each radiant, the molecule at the center is the predicted potential regulator and the outer circular molecules are the dataset’s DEGs. DEGs, differentially expressed genes; PE skin, peripheral edge of the lesional skin.
Ijms 23 04983 g005
Figure 6. The two most useful regulator effect networks. (A) The highest consistency score network. (B) The second-highest consistency score network. Each network shows how predicted upstream regulators (top tier) control their dataset’s DEGs (middle tier) and the downstream functions of the dataset’s DEGs (bottom tier). DEGs, differentially expressed genes; PE skin, peripheral edge of the lesional skin.
Figure 6. The two most useful regulator effect networks. (A) The highest consistency score network. (B) The second-highest consistency score network. Each network shows how predicted upstream regulators (top tier) control their dataset’s DEGs (middle tier) and the downstream functions of the dataset’s DEGs (bottom tier). DEGs, differentially expressed genes; PE skin, peripheral edge of the lesional skin.
Ijms 23 04983 g006
Figure 7. Quantitative multiplex Real-Time PCR analysis of DEGs and upstream regulators. The mRNA level of each gene was calculated relative to GAPDH, the housekeeping gene. (AJ) The relative RNA expression of selected DEGs. (KM) The relative mRNA expression levels of selected upstream regulators. (N) The fold change of significantly altered genes. All graphs were designated by geometric mean and standard deviation. The significance was determined using *; * p < 0.05, ** p < 0.01 and *** p < 0.001. n ≥ 10 for each gene. CXCL8, C-X-C motif chemokine ligand 8; DEFB103A/B, defensin beta 103A/B; DEGs, differentially expressed genes; IL17C, interleukin 17C; NFKBIZ, nuclear factor kappa B inhibitor zeta; PE skin, peripheral edge of lesional skin; PLSCR1, phospholipid scramblase 1; ROCK2, rho-associated kinase 2; S100A7, S100 calcium binding protein A7; SERPINB4, serpin family B member 4; UN skin, uninvolved skin.
Figure 7. Quantitative multiplex Real-Time PCR analysis of DEGs and upstream regulators. The mRNA level of each gene was calculated relative to GAPDH, the housekeeping gene. (AJ) The relative RNA expression of selected DEGs. (KM) The relative mRNA expression levels of selected upstream regulators. (N) The fold change of significantly altered genes. All graphs were designated by geometric mean and standard deviation. The significance was determined using *; * p < 0.05, ** p < 0.01 and *** p < 0.001. n ≥ 10 for each gene. CXCL8, C-X-C motif chemokine ligand 8; DEFB103A/B, defensin beta 103A/B; DEGs, differentially expressed genes; IL17C, interleukin 17C; NFKBIZ, nuclear factor kappa B inhibitor zeta; PE skin, peripheral edge of lesional skin; PLSCR1, phospholipid scramblase 1; ROCK2, rho-associated kinase 2; S100A7, S100 calcium binding protein A7; SERPINB4, serpin family B member 4; UN skin, uninvolved skin.
Ijms 23 04983 g007
Figure 8. The proposed pathogenesis network at the peripheral edge (PE) of the lesional skin. Injured keratinocytes secrete interleukin (IL)-8, IL-36α, and IL-36γ. IL-36 activates plasmacytoid dendritic cells (pDCs) to secrete interferon-alpha (IFN)-α in a Toll-like receptor (TLR)-9-dependent manner (IL-36/TLR-9 axis). IL-36α and IL-36γ activate myeloid dendritic cells (mDCs) to mature and secrete IL-6, which drives the differentiation of T helper (h) cells (IL-6/Th17 axis). ROCK2 (rho-associated kinase 2) regulates CD4+ T cells to secrete IL-17. Note that inflammation in the PE skin is potentiated by IL-36α, IL-36γ, IL-17C, IL-8, S100A7, S100A8, S100A9, S100A15, SERPINB4, and hBD-2 and that a self-sustaining circuit between the innate and adaptive immune response in the PE skin is mediated through IL-36α and IL-36γ.
Figure 8. The proposed pathogenesis network at the peripheral edge (PE) of the lesional skin. Injured keratinocytes secrete interleukin (IL)-8, IL-36α, and IL-36γ. IL-36 activates plasmacytoid dendritic cells (pDCs) to secrete interferon-alpha (IFN)-α in a Toll-like receptor (TLR)-9-dependent manner (IL-36/TLR-9 axis). IL-36α and IL-36γ activate myeloid dendritic cells (mDCs) to mature and secrete IL-6, which drives the differentiation of T helper (h) cells (IL-6/Th17 axis). ROCK2 (rho-associated kinase 2) regulates CD4+ T cells to secrete IL-17. Note that inflammation in the PE skin is potentiated by IL-36α, IL-36γ, IL-17C, IL-8, S100A7, S100A8, S100A9, S100A15, SERPINB4, and hBD-2 and that a self-sustaining circuit between the innate and adaptive immune response in the PE skin is mediated through IL-36α and IL-36γ.
Ijms 23 04983 g008
Figure 9. Peripheral edge of the lesional skin. The biopsy landmark was identified by the well-demarcated edge of the plaque. The piece of biopsy included equal sides of normal phenotypic skin and lesional skin.
Figure 9. Peripheral edge of the lesional skin. The biopsy landmark was identified by the well-demarcated edge of the plaque. The piece of biopsy included equal sides of normal phenotypic skin and lesional skin.
Ijms 23 04983 g009
Table 1. The top-50 differentially expressed genes of peripheral edge of the lesional skin compared with uninvolved skin.
Table 1. The top-50 differentially expressed genes of peripheral edge of the lesional skin compared with uninvolved skin.
SymbolMolecule Type *Log 2FCp-ValueSymbolMolecule Type *Log 2FCp-Value
Upregulated DEGsDownregulated DEGs
1A2Mtransporter15.9828.600 × 10−6ATP5F1Atransporter−27.3437.380 × 10−6
2H3C2other15.3679.340 × 10−5SEPTIN5enzyme−25.5271.410 × 10−3
3PLECother15.2951.360 × 10−5KHNYNother−24.2543.700 × 10−3
4ITSN2other14.9683.200 × 10−4PLD3enzyme−23.0105.630 × 10−3
5UGGT1enzyme14.9061.080 × 10−4FHL1other−21.9839.840 × 10−3
6STT3Aenzyme14.7453.460 × 10−3HEXAenzyme−20.8268.320 × 10−3
7HERC6enzyme14.7241.120 × 10−3DNM2enzyme−19.5847.500 × 10−3
8DEFB4A/DEFB4Bother14.6901.400 × 10−4TLE5transcription regulator−18.6133.990 × 10−3
9CAND1transcription regulator14.6548.540 × 10−3MYH14enzyme−16.8461.460 × 10−3
10HRNRother14.5971.340 × 10−5GSNother−16.8121.490 × 10−3
11IL36Acytokine14.5721.340 × 10−4SPSB3other−16.5931.100 × 10−3
12SMARCA4transcription regulator14.5172.630 × 10−5LCE1Cother−16.1531.620 × 10−3
13CYFIP1translation regulator14.5002.820 × 10−3LAMTOR4other−16.0589.610 × 10−3
14HIPK3kinase14.3656.200 × 10−3SYT8transporter−16.0149.260 × 10−3
15PLA2G4E-AS1other14.2748.590 × 10−5LTBP4growth factor−15.9852.600 × 10−3
16TRIP12enzyme14.1464.660 × 10−3ITGB4transmembrane receptor−15.7173.270 × 10−3
17CLIC1ion channel14.0914.110 × 10−5DCTN1other−15.5913.700 × 10−3
18CD164other14.0691.280 × 10−3RTN4other−15.5323.930 × 10−3
19BZW1translation regulator14.0256.550 × 10−3CAMSAP3other−15.5173.990 × 10−3
20NRBP1kinase13.9391.280 × 10−4NFIXtranscription regulator−15.2665.070 × 10−3
21PSMD12other13.6978.810 × 10−5CRATenzyme−15.2264.210 × 10−3
22EPB41L3other13.6622.000 × 10−3RAB40Cenzyme−15.1565.860 × 10−3
23SCYL2other13.6521.490 × 10−4CCL27other−15.0213.820 × 10−3
24H2BC10other13.6101.870 × 10−4SIK1/SIK1Bkinase−14.9946.560 × 10−3
25WDR75other13.5541.180 × 10−3CAPGother−14.9702.080 × 10−3
26H4C13other13.5501.080 × 10−3EPN2other−14.8797.780 × 10−3
27EXT2enzyme13.4324.740 × 10−4ARAP1other−14.8717.950 × 10−3
28H3C7other13.4051.480 × 10−3PDCD4other−14.8587.870 × 10−3
29PLEKHG5other13.3192.030 × 10−4XPO6other−14.7359.330 × 10−3
30DMXL1other13.3021.930 × 10−3PODNother−14.7059.560 × 10−3
31GBP6enzyme13.2782.320 × 10−3TUBGCP2peptidase−14.6949.730 × 10−3
32H3C11other13.0994.590 × 10−4CIRBPtranslation regulator−14.6589.860 × 10−3
33H2AC14other13.0924.200 × 10−6TSPAN14other−14.6409.940 × 10−3
34SMC6other12.9172.030 × 10−4ATP2B4transporter−14.5061.540 × 10−3
35H2AC21other12.8477.450 × 10−4CNOT1other−14.4148.590 × 10−3
36H2AC4other12.7935.310 × 10−4WBP2transcription regulator−14.4033.680 × 10−3
37CD59other12.7929.370 × 10−4MAN2C1enzyme−14.3242.770 × 10−3
38TMEM14EPother12.7536.240 × 10−4POLD2enzyme−14.2936.080 × 10−3
39GLB1enzyme12.7171.270 × 10−3MPSTenzyme−14.2375.410 × 10−3
40CTNNBL1other12.7032.090 × 10−3ACSS2enzyme−14.0294.050 × 10−3
41GFPT1enzyme12.6373.910 × 10−3AP2A1transporter−13.9444.740 × 10−3
42H3C8other12.5991.300 × 10−4SCRIBother−13.9006.080 × 10−3
43ZDHHC6enzyme12.5654.940 × 10−3TAP2transporter−13.8088.590 × 10−3
44IKBKBkinase12.4372.280 × 10−3PHYHIPother−13.7836.750 × 10−3
45GANABenzyme12.4122.170 × 10−4PLCH2enzyme−13.786.560 × 10−3
46TMEM184Bother12.4101.410 × 10−3TMEM63Bion channel−13.6176.260 × 10−3
47SEPTIN11other12.3779.900 × 10−4TAX1BP1other−13.6003.020 × 10−3
48RMDN3other12.3223.100 × 10−3PKMkinase−13.5945.310 × 10−3
49LPPother12.2925.600 × 10−3EPHX3enzyme−13.4728.680 × 10−3
50H2BC13other12.2922.600 × 10−3PTBP1enzyme−13.1678.660 × 10−3
Based on QIAGEN Knowledge Base as of October 2021, the blue molecules have been previously mentioned in psoriasis, either their mRNAs or their proteins. * Molecule types are classified into complex, cytokine, enzyme, fusion gene/product, G-protein coupled receptor, group, growth factor, ion channel, kinase, ligand-dependent nuclear receptor, mature microRNA, peptidase, phosphatase, transcription regulator, translation regulator, transmembrane receptor and transporter. None of these were classified as other. FC, fold change.
Table 2. The potential differentially expressed genes (DEGs).
Table 2. The potential differentially expressed genes (DEGs).
Gene SymbolLog2FCp-ValueGene SymbolLog2FCp-Value
UpregulationDownregulation
Anti-microbial peptides (AMPs)S100 proteins
DEFB4A/DEFB4B14.6901.400 × 10−4S100A16−5.3349.920 × 10−3
DEFB103A/DEFB103B7.6242.290 × 10−4S100A10−4.1002.780 × 10−5
S100A7A (S100A15) 9.5591.180 × 10−5S100A4−3.8191.450 × 10−4
S100A84.9442.590 × 10−4Cytokines
S100A94.7223.440 × 10−5FBRS−3.2042.140 × 10−3
S100A73.0023.930 × 10−5SLURP1−1.9753.420 × 10−3
LCN24.1828.240 × 10−4TIMP1−1.8574.290 × 10−3
Cytokines MIF−1.8424.860 × 10−3
IL36A14.5721.340 × 10−4CXCL14−1.7521.650 × 10−3
CXCL810.6816.580 × 10−4Transmembrane receptor
VAV39.1781.290 × 10−4ITGB4−15.7173.270 × 10−3
NAMPT3.9198.090 × 10−4Cell adhesion molecules
IL36G3.8972.090 × 10−3GJB3−5.4549.220 × 10−3
C10orf992.5672.490 × 10−3GJB5−2.3179.110 × 10−3
Transmembrane receptors LAMA5−2.1257.100 × 10−3
GFRA19.0235.660 × 10−3Keratin and late cornified envelope
TNFRSF1A3.3171.180 × 10−3KRT2−5.4849.530 × 10−4
ITGB12.0308.210 × 10−3KRT1−3.7003.250 × 10−4
Cell adhesion molecules KRT10−3.5211.510 × 10−4
DSC26.5951.260 × 10−4KRT77−3.4164.850 × 10−4
DSC12.4786.390 × 10−3LCE1C−16.1531.620 × 10−3
DSG33.5991.590 × 10−3LCE2C/LCE2D−5.3669.450 × 10−4
DSG12.4567.280 × 10−3LCE2B−4.6896.280 × 10−5
GJB25.4683.430 × 10−6LCE1B−3.8718.320 × 10−5
GJB65.1354.940 × 10−3LCE6A−3.5635.910 × 10−4
PLEC15.2951.360 × 10−5LCE2A−3.4671.230 × 10−4
LAMA21.6268.630 × 10−3LCE1D−3.3719.840 × 10−4
Keratin and late cornified envelopeLCE1A−3.1321.310 × 10−3
KRT16P18.1208.130 × 10−3LCE1F−2.9561.890 × 10−3
KRT6A5.7403.170 × 10−4Peptidase
KRT6C5.3101.170 × 10−4KLK11−5.9922.640 × 10−5
KRT164.6758.870 × 10−5Other interesting DEGs
KRT6B3.3842.000 × 10−3APOE−3.3006.710 × 10−3
LCE3A7.3832.430 × 10−5CALML5−2.1536.690 × 10−4
LCE3C5.4511.000 × 10−3CAV1−2.1215.430 × 10−3
LCE3E2.6522.210 × 10−3CCL27−15.0213.820 × 10−3
Peptidase CCND1−1.8274.370 × 10−3
PLAT4.3885.750 × 10−4CD81−3.6031.190 × 10−4
KLK74.2247.300 × 10−4COL16A1−2.6587.960 × 10−4
KLK63.9272.030 × 10−3COL1A1−2.6361.870 × 10−4
ADAM172.1463.400 × 10−3COL6A1−3.5478.270 × 10−5
Protease inhibitor COL6A2−3.3573.740 × 10−3
A2M15.9828.600 × 10−6COL7A1−1.7661.970 × 10−3
A2ML11.7037.220 × 10−3FBL−3.1224.290 × 10−3
SPINK51.5744.340 × 10−3GSN−16.8121.490 × 10−3
Other interesting DEGs GPX4−3.1374.430 × 10−3
CD5912.7929.370 × 10−3SERPINF1−2.1153.820 × 10−3
IKBKB12.4372.280 × 10−3IL20RB−1.7882.680 × 10−3
EDNRB11.8155.820 × 10−3IGFBP6−5.3991.260 × 10−3
TNIP310.7946.580 × 10−3IGFBP4−3.2481.140 × 10−4
SERPINB410.6917.430 × 10−3IGFBP7−1.5746.940 × 10−3
SERPINB37.6303.290 × 10−5ITGB5−2.0777.180 × 10−3
SERPINB93.9606.270 × 10−3LGALS7/LGALS7B−4.8808.880 × 10−5
XDH5.5656.140 × 10−3LGALS1−3.6811.030 × 10−3
PI35.1395.430 × 10−5LORICRIN−9.2448.540 × 10−4
STAT13.3843.940 × 10−3MMP28−11.5747.540 × 10−3
STAT31.6178.030 × 10−3NUPR1−2.1833.910 × 10−3
DDX212.5202.580 × 10−3PECAM1−2.5172.980 × 10−3
PLSCR12.3523.520 × 10−3PSORS1C2−3.3738.360 × 10−3
RICTOR2.3388.330 × 10−3TIMP2−3.0923.820 × 10−4
FLG2.1356.550 × 10−4TNXB−2.2382.460 × 10−3
KDM5A1.9052.380 × 10−3TYK2−1.9155.780 × 10−3
HLA-A1.8251.170 × 10−3VEGFB−3.3826.940 × 10−3
IGFBP31.6001.200 × 10−3WNT7B−3.5744.310 × 10−3
These DEGs were selected due to their biofunction and psoriasis association. Based on QIAGEN Knowledge Base as of October 2021, the bold molecules have not been previously mentioned in psoriasis. FC, fold change.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Boonpethkaew, S.; Meephansan, J.; Jumlongpim, O.; Tangtanatakul, P.; Soonthornchai, W.; Wongpiyabovorn, J.; Vipanurat, R.; Komine, M. Transcriptomic Profiling of Peripheral Edge of Lesions to Elucidate the Pathogenesis of Psoriasis Vulgaris. Int. J. Mol. Sci. 2022, 23, 4983. https://doi.org/10.3390/ijms23094983

AMA Style

Boonpethkaew S, Meephansan J, Jumlongpim O, Tangtanatakul P, Soonthornchai W, Wongpiyabovorn J, Vipanurat R, Komine M. Transcriptomic Profiling of Peripheral Edge of Lesions to Elucidate the Pathogenesis of Psoriasis Vulgaris. International Journal of Molecular Sciences. 2022; 23(9):4983. https://doi.org/10.3390/ijms23094983

Chicago/Turabian Style

Boonpethkaew, Suphagan, Jitlada Meephansan, Onjira Jumlongpim, Pattarin Tangtanatakul, Wipasiri Soonthornchai, Jongkonnee Wongpiyabovorn, Ratchanee Vipanurat, and Mayumi Komine. 2022. "Transcriptomic Profiling of Peripheral Edge of Lesions to Elucidate the Pathogenesis of Psoriasis Vulgaris" International Journal of Molecular Sciences 23, no. 9: 4983. https://doi.org/10.3390/ijms23094983

APA Style

Boonpethkaew, S., Meephansan, J., Jumlongpim, O., Tangtanatakul, P., Soonthornchai, W., Wongpiyabovorn, J., Vipanurat, R., & Komine, M. (2022). Transcriptomic Profiling of Peripheral Edge of Lesions to Elucidate the Pathogenesis of Psoriasis Vulgaris. International Journal of Molecular Sciences, 23(9), 4983. https://doi.org/10.3390/ijms23094983

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop