Oxidative Stress in Intervertebral Disc Degeneration: New Insights from Bioinformatic Strategies

Zhao, Yongzhao; Xiang, Qian; Lin, Jialiang; Jiang, Shuai; Li, Weishi

doi:https://doi.org/10.1155/2022/2239770

Oxidative Medicine and Cellular Longevity

On this page

Abstract Background Materials and Methods Results Discussion Abbreviations Data Availability Ethical Approval Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Special Issue

Oxidative Stress in Intervertebral Disc Degeneration and its Related Therapeutics 2021

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 2239770 | https://doi.org/10.1155/2022/2239770

Oxidative Stress in Intervertebral Disc Degeneration: New Insights from Bioinformatic Strategies

Yongzhao Zhao,^1,2,3Qian Xiang,^1,2,3Jialiang Lin,^1,2,3Shuai Jiang,^1,2,3and Weishi Li^1,2,3

Academic Editor: Sidong Yang

Received06 Jan 2022

Accepted10 Mar 2022

Published31 Mar 2022

Abstract

Oxidative stress has been proved to play important roles in the development of intervertebral disc degeneration (IDD); however, the underlying mechanism remains obscure to date. The aim of this study was to elucidate the vital roles of oxidative stress-related genes in the development of IDD using strict bioinformatic algorithms. The microarray data relevant to the IDD was downloaded from Gene Expression Omnibus database for further analysis. A series of bioinformatic strategies were used to determine the oxidative stress-related and IDD-related genes (OSIDDRGs), perform the function enrichment analysis and protein-protein interaction analysis, construct the lncRNA-miRNA-mRNA regulatory network, and investigate the potential relationship of oxidative stress to immunity abnormality and autophagy in IDD. We observed a significantly different status of oxidative stress between normal intervertebral disc tissues and IDD tissues. A total of 72 OSIDDRGs were screened out for the further function enrichment analysis, and 10 hub OSIDDRGs were selected to construct the lncRNA-miRNA-mRNA regulatory network. There was a very close association of oxidative stress with immunity abnormality and autophagy in IDD. Taken together, our findings can provide new insights into the mechanism research of oxidative stress in the development of IDD and offer new potential targets for the treatment strategies.

1. Background

It is estimated that approximately 84% adults will suffer from the low back pain (LBP) during their lifetimes [1, 2]. Intervertebral disc (IVD) degeneration (IDD) is the principal contributor to the LBP, but the precise pathogenesis remains obscure [3]. Several mechanisms have been proved to be involved in the initiation or progression of IDD, such as mechanical stress, immunity abnormality, metabolic disorders, and oxidative stress [3–8]. Especially, oxidative stress has gotten increased attention in the pathogenesis of IDD recently [9–11].

Oxidative stress refers to an imbalance between oxidants and antioxidants in the cells and tissues, which will result in the accumulation of reactive oxygen species [6]. Recently, increasing evidence shows that oxidative stress may play vital roles in IDD by inducing the premature senescence and promoting a catabolic phenotype in human nucleus pulposus (NP) cells [9–14]. Hu et al. found that exosomes derived from bone mesenchymal stem cells could relive the compression-induced apoptosis of NP cells through inhibiting the oxidative stress [10]. Seol et al. reported that the reduction of oxidative stress by amobarbital was an effective way to prevent the IDD progression [13]. Despite a growing number of relevant studies, the underlying mechanism of oxidative stress in IDD has not been fully elucidated up to now.

Long noncoding RNA (lncRNA) refers to the noncoding RNA whose length exceeds 200 nucleotides, and it can interact with the DNA, RNA, and protein [15]. LncRNA is involved in many diseases via functioning as a competing endogenous RNA (ceRNA) to suppress the microRNA (miRNA) functions and upregulate the targeted gene expression [16–18]. However, limited research has been performed to elucidate the ceRNA network of lncRNA in the regulation of oxidative stress in IDD.

Immunity abnormality is a very crucial contributor to the development of IDD [3, 19]. The proinflammatory cytokines secreted by the NP cells and the infiltrating immune cells (IICs) are involved in the progression of IDD [5]. Autophagy is an important intracellular degradation process to maintain the homeostasis of intracellular environment by removal of damaged organelles and nonfunction proteins and cyclic utilization of degraded components [20]. The activation of autophagy may be an important protection mechanism to relieve the damage of oxidative stress and thereby prevent the progression of IDD [21, 22]. However, the underlying association of oxidative stress with immunity abnormality and autophagy in the development of IDD has not been fully understood.

In this study, we used a series of mature and recognized bioinformatic strategies to determine the following items: (1) The oxidative stress-related genes (OSRGs) involved in the pathogenesis of IDD, their main biological functions, and the potential ceRNA regulatory network; and (2) the possible relationship of oxidative stress to immunity abnormality and autophagy in IDD.

2. Materials and Methods

2.1. Data Collection and Processing

This research has been approved by the institutional research ethic committee of Peking University Third Hospital, and the informed consent was exempted because all data was obtained from the common database. The flow chart of this research was shown in Figure 1. Gene expression data of mRNA and miRNA was downloaded from Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/). The dataset used for this research was selected according to following criterions: (a) species: Homo sapiens; (b) sequencing methods: microarray data or next-generation sequencing data; and (c) samples: the nucleus pulposus in normal IVD tissues and IDD tissues. The following dataset was excluded from this research: expression data from other tissues (e.g., blood sample) and IVD tissues from subjects with scoliosis. Finally, the mRNA expression dataset GSE56081 and miRNA expression dataset GSE63492 were selected and downloaded for further analysis. The expression files of these two datasets were from 5 IDD tissues and 5 normal IVD tissues by the microarray, and the details of these tissue samples were listed in the Supplementary Table 1. The OSRG list was obtained from Gene Set GOBP_RESPONSE_TO_OXIDATIVE_STRESS in Molecular Signatures database (http://www.gsea-msigdb.org/gsea/msigdb/index.jsp) [23] (Supplementary Table 2). The immunity-related gene list was downloaded from the immunology database and Analysis Portal (ImmPort) database (https://www.immport.org/home) [24] (Supplementary Table 3). The autophagy-related gene list was downloaded from the Human Autophagy (HADb) database (http://www.autophagy.lu/index.html) (Supplementary Table 4) [25].

2.2. Determination of Oxidative Stress in IDD

To investigate the potential role of oxidative stress in the pathogenesis of IDD, the oxidative stress score for each tissue sample was calculated using the single sample gene set enrichment analysis (ssGSEA), which was a computational approach to explore whether a priori defined set of genes has statistical significance and concordant differences in two biological conditions for a single sample [26–28].

2.3. Identification of Oxidative Stress-Related and IDD-Related Genes (OSIDDRGs) and IDD-Related miRNAs (IDDRmiRNAs)

IDD-related genes were obtained using the R package limma package, with the criterion of adjusted and . The OSIDDRGs were obtained with the intersection of the IDD-related genes and OSRGs using the Venn diagram. IDDRmiRNAs were obtained using the R package limma package, with the standard of and . The volcano plot and heat map were plotted using the R package ggplot2.

2.4. Functional Enrichment Analysis of OSIDDRGs

Gene Ontology (GO) analysis of OSIDDRGs was conducted to interpret their biological process, cell component, and molecular function. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was performed to explore the related signaling pathways of OSIDDRGs. The GO and KEGG analyses were both conducted using the DAVID database (https://david.ncifcrf.gov/) [29], and specific items of GO and KEGG analyses with p <0.05 were selected for further visualization using the R package ggplot2.

2.5. Protein-Protein Interaction (PPI) Analysis of OSIDDRGs

The STRING database (https://cn.string-db.org/) was a widely used tool to explore the association networks of functional proteins [30]. PPI analysis was conducted using the STRING database by imputing the OSIDDRGs into the multiple protein section, and protein pairs with were selected to construct the PPI network by Cytoscape software (https://cytoscape.org/). The PPI score was calculated using the MCC method by the cytoHubba plug-in, and top 10 OSIDDRGs ranked by the PPI score were selected for further analysis.

2.6. Construction of ceRNA Regulatory Network of lncRNA-miRNA-mRNA

For each OSIDDRG, the correlation analysis between IDDRmiRNAs and hub OSIDDRG was performed using the Pearson test, and significantly related IDDRmiRNAs with and were selected for further intersection with predicted miRNAs. The predicted miRNAs for 10 hub OSIDDRGs were obtained using TargetScanHuman database (http://www.targetscan.org/vert_80/) [31]. The intersection of significantly related IDDRmiRNAs and predicted miRNAs was performed using the Venn diagram to obtain the candidate miRNAs for the construction of miRNA-mRNA pairs. Then, the upstream targeted lncRNAs for candidate miRNAs were predicted to construct the lncRNA-miRNA pairs using the Predicted model in the LncBase Predicted database (http://carolina.imis.athenainnovation.gr/diana_tools/web/index.php?r=lncbasev2/index-predicted) [32]. Ultimately, the ceRNA regulatory network of lncRNA-miRNA-mRNA was constructed and presented by the Cytoscape software.

2.7. Correlation Analysis between Oxidative Stress and Immunity Abnormality

The immunity score for each tissue sample was calculated using the ssGSEA algorithm based on the immunity-related gene list [5] and compared between normal IVD tissues and IDD tissues. The immune infiltration analysis was performed using the CIBERSORT database (http://cibersort.stanford.edu/) [33]. The difference of IICs was compared between normal IVD tissues and IDD tissues, and IICs with were considered as IDD-related IICs (IDDRIICs). The correlation analysis between oxidative stress score and immunity score and the correlation analysis between 10 hub OSIDDRGs and 8 types of IDDRIICs were conducted to explore the potential relationship between oxidative stress and immunity abnormality in IDD using the Pearson test. To further elucidate the interaction between oxidative stress and immunity abnormality in the development of IDD, the miRNA-mRNA-IDDRIIC regulatory network was built by integrating the IDDRmiRNA-OSIDDRG pairs with and with OSIDDRG-IDDRIIC pairs with and . Ultimately, the miRNA-mRNA-IDDRIIC regulatory network was constructed and presented using the Sankey diagram.

2.8. Correlation Analysis between Oxidative Stress and Autophagy

The autophagy score for each tissue sample was calculated using the ssGSEA algorithm based on the autophagy-related gene list [26] and compared between normal IVD tissues and IDD tissues. To obtain the autophagy-related and IDD-related genes (ATGIDDRGs), the intersection between autophagy-related genes and IDD-related genes was conducted using the Venn diagram. Then, the PPI analysis of ATGIDDRGs was conducted to obtain the top 10 hub ATGIDDRGs using the MCC method by Cytoscape software. Furthermore, the correlation analysis between hub OSIDDRGs and hub ATGIDDRGs was conducted.

2.9. Statistical Analysis

Any analysis in this study was conducted by the R software 4.1.2. The ssGSEA score and percentage of IICs between normal IVD tissues and IDD tissues were compared using the Student’s test and indicating there was a statistical difference between two groups. Correlation analysis in this study was performed using the Pearson test, and the correlated pair with and was considered as statistically significant.

3. Results

3.1. Confirmation of Differential Oxidative Stress Score between Normal IVD Tissues and IDD Tissues

The oxidative stress score was calculated by ssGSEA algorithm based on the expression of OSRGs, and there was a significant difference of oxidative stress score between normal IVD tissues and IDD tissues () (Figure 2(a)). PCA analysis showed the expression of OSRGs could clearly distinguish the normal IVD tissues and IDD tissues (Figure 2(b)). Therefore, oxidative stress might play vital roles in the pathophysiology of IDD.

(a)

(b)

3.2. Determination of IDD-Related Genes and OSIDDRGs

As shown in the volcano plot (Figure 3(a)), a total of 2,269 IDD-related genes were identified with the criterion of and adjusted . The heat map showed the cluster analysis of these IDD-related genes could clearly distinguish the normal IVD tissues and IDD tissues (Figure 3(b)). A total of 72 OSIDDRGs were obtained by the intersection between OSRGs and IDD-related genes (Figure 3(c)), and 47 of them were upregulated and 25 were downregulated (Figure 3(d)) (Supplementary Table 5).

(a)

(b)

(c)

(d)

3.3. Function Analysis and PPI Analysis of OSIDDRGs

The GO analysis, which consisted of biological process, cell component, and molecular function, was conducted to explore the biological functions of OSIDDRGs. The biological process of these OSIDDRGs was mainly enriched at the response to oxidative stress, oxidation-reduction process, positive regulation of transcription of DNA-templated, apoptotic process, response to reactive oxygen species, aging, response to hypoxia, angiogenesis, and cellular response to interleukin-1 (Figure 4(a)). These OSIDDRGs elicited the biological functions mainly at the mitochondrion (Figure 4(b)), and the most common molecular function of these OSIDDRGs was the protein binding (Figure 4(c)). KEGG analysis showed these OSIDDRGs were mainly involved in the following biological pathways: TNF signaling pathway, estrogen signaling pathway, MAPK signaling pathway, GnRH signaling pathway, HIF-1 signaling pathway, NOD-like receptor signaling pathway, FoxO signaling pathway, and Toll-like receptor signaling pathway (Figure 4(d)).

PPI analysis of 72 OSIDDRGs was conducted with a network of 77 nodes and 279 edges (PPI enrichment ) (Figure 4(e)). The top 10 hub OSIDDRGs were selected, including IL6, PRDX1, MCL1, HMOX1, TXNRD1, MAPK1, HIF1A, FOXO1, JUN, and JAK2 (Figure 4(f)). The correlation analysis among top 10 hub OSIDDRGs was conducted (Figure 4(g)), and JUN-FOXO1 pair was the most positively correlated pair and IL6-MAPK1 pair (, ) (Figure 4(h)) was the most negatively correlated pair (, ) (Figure 4(i)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(i)

3.4. Construction of lncRNA-miRNA-mRNA Regulatory Network

A total of 47 IDDRmiRNAs were obtained with the standard of and (Figure 5(a)), and the heat map showed there was a different expression level of IDDRmiRNAs between the normal IVD tissues and IDD tissues (Figure 5(b)). For each hub OSIDDRG, the correlation analysis between the OSIDDRGs and IDDRmiRNAs was constructed. The significantly related IDDRmiRNAs ( and ) were further used to take the intersection with the predicted miRNAs using TargetScan Human database to obtain the candidate IDDRmiRNAs (Supplementary Figure 1). Ultimately, 13 miRNA-mRNA pairs consisted of 9 IDDRmiRNAs and 6 hub OSIDDRGs were obtained (Figure 5(c)). Then, 63 targeted lncRNAs for above-mentioned 9 IDDRmiRNAs were predicted using the LncBase database, and 90 lncRNA-miRNA pairs were obtained (Supplementary Table 6). Ultimately, the ceRNA regulatory network consisting of 63 lncRNAs, 9 IDDRmiRNAs, and 6 hub OSIDDRGs was constructed (Figure 5(d)).

(a)

(b)

(c)

(d)

3.5. Correlation Analysis between Oxidative Stress and Immunity Abnormality in IDD

There was a significant difference of immunity score between normal IVD tissues and IDD tissues, which indicated that immunity abnormality played important roles in the pathogenesis of IDD (Figure 6(a)). The percentage of IICs for each sample was shown in Figure 6(b), and 8 types of IICs significantly differed between the two groups, namely, IDDRIICs (Figure 6(c)). The correlation analysis among 8 types of IDDRIICs was performed, and 7 significantly correlated pairs (, ) were observed. Neutrophils-B_cells_memory was the most negatively correlated pair (, ), and dendritic_cells_activated-NK_cells_activated was the most positively correlated pair (, ) (Figure 6(d)). Moreover, there was a significant relationship between oxidative stress score and immunity score (, ), indicating there might be a potential link between oxidative stress and immunity abnormity in the development of IDD (Figure 6(e)). To further explore the association of oxidative stress with immunity abnormity in IDD, the correlation analysis between hub OSIDDRGs and IDDRIICs was conducted, 37 significantly correlated pairs (, ) were observed, and the neutrophils-MAPK1 pair was the most positively correlated pair (, ) and B_cells_memory-JAK2 was the most negatively correlated pair (, ) (Figure 6(f)). Combining the 13 miRNA-mRNA pairs with the 37 OSIDDRG-IDDRIIC pairs, a miRNA-mRNA-IDDRIIC regulatory work, which contained 9 IDDRmiRNAs, 6 hub OSIDDRGs, and 8 IDDRIICs, was constructed (Figure 6(g)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

Figure 6

Correlation analysis between oxidative stress and immunity abnormality in IDD. (a) comparison of immunity score between normal IVD tissues and IDD tissues; (b) distribution of IICs in normal IVD tissues and IDD tissues; (c) comparison of IICs between normal IVD tissues and IDD tissues; (d) correlation analysis among 8 types of IDDRIICs; (e) correlation analysis between oxidative stress score and immunity score; (f) correlation analysis between top 10 hub OSIDDRGs and 8 types of IDDRIICs; and (g) the miRNA-mRNA-IDDRIIC regulatory work.

3.6. Correlation Analysis between Oxidative Stress and Autophagy in IDD

An obvious difference of autophagy score was observed between the two groups (9) (Figure 7(a)), which manifested autophagy was involved in the development of IDD. There was a significant relationship between oxidative stress score and autophagy score (, ) (Figure 7(b)), which indicated there might be a potential link between oxidative stress and autophagy in the pathogenesis of IDD. Eighteen ATGIDDRGs were obtained with the intersection between IDD-related genes and autophagy-related genes (Figure 7(c)), and top 10 hub ATGIDDRGs were selected using PPI analysis (Figure 7(d)). Then, the correlation analysis between hub OSIDDRGs and hub ATGIDDRGs was performed and 74 OSIDDRG-ATGIDDRG pairs were observed (Figure 7(e)). The TXNRD1-HSPA5 pair was the most positively correlated pair (, ) (Figure 7(f)), and MAPK1- CDKN1A pair was the most negatively correlated pair (, ) (Figure 7(g)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

4. Discussion

IDD has become a leading cause for LBP, which results in huge economic burden for the patients and medical care system [1–3]. Despite abundant researches in recent years, the definite etiology for IDD remains obscure to date [3, 5, 22]. Recently, accumulating evidence showed that oxidative stress might play important roles in the pathogenesis of IDD, but the underlying mechanism was still unclear [12, 13, 34]. In this study, for the first time, we determined the different status of oxidative stress between normal IVD tissues and IDD tissues using the ssGESA algorithm, which indicated oxidative stress had a very vital role in the IDD. Furthermore, we screened out 72 OSIDDRGs, identified their relevant biological functions, and constructed the ceRNA regulatory network using a series of bioinformatic strategies. Moreover, we also preliminarily elucidated the association of oxidative stress with immunity abnormity and autophagy in IDD. Our findings provided new insights into the researches about the oxidative stress in the development of IDD.

Through several strict bioinformatic algorithms, 10 hub OSIDDRGs were identified, including IL6, PRDX1, MCL1, HMOX1, TXNRD1, MAPK1, HIF1A, FOXO1, JUN, and JAK2. The IL6 was a famous inflammatory cytokine and could promote the initiation and progression of IDD [35]. The HMOX1 could alleviate the senescence of NP cells through inducing the autophagy [36]. The MAPK1 has been reported to be upregulated by EZH2 and promote the IDD via suppressing the miR-129-5p [37]. The HIF1A could alleviate the compression-induced apoptosis of NP-derived stem cells trough upregulating the autophagy via the HIF1A-BNIP3-ATG7 axis [38]. The FOXO1 was a vital regulator of IVD homeostasis through the direct regulation of autophagy, adaptation to hypoxia, and resistance to oxidative stress [39]. The c-Jun was found to improve the expression of TGF-β and promote the NP cell proliferation via reducing the apoptosis and inflammatory response [40]. The JAK2/STAT signaling pathway has been testified to protect the IVD from oxidative stress induced degeneration [41]. To our knowledge, the biological roles of PRDX1, MCL1, and TXNRD1 in IDD have not been investigated up to now. The PRDX1 has been verified to be involved in the regulation of NK cells and protective autophagy in hepatocellular carcinoma [42]. The MCL1 was a famous antiapoptotic protein and could stimulate the progression and drug resistance of thymic carcinoma [43]. The TXNRD1 was a vital ferroptosis-related gene and has been approved to promote the invasion, progression, and metastasis of hepatocellular carcinoma [44]. Based on previous studies [42–44], we made the assumption that PRDX1, MCL1, and TXNRD1 might participate in IDD progression via the regulation of oxidative stress, cell apoptosis, and immunity status, but future studies should be conducted to interpret the precise underlying mechanism. Moreover, to further explore the biological functions of OSIDDRGs, we performed the function enrichment analysis and identified several crucial KEGG pathways which might be involved in IDD. In fact, several of these identified signaling pathways have been verified in IDD, including TNF signaling pathway [45], estrogen signaling pathway [7], MAPK signaling pathway [46], HIF-1 signaling pathway [47], Toll-like receptor signaling pathway [48], and glutathione metabolism [49]. However, how the OSIDDRGs participate in IDD via these signaling pathways has not been fully understood and needs the further investigation.

LncRNA can affect the gene expression by competitively binding miRNAs, namely, ceRNA mechanism, which has been proved to play important and complicated roles in the initiation or progression of human diseases, including IDD [15, 50, 51]. Lan et al. reviewed the latest researches about the interplay between miRNAs and autophagy in the IDD and determined the vital role of miRNA-autophagy pathways in the development of IDD [52]. Similarly, in Jiang et al. study, the potential molecular mechanisms of miRNAs, lncRNAs, and circular RNAs in the progression of IDD were described based on the latest literatures, which suggested that noncoding RNAs could serve as potential targets for the IDD treatment [53]. Moreover, a bioinformatic analysis conducted by Wang et al. showed that miRNAs could function as novel targets for preventing and treating IDD by regulating their target genes [54]. In a similar manner, a meta-analysis conducted by Sherafatian et al. also observed the vital potential roles of abnormal miRNAs in the development of IDD [55]. In the current study, we constructed a strict ceRNA regulatory network of lncRNA-miRNA-mRNA, which was composed of 63 lncRNAs, 9 candidate miRNAs, and 6 hub OSIDDRGs, such as TMEM191C-miRNA-1184-FOXO1 axis and C1orf143-miRNA-4747-5p-TXNRD1. These ceRNA pathways shed new lights on the investigation of oxidative stress in IDD, which are worth our validating and investigating in the future.

Immunity abnormality is another crucial contributor to the development of IDD, and many inflammatory cytokines and IICs have been identified to aggravate the IDD [3, 5, 19]. In our study, we confirmed the differential immunity status between normal IVD tissues and IDD tissues using the ssGSEA algorithm. Subsequently, we identified 8 types of IDDRIICs and performed the correlation analysis among them. We discovered that neutrophils-B_cell_memory pair and dendritic_cells_activated- NK_cells_activated pair were the most negatively and positively correlated pair, respectively. These significantly correlated pairs were worth well of further investigation because few studies have focused on the interaction among IICs in IDD up to now. To further elucidate the potential mechanism of the link between oxidative stress and immunity abnormality in IDD, like the ceRNA regulatory network of lncRNAs, we constructed the miRNA-mRNA-IDDRIICs regulatory network, such as miR-1184- FOXO1-neutrophils axis and miR-4533-MCL1-B_cells_memory axis. We hypothesized these miRNA-mRNA pairs might serve vital roles by changing the percentage of IDDRIICs in the pathogenesis of IDD, which deserved the special attention in future researches.

Previous studies have indicated that oxidative stress might promote the IDD via the regulation of autophagy [22, 34, 56]. Tang et al. used the hydrogen-peroxide to induce the oxidative stress in IDD, and they observed an obviously increased expression of LC3-II protein, indicating that oxidative stress could promote the formation of autophagosome in NP cells [34]. Similarly, Xu et al. study showed that a treatment of 400 μM tert-butyl-hydroperoxide for 6 h could also induce autophagy in NP cells by increasing the expression of LC3-II and Beclin-1 [57]. In our study, there was a significant correlation between oxidative stress score and autophagy score, suggesting a potential interaction between oxidative stress and autophagy in the pathomechanism of IDD. We identified 10 hub ATGIDDRGs through the PPI analysis and then attempted to interpret their relationship with OSIDDRGs using the Pearson analysis. Interestingly, most of hub OSIDDRGs had a significant association with hub ATGIDDRGs, indicating a close relationship between oxidative stress and autophagy in IDD. For instance, MAPK1-CDKNA1 pair was the most negatively correlated pair, which has not been fully investigated to date and deserved our major attention in the future.

Several limitations should be considered when interpreting our results. First, different from tumors, only few sequencing files of IDD were uploaded into public databases, which might be attributed to the rareness of control normal IVD tissues. As a result, the sample size of this study was relatively small, which might affect the stringency of our findings. Second, although the bioinformatic strategies used were strict and widely recognized, all findings in this study were obtained only based on the microarray data from GEO database, which needed further experimental verification. Third, although we tried to elucidate the association of oxidative stress with immunity abnormality or autophagy in IDD, the analysis performed in the current study was relatively thin, and in-depth scientific researches were needed. Forth, the retrospective design of this research would inevitably induce the bias, which might affect the reliability of results; therefore, prospective researches with the strict design were necessary to confirm our findings in the future.

Taken together, we found that normal IVD tissues and IDD tissues could be distinctly distinguished via the expression of OSRGs, and there was a significantly different oxidative stress score between the two groups, which indicated that oxidative stress might play important roles in the development of IDD. Then, through a series of strict bioinformatic strategies, the biological functions of OSIDDRGs were elucidated and 10 hub OSIDDRGs were identified (e.g., IL6, PRDX1, and MCL1). Moreover, the ceRNA network of lncRNA-miRNA-mRNA was constructed to further interpret the potential mechanism of oxidative stress in IDD. Furthermore, the relationship of oxidative stress to immunity abnormality and autophagy in the development of IDD was preliminarily evaluated. Our findings could provide new insights into the oxidative stress in IDD and deliver potential targets for the treatment of IDD.

Abbreviations

LBP:	Low back pain
IVD:	Intervertebral disc
IDD:	Intervertebral disc degeneration
NP:	Nucleus pulposus
LncRNA:	Long noncoding RNA
miRNA:	microRNA
ceRNA:	Competing endogenous RNA
IICs:	Infiltrating immune cells
OSRGs:	Oxidative stress-related genes
GEO:	Gene Expression Omnibus
ImmPort:	Immunology database and Analysis Portal database
HADb:	Human autophagy database
ssGSEA:	Single sample gene set enrichment analysis
OSIDDRGs:	Oxidative stress-related and IDD-related genes
IDDRmiRNAs:	IDD-related miRNAs
GO:	Gene Ontology
KEGG:	Kyoto Encyclopedia of Genes and Genomes analysis
PPI:	Protein-protein interaction
IDDRIICs:	IDD-related IICs.

Data Availability

The original contributions presented in the study are included in the article/supplementary material, and further inquiries can be directed to the corresponding author/s.

Ethical Approval

This study has been approved by the Ethics Committee of Peking University Third Hospital. The written informed consent was not necessary because all data was extracted from published studies.

Conflicts of Interest

The author(s) declare(s) no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Authors’ Contributions

Yongzhao Zhao, Qian Xiang, and Jialiang Lin contributed equally to this research and were listed as the co-first authors.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant No. 82172480).

Supplementary Materials

Supplementary Table 1: detailed information of ten tissues samples. Supplementary Table 2: the list of oxidative stress-related genes. Supplementary Table 3: the list of immunity-related genes. Supplementary Table 4: the list of autophagy-related genes. Supplementary Table 5: the details of 72 OSIDDRGs. Supplementary Table 6: the predicted lncRNAs for miRNAs using LncBase database. Supplementary Figure 1: the intersection between significantly related IDDRmiRNAs and predicted miRNAs for each hub OSIDDRGs. (a) IL6; (b) PRDX1; (c) MCL1; (d) HMOX1; (e) TXNRD1; (f) MAPK1; (g) HIF1A; (h) FOXO1; (i) JUN; (j) JAK2. (Supplementary Materials)

References

R. Buchbinder, M. van Tulder, B. Öberg et al., “Low back pain: a call for action,” Lancet, vol. 391, no. 10137, pp. 2384–2388, 2018.
View at: Publisher Site | Google Scholar
T. Vos, S. S. Lim, C. Abbafati et al., “Global burden of 369 diseases and injuries in 204 countries and territories, 1990-2019: A systematic analysis for the global burden of disease study 2019,” Lancet, vol. 396, no. 10258, pp. 1204–1222, 2020.
View at: Google Scholar
V. Francisco, J. Pino, M. González-Gay et al., “A new immunometabolic perspective of intervertebral disc degeneration,” Nature Reviews Rheumatology, vol. 18, no. 1, pp. 47–60, 2022.
View at: Publisher Site | Google Scholar
Y. Zhang, B. Yang, J. Wang et al., “Cell senescence: a nonnegligible cell state under survival stress in pathology of intervertebral disc degeneration,” Oxidative Medicine and Cellular Longevity, vol. 2020, Article ID 9503562, 12 pages, 2020.
View at: Publisher Site | Google Scholar
M. V. Risbud and I. M. Shapiro, “Role of cytokines in intervertebral disc degeneration: pain and disc content,” Nature Reviews Rheumatology, vol. 10, no. 1, pp. 44–56, 2014.
View at: Publisher Site | Google Scholar
J. C. Iatridis, K. Godburn, K. Wuertz, M. Alini, and P. J. Roughley, “Region-dependent aggrecan degradation patterns in the rat intervertebral disc are affected by mechanical loading in vivo,” Spine, vol. 36, no. 3, pp. 203–209, 2011.
View at: Google Scholar
S. Yang, F. Zhang, J. Ma, and W. Ding, “Intervertebral disc ageing and degeneration: the antiapoptotic effect of oestrogen,” Ageing Research Reviews, vol. 57, article 100978, 2020.
View at: Publisher Site | Google Scholar
Q. Xiang, L. Kang, J. Wang et al., “CircRNA-CIDN mitigated compression loading-induced damage in human nucleus pulposus cells via miR-34a-5p/SIRT1 axis,” eBioMedicine, vol. 53, pp. 102679–102679, 2020.
View at: Publisher Site | Google Scholar
X. Bao, Z. Wang, Q. Jia et al., “HIF-1α-mediated miR-623 regulates apoptosis and inflammatory responses of nucleus pulposus induced by oxidative stress via targeting TXNIP,” Oxidative Medicine and Cellular Longevity, vol. 2021, Article ID 6389568, 17 pages, 2021.
View at: Publisher Site | Google Scholar
Y. Hu, R. Tao, L. Wang et al., “Exosomes Derived from Bone Mesenchymal Stem Cells Alleviate Compression- Induced Nucleus Pulposus Cell Apoptosis by Inhibiting Oxidative Stress,” Oxidative Medicine and Cellular Longevity, vol. 2021, Article ID 2310025, 12 pages, 2021.
View at: Publisher Site | Google Scholar
F. Li, X. Sun, B. Zheng et al., “Arginase II promotes intervertebral disc degeneration through exacerbating senescence and apoptosis caused by oxidative stress and inflammation via the NF-κB pathway,” Frontiers in Cell and Development Biology, vol. 9, article 737809, 2021.
View at: Publisher Site | Google Scholar
A. Dimozi, E. Mavrogonatou, A. Sklirou, and D. Kletsas, “Oxidative stress inhibits the proliferation, induces premature senescence and promotes a catabolic phenotype in human nucleus pulposus intervertebral disc cells,” European Cells & Materials, vol. 30, no. 89-102, pp. 89–103, 2015.
View at: Publisher Site | Google Scholar
D. Seol, M. C. Coleman, J. A. Martin et al., “Targeting oxidative stress with amobarbital to prevent intervertebral disc degeneration: part I. in vitro and ex vivo studies,” The Spine Journal, vol. 21, no. 6, pp. 1021–1030, 2021.
View at: Publisher Site | Google Scholar
G. Cao, S. Yang, J. Cao et al., “The role of oxidative stress in intervertebral disc degeneration,” Oxidative Medicine and Cellular Longevity, vol. 2022, Article ID 2166817, 16 pages, 2022.
View at: Publisher Site | Google Scholar
P. McDonel and M. Guttman, “Approaches for understanding the mechanisms of long noncoding RNA regulation of gene expression,” Cold Spring Harbor Perspectives in Biology, vol. 11, no. 12, 2019.
View at: Publisher Site | Google Scholar
H. Sabaie, N. Amirinejad, M. R. Asadi et al., “Molecular insight into the therapeutic potential of long non-coding RNA-associated competing endogenous RNA axes in Alzheimer's disease: a systematic scoping review,” Frontiers in Aging Neuroscience, vol. 13, article 742242, 2021.
View at: Publisher Site | Google Scholar
H. Sabaie, Z. Salkhordeh, M. R. Asadi et al., “Long non-coding RNA- associated competing endogenous RNA axes in T-cells in multiple sclerosis,” Frontiers in Immunology, vol. 12, article 770679, 2021.
View at: Publisher Site | Google Scholar
H. Tang, S. Yuan, T. Chen, and P. Ji, “Development of an immune-related lncRNA–miRNA–mRNA network based on competing endogenous RNA in periodontitis,” Journal of Clinical Periodontology, vol. 48, no. 11, pp. 1470–1479, 2021.
View at: Publisher Site | Google Scholar
M. Molinos, C. R. Almeida, J. Caldeira, C. Cunha, R. M. Gonçalves, and M. A. Barbosa, “Inflammation in intervertebral disc degeneration and regeneration,” J R Soc Interface, vol. 12, no. 104, p. 20141191, 2015.
View at: Publisher Site | Google Scholar
A. M. Cuervo, “Autophagy and aging: keeping that old broom working,” Trends in Genetics, vol. 24, no. 12, pp. 604–612, 2008.
View at: Publisher Site | Google Scholar
L. P. Nan, F. Wang, Y. Liu et al., “6-gingerol protects nucleus pulposus-derived mesenchymal stem cells from oxidative injury by activating autophagy,” World J Stem Cells, vol. 12, no. 12, pp. 1603–1622, 2020.
View at: Publisher Site | Google Scholar
R. Kritschil, M. Scott, G. Sowa, and N. Vo, “Role of autophagy in intervertebral disc degeneration,” Journal of Cellular Physiology, vol. 237, no. 2, pp. 1266–1284, 2022.
View at: Publisher Site | Google Scholar
A. Liberzon, C. Birger, H. Thorvaldsdóttir, M. Ghandi, J. P. Mesirov, and P. Tamayo, “The molecular signatures database hallmark gene set collection,” Cell Systems, vol. 1, no. 6, pp. 417–425, 2015.
View at: Publisher Site | Google Scholar
S. Bhattacharya, P. Dunn, C. G. Thomas et al., “ImmPort, toward repurposing of open access immunological assay data for translational and clinical research,” Sci Data, vol. 5, no. 1, article 180015, 2018.
View at: Publisher Site | Google Scholar
E. Moussay, T. Kaoma, J. Baginska et al., “The acquisition of resistance to TNFα in breast cancer cells is associated with constitutive activation of autophagy as revealed by a transcriptome analysis using a custom microarray,” Autophagy, vol. 7, no. 7, pp. 760–770, 2011.
View at: Publisher Site | Google Scholar
D. A. Barbie, P. Tamayo, J. S. Boehm et al., “Systematic RNA interference reveals that oncogenic _KRAS_ -driven cancers require TBK1,” Nature, vol. 462, no. 7269, pp. 108–112, 2009.
View at: Publisher Site | Google Scholar
C. Zhu, H. Xiao, X. Jiang, R. Tong, and J. Guan, “Prognostic biomarker DDOST and its correlation with immune infiltrates in hepatocellular carcinoma,” Frontiers in Genetics, vol. 12, pp. 819520–819520, 2022.
View at: Publisher Site | Google Scholar
X. Lin, J. Deng, H. Deng et al., “Comprehensive analysis of the immune microenvironment in checkpoint inhibitor pneumonitis,” Frontiers in Immunology, vol. 12, article 818492, 2022.
View at: Publisher Site | Google Scholar
W. Huang da, 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. Szklarczyk, A. L. Gable, K. C. Nastou et al., “The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets,” Nucleic Acids Research, vol. 49, no. D1, pp. D605–d612, 2021.
View at: Publisher Site | Google Scholar
S. E. McGeary, K. S. Lin, C. Y. Shi et al., “The biochemical basis of microRNA targeting efficacy,” Science, vol. 366, no. 6472, 2019.
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
A. M. Newman, C. L. Liu, M. R. Green et al., “Robust enumeration of cell subsets from tissue expression profiles,” Nature Methods, vol. 12, no. 5, pp. 453–457, 2015.
View at: Publisher Site | Google Scholar
Z. Tang, B. Hu, F. Zang, J. Wang, X. Zhang, and H. Chen, “Nrf2 drives oxidative stress-induced autophagy in nucleus pulposus cells via a Keap1/Nrf2/p62 feedback loop to protect intervertebral disc from degeneration,” Cell Death & Disease, vol. 10, no. 7, p. 510, 2019.
View at: Publisher Site | Google Scholar
J. Chen, Z. Mei, B. Huang et al., “IL-6/YAP1/β-catenin signaling is involved in intervertebral disc degeneration,” Journal of Cellular Physiology, vol. 234, no. 5, pp. 5964–5971, 2019.
View at: Publisher Site | Google Scholar
W. Yi, H. Lan, Y. Wen et al., “HO-1 overexpression alleviates senescence by inducing autophagy via the mitochondrial route in human nucleus pulposus cells,” Journal of Cellular Physiology, vol. 235, no. 11, pp. 8402–8415, 2020.
View at: Publisher Site | Google Scholar
M. Zhou, S. J. He, W. Liu et al., “EZH2 upregulates the expression of MAPK1 to promote intervertebral disc degeneration via suppression of miR-129-5p,” The Journal of Gene Medicine, vol. 24, no. 3, article e3395, 2022.
View at: Publisher Site | Google Scholar
R. He, Z. Wang, M. Cui et al., “HIF1A alleviates compression-induced apoptosis of nucleus pulposus derived stem cells via upregulating autophagy,” Autophagy, vol. 17, no. 11, pp. 3338–3360, 2021.
View at: Publisher Site | Google Scholar
O. Alvarez-Garcia, T. Matsuzaki, M. Olmer et al., “FOXO are required for intervertebral disk homeostasis during aging and their deficiency promotes disk degeneration,” Aging Cell, vol. 17, no. 5, article e12800, 2018.
View at: Publisher Site | Google Scholar
M. Lei, K. Wang, S. Li et al., “The c-Jun signaling pathway has a protective effect on nucleus pulposus cells in patients with intervertebral disc degeneration,” Experimental and Therapeutic Medicine, vol. 20, no. 5, p. 123, 2020.
View at: Publisher Site | Google Scholar
S. Dai, T. Liang, X. Shi, Z. Luo, and H. Yang, “Salvianolic acid B protects intervertebral discs from oxidative stress-induced degeneration via activation of the JAK2/STAT3 signaling pathway,” Oxidative Medicine and Cellular Longevity, vol. 2021, Article ID 6672978, 13 pages, 2021.
View at: Publisher Site | Google Scholar
Y. An, J. Jiang, L. Zhou et al., “Peroxiredoxin 1 is essential for natamycin-triggered apoptosis and protective autophagy in hepatocellular carcinoma,” Cancer Letters, vol. 521, pp. 210–223, 2021.
View at: Publisher Site | Google Scholar
D. Müller, P. Mazzeo, R. Koch et al., “Functional apoptosis profiling identifies MCL-1 and BCL-xL as prognostic markers and therapeutic targets in advanced thymomas and thymic carcinomas,” BMC Medicine, vol. 19, no. 1, p. 300, 2021.
View at: Publisher Site | Google Scholar
M. R. McLoughlin, D. J. Orlicky, J. R. Prigge et al., “TrxR1, Gsr, and oxidative stress determine hepatocellular carcinoma malignancy,” Proceedings of the National Academy of Sciences of the United States of America, vol. 116, no. 23, pp. 11408–11417, 2019.
View at: Publisher Site | Google Scholar
K. Liu, J. Wei, G. Li et al., “Fexofenadine protects against intervertebral disc degeneration through TNF signaling,” Frontiers in Cell and Development Biology, vol. 9, article 687024, 2021.
View at: Publisher Site | Google Scholar
H. J. Zhang, H. Y. Liao, D. Y. Bai, Z. Q. Wang, and X. W. Xie, “MAPK /ERK signaling pathway: a potential target for the treatment of intervertebral disc degeneration,” Biomedicine & Pharmacotherapy, vol. 143, article 112170, 2021.
View at: Publisher Site | Google Scholar
Z. Xiong, J. Ding, J. Zhou, S. Yao, J. Zheng, and X. Guo, “Correlation between the HIF-1α/Notch signaling pathway and Modic changes in nucleus pulposus cells isolated from patients with low back pain,” BMC Musculoskeletal Disorders, vol. 21, no. 1, p. 500, 2020.
View at: Publisher Site | Google Scholar
Z. Li, X. Wang, H. Pan et al., “Resistin promotes CCL4 expression through toll-like receptor-4 and activation of the p38-MAPK and NF-κB signaling pathways: implications for intervertebral disc degeneration,” Osteoarthritis and Cartilage, vol. 25, no. 2, pp. 341–350, 2017.
View at: Publisher Site | Google Scholar
D. Yang, D. Wang, A. Shimer, F. H. Shen, X. Li, and X. Yang, “Glutathione protects human nucleus pulposus cells from cell apoptosis and inhibition of matrix synthesis,” Connective Tissue Research, vol. 55, no. 2, pp. 132–139, 2014.
View at: Publisher Site | Google Scholar
D. W. Thomson and M. E. Dinger, “Endogenous microRNA sponges: evidence and controversy,” Nature Reviews. Genetics, vol. 17, no. 5, pp. 272–283, 2016.
View at: Publisher Site | Google Scholar
X. Fan, G. Chen, F. Ma et al., “An lncRNA-miRNA-mRNA-ceRNA network regulates intervertebral disc degeneration: a bioinformatics study based on the dataset analysis,” General Physiology and Biophysics, vol. 40, no. 4, pp. 317–327, 2021.
View at: Publisher Site | Google Scholar
T. Lan, H. Shiyu, Z. Shen, B. Yan, and J. Chen, “New insights into the interplay between miRNAs and autophagy in the aging of intervertebral discs,” Ageing Research Reviews, vol. 65, article 101227, 2021.
View at: Publisher Site | Google Scholar
J. Jiang, Y. Sun, G. Xu, H. Wang, and L. Wang, “The role of miRNA, lncRNA and circRNA in the development of intervertebral disk degeneration (review),” Experimental and Therapeutic Medicine, vol. 21, no. 6, p. 555, 2021.
View at: Publisher Site | Google Scholar
X. Q. Wang, W. Z. Tu, J. B. Guo et al., “A bioinformatic analysis of microRNAs' role in human intervertebral disc degeneration,” Pain Medicine, vol. 20, no. 12, pp. 2459–2471, 2019.
View at: Publisher Site | Google Scholar
M. Sherafatian, H. R. Abdollahpour, F. Ghaffarpasand, S. Yaghmaei, M. Azadegan, and M. Heidari, “MicroRNA expression profiles, target genes, and pathways in intervertebral disk degeneration: a meta-analysis of 3 microarray studies,” World Neurosurgery, vol. 126, pp. 389–397, 2019.
View at: Publisher Site | Google Scholar
G. Zheng, Z. Pan, Y. Zhan et al., “TFEB protects nucleus pulposus cells against apoptosis and senescence via restoring autophagic flux,” Osteoarthritis and Cartilage, vol. 27, no. 2, pp. 347–357, 2019.
View at: Publisher Site | Google Scholar
W.-N. Xu, H.-L. Zheng, R.-Z. Yang et al., “Mitochondrial NDUFA4L2 attenuates the apoptosis of nucleus pulposus cells induced by oxidative stress via the inhibition of mitophagy,” Experimental & Molecular Medicine, vol. 51, no. 11, pp. 1–16, 2019.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Yongzhao Zhao 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

512

Downloads

575

Citations