- 1Graduate School of Tianjin Medical University, Tianjin, China
- 2Graduate School of Tianjin University of Traditional Chinese Medicine, Tianjin, China
- 3The Second Affiliated Hospital of Guizhou Medical University, Kaili, Guizhou, China
- 4Tianjin Academy of Traditional Chinese Medicine Affiliated Hospital, Tianjin, China
Psoriasis is a chronic inflammatory skin disease, the etiology of which has not been fully elucidated, in which CD8+ T cells play an important role in the pathogenesis of psoriasis. However, there is a lack of in-depth studies on the molecular characterization of different CD8+ T cell subtypes and their role in the pathogenesis of psoriasis. This study aims to further expound the pathogenesy of psoriasis at the single-cell level and to explore new ideas for clinical diagnosis and new therapeutic targets. Our study identified a unique subpopulation of CD8+ T cells highly infiltrated in psoriasis lesions. Subsequently, we analyzed the hub genes of the psoriasis-specific CD8+ T cell subpopulation using hdWGCNA and constructed a machine-learning prediction model, which demonstrated good efficacy. The model interpretation showed the influence of each independent variable in the model decision. Finally, we deployed the machine learning model to an online website to facilitate its clinical transformation.
1 Introduction
Psoriasis is a usual chronic inflammatory skin disease with scaly erythema or plaques as the main clinical manifestations, long duration and easy to recur (Rapalli et al., 2020). The global prevalence of psoriasis is about 3% and is increasing year by year (Scher et al., 2019). Currently, the diagnosis of psoriasis is mainly based on the characteristics of skin lesions and skin imaging, while histopathology is needed to assist in the diagnosis of some atypical lesions, but due to the traumatizing nature of biopsy, more portable auxiliary diagnostic methods are needed in clinical working. At present, psoriasis is still incurable. Mild psoriasis is mainly treated with topical therapy, while moderate to severe psoriasis needs to be treated with combined systemic medication, and patients who have poor response to traditional therapeutic drugs need to choose biological agents or small molecule drug therapy. The pathogenesis of psoriasis is complex, and it is currently believed that it is caused by immune abnormalities induced by genetic and environmental factors (Liu et al., 2024), in which a variety of cells and the cytokines they release play an considerable role in the pathogenesy of psoriasis (Hawkes et al., 2017), but the current treatment methods are only aimed at some of the links, and thus have limited efficacy, and an in-depth understanding of the cellular and molecular pathogenesis of psoriasis can help to find new therapeutic targets to optimize the combination of treatment options. This is of great clinical significance.
Through the application of single-cell sequencing technology, we are granted the opportunity to delve deeply into the immunological microenvironment within psoriatic tissues. Meticulous analyses of cellular components in both diseased and normal tissues have revealed specific subpopulations exhibiting heightened expression within psoriatic tissues. Through cell communication analyses, we have unveiled the molecular mechanisms through which these subpopulations exert their influence within the tissue microenvironment, offering a novel perspective for a nuanced comprehension of the finely tuned dynamics of the immune system during disease progression. This not only provides unprecedented opportunities for personalized therapies but also establishes the groundwork for more precise interventions.
Ultimately, we introduce the power of machine learning, employing various algorithms to construct a diagnostic model for psoriasis based on single-cell-level information within bulk RNA sequencing (The flowchart is shown in Figure 1). This model not only accurately distinguishes psoriasis patients but also serves as an intelligent tool for clinicians (https://hesijia.shinyapps.io/Psoriasis/), facilitating early diagnosis and more personalized treatment approaches. This comprehensive research methodology opens new avenues for in-depth exploration of psoriasis, contributing positively to the enhancement of patients’ quality of life and treatment outcomes.
2 Materials and methods
2.1 Data acquisition
We downloaded the scRNA-seq dataset GSE151177 (13 human psoriasis skin and 5 healthy volunteer skin) from the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/) (Kim et al., 2021). Bulk RNA-seq dataset GSE54456 (92 psoriatic and 82 normal punch biopsies) were obtained from the GEO database (Li et al., 2014).
2.2 scRNA-seq data analysis
R package Seurat was utilized for conducting the scRNA-seq analysis. Cell clustering was implemented by the “FindNeighbors” and “FindCluster” functions of the Seurat package. We manually annotated cell types using the original literature. Intercellular communication networks with annotated scRNA-seq data were analyzed using the R package CellChat and examined for ligand-receptor-mediated interplays between diverse cell types. Assessing major signal inputs and outputs in cell subpopulations by CellChatDB.Human on cellchat. The functional enrichment analysis of cell types was performed by R package clusterProfiler.
2.3 hdWGCNA analysis
The hdWGCNA workflow in R was implemented using the “hdWGCNA” package (Morabito et al., 2021; Morabito et al., 2023). Briefly, the hdWGCNA pipeline consists of the following steps: 1) data preprocessing, 2) gene network construction, 3) module identification, 4) module preservation analysis. The first step is to preprocess the gene expression data to get rid of noise and batch effects. The next step is to construct the gene co-expression network according to the pairwise correlation between genes. The third step involves identifying modules or gene clusters of highly related genes and further calculate module genes. The fourth step is to perform a module preservation analysis to evaluate the stability of the identified modules.
2.4 Machine learning and deep learning
We applied univariate logistic regression to recognize key diagnostic genes in IS CD8+ T cells. In addition, the variables were filtered using the Least Absolute Shrinkage and Selection Operator Logistic Regression (LASSOLR) algorithm, Boruta and Random Forest (RF) (Kang et al., 2021; Younis et al., 2021). Machine learning models were developed by the tidymodels package (https://www.tidymodels.org/). Nine machine learning models were used to select the optimum model, including, Logistic Regression, Enet, k-nearest neighbor (KNN), RF, SVM, multilayer perceptron (MLP), LightGBM, XGBoost and decision tree (DT). All models were evaluated on the test data. SHapley Additive exPlanations (SHAP) was used for model interpretation (Wang et al., 2021). A convolutional neural network (CNN) based deep learning model was developed using the Keras and TensorFlow frameworks.
2.5 Statistical analysis
R version 4.2.2. was used for analysis. The statistical significance level was set at p < 0.05. Pearson correlation analysis was used to reveal the correlation between the 6 hub genes and 10 types of immune cells.
3 Results
3.1 Single-cell RNA sequencing analysis identifies CD8+ T cells upregulated in psoriasis skin lesions
To explore the differences in cellular composition between psoriatic lesions and normal skin tissues, we used GSE151177 to analyze the cellular component of the two. UAMP analysis identified nine major cell populations, and after integration through harmony, dimensionality reduction, clustering and cell type annotation, we found that the alteration in the proportion of CD8+ T cells was most significant in psoriatic lesions (Figures 2A, B). To further understand the functions of various cell types in psoriasis lesions, we performed gene enrichment analysis, which showed that CD8+ T cell-related genes were mainly involved in positive regulation of leukocyte activation, T cell receptor signaling pathway, αβT cell activation and other processes (Figure 2C). Keratinocyte (KC) -related genes located in different locations of the epidermis were involved in different response processes, for example, KC-related genes located in the basal layer were mainly involved in receptor-mediated endocytosis, KC-related genes located in the stratum spinosum were mainly involved in lipopolysaccharide response, and KC-related genes located in the stratum granulosum were mainly involved in protein refolding process. With the gradual migration of KC to the stratum corneum, KC-related genes play an important role in keratinization and epidermal development. In addition, mature dendritic cell (DC) -related genes were mainly involved in epidermal development and keratinization, while semi-mature DCs mainly mediated phagocytosis and leukocyte migration. Macrophage-related genes are mainly involved in antigen processing and presentation.
Figure 2. Heterogeneity of cellular composition in normal and psoriatic tissues revealed by single-cell RNA sequencing analysis according to dataset GSE151177. (A) UMAP scatter plot showing increased proportion of CD8+ T cells in psoriatic lesions. (B) Stacked bar graphs showing changes in the proportion of each cell type. (C) Functional enrichment plot showing the most significant enrichment terms for each cell type.
3.2 Identification of CD8+ T cell subpopulations in psoriatic lesions and analysis of cellular communication
To further elucidate the variations in the composition of CD8+ T cell subpopulations in psoriatic lesions, we further identified ten different CD8+ T cell subsets by dimensionality reduction and clustering, which showed that not all CD8+ T cells were upregulated in psoriatic lesion samples, there was a significant increase in CD8+ T cell subsets 0, 3, 4, 6, and 9, as compared to the normal skin tissue samples (Figures 3A, B), which resulted in an elevated proportion of CD8+ T cells. Therefore, we selected the five subpopulations that were specifically upregulated in psoriatic lesions for further analysis and defined them as the psoriasis-specific CD8+ T cell subpopulation (IS CD8+ T cells), and defined the other CD8+ T cell subpopulations as the non-psoriasis-specific CD8+ T cell subpopulation (Other CD8+ T cells). To further explore the cell-cell communication network between IS CD8+ T cells and other cell types, we performed intercellular communication analysis and found that IS CD8+ T cells had even more and closer connections with other cell types, particularly with melanocytes (MC) and Other CD8+ T cells (Figure 3C). Moreover, among the ten cell types, MC emitted the most signals while IS CD8+ T cells received the most signals (Figure 3D). IS CD8+ T cells emitted more IL16 signals compared to Other CD8+ T cells and this signal was mostly received by Semimature DC (Figures 3E, F). Based on the above results, we speculated that such enhanced signaling might enhance the antigen presenting ability of semimature DC cells and the ability to activate CD4+ T cells, thus promoting the immune response. In addition, IS CD8+ T cells received more TNF signals from KC located in the stratum granulosum than Other CD8+ T cells (Figure 3G). From this, we inferred that specifically receiving too much TNF signaling and sending too much IL16 signaling may make IS CD8+ T cells a pathogenic subpopulation in psoriasis. Therefore, targeted intervention of IL16 signaling emitted by IS CD8+ T cells, as well as interfering with the reception of TNF signaling by IS CD8+ T cells, may be potential therapeutic targets for psoriasis.
Figure 3. CD8+ T cell subset analysis and cell communication analysis (A) UMAP scatter plot showing the distribution of CD8+ T cell subtypes in psoriatic lesions and normal skin. (B) Stacked bar graph showing the proportion of CD8+ T cell subtypes in psoriatic lesions and normal skin. (C) The number and strength of intercellular communication between various cell subtypes. The colors of the bubbles and lines in the graph indicate different cellular sources. (D) Scatter plot showing the strength of input and output interactions of various cell subtypes. (E) Signaling role analysis on the aggregated cell-cell communication network of all signaling pathways among various cell subtypes. (F) Role of various cell subtypes in the IL-16 signaling pathway network. (G) Role of Various cell subtypes in the TNF signaling pathway network.
3.3 hdWGCNA shows that IS CD8+ T cells are characterized by brown, yellow, green, red, and turquoise modules
We further investigated the function and characteristics of IS CD8+ T cells using hdWGCNA. We chose a power value of 9 to construct a scale-free network, and 8 gene modules were generated accordingly (Figures 4A, B). Among these modules, brown, yellow, green, red, and turquoise modules were highly expressed in subpopulations 0, 3, 4, 6, and 9 (Figure 4C). We then screened the top 25 genes most associated with each color module (Figure 4D), for a total of 125 genes, for use to develop machine learning diagnostic models described below. To obtain a more comprehensive understanding of these 125 genes, we performed gene enrichment analysis with Metascape (https://metascape.org/). The results showed that these genes were mainly involved in cell biosynthesis, IL12 pathway, cytokine signaling and immune response (Figure 4E).
Figure 4. High dimensional weighted gene co-expression network analysis (hdWGCNA) further characterizes IS CD8+ T cells. (A) Power value equal to 9 when the network reached a scale-free distribution. (B) 8 modules were identified by hdWGCNA clustering tree. (C) Feather plots depict the scores of the corresponding modules in IS CD8+ T cells. (D) Network plots show the top 25 genes most associated with the five color modules. (E) The enrichment analysis diagram shows the most significant enrichment items of 125 genes.
3.4 Development of a machine learning diagnostic models on account of psoriasis-associated CD8+ T cell hub genes
As mentioned, we have identified five modules, including brown, yellow, red, green and turquoise. The top 25 genes in each module were selected, resulting in a total of 125 hub genes, which represent psoriasis-associated CD8+ T cell subsets hub genes. A total of 110 hub genes were included in the bulk RNA-seq cohort. Then, univariate analysis was performed on the 110 hub genes and 95 genes with p < 0.05 remained for the following analysis. Finally, three machine learning algorithms including Boruta, least absolute shrinkage and selection operator (LASSO) and RF were applied to feature selection (Figure 5A). Six features including GZMB, GNAS, GBP5, FOXP3, LSP1 and CD81 were obtained by interacting the three machine learning algorithms.
Figure 5. Development of machine learning models. (A) A three-step pipeline for hub genes selection. Step1: identification of five psoriasis-associated module genes. Step2: univariate analysis to screen the psoriasis-associated genes. Step3: Three machine learning algorithms for further selection of the psoriasis-associated genes. (B) Evaluation of robustness of the nine machine learning algorithms on test data. (C) The absolute value of the SHAP value of each feature. The higher SHAP value represents higher importance in predicting the presence of psoriasis. (D) Scatter plot of gene expression levels versus SHAP values for each feature in each sample.
Nine machine learning algorithms were used to fit models based on the six features, including Logistic Regression, Enet, KNN, RF, SVM, multilayer perceptron (MLP), LightGBM, XGBoost and decision tree (DT). We randomly divided the bulk RNA-seq data into training set and test set in a ratio of 7:3. Nine machine learning models were trained on the training set and comprehensively evaluated on the test set (Figure 5B). We noticed that nearly all models exhibit great performance on test set. To better understand the contribution of features to model predictions, we selected the logistic regression model to conduct model explanatory analysis using SHAP. Among the six features, GZMB plays the most important role on influencing the prediction results (Figure 5C). The relationship between the expression of each feature and the SHAP value is shown in Figure 5D. A SHAP value greater than 0 means that it has a positive contribution to predicting the outcome as psoriasis. As the expression level of GZMB increases, the SHAP value gradually increases, indicating that the high expression of GZMB in CD8+ T cells drives the pathogenesis of psoriasis. Similarly, decreased expression level of GNAS in CD8+ T cells also promotes the development of psoriasis (Figure 5D). In summary, we identified 6 markers and constructed machine learning models to improve psoriasis diagnosis. Additionally, we also provided novel targets for CD8+ T cell-targeted therapy in psoriasis patients. For example, the use of GZMB knockout CAR-T cells in patients may delay psoriasis progression.
3.5 Establishment and validation of the novel convolutional neural network deep learning model
Given that interactions exist between hub genes and other immune cells in the immune microenvironment, and these interactions may be associated with the progression of psoriasis. We developed a novel gene-immune image-based convolutional neural network (CNN) model. The advantage of the model is independent of the batch effect of different datasets. We constructed a heatmap for each patient, and the value of each square was the ratio of that patient’s gene expression to the infiltration of specific immune cells (Figure 6A). Based on the heatmaps, the CNN was trained using 100 epochs in the training set and tested in the testing set (Figure 6B). The CNN model performed well in both the training and testing sets (training AUC = 1, sensitivity = 100%, specificity = 100%, testing AUC = 1, sensitivity = 100%, specificity = 100%), suggesting good diagnostic performance (Figure 6C).
Figure 6. Establishment and validation of the novel CNN deep learning model. (A) Image of a patient used for training CNN model. (B) Training process of the CNN model. (C) Performance of the CNN model in the training and testing data.
3.6 Development of an online web tool to deploy the machine learning model
To make it easier for patients and clinicians to use the psoriasis diagnostic model we developed, we deployed the machine learning model on an online website (https://hesijia.shinyapps.io/Psoriasis/) (Figure 7).
4 Discussion
As a chronic inflammatory skin disease mediated by T-lymphocytes under the background of polygenic inheritance, psoriasis is stubborn and prone to recurrent episodes, and can lead to a variety of complications such as metabolic syndrome, cardiovascular disease, and arthritis (Griffiths et al., 2021), which critically affects the living quality of patients and their physical and mental health. The etiology of psoriasis has not been fully elucidated, among which genetic, environmental and immune factors play a significant role in the occurrence and development of psoriasis (Armstrong and Read, 2020). In recent years, the treatment of psoriasis with biologics has made great breakthroughs, but it still cannot meet the treatment needs of all psoriasis patients, suggesting that there are still blind spots in the underlying pathogenesis of psoriasis, so it is particularly important to gain a deeper understanding of the complex molecular mechanisms of psoriasis. In recent years, the development of single-cell RNA sequencing technology and machine learning has made it possible to utilize molecular genetic information to improve the precise diagnosis of psoriasis. In our study, we used single-cell RNA sequencing technology and hdWGCNA to identify psoriasis-specific CD8+ T cell subpopulations and screened psoriasis-specific CD8+ T cell signature genes, which may be used in the future to assist in the diagnosis of psoriasis and provide a theoretical basis for an in-depth understanding of the immune microenvironment of psoriasis.
CD8+ T cells play a crucial role in autoimmune diseases, anti-tumor immunity and anti-infection immunity. Previous researchs have revealed that CD8+ T cells infiltrate the skin lesions and blood of psoriasis patients and CD8+ T cells can produce IFN-γ, IL-17 and IL-22 (Hijnen et al., 2013), suggesting that CD8+ T cells play an important role in the development of psoriasis, but their role in psoriasis has not yet been fully elucidated, and very few studies have explored the specific subtypes and corresponding functions of CD8+ T cells. We identified a distinct subset of CD8+ T cells based on single-cell RNA sequencing, termed the psoriasis-associated CD8+ T cell subset, which is highly infiltrated in psoriatic lesions. We identified GZMB, GNAS, GBP5, FOXP3, LSP1 and CD81 as characteristic genes of psoriasis-associated CD8+ T cells, among which, Granzyme B (GzmB), a serine protease, is produced by natural killer (NK) cells and CD8+ T cells, and is an vital mediator of skin injury, inflammation and repair (Turner et al., 2019). Studies have shown that the expression of GZMB is elevated in skin lesions and plasma of psoriasis patients (Yawalkar et al., 2001; Kamata et al., 2016), suggesting that GzmB may be involved in CD8+ T cell-mediated cell damage and inflammation. The Gαs protein encoded by GNAS gene is a G protein-coupled receptor, which is involved in various signaling pathways, including cell growth, differentiation and apoptosis. The role of GNAS in psoriasis remains unclear. We speculate that GNAS may indirectly affect the development of psoriasis due to its role in immune response and cell proliferation. Guanylate binding protein 5 (GBP5), a member of the guanosine triphosphatase family, is involved in immune and inflammatory responses, and its expression is induced by type I and type II interferons (Kutsch and Coers, 2021), suggesting a potential role for the GBP5 in inflammatory skin diseases. We speculate that GBP5 may be involved in the pathogenesis of psoriasis by regulating inflammatory signaling pathways, such as IFN-γ signaling pathway. Forkhead box protein-3 (FOXP3) is a marker of regulatory T cells (Tregs), especially CD4+ Tregs, however, FOXP3 is also expressed in CD8+ T cells (Niederlova et al., 2021) and FOXP3 expression is increased in psoriasis lesions (Fujimura et al., 2008). Our results suggest that the high expression of GZBM, GBP5 and FOXP3 in CD8+ T cells drives the pathogenesis of psoriasis, therefore, inhibiting the expression of GZBM, GBP5 and FOXP3 is expected to be a new strategy to alleviate the progression of psoriasis. Lymphocyte Specific Protein 1 (LSP1) is mainly involved in cytoskeleton construction and cell migration, so we speculated that LSP1 may affect the migration and localization of inflammatory cells such as T cells and indirectly participate in the pathological process of psoriasis. CD81 is a transmembrane protein that can affect the proliferation and activation of a variety of immune cells. CD81 is highly expressed in both normal and psoriatic tissues, but when CD81 expression is excessively elevated, it promotes the development of psoriasis, which we hypothesize may be related to CD81-mediated hyper- or over-immunization, and this suggests that the immunohyperplasia of CD81 in CD8+ T cell subsets may mediate the occurrence of psoriasis, providing new clues for targeted therapy of psoriasis.
Currently, the diagnosis of psoriasis is primarily based on skin lesion characteristics and skin imaging, and there are fewer studies on molecular markers for the early identification of psoriasis patients. In recent years, machine learning has shown great potential in clinical image recognition and disease classification (Anwar et al., 2018), and in the field of dermatology, machine learning algorithms have been used in a variety of diseases such as melanoma, basal cell carcinoma, and onychomycosis, etc. (Haenssle et al., 2018; Harangi, 2018; Marchetti et al., 2018; Zhang et al., 2018), which makes it possible to assist in the diagnosis of diseases by using machine learning methods, of which the first results have been seen in a number of aspects of psoriasis diagnosis, treatment and management (Patrick et al., 2018; Yu et al., 2020; Hammad et al., 2023). Verma et al. (2019) used a machine learning algorithm to predict skin diseases such as psoriasis and achieved a maximum accuracy of 98.64%. Kim et al. (2019) applied machine learning to spectral classification and diagnosis of psoriasis, demonstrating 95% accuracy, 95% sensitivity, and 100% specificity. Dash et al. (2019) constructed a Deep Convolutional Neural Network (DCNN) model to segment psoriasis lesions to assist in the diagnosis of psoriasis. The accuracy of the model was 94.8%, the sensitivity was 89.6%, and the specificity was 97.6%. Due to the diversity and large differences of psoriasis skin lesions, our study obtains six characteristic genes through three machine learning algorithms to predict the occurrence of psoriasis and assist clinical diagnosis and treatment, and shows 100% accuracy, which is higher than previous machine learning image recognition and other methods, and is more conducive to the disease management of psoriasis.
In conclusion, our study identified a psoriasis-specific CD8+ T cell subset and its characteristic genes, and established a machine-learning prediction model that provides new ideas for future psoriasis diagnosis. This study is unique in that it delves into the pathogenesis of psoriasis at the single-cell level, identifying a new cell subset specific to psoriasis. Through advanced bioinformatics methods and machine learning/deep learning techniques, it provides new big data-driven tools for the diagnosis of psoriasis, as well as potential new targets for clinical treatment.
Data availability statement
The data presented in the study are deposited in the Gene Expression Omnibus (GEO) repository, accession number GSE151177 and GSE54456.
Author contributions
SH: Conceptualization, Investigation, Writing–original draft, Writing–review and editing. LL: Data curation, Methodology, Writing–original draft, Writing–review and editing. XL: Formal Analysis, Writing–original draft. MG: Supervision, Writing–review and editing. MC: Project administration, Writing–original draft. JZ: Funding acquisition, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
We thank each of the authors for their hard work in the article.
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.
References
Anwar, S. M., Majid, M., Qayyum, A., Awais, M., Alnowami, M., and Khan, M. K. (2018). Medical image analysis using convolutional neural networks: a review. J. Med. Syst. 42 (11), 226. doi:10.1007/s10916-018-1088-1
Armstrong, A. W., and Read, C. (2020). Pathophysiology, clinical presentation, and treatment of psoriasis: a review. JAMA 323 (19), 1945–1960. doi:10.1001/jama.2020.4006
Dash, M., Londhe, N. D., Ghosh, S., Semwal, A., and Sonawane, R. S. (2019). PsLSNet: automated psoriasis skin lesion segmentation using modified U-Net-based fully convolutional network. Biomed. Signal Process. Control 52, 226–237. doi:10.1016/j.bspc.2019.04.002
Fujimura, T., Okuyama, R., Ito, Y., and Aiba, S. (2008). Profiles of Foxp3+ regulatory T cells in eczematous dermatitis, psoriasis vulgaris and mycosis fungoides. Br. J. Dermatol. 158 (6), 1256–1263. doi:10.1111/j.1365-2133.2008.08504.x
Griffiths, C., Armstrong, A., Gudjonsson, J., and Barker, J. (2021). Psoriasis. Lancet 397 (10281), 1301–1315. doi:10.1016/S0140-6736(20)32549-6
Haenssle, H. A., Fink, C., Schneiderbauer, R., Toberer, F., Buhl, T., Blum, A., et al. (2018). Man against machine: diagnostic performance of a deep learning convolutional neural network for dermoscopic melanoma recognition in comparison to 58 dermatologists. Ann. Oncol. 29 (8), 1836–1842. doi:10.1093/annonc/mdy166
Hammad, M., Pławiak, P., ElAffendi, M., El-Latif, A. A. A., and Latif, A. A. A. (2023). Enhanced deep learning approach for accurate eczema and psoriasis skin detection. Sensors 23 (16), 7295. doi:10.3390/s23167295
Harangi, B. (2018). Skin lesion classification with ensembles of deep convolutional neural networks. J. Biomed. Inf. 86, 25–32. doi:10.1016/j.jbi.2018.08.006
Hawkes, J. E., Chan, T. C., and Krueger, J. G. (2017). Psoriasis pathogenesis and the development of novel targeted immune therapies. J. Allergy Clin. Immunol. 140 (3), 645–653. doi:10.1016/j.jaci.2017.07.004
Hijnen, D., Knol, E. F., Gent, Y. Y., Giovannone, B., Beijn, S. J., Kupper, T. S., et al. (2013). CD8+ T cells in the lesional skin of atopic dermatitis and psoriasis patients are an important source of IFN-γ, IL-13, IL-17, and IL-22. J. Invest. Dermatol. 133 (4), 973–979. doi:10.1038/jid.2012.456
Kamata, Y., Kimura, U., Matsuda, H., Tengara, S., Kamo, A., Umehara, Y., et al. (2016). Relationships among plasma granzyme B level, pruritus and dermatitis in patients with atopic dermatitis. J. Dermatol. Sci. 84 (3), 266–271. doi:10.1016/j.jdermsci.2016.09.009
Kang, J., Choi, Y. J., Kim, I. K., Lee, H. S., Kim, H., Baik, S. H., et al. (2021). LASSO-based machine learning algorithm for prediction of lymph node metastasis in T1 colorectal cancer. Cancer Res. Treat. 53 (3), 773–783. doi:10.4143/crt.2020.974
Kim, J., Lee, J., Kim, H. J., Kameyama, N., Nazarian, R., Der, E., et al. (2021). Single-cell transcriptomics applied to emigrating cells from psoriasis elucidate pathogenic versus regulatory immune cell subsets. J. Allergy Clin. Immunol. 148 (5), 1281–1292. doi:10.1016/j.jaci.2021.04.021
Kim, S., Kim, J., Hwang, M., Kim, M., Jo, S. J., Je, M., et al. (2019). Smartphone-based multispectral imaging and machine-learning based analysis for discrimination between seborrheic dermatitis and psoriasis on the scalp. Biomed. Opt. Express 10 (2), 879–891. doi:10.1364/BOE.10.000879
Kutsch, M., and Coers, J. (2021). Human guanylate binding proteins: nanomachines orchestrating host defense. FEBS J. 288 (20), 5826–5849. doi:10.1111/febs.15662
Li, B., Tsoi, L. C., Swindell, W. R., Gudjonsson, J. E., Tejasvi, T., Johnston, A., et al. (2014). Transcriptome analysis of psoriasis in a large case–control sample: RNA-seq provides insights into disease mechanisms. J. Invest. Dermatol. 134 (7), 1828–1838. doi:10.1038/jid.2014.28
Liu, S., He, M., Jiang, J., Duan, X., Chai, B., Zhang, J., et al. (2024). Triggers for the onset and recurrence of psoriasis: a review and update. Cell Commun. Signal. 22 (1), 108. doi:10.1186/s12964-023-01381-0
Marchetti, M. A., Codella, N. C., Dusza, S. W., Gutman, D. A., Helba, B., Kalloo, A., et al. (2018). Results of the 2016 international skin imaging collaboration international symposium on biomedical imaging challenge: comparison of the accuracy of computer algorithms to dermatologists for the diagnosis of melanoma from dermoscopic images. J. Am. Acad. Dermatol. 78 (2), 270–277.e1. doi:10.1016/j.jaad.2017.08.016
Morabito, S., Miyoshi, E., Michael, N., Shahin, S., Martini, A. C., Head, E., et al. (2021). Single-nucleus chromatin accessibility and transcriptomic characterization of Alzheimer’s disease. Nat. Genet. 53 (8), 1143–1155. doi:10.1038/s41588-021-00894-z
Morabito, S., Reese, F., Rahimzadeh, N., Miyoshi, E., and Swarup, V. (2023). hdWGCNA identifies co-expression networks in high-dimensional transcriptomics data. Cell Rep. Methods 3 (6), 100498. doi:10.1016/j.crmeth.2023.100498
Niederlova, V., Tsyklauri, O., Chadimova, T., and Stepanek, O. (2021). CD8+ Tregs revisited: a heterogeneous population with different phenotypes and properties. Eur. J. Immunol. 51 (3), 512–530. doi:10.1002/eji.202048614
Patrick, M. T., Stuart, P. E., Raja, K., Gudjonsson, J. E., Tejasvi, T., Yang, J., et al. (2018). Genetic signature to provide robust risk assessment of psoriatic arthritis development in psoriasis patients. Nat. Commun. 9 (1), 4178. doi:10.1038/s41467-018-06672-6
Rapalli, V. K., Waghule, T., Gorantla, S., Dubey, S. K., Saha, R. N., and Singhvi, G. (2020). Psoriasis: pathological mechanisms, current pharmacological therapies, and emerging drug delivery systems. Drug Discov. Today 25 (12), 2212–2226. doi:10.1016/j.drudis.2020.09.023
Scher, J. U., Ogdie, A., Merola, J. F., and Ritchlin, C. (2019). Preventing psoriatic arthritis: focusing on patients with psoriasis at increased risk of transition. Nat. Rev. Rheumatol. 15 (3), 153–166. doi:10.1038/s41584-019-0175-0
Turner, C. T., Lim, D., and Granville, D. J. (2019). Granzyme B in skin inflammation and disease. Matrix Biol. 75, 126–140. doi:10.1016/j.matbio.2017.12.005
Verma, A. K., Pal, S., and Kumar, S. (2019). Classification of skin disease using ensemble data mining techniques. Asian pac. J. Cancer Prev. 20 (6), 1887–1894. doi:10.31557/APJCP.2019.20.6.1887
Wang, K., Tian, J., Zheng, C., Yang, H., Ren, J., Liu, Y., et al. (2021). Interpretable prediction of 3-year all-cause mortality in patients with heart failure caused by coronary heart disease based on machine learning and SHAP. Comput. Biol. Med. 137, 104813. doi:10.1016/j.compbiomed.2021.104813
Yawalkar, N., Schmid, S., Braathen, L., and Pichler, W. (2001). Perforin and granzyme B may contribute to skin inflammation in atopic dermatitis and psoriasis. Br. J. Dermatol. 144 (6), 1133–1139. doi:10.1046/j.1365-2133.2001.04222.x
Younis, H., Anwar, M. W., Khan, M. U. G., Sikandar, A., and Bajwa, U. I. (2021). A new sequential forward feature selection (SFFS) algorithm for mining best topological and biological features to predict protein complexes from protein–protein interaction networks (PPINs). Interdiscip. Sci. 13 (3), 371–388. doi:10.1007/s12539-021-00433-8
Yu, K., Syed, M. N., Bernardis, E., and Gelfand, J. M. (2020). Machine learning applications in the evaluation and management of psoriasis: a systematic review. J. Psoriasis Psoriatic Arthritis 5 (4), 147–159. doi:10.1177/2475530320950267
Keywords: single-cell RNA sequencing, hdWGCNA, machine learning, psoriasis, CD8+ T cells, predictive modeling
Citation: He S, Liu L, Long X, Ge M, Cai M and Zhang J (2024) Single-cell analysis and machine learning identify psoriasis-associated CD8+ T cells serve as biomarker for psoriasis. Front. Genet. 15:1387875. doi: 10.3389/fgene.2024.1387875
Received: 18 February 2024; Accepted: 24 May 2024;
Published: 10 June 2024.
Edited by:
Deborah Giordano, National Research Council (CNR), ItalyReviewed by:
Francesco Bardozzo, University of Salerno, ItalyJingtao Qiu, Stanford University, United States
Daniela Evangelista, National Research Council (CNR), Italy
Copyright © 2024 He, Liu, Long, Ge, Cai and Zhang. 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: Junling Zhang, emhhbmdqdW5saW5ndGpAMTI2LmNvbQ==
†These authors have contributed equally to this work