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.

Keywords

Intervertebral disc degeneration circular RNA Gene Expression Omnibus bioinformatics microRNA sponge

Introduction

Intervertebral disc degeneration (IDD) is an orthopedic disease that affects 97% of individuals 50 years or older worldwide.^1–3 IDD affects the stability of the lumbar region, leading to disc herniation and cervical spondylosis.^4–6 Numerous factors have been associated with IDD, including nutrition, genetics, aging, and mechanical loading.⁷ The pathological basis of IDD; that is, excessive apoptosis, leads to changes in cell composition and a decrease in active cells.⁸ Predominant lower back pain as a result of IDD is difficult to reverse because of changes inside the cells and the cell matrix.⁹ Hence, it is imperative to find new methods for IDD treatment.¹⁰ IDD causes annulus fibrosus of endplate cartilage at the top and bottom of the vertebral disc and affects the nucleus pulposus at the center of the disc. A routine spine surgery technique, spine fusion, is widely used to treat cervical spine instability, intervertebral disc injury, spinal deformity, and lumbar spine degeneration.¹¹ The advent of new technologies in molecular biology such as gene therapy may help to make the treatment of IDD more efficient and effective.¹²

RNA molecules that are covalently closed and have no 5′ cap structure or 3′ polyadenylated tail are called circular RNAs (circRNAs).¹³ The rapid growth of RNA sequencing and bioinformatics has led to the discovery of a large number of circRNA structures with known functions. Four types of circRNAs have been characterized: exon circRNAs, intron circRNAs, exon–intron circRNAs, and intergenic circRNAs.¹⁴ These circRNAs play crucial roles in many pathological and biological processes, including apoptosis, cell cycle, metastasis, invasion and migration, and cell proliferation.¹⁵ CircRNAs act as microRNA (miRNA) sponges that mediate downstream gene expression.¹⁶ Hence, circRNAs continuously provide a cushion to highly fluctuating miRNA effects, and circRNA–miRNA–mRNA networks are activated by overall gene expression that is directed to a specific disease process by the altered expression of circRNAs.¹⁷ Various types of circRNAs, such as circSEMA4B, VMA21, and circRNA-104670, have been shown to play crucial roles in IDD. Several diseases have been linked to the critical function of the circRNA–miRNA–mRNA axis.^18–20

Bioinformatics analysis and DNA sequencing methods have produced large amounts of data and led to the discovery of genomic mutations related to cancer.²¹ Critical roles that circRNAs, miRNAs, and mRNAs play in IDD also have been revealed.^22,23 Accurate predictions of drug therapy mechanisms together with rapid screening of drug targets have become more realistic with the development of bioinformatics tools.

In this study, we used the microarray expression profile datasets GSE67566 and GSE116726 from the Gene Expression Omnibus (GEO) to obtain circRNAs, miRNAs, and mRNAs related to IDD. Gene ontology (GO) enrichment analysis was applied to the differentially expressed genes (DEGs). A novel circRNA–miRNA–mRNA network that may play an important role in disc degeneration was constructed. Hub genes from among the common DEGs were detected by constructing a protein–protein interaction (PPI) network.

Material and methods

Ethical statement

Ethical approval for this study was deemed unnecessary because only data from bioinformatics databases were used.

Microarray expression data

We selected two GEO microarray datasets (GSE67566 and GSE116726) to screen for differentially expressed circRNAs (DECs), differentially expressed miRNAs (DEMs), and DEGs in IDD tissues compared with normal tissues.^24–26 The Agilent-070156 Human miRNA (GPL20712) platform was used for GSE116726, and the Agilent-069978 Arraystar Human CircRNA microarray V1 (GPL19978) platform (Agilent, Santa Clara, CA, USA) was used for GSE67566. The two datasets had eight control (normal) samples and eight IDD samples; GSE116726 had three control and three IDD samples, and GSE67566 had five control and five IDD samples (Table 1).

Table 1.

Details of the GSE67566 and GSE116726 microarray expression profile datasets.

Dataset	Organism	Platform	Last update date	IDD samples	Control samples
GSE67566	Homo sapiens	GPL19978 (Agilent-069978 Arraystar Human CircRNA microarray V1)	5 April 2017	5	5
GSE116726	Homo sapiens	GPL20712 (Agilent-070156 Human miRNA)	11 December 2018	3	3

IDD, intervertebral disc degeneration.

Data analysis to detect DEGs in the IDD tissues

The Bioconductor affy package (https://www.bioconductor.org/packages/release/bioc/html/affy.html) was used to analyze the annotation files for log₂ conversion, standardization, and background correction of the raw data, and the limma package in R (https://www.R-project.org/) was used to screen the DEGs in the IDD tissues. P-values <0.05 and |logFC| >1 (where FC is fold change) were considered significant.

Prediction of miRNA target genes

The top highly upregulated circRNAs were used for further investigation. To predict potential interactions between the miRNAs and circRNAs, we used the CircInteractome database (https://circinteractome.nia.nih.gov/);²⁷ to predict the target genes of the identified miRNAs, we used starBase v3.020 (http://starbase.sysu.edu.cn/). Target genes were identified by the miRNA target prediction programs TargetScan (http://www.targetscan.org), miRDB (http://mirdb.org/), and miRTarBase (http://mirtarbase.mbc.nctu.edu.tw/php/index.php).^28–31

Venn plot

We used the Venn diagram website (http://bioinformatics.psb.ugent.be/webtools/Venn/) to calculate and draw a diagram of overlapping genes. The result of target miRNA prediction was based on top upregulated DEC prediction and DEM from GSE116726. The overlapping genes from TargetScan, miRDB, and miRTarBase were considered miRNA target genes.

Gene Ontology and PPI analysis

The GO enrichment analysis was conducted using Metascape.³² GO enrichment analysis has been applied extensively to gene products or to specific gene annotations.³³ To better understand the relationships among the enriched terms, we constructed a network plot by connecting terms that had similarities >0.3 to form edges. Twenty clusters were formed and the terms with the lowest P-values were selected, with a cap of 250 terms and no more than 15 terms per cluster. Cytoscape was used to visualize the network.³⁴ The PPI analysis for proteins encoded by each listed gene was conducted using the OmniPath (https://omnipathdb.org/), InWeb_IM (https://inbio-discover.intomics.com/map.html#search), and BioGrid (https://wiki.thebiogrid.org/doku.php/statistics/) databases. The subset of proteins that formed physical interactions with at least one other protein was contained within the resultant network. We used the molecular complex detection (MCODE) plugin in Cytoscape to identify the components of the densely connected network.^35–37

Results

Identification of DEGs in IDD tissues

After processing the expression profile data and standardizing the gene expression, we screened the DEGs in the GSE67566 and GSE116726 datasets using bioinformatics tools (Table 1). A total of 970 DEMs were detected in GSE116726 and 105 DECs were detected in GSE67566 using the limma package. Volcano plots of the DEGs from each dataset are shown in Figure 1A and B.

Figure 1.

Differentially expressed genes (DEGs) and microRNAs (miRNAs) in intervertebral disc degeneration (IDD) tissues. (a, b) Volcano plots of DEGs between IDD and normal tissues in the (A) GSE67566 and (b) GSE116726 datasets. Red dots indicate significantly upregulated genes in IDD; green dots indicate significantly downregulated genes in IDD; black dots indicate genes that are not differentially expressed. P < 0.05 and |log2 fold change| >1 were considered significant. (c) Venn diagram of the 970 differentially expressed miRNAs detected from GSE116726 and their overlap with miRNAs from CircInteractome. (d) Venn diagram of the miRNA target genes predicted by TargetScan, miRDB, and miRTarBase.

Circ-0040039 functions as a sponge of several miRNAs

The ability of circRNAs to sponge miRNAs and inhibit their activities has been shown in numerous studies.^38–41 Analysis using the CircInteractome database (https://circinteractome.nia.nih.gov) identified two possible binding site types between hsa-circ-0040039 and its predicted target miRNAs. Twenty-three target miRNAs were identified; hsa-miR-370, hsa-miR-545, and hsa-miR-665 had both types of binding sites, 16 had one type of binding site and 4 had the other type of binding site (Table 2).

Table 2.

Predicted microRNA (miRNA) targets of hsa-circ-0040039.

CircRNA Mirbase ID	CircRNA (top) – miRNA (bottom) pairing
hsa_circ_0040039 (5′…3′)	GGUUGGUUCAUUCUGGCUCCGGA
	ǁ\| ǁǁǁ\|
hsa-miR-1203 (3′…5′)	CUCGACGUAGGACCGAGGCCC
hsa_circ_0040039 (5′…3′)	UCAUGGCUGUGACUGAGAAGGAU
	ǁ\| ǁǁǁ
hsa-miR-1248 (3′…5′)	AAAUCGUGUCACGAAUAUGU-UCUUCCA
hsa_circ_0040039 (5′…3′)	AUGGUGGAAGACAGCAAUGGAGA
	ǁǁǁ\|
hsa-miR-136 (3′…5′)	AGGUAGUAGUUUUGUUUACCUCA
hsa_circ_0040039 (5′…3′)	CAGGACAGGCUCUCGACAGGGCA
	ǁǁǁ\|
hsa-miR-146b-3p (3′…5′)	GGUCUUGACUCAGGUGUCCCGU
hsa_circ_0040039 (5′…3′)	AGUAUCAGAUCUGCCGUGGGAAG
	ǁǁǁ
hsa-miR-188-3p (3′…5′)	ACGUUUGGGACGUACACCCUC
hsa_circ_0040039 (5′…3′)	AUCCCUCUCAAAAUGUGCUUUGC
	ǁǁǁ\|
hsa-miR-330-3p (3′…5′)	AGAGACGUCCGGCACACGAAACG
hsa_circ_0040039 (5′…3′)	AUUCUGGCUCCGGAUGUCGAUCC
	ǁǁǁ\|
hsa-miR-369-5p (3′…5′)	CGCUUAUAUUGUGCCAGCUAGA
hsa_circ_0040039 (5′…3′)	CUACAUUCUCCUGAUAGCAGGAA
	ǁ\| ǁǁǁ
hsa-miR-370 (3′…5′)	UGGUCCAAGGUGGGG-UCGUCCG
hsa_circ_0040039 (5′…3′)	UUGGGGCAACCAGUACAGCAGGA
	ǁǁǁ\|
hsa-miR-370 (3′…5′)	UGGUCCAAGGUGGGGUCGUCCG
hsa_circ_0040039 (5′…3′)	CAUGCUGCUGCUGAGCUGAUCAA
	ǁǁǁ
hsa-miR-383 (3′…5′)	UCGGUGUUAGUGGAAGACUAGA
hsa_circ_0040039 (5′…3′)	UGCUGCUGAGCUGAUCAAGGAAG
	ǁǁǁ
hsa-miR-502-5p (3′…5′)	AUCGUGGGUCUAUCGUUCCUA
hsa_circ_0040039 (5′…3′)	CUCAAAAUGUGCUUUGCUGCUAG
	ǁǁǁ
hsa-miR-503 (3′…5′)	GACGUCUUGACAAGGGCGACGAU
hsa_circ_0040039 (5′…3′)	AAGUAACACCAUAUAUCAAGAAG
	ǁǁǁ
hsa-miR-526b (3′…5′)	UGUCUUUCACGAAGGGAGUUCUC
hsa_circ_0040039 (5′…3′)	CCUCUCAAAAUGUGC-UUUGCUGC
	ǁ\| ǁǁǁ\|
hsa-miR-545 (3′…5′)	CGUGUGUUAUUUACAAACGACU
hsa_circ_0040039 (5′…3′)	UGUGACUGAGAAGGAUUUGCUGC
	ǁǁǁ\|
hsa-miR-545 (3′…5′)	CGUGUGUUAUUUACAAACGACU
hsa_circ_0040039 (5′…3′)	GGAUCUGACCUUACAUUUGCUAC
	ǁǁǁ
hsa-miR-548p (3′…5′)	UUUCAUUGACGUCAAAAACGAU
hsa_circ_0040039 (5′…3′)	AAGAUACAGCCACAG---CACACUCC
	ǁǁ ǁǁǁ\|
hsa-miR-574-5p (3′…5′)	UGUGUGAGUGUGUGUGUGUGAGU
hsa_circ_0040039 (5′…3′)	AAACUAGAUGGUGGAAGACAGCA
	ǁǁǁ\|
hsa-miR-626 (3′…5′)	UUCUGUAAAAGUCUGUCGA
hsa_circ_0040039 (5′…3′)	GUGGGAAGGUGCAGCCCCCCAGU
	ǁǁǁ\|
hsa-miR-637 (3′…5′)	UGCGUCUCGGGCUUUCGGGGGUCA
hsa_circ_0040039 (5′…3′)	UCUGGCUCCGGAUGUCGAUCCCC
	ǁ\| ǁǁǁ\|
hsa-miR-638 (3′…5′)	UCCGGCGGUGGGCGGGC-GCUAGGGA
hsa_circ_0040039 (5′…3′)	CUCAAAAUGUGCUUUGCUGCUAG
	ǁǁǁ
hsa-miR-646 (3′…5′)	CGGAGUCUCCGUCGACGAA
hsa_circ_0040039 (5′…3′)	AAAGAUACAGCCACAGCACACUC
	ǁǁǁ\|
hsa-miR-648 (3′…5′)	UGGUCACGGGACGUGUGAA
hsa_circ_0040039 (5′…3′)	UAGUAUCAGAUCUGCCGUGGGAA
	ǁǁǁ
hsa-miR-662 (3′…5′)	GACGACCCGGUGUUGCACCCU
hsa_circ_0040039 (5′…3′)	UACAGCCACAGCACACUCCUGGU
	ǁǁǁ\|
hsa-miR-665 (3′…5′)	UCCCCGGAGUCGGAGGACCA
hsa_circ_0040039 (5′…3′)	CAUCGGGAUCUGUCAUCCUGGAC
	ǁǁǁ
hsa-miR-665 (3′…5′)	UCCCCGGAGUCGGAGGACCA
hsa_circ_0040039 (5′…3′)	AGGACAGGCUCUCGA--CAGGGCAU
	ǁǁ ǁǁǁ\|
hsa-miR-874 (3′…5′)	AGCCAGGGAGCCCGGUCCCGUC

Binding sites: hsa-miR-370, hsa-miR-545, and hsa-miR-665 had both types of binding sites, 16 had one type of binding site and 4 had the other type of binding site.

Venn plots

Among the 970 DEMs that were detected from GSE116726, five overlapped with miRNAs from the CircInteractome database, as shown in the Venn plot in Figure 1C. The five overlapping DEMs were hsa-miR-648, hsa-miR-638, hsa-miR-502-5p, hsa-miR-574-5p, and hsa-miR-662. Among the predicted miRNA target genes, 46 genes were predicted by all three prediction programs (TargetScan, miRDB, and miRTarBase; Figure 1D).

Gene Ontology and PPI

To investigate the biological functions of the DEGs associated with IDD, GO enrichment analysis was performed using Metascape. The results indicated that the target genes were significantly enriched in vesicle targeting, negative regulation of cell cycle G₂/M phase transition, cell junction assembly, negative regulation of cell cycle, negative regulation of cellular component organization, synapse assembly, and skeletal system morphogenesis (Figure 2). The PPI analysis showed that RAB1A, RAB1B, and MCFD2 play important roles in IDD (Figure 3).

Figure 2.

Gene Ontology (GO) enrichment analysis of the differentially expressed target genes associated with intervertebral disc degeneration tissues. (a) Heatmap of enriched GO terms assigned to the target genes. (b) Network of genes with enriched GO terms.

Figure 3.

Protein–protein interaction of three proteins—two members of the RAS oncogene family (RAB1A, RAB1B) and multiple coagulation factor deficiency (MCFD2)—associated with intervertebral disc degeneration.

Discussion

IDD is considered a global threat to human health because of the significant associated healthcare costs and because it is a gateway to other disc-related diseases.⁴² IDD is caused by a variety of traumatic, mental, and nutritional factors.⁴³ As well as the known behavioral and environmental factors, genetic factors have now been found to be associated with risk of IDD.^44–47 Traditional treatments for disc degeneration include surgery and preventive treatments¹⁰. The recent development of biological techniques to repair IDD has provided new avenues of research and potential new treatments for IDD. It has been suggested that Wnt signaling may offer insight into targets for pharmacologic or nonpharmacologic therapeutics for IDD.⁴⁸ More biomarkers and related signaling pathways must be identified to further improve the treatment of IDD.

In this study, we downloaded microarray expression data from GSE67566 and GSE116726 and used it to detect differences between IDD and normal tissues. The identified DECs and DEMs were used to create a circRNA–miRNA–mRNA network. GO enrichment analysis showed that the mRNAs in the network were significantly enriched in vesicle targeting, negative regulation of cell cycle G₂/M phase transition, cell junction assembly, negative regulation of cell cycle, negative regulation of cellular component organization, synapse assembly, and skeletal system morphogenesis. IDD has been shown to be closely associated with apoptosis of nucleus pulposus cells.⁴⁹ The pathology of IDD is influenced by decreases in the composition and synthesis of the extracellular matrix because of a reduction in the activity of the intervertebral disc cells due to increased cell apoptosis.⁵⁰ The PPI showed that RAB1A, RAB1B, and MCFD2 played important roles in disc degeneration. RAB1A and RAB1B are members of the RAS oncogene family and are associated with neoplasms, infections, and kidney diseases.^51–54 MCFD2, RAB1B, and RAB1A were enriched mainly in vesicle targeting, which is associated with phenotype-keeping of intervertebral disc chondroid cells and enhanced proliferation of notochord cells.⁵⁵ Further studies are necessary to elucidate the actions of the circRNA–miRNA–mRNA network in IDD to understand their potential relevance to IDD.

In summary, we identified a novel circRNA–miRNA–mRNA network related to the progression and initiation of IDD through a comprehensive bioinformatics analysis of two GEO microarray datasets of IDD. The identified hub genes could have potential as biomarkers in the treatment, prognosis, and diagnosis of IDD. Their crucial role in the pathogenesis of IDD was confirmed further through a series of detailed analyses. However, the specific role of hsa-circ-0040039 in the pathogenesis of IDD needs to be further studied. Overall, our results shed light on the potential integration of multiple biomarkers for clinical practice in the diagnosis and treatment of IDD.

Supplemental Material

sj-pdf-1-imr-10.1177_0300060520960983 - Supplemental material for A novel circRNA–miRNA–mRNA network reveals hsa-circ-0040039 as a biomarker for intervertebral disc degeneration

Supplemental material, sj-pdf-1-imr-10.1177_0300060520960983 for A novel circRNA–miRNA–mRNA network reveals hsa-circ-0040039 as a biomarker for intervertebral disc degeneration by Jianhua Tang, Chenlin Zhang, Shengru Wang and Jianfeng Chen in Journal of International Medical Research

Footnotes

Acknowledgment

We thank Margaret Biswas, PhD, from Liwen Bianji, Edanz Group China (), for editing the English text of a draft of this manuscript.

Declaration of conflicting interest

The authors declare that there is no conflict of interest.

Funding

The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the Scientific Research Projects of Wuxi Health Planning Commission (grant number Z201703Z201703); and the Training Program of Innovation and Entrepreneurship for College Students in Jiangsu (grant numbers 201913993001Y, 201913993021H).

References

Wang

, et al. Association between GDF5 single nucleotide polymorphism rs143383 and lumbar disc degeneration. Exp Ther Med 2018; 16: 1900–1904.

Elsaadany

Winters

Adams

, et al. Equiaxial strain modulates adipose-derived stem cell differentiation within 3D biphasic scaffolds towards annulus fibrosus. Sci Rep 2017; 7: 12868–12868.

Wismer

Grad

Fortunato

, et al. Biodegradable electrospun scaffolds for annulus fibrosus tissue engineering: effect of scaffold structure and composition on annulus fibrosus cells in vitro. Tissue Eng Part A 2014; 20: 672–682.

Zhang

Yang

Zhou

, et al. Knockdown of miR-222 inhibits inflammation and the apoptosis of LPS-stimulated human intervertebral disc nucleus pulposus cells. Int J Mol Med 2019; 44: 1357–1365.

Chan

and Gantenbein-Ritter

Intervertebral disc regeneration or repair with biomaterials and stem cell therapy–feasible or fiction?

Swiss Med Wkly 2012; 142: w13598.

Park

Lee

Park

, et al. Mitochondrial involvement in Fas-mediated apoptosis of human lumbar disc cells. J Bone Joint Surg Am 2005; 87: 1338–1342.

Gao

, et al. Degree of endplate chondrocyte degeneration in different tension regions during mechanical stimulation. Mol Med Rep 2018; 17: 4415–4421.

Ren

Guo

, et al. Aloperine attenuates hydrogen peroxide-induced injury via anti-apoptotic activity and suppression of the nuclear factor-κB signaling pathway. Exp Ther Med 2017; 13: 315–320.

Xiang

Lin

Shen

, et al. Construction and assessment of bio-engineered intervertebral discs. Exp Ther Med 2017; 14: 1929–1934.

10.

Huang

Cheng

Zheng

, et al. Vitamin D retards intervertebral disc degeneration through inactivation of the NF-κB pathway in mice. Am J Transl Res 2019; 11: 2496–2506.

11.

Lei

, et al. Application of a novel porous tantalum implant in rabbit anterior lumbar spine fusion model: in vitro and in vivo experiments. Chin Med J (Engl) 2019; 132: 51–62.

12.

Liu

, et al. Lentivirus-mediated TGF-β3, CTGF and TIMP1 gene transduction as a gene therapy for intervertebral disc degeneration in an in vivo rabbit model. Exp Ther Med 2016; 11: 1399–1404.

13.

Jeck

Sorrentino

Wang

, et al. Circular RNAs are abundant, conserved, and associated with ALU repeats. RNA 2013; 19: 141–157.

14.

Yang

and Huang

Circular RNAs in nasopharyngeal carcinoma. Clin Chim Acta 2020; 508: 240–248.

15.

Yuan

Zhou

Wang

, et al. Circular RNA Cdr1as sensitizes bladder cancer to cisplatin by upregulating APAF1 expression through miR-1270 inhibition. Mol Oncol 2019; 13: 1559–1576.

16.

Zheng

Qian

Jiang

, et al. CircRNA LRP6 promotes the development of osteosarcoma via negatively regulating KLF2 and APC levels. Am J Transl Res 2019; 11: 4126–4138.

17.

Lee

Moon

Jeon

, et al. Possible epigenetic regulatory effect of dysregulated circular RNAs in Alzheimer's disease model. Sci Rep 2019; 9: 11956–11956.

18.

Wang

Guo

Cheng

, et al. An integrated analysis of the circRNA-miRNA-mRNA network reveals novel insights into potential mechanisms of cell proliferation during liver regeneration. Artif Cells Nanomed Biotechnol 2019; 47: 3873–3884.

19.

Zhao

Zou

Han

, et al. Circlular RNA BARD1 (hsa_circ_0001098) overexpression in breast cancer cells with TCDD treatment could promote cell apoptosis via miR-3942/BARD1 axis. Cell Cycle 2018; 17: 2731–2744.

20.

Mahmoudi

Kiltschewskij

Fitzsimmons

, et al. Depolarization-associated CircRNA regulate neural gene expression and in some cases may function as templates for translation. Cells 2019; 9: 25.

21.

Sun

Zhu

, et al. Identification of 12 cancer types through genome deep learning. Sci Rep 2019; 9: 17256–17256.

22.

Zhao

, et al. Characterization of microRNA expression profiles in patients with intervertebral disc degeneration. Int J Mol Med 2014; 33: 43–50.

23.

Cheng

Zhang

, et al. Circular RNA VMA21 protects against intervertebral disc degeneration through targeting miR-200c and X linked inhibitor-of-apoptosis protein. Ann Rheum Dis 2018; 77: 770–779.

24.

Jiang

Zhang

, et al. Preclinical development of a microRNA-based therapy for intervertebral disc degeneration. Nat Commun 2018; 9: 5051.

25.

Liu

Che

Xie

, et al. Noncoding RNAs in human intervertebral disc degeneration: An integrated microarray study. Genom Data 2015; 5: 80–81.

26.

Lan

Liu

Pei

, et al. Landscape of RNAs in human lumbar disc degeneration. Oncotarget 2016; 7: 63166–63176.

27.

Dudekula

Panda

Grammatikakis

, et al. CircInteractome: A web tool for exploring circular RNAs and their interacting proteins and microRNAs. RNA Biol 2016; 13: 34–42.

28.

Liu

Zhou

, et al. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res 2014; 42: D92–D97.

29.

Liu

and Wang

Prediction of functional microRNA targets by integrative modeling of microRNA binding and target expression data.

Genome Biol 2019; 20: 18.

30.

Wong

and Wang

miRDB: an online resource for microRNA target prediction and functional annotations.

Nucleic Acids Res 2015; 43: D146–D152.

31.

Chou

Shrestha

Yang

, et al. miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res 2018; 46: D296–D302.

32.

Zhou

Pache

, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun 2019; 10: 1523.

33.

Gene Ontology Consortium. The Gene Ontology (GO) project in 2006. Nucleic Acids Res 2006; 34: D322–D326.

34.

Shannon

Markiel

Ozier

, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res 2003; 13: 2498–2504.

35.

Stark

Breitkreutz

Reguly

, et al. BioGRID: a general repository for interaction datasets. Nucleic Acids Res 2006; 34: D535–D539.

36.

Wernersson

Hansen

, et al. A scored human protein-protein interaction network to catalyze genomic interpretation. Nat Methods 2017; 14: 61–64.

37.

Bader

and Hogue

CW.

An automated method for finding molecular complexes in large protein interaction networks.

BMC Bioinformatics 2003; 4: 2.

38.

Shi

Wang

Feng

, et al. circRNAs and exosomes: a mysterious frontier for human cancer. Mol Ther Nucleic Acids 2019; 19: 384–392.

39.

Naeli

Pourhanifeh

Karimzadeh

, et al. Circular RNAs and gastrointestinal cancers: Epigenetic regulators with a prognostic and therapeutic role. Crit Rev Oncol Hematol 2019; 145: 102854.

40.

Jiang

Shi

, et al. CircRNA: a rising star in gastric cancer. Cell Mol Life Sci 2019.

41.

Liu

Luo

, et al. Circular RNAs: The star molecules in cancer. Mol Aspects Med 2019; 70: 141–152.

42.

Deng

Zhao

Kang

, et al. Elevated interleukin-6 expression levels are associated with intervertebral disc degeneration. Exp Ther Med 2016; 11: 1425–1432.

43.

Xiong

Zhou

Figini

, et al. Multi-parameter evaluation of lumbar intervertebral disc degeneration using quantitative magnetic resonance imaging techniques. Am J Transl Res 2018; 10: 444–454.

44.

Bai

Qiu

, et al. Association between TRAIL gene polymorphisms and the susceptibility and severity of lumbar disc degeneration. Int J Clin Exp Pathol 2015; 8: 7415–7420.

45.

Ala-Kokko

Genetic risk factors for lumbar disc disease.

Ann Med 2002; 34: 42–47.

46.

Kalb

Martirosyan

Kalani

, et al. Genetics of the degenerated intervertebral disc. World Neurosurg 2012; 77: 491–501.

47.

Mayer

Iatridis

Chan

, et al. Genetic polymorphisms associated with intervertebral disc degeneration. Spine J 2013; 13: 299–317.

48.

Holguin

and Silva

MJ.

In-vivo nucleus pulposus-specific regulation of adult murine intervertebral disc degeneration via Wnt/beta-catenin signaling.

Sci Rep 2018; 8: 11191.

49.

Zhan

Wang

, et al. Expression of miR-625 and Fas in cervical vertebral cartilage endplate. Exp Ther Med 2018; 15: 513–519.

50.

Kermani

Hoboubati

Esmaeili-Mahani

, et al. Induction of intervertebral disc cell apoptosis and degeneration by chronic unpredictable stress. J Neurosurg Spine 2014; 20: 578–584.

51.

Yang

, et al. Aberrant amino acid signaling promotes growth and metastasis of hepatocellular carcinomas through Rab1A-dependent activation of mTORC1 by Rab1A. Oncotarget 2015; 6: 20813–20828.

52.

Thomas

Zhang

Wei

, et al. Rab1A is an mTORC1 activator and a colorectal oncogene. Cancer Cell 2014; 26: 754–769.

53.

Nikoshkov

Broliden

Attarha

, et al. Expression pattern of the PRDX2, RAB1A, RAB1B, RAB5A and RAB25 genes in normal and cancer cervical tissues. Int J Oncol 2015; 46: 107–112.

54.

Tanaka

Mun

Harada

, et al. Hsc70 contributes to cancer cell survival by preventing Rab1A degradation under stress conditions. PLoS One 2014; 9: e96785.

55.

Zhao

Liu

Feng

, et al. [ Notochord cells enhance proliferation and phenotype-keeping of intervertebral disc chondroid cells]. Zhongguo Xiu Fu Chong Jian Wai Ke Za Zhi 2008; 22: 939–943.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

1.86 MB