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

Hereditas ◽  
2021 ◽  
Vol 158 (1) ◽  
Author(s):  
Yanzhi Ge ◽  
Zuxiang Chen ◽  
Yanbin Fu ◽  
Xiujuan Xiao ◽  
Haipeng Xu ◽  
...  

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.

2021 ◽  
Author(s):  
Li Tao ◽  
ChaoLiang Xiong ◽  
Li Xue

Abstract Background: Rheumatoid arthritis (RA) is a chronic autoimmune disease characterized by synovitis and subsequent destruction of cartilage and bone. This study aimed to explore RA-related gene markers and the underlying molecular mechanism.Material and Methods: The expression profiles of GSE77298, GSE55235 and GSE12021 were obtained from the Gene Expression Omnibus database. Then, the differential gene expression analysis was conducted between GSE77298 and GSE55235 datasets. Limma package and a Venn diagram were utilized to screen the overlapping differentially expressed genes (DEGs), and Functional enrichment and pathway analysis were performed by using DAVID database. Subsequently, a protein-protein interaction (PPI) network was established, and candidate hub genes were recognized by using STRING and Cytoscape software. Finally, another dataset (GSE12021) was used for the validation of diagnostic value of the candidate hub genes and to identify real hub genes by using receiver operating characteristic (ROC) curves.Results: A total of 385 DEGs were detected, which include 19 downregulated genes and 366 upregulated genes. GO and KEGG pathway analysis showed that DEGs was mainly enriched in various immune and inflammatory response-related functions and pathways. The PPI network was composed of 374 nodes and 767 edges. A total of 8 real hub genes (HLA-DRA, HLA-DRB1, LCK, VAV1, HLA-DPA1, HLA-DPB1, C3AR1 and CD3D) which displayed an excellent diagnostic value for RA were identified.Conclusion: these findings may provide novel and reliable biomarkers for RA, which have some interesting implications for early diagnosis, prognosis and targeted therapy.


2021 ◽  
Author(s):  
Shuo Xu ◽  
Dingsheng Liu ◽  
Mingming Cui ◽  
Yao Zhang ◽  
Yu Zhang ◽  
...  

Abstract Background: Colon adenocarcinoma (COAD) is among the most common digestive system malignancies worldwide. Although some molecular analyses of colon cancer were previously performed, the pathogenesis and gene signatures remain unclear. This study explored the genetic characteristics and molecular mechanisms underlying colon cancer development. Methods: Three gene expression data sets were obtained from the Gene Expression Omnibus (GEO) database. GEO2R was used to determine differentially expressed genes (DEGs) between COAD and normal tissues. Then, the intersection of the data sets was obtained. Metascape was used to perform the functional enrichment analyses. Next, STRING was used to build protein-protein interaction (PPI) networks. Hub genes were identified and analysed using Cytoscape. Next, survival analysis of the hub genes was performed. To explore the early diagnostic value of the hub genes, UALCAN was used to analyse expression characteristics. Finally, alterations in the hub genes were predicted and analysed by cBioPortal. Results: Altogether 436 DEGs were detected. The DEGs were mainly enriched in cell cycle phase transition, nuclear division, meiotic nuclear division, cytokinesis. Based on PPI networks, 20 hub genes were selected. Among them, 6 hub genes (CCNB1, CCNA2, AURKA, NCAPG, DLGAP5, and CENPE) showed significant prognostic value in colon cancer (p<0.05), while 5 hub genes (CDK1, CCNB1, CCNA2, MAD2L1, and DLGAP5) were associated with early colon cancer diagnosis. Conclusions: In conclusion, integrated bioinformatics analysis was used to identify hub genes that reveal the potential mechanism of carcinogenesis and progression of colon cancer. The hub genes might be novel biomarkers for early diagnosis, treatment, and prognosis of colon cancer.


2021 ◽  
Vol 80 (Suppl 1) ◽  
pp. 111-112
Author(s):  
J. Qiao ◽  
S. X. Zhang ◽  
H. Wang ◽  
J. Q. Zhang ◽  
M. T. Qiu ◽  
...  

Background:Rheumatoid arthritis (RA) is an aggressive immune-mediated joint disease characterized by synovial proliferation and inflammation, cartilage destruction, and joint destruction1. Despite efforts to characterize the disease subsets and to predict the differential prognosis in RA patients, disease heterogeneity is not adequately translated into the current clinical subclassification2.Objectives:To develop and validate an integrative system approach for stratifying patients with RA according to disease status and whole-genome gene expression data.Methods:An RNA sequencing dataset of synovial tissues from 124 RA patients (including 57 patients with early RA, 95 with established RA) and 15 healthy controls (HC) was imported from the Gene Expression Omnibus (GEO) database (GSE89408) by software package R (version 4.0.3). After filtrating of differentially expressed genes (DEGs) between RA and HC, non-negative matrix factorization, functional enrichment, and immune cell infiltration were applied to illustrate the landscapes of these patients for classification. Clinical features (age, gender, and auto-antibodies) were also compared to discover the signatures of these classifications.Results:A matrix of 576 DEGs from RA samples was classified into 5 subtypes (early/C1–C3, established/C4-C5) with distinct molecular and cellular signatures and two sub-groups (S1 and S2) (Figure 1A-1D). New-onset patients (early C2) and established C4 patients were named as S1, they shared similar gene signatures mainly characterized by prominent immune cells and proinflammatory signatures, and enriched in the chemokine-mediated signaling pathway, lymphocyte activation, response to bacterium and Primary immunodeficiency. S2(C1, C3 and C5) were more occupied by synovial fibroblasts of destructive phenotype. They were mainly enriched in the response to external factors and PPAR signaling pathway (Figure 1E-1H). Interestingly, combined with clinical information, S1 and S2 had no significance in age and gender (P > 0.05). But patients in S1 had a stronger association with the presence of anti-citrullinated protein antibodies (ACPA) (P < 0.05) (Figure 1I-1J).Conclusion:We successfully deconvoluted RA synovial tissues into pathobiological discrete subsets using an unsupervised machine learning method and described their distinct molecular and cellular characteristics. These results provide important insights into divergent and shared mechanistic features of RA and serve as a template for future studies to guide drug tar-get discovery by synovial molecular signatures and de-sign stratified approaches for patients with RA.References:[1]Smolen JS, Aletaha D, McInnes IB. Rheumatoid arthritis. Lancet 2016;388(10055):2023-38. doi: 10.1016/S0140-6736(16)30173-8 [published Online First: 2016/10/30][2]Jung SM, Park KS, Kim KJ. Deep phenotyping of synovial molecular signatures by integrative systems analysis in rheumatoid arthritis. Rheumatology (Oxford) 2020 doi: 10.1093/rheumatology/keaa751 [published Online First: 2020/11/25]Acknowledgements:This project was supported by National Science Foundation of China (82001740), Open Fund from the Key Laboratory of Cellular Physiology (Shanxi Medical University) (KLCP2019) and Innovation Plan for Postgraduate Education in Shanxi Province (2020BY078).Disclosure of Interests:None declared


2021 ◽  
Vol 11 ◽  
Author(s):  
Fu Chen ◽  
Yong Zhou ◽  
Zhiyuan Wu ◽  
Yunze Li ◽  
Wenlong Zhou ◽  
...  

BackgroundAs the incidence of nonalcoholic fatty liver disease (NAFLD) increases globally, nonalcoholic steatohepatitis (NASH) has become the second common cause of liver transplantation for liver diseases. Recent evidence shows that Roux-en-Y gastric bypass (RYGB) surgery obviously alleviates NASH. However, the mechanism underlying RYGB induced NASH improvement is still elusive.MethodsWe obtained datasets, including hepatic gene expression data and histologic NASH status, at baseline and 1 year after RYGB surgery. Differentially expressed genes (DEGs) were identified comparing gene expression before and after RYGB surgery in each dataset. Common DEGs were obtained between both datasets and further subjected to functional and pathway enrichment analysis. Protein–protein interaction (PPI) network was constructed, and key modules and hub genes were also identified.ResultsIn the present study, GSE106737 and GSE83452 datasets were included. One hundred thirty common DEGs (29 up-regulated and 101 down-regulated) were identified between GSE106737 and GSE83452 datasets. KEGG analysis showed that mineral absorption, IL-17 signaling pathway, osteoclast differentiation, and TNF signaling pathway were significantly enriched. Based on the PPI network, IGF1, JUN, FOS, LDLR, TYROBP, DUSP1, CXCR4, ATF3, CXCL2, EGR1, SAA1, CTSS, and PPARA were identified as hub genes, and three functional modules were also extracted.ConclusionThis study identifies the global gene expression change in the liver of NASH patients before and after RYGB surgery in a bioinformatic method. Our findings will contribute to the understanding of molecular biological changes underlying NASH improvement after RYGB surgery.


2021 ◽  
Author(s):  
Jiaying Shi ◽  
Shule Wang ◽  
Jingfei Zhang ◽  
Xueli Chang ◽  
Juan Wang ◽  
...  

Abstract Background: Immune-mediated necrotizing myopathy (IMNM) is a type of autoimmune myopathy with limited therapeutic measures. This study aims to elucidate the potential biomarkers and investigate the underlying mechanisms in IMNM. Materials and Methods: Microarray datasets in GSE128470 and GSE39454 were obtained from the Gene Expression Omnibus. Differentially expressed genes (DEGs) were filtrated by limma package in R statistical software. Functional enrichment analyses were performed using DAVID online tools. STRING database was used to construct protein‑protein interaction (PPI) networks. The module analysis and hub genes validation were performed using Cytoscape software. Results: Integrated analysis of two databases revealed 160 co-expressed DEGs in IMNM, including 80 downregulated genes and 80 upregulated genes. GO enrichment analysis revealed that sarcomere is the most significantly enriched GO term within the DEGs. KEGG pathway enrichment analysis revealed significant enrichment pathways in cancer. A PPI network consisting of 115 nodes and 205 edges were constructed and top 20 hub genes were identified. Two key modules from the network were identified. Eight hub genes in module 1 (MYH3, MYH7B, MYH8, MYL5, MYBPH, ACTC1, YNNT 2 and MYOG) are tightly associated with skeletal muscle construction. Seven hub genes in module 2 (C1QA, TYROBP, MS4A6A, RNASE6, FCGR2A, FCER1G and LAPTM5) mainly take part in immune response. Conclusions: Our research indicated that cancer-related pathways, skeletal muscle construction pathways and immune-mediated pathways might participate in the development of IMNM. Identified hub genes may serve as potential biomarkers or targets for early diagnosis.


2021 ◽  
Vol 16 (1) ◽  
Author(s):  
Liuxun Li ◽  
Xiaokang Du ◽  
Haiqian Ling ◽  
Yuhang Li ◽  
Xuemin Wu ◽  
...  

Abstract Background Sciatic nerve injury (SNI), which frequently occurs under the traumatic hip and hip fracture dislocation, induces serious complications such as motor and sensory loss, muscle atrophy, or even disabling. The present work aimed to determine the regulating factors and gene network related to the SNI pathology. Methods Sciatic nerve injury dataset GSE18803 with 24 samples was divided into adult group and neonate group. Weighted gene co-expression network analysis (WGCNA) was carried out to identify modules associated with SNI in the two groups. Moreover, differentially expressed genes (DEGs) were determined from every group, separately. Subsequently, co-expression network and protein–protein interaction (PPI) network were overlapped to identify hub genes, while functional enrichment and Reactome analysis were used for a comprehensive analysis of potential pathways. GSE30165 was used as the test set for investigating the hub gene involvement within SNI. Gene set enrichment analysis (GSEA) was performed separately using difference between samples and gene expression level as phenotype label to further prove SNI-related signaling pathways. In addition, immune infiltration analysis was accomplished by CIBERSORT. Finally, Drug–Gene Interaction database (DGIdb) was employed for predicting the possible therapeutic agents. Results 14 SNI status modules and 97 DEGs were identified in adult group, while 15 modules and 21 DEGs in neonate group. A total of 12 hub genes was overlapping from co-expression and PPI network. After the results from both test and training sets were overlapped, we verified that the ten real hub genes showed remarkably up-regulation within SNI. According to functional enrichment of hub genes, the above genes participated in the immune effector process, inflammatory responses, the antigen processing and presentation, and the phagocytosis. GSEA also supported that gene sets with the highest significance were mostly related to the cytokine–cytokine receptor interaction. Analysis of hub genes possible related signaling pathways using gene expression level as phenotype label revealed an enrichment involved in Lysosome, Chemokine signaling pathway, and Neurotrophin signaling pathway. Immune infiltration analysis showed that Macrophages M2 and Regulatory T cells may participate in the development of SNI. At last, 25 drugs were screened from DGIdb to improve SNI treatment. Conclusions The gene expression network is determined in the present work based on the related regulating factors within SNI, which sheds more light on SNI pathology and offers the possible biomarkers and therapeutic targets in subsequent research.


Open Medicine ◽  
2021 ◽  
Vol 17 (1) ◽  
pp. 96-112
Author(s):  
Chengrui Li ◽  
Yufeng Wan ◽  
Weijun Deng ◽  
Fan Fei ◽  
Linlin Wang ◽  
...  

Abstract Lung adenocarcinoma (LUAD) is the most common subtype of non-small cell lung cancer associated with an unstable prognosis. Thus, there is an urgent demand for the identification of novel diagnostic and prognostic biomarkers as well as targeted drugs for LUAD. The present study aimed to identify potential new biomarkers associated with the pathogenesis and prognosis of LUAD. Three microarray datasets (GSE10072, GSE31210, and GSE40791) from the Gene Expression Omnibus database were integrated to identify the differentially expressed genes (DEGs) in normal and LUAD samples using the limma package. Bioinformatics tools were used to perform functional and signaling pathway enrichment analyses for the DEGs. The expression and prognostic values of the hub genes were further evaluated by Gene Expression Profiling Interactive Analysis and real-time quantitative polymerase chain reaction. Furthermore, we mined the “Connectivity Map” (CMap) to explore candidate small molecules that can reverse the tumoral of LUAD based on the DEGs. A total of 505 DEGs were identified, which included 337 downregulated and 168 upregulated genes. The PPI network was established with 1,860 interactions and 373 nodes. The most significant pathway and functional enrichment associated with the genes were cell adhesion and extracellular matrix-receptor interaction, respectively. Seven DEGs with high connectivity degrees (ZWINT, RRM2, NDC80, KIF4A, CEP55, CENPU, and CENPF) that were significantly associated with worse survival were chosen as hub genes. Lastly, top 20 most important small molecules which reverses the LUAD gene expressions were identified. The findings contribute to revealing the molecular mechanisms of the initiation and progression of LUAD and provide new insights for integrating multiple biomarkers in clinical practice.


2020 ◽  
Vol 40 (12) ◽  
Author(s):  
Bin Zuo ◽  
JunFeng Zhu ◽  
Fei Xiao ◽  
ChengLong Wang ◽  
Yun Shen ◽  
...  

Abstract Background: 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. Methods: 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. Results: 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. Conclusion: 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.


2020 ◽  
Vol 23 (8) ◽  
pp. 805-813
Author(s):  
Ai Jiang ◽  
Peng Xu ◽  
Zhenda Zhao ◽  
Qizhao Tan ◽  
Shang Sun ◽  
...  

Background: Osteoarthritis (OA) is a joint disease that leads to a high disability rate and a low quality of life. With the development of modern molecular biology techniques, some key genes and diagnostic markers have been reported. However, the etiology and pathogenesis of OA are still unknown. Objective: To develop a gene signature in OA. Method: In this study, five microarray data sets were integrated to conduct a comprehensive network and pathway analysis of the biological functions of OA related genes, which can provide valuable information and further explore the etiology and pathogenesis of OA. Results and Discussion: Differential expression analysis identified 180 genes with significantly expressed expression in OA. Functional enrichment analysis showed that the up-regulated genes were associated with rheumatoid arthritis (p < 0.01). Down-regulated genes regulate the biological processes of negative regulation of kinase activity and some signaling pathways such as MAPK signaling pathway (p < 0.001) and IL-17 signaling pathway (p < 0.001). In addition, the OA specific protein-protein interaction (PPI) network was constructed based on the differentially expressed genes. The analysis of network topological attributes showed that differentially upregulated VEGFA, MYC, ATF3 and JUN genes were hub genes of the network, which may influence the occurrence and development of OA through regulating cell cycle or apoptosis, and were potential biomarkers of OA. Finally, the support vector machine (SVM) method was used to establish the diagnosis model of OA, which not only had excellent predictive power in internal and external data sets (AUC > 0.9), but also had high predictive performance in different chip platforms (AUC > 0.9) and also had effective ability in blood samples (AUC > 0.8). Conclusion: The 4-genes diagnostic model may be of great help to the early diagnosis and prediction of OA.


2021 ◽  
Vol 80 (Suppl 1) ◽  
pp. 460.1-460
Author(s):  
L. Cheng ◽  
S. X. Zhang ◽  
S. Song ◽  
C. Zheng ◽  
X. Sun ◽  
...  

Background:Rheumatoid arthritis (RA) is a chronic, inflammatory synovitis based systemic disease of unknown etiology1. The genes and pathways in the inflamed synovium of RA patients are poorly understood.Objectives:This study aims to identify differentially expressed genes (DEGs) associated with the progression of synovitis in RA using bioinformatics analysis and explore its pathogenesis2.Methods:RA expression profile microarray data GSE89408 were acquired from the public gene chip database (GEO), including 152 synovial tissue samples from RA and 28 healthy synovial tissue samples. The DEGs of RA synovial tissues were screened by adopting the R software. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed. Protein-protein interaction (PPI) networks were assembled with Cytoscape software.Results:A total of 654 DEGs (268 up-regulated genes and 386 down-regulated genes) were obtained by the differential analysis. The GO enrichment results showed that the up-regulated genes were significantly enriched in the biological processes of myeloid leukocyte activation, cellular response to interferon-gamma and immune response-regulating signaling pathway, and the down-regulated genes were significantly enriched in the biological processes of extracellular matrix, retinoid metabolic process and regulation of lipid metabolic process. The KEGG annotation showed the up-regulated genes mainly participated in the staphylococcus aureus infection, chemokine signaling pathway, lysosome signaling pathway and the down-regulated genes mainly participated in the PPAR signaling pathway, AMPK signaling pathway, ECM-receptor interaction and so on. The 9 hub genes (PTPRC, TLR2, tyrobp, CTSS, CCL2, CCR5, B2M, fcgr1a and PPBP) were obtained based on the String database model by using the Cytoscape software and cytoHubba plugin3.Conclusion:The findings identified the molecular mechanisms and the key hub genes of pathogenesis and progression of RA.References:[1]Xiong Y, Mi BB, Liu MF, et al. Bioinformatics Analysis and Identification of Genes and Molecular Pathways Involved in Synovial Inflammation in Rheumatoid Arthritis. Med Sci Monit 2019;25:2246-56. doi: 10.12659/MSM.915451 [published Online First: 2019/03/28][2]Mun S, Lee J, Park A, et al. Proteomics Approach for the Discovery of Rheumatoid Arthritis Biomarkers Using Mass Spectrometry. Int J Mol Sci 2019;20(18) doi: 10.3390/ijms20184368 [published Online First: 2019/09/08][3]Zhu N, Hou J, Wu Y, et al. Identification of key genes in rheumatoid arthritis and osteoarthritis based on bioinformatics analysis. Medicine (Baltimore) 2018;97(22):e10997. doi: 10.1097/MD.0000000000010997 [published Online First: 2018/06/01]Acknowledgements:This project was supported by National Science Foundation of China (82001740), Open Fund from the Key Laboratory of Cellular Physiology (Shanxi Medical University) (KLCP2019) and Innovation Plan for Postgraduate Education in Shanxi Province (2020BY078).Disclosure of Interests:None declared


Sign in / Sign up

Export Citation Format

Share Document