Identification of Key lncRNAs, circRNAs, and mRNAs in Osteoarthritis via Bioinformatics Analysis

Osteoarthritis (OA) is a common degenerative joint disorder that adversely affects the quality of life of patients. Identification of novel diagnostic biomarkers is pivotal for the early detection and prevention of OA. Dataset GSE185059 was selected from Gene Expression Omnibus database to obtain differentially expressed lncRNAs (DE-lncRNAs), mRNAs (DE-mRNAs), and circRNAs (DE-circRNAs) between OA and normal samples. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses as well as protein-protein interaction (PPI) network construction of DE-mRNAs were conducted. Hub genes were identified from PPI networks and validated by RT-qPCR. starBase database was utilized for predicting miRNAs binding with hub genes, selected DE-lncRNAs and DE-circRNAs, respectively. The competing endogenous RNA (ceRNA) networks were constructed. A total of 818 DE-mRNAs, 191 DE-lncRNAs, and 2053 DE-circRNAs were identified. The DE-mRNAs were significantly enriched in several inflammation-related GO terms and KEGG pathways such as positive regulation of cell-cell adhesion, TNF-alpha signaling pathway and NF-kappa B signaling pathway. Thirteen hub genes were identified, which were CFTR, GART, SMAD2, NCK1, TJP1, UBE2D1, EFTUD2, PRKACB, IL10, SNRPG, CHD4, RPS24, and SRSF6. OA-related DE-lncRNA/circRNA-miRNA-hub gene networks were constructed. We identified 13 hub genes and constructed the ceRNA networks related to OA, providing a theoretical basis for further research.