Identification of diagnostic genes and drug prediction in metabolic syndrome-associated rheumatoid arthritis by integrated bioinformatics analysis, machine learning, and molecular docking.

Interactions between the immune and metabolic systems may play a crucial role in the pathogenesis of metabolic syndrome-associated rheumatoid arthritis (MetS-RA). The purpose of this study was to discover candidate biomarkers for the diagnosis of RA patients who also had MetS. Three RA datasets and one MetS dataset were obtained from the Gene Expression Omnibus (GEO) database. Differential expression analysis, weighted gene co-expression network analysis (WGCNA), and machine learning algorithms including Least Absolute Shrinkage and Selection Operator (LASSO) regression and Random Forest (RF) were employed to identify hub genes in MetS-RA. Enrichment analysis was used to explore underlying common pathways between MetS and RA. Receiver operating characteristic curves were applied to assess the diagnostic performance of nomogram constructed based on hub genes. Protein-protein interaction, Connectivity Map (CMap) analyses, and molecular docking were utilized to predict the potential small molecule compounds for MetS-RA treatment. qRT-PCR was used to verify the expression of hub genes in fibroblast-like synoviocytes (FLS) of MetS-RA. The effects of small molecule compounds on the function of RA-FLS were evaluated by wound-healing assays and angiogenesis experiments. The CIBERSORT algorithm was used to explore immune cell infiltration in MetS and RA. MetS-RA key genes were mainly enriched in immune cell-related signaling pathways and immune-related processes. Two hub genes (TYK2 and TRAF2) were selected as candidate biomarkers for developing nomogram with ideal diagnostic performance through machine learning and proved to have a high diagnostic value (area under the curve, TYK2, 0.92; TRAF2, 0.90). qRT-PCR results showed that the expression of TYK2 and TRAF2 in MetS-RA-FLS was significantly higher than that in non-MetS-RA-FLS (nMetS-RA-FLS). The combination of CMap analysis and molecular docking predicted camptothecin (CPT) as a potential drug for MetS-RA treatment. In vitro validation, CPT was observed to suppress the cell migration capacity and angiogenesis capacity of MetS-RA-FLS. Immune cell infiltration results revealed immune dysregulation in MetS and RA. Two hub genes were identified in MetS-RA, a nomogram for the diagnosis of RA and MetS was established based on them, and a potential therapeutic small molecule compound for MetS-RA was predicted, which offered a novel research perspective for future serum-based diagnosis and therapeutic intervention of MetS-RA.

Huang Y ，Yue S ，Qiao J ，Dong Y ，Liu Y ，Zhang M ，Zhang C ，Chen C ，Tang Y ，Zheng J ... -  《Frontiers in Immunology》

Machine learning and weighted gene co-expression network analysis identify a three-gene signature to diagnose rheumatoid arthritis.

Rheumatoid arthritis (RA) is a systemic immune-related disease characterized by synovial inflammation and destruction of joint cartilage. The pathogenesis of RA remains unclear, and diagnostic markers with high sensitivity and specificity are needed urgently. This study aims to identify potential biomarkers in the synovium for diagnosing RA and to investigate their association with immune infiltration. We downloaded four datasets containing 51 RA and 36 healthy synovium samples from the Gene Expression Omnibus database. Differentially expressed genes were identified using R. Then, various enrichment analyses were conducted. Subsequently, weighted gene co-expression network analysis (WGCNA), random forest (RF), support vector machine-recursive feature elimination (SVM-RFE), and least absolute shrinkage and selection operator (LASSO) were used to identify the hub genes for RA diagnosis. Receiver operating characteristic curves and nomogram models were used to validate the specificity and sensitivity of hub genes. Additionally, we analyzed the infiltration levels of 28 immune cells in the expression profile and their relationship with the hub genes using single-sample gene set enrichment analysis. Three hub genes, namely, ribonucleotide reductase regulatory subunit M2 (RRM2), DLG-associated protein 5 (DLGAP5), and kinesin family member 11 (KIF11), were identified through WGCNA, LASSO, SVM-RFE, and RF algorithms. These hub genes correlated strongly with T cells, natural killer cells, and macrophage cells as indicated by immune cell infiltration analysis. RRM2, DLGAP5, and KIF11 could serve as potential diagnostic indicators and treatment targets for RA. The infiltration of immune cells offers additional insights into the underlying mechanisms involved in the progression of RA.

Wu YK ，Liu CD ，Liu C ，Wu J ，Xie ZG ... -  《Frontiers in Immunology》

Identification of Immune-Associated Genes in Diagnosing Aortic Valve Calcification With Metabolic Syndrome by Integrated Bioinformatics Analysis and Machine Learning.

Immune system dysregulation plays a critical role in aortic valve calcification (AVC) and metabolic syndrome (MS) pathogenesis. The study aimed to identify pivotal diagnostic candidate genes for AVC patients with MS. We obtained three AVC and one MS dataset from the gene expression omnibus (GEO) database. Identification of differentially expressed genes (DEGs) and module gene via Limma and weighted gene co-expression network analysis (WGCNA), functional enrichment analysis, protein-protein interaction (PPI) network construction, and machine learning algorithms (least absolute shrinkage and selection operator (LASSO) regression and random forest) were used to identify candidate immune-associated hub genes for diagnosing AVC with MS. To assess the diagnostic value, the nomogram and receiver operating characteristic (ROC) curve were developed. Finally, immune cell infiltration was created to investigate immune cell dysregulation in AVC. The merged AVC dataset included 587 DEGs, and 1,438 module genes were screened out in MS. MS DEGs were primarily enriched in immune regulation. The intersection of DEGs for AVC and module genes for MS was 50, which were mainly enriched in the immune system as well. Following the development of the PPI network, 26 node genes were filtered, and five candidate hub genes were chosen for nomogram building and diagnostic value evaluation after machine learning. The nomogram and all five candidate hub genes had high diagnostic values (area under the curve from 0.732 to 0.982). Various dysregulated immune cells were observed as well. Five immune-associated candidate hub genes (BEX2, SPRY2, CXCL16, ITGAL, and MORF4L2) were identified, and the nomogram was constructed for AVC with MS diagnosis. Our study could provide potential peripheral blood diagnostic candidate genes for AVC in MS patients.

Zhou Y ，Shi W ，Zhao D ，Xiao S ，Wang K ，Wang J ... -  《Frontiers in Immunology》

Identifying hub circadian rhythm biomarkers and immune cell infiltration in rheumatoid arthritis.

Rheumatoid arthritis (RA) is a chronic systemic autoimmune disease with symptoms characterized by typical circadian rhythmic changes. This study aimed to identify the hub circadian rhythm genes (CRGs) in RA and explore their association with immune cell infiltration and pathogenesis of RA. The differentially expressed CRGs (DECRGs) between RA and normal control samples were screened from Datasets GSE12021 and GSE55235. Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, and Gene Set Enrichment Analysis were used to explore the potential functional mechanisms of DECRGs in RA. Weighted Gene Co-expression Network Analysis and Least Absolute Shrinkage and Selection Operator regression analysis were performed to identify hub CRGs of RA. CIBERSORT was conducted to compare the infiltration level of immune cells in RA and control synovial tissue and their relationship with hub genes. In addition, the diagnostic value of hub biomarkers was evaluated by the area under the receiver operator characteristic curve. Further, a nomogram prediction model was constructed and its significance for clinical decision-making was evaluated. The green module was identified as the hub module associated with RA. Four hub CRGs (EGR1, FOSL2, GADD45B, and NFIL3) were identified and showed that they had the highest specificity and sensitivity for RA diagnosis, respectively. The expression levels and diagnostic values of these genes were externally validated in the dataset GSE55457. A nomogram prediction model based on the four hub CRGs was constructed and proved to have a certain clinical decision value. Additionally, the correlation analysis of immune cells with hub genes showed that all hub genes were significantly positively correlated with activated mast cells, resting memory CD4+ T cells, and monocytes. Whereas, all hub genes were negatively correlated with plasma cells, CD8+ T cells, and activated memory CD4+ T cells. Meanwhile, FOSL2 and GADD45B were negatively correlated with Tfh cells. Four hub CRGs were identified and showed excellent diagnostic value for RA. These genes may be involved in the pathological process of RA by disrupting the rhythmic oscillations of cytokines through immune-related pathways and could be considered molecular targets for future chronotherapy against RA.

Wen P ，Ma T ，Zhang B ，Hao L ，Wang Y ，Guo J ，Song W ，Wang J ，Zhang Y ... -  《Frontiers in Immunology》

Identification of immune-related genes in diagnosing atherosclerosis with rheumatoid arthritis through bioinformatics analysis and machine learning.

Increasing evidence has proven that rheumatoid arthritis (RA) can aggravate atherosclerosis (AS), and we aimed to explore potential diagnostic genes for patients with AS and RA. We obtained the data from public databases, including Gene Expression Omnibus (GEO) and STRING, and obtained the differentially expressed genes (DEGs) and module genes with Limma and weighted gene co-expression network analysis (WGCNA). Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analysis, the protein-protein interaction (PPI) network, and machine learning algorithms [least absolute shrinkage and selection operator (LASSO) regression and random forest] were performed to explore the immune-related hub genes. We used a nomogram and receiver operating characteristic (ROC) curve to assess the diagnostic efficacy, which has been validated with GSE55235 and GSE73754. Finally, immune infiltration was developed in AS. The AS dataset included 5,322 DEGs, while there were 1,439 DEGs and 206 module genes in RA. The intersection of DEGs for AS and crucial genes for RA was 53, which were involved in immunity. After the PPI network and machine learning construction, six hub genes were used for the construction of a nomogram and for diagnostic efficacy assessment, which showed great diagnostic value (area under the curve from 0.723 to 1). Immune infiltration also revealed the disorder of immunocytes. Six immune-related hub genes (NFIL3, EED, GRK2, MAP3K11, RMI1, and TPST1) were recognized, and the nomogram was developed for AS with RA diagnosis.

Liu F ，Huang Y ，Liu F ，Wang H ... -  《Frontiers in Immunology》