Dataset Information

Identification and validation of hub genes of synovial tissue for patients with osteoarthritis and rheumatoid arthritis.

ABSTRACT:

Background

Osteoarthritis (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.

Methods

Microarray 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).

Results

A 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).

Conclusion

This 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.

SUBMITTER: Ge Y

PROVIDER: S-EPMC8480049 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundRheumatoid 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.MethodsThree 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.ResultsA 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).ConclusionsThrough 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.

Project description:BackgroundOsteoarthritis (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.ObjectiveThis 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.MethodsGene 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.ResultsWe 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.ConclusionsWe 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.

Project description:Background: Rheumatoid arthritis (RA) is a common autoimmune disease characterized by progressive, destructive polyarthritis. However, the cause and underlying molecular events of RA are not clear. Here, we applied integrated bioinformatics to identify tissue-specific expressed hub genes involved in RA and reveal potential targeted drugs. Methods: Three expression profiles of human microarray datasets involving fibroblast-like synoviocytes (FLS) were downloaded from the Gene Expression Omnibus (GEO) database, the differentially expressed mRNAs (DEGs), miRNAs (DEMs), and lncRNAs (DELs) between normal and RA synovial samples were screened using GEO2R tool. BioGPS was used to identified tissue-specific expressed genes. Functional and pathway enrichment analyses were performed for common DEGs using the DAVID database, and the protein-protein interaction (PPI) network of common DEGs was constructed to recognize hub genes by the STRING database. Based on receiver operating characteristic (ROC) curve, we further investigated the prognostic values of tissue-specific expressed hub genes in RA patients. Connectivity Map (CMap) was run to identify novel anti-RA potential drugs. The DEM-DEG pairs and ceRNA network containing key DEMs were established by Cytoscape. Results: We obtain a total of 418 DEGs, 23 DEMs and 49 DELs. 64 DEGs were verified as tissue-specific expressed genes, most derive from the hematologic/immune system (20/64, 31.25%). GO term and KEGG pathway enrichment analysis showed that DEGs focused primarily on immune-related biological process and NF-κB pathway. 10 hub genes were generated via using MCODE plugin. Among them, SPAG5, CUX2, and THEMIS2 were identified as tissue-specific expressed hub genes, these 3 tissue-specific expressed hub genes have superior diagnostic value in the RA samples compared with osteoarthritis (OA) samples. 5 compounds (troleandomycin, levodopa, trichostatin A, LY-294002, and levamisole) rank among the top five in connectivity score. In addition, 5 miRNAs were identified to be key DEMs, the lncRNA-miRNA-mRNA network with five key DEMs was formed. The networks containing tissue-specific expressed hub genes are as follows: ARAP1-AS2/miR-20b-3p/TRIM3, ARAP1-AS2/miR-30c-3p/FRZB. Conclusion: This study indicates that screening for identify tissue-specific expressed hub genes and ceRNA network in RA using integrated bioinformatics analyses could help us understand the mechanism of development of RA. Besides, SPAG5 and THEMIS2 might be candidate biomarkers for diagnosis of RA. LY-294002, trichostatin A, and troleandomycin may be potential drugs for RA.

Project description:BACKGROUND:Rheumatoid arthritis and osteoarthritis are two common musculoskeletal disorders that affect the joints. Despite high prevalence rates, etiological factors involved in these disorders remain largely unknown. Dissecting the molecular aspects of these disorders will significantly contribute to improving their diagnosis and clinical management. In order to identify proteins that are differentially expressed between these two conditions, a quantitative proteomic profiling of synovial fluid obtained from rheumatoid arthritis and osteoarthritis patients was carried out by using iTRAQ labeling followed by high resolution mass spectrometry analysis. RESULTS:We have identified 575 proteins out of which 135 proteins were found to be differentially expressed by ?3-fold in the synovial fluid of rheumatoid arthritis and osteoarthritis patients. Proteins not previously reported to be associated with rheumatoid arthritis including, coronin-1A (CORO1A), fibrinogen like-2 (FGL2), and macrophage capping protein (CAPG) were found to be upregulated in rheumatoid arthritis. Proteins such as CD5 molecule-like protein (CD5L), soluble scavenger receptor cysteine-rich domain-containing protein (SSC5D), and TTK protein kinase (TTK) were found to be upregulated in the synovial fluid of osteoarthritis patients. We confirmed the upregulation of CAPG in rheumatoid arthritis synovial fluid by multiple reaction monitoring assay as well as by Western blot. Pathway analysis of differentially expressed proteins revealed a significant enrichment of genes involved in glycolytic pathway in rheumatoid arthritis. CONCLUSIONS:We report here the largest identification of proteins from the synovial fluid of rheumatoid arthritis and osteoarthritis patients using a quantitative proteomics approach. The novel proteins identified from our study needs to be explored further for their role in the disease pathogenesis of rheumatoid arthritis and osteoarthritis.Sartaj Ahmad and Raja Sekhar Nirujogi contributed equally to this article.

Project description:Rheumatoid arthritis (RA) is an autoimmune disease characterized by systemic inflammation, especially synovitis, leading to joint damage. It is important to explore potential biomarkers and therapeutic targets to improve the clinical treatment of RA. However, the potential underlying mechanisms of action of available treatments for RA have not yet been fully elucidated. The present study investigated the potential biomarkers of RA and identified specific targets for therapeutic intervention. A comprehensive analysis was performed using mRNA files downloaded from the Gene Expression Omnibus. Differences in gene expression were analyzed and compared between the normal and RA groups. In addition, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed on differentially expressed genes (DEGs). A protein-protein interaction network, Molecular Complex Detection and cytoHubba network were evaluated to identify hub genes. Finally, using an experimental RA rat model induced by Freund's complete adjuvant (FCA), the expression of potential biomarkers or target genes in RA were verified through reverse transcription-quantitative PCR. The results of the mRNA dataset processing revealed 195 DEGs in patients with RA when compared with the healthy controls. Moreover, 10 hub genes were identified in patients with RA and four candidate mRNAs were identified, as follows: Discs large homolog-associated protein 5 (DLGAP5), kinesin family member 20A (KIF20A), maternal embryonic leucine zipper kinase (MELK) and nuclear division cycle 80 (NDC80). Finally, the bioinformatics analysis results were validated by quantifying the expression of the DLGAP5, KIF20A, MELK and NDC80 genes in the FCA-induced experimental RA rat model. The findings of the present study suggested that the treatment of RA may be successful through the inhibition of DLGAP5, KIF20A, MELK and NDC80 expression. Therefore, the targeting of these genes may result in more effective treatments for patients with RA.