Identification of Critical Biomarkers and Immune Infiltration in Rheumatoid Arthritis Based on WGCNA and LASSO Algorithm.

Rheumatoid arthritis(RA) is the most common inflammatory arthritis, and a significant cause of morbidity and mortality. RA patients' synovial inflammation contains a variety of genes and signalling pathways that are poorly understood. It was the goal of this research to discover the major biomarkers related to the course of RA and how they connect to immune cell infiltration. The Gene Expression Omnibus was used to download gene microarray data. Differential expression analysis, weighted gene co-expression network analysis (WGCNA), and least absolute shrinkage and selection operator (LASSO) regression were used to identify hub markers for RA. Single-sample GSEA was used to examine the infiltration levels of 28 immune cells and their connection to hub gene markers. The hub genes' expression in RA-HFLS and HFLS cells was verified by RT-PCR. The CCK-8 assay was applied to determine the roles of hub genes in RA. In this study, we identified 21 differentially expressed genes (DEGs) in RA. WGCNA yielded two co-expression modules, one of which exhibited the strongest connection with RA. Using a combination of differential genes, a total of 6 intersecting genes was discovered. Six hub genes were identified as possible biomarkers for RA after a lasso analysis was performed on the data. Three hub genes, CKS2, CSTA, and LY96, were found to have high diagnostic value using ROC curve analysis. They were shown to be closely related to the concentrations of several immune cells. RT-PCR confirmed that the expressions of CKS2, CSTA and LY96 were distinctly upregulated in RA-HFLS cells compared with HFLS cells. More importantly, knockdown of CKS2 suppressed the proliferation of RA-HFLS cells. Overall, to help diagnose and treat RA, it's expected that CKS2, CSTA, and LY96 will be available, and the aforementioned infiltration of immune cells may have a significant impact on the onset and progression of the disease.

Jiang F ，Zhou H ，Shen H 《Frontiers in Immunology》

Screening and identification of potential hub genes and immune cell infiltration in the synovial tissue of rheumatoid arthritis by bioinformatic approach.

Rheumatoid arthritis (RA) is an autoimmune disease that affects individuals of all ages. The basic pathological manifestations are synovial inflammation, pannus formation, and erosion of articular cartilage, bone destruction will eventually lead to joint deformities and loss of function. However, the specific molecular mechanisms of synovitis tissue in RA are still unclear. Therefore, this study aimed to screen and explore the potential hub genes and immune cell infiltration in RA. Three microarray datasets (GSE12021, GSE55457, and GSE55235), from the Gene Expression Omnibus (GEO) database, have been analyzed to explore the potential hub genes and immune cell infiltration in RA. First, the LIMMA package was used to screen the differentially expression genes (DEGs) after removing the batch effect. Then the clusterProfiler package was used to perform functional enrichment analyses. Second, through weighted coexpression network analysis (WGCNA), the key module was identified in the coexpression network of the gene set. Third, the protein-protein interaction (PPI) network was constructed through STRING website and the module analysis was performed using Cytoscape software. Fourth, the CIBERSORT and ssGSEA algorithm were used to analyze the immune status of RA and healthy synovial tissue, and the associations between immune cell infiltration and RA-related diagnostic biomarkers were evaluated. Fifth, we used the quantitative reverse transcription-polymerase chain reaction (qRT-PCR) to validate the expression levels of the hub genes, and ROC curve analysis of hub genes for discriminating between RA and healthy tissue. Finally, the gene-drug interaction network was constructed using DrugCentral database, and identification of drug molecules based on hub genes using the Drug Signature Database (DSigDB) by Enrichr. A total of 679 DEGs were identified, containing 270 downregulated genes and 409 upregulated genes. DEGs were primarily enriched in immune response and chemokine signaling pathways, according to functional enrichment analysis of DEGs. WGCNA explored the co-expression network of the gene set and identified key modules, the blue module was selected as the key module associated with RA. Seven hub genes are identified when PPI network and WGCNA core modules are intersected. Immune infiltration analysis using CIBERSORT and ssGSEA algorithms revealed that multiple types of immune infiltration were found to be upregulated in RA tissue compared to normal tissue. Furthermore, the levels of 7 hub genes were closely related to the relative proportions of multiple immune cells in RA. The results of the qRT-PCR demonstrated that the relative expression levels of 6 hub genes (CD27, LCK, CD2, GZMB, IL7R, and IL2RG) were up-regulated in RA synovial tissue, compared with normal tissue. Simultaneously, ROC curves indicated that the above 6 hub genes had strong biomarker potential for RA (AUC >0.8). Through bioinformatics analysis and qRT-PCR experiment, our study ultimately discovered 6 hub genes (CD27, LCK, CD2, GZMB, IL7R, and IL2RG) that closely related to RA. These findings may provide valuable direction for future RA clinical diagnosis, treatment, and associated research.

Feng ZW ，Tang YC ，Sheng XY ，Wang SH ，Wang YB ，Liu ZC ，Liu JM ，Geng B ，Xia YY ... -  《Heliyon》

CKS2 and S100A12: Two Novel Diagnostic Biomarkers for Rheumatoid Arthritis.

Rheumatoid arthritis (RA) is a chronic systematicness autoimmunity disease with joint inflammation. RA etiology is still unknown. Early and exact diagnosing is still hard to reach. In the paper, we purposed to discover novel diagnosis biological marker for RA. Two open, usable gene expression profiles of human RA as well as controlled specimens (dataset GSE17755 as well as GSE93272) were downloaded from the GEO database. Differentially expressed genes (DEGs) were screened between 331 RA and 88 control samples. Functional enrichment analysis was applied to explore the possible function of DEGs. Expression levels as well as diagnosis values of biological marker in RA were further verified in our cohort by the use of RT-PCR and ROC assays. We identified 13 DEGs between RA samples and control samples. 13 DEGs were remarkably abundant in NF-kappa B signal pathway. Among the 13 DEGs, CKS2, S100A12, LY96, and ANXA3 exhibited a strong diagnostic ability in screening RA specimens from normal specimens using all AUC > 0.8. Moreover, we confirmed that the expression of CKS2 and S100A12 was distinctly upregulated in RA specimens contrasted to normal specimens. Overall, serum CKS2 and S100A12 could be used as novel diagnosis biological markers for RA patients.

Wang ZW ，Zhang LJ ，Zhuang Y ，Lv ZF ，Tan ZM ... -  《-》

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 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》