B cell subsets-related biomarkers and molecular pathways for systemic lupus erythematosus by transcriptomics analyses

Systemic lupus erythematosus (SLE), an autoimmune disease, is characterised by B-cell abnormalities and a loss of tolerance that can produce autoantibody. However, the imperative genes and molecular pathways involved in the change of B cell populations remain unclear. The expression of B cell subsets between SLE and healthy controls (HCs) was detected based on micro-array transcriptome data. The Weighted Gene Co-Expression Network Analysis (WGCNA) further revealed the co-expression modules of naïve and memory B cells. Whereafter, we performed the functional enrichment analysis, Protein-protein interaction (PPI) networks construction and feature selection to screen hub genes. Ultimately, we recruited SLE patients and HCs from the Second Hospital of Shanxi Medical University and further verified these genes in transcriptome sequencing samples. Total of 1087 SLE patients and 86 HCs constituted in the study. Compared to HCs, the levels of peripheral naïve B cells of SLE patients decreased, while memory B cells increased. WGCNA identified two modules with the highest correlation for the subsequent analysis. The purple module was primarily in connection with naïve B cells, and the GO analysis indicated that these genes were mainly abundant in B cell activation. The blue module relevant to memory B cells was most significantly enriched in the “defence response to virus” correlation pathway. Then we screened six hub genes by PPI and feature selection. Finally, four biomarkers (IFI27, IFITM1, MX2, IRF7) were identified by transcriptome sequencing verification. Our study identified hub genes and key pathways associated with the naïve and memory B cells respectively, which may offer novel insights into the behaviours of B cells and the pathogenesis of SLE.