Identification of potential biomarkers and immunoregulatory mechanisms of rheumatoid arthritis based on multichip co-analysis of GEO database

doi:10.12122/j.issn.1673-4254.2024.06.10

Abstract

Abstract:

Objective To identify the biomarkers for early rheumatoid arthritis (RA) diagnosis and explore the possible immune regulatory mechanisms. Methods The differentially expressed genesin RA were screened and functionally annotated using the limma, RRA, batch correction, and clusterProfiler. The protein-protein interaction network was retrieved from the STRING database, and Cytoscape 3.8.0 and GeneMANIA were used to select the key genes and predicting their interaction mechanisms. ROC curves was used to validate the accuracy of diagnostic models based on the key genes. The disease-specific immune cells were selected via machine learning, and their correlation with the key genes were analyzed using Corrplot package. Biological functions of the key genes were explored using GSEA method. The expression of STAT1 was investigated in the synovial tissue of rats with collagen-induced arthritis (CIA). Results We identified 9 core key genes in RA (CD3G, CD8A, SYK, LCK, IL2RG, STAT1, CCR5, ITGB2, and ITGAL), which regulate synovial inflammation primarily through cytokines-related pathways. ROC curve analysis showed a high predictive accuracy of the 9 core genes, among which STAT1 had the highest AUC (0.909). Correlation analysis revealed strong correlations of CD3G, ITGAL, LCK, CD8A, and STAT1 with disease-specific immune cells, and STAT1 showed the strongest correlation with M1-type macrophages (R=0.68, P=2.9e-08). The synovial tissues of the ankle joints of CIA rats showed high expressions of STAT1 and p-STAT1 with significant differential expression of STAT1 between the nucleus and the cytoplasm of the synovial fibroblasts. The protein expressions of p-STAT1 and STAT1 in the cell nuclei were significantly reduced after treatment. Conclusion CD3G, CD8A, SYK, LCK, IL2RG, STAT1, CCR5, ITGB2, and ITGAL may serve as biomarkers for early diagnosis of RA. Gene-immune cell pathways such as CD3G/CD8A/LCK-γδ T cells, ITGAL-Tfh cells, and STAT1-M1-type macrophages may be closely related with the development of RA.

Key words: rheumatoid arthritis, biomarkers, immune regulatory, bionformatic, machine learning

Lili CHEN, Tianyu WU, Ming ZHANG, Zixia DING, Yan ZHANG, Yiqing YANG, Jiaqian ZHENG, Xiaonan ZHANG. Identification of potential biomarkers and immunoregulatory mechanisms of rheumatoid arthritis based on multichip co-analysis of GEO database[J]. Journal of Southern Medical University, 2024, 44(6): 1098-1108.

Figures/Tables 8

Fig.1 Differential gene expressions in different datasets of RA. A-E: RA-related differentially expressed genes in the GSE1919, GSE12021, GSE55235, GSE55457 and GSE77298 datasets.

Fig.2 GO and KEGG enrichment analysis of the differentially expressed genes (DEGs). A: Circle diagram for GO enrichment analysis of the DEGs. B: Sankey diagrams for KEGG enrichment analysis of the DEGs.

Fig.3 Screening of the differentially expressed core genes and their expressions. A: Results of the core gene screening. B: Complex Heatmap of the core genes. C-K: Expression of the core genes in control, early RA and advanced RA in the validation dataset (C: SYK; D: STAT1; E: LCK; F: IL2RG; G: ITGB2; H: ITGAL; I: CD3G; J: CCR5; K: CD8A).

Fig.4 Core gene function prediction and receiver operating characteristic (ROC) curve results. A:Analysis of GeneMANIA results for core genes. B-J: Area under ROC for the core genes (B: CD8A; C: CD3G; D: IL2RG; E: ITGB2; F: ITGAL; G: LCK; H: SYK; I: CCR5; J: STAT1).

Fig.5 Disease-specific immune cell expression and screening. A: Immune cell infiltration in RA. B: Differential expression of immune cells in RA. C: Correlation Heatmap of the immune cells. D: Intersections of disease-specific cells obtained by screening with Lasso regression and CIBERSORT methods.

Fig.6 Correlation analysis of the core genes with the disease-specific immune cells. A: Heatmap of correlations between core genes and disease-specific immune cells. B-F: Correlation analysis of the core genes with disease-specific immune cells (B: CD8A with γδ T cell; C: CD3G with γδ T cell; D: ITGAL with Tfh cell; E: LCK with γδ T cell; F: STAT1 with macrophages M1).

Fig.7 Expression of STAT1 in ankle joints of CIA rats. A: GSEA analysis of biological processes potentially mediated by upregulated STAT1. B: HE staining of the knee joints of CIA rats. C: Safranine O-Fast Green staining of the knee joints. D: STAT1 expression in the synovial tissue of CIA rats.

Fig.8 Expression of p-STAT1 in the synovial tissue and fibroblast-like synoviocytes of CIA rats. A: p-STAT1 expression in the synovial tissue of CIA rats. B: Immunoblotting showing expressions of STAT1/p-STAT1 in the cytoplasm and nucleus of the synovial tissue cells in different groups. C, D: Western blotting of p-STAT and STAT1 proteins in cytoplasm and nuclei of synovial fibroblasts in the joint tissue of CIA rats. *P<0.05,**P<0.01,***P<0.001, ****P<0.0001.

References 39

1	Alivernini S, Firestein GS, McInnes IB. The pathogenesis of rheumatoid arthritis[J]. Immunity, 2022, 55(12): 2255-70.
2	Xu Q, Ni JJ, Han BX, et al. Causal relationship between gut microbiota and autoimmune diseases: a two-sample Mendelian randomization study[J]. Front Immunol, 2021, 12: 746998.
3	Kronzer VL, Davis JM. Etiologies of rheumatoid arthritis: update on mucosal, genetic, and cellular pathogenesis[J]. Curr Rheumatol Rep, 2021, 23(4): 21.
4	Chiang PH, Ju PC, Chiang YC, et al. Correspondence on ‘Incidence trend of five common musculoskeletal disorders from 1990 to 2017 at the global, regional and national level: results from the global burden of disease study 2017'[J]. Ann Rheum Dis, 2023, 82(2): e46.
5	Guan SY, Zheng JX, Sam NB, et al. Global burden and risk factors of musculoskeletal disorders among adolescents and young adults in 204 countries and territories, 1990-2019[J]. Autoimmun Rev, 2023, 22(8): 103361.
6	Lenti MV, Rossi CM, Melazzini F, et al. Seronegative autoimmune diseases: a challenging diagnosis[J]. Autoimmun Rev, 2022, 21(9): 103143.
7	Tidblad L, Westerlind H, Delcoigne B, et al. Comorbidities and treatment patterns in early rheumatoid arthritis: a nationwide Swedish study[J]. RMD Open, 2022, 8(2): e002700.
8	Ben Mrid R, Bouchmaa N, Ainani H, et al. Anti-rheumatoid drugs advancements: new insights into the molecular treatment of rheumatoid arthritis[J]. Biomed Pharmacother, 2022, 151: 113126.
9	Tan Y, Buch MH. 'Difficult to treat' rheumatoid arthritis: current position and considerations for next steps[J]. RMD Open, 2022, 8(2): e002387.
10	Clough E, Barrett T. The gene expression omnibus database[J]. Methods Mol Biol, 2016, 1418: 93-110.
11	Zhao T, Xie Z, Xi Y, et al. How to model rheumatoid arthritis in animals: from rodents to non-human Primates[J]. Front Immunol, 2022, 13: 887460.
12	Zhang X, Zhang X, Wang X, et al. Efficient delivery of triptolide plus a miR-30-5p inhibitor through the use of near infrared laser responsive or CADY modified MSNs for efficacy in rheumatoid arthritis therapeutics[J]. Front Bioeng Biotechnol, 2020, 8: 170.
13	Rincón-Riveros A, Morales D, Rodríguez JA, et al. Bioinformatic tools for the analysis and prediction of ncRNA interactions[J]. Int J Mol Sci, 2021, 22(21): 11397.
14	Jang S, Kwon EJ, Lee JJ. Rheumatoid arthritis: pathogenic roles of diverse immune cells[J]. Int J Mol Sci, 2022, 23(2): 905.
15	Yu R, Zhang J, Zhuo Y, et al. Identification of diagnostic signatures and immune cell infiltration characteristics in rheumatoid arthritis by integrating bioinformatic analysis and machine-learning strategies[J]. Front Immunol, 2021, 12: 724934.
16	Cheng Q, Chen X, Wu H, et al. Three hematologic/immune system-specific expressed genes are considered as the potential biomarkers for the diagnosis of early rheumatoid arthritis through bioinformatics analysis[J]. J Transl Med, 2021, 19(1): 18.
17	Firestein GS, McInnes IB. Immunopathogenesis of rheumatoid arthritis[J]. Immunity, 2017, 46(2): 183-96.
18	Wu D, Luo Y, Li T, et al. Systemic complications of rheumatoid arthritis: focus on pathogenesis and treatment[J]. Front Immunol, 2022, 13: 1051082.
19	Matsumoto H, Fujita Y, Onizawa M, et al. Increased CEACAM1 expression on peripheral blood neutrophils in patients with rheumatoid arthritis[J]. Front Immunol, 2022, 13: 978435.
20	Li Z, Xu M, Li R, et al. Identification of biomarkers associated with synovitis in rheumatoid arthritis by bioinformatics analyses[J]. Biosci Rep, 2020, 40(9): BSR20201713.
21	Zhou S, Lu H, Xiong M. Identifying immune cell infiltration and effective diagnostic biomarkers in rheumatoid arthritis by bioinformatics analysis[J]. Front Immunol, 2021, 12: 726747.
22	Vantourout P, Hayday A. Six-of-the-best: unique contributions of γδ T cells to immunology[J]. Nat Rev Immunol, 2013, 13: 88-100.
23	Hayday AC. γδ T cell update: adaptate orchestrators of immune surveillance[J]. J Immunol, 2019, 203(2): 311-20.
24	Mo WX, Yin SS, Chen H, et al. Chemotaxis of Vδ2 T cells to the joints contributes to the pathogenesis of rheumatoid arthritis[J]. Ann Rheum Dis, 2017, 76(12): 2075-84.
25	Jacobs MR, Haynes BF. Increase in TCR gamma delta T lymphocytes in synovia from rheumatoid arthritis patients with active synovitis[J]. J Clin Immunol, 1992, 12(2): 130-8.
26	Pascual-García S, Martínez-Peinado P, López-Jaén AB, et al. Analysis of novel immunological biomarkers related to rheumatoid arthritis disease severity[J]. Int J Mol Sci, 2023, 24(15): 12351.
27	Lu J, Wu J, Xia XL, et al. Follicular helper T cells: potential therapeutic targets in rheumatoid arthritis[J]. Cell Mol Life Sci, 2021, 78(12): 5095-106.
28	Zhao S, Wang Y, Liang Y, et al. MicroRNA-126 regulates DNA methylation in CD4⁺ T cells and contributes to systemic lupus erythematosus by targeting DNA methyltransferase 1[J]. Arthritis Rheum, 2011, 63(5): 1376-86.
29	Cutolo M, Campitiello R, Gotelli E, et al. The role of M1/M2 macrophage polarization in rheumatoid arthritis synovitis[J]. Front Immunol, 2022, 13: 867260.
30	Tardito S, Martinelli G, Soldano S, et al. Macrophage M1/M2 polarization and rheumatoid arthritis: a systematic review[J]. Autoimmun Rev, 2019, 18(11): 102397.
31	Ye Q, Luo F, Yan T. Transcription factor KLF4 regulated STAT1 to promote M1 polarization of macrophages in rheumatoid arthritis[J]. Aging: Albany NY, 2022, 14(14): 5669-80.
32	Guo DG, Lv JH, Chen X, et al. Study of miRNA interactome in active rheumatoid arthritis patients reveals key pathogenic roles of dysbiosis in the infection-immune network[J]. Rheumatology: Oxford, 2021, 60(3): 1512-22.
33	Zaiss MM, Wu HJ J, Mauro D, et al. The gut-joint axis in rheumatoid arthritis[J]. Nat Rev Rheumatol, 2021, 17: 224-37.
34	Kondo N, Kuroda T, Kobayashi D. Cytokine networks in the pathogenesis of rheumatoid arthritis[J]. Int J Mol Sci, 2021, 22(20): 10922.
35	Simon LS, Taylor PC, Choy EH, et al. The Jak/STAT pathway: a focus on pain in rheumatoid arthritis[J]. Semin Arthritis Rheum, 2021, 51(1): 278-84.
36	Banerjee S, Biehl A, Gadina M, et al. JAK-STAT signaling as a target for inflammatory and autoimmune diseases: current and future prospects[J]. Drugs, 2017, 77(5): 521-46.
37	Hu L, Liu RJ, Zhang LL. Advance in bone destruction participated by JAK/STAT in rheumatoid arthritis and therapeutic effect of JAK/STAT inhibitors[J]. Int Immunopharmacol, 2022, 111: 109095.
38	Karonitsch T, Saferding V, Kieler M, et al. Amino acids fueling fibroblast-like synoviocyte activation and arthritis by regulating chemokine expression and leukocyte migration[J]. Arthritis Rheumatol, 2024, 76(4): 531-40.
39	Mai YP, Yu XT, Gao T, et al. Autoantigenic peptide and immunomodulator codelivery system for rheumatoid arthritis treatment by reestablishing immune tolerance[J]. ACS Appl Mater Interfaces, 2024: Online ahead of print.

[1]	Weiyi LI, Lu JIANG, Zongxing ZHANG, Dan CHEN, Zhuoma BAO, Li HUANG, Lin YUAN. Qianggu Kangshu Formula attenuates osteoclast differentiation in rheumatoid arthritis by inhibiting the HIF-1α/BNIP3 autophagy signaling pathway [J]. Journal of Southern Medical University, 2025, 45(7): 1389-1396.
[2]	Meimei CHEN, Yang WANG, Huangwei LEI, Fei ZHANG, Ruina HUANG, Zhaoyang YANG. Construction of recognition models for subthreshold depression based on multiple machine learning algorithms and vocal emotional characteristics [J]. Journal of Southern Medical University, 2025, 45(4): 711-717.
[3]	Zhi GAO, Ao WU, Zhongxiang HU, Peiyang SUN. Bioinformatics analysis of oxidative stress and immune infiltration in rheumatoid arthritis [J]. Journal of Southern Medical University, 2025, 45(4): 862-870.
[4]	Zhoufang CAO, Yuan WANG, Mengna WANG, Yue SUN, Feifei LIU. LINC00837/miR-671-5p/SERPINE2 functional axis promotes pathological processes of fibroblast-like synovial cells in rheumatoid arthritis [J]. Journal of Southern Medical University, 2025, 45(2): 371-378.
[5]	Junping ZHAN, Shuo HUANG, Qingliang MENG, Wei FAN, Huimin GU, Jiakang CUI, Huilian WANG. Buyang Huanwu Decoction reduces mitochondrial autophagy in rheumatoid arthritis synovial fibroblasts in hypoxic culture by inhibiting the BNIP3-PI3K/Akt pathway [J]. Journal of Southern Medical University, 2025, 45(1): 35-42.
[6]	Rui YANG, Yuqin SHU, Huijie WEN, Xi CAI, Zhen WANG, Chen ZHANG, Yang XIANG, Hao WU. Pterocarya hupehensis Skan total flavones ameliorate rheumatoid arthritis in rats by suppressing formation of neutrophil extracellular traps [J]. Journal of Southern Medical University, 2024, 44(9): 1645-1652.
[7]	Shan XIANG, Zongxing ZHANG, Lu JIANG, Daozhong LIU, Weiyi LI, Zhuoma BAO, Rui TIAN, Dan CHENG, Lin YUAN. Tujia medicine Toddalia asiatica improves synovial pannus in rats with collagen-induced arthritis through the PI3K/Akt signaling pathway [J]. Journal of Southern Medical University, 2024, 44(8): 1582-1588.
[8]	Caiyu SHEN, Shuai WANG, Ruiying ZHOU, Yuhe WANG, Qin GAO, Xingzhi CHEN, Shu YANG. Prediction of risk of in-hospital death in patients with chronic heart failure complicated by lung infections using interpretable machine learning [J]. Journal of Southern Medical University, 2024, 44(6): 1141-1148.
[9]	Zhiwei ZUO, Qingliang MENG, Jiakang CUI, Kelei GUO, Hua BIAN. An artificial neural network diagnostic model for scleroderma and immune cell infiltration analysis based on mitochondria-associated genes [J]. Journal of Southern Medical University, 2024, 44(5): 920-929.
[10]	LI Yunfei, YANG Jingyi, ZHANG Ying, ZHANG Caixia, WEI Yuxiang, WANG Yiying, WU Ning, SUN Jianfei, WU Zunqiu. The Miao medicine Sidaxue alleviates rheumatoid arthritis in rats possibly by downregulating matrix metalloproteinases [J]. Journal of Southern Medical University, 2024, 44(4): 739-747.
[11]	Wen ZHAO, Hejing RUAN, Siyuan WANG, Yuzhe CHENG, Miao LEI, Jiufa ZHAO, Chuanmiao LIU. Inhibiting Yes-associated protein alleviates CCl₄ liver fibrosis in mice by reducing epithelial mesenchymal transition [J]. Journal of Southern Medical University, 2024, 44(10): 1839-1849.
[12]	Xiaoyin HUANG, Fenglian CHEN, Yu ZHANG, Shujun LIANG. A predictive model for survival outcomes of glioma patients based on multi-parametric, multi-regional MRI radiomics features and clinical features [J]. Journal of Southern Medical University, 2024, 44(10): 2004-2014.
[13]	HE Huishan, GUO Erjia, MENG Wenyi, WANG Yu, WANG Wen, HE Wenle, WU Yuankui, YANG Wei. Predicting cerebral glioma enhancement pattern using a machine learning-based magnetic resonance imaging radiomics model [J]. Journal of Southern Medical University, 2024, 44(1): 194-200.
[14]	WANG Fanfan, LIU Jian, FANG Yanyan, WEN Jianting, HE Mingyu, HAN Qi, LI Xu. Traditional Chinese medicine may reduce the risk of readmission in patients with rheumatoid arthritis complicated with elevated platelet count: a matched cohort study [J]. Journal of Southern Medical University, 2023, 43(9): 1548-1557.
[15]	LUO Xiao, CHENG Yi, WU Cheng, HE Jia. An interpretable machine learning-based prediction model for risk of death for patients with ischemic stroke in intensive care unit [J]. Journal of Southern Medical University, 2023, 43(7): 1241-1247.