A novel circRNA–miRNA–mRNA network reveals hsa-circ-0040039 as a biomarker for intervertebral disc degeneration

Abstract

Objective Alterations in the structure and function of intervertebral discs by multifaceted chronic processes can result in intervertebral disc degeneration (IDD). The mechanisms involved in IDD are still unknown. Methods We investigated the possible mechanisms underlying IDD using a bioinformatics analysis of publicly available microarray expression datasets and built a circular RNA–microRNA–mRNA (circRNA–miRNA–mRNA) network based on the results. Datasets GSE67566 and GSE116726 were downloaded from the Gene Expression Omnibus (GEO) and analyzed using the limma package in R. The CircInteractome database was used to detect miRNAs related to circRNA, and TargetScan, miRDB, and miRTarBase were used to predict target mRNAs. Key target genes were annotated using Gene Ontology terms. Results The circRNA hsa-circ-0040039 was found to have the top log fold-change score. Analysis using Metascape showed that the associated genes were enriched mainly in the cell cycle. The Cytoscape plugin MCODE predicted that two members of the RAS oncogene family—RAB1A and RAB1B—and multiple coagulation factor deficiency (MCFD2) may play key roles in IDD. Conclusion Our results suggested that hsa-circ-0040039 and the related network may be potential biomarkers for IDD.