Differentially-expressed mRNAs, microRNAs and long noncoding RNAs in intervertebral disc degeneration identified by RNA-sequencing

Abstract

ABSTRACT The underlying molecular mechanisms of intervertebral disc degeneration (IDD) remain unclear. This study aimed to identify the crucial molecules and explore the function of noncoding RNAs and related pathways in IDD. We randomly selected three samples each from an IDD and a spinal cord injury group (control) for RNA-sequencing. We identified 463 differentially-expressed long noncoding RNAs (lncRNAs), 47 differentially-expressed microRNAs (miRNAs), and 1,334 differentially-expressed mRNAs in IDD. Three hundred fifty-eight lncRNAs as cis-regulators could potentially target 865 genes. Protein–protein interaction (PPI) network analysis confirmed that IL-6, VEGFA, IGF1, MMP9, CXCL8, FGF2, IL1B, CCND1, ITGAM, PTPRC, FOS and PTGS2 were hub genes. We built a competing endogenous RNA (ceRNA) network and identified lncRNA XIST–hsa-miR-4775–PLA2G7 and lncRNA XIST–hsa-miR-424-5p–AMOT/TGFBR3 ceRNA axes. Quantitative real-time PCR (qRT-PCR) was implemented in 15 IDD samples and 15 controls to validate differentially-expressed genes in ceRNA axes. From the ceRNA network, gene ontology (GO) enrichment analysis indicated that noncoding RNAs were associated with several biological processes, including extracellular matrix organization, extracellular structure organization, leukocyte migration, and mesenchyme development. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that noncoding RNAs were associated with several pathways including the AGE-RAGE signaling pathway, PI3K-Akt signaling pathway, axon guidance, and osteoclast differentiation. These results indicate that some specific noncoding RNAs and ceRNA axes may be vital during the development of IDD, and may have potential as alternative diagnostic biomarkers as well as novel therapeutic strategies for IDD.