- 1Department of Rheumatology and Immunology, Shengjing Hospital, China Medical University, Shenyang, China
- 2Computational Systems Biology Laboratory, Department of Bioinformatics, Shantou University Medical College (SUMC), Shantou, China
Systemic lupus erythematosus (SLE) is a complex autoimmune disease that affects several organs and causes variable clinical symptoms. Exploring new insights on genetic factors may help reveal SLE etiology and improve the survival of SLE patients. The current study is designed to identify key genes involved in SLE and develop potential diagnostic biomarkers for SLE in clinical practice. Expression data of all genes of SLE and control samples in GSE65391 and GSE72509 datasets were downloaded from the Gene Expression Omnibus (GEO) database. A total of 11 accurate differentially expressed genes (DEGs) were identified by the “limma” and “RobustRankAggreg” R package. All these genes were functionally associated with several immune-related biological processes and a single KEGG (Kyoto Encyclopedia of Genes and Genome) pathway of necroptosis. The PPI analysis showed that IFI44, IFI44L, EIF2AK2, IFIT3, IFITM3, ZBP1, TRIM22, PRIC285, XAF1, and PARP9 could interact with each other. In addition, the expression patterns of these DEGs were found to be consistent in GSE39088. Moreover, Receiver operating characteristic (ROC) curves analysis indicated that all these DEGs could serve as potential diagnostic biomarkers according to the area under the ROC curve (AUC) values. Furthermore, we constructed the transcription factor (TF)-diagnostic biomarker-microRNA (miRNA) network composed of 278 nodes and 405 edges, and a drug-diagnostic biomarker network consisting of 218 nodes and 459 edges. To investigate the relationship between diagnostic biomarkers and the immune system, we evaluated the immune infiltration landscape of SLE and control samples from GSE6539. Finally, using a variety of machine learning methods, IFI44 was determined to be the optimal diagnostic biomarker of SLE and then verified by quantitative real-time PCR (qRT-PCR) in an independent cohort. Our findings may benefit the diagnosis of patients with SLE and guide in developing novel targeted therapy in treating SLE patients.
Introduction
Systemic lupus erythematosus (SLE) is a chronic autoimmune disease mediated by autoimmune reactions and is characterized by autoimmune inflammation (Bakshi et al., 2018). The serum of SLE patients contains a variety of autoantibodies represented by antinuclear antibodies, which combined with the corresponding autoantigens in the body to form immune complexes and then deposited in the glomerulus, liver, joints, skin, and other parts, causing multiple symptoms and manifestations (Ahlin et al., 2012). Currently, the pathogenesis and etiology of SLE have not been fully elucidated, and biomarkers available in clinical practice are still limited, mainly including anti-dsDNA antibodies, complement molecules, and white blood cell counts (Piga and Arnaud, 2021). Because of the diversity of symptoms, it can sometimes be difficult to distinguish SLE from other diseases with similar symptoms, such as rheumatoid arthritis (RA) and myositis. However, early diagnosis and timely intervention can help reduce SLE recurrence and hospitalization rates and the accumulation of chronic organ damage (Mak et al., 2013; Piga and Arnaud, 2021). Therefore, the identification of reliable biomarkers and revealing the underlying molecular mechanisms are essential for better diagnosis and effective treatment of SLE.
In recent years, numerous biomarkers in the different processes of autoimmune diseases have been identified through comprehensive bioinformatics analyses, including SLE (Zhao et al., 2021), RA (Cheng et al., 2021) and ulcerative Colitis (Chen et al., 2020), laying the foundations for exploring the potential molecular mechanisms in autoimmune diseases. Meanwhile, with the rapid development of artificial intelligence (AI), machine learning algorithm, as an important branch, has been widely used in the diagnostic classification and prognostic prediction of diseases. For example, machine learning methods were used to identify key prognostic molecules in esophageal squamous cell carcinoma (Li et al., 2021). Potential diagnostic biomarkers of acute myocardial infarction were identified by the least absolute shrinkage and selection operator (LASSO) regression model and support vector machine recursive feature elimination (SVM-RFE) (Zhao et al., 2020). Machine learning algorithms are generally divided into weak classifier algorithms and strong classifier algorithms. For example, logistic regression (LR), support vector machine (SVM), and artificial neural network (ANN) are weak classifier algorithms, and random forests (RF) and extreme gradient enhancement (XGBoost) are strong classifier algorithms. A strong classifier can be composed of more than one weak classifier. As proposed by Stafford et al., depending on the high classification performance of clinical and genomic data, RF and SVM were most frequently utilized in the diagnosis of autoimmune diseases (Stafford et al., 2020).
Therefore, based on bioinformatic analyses and machine learning algorithms, the present study was aimed to identify potential diagnostic biomarkers in SLE and construct the molecular regulatory networks related to diagnostic biomarkers, laying a foundation for in-depth exploration of molecular mechanisms of SLE.
Materials and Methods
Data Sources
The expression profiles of 924 SLE and 48 control samples in GSE65391, 99 SLE and 18 control samples in GSE72509, and 78 SLE and 46 control samples in GSE39088 were downloaded from the GEO database (https://www.ncbi.nlm.nih.gov/geo/). GSE65391 and GSE39088 were used to identify potential diagnostic biomarkers in SLE, while GSE39088 was used to test the reliability of these diagnostic biomarkers.
Screening and Functional Analysis of DEGs
DEGs between SLE and control samples in GSE65391 and GSE39088 datasets were identified by the “limma” R package with the threshold |log2FC| >1 and p-value < 0.05. The general views of DEGs were shown as volcano plots. The expressions of DEGs were visualized in the heatmap. The volcano plots and heatmaps were generated by the “ggplot2” R package. “RobustRankAggreg” R package was performed to screen accurate DEGs from GSE65391 and GSE39088 (Zhou et al., 2021). Then the functions of screened DEGs were analyzed by the “clusterProfiler” R package. p-value < 0.05 was considered as significantly enriched.
Identification of Diagnostical Biomarkers and Prediction of Regulators and Drugs
The screened DEGs were submitted into the STRING database (https://string-db.org), and then a PPI network was constructed by setting the confidence as 0.4. The correlations among DEGs in the PPI network were evaluated by the Pearson method and visualized in the heatmap. Then the expressions of correlated DEGs (Cor >0.7) were tested in GSE39088. Thereafter, the performance of these DEGs in distinguishing SLE and control samples were evaluated by ROC curves analyses, and DEGs with AUC greater than 0.7 were identified as diagnostic biomarkers of SLE. Furthermore, the reliability of these biomarkers was tested in GSE39088. The miRWalk database (http://mirwalk.umm.uni-heidelberg.de/) was used to predict the miRNAs targeting these diagnostic biomarkers. The ChEA3 database (https://amp.pharm.mssm.edu/chea3/) was used to predict the TFs targeting diagnostic biomarkers. Then the miRNA-diagnostic biomarker and TF-diagnostic biomarker pairs were integrated into a miRNA-diagnostic biomarker-TF regulatory network and visualized by Cytoscape software. The CTD database (http://ctdbase.org/) was used to search for drugs targeting these diagnostic biomarkers, and the PubChem database (http://www.pubchem.ncbi.nlm.gov) was used to display the chemical structures of several drugs, and a drug-gene network was constructed and visualized.
Evaluation of Immune Cell Infiltration
The immune infiltrations of 28 types of immune cells in SLE and control samples from GSE65391 were evaluated by the ssGSEA method (Ye et al., 2019). Differentially infiltrated immune cells between SLE and control samples were identified by the Wilcoxon test using p-value < 0.05 as the cutoff. Furthermore, the correlations between differentially infiltrated immune cells and diagnostic biomarkers were calculated by the Spearman method and shown in the heatmap.
Machine Learning
Machine learning methods, including LR, RF, XGBoost, SVM, and ANN (Li et al., 2021), were performed by the “glmnet”, “randomForest”, “xgboost”, “e1071”, and “neuralnet” R packages to develop classifiers for diagnostic classification, respectively. For each machine learning algorithm, 1,023 models representing all combinations of 10 identified biomarkers were established, and AUCs of the models were calculated. Among all classifiers, the top 100 models with the highest AUC values were selected, and the occurrence frequencies of each diagnostic biomarker were counted. The top five diagnostic biomarkers with the highest occurrence frequency in every classifier were extracted and overlapped to identify the optimal diagnostic biomarkers.
Collection of Clinical Characteristics
A total of 26 SLE patients and 20 sex- and age-matched healthy controls were recruited from Shengjing Hospital of China Medical University. The clinical characteristics of subjects, such as age, sex, course of the disease, and clinical and laboratory indices, were obtained from electronic medical records. The diagnosis of SLE was followed according to the European League Against Rheumatism (EULAR)/American College of Rheumatology (ACR) 2019 criteria and the SLE disease activity was assessed according to the systemic lupus erythematosus disease activity index 2000 (SLEDAI-2K) (Gladman et al., 2002). Patients with malignant tumors, pathogen infection, and other autoimmune diseases, such as RA and systemic sclerosis, were excluded. The research protocol was approved by the Medical Ethics Committees of the Shengjing Hospital of China Medical University. All experiments were conducted by the principles and regulations formulated by the ethics committee.
RNA Extraction and qRT-PCR
Blood samples from each subject were collected in ethylenediaminetetraacetic (EDTA) tubes. Peripheral blood mononuclear cells (PBMCs) were obtained by density gradient centrifugation (Solarbio Life Sciences, Beijing, China). Total cellular RNA was extracted from PBMCs using RNA Extraction Kit (Omega, Guangzhou, China). 200 ng RNA per sample was submitted for reverse transcription using Evo M-MLV RT Kit (Accurate Biotechnology, Changsha, China) following the manufacturer’s instructions. The qPCR cycle was conducted using SYBR Green Premix Pro Taq HS qPCR Kit (Accurate Biotechnology, Changsha, China) on Light Cycler 480 real-time PCR instrument (Roche, Basel, Switzerland). The 2−ΔΔCt method was used to calculate the relative expression of mRNA. β-actin was used as the internal control for normalization. The gene-specific primers are available in Supplementary Table S1.
Statistical Analysis
All the statistical analyses were performed with R software (version 4.1.0). The student’s t-test was performed to compare gene expressions between different groups. ROC curve analysis was used to evaluate the performance of biomarkers for diagnosing SLE. p < 0.05 was considered statistically significant.
Results
DEGs Involved in SLE
A total of 161 DEGs, including 124 up-regulated and 37 down-regulated genes in SLE samples relative to control samples, were identified in GSE65391 (Supplementary Table S2; Figure 1A), and the expression levels were shown as a heatmap (Figure 1B). Meantime, a total of 125 DEGs, including 116 up-regulated and nine down-regulated genes in SLE, were identified in GSE72509 (Supplementary Table S3; Figure 1C), and the expression levels were also shown as a heatmap (Figure 1D). By RobustRankAggreg method, IFI44, IFI44L, EIF2AK2, IFIT3, IFITM3, ZBP1, TRIM22, PRIC285, XAF1, PARP9, and ODF3B were screened as accurate DEGs (Supplementary Table S4). All the expressions of them were up-regulated in SLE samples (Figure 2A). The top 10 biological processes into which these DEGs were significantly enriched were associated with immunity, such as response to virus, type Ⅰ interferon signaling pathway, cellular response to type Ⅰ interferon, and positive regulation of cytokine-mediated signaling pathway (Figure 2B). In addition, these DEGs were markedly associated with a KEGG pathway of necroptosis (Figure 2C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/865559/fgene-13-865559-HTML/image_m/fgene-13-865559-g001.jpg)
FIGURE1. Volcano plots and heatmaps of DEGs. Each colored dot represents a DEG (|log2FC| >1 and p-value < 0.05). The blue dots represent the downregulated genes and the red dots represent the up-regulated genes. (A) Volcano map of the 161 EDGs identified in GSE65391. (B) Heatmap of the 161 EDGs identified in GSE65391. (C) Volcano map of the 125 EDGs identified in GSE72509. (D) Heatmap of the 125 EDGs identified in GSE72509.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/865559/fgene-13-865559-HTML/image_m/fgene-13-865559-g002.jpg)
FIGURE 2. The expression heatmap and enrichment analysis of accurate DEGs. (A) The expression heatmap of 11 accurate DEGs in GSE65391 and GSE72509. (B) The top 10 enriched GO terms. (C) The enriched pathway.
Ten Diagnostic Biomarkers Identified in SLE
Next, we constructed a PPI network of the DEGs, including IFI44, IFI44L, EIF2AK2, IFIT3, IFITM3, ZBP1, TRIM22, PRIC285 (also named as HELZ2), XAF1, and PARP9 (Figure 3A). These 10 genes had strong positive correlations with each other (Cor >0.7, Figure 3B), so we constructed a correlation network according to their correlations (Supplementary Table S5; Figure 3C). Moreover, the expression patterns of these 10 genes were validated in the GSE39088 dataset and were all up-regulated in SLE samples (Figure 3D). To identify their performance in distinguishing SLE and control samples, we plotted ROC curves in GSE65391 GSE72509 and GSE39088 datasets. The AUCs of the 10 genes were higher than 0.9 in GSE65391 (Figure 4A), higher than 0.85 in GSE72509 (Figure 4B), and higher than 0.7 in GSE39088, indicating that these genes had high accuracy and reliability in distinguishing between SLE and control samples. Thus, these 10 genes were identified as diagnostic biomarkers in SLE.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/865559/fgene-13-865559-HTML/image_m/fgene-13-865559-g003.jpg)
FIGURE 3. Interactions among the ten diagnostic biomarkers. (A) The PPI network of the ten biomarkers. (B) The correlation heatmap of the ten biomarkers. (C) The correlation network of the ten biomarkers. (D) The boxplot of the expression levels of ten biomarkers validated in GSE39088.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/865559/fgene-13-865559-HTML/image_m/fgene-13-865559-g004.jpg)
FIGURE4. The diagnostic performance of the ten biomarkers. (A)Training set: GSE65391. (B) Validation set: GSE72509. (C) Test set: GSE39088.
Construction of the TF-Diagnostic Biomarker-miRNA Network
Thereafter, we investigated the miRNAs and TFs that regulate the expression of the identified biomarkers. By miRWalk database, 176 miRNAs targeting IFI44, IFI44L, EIF2AK2, IFIT3, IFITM3, ZBP1, TRIM22, XAF1, and PARP9 were predicted, and a miRNA-biomarker network composed of 185 nodes and 179 edges was constructed (Supplementary Figure S1). In addition, 93 TFs binding with and regulating the expressions of IFI44, IFI44L, EIF2AK2, IFIT3, IFITM3, ZBP1, TRIM22, XAF1, and PARP9 were obtained from the ChEA3 database, and a TF-biomarker network composed of 102 nodes and 226 edges was constructed (Supplementary Figure S2). After integration, a TF-biomarker-miRNA regulatory network was visualized by Cytoscape, including 278 nodes and 405 edges (Figure 5). Furthermore, we also predicted the drugs targeting the biomarkers by CTD (Supplementary Table S6), extracted the drug-biomarker relation pairs, and constructed a drug-gene network composed of 218 nodes and 459 edges, including nine biomarkers and 209 drugs (Figure 6A). The chemical structures of several drugs, including (+)-JQ1 compound, acetaminophen, Benzo(a) pyrene, Estradiol, and Valproic Acid are displayed in Figure 6B.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/865559/fgene-13-865559-HTML/image_m/fgene-13-865559-g005.jpg)
FIGURE 5. TF-biomarker-miRNA regulatory network. The orange ellipses represent the biomarkers, the purple round rectangles represent the miRNAs and the green diamonds represent the TFs.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/865559/fgene-13-865559-HTML/image_m/fgene-13-865559-g006.jpg)
FIGURE 6. Predicting the drugs targeting the diagnostic biomarkers. (A) Drug-gene network. (B) Chemical structures of few of the drugs, from left to right, are (+)-JQ1 compound, acetaminophen, Benzo(A)pyrene, Estradiol and Valproic Acid.
Immune Cell Infiltration Results
Considering that SLE is an autoimmune disease and immune cells play important roles in affecting disease progression, we analyzed and compared the immune infiltration in SLE and control samples. The infiltration of 28 types of immune cells in each sample was calculated by the ssGSEA algorithm (Supplementary Table S7) and the results are displayed in the heatmap (Figure 7A). We found that the infiltrations of most immune cells were significantly different between SLE and control samples (Figure 7B). Moreover, we found that the expression of biomarker genes was positively correlated with activated dendritic cells, central memory CD8 T cells, gamma delta T cells, neutrophils, and type 2 T helper cells, but negatively correlated with activated B cells, activated CD8 T cells, CD56bright natural killer cells, CD56dim natural killer cells and central memory CD4 T cells (Figure 7C), indicating that these biomarkers may regulate SLE via interplay with the immune environment.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/865559/fgene-13-865559-HTML/image_m/fgene-13-865559-g007.jpg)
FIGURE 7. The relationship between diagnostic biomarkers and immune cell infiltration. (A) The heatmap of the infiltration proportion of 28 types of immune cells. (B)The boxplot of the infiltration proportion of 28 types of immune cells in SLE and control samples. (C) The heatmap of the correlations between diagnostic biomarkers and immune cells.
Identification of Optimal Biomarkers in SLE
To detect the optimal SLE biomarkers, we performed machine learning analyses, in which the importance of these 10 biomarkers was weighted by their occurrence frequencies in the top 100 models (Figure 8A). The top five important biomarkers in each machine learning method were shown in Table 1, and the only intersecting biomarker was IFI44 (Figure 8B), indicating that IFI44 was the optimal SLE biomarker. Then the IFI44 related network composed of 76 nodes and 75 edges was extracted, including one miRNA, 18 TFs, and 56 drugs (Figure 8C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/865559/fgene-13-865559-HTML/image_m/fgene-13-865559-g008.jpg)
FIGURE 8. Identification of the optimal diagnostic biomarkers. (A) The occurrence frequencies of ten diagnostic biomarkers in top 100 models. (B) Venn diagram of top five important biomarkers in each machine learning analysis. (C) The IFI44 related network.
To further explore whether IFI44 could be a reliable biomarker, the expression levels of IFI44 were identified by qRT-PCR in an independent cohort of 26 SLE patients and 20 healthy controls. The main clinical features of patients and controls are summarized in Table 2. Results showed that compared to that in healthy controls, the expression of IFI44 in SLE patients’ PBMCs was significantly higher no matter whether lupus nephritis (LN) was present (Figures 9A, B). ROC curve analysis showed that AUC for IFI44 was 0.850 when distinguishing SLE patients from healthy controls, and the diagnostic sensitivity and specificity were 0.923 and 0.850, respectively (Figure 9C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/865559/fgene-13-865559-HTML/image_m/fgene-13-865559-g009.jpg)
FIGURE 9. Validation of IFI44 as SLE diagnosis biomarker in an independent cohort. (A) Expression of IFI44 in SLE patients and healthy controls. (B) Expression of IFI44 in SLE patients with or without LN and healthy controls. (C) ROC curves of IFI44 for SLE diagnosis. ***p < 0.001, ****p < 0.0001.
Discussion
In this study, two gene expression profile datasets were integrated and analyzed by multiple bioinformatic approaches. 11 DEGs between SLE and control samples were identified and analyzed by GO and KEGG. The results of GO analysis in the biological process have mainly enriched the response of type I and type II interferon (IFN) and the regulation of cytokines, while KEGG pathway analysis showed these DEGs were involved in the pathway of necroptosis. There was sufficient evidence to support that the impaired expression of type I IFN and its related genes were widely involved in the pathology of SLE (Postal et al., 2020), and the activity of type I IFN is related to the level of circulating type II IFN (Oke et al., 2019). The dysregulated secretion of cytokines and associated impairment of immune regulation is a key factor influencing the symptoms and disease activity in SLE patients (Howe and Leung, 2019). Necroptosis, a specialized programmed cell death, is a regulated mode of necrotizing cell death mediated by the RIP1 and RIP3 kinases, the hyperactivation of which leads to autoimmunity (O'Donnell et al., 2018). It has been reported that elevated IFN signaling in SLE increased necroptosis, leading to tissue damage (Sarhan et al., 2019). Further, necroptosis was also associated with B cell reduction in SLE patients (Fan et al., 2014).
Through PPI analysis, we found that IFI44, IFI44L, EIF2AK2, IFIT3, IFITM3, ZBP1, TRIM22, PRIC285, XAF1, and PARP9, interacted with each other and had strong positive correlations, and all of them have high diagnostic efficiency in different datasets. Therefore, these 10 genes could be used as diagnostic biomarkers for SLE. The expression of IFI44 is induced by interferon, which may reflect the whole blood interferon signature in SLE (Strauß et al., 2017). DNA methylation is often inversely proportional to the transcriptional activity of genes, and PARP9 and IFI44L showed marked hypomethylation in a variety of immune cells in SLE patients (Ulff-Møller et al., 2018), including CD4+ T cells, monocytes, granulocytes, and B cells. Hypomethylation of the IFI44L promoter region has excellent sensitivity and specificity for diagnosing SLE and distinguishing it from other autoimmune diseases (Zhao et al., 2016). EIF2AK2 is highly expressed in SLE and selectively modulates immune responses and transcription of SLE-related histone genes by targeting TFs (Ge et al., 2021). IFIT3 promotes the production of type I IFN and other pro-inflammatory cytokines in SLE patients by positively regulating the cGAS-STING signaling pathway, aggravating the symptoms of SLE (Wang et al., 2018). IFITM3 is an interferon-induced transmembrane protein whose role in SLE is unclear, but it can inhibit the production of IL-6 (Stacey et al., 2017) and regulate the differentiation of T helper cells (Yánez et al., 2020), so it may help to regulate the inflammatory response and immune regulation in SLE. ZBP1 is closely related to necroptosis. Stimulated by IFN, ZBP1 protein can interact with RIPK3 to initiate RIPK3-dependent necroptosis (Yang et al., 2020). TRIM22 is a viral restriction factor that may play a role in certain autoimmune diseases such as multiple sclerosis (Jefferies et al., 2011). PRIC285 is a transcriptional coactivator involved in PPAR-γ signaling (Fairfax et al., 2012), and PPAR-γ can inhibit the activation of macrophages and regulate their differentiation, improving SLE symptoms (Kiss et al., 2013). XAF1 can increase p53 transcriptional activity (Pinto et al., 2020), and p53-dependent apoptosis has been implicated in the pathogenesis and disease activity of SLE (Chen et al., 2021).
The organ damage caused by SLE is attributed to the deposition of immune complexes on the one hand and the infiltration of activated immune cells on the other hand (Apostolidis et al., 2011). Therefore, we compared the immune cell infiltration between SLE and control samples in GSE65391 and found that the proportions of 18 immune cells were significantly different between the two groups, of which 10 were significantly associated with the screened diagnostic genes. Among them, the proportion of central memory CD8+ T cells was significantly high, and there was a strong positive correlation with diagnostic genes, which may be associated with the maintenance of chronic inflammation (Liu et al., 2007). It is generally believed that the Th2-dominated Th1/Th2 imbalance and the pathogenesis of SLE are intimately connected (Jiang et al., 2021), and our results showed that the diagnostic biomarkers have a significant positive correlation with the increased frequency of Th2, but not with Th1, which supported the reliability of the diagnostic biomarkers screened by us.
In addition, we constructed the miRNA-diagnostic biomarker-TF network to explore the regulatory mechanisms of the selected genes. MicroRNA (miRNA) is one of the main epigenetic regulators of SLE-related genes. A considerable amount of research progress has been made in the development of biomarkers and therapeutic methods based on miRNA (Hong et al., 2020). The network indicated that hsa-miR-6799-5p and hsa-miR-6759-5p, two tumor suppressor-related miRNAs, could interact with EIF2AK2 and IFI44L simultaneously. Hsa-miR-6759-5p can regulate the PI3K/AKT pathway that plays an important role in chronic inflammation by targeting IGF2 (Liu et al., 2020), while the molecular function of hsa-miR-6799-5p is still unclear. We also used the CTD database to predict the drugs associated with the diagnostic genes and established a drug-gene network, which can provide a reference for constructing new treatment options or mining potential pathogenic factors for SLE. For example, the (+)-JQ1 compound simultaneously targeted seven diagnostic genes in the network, and in vitro treatment of CD4+ T cells from SLE patients with JQ1 has been reported to reverse immune dysregulation and reduce inflammatory cytokines such as IFN-γ and IL-21 (Gao et al., 2018), suggesting that it could be a potential SLE therapeutic drug. The notorious carcinogen Benzo(a)pyrene also affects multiple diagnostic genes, but it is also an immunomodulator that can act as a ligand for aryl hydrocarbon receptors to alleviate arthritis symptoms in certain autoimmune diseases such as RA (Hui and Dai, 2020).
In this study, the RobustRankAggreg algorithm was used to evaluate the expression consistency of diagnostic genes in multiple datasets, and multiple machine learning algorithms were used to evaluate the contribution of different diagnostic genes to distinguish disease and control samples, to make the identification of biomarkers more accurate. Finally, the results showed that IFI44 had the highest contribution, suggesting that it may be the optimal SLE diagnostic biomarker. IFI44 is a type I IFN signature gene, which was hypomethylated in SLE patients (Joseph et al., 2019) and negatively regulated the innate immune response induced by the virus (DeDiego et al., 2019). Therefore, it may be related to the immune imbalance of autoimmune diseases. IFI44 has been considered as a key diagnostic biomarker in various diseases, including Sjogren’s syndrome (Xu et al., 2021) and psoriasis (Wang et al., 2020). A recent study showed that IFI44 can serve as a key biomarker for LN from IgA nephritis and healthy controls, and was associated with LN disease activity (Shen et al., 2021), suggesting that IFI44 was not only involved in the damage of immune complexes to the kidney but also closely related to the pathogenesis of SLE. Lupus nephritis is a frequent and severe complication of SLE, occurs in about 40% of SLE patients, which often indicates a poor prognosis (Gasparotto et al., 2020). The subjects included in our study were not limited to lupus nephritis, so it can better reflect the value of IFI44 in the diagnosis of SLE. The qRT-PCR results showed the up-regulation of IFI44 differed significantly between SLE patients with or without LN and healthy controls, indicating that IFI44 might be a reliable SLE diagnostic biomarker.
However, our study has certain limitations. Firstly, the samples of one of the datasets we used, GSE65391, were mainly from pediatric patients, and it is difficult to determine whether the age factor affected the research results. Secondly, we need further experiments to verify our findings, such as validating in a larger scale and rigorous trial, evaluating the expression of the other biomarkers we identified, and comparing the expression of IFI44 in various autoimmune diseases.
Conclusion
In conclusion, we found ten potential diagnostic biomarkers (IFI44, IFI44L, EIF2AK2, IFIT3, IFITM3, ZBP1, TRIM22, and PRIC285) for SLE by integrating bioinformatics methods, and discovered the potential of IFI44 as an optimal biomarker by five machine learning algorithms. The qRT-PCR and ROC curve analysis were performed to validate the diagnostic performance of IFI44 in an independent cohort. Immune cell infiltration showed the proportion of central memory CD8+ T cells was significantly high and positively correlated with selected biomarkers in SLE patients. The construction of miRNA-diagnostic biomarker-TF regulatory network and drug-gene network provides ideas for further exploring the pathogenesis at the genetic level and treatment of SLE.
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
ZJ and MS analyzed and interpreted data, analyzed the conclusions, and ZJ collected samples, performed experiments, and wrote the manuscript. XD and ZP reviewed the conclusions and revised the manuscript. DL and ZJ presented the idea and design. All authors contributed to the article and approved the submitted version.
Funding
The present study was supported by the National Natural Science Foundation of China (grant no. 81501408).
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/fgene.2022.865559/full#supplementary-material
References
Ahlin, E., Mathsson, L., Eloranta, M.-L., Jonsdottir, T., Gunnarsson, I., Rönnblom, L., et al. (2012). Autoantibodies Associated with RNA Are More Enriched Than Anti-dsDNA Antibodies in Circulating Immune Complexes in SLE. Lupus 21, 586–595. doi:10.1177/0961203311434938
Apostolidis, S. A., Lieberman, L. A., Kis-Toth, K., Crispín, J. C., and Tsokos, G. C. (2011). The Dysregulation of Cytokine Networks in Systemic Lupus Erythematosus. J. Interferon Cytokine Res. 31, 769–779. doi:10.1089/jir.2011.0029
Bakshi, J., Segura, B. T., Wincup, C., and Rahman, A. (2018). Unmet Needs in the Pathogenesis and Treatment of Systemic Lupus Erythematosus. Clinic Rev. Allerg Immunol. 55, 352–367. doi:10.1007/s12016-017-8640-5
Chen, Y., Li, H., Lai, L., Feng, Q., and Shen, J. (2020). Identification of Common Differentially Expressed Genes and Potential Therapeutic Targets in Ulcerative Colitis and Rheumatoid Arthritis. Front. Genet. 11, 572194. doi:10.3389/fgene.2020.572194
Chen, Y. C., Chou, Y. C., Hsieh, Y. T., Kuo, P. Y., Yang, M. L., Chong, H. E., et al. (2021). Targeting Intra-Pulmonary P53-Dependent Long Non-coding RNA Expression as a Therapeutic Intervention for Systemic Lupus Erythematosus-Associated Diffuse Alveolar Hemorrhage. Int. J. Mol. Sci. 22, 6948. doi:10.3390/ijms22136948
Cheng, Q., Chen, X., Wu, H., and Du, Y. (2021). Three Hematologic/immune System-specific Expressed Genes Are Considered as the Potential Biomarkers for the Diagnosis of Early Rheumatoid Arthritis through Bioinformatics Analysis. J. Transl Med. 19, 18. doi:10.1186/s12967-020-02689-y
DeDiego, M. L., Nogales, A., Martinez-Sobrido, L., and Topham, D. J. (2019). Interferon-Induced Protein 44 Interacts with Cellular FK506-Binding Protein 5, Negatively Regulates Host Antiviral Responses, and Supports Virus Replication. mBio 10, e01839–19. doi:10.1128/mBio.01839-19
Fairfax, B. P., Makino, S., Radhakrishnan, J., Plant, K., Leslie, S., Dilthey, A., et al. (2012). Genetics of Gene Expression in Primary Immune Cells Identifies Cell Type-specific Master Regulators and Roles of HLA Alleles. Nat. Genet. 44, 502–510. doi:10.1038/ng.2205
Fan, H., Liu, F., Dong, G., Ren, D., Xu, Y., Dou, J., et al. (2014). Activation-induced Necroptosis Contributes to B-Cell Lymphopenia in Active Systemic Lupus Erythematosus. Cell Death Dis 5, e1416. doi:10.1038/cddis.2014.375
Gao, X., Gao, K., Wu, J., and Zhao, M. (2018). Effect of JQ1 on Expression of Autoimmune-Related Genes in CD4+T Cells of Systemic Lupus Erythematosus. Zhong Nan Da Xue Xue Bao Yi Xue Ban 43, 704–710. doi:10.11817/j.issn.1672-7347.2018.07.002
Gasparotto, M., Gatto, M., Binda, V., Doria, A., and Moroni, G. (2020). Lupus Nephritis: Clinical Presentations and Outcomes in the 21st century. Rheumatology (Oxford) 59, v39–v51. doi:10.1093/rheumatology/keaa381
Ge, L., Zhang, Y., Zhao, X., Wang, J., Zhang, Y., Wang, Q., et al. (2021). EIF2AK2 Selectively Regulates the Gene Transcription in Immune Response and Histones Associated with Systemic Lupus Erythematosus. Mol. Immunol. 132, 132–141. doi:10.1016/j.molimm.2021.01.030
Gladman, D. D., Ibañez, D., and Urowitz, M. B. (2002). Systemic Lupus Erythematosus Disease Activity index 2000. J. Rheumatol. 29, 288–291.
Hong, S.-M., Liu, C., Yin, Z., Wu, L., Qu, B., and Shen, N. (2020). MicroRNAs in Systemic Lupus Erythematosus: a Perspective on the Path from Biological Discoveries to Clinical Practice. Curr. Rheumatol. Rep. 22, 17. doi:10.1007/s11926-020-00895-7
Howe, H. S., and Leung, B. P. L. (2019). Anti-Cytokine Autoantibodies in Systemic Lupus Erythematosus. Cells 9, 72. doi:10.3390/cells9010072
Hui, W., and Dai, Y. (2020). Therapeutic Potential of Aryl Hydrocarbon Receptor Ligands Derived from Natural Products in Rheumatoid Arthritis. Basic Clin. Pharmacol. Toxicol. 126, 469–474. doi:10.1111/bcpt.13372
Jefferies, C., Wynne, C., and Higgs, R. (2011). Antiviral TRIMs: Friend or Foe in Autoimmune and Autoinflammatory Disease? Nat. Rev. Immunol. 11, 617–625. doi:10.1038/nri3043
Jiang, Y., Zhao, Y., and Mo, X. (2021). Expression of lncRNA NEAT1 in Peripheral Blood Mononuclear Cells of Patients with Systemic Lupus Erythematosus and its Correlation with Th1/Th2 Balance. Int. J. Clin. Exp. Pathol. 14, 646–652.
Joseph, S., George, N. I., Green-Knox, B., Treadwell, E. L., Word, B., Yim, S., et al. (2019). Epigenome-wide Association Study of Peripheral Blood Mononuclear Cells in Systemic Lupus Erythematosus: Identifying DNA Methylation Signatures Associated with Interferon-Related Genes Based on Ethnicity and SLEDAI. J. Autoimmun. 96, 147–157. doi:10.1016/j.jaut.2018.09.007
Kiss, M., Czimmerer, Z., and Nagy, L. (2013). The Role of Lipid-Activated Nuclear Receptors in Shaping Macrophage and Dendritic Cell Function: From Physiology to Pathology. J. Allergy Clin. Immunol. 132, 264–286. doi:10.1016/j.jaci.2013.05.044
Li, M.-X., Sun, X.-M., Cheng, W.-G., Ruan, H.-J., Liu, K., Chen, P., et al. (2021). Using a Machine Learning Approach to Identify Key Prognostic Molecules for Esophageal Squamous Cell Carcinoma. BMC Cancer 21, 906. doi:10.1186/s12885-021-08647-1
Liu, G.-Z., Fang, L.-B., Hjelmström, P., and Gao, X.-G. (2007). Increased CD8+ central Memory T Cells in Patients with Multiple Sclerosis. Mult. Scler. 13, 149–155. doi:10.1177/1352458506069246
Liu, G., Guo, W., Chen, G., Li, W., Cui, Y., Qin, J., et al. (2020). Lnc-MCEI Mediated the Chemosensitivity of Esophageal Squamous Cell Carcinoma via miR-6759-5p to Competitively Regulate IGF2. Int. J. Biol. Sci. 16, 2938–2950. doi:10.7150/ijbs.47051
Mak, A., Isenberg, D. A., and Lau, C.-S. (2013). Global Trends, Potential Mechanisms and Early Detection of Organ Damage in SLE. Nat. Rev. Rheumatol. 9, 301–310. doi:10.1038/nrrheum.2012.208
O'Donnell, J. A., Lehman, J., Roderick, J. E., Martinez-Marin, D., Zelic, M., Doran, C., et al. (2018). Dendritic Cell RIPK1 Maintains Immune Homeostasis by Preventing Inflammation and Autoimmunity. J. Immunol. 200, 737–748. doi:10.4049/jimmunol.1701229
Oke, V., Gunnarsson, I., Dorschner, J., Eketjäll, S., Zickert, A., Niewold, T. B., et al. (2019). High Levels of Circulating Interferons Type I, Type II and Type III Associate with Distinct Clinical Features of Active Systemic Lupus Erythematosus. Arthritis Res. Ther. 21, 107. doi:10.1186/s13075-019-1878-y
Piga, M., and Arnaud, L. (2021). The Main Challenges in Systemic Lupus Erythematosus: Where Do We Stand? J. Clin. Med. 10, 243. doi:10.3390/jcm10020243
Pinto, E. M., Figueiredo, B. C., Chen, W., Galvao, H. C. R., Formiga, M. N., Fragoso, M. C. B. V., et al. (2020). XAF1 as a Modifier of P53 Function and Cancer Susceptibility. Sci. Adv. 6, eaba3231. doi:10.1126/sciadv.aba3231
Postal, M., Vivaldo, J. F., Fernandez-Ruiz, R., Paredes, J. L., Appenzeller, S., and Niewold, T. B. (2020). Type I Interferon in the Pathogenesis of Systemic Lupus Erythematosus. Curr. Opin. Immunol. 67, 87–94. doi:10.1016/j.coi.2020.10.014
Sarhan, J., Liu, B. C., Muendlein, H. I., Weindel, C. G., Smirnova, I., Tang, A. Y., et al. (2019). Constitutive Interferon Signaling Maintains Critical Threshold of MLKL Expression to License Necroptosis. Cell Death Differ. 26, 332–347. doi:10.1038/s41418-018-0122-7
Shen, L., Lan, L., Zhu, T., Chen, H., Gu, H., Wang, C., et al. (2021). Identification and Validation of IFI44 as Key Biomarker in Lupus Nephritis. Front. Med. 8, 762848. doi:10.3389/fmed.2021.762848
Stacey, M. A., Clare, S., Clement, M., Marsden, M., Abdul-Karim, J., Kane, L., et al. (2017). The Antiviral Restriction Factor IFN-Induced Transmembrane Protein 3 Prevents Cytokine-Driven CMV Pathogenesis. J. Clin. Invest. 127, 1463–1474. doi:10.1172/jci84889
Stafford, I. S., Kellermann, M., Mossotto, E., Beattie, R. M., Macarthur, B. D., and Ennis, S. (2020). A Systematic Review of the Applications of Artificial Intelligence and Machine Learning in Autoimmune Diseases. Npj Digit. Med. 3, 30. doi:10.1038/s41746-020-0229-3
Strauß, R., Rose, T., Flint, S. M., Klotsche, J., Häupl, T., Peck-Radosavljevic, M., et al. (2017). Type I Interferon as a Biomarker in Autoimmunity and Viral Infection: a Leukocyte Subset-specific Analysis Unveils Hidden Diagnostic Options. J. Mol. Med. (Berl) 95, 753–765. doi:10.1007/s00109-017-1515-7
Ulff-Møller, C. J., Asmar, F., Liu, Y., Svendsen, A. J., Busato, F., Grønbaek, K., et al. (2018). Twin DNA Methylation Profiling Reveals Flare-Dependent Interferon Signature and B Cell Promoter Hypermethylation in Systemic Lupus Erythematosus. Arthritis Rheumatol. 70, 878–890. doi:10.1002/art.40422
Wang, J., Dai, M., Cui, Y., Hou, G., Deng, J., Gao, X., et al. (2018). Association of Abnormal Elevations in IFIT 3 with Overactive Cyclic GMP ‐ AMP Synthase/Stimulator of Interferon Genes Signaling in Human Systemic Lupus Erythematosus Monocytes. Arthritis Rheumatol. 70, 2036–2045. doi:10.1002/art.40576
Wang, X., Liu, X., Liu, N., and Chen, H. (2020). Prediction of Crucial Epigenetically-associated, Differentially Expressed Genes by Integrated Bioinformatics Analysis and the identification of S100A9 as a Novel Biomarker in Psoriasis. Int. J. Mol. Med. 45, 93–102. doi:10.3892/ijmm.2019.4392
Xu, H., Chen, J., Wang, Y., Wu, Y., and Liang, Y. (2021). SELL and IFI44 as Potential Biomarkers of Sjögren's Syndrome and Their Correlation with Immune Cell Infiltration. Genes Genet. Syst. 96, 71–80. doi:10.1266/ggs.20-00053
Yánez, D. C., Ross, S., and Crompton, T. (2020). The IFITM Protein Family in Adaptive Immunity. Immunology 159, 365–372. doi:10.1111/imm.13163
Yang, D., Liang, Y., Zhao, S., Ding, Y., Zhuang, Q., Shi, Q., et al. (2020). ZBP1 Mediates Interferon-Induced Necroptosis. Cell Mol. Immunol. 17, 356–368. doi:10.1038/s41423-019-0237-x
Ye, L., Zhang, T., Kang, Z., Guo, G., Sun, Y., Lin, K., et al. (2019). Tumor-Infiltrating Immune Cells Act as a Marker for Prognosis in Colorectal Cancer. Front. Immunol. 10, 2368. doi:10.3389/fimmu.2019.02368
Zhao, M., Zhou, Y., Zhu, B., Wan, M., Jiang, T., Tan, Q., et al. (2016). IFI44L Promoter Methylation as a Blood Biomarker for Systemic Lupus Erythematosus. Ann. Rheum. Dis. 75, 1998–2006. doi:10.1136/annrheumdis-2015-208410
Zhao, E., Xie, H., and Zhang, Y. (2020). Predicting Diagnostic Gene Biomarkers Associated with Immune Infiltration in Patients with Acute Myocardial Infarction. Front. Cardiovasc. Med. 7, 586871. doi:10.3389/fcvm.2020.586871
Zhao, X., Zhang, L., Wang, J., Zhang, M., Song, Z., Ni, B., et al. (2021). Identification of Key Biomarkers and Immune Infiltration in Systemic Lupus Erythematosus by Integrated Bioinformatics Analysis. J. Transl Med. 19, 35. doi:10.1186/s12967-020-02698-x
Keywords: systemic lupus erythematosus, machine learning, integrated bioinformatics, diagnostic biomarkers, immune infiltration
Citation: Jiang Z, Shao M, Dai X, Pan Z and Liu D (2022) Identification of Diagnostic Biomarkers in Systemic Lupus Erythematosus Based on Bioinformatics Analysis and Machine Learning. Front. Genet. 13:865559. doi: 10.3389/fgene.2022.865559
Received: 30 January 2022; Accepted: 01 April 2022;
Published: 14 April 2022.
Edited by:
Shicheng Guo, University of Wisconsin-Madison, United StatesReviewed by:
Zhiqiu Xia, University of Nebraska Medical Center, United StatesFeng Cheng, University of California, San Francisco, United States
Songbo Wei, University of California, Riverside, United States
Shichang Li, Beckman Research Institute at City of Hope, United States
Copyright © 2022 Jiang, Shao, Dai, Pan and Liu. 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: Dongmei Liu, 18940256961@163.com