A correction has been applied to this article in:
Corrigendum: A microarray data analysis investigating the pathogenesis and potential biomarkers of autophagy and ferroptosis in intervertebral disc degeneration
Wenhao Kuang
Cong Jiang
Cheng Yu
Jinwei Hu
Yang Duan
Zhong Chen- Department of Spinal Surgery, Zhujiang Hospital, Southern Medical University, Guangzhou, China
Background: Intervertebral disc degeneration (IDD) entails complex pathological changes and causes lower back pain (LBP). However, there is still a lack of understanding of the mechanisms involved in IDD, particularly regarding the roles of autophagy and ferroptosis. The current study used microarray data to investigate the pathogenesis of IDD and potential biomarkers related to autophagy and ferroptosis in IDD.
Methods: Differentially expressed genes (DEGs) were identified by analyzing the mRNA and miRNA expression profiles of IDD patients from the Gene Expression Omnibus (GEO). The protein-protein interaction network, Gene Ontology (GO) enrichment, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, and gene set enrichment analysis (GSEA) were utilized. The Human Autophagy Database (HADb) and Ferroptosis Database were used in conjunction with hub genes to identify autophagy- and ferroptosis-related genes. The Transcription Factor -hub gene-miRNA network was constructed. Lastly, the expression of DEGs in normal and degenerated nucleus pulposus cells (NPCs) was investigated via the quantitative reverse transcription polymerase chain reaction (qRT-PCR).
Results: A total of 362 DEGs associated with IDD were identified. GO and KEGG analyses indicated that oxidative stress, extracellular matrix, PI3K-AKT signaling pathway, and ferroptosis were key factors in IDD occurrence. GSEA indicated that IDD was associated with changes in autophagy, iron ion homeostasis, extracellular matrix, and oxidative stress. Eighty-nine hub genes were obtained, including five that were autophagy-related and three that were ferroptosis-related. Of these, TP53 and SESN2 were the intersections of autophagy- and ferroptosis-related genes. In qRT-PCR analysis, CANX, SLC38A1, and TP53 were downregulated in degenerative NPCs, whereas GNAI3, SESN2, and VAMP3 were upregulated.
Conclusion: The current study revealed aspects of autophagy- and ferroptosis-related genes involved in IDD pathogenesis, warranting further investigation.
Introduction
Lower back pain (LBP) affects people of all cultures and ages. It reportedly affects 700 million patients worldwide physically, psychologically, and economically and reduces the productivity of society by causing massive economic losses (Francisco et al., 2022). Intervertebral disc degeneration (IDD) is thought to be the most important component of LBP (Zhao et al., 2021). Intervertebral discs (IVDs) are composed of fibrous tissue and cartilage and have an avascular structure that contains three layers: the annulus fibrosus with its outer and inner parts, the central nucleus pulposus, and the endplate. The normal function of an IVD depends on the cells and stroma (Kos et al., 2019). Multiple studies have indicated that various pathological factors, such as extracellular matrix (ECM) degradation, oxidative stress, inflammatory responses, mitochondrial dysfunction, nutritional deficiency, abnormal mechanical loading, and epigenetic changes, are involved in IDD progression (Zhang G. Z. et al., 2020; Zhang et al., 2021; Liang et al., 2022). However, the regulatory mechanisms involved in IDD are complex; thus, many of them remain to be elucidated.
Autophagy is a process in which cells self-digest their own components. It nourishes cells during fasting, clearing them of excess or damaged organelles, misfolded proteins, and invading microorganisms (Levine and Kroemer 2008). In recent years, numerous studies have indicated that disordered autophagy may lead to IDD (Gong and Zhang 2021; Madhu et al., 2021). Studies investigating autophagy in IDD are still at a very early stage, and the physiological role of autophagy in normal IVDs and specific mechanisms involved in autophagy regulation in IVD cells remain unclear (Kritschil et al., 2022). More studies on relevant functions and mechanisms are needed to clarify the role of autophagy in IDD.
The term ferroptosis, coined in 2012 (Dixon et al., 2012), describes a unique form of cell death. It is driven by iron-dependent phospholipid peroxidation and regulated by multiple cellular metabolic events, such as redox homeostasis, mitochondrial activity, and metabolism of amino acids, lipids, and sugars, as well as many disease-related signaling pathways. Ferroptosis may be involved in many organ damage and degenerative pathologies (Jiang et al., 2021). Previous studies have suggested that ferroptosis might play a role in IDD (Zhang GZ. et al., 2020; Lu et al., 2021).
The current study used microarray data analysis to investigate IDD pathogenesis and potential biomarkers related to autophagy or ferroptosis in IDD.
Materials and methods
Microarray datasets
Four datasets were retrieved from the Gene Expression Omnibus (GEO) database (Edgar et al., 2002; Barrett et al., 2012) (http://www.ncbi.nlm.nih.gov/geo): two mRNA expression datasets [GSE23130 (Gruber et al., 2012) and GSE15227 (Gruber et al., 2009)] and two miRNA expression datasets [GSE63492 (Liu et al., 2015) and GSE116726 (Ji et al., 2018)]. GSE15227 contained 12 healthy samples and three IDD samples, GSE23130 contained 15 healthy samples and eight IDD samples, GSE116726 contained three healthy samples and three IDD samples, and GSE63492 contained five healthy samples and five IDD samples.
Identification of differentially expressed genes and hub genes
The GEO2R network tool (https://www.ncbi.nlm.nih.gov/geo/geo2r/) is based on the R software “LIMMA” and GEO query (Barrett et al., 2012), and it was used in the current study to identify differentially expressed mRNAs and miRNAs in healthy and IDD samples. A log2 fold change of >1 and p of <.05 were used as the criteria to identify differentially expressed genes (DEGs). The intersections of the DEGs in the two mRNA datasets and two miRNA datasets were used to obtain differentially expressed mRNAs and miRNAs. Target genes were predicted from the differentially expressed miRNAs using MIRNET (https://www.mirnet.ca/), and hub genes were identified after crossing the target genes with the differentially expressed mRNAs. VEEN diagrams were then constructed to detect and visualize the intersections of these databases.
Gene set enrichment analysis
“c2.all.v7.4.symbols” and “c5.all.v7.4.symbols” gene sets in the MSigDB database (Liberzon et al., 2015) were used as reference gene sets. Then, gene set enrichment analysis (GSEA) was performed (Subramanian et al., 2005) on GSE23130 and GSE15227.
Protein-protein interaction network construction
Genes with a minimum required interaction score of >.4 were selected to construct a differentially expressed mRNA-associated protein-protein interaction (PPI) network using the STRING database. The DEGs were functionally annotated using the Cytoscape 3.8.0 plug-in “clueGO” (Bindea et al., 2009).
Functional enrichment analyses of hub genes
Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses of the hub genes were performed using the R software “clusterProfiler” (Yu et al., 2012), and the results were visualized with the R software “ggplot2”.
Transcription factor-hub gene-miRNA network construction
Network Analyst (https://www.networkanalyst.ca/) was used to predict transcription factors that could regulate hub genes to construct the transcription factor-hub gene-miRNA network, and this network was visualized using Cytoscape (v.3.8.0).
Autophagy-related and ferroptosis-related genes in IDD
Autophagy-related genes were obtained from the Human Autophagy Database (HADb) (http://www.autophagy.lu/index.html) and then intersected with the hub genes to obtain the autophagy-related genes in IDD. Similarly, driver, suppressor, and marker data were obtained from the Ferroptosis Database (FerrDb) (http://www.zhounan.org/ferrdb), with removed duplication, and then intersected to obtain ferroptosis-related genes for subsequent analysis.
Isolation and culture of nucleus pulposus cells
Normal human nucleus pulposus tissue was collected, washed with phosphate-buffered saline (PBS), and then cut into pieces. After digestion with .2% collagenase type II (Thermo Fisher Scientific, United States) for 2 h, followed by .25% digestion with trypsin for .5 h, the samples were washed again with PBS. After centrifugation, they were transferred to culture flasks with Dulbecco’s Modified Eagle’s Medium)/F12 medium containing 10% fetal bovine serum (Thermo Fisher Science, United States) and 1% penicillin-streptomycin (Thermo Fisher Science, United States) at 37°C in a 5% CO2 incubator for cell isolation and cultivation. The cell culture medium was replaced on the fifth day and every two or three days thereafter. After 80%–90% confluency was reached, nucleus pulposus cells (NPCs) were passaged. Third-passage NPCs were used in subsequent experiments. The experimental protocol was approved by the Ethics Committee of Zhujiang Hospital of Southern Medical University.
Quantitative reverse transcription-PCR
NPCs were divided into normal and degeneration groups. In the normal group, the cell culture medium was replaced every 48–72 h. In the degeneration group, the NPCs were treated with 100 ng/ml tumor necrosis factor (TNF)-α for 24 h. Total RNA was extracted and reverse-transcribed via a kit in accordance with the manufacturer’s instructions (Accurate Biotechnology, China). Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was performed using the CFX Connection Real-Time System (Bio-Rad, United States). The 2−ΔΔCt method was used to evaluate gene expression relative to GAPDH expression. Sequence fragments of RNAs are shown in Table 3.
Statistical analysis
Data were expressed as mean ± the standard deviation. Gene expression levels of clinical samples were compared using the independent sample t-test. Benjamini–Hochberg procedure was used to reduce the false-positive rate of multiple tests. A p of <.05 was considered statistically significant. All statistical analyses were conducted using R programming (https://www.r-project.org/, version 4.0.2).
Results
Identification of DEGs and hub genes
A total of 1,817 upregulated and 1,689 downregulated mRNA and 470 upregulated and 538 downregulated miRNA were identified. Of these, 362 overlapping mRNA were defined as differentially expressed mRNA, and 23 overlapping miRNAs were defined as differentially expressed miRNAs (Figures 1A,B). A total of 2773 target mRNAs were predicted from these 23. Eighty-nine hub genes were identified by intersecting the 2773 target mRNAs with 362 differentially expressed mRNAs (Figure 1C; Table 1).
FIGURE 1
FIGURE 1. Differential expression of mRNAs, miRNAs, and hub genes in IDD. (A) Venn diagram of differentially expressed mRNAs in the two mRNA datasets. (B) Venn diagrams of differentially expressed miRNAs in two miRNA datasets. (C) Venn diagram of hub genes in IDD.
TABLE 1
TABLE 1. 89 hub genes.
PPI network analysis
The PPI network consisted of 361 nodes and 1595 edges (Figure 2A). The average node degree was 8.84, and the local clustering coefficient was .397. The three proteins with the highest node degrees were ACTB (90 nodal degrees), TP53 (69 nodal degrees), and FN1 (60 nodal degrees).
FIGURE 2
FIGURE 2. mRNA-related PPI networks in IDD. (A) PPI network of differentially expressed mRNAs in IDD. (B) clueGO-GO analysis of genes in PPIs that had differentially expressed mRNAs associated with IDD. (C) clueGO-KEGG analysis of genes in PPIs that had differentially expressed mRNAs associated with IDD.
GO and KEGG enrichment analyses of DEGs
GO and KEGG analyses of the DEGs obtained from the PPI network were performed to verify the functions of differentially expressed mRNAs in IDD. GO analysis results indicated changes in responses to oxidative stress, reactive oxygen species (ROS) metabolic processes, transforming growth factor beta (TGF-β) receptor signaling pathway, TGF-β, growth factor binding, collagen-containing ECM, ECM organization, mitochondrial respirasome, apoptotic signaling pathway, intrinsic apoptotic signaling pathway in response to oxidative stress, and lysosomes (Figure 2B). KEGG pathway enrichment analysis indicated that the DEGs were significantly enriched in ferroptosis, oxidative phosphorylation, PI3K-Akt signaling pathway, lysosomes, and ECM-receptor interaction (Figure 2C).
GSEA
Genes in GSE15227 were significantly involved in five biological functions: regulation of autophagy, autolysosomes, iron ion transport, responses to oxidative stress, and collagen-containing ECM (Figure 3A). Genes in GSE23130 were also significantly involved in five biological functions: autophagosomes, processes utilizing autophagic mechanisms, iron ion homeostasis, responses to oxidative stress, and collagen-containing ECM (Figure 3B).
FIGURE 3
FIGURE 3. GSEA of the two datasets. (A) GSEA of genes in the GSE15227 dataset. (B) GSEA of genes in the GSE23130 dataset.
GO and KEGG enrichment analyses of hub genes
In GO analysis, hub genes were mainly enriched in ROS metabolic processes and collagen-containing ECM. In KEGG pathway analysis, hub genes were mainly enriched in the PI3K-Akt signaling pathway, ECM-receptor interaction, and cellular senescence (Figure 4; Table 2).
FIGURE 4
FIGURE 4. GO and KEGG enrichment analyses of hub genes.
TABLE 2
TABLE 2. The GO and KEGG enrichment analyses of hub genes in IDD.
TF-hub Gene-miRNA network
A TF-hub gene-miRNA network containing 89 mRNAs, 22 miRNAs, and 81 TFs was generated. The top targeted miRNA hub gene was TP53, which was modulated by seven miRNAs. The top targeted TF hub gene was GPBP1, which was modulated by 20 TFs. The TF that regulated the most hub genes was FOXC1, which regulated 56 hub genes. The DEM that regulated the most hub genes, 22, was hsa-mir-185-5p (Figure 5A). Six genes were obtained after the 81 predicted TFs were intersected with the 89 hub genes: TP53, DAZAP2, HLA-B, MSN, UBLCP1, and PPT1. This indicated that they are simultaneously regulated by TFs and differentially expressed miRNAs and also regulate some hub genes (Figure 5B). TP53 was evidently a key juncture of the network (Figure 5C).
FIGURE 5
FIGURE 5. TF-hub gene-miRNA network. (A) Red nodes represent TFs, green nodes represent miRNAs, and yellow nodes represent hub genes. (B) Venn diagram of hub genes acting as TFs in IDD. (C) Yellow nodes represent hub genes, and yellow nodes with red borders represent hub genes acting as transcription factors.
Autophagy-related and ferroptosis-related hub genes in IDD
Two hundred and twenty-two autophagy-related genes and 487 ferroptosis-related genes were obtained. Intersecting 222 autophagy-related genes with 89 hub genes resulted in the identification of five autophagy-related hub genes in IDD: CANX, GNAI3, SESN2, TP53, and VAMP3. Intersecting 487 ferroptosis-related genes with 89 hub genes resulted in the identification of three ferroptosis-related hub genes in IDD: TP53, SLC38A1, and SESN2. TP53 and SESN2 were simultaneously associated with both autophagy and ferroptosis (Figure 6A). The related TF-hub gene-miRNA network was then constructed (Figure 6B).
FIGURE 6
FIGURE 6. Autophagy-related and ferroptosis-related hub genes in IDD and associated TF-hub gene-miRNA network. (A) Autophagy- and ferroptosis-related hub genes in IDD, indicating that there are two genes simultaneously associated with both autophagy and ferroptosis. (B) Red rectangles represent TFs, green rectangles represent miRNAs, yellow and blue rectangles represent autophagy and ferroptosis genes, respectively, and yellow rectangles with a blue frame represent genes associated with both autophagy and ferroptosis.
Expression levels of hub genes
CANX, SLC38A1, and TP53 were downregulated in degenerative NPCs compared to normal NPCs, whereas GNAI3, SESN2, and VAMP3 were upregulated (Figure 7) (Table 3).
FIGURE 7
FIGURE 7. Expression levels of hub genes in normal NPCs and degenerated NPCs. *p < 0.05; **p < 0.01; ***p < 0.001; ns, not significant.
TABLE 3
TABLE 3. Primer sequences.
Discussion
IDD is an important cause of LBP and involves complex pathophysiological changes and molecular mechanisms. Research on the occurrence and developmental mechanisms involved in IDD has been limited, particularly regarding mechanisms relating to autophagy and ferroptosis. In the current study, microarray and bioinformatics analyses were used to investigate possible molecular mechanisms involved in IDD occurrence and development, and potential biomarkers were predicted.
Both GO and KEGG enrichment analyses of DEGs and hub genes highlighted the importance of oxidative stress, ECM, and PI3K-AKT signaling pathway, and all three of these were related to autophagy and ferroptosis. Activation of the PI3K-AKT signaling pathway can increase ECM content, prevent apoptosis, promote cell proliferation, induce or prevent autophagy, mitigate oxidative damage, and contribute to adaptation to a hypoxic microenvironment, thereby preventing IDD (Ouyang et al., 2017). Exosomes of the cartilage endplate activate the PI3K/AKT/autophagy pathway, inhibit NPC apoptosis, and slow the IDD process (Luo et al., 2021). Islet amyloid polypeptide regulates gene expression responsible for ECM anabolism and catabolism and controls the interaction between apoptosis and NPC autophagy via the PI3K/AKT/mTOR and p38/JNK MAPK signaling pathways (Wu et al., 2017). XIST may act as a competitive endogenous RNA of MIRI-19, and XIST/MIRI-19/PTEN crosstalk reduces NPC viability, promotes NPC degeneration, and induces NPC autophagy in IVD tissues via the PI3K/Akt signaling pathway inactivation, slowing IDD development (Chen W et al., 2021). Bone marrow mesenchymal stem cell-derived extracellular vesicles harboring circ_0072464 inhibit NPC ferroptosis by inhibiting miR-431 and upregulating NRF2, which promotes NPC matrix synthesis and proliferation, alleviating IDD (Yu et al., 2022). Under oxidative stress, ferroprotein imbalance leads to NP intercellular iron overload, resulting in NPC ferroptosis and accelerating IDD progression in vivo (Lu et al., 2021). Homocysteine promotes IDD development by upregulating oxidative stress and ferroptosis in the nucleus pulposus by enhancing GPX4 methylation (Zhang G. Z. et al., 2020). Previous studies have indicated that oxidative stress, ECM, and PI3K-Akt signaling pathway might be clinical biomarkers or therapeutic indicators in IDD, which is consistent with the current study. Based on the present analysis and previous studies, we suggest that the roles of autophagy and ferroptosis in IDD warrant attention.
TP53 encodes a tumor suppressor protein containing transcriptional activation, DNA binding, and oligomerization domains. The p53-mediated transcriptional suppression of SLC7A11 promotes ferroptosis in cancer cells. p53 R273H and R175H inhibit the expression of SLC7A11 by not binding to DNA and by inhibiting the activity of other transcription factors, suggesting that the expression of this core ferroptosis regulator is regulated by an integrated network of transcription factors (Chen X et al., 2021). p53 can enhance ferroptosis by enhancing the expression of SAT1 (spermidine/spermine N1-acetyltransferase 1) and GLS2 (glutaminase 2) or by inhibiting that of SLC7A11 (solute carrier family seven member 11). On the other hand, p53 suppresses ferroptosis through the direct inhibition of DPP4 (dipeptidyl peptidase 4) activity or by the promotion of CDKN1A/p21 (cyclin-dependent kinase inhibitor 1A) expression (Kang et al., 2019). The loss of TP53 prevents nuclear accumulation of DPP4, which promotes plasma membrane-associated DPP4-dependent lipid peroxidation, leading to ferroptosis (Xie et al., 2017). Previous studies have demonstrated the central role of TP53 in the ferroptosis regulatory network, controlling ferroptosis by directly regulating the expression of downstream ferroptosis-related target genes or regulating the expression of other transcription factors to control the expression of death-related genes. Our results indicated that TP53 might regulate downstream ferroptosis-related target gene SESN2 to regulate ferroptosis, while TP53 might be regulated by other transcription factors (such as DAZAP2, HLA-B, MSN, UBLCP1, and PPT1) to form a complex and integrated regulatory network of transcription factors, which was consistent with previous studies.
SESN2 encodes a member of the sestrin family of PA26-related proteins. SESN2 can play an antioxidant role in sepsis, downregulating the ATF4-CHOP-CHAC1 signaling pathway and inhibiting ferroptosis in dendritic cells (Li et al., 2021). Ferroptosis-mediated SESN2 induction protects against iron overload and ferroptosis-induced liver injury, relying on NRF2 (Park et al., 2019). Based on existing studies, we believe that TP53 may be a central gene in the regulatory network of ferroptosis in IDD, and it may play an important role in IDD by regulating SESN2 or a complex transcription factor regulatory network to control ferroptosis.
SLC38A1 functions as a sodium-dependent amino acid transporter. As a competitive endogenous RNA, long non-coding RNA ZFAS1 promotes the conversion of lung fibroblasts to myofibroblasts and ferroptosis via the MIR-150-5P/SLC38A1 axis (Yang et al., 2020). OGFRP1 regulates lung cancer cell proliferation and ferroptosis by inhibiting miR-299-3p to enhance SLC38A1 expression (Liu et al., 2022). SLC7A11 has been confirmed to be regulated by TP53. SLC38A1 and SLC7A11 belong to different family members of solute carrier (SLC) and participate in the regulation of ferroptosis. Whether TP53 and SLC38A1 have a regulatory relationship deserves further attention and research.
TP53 is reported to be involved in autophagy. DNA damage-regulated autophagy modulator and TP53-targeted autophagy-related genes, including unc-51-like autophagy activating kinase one and autophagy-related 7, are also essential for TP53-mediated apoptosis. Cytoplasmic TP53 inhibits AMP-dependent kinase, which, in turn, activates MTOR, leading to the inhibition of autophagy, whereas nuclear TP53 has the opposite effect. Mild stressors activate TP53, which, in turn, upregulates the level of basal autophagy without activation of MTOR-dependent induced autophagy. TP53/p53-FBXO22-TFEB controls basal autophagy, regulating hormesis (Suzuki et al., 2021). Increased levels of p53 acetylation suppress renal tubular epithelial cells (RTEC) autophagy after sepsis. Autophagy activation induced by p53 following deacetylation by Sirt1 reduces sepsis-associated acute kidney injury (SAKI) (Sun et al., 2021). UBE2T upregulates autophagy in NSCLC cells by activating the p53/AMPK/mTOR signaling pathway (Zhu et al., 2021).
CANX encodes a member of the calnexin family of molecular chaperones. Blocked mitosis in hypoxic cells may be due to the fact that after CANX knockout, the number of elongated mitochondria is increased, and the colocalization of autophagosomes and mitochondria is reduced (Wu et al., 2016). The endoplasmic reticulum (ER)-resident lectin chaperone Calnexin (CANX) and ER-phagy receptor FAM134B are required for autophagy-mediated quality control of endogenous procollagens (Forrester et al., 2019). VAMP3 is a member of the vesicle-associated membrane protein (VAMP)/synaptobrevin family. VAMP3 is required for the fusion between autophagosomes and multivesicular bodies (MVBs) in K562 cells (Fader et al., 2009). The STX6-VTI1B-VAMP3 complex facilitates xenophagy by regulating the fusion between recycling endosomes and autophagosomes (Nozawa et al., 2017). TP53 is not only an autophagy-related gene but also an important and potent ferroptosis-related gene. Revealing its role in IDD autophagy and regulatory axis of TP53-CANX/TP53-VAMP3 may inspire new ideas for the research on the regulatory network of IDD autophagy.
GNAI3 belongs to the G-alpha family. GNAI3 reduces autophagy flux and inhibits the autophagic degradation of hepatitis B virus antigen by enhancing the insulin-induced AKT-MTOR pathway (Wang et al., 2022). The activator of G-protein signaling 3 (AGS3), which is directly bound to light chain 3 (LC3), recruits Gαi3 to the LC3-positive membrane upon starvation and enhances autophagy by inhibiting G protein. Gα-interacting vesicle-associated protein (GIV) destroys the Gαi3-AGS three complex under the stimulation of growth factors, releases Gαi3 from the LC3-positive membrane, and inhibits autophagy by activating G protein (Garcia-Marcos et al., 2011). At present, there is no research involving GNAI3 and its encoded Gαi3 in the field of IDD. Elucidation of the role of GNAI3 in IDD may contribute to a further understanding of autophagy in IDD.
The present study had some limitations. First, the data used were downloaded from public databases, and it was not possible to evaluate whether some of those data contained errors. Second, only IDD tissue samples were analyzed, but a comprehensive analysis of venous blood samples and IDD tissue may be more helpful with respect to providing a general understanding of the potential mechanisms involved in IDD, which would be an improvement for future research. Lastly, a large amount of evidence is still needed to support our hypothesis. In subsequent studies, we will use various experiments, such as western blotting, immunohistochemistry assays, and dual luciferase reporter gene assays, to verify our conclusions and further clarify the mechanisms of IDD autophagy and ferroptosis.
In conclusion, the purpose of this study was to investigate the molecular mechanisms involved in IDD progression via comprehensive bioinformatics analysis to complement the understanding of autophagy and ferroptosis in IDD. This study might lead to the discovery of new IDD biomarkers, and more research is required to validate our results.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving human participants were reviewed and approved by the Ethics Committee of Zhujiang Hospital of Southern Medical University. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.
Author contributions
WK and CJ: data analysis, validation, and manuscript writing. CY, JH, and YD: the extraction and culture of nucleus pulposus cells and qRT-PCR experiments. ZC: the design of the subject. All authors read and approved the final manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.1090467/full#supplementary-material
References
Barrett, T., Wilhite, S. E., Ledoux, P., Evangelista, C., Kim, I. F., Tomashevsky, M., et al. (2012). NCBI GEO: Archive for functional genomics data sets—update. Nucleic Acids Res. 41 (D1), D991–D995. doi:10.1093/nar/gks1193
Bindea, G., Mlecnik, B., Hubert, H., Charoentong, P., Tosolini, M., Amos, K., et al. (2009). ClueGO: A Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics 25 (8), 1091–1093. doi:10.1093/bioinformatics/btp101
Chen W, W., Li, S., and Zhang, F. (2021). Role of lncRNA XIST/microRNA-19/PTEN network in autophagy of nucleus pulposus cells in intervertebral disc degeneration via the PI3K/Akt signaling pathway. Cell cycleGeorget. Tex.) 20 (17), 1629–1641. doi:10.1080/15384101.2021.1924450
Chen X, X., Kang, R., Kroemer, G., and Tang, D. (2021). Broadening horizons: The role of ferroptosis in cancer. Nat. Rev. Clin. Oncol. 18 (5), 280–296. doi:10.1038/s41571-020-00462-0
Dixon, S. J., Lemberg, K. M., Lamprecht, M. R., Skouta, R., Zaitsev, E. M., Gleason, C. E., et al. (2012). Ferroptosis: An iron-dependent form of nonapoptotic cell death. Cell 149 (5), 1060–1072. doi:10.1016/j.cell.2012.03.042
Edgar, R., Domrachev, M., and Lash, A. E. (2002). Gene expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 30 (1), 207–210. doi:10.1093/nar/30.1.207
Fader, C. M., Germán Sánchez, D., Mestre, M. B., and Colombo, M. I. (2009). TI-VAMP/VAMP7 and VAMP3/cellubrevin: Two v-SNARE proteins involved in specific steps of the autophagy/multivesicular body pathways. Biochimica Biophysica Acta - Mol. Cell Res. 1793 (12), 1901–1916. doi:10.1016/j.bbamcr.2009.09.011
Forrester, A., De Leonibus, C., Grumati, P., Fasana, E., Piemontese, M., Staiano, L., et al. (2019). A selective ER-phagy exerts procollagen quality control via a Calnexin-FAM134B complex. EMBO J. 38 (2), e99847. doi:10.15252/embj.201899847
Francisco, V., Pino, J., Ángel González-Gay, M., Lago, F., Karppinen, J., Tervonen, O., et al. (2022). A new immunometabolic perspective of intervertebral disc degeneration. Nat. Rev. Rheumatol. 18 (1), 47–60. doi:10.1038/s41584-021-00713-z
Garcia-Marcos, M., Ear, J., Farquhar, M. G., Ghosh, P., and Yoshimori, T. (2011). A GDI (AGS3) and a GEF (GIV) regulate autophagy by balancing G protein activity and growth factor signals. Mol. Biol. Cell 22 (5), 673–686. doi:10.1091/mbc.e10-08-0738
Gong, C-Y., and Zhang, H-H. (2021). Autophagy as a potential therapeutic target in intervertebral disc degeneration. Life Sci. 273, 119266. doi:10.1016/j.lfs.2021.119266
Gruber, H. E., Hoelscher, G., Loeffler, B., Yin, C., Ingram, J. A., Halligan, W., et al. (2009). Prostaglandin E1 and misoprostol increase epidermal growth factor production in 3D-cultured human annulus cells. Spine J. 9 (9), 760–766. doi:10.1016/j.spinee.2009.04.024
Gruber, H. E., Hoelscher, G. L., Ingram, J. A., and Hanley, E. N. (2012). Genome-wide analysis of pain-nerve- and neurotrophin -related gene expression in the degenerating human annulus. Mol. Pain 8, 63. doi:10.1186/1744-8069-8-63
Ji, M-L., Jiang, H., Zhang, X-J., Shi, P-L., Li, C., Wu, H., et al. (2018). Preclinical development of a microRNA-based therapy for intervertebral disc degeneration. Nat. Commun. 9, 5051. doi:10.1038/s41467-018-07360-1
Jiang, X., Stockwell, B. R., and Conrad, M. (2021). Ferroptosis: Mechanisms, biology and role in disease. Nat. Rev. Mol. Cell Biol. 22 (4), 266–282. doi:10.1038/s41580-020-00324-8
Kang, R., Guido, K., and Tang, D. (2019). The tumor suppressor protein p53 and the ferroptosis network. Free Radic. Biol. Med. 133, 162–168. doi:10.1016/j.freeradbiomed.2018.05.074
Kos, N., Gradisnik, L., and Velnar, T. (2019). A brief review of the degenerative intervertebral disc disease. Med. Arch. 73 (6), 421–424. doi:10.5455/medarh.2019.73.421-424
Kritschil, R., Scott, M., Sowa, G., and Vo, N. (2022). Role of autophagy in intervertebral disc degeneration. J. Cell. Physiology 237 (2), 1266–1284. doi:10.1002/jcp.30631
Levine, B., and Kroemer, G. (2008). Autophagy in the pathogenesis of disease. Cell 132 (1), 27–42. doi:10.1016/j.cell.2007.12.018
Li, J. Y., Ren, C., Wang, L. X., Yao, R. Q., Dong, N., Wu, Y., et al. (2021). Sestrin2 protects dendrite cells against ferroptosis induced by sepsis. Cell Death Dis. 12 (9), 834. doi:10.1038/s41419-021-04122-8
Liang, Huaizhen, Luo, Rongjin, Li, Gaocai, Zhang, Weifeng, Song, Yu, and Yang, Cao (2022). The proteolysis of ECM in intervertebral disc degeneration. Int. J. Mol. Sci. 23 (3), 1715. doi:10.3390/ijms23031715
Liberzon, A., Birger, C., Thorvaldsdóttir, H., Ghandi, M., Mesirov, J. P., and Tamayo, P. (2015). The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst. 1 (6), 417–425. doi:10.1016/j.cels.2015.12.004
Liu, X., Lu, C., Xie, Y-K., Hu, Q-J., Ma, C-J., Pei, Y-J., et al. (2015). Noncoding RNAs in human intervertebral disc degeneration: An integrated microarray study. Genomics Data 5, 80–81. doi:10.1016/j.gdata.2015.05.027
Liu, L., Su, S., Ye, D., Yu, Z., Lu, W., and Li, X. (2022). Long non-coding RNA OGFRP1 regulates cell proliferation and ferroptosis by miR-299-3p/SLC38A1 axis in lung cancer. Anti-Cancer Drugs 33 (9), 826–839. doi:10.1097/CAD.0000000000001328
Lu, S., Song, Y., Luo, R., Li, S., Li, G., Wang, K., et al. (2021). Ferroportin-dependent iron homeostasis protects against oxidative stress-induced nucleus pulposus cell ferroptosis and ameliorates intervertebral disc segeneration in vivo. Oxid. Med. Cell. Longev. 2021, 1–18. doi:10.1155/2021/6670497
Luo, L., Jian, X., Sun, H., Qin, J., Wang, Y., Zhang, J., et al. (2021). Cartilage endplate stem cells inhibit intervertebral disc degeneration by releasing exosomes to nucleus pulposus cells to activate Akt/autophagy. Stem Cells 39 (4), 467–481. doi:10.1002/stem.3322
Madhu, V., Guntur, A. R., and Risbud, M. V. (2021). Role of autophagy in intervertebral disc and cartilage function: Implications in health and disease. Matrix Biol. 100-101, 207–220. doi:10.1016/j.matbio.2020.12.002
Nozawa, T., Minowa-Nozawa, A., Aikawa, C., and Nakagawa, I. (2017). The STX6-VTI1B-VAMP3 complex facilitates xenophagy by regulating the fusion between recycling endosomes and autophagosomes. Autophagy 13 (1), 57–69. doi:10.1080/15548627.2016.1241924
Ouyang, Z. H., Wang, W. J., Yan, Y. G., Wang, B., and Lv, G. H. (2017). The PI3K/akt pathway: A critical player in intervertebral disc degeneration. Oncotarget 8 (34), 57870–57881. doi:10.18632/oncotarget.18628
Park, S. J., Kim, J. H., Yang, J. W., Han, Y., Ku, S. K., Cho, I. J., et al. (2019). Protective effect of sestrin2 against iron overload and ferroptosis-induced liver injury. Toxicol. Appl. Pharmacol. 379, 114665. doi:10.1016/j.taap.2019.114665
Subramanian, A., Tamayo, P., Mootha, V. K., Mukherjee, S., Ebert, B. L., Gillette, M. A., et al. (2005). Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. 102 (43), 15545–15550. doi:10.1073/pnas.0506580102
Sun, M., Li, J., Mao, L., Wu, J., Deng, Z., He, M., et al. (2021). p53 deacetylation alleviates sepsis-induced acute kidney injury by promoting autophagy. Front. Immunol. 12, 685523. doi:10.3389/fimmu.2021.685523
Suzuki, N., Johmura, Y., Wang, T. W., Migita, T., Wu, W., Noguchi, R., et al. (2021). TP53/p53-FBXO22-TFEB controls basal autophagy to govern hormesis. Autophagy 17 (11), 3776–3793. doi:10.1080/15548627.2021.1897961
Wang, X., Wei, Z., Lan, T., He, Y., Cheng, B., Li, R., et al. (2022). CCDC88A/GIV promotes HBV replication and progeny secretion via enhancing endosomal trafficking and blocking autophagic degradation. Autophagy 18 (2), 357–374. doi:10.1080/15548627.2021.1934271
Wu, W., Li, W., Chen, H., Jiang, L., Zhu, R., and Feng, D. (2016). FUNDC1 is a novel mitochondrial-associated-membrane (MAM) protein required for hypoxia-induced mitochondrial fission and mitophagy. Autophagy 12 (9), 1675–1676. doi:10.1080/15548627.2016.1193656
Wu, X., Song, Y., Liu, W., Wang, K., Gao, Y., Li, S., et al. (2017). IAPP modulates cellular autophagy, apoptosis, and extracellular matrix metabolism in human intervertebral disc cells. Cell Death Discov. 3, 16107. doi:10.1038/cddiscovery.2016.107
Xie, Y., Zhu, S., Song, X., Sun, X., Fan, Y., Liu, J., et al. (2017). The tumor suppressor p53 limits ferroptosis by blocking DPP4 activity. Cell Rep. 20 (7), 1692–1704. doi:10.1016/j.celrep.2017.07.055
Yang, Y., Tai, W., Lu, N., Li, T., Liu, Y., Wu, W., et al. (2020). lncRNA ZFAS1 promotes lung fibroblast-to-myofibroblast transition and ferroptosis via functioning as a ceRNA through miR-150-5p/SLC38A1 axis. Aging (Albany NY) 12 (10), 9085–9102. doi:10.18632/aging.103176
Yu, G., Wang, L-G., Han, Y., and He, Q-Y. (2012). clusterProfiler: an R Package for comparing biological themes among gene clusters. OMICS A J. Integr. Biol. 16 (5), 284–287. doi:10.1089/omi.2011.0118
Yu, X., Xu, H., Liu, Q., Wang, Y., Wang, S., Lu, R., et al. (2022). circ_0072464 shuttled by bone mesenchymal stem cell-secreted extracellular vesicles inhibits nucleus pulposus cell ferroptosis to relieve intervertebral disc degeneration. Oxidative Med. Cell. Longev. 2022, 1–21. doi:10.1155/2022/2948090
Zhang, G. Z., Liu, M. Q., Chen, H. W., Wu, Z. L., Gao, Y. C., Ma, J., et al. (2021). NF-κB signalling pathways in nucleus pulposus cell function and intervertebral disc degeneration. Cell Prolif. 54, e13057. doi:10.1111/cpr.13057
Zhang, G. Z., Deng, Y. J., Xie, Q. Q., Ren, E. H., Ma, Z. J., He, X. G., et al. (2020a). Sirtuins and intervertebral disc degeneration: Roles in inflammation, oxidative stress, and mitochondrial function. Clin. Chim. Acta 508, 33–42. doi:10.1016/j.cca.2020.04.016
Zhang, G. Z., Liu, M. Q., Chen, H. W., Wu, Z. L., Gao, Y. C., Ma, Z. J., et al. (2020b). Homocysteine induces oxidative stress and ferroptosis of nucleus pulposus via enhancing methylation of GPX4. Free Radic. Biol. Med. 160, 552–565. doi:10.1016/j.freeradbiomed.2020.08.029
Zhao, K., Ran, A., Xiang, Q., Li, G., Wang, K., Song, Y., et al. (2021). Acid-sensing ion channels regulate nucleus pulposus cell inflammation and pyroptosis via the NLRP3 inflammasome in intervertebral disc degeneration. Cell Prolif. 54, e12941. doi:10.1111/cpr.12941
Keywords: autophagy, DEGs, ferroptosis, hub genes, integrated bioinformatics, intervertebral disc degeneration
Citation: Kuang W, Jiang C, Yu C, Hu J, Duan Y and Chen Z (2023) A microarray data analysis investigating the pathogenesis and potential biomarkers of autophagy and ferroptosis in intervertebral disc degeneration. Front. Genet. 13:1090467. doi: 10.3389/fgene.2022.1090467
Received: 05 November 2022; Accepted: 13 December 2022;
Published: 04 January 2023.
Edited by:
Zili Zhang, Nanjing University of Chinese Medicine, ChinaReviewed by:
Yang Wu, Nanjing Medical University, ChinaChun Zhang, Ludwig Maximilian University of Munich, Germany
Zhuang Yanshuang, Taizhou Traditional Chinese Medicine Hospital, China
Copyright © 2023 Kuang, Jiang, Yu, Hu, Duan and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zhong Chen, Y2hlbnpob25nQHNtdS5lZHUuY24=
†These authors have contributed equally to this work