POTENTIAL KEY BIOMARKERS OF OSTEOARTHRITIS BY INTEGRATIVE GENOME-WIDE GENE EXPRESSION PROFILING ANALYSIS
Rongqiang Zhang,Aimin Yang,Xiaomei Ren,Jie Zhang,Xiaoli Yang,Qiling Liu,Na Sun,Puwei Yuan,Yongmin Xiong
DOI: https://doi.org/10.1016/j.joca.2018.02.383
IF: 7.507
2018-01-01
Osteoarthritis and Cartilage
Abstract:Purpose: Osteoarthritis (OA), a disorder of movable joints characterized by cell stress and extracellular matrix degradation, is initiated by micro- and macro-injuries that activate maladaptive repair responses including pro-inflammatory pathways of innate immunity. Some key genes and diagnostic markers had been identified using microarrays, however, they are not able to entirely reveal the complex pathogenesis of OA. To better understand its mechanism, we conducted an updated integrative analysis to identify the key biomarkers of synovial tissue from OA patients. Methods: We screened two new expression profiling datasets (GSE 55235, GSE 55457) with the following inclusion criteria: detection using the same platform, samples from the synovial tissue of OA patients and healthy controls and sample size >10 from the NCBI GEO database. Normalization was performed to minimize inconsistency and heterogeneity of the raw data. Differentially expressed genes (DEGs, p<0.05, FDR<0.05, Fold Change >2) between OA patients and controls from the 2 datasets were explored using R v3.4.0 and GENEVESTIGATOR. Gene Ontology (GO) and KEGG pathway enrichment analyses of the DEGs were performed using Cytoscape 3.4.0. The protein-protein interaction (PPI) network of the DEGs was obtained from STRING database v9.05. Results: 793 and 275 DEGs from GSE 55235 and GSE 55457 were obtained, respectively. Among them, 44 DEGs including 26 up-regulated genes (SCRG1, CX3CR1, PTN, NAP1L3, etc.) and 18 down-regulated genes (ELF1, ZNF24, AKAP13, RLF, etc.), exhibited identical expression trends in the 2 datasets (Fig. 1). The DEGs were mainly involved in osteoclast differentiation and the signaling pathways of Toll-like receptor, TNF and PI3K-Akt. JUN was the core gene of the PPI network of 44 DEGs (Fig. 2A) and expression of JUN in OA patients was clearly decreased (Fig. 2B). The receiver operating characteristic curve (ROC) of JUN for OA diagnosis was shown in Fig. 2C. The area under the curve (AUC) was 0.728 (95% CI: 0.566-0.889, P=0.014), with a high specificity of 100% and a fair sensitivity of 50%. Conclusions: Several apoptosis and inflammation related pathways were elevated in OA, suggesting that apoptosis may play an important role in the occurrence of OA. JUN (mainly relate to malignancies and inflammation) may be a key biomarker of OA. Our results provide some insights into the molecular mechanisms underlying OA, which could aid the diagnosis and treatment of this disease. Rigorously-designed experiments are required to validate these findings in further studies.View Large Image Figure ViewerDownload Hi-res image Download (PPT)