Dataset Information

Major ceRNA regulation and key metabolic signature analysis of intervertebral disc degeneration.

ABSTRACT:

Background and objective

Intervertebral disc degeneration (IDD) is a complex multifactorial and irreversible pathological process. In IDD, multiple competing endogenous RNAs (ceRNA, including mRNA, lncRNA, and pseudogenes) can compete to bind with miRNAs. However, the potential metabolic signatures in nucleus pulposus (NP) cells remain poorly understood. This study investigated key metabolic genes and the ceRNA regulatory mechanisms in the pathogenesis of IDD based on microarray datasets.

Methods

We retrieved and downloaded four independent IDD microarray datasets from the Gene Expression Omnibus. Combining the predicted interactions from online databases (miRcode, miRDB, miRTarBase, and TargetScan), differentially expressed lncRNAs (DElncRNAs), miRNAs (DEmiRNAs), and mRNAs (DEmRNAs) were identified. A ceRNA network was constructed and annotated using GO and KEGG pathway enrichment analyses. Moreover, we searched the online metabolic gene set and used support vector machine (SVM) to find the critical metabolic DEmRNA(s) and other DERNAs. Differential gene expression was validated with a merged dataset.

Results

A total of 45 DEmRNAs, 36 DElncRNAs, and only one DEmiRNA (miR-338-3p) were identified in the IDD microarray datasets. GO and KEGG pathway enrichment analyses revealed that the DEmRNAs were predominantly enriched in the PI3K-Akt signaling pathway, MAPK signaling pathway, IL-17 signaling pathway, apoptosis, and cellular response to oxidative stress. Based on SVM screening, 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase (PFK/FBPase) 2 is the critical metabolic gene with lower expression in IDD, and AC063977.6 is the key lncRNA with lower expression in IDD. The ceRNA hypothesis suggests that AC063977.6, miR-338-3p (high expression), and PFKFB2 are dysregulated as an axis in IDD.

Conclusions

The results suggest that lncRNA AC063977.6 correlate with PFKFB2, the vital metabolic signature gene, via targeting miR-338-3p during IDD pathogenesis. The current study may shed light on unraveling the pathogenesis of IDD.

SUBMITTER: Cao S

PROVIDER: S-EPMC7937257 | biostudies-literature | 2021 Mar

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Major ceRNA regulation and key metabolic signature analysis of intervertebral disc degeneration.

Cao Shuai S Li Jie J Yang Kai K Li Haopeng H

BMC musculoskeletal disorders 20210306 1

<h4>Background and objective</h4>Intervertebral disc degeneration (IDD) is a complex multifactorial and irreversible pathological process. In IDD, multiple competing endogenous RNAs (ceRNA, including mRNA, lncRNA, and pseudogenes) can compete to bind with miRNAs. However, the potential metabolic signatures in nucleus pulposus (NP) cells remain poorly understood. This study investigated key metabolic genes and the ceRNA regulatory mechanisms in the pathogenesis of IDD based on microarray datasets ...[more]

PMID: 33676464

Similar Datasets

Project description:Accumulating evidence has indicated that noncoding RNAs are involved in intervertebral disc degeneration (IDD); however, the competing endogenous RNA (ceRNA)‑mediated regulatory mechanisms in IDD remain rarely reported. The present study aimed to comprehensively investigate the alterations in expression levels of circular RNA (circRNA), long noncoding RNA (lncRNA), microRNA (miRNA/miR) and mRNA in the nucleus pulposus (NP) of patients with IDD. In addition, crucial lncRNA/circRNA‑miRNA‑mRNA ceRNA interaction axes were screened using the GSE67567 microarray dataset obtained from the Gene Expression Omnibus database. After data preprocessing, differentially expressed circRNAs (DECs), lncRNAs (DELs), miRNAs (DEMs) or genes (DEGs) between IDD and normal controls were identified using the Linear Models for Microarray data method. A protein‑protein interaction (PPI) network was constructed for DEGs based on protein databases, followed by module analysis. The ceRNA network was constructed based on the interaction between miRNAs and mRNAs, and lncRNAs/circRNAs and miRNAs. The underlying functions of mRNAs were predicted using the Database for Annotation, Visualization and Integrated Discovery database. The present study identified 636 DECs, 115 DELs, 84 DEMs and 1,040 DEGs between patients with IDD and control individuals. PPI network analysis demonstrated that Fos proto‑oncogene, AP‑1 transcription factor subunit (FOS), mitogen‑activated protein kinase 1 (MAPK1), hypoxia inducible factor 1 subunit α (HIF1A) and transforming growth factor β1 (TGFB1) were hub genes and enriched in modules. Metastasis‑associated lung adenocarcinoma transcript 1 (MALAT1)/hsa_circRNA_102348‑hsa‑miR‑185‑5p‑TGFB1/FOS, MALAT1‑hsa‑miR‑155‑5p‑HIF1A, hsa_circRNA_102399‑hsa‑miR‑302a‑3p‑HIF1A, MALAT1‑hsa‑miR‑519d‑3p‑MAPK1 and hsa_circRNA_100086‑hsa‑miR‑509‑3p‑MAPK1 ceRNA axes were obtained by constructing the ceRNA networks. In conclusion, these identified ceRNA interaction axes may be crucial targets for the treatment of IDD.

Dataset Information

Major ceRNA regulation and key metabolic signature analysis of intervertebral disc degeneration.

Background and objective

Methods

Results

Conclusions

Publications

Major ceRNA regulation and key metabolic signature analysis of intervertebral disc degeneration.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets