Long noncoding RNA expression profile reveals lncRNAs signature associated with extracellular matrix degradation in kashin-beck disease

Kashin-Beck disease (KBD) is a deformative, endemic osteochondropathy involving degeneration and necrosis of growth plates and articular cartilage. The pathogenesis of KBD is related to gene expression and regulation mechanisms, but long noncoding RNAs (lncRNAs) in KBD have not been investigated. In this study, we identified 316 up-regulated and 631 down-regulated lncRNAs (≥ 2-fold change) in KBD chondrocytes using microarray analysis, of which more than three-quarters were intergenic lncRNAs and antisense lncRNAs. We also identified 232 up-regulated and 427 down-regulated mRNAs (≥ 2-fold change). A lncRNA-mRNA correlation analysis combined 343 lncRNAs and 292 mRNAs to form 509 coding-noncoding gene co-expression networks (CNC networks). Eleven lncRNAs were predicted to have cis-regulated target genes, including NAV2 (neuron navigator 2), TOX (thymocyte selection-associated high mobility group box), LAMA4 (laminin, alpha 4), and DEPTOR (DEP domain containing mTOR-interacting protein). The differentially expressed mRNAs in KBD significantly contribute to biological events associated with the extracellular matrix. Meanwhile, 34 mRNAs and 55 co-expressed lncRNAs constituted a network that influences the extracellular matrix. In the network, FBLN1 and LAMA 4 were the core genes with the highest significance. These novel findings indicate that lncRNAs may play a role in extracellular matrix destruction in KBD.