Corrigendum: IL36G is associated with cutaneous antiviral competence in psoriasis
- 1Department of Dermatology and Venereology, First Affiliated Hospital of Kunming Medical University, Kunming, China
- 2Hubei Provincial Key Laboratory of Occurrence and Intervention of Kidney Diseases, Medical College, Hubei Polytechnic University, Huangshi, China
- 3Department of Dermatology, Huangshi Central Hospital, Affiliated Hospital of Hubei Polytechnic University, Edong Healthcare Group, Huangshi, China
- 4NHC Key Laboratory of Drug Addiction Medicine, First Affiliated Hospital of Kunming Medical University, Kunming Medical University, Kunming, China
- 5Scientific Research Laboratory Center, First Affiliated Hospital of Kunming Medical University, Kunming, China
Background: Psoriasis is a common inflammatory skin disease that has a great impact on patients’ physical and mental health. However, the causes and underlying molecular mechanisms of psoriasis are still largely unknown.
Methods: The expression profiles of genes from psoriatic lesion samples and skin samples from healthy controls were integrated via the sva software package, and differentially expressed genes (DEGs) between psoriasis and healthy skin were screened by the limma package. Furthermore, GO and KEGG pathway enrichments for the DEGs were performed using the Clusterprofiler package. Protein–protein interaction (PPI) networks for the DEGs were then constructed to identify hub genes. scGESA analysis was performed on a single-cell RNA sequencing dataset via irGSEA. In order to find the cytokines correlated with the hub genes expression, single cell weighted gene co-expression network analyses (scWGCNA) were utilized to build a gene co-expression network. Furthermore, the featured genes of psoriasis found in suprabasal keratinocytes were intersected with hub genes. We then analyzed the expression of the intersection genes and cytokines in the integrated dataset. After that, we used other datasets to reveal the changes in the intersection genes’ expression levels during biological therapy. The relationship between intersection genes and PASI scores was also explored.
Results: We identified 148 DEGs between psoriatic and healthy samples. GO and KEGG pathway enrichment analysis suggested that DEGs are mainly involved in the defense response to other organisms. The PPI network showed that 11 antiviral proteins (AVPs) were hub genes. scGSEA analysis in the single-cell transcriptome dataset showed that those hub genes are highly expressed in keratinocytes, especially in suprabasal keratinocytes. ISG15, MX1, IFI44L, and IFI27 were the characteristic genes of psoriasis in suprabasal keratinocytes. scWGCNA showed that three cytokines—IL36G, MIF, and IL17RA—were co-expressed in the turquoise module. Only interleukin-36 gamma (IL36G) was positively correlated with AVPs in the integrated dataset. IL36G and AVPs were found co-expressed in a substantial number of suprabasal keratinocytes. Furthermore, we found that the expression levels of IL36G and the 4 AVPs showed positive correlation with PASI score in patients with psoriasis, and that these levels decreased significantly during treatment with biological therapies, but not with methotrexate.
Conclusion: IL36G and antiviral proteins may be closely related with the pathogenesis of psoriasis, and they may represent new candidate molecular markers for the occurrence and severity of psoriasis.
Introduction
Psoriasis is a common immune-mediated inflammatory dermatosis, prevalent worldwide, that can have serious effects on patients’ physical and mental health (1). The pathogenesis of psoriasis depends on genetic susceptibility, infection, immune abnormalities, and psychoneurotic factors, resulting in chronicity and leaving patients prone to the recurrence of erythematous plaques with adherent silvery scales (2). The main pathological characteristics of psoriatic lesions are hyperproliferation and aberrant differentiation of epidermal keratinocytes. Even though much is known about this condition, the multifactorial nature of its pathogenesis has not been studied thoroughly. In particular, the extent of the role infection plays in the pathogenesis of psoriasis remains controversial.
A previous study found that persistent subclinical streptococcal and staphylococcal infections may be responsible for not only recurrent acute guttate psoriasis, but also new episodes of chronic plaque psoriasis (3). Recent studies have demonstrated that the pathophysiology of psoriasis is closely associated with antimicrobial peptides (AMPs), including cathelicidin (LL-37), human β-defensins, S100 proteins, and lipocalin 2 (4–6). AMPs are an essential part of the innate immune system’s defense against pathogen-based infections (7). Furthermore, it has been confirmed that infection is involved in the pathological mechanism of psoriasis. Levels of keratinocyte-derived antiviral proteins (AVPs) of 2–5-oligoadenylate synthase 2 (OAS2) were found to be significantly elevated in both the epidermis and serum of psoriasis patients (8). Gao et al. found that AVPs, such as ISG15, MX1, OAS2, OASL, and OAS3, were overexpressed in cultured human keratinocyte cells (HaCaT) stimulated by tumor necrosis factor-alpha (TNF-α) in a psoriasis cell model (9). Moreover, ISG15, RSAD2, IRF7, MX2, and TRIM22 have been found to be overexpressed in psoriatic skin. However, another study found no increased expression of AVPs in skin affected by atopic dermatitis, a chronic inflammatory disease similar to psoriasis (10). It is therefore of interest to determine the link between skin antiviral phenotype and the pathogenesis of psoriasis. Therefore, it is imperative to explore the roles of AVPs in the occurrence of psoriasis.
Herein, we merged five mRNA microarray datasets from the gene expression omnibus (GEO) database to screen for differentially expressed genes (DEGs) between psoriasis and healthy samples. Subsequently, enrichment analysis and protein–protein interaction (PPI) network analysis were performed to screen for the hub gene. We then explored the underlying mechanism of hub genes in a single-cell transcriptome. The workflow for this study is shown in Figure S1.
Material and methods
Data sources and preprocessing
We obtained the psoriasis transcriptome dataset from GEO, which included the GSE13355 dataset (58 psoriatic lesion samples and 58 non-lesion samples), the GSE30999 dataset (85 psoriatic lesion samples and 85 non-lesion samples), the GSE34248 dataset (14 psoriatic lesion samples and 14 non-lesion samples), the GSE41662 dataset (24 psoriatic lesion samples and 24 non-lesion samples), the GSE14905 dataset (33 psoriatic lesion samples and 28 non-lesion samples), and the GSE162183 dataset (including the lesions of 3 patients and similar regions from 3 healthy donors). Besides these datasets, the GSE85034 (30 psoriatic lesion tissues before and after adalimumab or methotrexate (MTX) treatment), GSE51440 (39 psoriatic lesion tissues before and after guselkumab treatment), GSE117468 (41 psoriatic lesion tissues before and after brodalumab treatment and 44 psoriatic lesion tissues before and after ustekinumab treatment), GSE41664 (36 psoriatic lesion tissues before and after etanercept treatment), and GSE69967 (26 psoriatic lesion tissues before and after tofacitinib treatment) datasets are also included in our study (Table 1).
All patients whose samples were included in this study were diagnosed with psoriasis vulgaris. The list of cytokines was obtained from https://www.immport.org/resources/cytokineRegistry. During data integration, when a given gene expression was detected across all five datasets, the genes were preserved for further analysis. The expression profiles were integrated via the Combat function of the sva package (11) to remove the batch effects.
Analysis of differentially expressed genes
The DEGs between psoriasis and normal tissues were identified by the limma R software package. Moreover, |log2 fold change (FC) |≥ 2, with an adjusted p-value of < 0.05, was set as the threshold for differentially expressed genes selection.
GO and KEGG pathway enrichment analyses
The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed using the clusterProfiler software package (v 4.4.1) to explore the biological function of the DEGs. GO terms of the KEGG pathway with a false discovery rate (FDR) < 0.05 were considered to be significantly enriched.
Construction of the protein–protein interaction network and hub gene screening
The DEGs were imported into the STRING protein database (STRING11.5; https://string-db.org/) with a combined score >0.9 was adjusted to obtain the protein–protein interaction network. The plug-in CytoHubba (http://apps.cytoscape.org/cytohubba) was used to calculate the topological characteristic parameters of the nodes in the PPI network to screen out the hub genes. The degree method was used to explore the hub genes in the PPI network, and the screening condition was degrees ≥ 10.
Data processing for single-cell RNA sequencing
The GSE162183 expression matrices were obtained from the GEO database. Quality control and dimensionality reduction were performed by the Seurat software package (v4.1.1). For the initial QC step, we created Seurat objects for the normal and psoriasis groups and filtered out the cells that expressed < 200 genes. Genes expressed in fewer than 3 cells were also excluded. I gene expression profiles of the remaining cells were then normalized, and 2,000 highly variable genes from each sample were identified by the vst method. All genes were scaled, and the principal component analysis was conducted. The cells were clustered by unsupervised clustering (resolution = 0.5) and visualized by umap using the top 20 principal components. Cell-type annotation was performed using the singleR software package (v 1.8.1) and refined with manual annotation. Integration of all single-cell rank-based gene set enrichment analyses with the AVPs was performed by irGSEA (v 1.1.2), and the calculation of the enrichment scoring method was set as UCell.
Single-nucleus consensus weighted gene coexpression network analysis
Single-nucleus consensus weighted gene coexpression network analysis (scWGCNA) is a systematic biological method used to construct a scale-free network based on single-cell gene expression profiles. In the present study, scWGCNA was performed on the expressions and clinical prototypes found in the GSE162183 dataset. First, we constructed metacells with k = 200, then a signed similarity matrix was created. The signed similarity matrix was then raised to power β = 16, varying between cell types, to emphasize strong correlations and reduce the emphasis of weak correlations on an exponential scale. The resulting adjacency matrix was then transformed into a topological overlap matrix. Modules were defined using specific module-cutting parameters, including a minimum module size of 25 genes, a deepSplit score of 2, and a merge Cut Height of 0.6. Module eigengenes (MEs) were identified using the package’s ME function, and signed module membership (MM) was utilized to determine the correlation of MEs with individual genes. Gene significance (GS) was defined as the correlation between the gene expression level and the psoriasis. Modules with the strongest positive correlations with psoriasis were selected.
Statistical analysis
All statistical analyses were performed on R software (v4.1.2). Differences between groups were compared by the Mann-Whitney-Wilcoxon test. The correlation between immune infiltrating cells and core genes was measured by Pear’on’s chi-square test. An adjusted p-value of < 0.05 was considered statistically significant.
Results
Differential gene expression analysis
The GSE13355, GSE30999, GSE34248, GSE41662, and GSE14905 datasets were merged into an integrated dataset, which contained 214 psoriatic lesion samples (LS) and 209 non-lesion samples (NL). Before the adjustment, the clustering of samples was largely driven by batch effect, which is due to merging the data from different datasets (Figure 1A). After applying ComBat, the batch effect was mitigated (i.e., samples from different datasets were mixed together) (Figure 1B). The gene differential expression in the integrated dataset was analyzed via the limma package. Based on standard threshold values for |log2 fold change (FC) |≥ 2, with an adjusted p-value of < 0.05, 148 DEGs were identified, 125 of which were up-regulated and 23 of which were down-regulated (Figures 1C, D) (Table S1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/971071/fimmu-13-971071-HTML-r1/image_m/fimmu-13-971071-g001.jpg)
Figure 1 Analysis of the differentially expressed genes within the integrated dataset. (A) Principal component analysis (PCA) before batch effect adjustment. Samples from the different datasets cluster together. (B) PCA after batch effect adjustment. Samples from different datasets overlap. Red points: samples from GSE1335. Olive green points: samples from GSE14905. Spring green points: samples from GSE30999. Blue points: samples from GSE34248. Pink points: samples from GSE41662. (C) The DEG Volcano map shows genes in psoriatic lesion samples (LS) and non-lesion samples (NL). Up-regulated genes are in red; down-regulated genes are in green. (D) DEG expression heat map of skin tissue, with high expression in red and low expression in green.
Pathway enrichment analysis and protein–protein interaction analysis
To further investigate the biological functions of the DEGs, we performed a pathway enrichment analysis. The result of the GO enrichment analysis showed that DEGs were mainly related to defense responses to other organisms, to the antimicrobial humoral response, and to the response to virus (Figure 2A; Table S2). The KEGG enrichment analysis showed that DEGs were mainly associated with viral proteins interacting with cytokines and cytokine receptors, cytokine–cytokine receptor interactions, and the interleukin (IL)-17 signaling pathway (Figure 2B) (Table S3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/971071/fimmu-13-971071-HTML-r1/image_m/fimmu-13-971071-g002.jpg)
Figure 2 Pathway enrichment analyses and PPI network of DEGs. (A) Bar chart of differentially expressed genes’ GO enrichment analysis. (B) Bubble chart of differentially expressed genes’ KEGG enrichment analysis. The top 10 pathway enrichments are shown. (C) The PPI network of the DEGs. (D) Hub genes’ expression between psoriatic lesion samples (LS) and non-lesion samples (NL). p values were calculated using a Wilcoxon signed-rank test. Asterisks corresponding to p values<0.05 (****< 0.0001).
To screen for hub genes, we first constructed a protein–protein interaction (PPI) network with 32 nodes and 1347 interaction edges (Figure 2C). In the PPI network, a total of 11 genes (RSAD2, IFIT1, MX1, ISG15, IFI27, OAS3, OAS2, IFI6, IRF7, OAS1, and IFI44L) were identified as hub genes with degrees ≥ 10 (Table S4). All of those genes were AVPs, and were significantly upregulated in psoriatic lesions (Figure 2D).
GSEA of scRNA-seq
To investigate the expression of AVPs in psoriatic lesions at single-cell resolution, we collected tissue from 3 psoriasis patients and tissue of a similar region from 3 healthy donors and performed single-cell RNA sequencing (scRNA-seq). After standard data processing and quality control procedures, we obtained transcriptomic profiles for 24,234 cells. Principal cell clusters were classified using an unsupervised graph-based clustering strategy. Cells with similar profiles were annotated using the singleR package.
The cellular composition is shown in Figure 3A. Single-cell rank-based gene set enrichment analysis (GSEA) showed that the enrichment score was higher in keratinocytes (including suprabasal keratinocyte and basal keratinocyte) (Figure 3B). The AVP expression levels in cell populations of the psoriatic lesions and normal skin are shown in Figure S2. ISG15, MX1, IFI44L, and IFI27 were the characteristic psoriasis genes found in suprabasal keratinocytes (Table S5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/971071/fimmu-13-971071-HTML-r1/image_m/fimmu-13-971071-g003.jpg)
Figure 3 Overview of the single-cell landscape for psoriatic lesions. (A) A UMAP view and clustering analysis of combined single-cell transcriptome data from psoriasis and normal skin (n = 24234). Clusters are distinguished by different colors, with the general identity of each cell cluster shown on the right. (B) Density scatterplot of AVP gene set. UMAP, uniform manifold approximation and projection; AVP, antiviral protein.
Use of scWGCNA to screen the cytokines correlated with AVP expression
The scWGCNA package was used to merge 200 cells into a metacell, which was then used to screen module genes associated with psoriasis in the dataset (GSE162183). A dendrogram of samples with clinical traits were clustered using the average linkage method and Pearson’s correlation method (Figure 4A). We determined the soft threshold of 16 by calculating the scale-free model fit and mean connectivity (Figures 4B, C). Different module genes in the dynamic tree cut were re-clustered through a topological similarity strategy. Eventually, seven modules developed after merging, and the turquoise module was found to be strongly correlated with psoriasis (Figure 4D). We then analyzed the correlation of each module with two clinical traits (normal and psoriasis). We found that AVPs were present in the turquoise module (cor = 0.39, p = 2.5e−165) (Figure 4E) and therefore extracted cytokines from that module (Figure 4F; Table S6). We obtained three cytokines (MIF, IL17RA, and IL36G) and analyzed the relationship between those cytokines and the AVPs. Only interleukin-36 gamma (IL36G) was found to positively correlate with AVPs (Table S7).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/971071/fimmu-13-971071-HTML-r1/image_m/fimmu-13-971071-g004.jpg)
Figure 4 Construction of co-expression modules by the scWGCNA. (A) Clustering dendrogram of samples with trait heatmap. (B) The relationship between the scale-free fit index and various soft-thresholding powers. (C) The relationship between the mean connectivity and various soft-thresholding powers. (D) The cluster dendrogram of genes in GSE162183. Branches of the cluster dendrogram of the most connected genes gave rise to 7 gene co-expression modules. (E) Heatmap of the correlation between module eigengenes and clinical phenotype (normal and psoriasis). (F) Scatter plot describing the relationship between MM and GS in the turquoise module.
High levels of AVPs and IL36G in suprabasal keratinocytes
Of all the IL36G-expressing cells in the scRNA-seq data, IL36G-expressing cells were mainly found in suprabasal keratinocytes (Figure 5A). IL36G was upregulated in psoriatic skin compared with normal skin (Figure 5B). The IL36G-expressing suprabasal keratinocytes from normal skin did not express any AVPs (Figure 5C). The IL36G that was highly expressed in the suprabasal keratinocytes from the psoriatic lesions was also highly co-expressed with AVPs (Figure 5D). Suprabasal keratinocytes from psoriatic lesions expressed more AVPs and IL36G than suprabasal keratinocytes from the control skin (Figures 5E, F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/971071/fimmu-13-971071-HTML-r1/image_m/fimmu-13-971071-g005.jpg)
Figure 5 Expression of IL36G and AVPs in suprabasal keratinocytes of psoriatic lesions. (A) IL36G expression in total scRNA-seq data is visualized in low-dimensional space. (B) Violin plot of each cell cluster. (C) Co-expression of IL36G and AVPs in normal skin keratinocytes. (D) Co-expression of IL36G and AVPs in suprabasal keratinocytes of psoriatic lesions. (E, F) Cells with expression of ISG15, MX1, and IL36G (E) and IFI44L, IFI27, and IL36G (F) within suprabasal keratinocytes in normal and psoriatic skin are quantified by number of cells.
Expression of AVPs and IL36G in patients with psoriasis before and after therapy
The expression levels of the AVPs and IL36G in psoriatic lesions were significantly decreased during treatment, especially at 12 weeks or 16 weeks of treatment. The downward trend was even more pronounced with biologic therapies (adalimumab, ustekinumab and brodalumab, tofacitinib, etanercept) than with MTX (Figures 6A–F). We therefore conclude that biological therapy is superior to other drugs in reducing the level of AVPs in psoriasis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/971071/fimmu-13-971071-HTML-r1/image_m/fimmu-13-971071-g006.jpg)
Figure 6 Expression of AVPs and IL36G during treatment. (A) Expression of AVPs and IL36G during ustekinumab treatment. (B) Expression of AVPs and IL36G during brodalumab treatment. (C) Expression of AVPs and IL36G during tofacitinib treatment. (D) Expression of AVPs and IL36G during etanercept treatment. (E) Expression of AVPs and IL36G during adalimumab treatment. (F) Expression of AVPs and IL36G during MTX treatment. MTX, methotrexate. p values were calculated using a Wilcoxon signed-rank test. Asterisks corresponding to p values<0.05 (*< 0.05, **< 0.01, ***< 0.001, ****< 0.0001) and those with p values> 0.05 as not significant "ns" (stat_compare_means function).
The correlation between AVPs, IL36G expression levels, and psoriasis area severity index scores during therapy
In order to explore the association between the AVPs and IL36G expression levels in the context of psoriasis severity, we further investigated the relationship between levels of AVPs and IL36G and Psoriasis Area Severity Index (PASI) scores during psoriasis treatment. Levels of AVPs and IL36G were positively correlated with PASI scores during treatment with adalimumab, ustekinumab, brodalumab, MTX, and tofacitinib. These results therefore suggest that AVPs and IL36G play a role in the severity of psoriasis (Figures 7A–E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/971071/fimmu-13-971071-HTML-r1/image_m/fimmu-13-971071-g007.jpg)
Figure 7 (A–E) The relationship between the expression of three biomarkers and PASI scores during psoriasis treatment. (A) Ustekinumab treatment. (B) Brodalumab treatment. (C) Tofacitinib treatment. (D) Adalimumab treatment. (E) Methotrexate treatment.
Discussion
It is well accepted that infection with many different microorganisms or viruses, including β-haemolytic streptococci, Staphylococcus aureus, Po`rphyromonas gingivalis, Candida albicans, Chlamydia psittaci, HIV, and hepatitis C virus, is a common trigger or exacerbation factor for psoriasis (12). In particular, acute guttate psoriasis is generally thought to be caused by bacterial infection, such as streptococcal infection (13). Beyond the classical role of T helper cell 1 (Th1)/Th2 homeostasis and the IL17/IL23 axis, recently described AMPs have been found to play a role in the immunopathogenesis of psoriasis, further confirming the relationship between infection and psoriasis (4). Numerous high-throughput sequencing studies have suggested that the skin microbiome could play a role in the pathogenesis and therapeutic effect monitoring of psoriasis (14–16) and that AMPs could modify host microenvironments and regulate the colonization of microorganisms (17). It has been suggested that increased expression levels of AMPs could be involved in the pathogenesis of psoriasis (18). AMPs and AVPs are two important functional protein clusters mediating innate immune in psoriatic epidermis. However, the role of AVPs in the keratinocytes of psoriatic lesions remains unelucidated.
In this study, a series of bioinformatics analyses was performed on microarray and single-cell gene expression datasets. We found that the DEGs were mainly associated with a defense response to other organisms, with the type I interferon (IFN) signaling pathway, with the interactions between viral proteins and cytokines, and with the cytokine receptor and cytokine–cytokine receptor interaction pathway. The hub genes were AVPs (RSAD2, IFIT1, MX1, ISG15, IFI27, OAS3, OAS2, IFI6, IRF7, OAS1, and IFI44), which were highly expressed within the keratinocytes in psoriatic lesions. Moreover, IL36G may interplay with ISG15, MX1, IFI44L, and IFI27 in the execution of antiviral function within suprabasal keratinocytes.
Recent findings have demonstrated that the genes of AVPs were overexpressed in psoriatic lesions (8, 10, 19). Zhou et al. found that the keratinocyte-derived OAS2 expression level was positively correlated with PASI scores and decreased after therapy (8). Moreover, in TNF-α–induced HaCaT cells, ISG15 and MX1 were significantly up-regulated (9). Within this series of AVP genes, Mx proteins, such as MX1, are known for inhibiting negative-strand RNA viruses (4, 20).
Interferon-stimulated genes (ISGs) are cellular products that mediate the type I interferon response against a wide range of invading viruses. Only a few ISGs, such as ISG15, IRF7, IFIT1, IFI27, IFI6, and IFI44L, have antimicrobial activity (21). ISG15 (interferon-stimulated gene product 15) participates in many antiviral signaling pathways to directly promote viral clearance (22). IFI27 (interferon α–inducible protein 27) is involved in the proliferation of skin keratinocytes, both in imiquimod-induced psoriasis-like skin and in HaCaT cells (23). Although extensive research confirms that the ISG family are involved in various bacterial and viral infections, the novel functions of AVPs in psoriatic keratinocytes need to be thoroughly investigated.
Although T cells are currently considered to be the main driver of psoriasis, keratinocytes play an important role in regulating inflammation and relapse (24). Epidermal keratinocytes secrete antimicrobial peptides, LL-37, and cytokines like TNF-α, IL6, and IL36 to activate Th1/Th17 cells, which interact with each other to regulate the proliferation and differentiation of keratinocytes (25). IL36G is a pro-inflammatory cytokine whose role in the pathogenesis of psoriasis has been extensively described (26–29). More noteworthy is that our findings showed that AVP genes were abundantly co-expressed with IL36G on psoriatic keratinocytes. Previous research has suggested that IL36G may be an alarmin that signals viral infection (30). One study suggested that IL36G itself is a novel antiviral protein involved in the defense against influenza virus (28). Infectious agents may induce the primary skin barrier layers to produce large amounts of AVPs as a response. Another study found that Th17 cells derived IL29, but not type I IFN-mediated psoriatic keratinocyte is responsible for the expression of the AVPs (31). Whether upregulation of IL36G and AVPs expression in keratinocytes trigger psoriasis is currently unknown.
In addition, the expression levels of AVPs and IL36G were significantly positively correlated with PASI scores, and these levels decreased dramatically during biological treatment. These observations indicate that IL36G and AVPs may be related to the severity of psoriasis and could be used as an indicator for therapeutic efficacy. Research conducted by D’Erme et al. (26) also found that IL36G was closely associated with the disease activity of psoriasis. However, the exact production mechanism of IL36G and AVPs, and the regulatory relationship between them in the context of psoriasis, have still not been sufficiently investigated.
In summary, we demonstrated that the presence of IL36G may mediate the keratinocyte expression of AVPs such as MX1, ISG15, IFI27, and IFI44L in psoriasis vulgaris. The co-expression of AVP genes and IL36G was associated with psoriasis severity and therapeutic efficacy. Our results also suggest that infection factors and AVPs produced by keratinocytes may play an important role in the pathogenesis of psoriasis.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
Y-WL, H-MW, and Y-JC performed the data analysis and interpreted the results. NS and L-HY prepared the draft. R-JD, Y-QK, and Y-YL designed the research and revised the draft. All authors contributed to the article and approved the submitted version.
Funding
This work was supported in part by the National Natural Science Foundation of China (81860553), the Talent Introduction Project of Hubei Polytechnic University (21xjz33R, 21xjz34R), the “Ten-thousand Talents Program” of Yunnan Province (YNWR-MY-2018-039), the Project of AIDS Bureau of Yunnan Province, the Yunnan Province Clinical Center for Skin Immune Diseases (ZX2019-03-02), and the Yunnan Province Clinical Research Center for Skin Immune Diseases (2019ZF012).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.971071/full#supplementary-material
References
1. Griffiths CEM, Armstrong AW, Gudjonsson JE, Barker J. Psoriasis. Lancet (2021) 397(10281):1301–15. doi: 10.1016/S0140-6736(20)32549-6
2. Rendon A, Schakel K. Psoriasis pathogenesis and treatment. Int J Mol Sci (2019) 20(6):1475. doi: 10.3390/ijms20061475
3. Bartenjev I, Rogl Butina M, Potocnik M. Subclinical microbial infection in patients with chronic plaque psoriasis. Acta Derm Venereol Suppl (Stockh) (2000) 211):17–8. doi: 10.1080/00015550050500068
4. Ma JY, Shao S, Wang G. Antimicrobial peptides: Bridging innate and adaptive immunity in the pathogenesis of psoriasis. Chin Med J (Engl) (2020) 133(24):2966–75. doi: 10.1097/CM9.0000000000001240
5. Takahashi T, Gallo RL. The critical and multifunctional roles of antimicrobial peptides in dermatology. Dermatol Clin (2017) 35(1):39–50. doi: 10.1016/j.det.2016.07.006
6. Qian L, Chen W, Sun W, Li M, Zheng R, Qian Q, et al. Antimicrobial peptide ll-37 along with peptidoglycan drive monocyte polarization toward Cd14(High)Cd16(+) subset and may play a crucial role in the pathogenesis of psoriasis guttata. Am J Transl Res (2015) 7(6):1081–94.
7. Pahar B, Madonna S, Das A, Albanesi C, Girolomoni G. Immunomodulatory role of the antimicrobial ll-37 peptide in autoimmune diseases and viral infections. Vaccines (Basel) (2020) 8(3):517. doi: 10.3390/vaccines8030517
8. Zhou Y, Wang P, Yan BX, Chen XY, Landeck L, Wang ZY, et al. Quantitative proteomic profile of psoriatic epidermis identifies Oas2 as a novel biomarker for disease activity. Front Immunol (2020) 11:1432. doi: 10.3389/fimmu.2020.01432
9. Gao LJ, Shen J, Ren YN, Shi JY, Wang DP, Cao JM. Discovering novel hub genes and pathways associated with the pathogenesis of psoriasis. Dermatol Ther (2020) 33(6):e13993. doi: 10.1111/dth.13993
10. Raposo RA, Gupta R, Abdel-Mohsen M, Dimon M, Debbaneh M, Jiang W, et al. Antiviral gene expression in psoriasis. J Eur Acad Dermatol Venereol (2015) 29(10):1951–7. doi: 10.1111/jdv.13091
11. Leek JT, Johnson WE, Parker HS, Jaffe AE, Storey JD. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics (2012) 28(6):882–3. doi: 10.1093/bioinformatics/bts034
12. Rademaker M, Agnew K, Anagnostou N, Andrews M, Armour K, Baker C, et al. Psoriasis and infection. a clinical practice narrative. Australas J Dermatol (2019) 60(2):91–8. doi: 10.1111/ajd.12895
13. Garritsen FM, Kraag DE, de Graaf M. Guttate psoriasis triggered by perianal streptococcal infection. Clin Exp Dermatol (2017) 42(5):536–8. doi: 10.1111/ced.13129
14. Wang H, Chan HH, Ni MY, Lam WW, Chan WMM, Pang H. Bacteriophage of the skin microbiome in patients with psoriasis and healthy family controls. J Invest Dermatol (2020) 140(1):182–90 e5. doi: 10.1016/j.jid.2019.05.023
15. Langan EA, Kunstner A, Miodovnik M, Zillikens D, Thaci D, Baines JF, et al. Combined culture and metagenomic analyses reveal significant shifts in the composition of the cutaneous microbiome in psoriasis. Br J Dermatol (2019) 181(6):1254–64. doi: 10.1111/bjd.17989
16. Stehlikova Z, Kostovcik M, Kostovcikova K, Kverka M, Juzlova K, Rob F, et al. Dysbiosis of skin microbiota in psoriatic patients: Co-occurrence of fungal and bacterial communities. Front Microbiol (2019) 10:438. doi: 10.3389/fmicb.2019.00438
17. Heimlich DR, Harrison A, Mason KM. Host antimicrobial peptides in bacterial homeostasis and pathogenesis of disease. Antibiotics (Basel) (2014) 3(4):645–76. doi: 10.3390/antibiotics3040645
18. Morizane S, Gallo RL. Antimicrobial peptides in the pathogenesis of psoriasis. J Dermatol (2012) 39(3):225–30. doi: 10.1111/j.1346-8138.2011.01483.x
19. Zou A, Jian Q. Cxcl10 and its related key genes as potential biomarkers for psoriasis: Evidence from bioinformatics and real-time quantitative polymerase chain reaction. Med (Baltimore) (2021) 100(38):e27365. doi: 10.1097/MD.0000000000027365
20. Verhelst J, Hulpiau P, Saelens X. Mx proteins: Antiviral gatekeepers that restrain the uninvited. Microbiol Mol Biol Rev (2013) 77(4):551–66. doi: 10.1128/MMBR.00024-13
21. Schoggins JW, Wilson SJ, Panis M, Murphy MY, Jones CT, Bieniasz P, et al. A diverse range of gene products are effectors of the type I interferon antiviral response. Nature (2011) 472(7344):481–5. doi: 10.1038/nature09907
22. Freitas BT, Scholte FEM, Bergeron E, Pegan SD. How Isg15 combats viral infection. Virus Res (2020) 286:198036. doi: 10.1016/j.virusres.2020.198036
23. Hsieh WL, Huang YH, Wang TM, Ming YC, Tsai CN, Pang JH. Ifi27, a novel epidermal growth factor-stabilized protein, is functionally involved in proliferation and cell cycling of human epidermal keratinocytes. Cell Prolif (2015) 48(2):187–97. doi: 10.1111/cpr.12168
24. Ni X, Lai Y. Keratinocyte: A trigger or an executor of psoriasis? J Leukoc Biol (2020) 108(2):485–91. doi: 10.1002/JLB.5MR0120-439R
25. Benhadou F, Mintoff D, Del Marmol V. Psoriasis: Keratinocytes or immune cells - which is the trigger? Dermatology (2019) 235(2):91–100. doi: 10.1159/000495291
26. D’Erme AM, Wilsmann-Theis D, Wagenpfeil J, Holzel M, Ferring-Schmitt S, Sternberg S, et al. Il-36gamma (Il-1f9) is a biomarker for psoriasis skin lesions. J Invest Dermatol (2015) 135(4):1025–32. doi: 10.1038/jid.2014.532
27. Keermann M, Koks S, Reimann E, Prans E, Abram K, Kingo K. Transcriptional landscape of psoriasis identifies the involvement of Il36 and Il36rn. BMC Genomics (2015) 16:322. doi: 10.1186/s12864-015-1508-2
28. Ovesen SK, Schulze-Osthoff K, Iversen L, Johansen C. Ikbzeta 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(2):adv00386. doi: 10.2340/00015555-3749
29. Traks T, Keermann M, Prans E, Karelson M, Loite U, Koks G, et al. Polymorphisms in Il36g gene are associated with plaque psoriasis. BMC Med Genet (2019) 20(1):10. doi: 10.1186/s12881-018-0742-2
30. Lian LH, Milora KA, Manupipatpong KK, Jensen LE. The double-stranded rna analogue polyinosinic-polycytidylic acid induces keratinocyte pyroptosis and release of il-36gamma. J Invest Dermatol (2012) 132(5):1346–53. doi: 10.1038/jid.2011.482
Keywords: psoriasis, antiviral protein (AVP), IL36G, keratinocyte, treatment
Citation: Lu Y-W, Chen Y-J, Shi N, Yang L-H, Wang H-M, Dong R-J, Kuang Y-Q and Li Y-Y (2022) IL36G is associated with cutaneous antiviral competence in psoriasis. Front. Immunol. 13:971071. doi: 10.3389/fimmu.2022.971071
Received: 16 June 2022; Accepted: 22 August 2022;
Published: 12 September 2022.
Edited by:
Riichiro Abe, Niigata University, JapanReviewed by:
Sulev Kõks, Murdoch University, AustraliaKazumitsu Sugiura, Fujita Health University, Japan
Adam Reich, University of Rzeszow, Poland
Copyright © 2022 Lu, Chen, Shi, Yang, Wang, Dong, Kuang and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rong-Jing Dong, dongrongjing@hbpu.edu.cn; Yi-Qun Kuang, yq610433@hotmail.com; Yu-Ye Li, yyeli2000@126.com
†These authors have contributed equally to this work and share first authorship