Exploration of Potential Integrated Models of N6-Methyladenosine Immunity in Systemic Lupus Erythematosus by Bioinformatic Analyses

Systemic lupus erythematosus (SLE) is a prototypical systemic autoimmune disease of unknown etiology. The epigenetic regulation of N6-methyladenosine (m6A) modification in immunity is emerging. However, few studies have focused on SLE and m6A immune regulation. In this study, we aimed to explore a potential integrated model of m6A immunity in SLE. The models were constructed based on RNA-seq data of SLE. A consensus clustering algorithm was applied to reveal the m6A-immune signature using principal component analysis (PCA). Univariate and multivariate Cox regression analyses and Kaplan–Meier analysis were used to evaluate diagnostic differences between groups. The effects of m6A immune-related characteristics were investigated, including risk evaluation of m6A immune phenotype-related characteristics, immune cell infiltration profiles, diagnostic value, and enrichment pathways. CIBERSORT, ESTIMATE, and single-sample gene set enrichment analysis (ssGSEA) were used to evaluate the relative immune cell infiltrations (ICIs) of the samples. Conventional bioinformatics methods were used to identify key m6A regulators, pathways, gene modules, and the coexpression network of SLE. In summary, our study revealed that IGFBP3 (as a key m6A regulator) and two pivotal immune genes (CD14 and IDO1) may aid in the diagnosis and treatment of SLE. The potential integrated models of m6A immunity that we developed could guide clinical management and may contribute to the development of personalized immunotherapy strategies.