Abstract
Systemic lupus erythematosus (SLE) is an autoimmune disease involving multiple systems. Its recurrent episodes and fluctuating disease courses have a severe impact on patients. Biomarkers to predict disease prognosis and remission are still lacking in SLE. We downloaded the GSE50772 dataset from the Gene Expression Omnibus database and identified differentially expressed genes (DEGs) between SLE and healthy controls. Weighted gene co-expression network analysis was used to identify key gene modules and corresponding genes in SLE. The overlapped genes in DEGs and key modules are used as key genes for subsequent analysis. These key genes were analyzed using 3 machine learning algorithms, including the least absolute shrinkage and selection operator, support vector machine recursive elimination, and random forest algorithms. The overlapped genes were obtained as potential biomarkers for further analysis, investigating and validating the potential biomarkers' possible functions, regulatory mechanisms, diagnostic value, and expression levels. And finally studied the differences between groups in level of immune cell infiltration and explored the relationship between potential biomarkers and immunity. A total of 234 overlapped genes in DEGs and key modules are used as key genes for subsequent analysis. After taking the intersection of the key genes obtained by 3 algorithms, we got 4 potential biomarkers (ARID2, CYSTM1, DDIT3, and RNASE1) with high diagnostic values. Finally, further immune infiltration analysis showed differences in various immune cells in the SLE and healthy control samples. ARID2, CYSTM1, DDIT3, and RNASE1 can affect the immune function of SLE patients. ARID2, CYSTM1, DDIT3, and RNASE1 could be used as immune-related potential biomarkers and therapeutic or diagnostic targets for further research.