Prediction of a Potential Mechanism of Intervertebral Disc Degeneration Based on a Novel Competitive Endogenous RNA Network

Huang, Junshen; Li, Yuxi; Ye, Ziwei; Cheng, Ziying; Huang, Jiajun; Lu, Shixin; Su, Kaihui; Liang, Yuwei; Li, Ming; Huang, Lin

doi:https://doi.org/10.1155/2021/6618834

BioMed Research International

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusions Data Availability Conflicts of Interest Authors’ Contributions Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 6618834 | https://doi.org/10.1155/2021/6618834

Prediction of a Potential Mechanism of Intervertebral Disc Degeneration Based on a Novel Competitive Endogenous RNA Network

Junshen Huang,¹Yuxi Li,¹Ziwei Ye,²Ziying Cheng,¹Jiajun Huang,¹Shixin Lu,¹Kaihui Su,¹Yuwei Liang,¹Ming Li,¹and Lin Huang¹

Academic Editor: Guoxian Yu

Received14 Oct 2020

Revised19 Apr 2021

Accepted17 Jun 2021

Published30 Jun 2021

Abstract

Low back pain which resulted from intervertebral disc degeneration (IDD) is a common health problem that afflicts people all over the world. Due to the lack of an overall understanding of the molecular interactions involved in IDD, we hope to better understand the pathogenetic mechanisms that drive the degenerative process. The purpose of this study is to obtain mRNAs, miRNAs, lncRNAs, and circRNAs associated with IDD gained from public databases and to establish an interaction network. According to the results of microarray analysis and bioinformatics analysis from the contrast of IDD and normal nucleus pulposus tissues, a total of 49 mRNAs, 10 miRNAs, 30 lncRNAs, and 4 circRNAs were obtained and a lncRNA/circRNA–miRNA–mRNA interaction network was constructed. NEAT1–miR-5100–COL10A1 and miR663AHG/HEIH/hsa-circ-0003600–miR-4741–HAS2/HYAL1/LYVE1 might be potential interaction axes of the molecular mechanism in IDD. The increased expression of NEAT1 might inhibit miR-5100 and subsequently upregulate the expression of COL10A1, which leads to IDD, while the increased expression of miR663AHG/HEIH/hsa-circ-0003600 might inhibit miR-4741 and indirectly upregulate HAS2/HYAL1/LYVE1, and leads to the protection from IDD. More interaction axes are to be exploited to provide theoretical bases for further study on IDD.

1. Introduction

Low back pain (LBP) is a common health problem that afflicts people all over the world, and it is estimated that around 84% of the world’s population will experience low back pain in their lives [1]. Not only does it affect patients’ quality of life, but it also has a broad socioeconomic impact, with some recent studies showing that related costs in the United States exceed $100 billion a year [2, 3]. Importantly, intervertebral disc degeneration (IDD) is one of the major risk factors related to LBP [4]. It has been determined that 40% of chronic LBP is caused by IDD [5, 6], and the main factors affecting the pathophysiology of IDD include genetic susceptibility, old age, smoking, alcohol addiction, obesity, and diabetes [7, 8]. Currently, there is no specific treatment for IDD-induced LBP. The current treatment strategies of IDD-induced LBP focus only on pharmacological approaches to relieve pain or inflammation, followed by surgery as a last resort for severe disease stages. With the lack of long-term efficacy and increasing drug abuse, it is critical to develop a molecular targeted therapy for IDD, which may fundamentally prevent IDD or restore the function of the intervertebral disc to solve the problem [9, 10].

IDD is a chronic pathological process with multiple complex etiologies [11]. The intervertebral disc is composed of the inner-most glycosaminoglycan, notochord-derived nucleus pulpous (NP), fibrocartilaginous annulus fibrosus, and superior and inferior cartilaginous end plates. As NP ages and degenerates, it gradually loses its water-bound matrix and becomes fibrotic and its ability to absorb mechanical loads is also diminished, which is accompanied by the phenotype changes in nucleus pulposus cells (NPC) [12]. NPC regulates the balance of synthesis/catabolism by secreting specific biological factors such as CTGF, Shh, and TGF-1, and the secretion changes of these biological factors may also lead to IDD [13]. Understanding the changes of NPC in IDD is essential for the pathogenesis and treatment of the disease.

So far, there have been many reports on the potential influence of IDD genes, such as COL1A1, COL9A2, COL9A3, COL11A2, IL-6, AGC1, VDR, and MMP-3 [14–18]. However, increasing evidences also show that noncoding RNAs (ncRNAs) including miRNA, lncRNA, and circRNA can exert effects on coding RNAs and influence biological processes such as cell proliferation and apoptosis. miRNAs act by binding to complementary sequences in the 3-untranslated region (UTR) of their target mRNA to trigger transcriptional inhibition or mRNA degradation [19]. lncRNAs and circRNAs can competitively bind to miRNAs through their miRNA reaction elements, thus acting as competitive endogenous RNAs (ceRNAs) to regulate miRNA-target mRNAs’ expression levels [20]. More and more studies are revealing the influence of the molecular expression in NPC in the process of IDD disease. However, these results are relatively independent, which cannot provide us a more comprehensive understanding of the molecular expression in the NPC of IDD. Meanwhile, there are still many unknown genes to be explored for a more comprehensive explanation of IDD. There have been few reports on the construction of IDD-related lncRNA/circRNA–miRNA–mRNA networks so far. Therefore, the goal of this study is to filtrate differentially expressed RNAs by using multiple microarray datasets collected from public databases to further construct the regulatory mechanism network of lncRNA/circRNA–miRNA–mRNA in IDD. We aim to identify key molecules in the IDD process and their associated interaction axes, which may provide possible targets for further study.

2. Materials and Methods

An overall framework of this study, firstly, is shown in Figure 1. The logical framework is mainly divided into “data collection and screening,” “data analysis,” and “network construction”.

Figure 1

The logical framework is mainly divided into data collection and screening, data analysis, and network construction. The datasets of mRNA, miRNA, lncRNA, and circRNA correlated to IDD were collected from the public database, and DEMs, DEGs, DELs, and DECs with significant expression differences were screened out. Bioinformatics analysis tools and databases were used to conduct regulatory prediction analysis on DEGs, DELs, and DECs and obtain miRNAs paired with them. DEMs was then matched to establish a complete lncRNA/circRNA-miRNA-mRNA network. Then, the network was used to conduct ceRNA analysis and functional enrichment analysis. IDD: intervertebral disc degeneration; DEM: differentially expressed miRNA; DEG: differentially expressed mRNA; DEL: differentially expressed lncRNA; DEC: differentially expressed circRNA; ceRNA: competitive endogenous RNA.

2.1. Data Collection and Screening

Public functional genomics data was obtained through Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/). To obtain microarray expression datasets for IDD mRNA, miRNA, lncRNA, and circRNA, we used the search keywords “intervertebral Disc” and “Intervertebral Disc” and limited the organism to Homo sapiens. By reviewing the summaries, overall designs, and sample sources of these datasets manually, we required that the sample source should be NPC. And in the overall designs, we also required that there were no drug therapy or other interventions in the population to ensure that we could use the selected datasets to screen out the differentially expressed molecules of degenerated NPC.

2.2. Data Analysis

Firstly, the differential expression of each dataset was analyzed by using the GEO2R tool [21]. After analysis, the differentially expressed genes (DEGs), differentially expressed miRNAs (DEMs), differentially expressed lncRNAs (DELs), and differentially expressed circRNAs (DECs) between the patients and the control group were obtained. The statistical significance was set to , adj. value < 0.05 (raw value was selected if there were too few molecules screened under the condition of adj. value < 0.05). Then, the series matrix files of these datasets were downloaded to do heat map analysis by Morpheus (https://software.broadinstitute.org/morpheus). Intersection analysis was accomplished by the Venn drawing tool of Bioinformatics & Evolutionary Genomics (http://bioinformatics.psb.ugent.be/webtools/Venn/) to gain common differentially expressed RNAs. Besides, functional enrichment analyses on the common DEGs, including Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, were done by the Database for Annotation, Visualization, and Integrated Discovery (DAVID) v6.8 tool [22, 23].

2.3. Protein–Protein Interaction (PPI) Network Construction

In the construction of PPI relationship, we used STRING as the prediction tool. STRING is a known database for predicting protein–protein interactions, including direct physical and indirect functional connections. The interaction relationship is derived from computer prediction, knowledge transfer between organisms, and interactions aggregated from other databases [24]. STRING extract experimental data from numbers of databases like BIND, DIP, and GRID and extract curated data from Biocarta, BioCyc, Reactome, and so on, which guarantee the reliability of mRNA–mRNA relationships in the interaction network. We mapped the DEGs to STRING and validated experimentally only interactions with a combined , which was selected as significant. Hub genes could be found through PPI analysis of common DEGs with the STRING tool, and the PPI network of DEGs was constructed and visualized using Cytoscape software (version 3.6.1; https://www.cytoscape.org.).

2.4. miRNA–mRNA Relationship Construction

miRNA-mRNA interactions were established by TargetScan Human 7.2, which is a tool that predicts biological targets of miRNAs by searching for the presence of conserved 8mer, 7mer, and 6mer sites that match the seed region of each miRNA [25]. The target mRNAs associated with DEMs could be predicted using TargetScan Human 7.2. We entered the gene names of each DEMs at TargetScan Human 7.2 to search for their predictive binding mRNAs, and among the predicted results, we selected those with low total context score and high aggregate PCT which suggest high possibility. Then, the target mRNAs were overlapped with the origin DEGs to find the very DEGs that are paired with DEMs and construct the miRNA–mRNA relationship using Cytoscape software.

2.5. lncRNA–miRNA Relationship Construction

lncRNA–miRNA relationships were constructed by LncBase Predicted v.2, a tool developed by DIANA Lab which presents an extensive collection of miRNA–lncRNA interactions. This significantly enhanced database includes more than 70000 low- and high-throughput, (in) direct miRNA–lncRNA experimentally supported interactions, derived from manually curated publications and the analysis of 153 AGO CLIP-Seq libraries. LncBase hosts in silico-predicted miRNA targets on lncRNA, identified with the DIANA-microT algorithm [26]. We used LncBase to predict the target lncRNAs related to DEMs, and then, the target lncRNAs were overlapped with the origin DELs to screen the DELs that are paired with DEMs and to construct the lncRNA–miRNA network using Cytoscape software.

2.6. circRNA–miRNA Relationship Construction

We used two tools, circBase and RegRNA 2.0, to help construct the circRNA–miRNA relationship. circBase can explore public circRNA datasets and export FASTA files containing genomic sequences [27], while RegRNA 2.0 is a tool that could comprehensively identify the functional RNA motifs and sites in an input RNA sequence, as well as predicting miRNAs that are related to the input circRNA [28]. DEC sequences were obtained from circBase, and then, RegRNA 2.0 was used to predict miRNAs that are related to DECs. After these miRNAs were overlapped with DEMs, DECs paired with DEMs was screened to construct the circRNA–miRNA relationship. The constructed lncRNA–miRNA and circRNA–miRNA relationships were then correlated with the miRNA–mRNA network, and Cytoscape software was used to build the lncRNA/circRNA–miRNA–mRNA network.

3. Results

After searching, we obtained a total of 6 mRNA datasets (GSE34095, GSE70362, GSE112216, GSE114169, GSE118927, and GSE124272), 3 miRNA datasets (GSE19943, GSE63492, and GSE116726), 3 lncRNA datasets (GSE56081, GSE112216, and GSE118927), and 1 circRNA dataset, which is GSE67566. The basic information of each dataset is shown in Table 1. Then, by reviewing the summaries, overall designs, and sample sources of these datasets manually, we required that the samples of the selected data sets should be NPC and there were no drug therapy or other interventions in the population in the overall designs.

It is worth noting that Lan et al. considered that the discs in scoliosis were not strictly normal discs [29], so they reused the normal discs as the control and replaced GSE19943 with GSE63492. Therefore, we also excluded the inclusion of GSE34095. Finally, we selected GSE70362 as the mRNA dataset, GSE63492 and GSE116726 as the miRNA dataset, GSE56081 as the lncRNA dataset, and GSE67566 as the circRNA dataset. According to the preset threshold (adj. value/raw value < 0.05 and ), we obtained 80 DEGs from GSE70362 (including 32 upregulated and 48 downregulated), 95 DEMs from GSE63492 (including 57 upregulated and 38 downregulated), and 983 DEMs from GSE116726 (including 527 upregulated and 456 downregulated). Besides, 115 DELs were screened out from GSE56081 (containing 50 upregulated and 65 downregulated) and 636 DECs were gained from GSE67566 (consisting of 354 upregulated and 282 downregulated). The top 20 DEGs, DEMs, DELs, and DECs are presented in Tables 2–4. The hierarchical cluster heat maps (Figure 2(a)) indicate that the expression differences of these RNAs were obvious between IDD NPC and normal NPC. Using the Venn drawing tool, the intersection analysis of the two miRNA datasets was performed and 33 DEMs were obtained, as shown in Figure 2(b). The respective expressions of these 33 miRNAs in the two datasets are shown in Figure 2(c), and 14 DEMs with consistent expression were ultimately selected, among which 9 were upregulated and 5 were downregulated. The 80 DEGs were analyzed using DAVID v6.8 to predict their possible biological functions and mechanism pathways. GO enrichment analysis shows 37 terms in total, and Table 5 exhibits 20 of them in which values are under 0.05 (as the values of KEGG analysis results were all greater than 0.05, they were not shown in the table). These DEGs may be involved in biological functions including extracellular exosomes, extracellular space, cellular response to zinc ions, negative regulation of growth, virus receptor activity, and extracellular region. The intuitive display of the enrichment analysis results is shown in Figure 3.

(a)

(b)

(c)

DEM-binding target genes predicted by TargetScan Human 7.2 were compared with 80 DEGs for overlapping. Thus, 10 DEMs and 49 DEGs formed a miRNA–mRNA interaction relationship and GO enrichment analysis was also performed on these 49 DEGs. In Table 5, these DEGs may be involved in the biological processes of the extracellular region, extracellular exosome, transcriptional activator activity, cell inhibiting, transcription factor binding, etc.

The direct display of the enrichment analysis results is shown in Figure 3. A PPI network consisting of 49 nodes and 21 interaction pairs was built using the STRING tool. Among these DEGs, CYP1B1, NQO1, FOXF2, FOXQ1, and GATA6 form an interaction network together, COL10A1, IBSP, and DLX3 form an interaction axis, and CD1D and CD80 form another interaction axis. Figure 4 shows the protein-protein interaction (PPI) network built by Cytoscape software. Using LncBase Predicted v.2, a total of 5973 lncRNAs were predicted by 10 DEMs, and after eliminating the duplicate data, they were superimposed with 3073 DELs. Thus, the lncRNA–miRNA interaction network between 30 DELs and 8 DEMs was built (11 DELs with downregulated expression regulate 5 DEMs with upregulated expression, while 19 DELs with upregulated expression regulate 3 DEMs).

Through circBase and RegRNA 2.0 database, 505 miRNAs predicted to interact with DECs were obtained to overlap with 10 DEMs and eventually 4 down-expressed DECs and 3 up-expressed DEMs form an interaction relationship. Table 6 exhibits the lncRNA/circRNA–miRNA–mRNA relationship and Figure 5 shows the interaction network constructed by Cytoscape software.

4. Discussion

As a ubiquitous health problem that affects global health and socioeconomics, it is critical to further study the biological mechanisms behind IDD. Except for the potential IDD mRNA, noncoding RNAs including miRNA, lncRNA, and circRNA may also be the targets used to develop potential therapy strategies. We have reviewed some existing research on the mechanisms of IDD ceRNA. Wang et al. suggested that lncRNA TRPC7-AS1 regulates nucleus pulposus cellular senescence and ECM synthesis via competing with HPN for miR-4769-5p binding [30], and Yang et al. found that lncRNA SLC20A1 promotes extracellular matrix degradation in NP cells by targeting the miR-31-5p/MMP3 axis [31]. In addition, a research suggested that LINC00969 promotes the degeneration of the intervertebral disc by sponging miR-335-3p and regulating NLRP3 inflammasome activation [32]. Except for these researches on lncRNA, there are also literatures reporting the mechanism of circRNA, like circularRNA_104670 can directly bind to miR-17-3p and correct the negative regulation of miR-17-3p on MMP-2, thus inhibiting the apoptosis of NP cells [33]. And hsa_circ_4099 can act as a “sponge” by competitive binding of miR-616-5p, thus reversing the inhibitory effect of miR-616-5p on Sox9 [34]. Noncoding RNA plays an important role in the process of IDD. However, there is still a lack of a comprehensive understanding of the molecular interactions involved in IDD, as there may be other important functional regulation axes. A broad regulatory network of the disease is in need. Based on the existing data from public databases, our study integrated and analyzed the genes and ncRNAs and formed a molecular interaction network. The network we built consists of 49 DEGs, 10 DEMs, 30 DELs, and 4 DECs, as well as containing 375 potential lncRNA/circRNA–miRNA–mRNA axes. We look forward to finding the genes and ncRNAs that are most likely related to IDD through this network, in order to lay a foundation for further research in the future.

GO enrichment analysis results show that the metabolism of hyaluronan may be one of the important biological processes for IDD. The hyaluronan synthase activity, hyaluronan biosynthetic process, and hyaluronan catabolic process are shown repeatedly in the GO analysis result list. Hyaluronic acid (HA) is one of the major glycosaminoglycan components of the extracellular matrix and is thought to be involved in cell proliferation, migration, and differentiation. HA plays a crucial role in retaining water of the intervertebral disc, so it can provide flexibility and shock absorbance in the spine [35]. Seen from the GO enrichment analysis results, HYAL1, HAS2, and LYVE1 are common gene molecules in hyaluronan-related metabolism. We searched these three genes in the constructed PPI interaction network and found that they constitute the interaction axis of HYAL1–HAS2–LYVE1. These findings suggest that they may indeed participate in the metabolism of hyaluronan and play a common regulatory role in IDD. By reviewing literatures, we found that HAS2 is a member of the gene family that encodes HA synthase, which was elevated in a study of HA-based hydrogel-induced NPC amplification [36]. HYAL1 encodes lysosomal hyaluronidase, which degrades HA in cells. In one study, the gene expression of HYAL1 and the protein expression of HAYL2 were significantly increased in moderate/severe IDD samples compared with those without or with low IDD [37]. The protein encoded by LYVE1 acts as a receptor and binds to both soluble and fixed HA and may also play a role in lymphatic HA transport. It is also been reported in the literature that ingrowth of vascularized fibrous tissue was seen at the edge or within fragments of degenerate disc tissue. Some scattered small vessels were lined by LYVE1+ endothelial cell [38]. In Reactome (https://www.reactome.org), we found a pathway through which HAS2, HYAL1, and LYVE1 coparticipate in HA biological processes (Figure 6). This pathway indicates that the integral membrane dual-action glycosyltransferase proteins hyaluronan synthases 1–3 (HAS1-3) mediate the polymerization of glucuronic acid (GlcA) with N-acetylglucosamine (GlcNAc) to form HA. The resulting polymer has the arrangement [-4GlcA-1,3GlcNAc-]n [39–41]. LYVE1, HMMR, STAB2, and CD44 together form the HAR complex. High molecular weight HA is tethered to the cell surface by HA receptors and the GPI-linked hyaluronidase 2 (HYAL2) to form a HA : HAR : HYAL2 complex in the plasma membrane that localizes to caveolae, continuing to complete the degradation of HA [42–45]. In the acidic environment of the lysosome, hyaluronidase 1 (HYAL1) is able to hydrolyze large 50 disaccharide unit HA fragments to 2 disaccharide units [46]. These theories suggest that the prediction of interaction relationship in this study has certain effectiveness.

Figure 6

Schematic diagram of hyaluronic acid-related IDD molecular mechanism. The integral membrane dual-action glycosyltransferase proteins hyaluronan synthases 1–3 (HAS1-3) mediate the polymerization of glucuronic acid (GlcA) with N-acetylglucosamine (GlcNAc) to form HA. The resulting polymer has the arrangement [-4GlcA-1,3GlcNAc-]n. LYVE1, HMMR, STAB2, and CD44 together form the HAR complex. High molecular weight HA is tethered to the cell surface by HA receptors and the GPI-linked hyaluronidase 2 (HYAL2) to form a HA : HAR : HYAL2 complex in the plasma membrane that localizes to caveolae, continuing to complete the degradation of HA. In the acidic environment of the lysosome, hyaluronidase 1 (HYAL1) could hydrolyze large 50 disaccharide unit HA fragments to 2 disaccharide units. HAS2, HYAL1, and LYVE1 may be negatively regulated by miR-4741, while mi663ahg/HEIH (lncRNA) and hsa-circ-0003600 (circRNA) may competitively be bound on miR-4741 to act as ceRNAs, correcting the negative regulation of miR-4741 to prevent IDD.

In the interaction network constructed in this study, HAS2, HYAL1, and LYVE1 may all be negatively regulated by miR-4741. Meanwhile, mi663ahg/HEIH (lncRNA) and hsa-circ-0003600 (circRNA) may competitively be bound on miR-4741 to act as ceRNAs, correcting the negative regulation of miR-4741. Through literature searching, we found that the mechanism of miR663AHG/HEIH/hsa-circ-0003600 is still lacking in the study of IDD. The interaction axis of miR663AHG/HEIH/hsa-circ-0003600–miR-4741–HAS2/HYAL1/LYVE1 may be one of the potential mechanisms to be studied in IDD.

We also noticed that COL10A1 was the gene with the greatest difference in expression between degenerative NP and normal NP (Table 2), and the results show that it had an aberrant expression in degenerative NP, which is associated with some existing research conclusions. COL10A1 encodes type X collagen, which is specifically expressed and synthesized by hypertrophic chondrocytes, and it plays an important role in the process of cartilage ossification and may be related to matrix degradation, calcification, and vascular invasion [47]. Type X collagen can be localized in combination with advanced intervertebral disc degeneration, and the positive staining of type X collagen is more obvious in the intervertebral disc matrix at the late stage of IDD in elderly patients [48]. In some experiments on mice that have led to disc degeneration, abnormal expression of type X collagen, increased apoptosis, inappropriate intervertebral disc ossification, and elevated expression of COL10A1, Runx2, MMP-13, etc. were observed in degenerative NP [11, 49–51]. Our network analysis also has confirmed that COL10A plays an important role in the pathogenesis of IDD. This not only verifies the credibility of our analysis but also provides an additional theoretical basis for further study on COL10A.

In the lncRNA/circRNA–miRNA–mRNA interaction network constructed in this study, COL10A1 may be negatively regulated by miR-5100, while lncRNA LINC00475, SOX2-OT, NEAT1, SNHG12, GAS5, AGAP2-AS1, and TARID may act as ceRNA for competitively binding miR-5100 and thus correcting the negative regulation of miR-5100. Wang et al. verified the contribution of miR-5100 to osteoblast differentiation and mineralization through functional gain and loss experiments, but there is still a lack of studies on miR-5100 in hypertrophic chondrocytes, type X collagen, or IDD. For lncRNA that may negatively regulate miR-5100, the current studies of LINC00475, Sox2-OT, SNHG12, AGAP2-AS1, and TARID are focused on various types of tumours, such as glioma, lymphoma, ovarian cancer, and renal clear cell sarcoma [52–55], while Wang et al. suggested that the overexpression of lncRNA GAS5 might promote the apoptosis of NPC through downregulation of Bcl-2 and upregulation of caspase-3 [56]. It is worth noting that NEAT1 may be the ceRNA-regulating COL10A1. Ruan et al. suggest that NEAT1 has been demonstrated to participate in ECM remodelling and they found the that increased expression of NEAT1 decreased aggrecan and type II collagen levels and increased ADAMTS4 and MMP-13 levels [57]. Our analysis results showed that NEAT1 expression was increased in IDD patients (), further confirming that NEAT1 may play an important role in IDD. Some research has proved that a robust expression of type X collagen and MMP-13 mRNA and protein could occur by hypertrophic chondrocytes during IDD [51, 58]. This suggests that NEAT1 and COL10A1 are connected by type X collagen and MMP-13. The ceRNA network shows that both NEAT1 and COL10A1 can interact with miR-5100 to form the interaction axis. Thus, we build a schematic diagram of the IDD molecular mechanism related to NEAT1 and COL10A1 (Figure 7). Although it can be hypothesized that in IDD patients, the NEAT1–miR-5100–COL10A1 interaction axis may be unbalanced, leading to the occurrence of IDD; there is still a lack of evidence that NEAT1 can negatively regulate miR-5100 and thereby regulate COL10A1 to realize ECM remodelling. Further elucidation is needed to confirm this hypothesis in the future.

Figure 7

Schematic diagram of COL10A1-related IDD molecular mechanism. The red line represents positive regulation, the green line represents negative regulation, and the arrow represents the direction of regulation. The solid line represents the verified regulatory relationship, and the dashed line represents the regulatory hypothesis proposed in this study. COL10A1 may be negatively regulated by miR-5100, while lncRNA LINC00475, SOX2-OT, NEAT1, SNHG12, GAS5, AGAP2-AS1, and TARID may act as ceRNA for competitively binding miR-5100 and thus correcting the negative regulation of miR-5100, inducing IDD.

5. Conclusions

Through bioinformatics analysis of multiple datasets, we constructed a ceRNA interaction network of lncRNA/circRNA–miRNA–mRNA for IDD. The result of this study verifies that the process of cartilage ossification, matrix degradation, and calcification following the synthesis of type X collagen and the imbalance of hyaluronic acid metabolism are two explanations for the occurrence of IDD. NEAT1–miR-5100–COL10A1 and miR663AHG/HEIH/hsa-circ-0003600–miR-4741–HAS2/HYAL1/LYVE1 might be potential interaction axes of the molecular mechanism in IDD. The increased expression of NEAT1 might inhibit miR-5100 and subsequently upregulate the expression of COL10A1, which leads to IDD, while the increased expression of miR663AHG/HEIH/hsa-circ-0003600 might inhibit miR-4741 and indirectly upregulate HAS2/HYAL1/LYVE1, and leads to the protection from IDD. This novel ceRNA network could be used to predict potential mechanisms to reach a comprehensive understanding and could provide possible targets for the further study of IDD in the future.

Data Availability

All data generated or analyzed during this study are included in this published article. The datasets analyzed during the current study are available in the Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/).

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Authors’ Contributions

Junshen Huang and Yuxi Li contributed equally to this work.

Acknowledgments

The authors thank Prof. Qin for idea guidance and thank Prof. Saw for providing English language editing. This research was funded by the Natural Science Foundation of Guangdong Province, grant number 2020A1515010371, 2017A030310554, and 2017A030313652, and the APC was funded by the China Postdoctoral Science Foundation, grant number 2020TQ0376.

References

M. Gore, A. Sadosky, B. R. Stacey, K.-S. Tai, and D. Leslie, “The burden of chronic low back pain: clinical comorbidities, treatment patterns, and health care costs in usual care settings,” Spine, vol. 37, no. 11, pp. E668–E677, 2012.
View at: Publisher Site | Google Scholar
P. B. Polatin, R. K. Kinney, R. J. Gatchel, E. Lillo, and T. G. Mayer, “Psychiatric illness and chronic low-back pain. The mind and the spine--which goes first?” Spine, vol. 18, no. 1, pp. 66–71, 1993.
View at: Publisher Site | Google Scholar
X. Luo, R. Pietrobon, S. X. Sun, G. G. Liu, and L. Hey, “Estimates and patterns of direct health care expenditures among individuals with back pain in the United States,” Spine, vol. 29, no. 1, pp. 79–86, 2004.
View at: Publisher Site | Google Scholar
G. K. Lutz, M. Butzlaff, and U. Schultz-Venrath, “Looking back on back pain: trial and error of diagnoses in the 20th century,” Spine, vol. 28, no. 16, pp. 1899–1905, 2003.
View at: Publisher Site | Google Scholar
K. M. C. Cheung, J. Karppinen, D. Chan et al., “Prevalence and pattern of lumbar magnetic resonance imaging changes in a population study of one thousand forty-three individuals,” Spine, vol. 34, no. 9, pp. 934–940, 2009.
View at: Publisher Site | Google Scholar
G. B. Andersson, “Epidemiological features of chronic low-back pain,” The Lancet (London, England), vol. 354, no. 9178, pp. 581–585, 1999.
View at: Publisher Site | Google Scholar
H. E. Gruber and E. N. Hanley, “Do we need biomarkers for disc degeneration?” Biomarker Insights, vol. 1, pp. 131–133, 2007.
View at: Google Scholar
D. Samartzis, J. Karppinen, D. Chan, K. D. K. Luk, and K. M. C. Cheung, “The association of lumbar intervertebral disc degeneration on magnetic resonance imaging with body mass index in overweight and obese adults: a population-based study,” Arthritis and Rheumatism, vol. 64, no. 5, pp. 1488–1496, 2012.
View at: Publisher Site | Google Scholar
S. Okie, “A flood of opioids, a rising tide of deaths,” The New England Journal of Medicine, vol. 363, no. 21, pp. 1981–1985, 2010.
View at: Publisher Site | Google Scholar
P. H. Wu, H. S. Kim, and I.-T. Jang, “Intervertebral disc diseases part 2: a review of the current diagnostic and treatment strategies for intervertebral disc disease,” International Journal of Molecular Sciences, vol. 21, no. 6, p. 2135, 2020.
View at: Publisher Site | Google Scholar
H. Choi, S. Tessier, E. S. Silagi et al., “A novel mouse model of intervertebral disc degeneration shows altered cell fate and matrix homeostasis,” Matrix Biology, vol. 70, pp. 102–122, 2018.
View at: Publisher Site | Google Scholar
S. Sobajima, J. F. Kompel, J. S. Kim et al., “A slowly progressive and reproducible animal model of intervertebral disc degeneration characterized by mri, x-ray, and histology,” Spine, vol. 30, no. 1, pp. 15–24, 2005.
View at: Publisher Site | Google Scholar
L. Frapin, J. Clouet, V. Delplace, M. Fusellier, J. Guicheux, and C. Le Visage, “Lessons learned from intervertebral disc pathophysiology to guide rational design of sequential delivery systems for therapeutic biological factors,” Advanced Drug Delivery Reviews, vol. 149-150, pp. 49–71, 2019.
View at: Publisher Site | Google Scholar
R. V. T. Filho, G. de Matos Abe, and M. T. Daher, “Genetic influence in disc degeneration - systematic review of literature,” Revista Brasileira de Ortopedia, vol. 55, no. 2, pp. 131–138, 2020.
View at: Publisher Site | Google Scholar
L. Kalichman and D. J. Hunter, “The genetics of intervertebral disc degeneration.Associated genes,” Joint, Bone, Spine, vol. 75, no. 4, pp. 388–396, 2008.
View at: Publisher Site | Google Scholar
N. L. Martirosyan, A. A. Patel, A. Carotenuto et al., “Genetic alterations in intervertebral disc disease,” Frontiers in Surgery, vol. 3, p. 59, 2016.
View at: Publisher Site | Google Scholar
S. Kalb, N. L. Martirosyan, M. Y. S. Kalani, G. G. Broc, and N. Theodore, “Genetics of the degenerated intervertebral disc,” World Neurosurgery, vol. 77, no. 3-4, pp. 491–501, 2012.
View at: Publisher Site | Google Scholar
S. Hanaei, S. Abdollahzade, A. Khoshnevisan, C. K. Kepler, and N. Rezaei, “Genetic aspects of intervertebral disc degeneration,” Reviews in the Neurosciences, vol. 26, no. 5, pp. 581–606, 2015.
View at: Publisher Site | Google Scholar
J. Hausser and M. Zavolan, “Identification and consequences of miRNA-target interactions -- beyond repression of gene expression,” Nature Reviews. Genetics, vol. 15, no. 9, pp. 599–612, 2014.
View at: Publisher Site | Google Scholar
L. Salmena, L. Poliseno, Y. Tay, L. Kats, and P. P. Pandolfi, “A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language?” Cell, vol. 146, no. 3, pp. 353–358, 2011.
View at: Publisher Site | Google Scholar
S. Davis and P. S. Meltzer, “Geoquery: a bridge between the gene expression omnibus (geo) and bioconductor,” Bioinformatics (Oxford, England), vol. 23, no. 14, pp. 1846-1847, 2007.
View at: Publisher Site | Google Scholar
D. W. Huang, B. T. Sherman, and R. A. Lempicki, “Systematic and integrative analysis of large gene lists using David bioinformatics resources,” Nature Protocols, vol. 4, no. 1, pp. 44–57, 2009.
View at: Publisher Site | Google Scholar
D. W. Huang, B. T. Sherman, and R. A. Lempicki, “Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists,” Nucleic Acids Research, vol. 37, no. 1, pp. 1–13, 2009.
View at: Publisher Site | Google Scholar
D. Szklarczyk, A. L. Gable, D. Lyon et al., “String v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets,” Nucleic Acids Research, vol. 47, no. D1, pp. D607–D613, 2019.
View at: Publisher Site | Google Scholar
B. P. Lewis, C. B. Burge, and D. P. Bartel, “Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microrna targets,” Cell, vol. 120, no. 1, pp. 15–20, 2005.
View at: Publisher Site | Google Scholar
M. D. Paraskevopoulou, I. S. Vlachos, D. Karagkouni et al., “Diana-lncbase v2: indexing microrna targets on non-coding transcripts,” Nucleic Acids Research, vol. 44, no. D1, pp. D231–D238, 2016.
View at: Publisher Site | Google Scholar
P. Glažar, P. Papavasileiou, and N. Rajewsky, “Circbase: a database for circular rnas,” RNA (New York, N.Y.), vol. 20, no. 11, pp. 1666–1670, 2014.
View at: Publisher Site | Google Scholar
T.-H. Chang, H.-Y. Huang, J. B.-K. Hsu, S.-L. Weng, J.-T. Horng, and H.-D. Huang, “An enhanced computational platform for investigating the roles of regulatory rna and for identifying functional rna motifs,” BMC Bioinformatics, vol. 14, Supplement 2, 2013.
View at: Publisher Site | Google Scholar
P.-H. Lan, Z.-H. Liu, Y.-J. Pei et al., “Landscape of rnas in human lumbar disc degeneration,” Oncotarget, vol. 7, no. 39, pp. 63166–63176, 2016.
View at: Publisher Site | Google Scholar
X. Wang, D. Li, H. Wu et al., “Lncrna trpc7-as1 regulates nucleus pulposus cellular senescence and ecm synthesis via competing with hpn for mir-4769-5p binding,” Mechanisms of Ageing and Development, vol. 190, p. 111293, 2020.
View at: Publisher Site | Google Scholar
Y. Yang, Z. Zhong, Y. Zhao, K. Ren, and N. Li, “Lincrna-slc20a1 (slc20a1) promotes extracellular matrix degradation in nucleus pulposus cells in human intervertebral disc degeneration by targeting the mir-31-5p/mmp3 axis,” International Journal of Clinical and Experimental Pathology, vol. 12, no. 9, pp. 3632–3643, 2019.
View at: Google Scholar
L. Yu, Y. Hao, C. Xu, G. Zhu, and Y. Cai, “Linc00969 promotes the degeneration of intervertebral disk by sponging mir-335-3p and regulating nlrp3 inflammasome activation,” IUBMB Life, vol. 71, no. 5, pp. 611–618, 2019.
View at: Publisher Site | Google Scholar
J. Song, H.-L. Wang, K.-H. Song et al., “Circularrna_104670 plays a critical role in intervertebral disc degeneration by functioning as a cerna,” Experimental & Molecular Medicine, vol. 50, no. 8, pp. 1–12, 2018.
View at: Publisher Site | Google Scholar
H. Wang, P. He, H. Pan et al., “Circular rna circ-4099 is induced by tnf-α and regulates ecm synthesis by blocking mir-616-5p inhibition of sox9 in intervertebral disc degeneration,” Experimental & Molecular Medicine, vol. 50, no. 4, pp. 1–14, 2018.
View at: Publisher Site | Google Scholar
K. Benz, C. Stippich, C. Osswald et al., “Rheological and biological properties of a hydrogel support for cells intended for intervertebral disc repair,” BMC Musculoskeletal Disorders, vol. 13, no. 1, 2012.
View at: Publisher Site | Google Scholar
A. Peroglio, S. Grad, D. Mortisen et al., “Injectable thermoreversible hyaluronan-based hydrogels for nucleus pulposus cell encapsulation,” European Spine Journal, vol. 21, no. S6, pp. 839–849, 2011.
View at: Publisher Site | Google Scholar
O. Krupkova, H. Greutert, N. Boos, J. Lemcke, T. Liebscher, and K. Wuertz-Kozak, “Expression and activity of hyaluronidases HYAL-₁, HYAL-₂ and HYAL-₃ in the human intervertebral disc,” European Spine Journal, vol. 29, no. 3, pp. 605–615, 2019.
View at: Publisher Site | Google Scholar
K. Benz, C. Stippich, C. Osswald et al., “Rheological and biological properties of a hydrogel support for cells intended for intervertebral disc repair,” BMC Musculoskeletal Disorders, vol. 13, no. 1, p. 54, 2012.
View at: Publisher Site | Google Scholar
K. Watanabe and Y. Yamaguchi, “Molecular identification of a putative human hyaluronan synthase,” The Journal of Biological Chemistry, vol. 271, no. 38, pp. 22945–22948, 1996.
View at: Publisher Site | Google Scholar
A. P. Spicer, J. S. Olson, and J. A. McDonald, “Molecular cloning and characterization of a cDNA encoding the third putative mammalian hyaluronan synthase,” The Journal of Biological Chemistry, vol. 272, no. 14, pp. 8957–8961, 1997.
View at: Publisher Site | Google Scholar
A. M. Shyjan, P. Heldin, E. C. Butcher, T. Yoshino, and M. J. Briskin, “Functional cloning of the cDNA for a human hyaluronan synthase,” The Journal of Biological Chemistry, vol. 271, no. 38, pp. 23395–23399, 1996.
View at: Publisher Site | Google Scholar
E. N. Harris, S. V. Kyosseva, J. A. Weigel, and P. H. Weigel, “Expression, processing, and glycosaminoglycan binding activity of the recombinant human 315-kDa hyaluronic acid receptor for endocytosis (HARE),” The Journal of Biological Chemistry, vol. 282, no. 5, pp. 2785–2797, 2007.
View at: Publisher Site | Google Scholar
S. Banerji, J. Ni, S. X. Wang et al., “Lyve-1, a new homologue of the cd44 glycoprotein, is a lymph-specific receptor for hyaluronan,” The Journal of Cell Biology, vol. 144, no. 4, pp. 789–801, 1999.
View at: Publisher Site | Google Scholar
R. Asher and A. Bignami, “Hyaluronate binding and cd44 expression in human glioblastoma cells and astrocytes,” Experimental Cell Research, vol. 203, no. 1, pp. 80–90, 1992.
View at: Publisher Site | Google Scholar
C. Wang, J. Entwistle, G. Hou, Q. Li, and E. A. Turley, “The characterization of a human RHAMM cDNA: conservation of the hyaluronan- binding domains,” Gene, vol. 174, no. 2, pp. 299–306, 1996.
View at: Publisher Site | Google Scholar
M. J. Jedrzejas and R. Stern, “Structures of vertebrate hyaluronidases and their unique enzymatic mechanism of hydrolysis,” Proteins, vol. 61, no. 2, pp. 227–238, 2005.
View at: Publisher Site | Google Scholar
G. Shen, “The role of type x collagen in facilitating and regulating endochondral ossification of articular cartilage,” Orthodontics & Craniofacial Research, vol. 8, no. 1, pp. 11–17, 2005.
View at: Publisher Site | Google Scholar
N. Boos, A. G. Nerlich, I. Wiest, K. von der Mark, and M. Aebi, “Immunolocalization of type x collagen in human lumbar intervertebral discs during ageing and degeneration,” Histochemistry and Cell Biology, vol. 108, no. 6, pp. 471–480, 1997.
View at: Publisher Site | Google Scholar
Y.-J. Lao, T.-T. Xu, H.-T. Jin et al., “Accumulated spinal axial biomechanical loading induces degeneration in intervertebral disc of mice lumbar spine,” Orthopaedic Surgery, vol. 10, no. 1, pp. 56–63, 2018.
View at: Publisher Site | Google Scholar
Q.-Q. Liang, X.-J. Cui, Z.-J. Xi et al., “Prolonged upright posture induces degenerative changes in intervertebral discs of rat cervical spine,” Spine, vol. 36, no. 1, pp. E14–E19, 2011.
View at: Publisher Site | Google Scholar
J. Zieba, K. N. Forlenza, J. S. Khatra et al., “Tgfβ and bmp dependent cell fate changes due to loss of filamin b produces disc degeneration and progressive vertebral fusions,” PLoS Genetics, vol. 12, no. 3, article e1005936, 2016.
View at: Publisher Site | Google Scholar
L. Yu, S. Gui, Y. Liu et al., “Long intergenic non-protein coding rna 00475 silencing acts as a tumor suppressor in glioma under hypoxic condition by impairing microrna-449b-5p-dependent agap2 up-regulation,” Therapeutic Advances in Medical Oncology, vol. 12, 2020.
View at: Publisher Site | Google Scholar
L.-Y. Chen, X.-M. Zhang, B.-Q. Han, and H.-B. Dai, “Long noncoding rna snhg12 indicates the prognosis and accelerates tumorigenesis of diffuse large b-cell lymphoma through sponging micror-195,” Oncotargets and Therapy, vol. 13, pp. 5563–5574, 2020.
View at: Publisher Site | Google Scholar
Z. Tingting, L. Xiaojing, T. Xiaoyan, H. Keqin, and Q. Junjun, “The antisense long noncoding rna agap2-as1 regulates cell proliferation and metastasis in epithelial ovarian cancer,” Journal of Cancer, vol. 11, no. 18, pp. 5318–5328, 2020.
View at: Publisher Site | Google Scholar
S. L. Gooskens, S. Gadd, J. M. G. Auvil et al., “Tcf21 hypermethylation in genetically quiescent clear cell sarcoma of the kidney,” Oncotarget, vol. 6, no. 18, pp. 15828–15841, 2015.
View at: Publisher Site | Google Scholar
Y. Wang, Q. Song, X. Huang et al., “Long noncoding rna gas5 promotes apoptosis in primary nucleus pulposus cells derived from the human intervertebral disc via bcl-2 downregulation and caspase-3 upregulation,” Molecular Medicine Reports, vol. 19, pp. 2164–2172, 2019.
View at: Publisher Site | Google Scholar
Z. Ruan, H. Ma, J. Li, H. Liu, H. Jia, and F. Li, “The long non-coding rna neat1 contributes to extracellular matrix degradation in degenerative human nucleus pulposus cells,” Experimental Biology and Medicine (Maywood, N.J.), vol. 243, pp. 595–600, 2018.
View at: Publisher Site | Google Scholar
H. Jin, J. Shen, B. Wang, M. Wang, B. Shu, and D. Chen, “Tgf-β signaling plays an essential role in the growth and maintenance of intervertebral disc tissue,” FEBS Letters, vol. 585, no. 8, pp. 1209–1215, 2011.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Junshen Huang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

383

Downloads

1066

Citations