Identification of key genes in rheumatoid arthritis and osteoarthritis based on bioinformatics analysis.
Ontology highlight
ABSTRACT:
SUBMITTER: Zhu N
PROVIDER: S-EPMC6392928 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BACKGROUND:Rheumatoid arthritis (RA) is a chronic auto-inflammatory disorder of joints. The present study aimed to identify the key genes in RA for better understanding the underlying mechanisms of RA. METHODS:The integrated analysis of expression profiling was conducted to identify differentially expressed genes (DEGs) in RA. Moreover, functional annotation, protein-protein interaction (PPI) network and transcription factor (TF) regulatory network construction were applied for exploring the potential biological roles of DEGs in RA. In addition, the expression level of identified candidate DEGs was preliminarily detected in peripheral blood cells of RA patients in the GSE17755 dataset. Quantitative real-time polymerase chain reaction (qRT-PCR) was conducted to validate the expression levels of identified DEGs in RA. RESULTS:A total of 378 DEGs, including 202 up- and 176 down-regulated genes, were identified in synovial tissues of RA patients compared with healthy controls. DEGs were significantly enriched in axon guidance, RNA transport and MAPK signaling pathway. RBFOX2, LCK and SERBP1 were the hub proteins in the PPI network. In the TF-target gene network, RBFOX2, POU6F1, WIPF1 and PFKFB3 had the high connectivity with TFs. The expression status of 11 candidate DEGs was detected in GSE17755, the expression levels of MAT2A and NSA2 were significantly down-regulated and CD47 had the up-regulated tendency in peripheral blood cells of patients with RA compared with healthy individuals. qRT-PCR results of MAT2A, NSA2, CD47 were compatible with our bioinformatics analyses. DISCUSSION:Our study might provide valuable information for exploring the pathogenesis mechanism of RA and identifying the potential biomarkers for RA diagnosis.
| S-EPMC5356478 | biostudies-literature
Project description:We aimed to identify key genes associated with rheumatoid arthritis (RA).The microarray datasets of GSE1919, GSE12021, and GSE21959 (35 RA samples and 32 normal controls) were downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) in RA samples were identified using the t test in limma package. Functional enrichment analysis was performed using clusterProfiler package. A protein-protein interaction (PPI) network of selected DEGs was constructed based on the Human Protein Reference Database. Active modules were explored using the jActiveModules plug-in in the Cytoscape Network Modeling package.In total, 537 DEGs in RA samples were identified, including 241 upregulated and 296 downregulated genes. A total of 24,451 PPI pairs were collected, and 5 active modules were screened. Furthermore, 19 submodules were acquired from the 5 active modules. Discs large homolog 1 (DLG1) and related DEGs such as guanylate cyclase 1, soluble, alpha 2 (GUCY1A2), N-methyl d-aspartate receptor 2A subunit (GRIN2A), and potassium voltage-gated channel member 1 (KCNA1) were identified in 8 submodules. Plasminogen (PLG) and related DEGs such as chemokine (C-X-C motif) ligand 2 (CXCL2), laminin, alpha 3 (LAMA3), complement component 7 (C7), and coagulation factor X (F10) were identified in 4 submodules.Our results indicate that DLG1, GUCY1A2, GRIN2A, KCNA1, PLG, CXCL2, LAMA3, C7, and F10 may play key roles in the progression of RA and may serve as putative therapeutic targets for treating RA.
| S-EPMC5626145 | biostudies-literature
Project description:MicroRNAs (miRNAs) act as epigenetic markers and regulate the expression of their target genes, including those characterized as regulators in autoimmune diseases. Rheumatoid arthritis (RA) is one of the most common autoimmune diseases. The potential roles of miRNA-regulated genes in RA pathogenesis have greatly aroused the interest of clinicians and researchers in recent years. In the current study, RA-related miRNAs records were obtained from PubMed through conditional literature retrieval. After analyzing the selected records, miRNA targeted genes were predicted. We identified 14 RA-associated miRNAs, and their sub-analysis in 5 microarray or RNA sequencing (RNA-seq) datasets was performed. The microarray and RNA-seq data of RA were also downloaded from NCBI Gene Expression Omnibus (GEO) and Sequence Read Archive (SRA), analyzed, and annotated. Using a bioinformatics approach, we identified a series of differentially expressed genes (DEGs) by comparing studies on RA and the controls. The RA-related gene expression profile was thus obtained and the expression of miRNA-regulated genes was analyzed. After functional annotation analysis, we found GO molecular function (MF) terms significantly enriched in calcium ion binding (GO: 0005509). Moreover, some novel dysregulated target genes were identified in RA through integrated analysis of miRNA/mRNA expression. The result revealed that the expression of a number of genes, including ROR2, ABI3BP, SMOC2, etc., was not only affected by dysregulated miRNAs, but also altered in RA. Our findings indicate that there is a close association between negatively correlated mRNA/miRNA pairs and RA. These findings may be applied to identify genetic markers for RA diagnosis and treatment in the future.
| S-EPMC4567271 | biostudies-literature
Project description:The current study aimed to identify gene signatures during rheumatoid arthritis (RA) and osteoarthritis (OA), and used these to elucidate the underlying modular mechanisms. Using the Gene Expression Omnibus database, the present study obtained the GSE7669 mRNA expression microarray data from RA and OA synovial fibroblasts (n=6 each). The differentially expressed genes (DEGs) in RA synovial samples compared with OA samples were identified using the Linear Models for Microarray Analysis package. The Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed using the Database for Annotation Visualization and Integrated Discovery. A protein‑protein interaction network was constructed and the modules were further analyzed using the Molecular Complex Detection plugin of Cytoscape. A total of 181 DEGs were identified by comparing RA and OA synovial samples (96 up‑ and 85 downregulated genes). The significant DEGs in module 1, including collagen, type I, α 1 (COL1A1), COL3A1, COL4A1 and COL11A1, were predominantly enriched in the extracellular matrix (ECM)‑receptor interaction and focal adhesion pathways. Additionally, significant DEGs in module 2, including radical S‑adenosyl methionine domain containing 2 (RSAD2), 2'‑5'‑oligoadenylate synthetase 2 (OAS2), myxovirus (influenza virus) resistance 1 (MX1) and ISG15 ubiquitin‑like modifier (ISG15), were predominantly associated with immune function pathways. In conclusion, the present study indicated that RSAD2, OAS2, MX1 and ISG15 may be notable gene signatures in RA development via regulation of the immune response. COL3A1, COL4A1, COL1A1 and COL11A1 may be important gene signatures in OA development via involvement in the pathways of ECM-receptor interactions and focal adhesions.
| S-EPMC4940106 | biostudies-literature
Project description:Sex differences have been suggested to play critical roles in the pathophysiology of osteoarthritis (OA), resulting in sex-specific prevalence and incidence. However, their roles in the development of OA remain largely unknown. The aim of this study was to screen out key genes and pathways mediating biological differences between OA females after menopause and OA males. First, the gene expression data of GSE36700 and GSE55457 were downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) between sexes were identified using R software, respectively. The overlapping DEGs were obtained. Then, protein-protein interactive (PPI) network was constructed to further analyze interactions between the overlapping DEGs. Finally, enrichment analyses were separately performed using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes tools. In our results, a total of 278 overlapping DEGs were identified between OA females after menopause and OA males, including 219 upregulated and 59 downregulated genes. In the PPI network, seven hub genes were identified, including EGF, ERBB2, CDC42, PIK3R2, LCK, CBL, and STAT1. Functional enrichment analysis revealed that these genes were mainly enriched in PI3K-Akt signaling pathway, osteoclast differentiation, and focal adhesion. In conclusion, the results in the current study suggest that pathways of PI3K-Akt, osteoclast differentiation, and focal adhesion may play important roles in the development of OA females after menopause. EGFR, ERBB2, CDC42, and STAT1 may be key genes related to OA progression in postmenopausal women and may be promising therapeutic targets for OA.
| S-EPMC6925789 | biostudies-literature
Project description:The present study aimed to identify key genes involved in osteoarthritis (OA). Based on a bioinformatics analysis of five gene expression profiling datasets (GSE55457, GSE55235, GSE82107, GSE12021 and GSE1919), differentially expressed genes (DEGs) in OA were identified. Subsequently, a protein-protein interaction (PPI) network was constructed and its topological structure was analyzed. In addition, key genes in OA were identified following a principal component analysis (PCA) based on the DEGs in the PPI network. Finally, the functions and pathways enriched by these key genes were also analyzed. The PPI network consisted of 241 nodes and 576 interactives, including a total of 171 upregulated DEGs [e.g., aspartylglucosaminidase (AGA), CD58 and CD86] and a total of 70 downregulated DEGs (e.g., acetyl-CoA carboxylase ? and dihydropyrimidine dehydrogenase). The PPI network complied with an attribute of scale-free small-world network. After PCA, 47 key genes were identified, including ?-1,4-galactosyltransferase-1 (B4GALT1), AGA, CD58, CD86, ezrin, and eukaryotic translation initiation factor 4 ? 1 (EIF4G1). Subsequently, the 47 key genes were identified to be enriched in 13 Gene Ontology (GO) terms and 2 Kyoto Encyclopedia of Genes and Genomes pathways, with the GO terms involving B4GALT1 including positive regulation of developmental processes, protein amino acid terminal glycosylation and protein amino acid terminal N-glycosylation. In addition, B4GALT1 and EIF4G1 were confirmed to be downregulated in OA samples compared with healthy controls, but only EIF4G1 was determined to be significantly downregulated in OA samples, as determined via a meta-analysis of the 5 abovementioned datasets. In conclusion, B4GALT1 and EIF4G1 were indicated to have significant roles in OA, and B4GALT1 may be involved in positive regulation of developmental processes, protein amino acid terminal glycosylation and protein amino acid terminal N-glycosylation. The present study may enhance the current understanding of the molecular mechanisms of OA and provide novel therapeutic targets.
| S-EPMC7678638 | biostudies-literature
Project description:Objective Rheumatoid arthritis (RA) is one of the most prevalent inflammatory arthritis worldwide. However, the genes and pathways associated with macrophages from synovial fluids in RA patients still remain unclear. This study aims to screen and verify differentially expressed genes (DEGs) related to identifying candidate genes related to synovial macrophages in rheumatoid arthritis by bioinformatics analysis. Methods We searched the Gene Expression Omnibus (GEO) database, and GSE97779 and GSE10500 with synovial macrophages expression profiling from multiple RA microarray dataset were selected to conduct a systematic analysis. GSE97779 included nine macrophage samples from synovial fluids of RA patients and five macrophage samples from primary human blood of HC. GSE10500 included five macrophage samples from synovial fluids of RA patients and three macrophage samples from primary human blood of HC. Functional annotation of DEGs was performed, including Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Protein–protein interaction (PPI) network of DEGs was established using the STRING database. CytoHubba was used to identify hub genes. MCODE was used to determine gene clusters in the interactive network. Results There were 2638 DEGs (1425 upregulated genes and 1213 downregulated ones) and 889 DEGs (438 upregulated genes and 451 downregulated ones) selected from GSE97779 and GSE10500, respectively. Venn diagrams showed that 173 genes were upregulated and 106 downregulated in both two datasets. The top 10 hub genes, including FN1, VEGFA, HGF, SERPINA1, MMP9, PPBP, CD44, FPR2, IGF1, and ITGAM, were identified using the PPI network. Conclusion This study provides new insights for the potential biomarkers and the relevant molecular mechanisms in RA patients. Our findings suggest that the 10 candidate genes might be used in diagnosis, prognosis, and therapy of RA in the future. However, further studies are required to confirm the expression of these genes in synovial macrophages in RA and control specimen.
| S-EPMC8575484 | biostudies-literature
Project description:Rheumatoid arthritis (RA) and osteoarthritis (OA) are two major types of joint diseases. The present study aimed to identify hub genes involved in the pathogenesis and further explore the potential treatment targets of RA and OA. The gene expression profile of GSE12021 was downloaded from Gene Expression Omnibus (GEO). Total 31 samples (12 RA, 10 OA and 9 NC samples) were used. The differentially expressed genes (DEGs) in RA versus NC, OA versus NC and RA versus OA groups were screened using limma package. We also verified the DEGs in GSE55235 and GSE100786. Functional annotation and protein-protein interaction (PPI) network construction of OA- and RA-specific DEGs were performed. Finally, the candidate small molecules as potential drugs to treat RA and OA were predicted in CMap database. 165 up-regulated and 163 down-regulated DEGs between RA and NC samples, 73 up-regulated and 293 down-regulated DEGs between OA and NC samples, 92 up-regulated and 98 down-regulated DEGs between RA and OA samples were identified. Immune response and TNF signaling pathway were significantly enriched pathways for RA- and OA-specific DEGs, respectively. The hub genes were mainly associated with 'Primary immunodeficiency' (RA vs. NC group), 'Ribosome' (OA vs. NC group), and 'Chemokine signaling pathway' (RA vs. OA group). Arecoline and Cefamandole were the most promising small molecule to reverse the RA and OA gene expression. Our findings suggest new insights into the underlying pathogenesis of RA and OA, which may improve the diagnosis and treatment of these intractable chronic diseases.
| S-EPMC7744737 | biostudies-literature
Project description:Osteoarthritis (OA) is a chronic arthropathy that occurs in the middle‑aged and elderly population. The present study aimed to identify gene signature differences between synovial cells from OA synovial membrane with and without inflammation, and to explain the potential mechanisms involved. The differentially expressed genes (DEGs) between 12 synovial membrane with inflammation and 12 synovial membrane without inflammation from the dataset GSE46750 were identified using the Gene Expression Omnibus 2R. The DEGs were subjected to enrichment analysis, protein‑protein interaction (PPI) analysis and module analysis. The analysis results were compared with text‑mining results. A total of 174 DEGs were identified. Gene Ontology enrichment results demonstrated that functional molecules encoded by the DEGs primarily had extracellular location, molecular functions predominantly involving 'chemokine activity' and 'cytokine activity', and were associated with biological processes, including 'inflammatory response' and 'immune response'. The Kyoto Encyclopedia of Genes and Genomes results demonstrated that DEGS may function through pathways associated with 'rheumatoid arthritis', 'chemokine signaling pathway', 'complement and coagulation cascades', 'TNF signaling pathway', 'intestinal immune networks for IgA production', 'cytokine‑cytokine receptor interaction', 'allograft rejection', 'Toll‑like receptor signaling pathway' and 'antigen processing and presentation'. The top 10 hub genes [interleukin (IL)6, IL8, matrix metallopeptidase (MMP)9, colony stimulating factor 1 receptor, FOS proto‑oncogene, AP1 transcription factor subunit, insulin‑like growth factor 1, TYRO protein tyrosine kinase binding protein, MMP3, cluster of differentiation (CD)14 and CD163] and four gene modules were identified from the PPI network using Cytoscape. In addition, text‑mining was used to identify the commonly used drugs and their targets for the treatment of OA. It was initially verified whether the results of the present study were useful for the study of OA treatment targets and pathways. The present study provided insight for the molecular mechanisms of OA synovitis. The hub genes and associated pathways derived from analysis may be targets for OA treatment. IL8 and MMP9, which were validated by text‑mining, may be used as molecular targets for the OA treatment, while other hub genes require further validation.
| S-EPMC6236257 | biostudies-other
Project description:Discrimination of rheumatoid arthritis (RA) patients from patients with other inflammatory/degenerative joint diseases or healthy individuals purely on the basis of genes differentially expressed in high-throughput data has proven very difficult. Thus, the present study sought to achieve such discrimination by employing a novel unbiased approach using rule-based classifiers. Three multi-center genome-wide transcriptomic data sets (Affymetrix HG- U133 A/B) from a total of 79 individuals, including 20 healthy controls (control group - CG), as well as 26 osteoarthritis (OA) and 33 RA patients, were used to infer rule- based classifiers to discriminate the disease groups. The rules were ranked with respect to Kiendl’s statistical relevance index, and the resulting rule set was optimized by pruning. The rule sets were inferred separately from data of one of three centers and applied to the two remaining centers for validation. All rules from the optimized rule sets of all centers were used to analyze their biological relevance applying the software Pathway Studio. The optimized rule sets for the three centers contained a total of 29, 20, and 8 rules (including 10, 8, and 4 rules for ‘RA’), respectively. The mean sensitivity for the prediction of RA based on six center-to-center tests was 96% (range 90% to 100%), that for OA 86% (range 40% to 100%). The mean specificity for RA prediction was 94% (range 80% to 100%), that for OA 96% (range 83.3% to 100%). The average overall accuracy of the three different rule-based classifiers was 91% (range 80% to 100%). Unbiased analyses by Pathway Studio of the gene sets obtained by discrimination of RA from OA and CG with rule-based classifiers resulted in the identification of the pathogenetically and/or therapeutically relevant interferon-gamma and GM-CSF pathways. First-time application of rule-based classifiers for the discrimination of RA resulted in high performance, with means for all assessment parameters close to or higher than 90%. In addition, this unbiased, new approach resulted in the identification not only of pathways known to be critical to RA, but also of novel molecules such as serine/threonine kinase 10. Synovial tisssue from healthy joints, OA joints, and RA joints.
2014-02-21 | E-GEOD-55235 | biostudies-arrayexpress