Identification and functional analysis of long non-coding RNAs in the synovial membrane of osteoarthritis patients

Abstract

Osteoarthritis (OA), the most common chronic joint disease in the elderly, has become a significant economic burden for families and societies worldwide. Although treatments are continually improving, current drugs only target joint pain, with no effective therapies modifying OA progression. Long noncoding RNAs (lncRNAs), which have received increasing attention in recent years, are abnormally expressed in OA cartilage. In the present study, weighted coexpression network analysis (WGCNA) was applied to identify modules related to certain OA clinical traits. In total, 4404 coding genes and 161 lncRNAs were differentially expressed based on two OA expression profile data sets and normal control samples. Subsequently, 11 independent modules were acquired, and the green module, with a total of 49 hub genes, was identified as the most relevant to OA. These hub genes were validated using the GSE12021 data set. There was only one lncRNA among the hub genes, namely, NONHSAG034351. Thus, we further explored the function of NONHSAG034351-related genes in the network. Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that NONHSAG034351-associated genes are involved in the response to lipopolysaccharide, angiogenesis, tumour necrosis factor (TNF) signalling, and mitogen-activated protein kinase (MAPK) signalling pathways. In conclusion, we identified modules through WGCNA related to OA clinical traits. NONHSAG034351, the only hub-lncRNA, was downregulated in OA synovial tissue and might play a significant role in the pathological progression of this disease. Our findings have important clinical implications and could provide novel biomarkers that indicate the molecular mechanisms of OA and act as potential therapeutic targets. Significance of this study: Long noncoding RNAs (lncRNAs) have been reported to be abnormally expressed in osteoarthritis (OA), which is the most common chronic joint disease among the elderly. In the present study, we report the expression profiles of lncRNAs in OA and the identification of modules through WGCNA related to OA clinical traits. NONHSAG034351, the only hub-lncRNA identified to be downregulated in the synovial tissue of OA patients, might play a significant role in the pathological progression of OA. Furthermore, our findings provide novel biomarkers associated with the molecular mechanisms underlying OA pathogenesis, thus implying potential therapeutic targets with important clinical implications.