Dataset Information

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

ABSTRACT:

Background

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.

Objective

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.

Methods

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.

Results

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.

Conclusions

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.

SUBMITTER: Sun F

PROVIDER: S-EPMC8152926 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundOsteoarthritis (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.ObjectiveThis study was designed to identify key biomarkers and immune infiltration in the synovial tissue of osteoarthritis by bioinformatics analysis.Materials and methodsThe 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.ResultsA 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.ConclusionThe hub genes and the difference in immune infiltration in synovial tissue between osteoarthritis and normal controls might provide new insight for understanding OA development.

Project description:BackgroundRheumatoid arthritis (RA) is a chronic systemic autoimmune disorder characterized by inflammatory cell infiltration, leading to persistent synovitis and joint destruction. The pathogenesis of RA remains unclear. This study aims to explore the immune molecular mechanism of RA through bioinformatics analysis.MethodsFive microarray datasets and a high throughput sequencing dataset were downloaded. CIBERSORT algorithm was performed to evaluate immune cell infiltration in synovial tissues between RA and healthy control (HC). Wilcoxon test and Least Absolute Shrinkage and Selection Operator (LASSO) regression were conducted to identify the significantly different infiltrates of immune cells. Differentially expressed genes (DEGs) were screened by "Batch correction" and "RobustRankAggreg" methods. Functional correlation of DEGs were analyzed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Candidate biomarkers were identified by cytoHubba of Cytoscape, and their diagnostic effectiveness was predicted by Receiver Operator Characteristic Curve (ROC) analysis. The association of the identified biomarkers with infiltrating immune cells was explored using Spearman's rank correlation analysis in R software.ResultsTen significantly different types of immune cells between RA and HC were identified. A total of 202 DEGs were obtained by intersection of DEGs screened by two methods. The function of DEGs were significantly associated with immune cells. Five hub genes (CXCR4, CCL5, CD8A, CD247, and GZMA) were screened by R package "UpSet". CCL5+CXCR4 and GZMA+CD8A were verified to have the capability to diagnose RA and early RA with the most excellent specificity and sensitivity, respectively. The correlation between immune cells and biomarkers showed that CCL5 was positively correlated with M1 macrophages, CXCR4 was positively correlated with memory activated CD4+ T cells and follicular helper T (Tfh) cells, and GZMA was positively correlated with Tfh cells.ConclusionsCCL5, CXCR4, GZMA, and CD8A can be used as diagnostic biomarker for RA. GZMA-Tfh cells, CCL5-M1 macrophages, and CXCR4- memory activated CD4+ T cells/Tfh cells may participate in the occurrence and development of RA, especially GZMA-Tfh cells for the early pathogenesis of RA.

Project description:INTRODUCTION: MicroRNAs (miRNAs), endogenous small noncoding RNAs regulating the activities of target mRNAs and cellular processes, are present in human plasma in a stable form. In this study, we investigated whether miRNAs are also stably present in synovial fluids and whether plasma and synovial fluid miRNAs could be biomarkers of rheumatoid arthritis (RA) and osteoarthritis (OA). METHODS: We measured concentrations of miR-16, miR-132, miR-146a, miR-155 and miR-223 in synovial fluid from patients with RA and OA, and those in plasma from RA, OA and healthy controls (HCs) by quantitative reverse transcription-polymerase chain reaction. Furthermore, miRNAs in the conditioned medium of synovial tissues, monolayer fibroblast-like synoviocytes, and mononuclear cells were examined. Correlations between miRNAs and biomarkers or disease activities of RA were statistically examined. RESULTS: Synovial fluid miRNAs were present and as stable as plasma miRNAs for storage at -20 degrees C and freeze-thawing from -20 degrees C to 4 degrees C. In RA and OA, synovial fluid concentrations of miR-16, miR-132, miR-146a, and miR-223 were significantly lower than their plasma concentrations, and there were no correlation between plasma and synovial fluid miRNAs. Interestingly, synovial tissues, fibroblast-like synoviocytes, and mononuclear cells secreted miRNAs in distinct patterns. The expression patterns of miRNAs in synovial fluid of OA were similar to miRNAs secreted by synovial tissues. Synovial fluid miRNAs of RA were likely to originate from synovial tissues and infiltrating cells. Plasma miR-132 of HC was significantly higher than that of RA or OA with high diagnosability. Synovial fluid concentrations of miR-16, miR-146a miR-155 and miR-223 of RA were significantly higher than those of OA. Plasma miRNAs or ratio of synovial fluid miRNAs to plasma miRNAs, including miR-16 and miR-146a, significantly correlated with tender joint counts and 28-joint Disease Activity Score. CONCLUSIONS: Plasma miRNAs had distinct patterns from synovial fluid miRNAs, which appeared to originate from synovial tissue. Plasma miR-132 well differentiated HCs from patients with RA or OA, while synovial fluid miRNAs differentiated RA and OA. Furthermore, plasma miRNAs correlated with the disease activities of RA. Thus, synovial fluid and plasma miRNAs have potential as diagnostic biomarkers for RA and OA and as a tool for the analysis of their pathogenesis.

Project description:BackgroundOsteoarthritis (OA) and rheumatoid arthritis (RA) were two major joint diseases with similar clinical phenotypes. This study aimed to determine the mechanistic similarities and differences between OA and RA by integrated analysis of multiple gene expression data sets.MethodsMicroarray data sets of OA and RA were obtained from the Gene Expression Omnibus (GEO). By integrating multiple gene data sets, specific differentially expressed genes (DEGs) were identified. The Gene Ontology (GO) functional annotation, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and protein-protein interaction (PPI) network analysis of DEGs were conducted to determine hub genes and pathways. The "Cell Type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT)" algorithm was employed to evaluate the immune infiltration cells (IICs) profiles in OA and RA. Moreover, mouse models of RA and OA were established, and selected hub genes were verified in synovial tissues with quantitative polymerase chain reaction (qPCR).ResultsA total of 1116 DEGs were identified between OA and RA. GO functional enrichment analysis showed that DEGs were enriched in regulation of cell morphogenesis involved in differentiation, positive regulation of neuron differentiation, nuclear speck, RNA polymerase II transcription factor complex, protein serine/threonine kinase activity and proximal promoter sequence-specific DNA binding. KEGG pathway analysis showed that DEGs were enriched in EGFR tyrosine kinase inhibitor resistance, ubiquitin mediated proteolysis, FoxO signaling pathway and TGF-beta signaling pathway. Immune cell infiltration analysis identified 9 IICs with significantly different distributions between OA and RA samples. qPCR results showed that the expression levels of the hub genes (RPS6, RPS14, RPS25, RPL11, RPL27, SNRPE, EEF2 and RPL19) were significantly increased in OA samples compared to their counterparts in RA samples (P < 0.05).ConclusionThis large-scale gene analyses provided new insights for disease-associated genes, molecular mechanisms as well as IICs profiles in OA and RA, which may offer a new direction for distinguishing diagnosis and treatment between OA and RA.

Project description:Heart failure (HF) and osteoarthritis (OA) are medical conditions that can significantly impact daily activities. Evidence has shown that HF and OA may share some pathogenic mechanisms. However, the underlying genomic mechanisms remain unclear. This study aimed to explore the underlying molecular mechanism and identify diagnostic biomarkers for HF and OA. With the cutoff criteria of fold change (FC) > 1.3 and P < .05, 920, 1500, 2195, and 2164 differentially expressed genes (DEGs) were identified in GSE57338, GSE116250, GSE114007, and GSE169077, respectively. After making the intersection of DEGs, we obtained 90 upregulated DEGs and 51 downregulated DEGs in HF datasets and 115 upregulated DEGs and 75 downregulated DEGs in OA datasets. Afterward, we conducted genome ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, protein-protein interaction (PPI) networks, and hub genes screening based on DEGs. Then, 4 common DEGs (fibroblast activation protein alpha [FAP], secreted frizzled-related protein 4 (SFRP4), Thy-1 cell surface antigen (THY1), matrix remodeling associated 5 [MXRA5]) between HF and OA were screened and validated in GSE5406 and GSE113825 datasets, based on which we established the support vector machine (SVM) models. The combined area under the receiver operating characteristic curve (AUC) of THY1, FAP, SFRP4, and MXRA5 in the HF training and test sets reached 0.949 and 0.928. While in the OA training set and test set, the combined AUC of THY1, FAP, SFRP4, and MXRA5 reached 1 and 1, respectively. The analysis of immune cells in HF revealed high levels of dendritic cell (DC), B cells, natural killer T cell (NKT), Type 1 regulatory T cell (Tr1), cytotoxic T cell (Tc), exhausted T cell (Tex), and mucosal-associated invariant T cell (MAIT), while displaying lower levels of monocytes, macrophages, NK, CD4 + T, gamma delta T (γδ T), T helper type 1 (Th1), T helper type 2 (Th2), and effector memory T cell (Tem). Moreover, the 4 common DEGs were positively correlated with DCs and B cells and negatively correlated with γδ T. In OA patients, the abundance of monocyte, macrophage, CD4 + naïve, and natural T regulatory cell (nTreg) was higher, while the infiltration of CD8 + T, γδ T, CD8 + naïve, and MAIT was lower. The expression of THY1 and FAP was significantly correlated with macrophage, CD8 + T, nTreg, and CD8 + naïve. SFRP4 was correlated with monocyte, CD8 + T, γδ T, CD4 + naïve, nTreg, CD8 + naïve and MAIT. MXRA5 was correlated with macrophage, CD8 + T, nTreg and CD8 + naïve. FAP, THY1, MXRA5, and SFRP4 may be diagnostic biomarkers for both HF and OA, and their correlation with immune cell infiltrations suggests shared immune pathogenesis.

Project description:BackgroundDifferential diagnosis among psoriatic arthritis (PsA) and seronegative rheumatoid arthritis (Abneg RA) can be challenging particularly in the clinical setting of peripheral phenotype and autoantibodies seronegativity. The aim of the study was to identify synovial tissue (ST) biomarkers differentially expressed in PsA and Abneg RA and test their predictive value of therapeutic response.MethodsThirty-four PsA patients [12 DMARD naive and 22 non-responder to methotrexate (MTX-IR)] with peripheral joint involvement and 55 Abneg RA (27 DMARD naive and 28 MTX-IR) underwent US-guided ST biopsy and immunohistochemistry (IHC) for CD68+, CD3+, CD20+, CD21+, CD117+, and CD138+ cells. After study entry, each DMARD-naive patient started MTX therapy and was followed in an outpatient setting for at least 6 months to define the achievement of Minimal Disease Activity (PsA) and DAS remission (Abneg RA) status respectively. Each IR-MTX patient was treated according to EULAR recommendations.ResultsAt study entry, IHC analysis revealed that PsA patients had comparable levels of lining and sublining CD68+ and sublining CD21+, CD20+, and CD3+ cells than Abneg RA, despite the therapeutic regimen. Moreover, regardless of the therapeutic scheme, PsA patients showed higher IHC score of CD117+ cells (p = 0.0004 and p = 0.0005 for naive and MTX-IR patients respectively) compared to Abneg RA patients. Conversely, Abneg RA patients showed higher IHC score of CD138+ cells, irrespective to the therapeutic scheme (p = 0.04 and p = 0.002 for naive and MTX-IR patients respectively). Analyzing the response rate to the therapeutic scheme, naive PsA patients reaching MDA status at 6 months follow-up, showed, at the study entry, lower IHC score of CD3+ cells compared to PsA patients not reaching this outcome (p = 0.02); conversely, naive Abneg RA patients reaching DAS remission status at 6 months follow-up, showed, at the study entry, lower IHC score of sublining CD68+ cells compared to Abneg RA patients not reaching this outcome (p < 0.001).ConclusionsCD117+ and CD138+ cells are differentially distributed among PsA and Abneg RA. Histological analysis of ST may help to solve the clinical overlap between the two diseases and provides prognostic data about the therapy success.