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》

Identification of Disease-Specific Hub Biomarkers and Immune Infiltration in Osteoarthritis and Rheumatoid Arthritis Synovial Tissues by Bioinformatics Analysis.

Osteoarthritis (OA) and rheumatoid arthritis (RA) are well-known cause of joint disability. Although they have shown the analogous clinical features involving chronic synovitis that progresses to cartilage and bone destruction, the pathogenesis that initiates and perpetuates synovial lesions between RA and OA remains elusive. This study is aimed at identifying disease-specific hub genes, exploring immune cell infiltration, and elucidating the underlying mechanisms associated with RA and OA synovial lesion. Gene expression profiles (GSE55235, GSE55457, GSE55584, and GSE12021) were selected from Gene Expression Omnibus for analysis. Differentially expressed genes (DEGs) were identified by the "LIMMA" package in Bioconductor. The DEGs were identified by Gene Ontology (GO) and KEGG pathway analysis. A protein-protein interaction network was constructed to identify candidate hub genes by using STRING and Cytoscape. Hub genes were identified by validating from GSE12021. Furthermore, we employed the CIBERSORT website to assess immune cell infiltration between OA and RA. Finally, we explored the correlation between the levels of hub genes and relative proportion of immune cells in OA and RA. We identified 68 DEGs which were mainly enriched in immune response and chemokine signaling pathway. Six hub genes with a cutoff of AUC > 0.80 by ROC analysis and relative expression of P < 0.05 were identified successfully. Compared with OA, the RA synovial tissues consisted of a higher proportion of 7 immune cells, whereas 4 immune cells were found in relatively lower proportion (P < 0.05). In addition, the levels of 6 hub genes were closely associated with relative proportion of 11 immune cells in OA and RA. We used bioinformatics analysis to identify hub genes and explored immune cell infiltration of immune microenvironment in synovial tissues. Our results should offer insights into the underlying molecular mechanisms of synovial lesion and provide potential target for immune-based therapies of OA and RA.

Sun F ，Zhou JL ，Peng PJ ，Qiu C ，Cao JR ，Peng H ... -  《-》

Identification of key biomarkers and immune infiltration in the synovial tissue of osteoarthritis by bioinformatics analysis.

Osteoarthritis (OA) is the most common chronic degenerative joint disease and is mainly characterized by cartilage degeneration, subcartilage bone hyperplasia, osteophyte formation and joint space stenosis. Recent studies showed that synovitis might also be an important pathological change of OA. However, the molecular mechanisms of synovitis in OA are still not well understood. This study was designed to identify key biomarkers and immune infiltration in the synovial tissue of osteoarthritis by bioinformatics analysis. The gene expression profiles of GSE12021, GSE55235 and GSE55457 were downloaded from the GEO database. The differentially expressed genes (DEGs) were identified by the LIMMA package in Bioconductor, and functional enrichment analyses were performed. A protein-protein interaction network (PPI) was constructed, and module analysis was performed using STRING and Cytoscape. The CIBERSORT algorithm was used to analyze the immune infiltration of synovial tissue between OA and normal controls. A total of 106 differentially expressed genes, including 68 downregulated genes and 38 upregulated genes, were detected. The PPI network was assessed, and the most significant module containing 14 hub genes was identified. Gene Ontology analysis revealed that the hub genes were significantly enriched in immune cell chemotaxis and cytokine activity. KEGG pathway analysis showed that the hub genes were significantly enriched in the rheumatoid arthritis signaling pathway, IL-17 signaling pathway and cytokine-cytokine receptor interaction signaling pathway. The immune infiltration profiles varied significantly between osteoarthritis and normal controls. Compared with normal tissue, OA synovial tissue contained a higher proportion of memory B cells, naive CD4+ T cells, regulatory T cells, resting dendritic cells and resting mast cells, while naive CD4+ T cells, activated NK cells, activated mast cells and eosinophils contributed to a relatively lower portion (P > 0.05). Finally, the expression levels of 11 hub genes were confirmed by RT-PCR. The hub genes and the difference in immune infiltration in synovial tissue between osteoarthritis and normal controls might provide new insight for understanding OA development.

Cai W ，Li H ，Zhang Y ，Han G ... -  《PeerJ》

Identifying the Hub Genes and Immune Cell Infiltration in Synovial Tissue between Osteoarthritic and Rheumatoid Arthritic Patients by Bioinformatic Approach.

Osteoarthritis (OA) and rheumatoid arthritis (RA) are two common diseases that result in limb disability and a decrease in quality of life. The major symptoms of OA and RA are pain, swelling, stiffness, and malformation of joints, and each disease also has unique characteristics. To compare the pathological mechanisms of OA and RA via weighted correlation network analysis (WGCNA) and immune infiltration analysis and find potential diagnostic and pharmaceutical targets for the treatment of OA and RA. The gene expression profiles of ten OA and ten RA synovial tissue samples were downloaded from the Gene Expression Omnibus (GEO) database (GSE55235). After obtaining differentially expressed genes (DEGs) via GEO2R, WGCNA was conducted using an R package, and modules and genes that were highly correlated with OA and RA were identified. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment, and protein-protein interaction (PPI) network analyses were also conducted. Hub genes were identified using the Search Tool for the Retrieval of Interacting Genes (STRING) and Cytoscape software. Immune infiltration analysis was conducted using the Perl program and CIBERSORT software. Two hundred ninety-nine DEGs, 24 modules, 16 GO enrichment terms, 6 KEGG pathway enrichment terms, 10 hub genes (CXCL9, CXCL10, CXCR4, CD27, CD69, CD3D, IL7R, STAT1, RGS1, and ISG20), and 8 kinds of different infiltrating immune cells (plasma cells, CD8 T cells, activated memory CD4 T cells, T helper follicular cells, M1 macrophages, Tregs, resting mast cells, and neutrophils) were found to be involved in the different pathological mechanisms of OA and RA. Inflammation-associated genes were the top differentially expressed hub genes between OA and RA, and their expression was downregulated in OA. Genes associated with lipid metabolism may have upregulated expression in OA. In addition, immune cells that participate in the adaptive immune response play an important role in RA. OA mainly involves immune cells that are associated with the innate immune response.

Wang J ，Fan Q ，Yu T ，Zhang Y ... -  《-》

Identification of PPARG as key gene to link coronary atherosclerosis disease and rheumatoid arthritis via microarray data analysis.

Inflammation is the common pathogenesis of coronary atherosclerosis disease (CAD) and rheumatoid arthritis (RA). Although it is established that RA increases the risk of CAD, the underlining mechanism remained indefinite. This study seeks to explore the molecular mechanisms of RA linked CAD and identify potential target gene for early prediction of CAD in RA patients. The study utilized five raw datasets: GSE55235, GSE55457, GSE12021 for RA patients, and GSE42148 and GSE20680 for CAD patients. Gene Set Enrichment Analysis (GSEA) was used to investigate common signaling pathways associated with RA and CAD. Then, weighted gene co-expression network analysis (WGCNA) was performed on RA and CAD training datasets to identify gene modules related to single-sample GSEA (ssGSEA) scores. Overlapping module genes and differentially expressed genes (DEGs) were considered as co-susceptible genes for both diseases. Three hub genes were screened using a protein-protein interaction (PPI) network analysis via Cytoscape plug-ins. The signaling pathways, immune infiltration, and transcription factors associated with these hub genes were analyzed to explore the underlying mechanism connecting both diseases. Immunohistochemistry and qRT-PCR were conducted to validate the expression of the key candidate gene, PPARG, in macrophages of synovial tissue and arterial walls from RA and CAD patients. The study found that Fc-gamma receptor-mediated endocytosis is a common signaling pathway for both RA and CAD. A total of 25 genes were screened by WGCNA and DEGs, which are involved in inflammation-related ligand-receptor interactions, cytoskeleton, and endocytosis signaling pathways. The principal component analysis(PCA) and support vector machine (SVM) and receiver-operator characteristic (ROC) analysis demonstrate that 25 DEGs can effectively distinguish RA and CAD groups from normal groups. Three hub genes TUBB2A, FKBP5, and PPARG were further identified by the Cytoscape software. Both FKBP5 and PPARG were downregulated in synovial tissue of RA and upregulated in the peripheral blood of CAD patients and differential mRNAexpreesion between normal and disease groups in both diseases were validated by qRT-PCR.Association of PPARG with monocyte was demonstrated across both training and validation datasets in CAD. PPARG expression is observed in control synovial epithelial cells and foamy macrophages of arterial walls, but was decreased in synovial epithelium of RA patients. Its expression in foamy macrophages of atherosclerotic vascular walls exhibits a positive correlation (r = 0.6276, p = 0.0002) with CD68. Our findings suggest that PPARG may serve as a potentially predictive marker for CAD in RA patients, which provides new insights into the molecular mechanism underling RA linked CAD.

Zhang Z ，Chen Y ，Fu X ，Chen L ，Wang J ，Zheng Q ，Zhang S ，Zhu X ... -  《PLoS One》