Wei Guo1,2,3
Lin-Feng Wang
Xiao-Ming Li- 1Department of Orthopaedics, The Third Hospital of Hebei Medical University, Shijiazhuang, China
- 2Department of Orthopaedics, Hebei Province Cangzhou Hospital of Integrated Traditional Chinese Medicine-Western Medicine, Cangzhou, China
- 3Hebei Key Laboratory of Integrated Traditional and Western Medicine in Osteoarthrosis Research, Cangzhou, China
- 4Department of Breast Surgery, Hebei Province Cangzhou Hospital of Integrated Traditional Chinese Medicine-Western Medicine, Cangzhou, China
Intervertebral disc (IVD) degeneration and its inflammatory microenvironment can result in discogenic pain, which has been shown to stem from the nucleus pulposus (NP). Increasing evidence suggests that mitochondrial related genes are strictly connected to cell functionality and, importantly, it can regulate cell immune activity in response to damaged associated signals. Therefore, identification of mitochondria related genes might offer new diagnostic markers and therapeutic targets for IVD degeneration. In this study, we identified key genes involved in NP tissue immune cell infiltration during IVD degeneration by bioinformatic analysis. The key modules were screened by weighted gene co-expression network analysis (WCGNA). Characteristic genes were identified by random forest analysis. Then gene set enrichment analysis (GSEA) was used to explore the signaling pathways associated with the signature genes. Subsequently, CIBERSORT was used to classify the infiltration of immune cells. Function of the hub gene was confirmed by PCR, Western blotting and ELISA. Finally, we identified MFN2 as a crucial molecule in the process of NP cell pyroptosis and NLRP3 inflammasome activation. We speculate that the increased MFN2 expression in NP tissue along with the infiltration of CD8+ T cells, NK cell and neutrophils play important roles in the pathogenesis of IVD degeneration.
Introduction
Low back pain (LBP) is a major musculoskeletal disease with adversely affects in people of all ages and socioeconomic groups (Knezevic et al., 2021). It is the most common reason for medical consultation and the leading reason of disability worldwide (Knezevic et al., 2021). LBP can be caused by a variety of reasons, but intervertebral disc (IVD) degeneration has been indicated as the most important reason (Martin et al., 2008; Livshits et al., 2011). IVD is composed of nucleus pulposus (NP) located centrally, fibrous annulus in the periphery, and cartilaginous end plates that connect cranially and caudally (Colombier et al., 2014). The healthy NP tissue serves as the hydrogel-like core of the IVD, consisting mainly of NP cells and extracellular matrix (ECM) (Lyu et al., 2019). It is generally believed that NP cells are of great important for IVD function and ECM metabolism, which first exhibits degenerative changes during IVD degeneration (Lan et al., 2021). Therefore, a further study of the molecular mechanism of NP cell loss may provide a novel therapeutic target for improving the IVD degeneration.
Mitochondria is a closed, double-membrane organelle which is found in nearly all eukaryotes. Mitochondria generate adenosine triphosphate mainly via oxidative phosphorylation of substances. Mitochondrial membrane potential and intimal permeability are the main characteristics of healthy mitochondria. In order to maintain oxidative phosphorylation, various carbon substrates need to be metabolized through certain pathways, and finally recycle and polymerize with tricarboxylic acid (TCA) to produce reductive equivalent. Mitochondria contain an active calcium ion transport system that activates many enzymes associated with oxidative metabolic pathways. Previous study has linked mitochondria to cell damage and a wide range of age-related diseases (Giorgi et al., 2018). Cell death induced by mitochondria is an important mechanism leading to IVD degeneration (Sun et al., 2021). It has been confirmed that disruption of mitochondrial dynamics is also closely related to mitochondrial dysfunction and oxidative stress during IVD degeneration. Xu et al. (2019) found that the accumulation of progerin in human NP tissues was correlated with IVD degeneration progression, and further studies confirmed that progerin stimulation could shift mitochondrial dynamics to fission events by reducing the levels of mitochondrial fusion factors Opa1 and increasing the levels of mitochondrial fission factor Drp1 (Song et al., 2021). NP cell function loss occurs in the IVD degeneration progression stage due to cell senescence, death, inflammatory responses, and imbalances in anabolic and catabolic metabolism, which are also closely associated with mitochondrial damage (Wang et al., 2020a; Huang et al., 2020). Most studies on mitochondria-related genes focus on their regulation of mitochondrial function (Prasun, 2020; Sun et al., 2021). However, increasing evidence suggests that mitochondrial related genes are strictly connected to cell functionality and, importantly, it can regulate cell immune activity in response to damaged associated signals (Zhang et al., 2021; Deng et al., 2022). IVD was considered as an immunologically privileged organ that isolated NP tissue from the host immune system due to their unique structure. When the blood-NP barrier is breached, NP tissue triggers an immune response. This process plays an important role in IVD degeneration and result in a variety of subsequent pathological processes (Bridgen et al., 2017). Therefore, the relationship of mitochondrial related genes with NP cell inflammatory response in the process of IVD degeneration has become a research hotspot.
In this study, bioinformatics analysis based on transcriptome microarrays was used to identify mitochondria-associated genes in patients with IVD degeneration, providing new insights into the diagnosis and treatment of the disease. Moreover, we investigated the effects of Mitofusin-2 (MFN2) overexpression on human NP cells pyroptosis and NLRP3 inflammasome activation.
Materials and methods
Ethics statement
This study was approved by the ethics committees of Hebei Province Cangzhou Hospital of Integrated Traditional Chinese Medicine-Western Medicine (2020022). Human NP tissue samples were obtained from patients undergoing surgery at Hebei Province Cangzhou Hospital of Integrated Traditional Chinese Medicine-Western Medicine. Written informed consent was obtained from all patients for the use of their tissue specimens for research purpose.
Clinical specimens
Human lumbar degenerative NP specimens were obtained from 10 patients with IVD degeneration undergoing discectomy. The control samples were taken from 10 age- and sex-matched patients with fresh traumatic vertebral fracture undergoing decompressive surgery because of neurological deficits.
Data sources
In this study, two data sets were downloaded from the Gene Expression Omnibus (GEO) (http://www.ncbi.nlm.nih.gov/geo/), namely, GSE56081 [platform: GPL15314 Arraystar Human LncRNA microarray V2.0 (Agilent-033010 probe name version) and GSE165722 (GPL24676 Illumina NovaSeq) 6,000 (Homo sapiens)]. GSE56081 contained two groups of mRNA expression profiles, containing 5 NP tissue from IVD degeneration patients and five from controls. GSE165722 contained eight samples of Single-cell transcriptome profiling from NP tissue. A total of 1,136 mitochondria-related genes (Mito-Genes) were downloaded from the Mitocarta3.0 database (Rath et al., 2021).
Weighted gene coexpression network analysis (WGCNA) was used to identify IVD degeneration genes
For this study, a R package named “WGCNA” was used to construct the co-expressed clusters by analysis the gene expression level of GSE56081 (Langfelder and Horvath, 2008). Prior to the network construction, the raw data was converted into a recognizable format using R package “affy,” then the missing values were estimated using k nearest neighbor based approach, and then the median method was used for normalization (Troyanskaya et al., 2001). The pickSoftThreshold function of WGCNA package is used to calculate the soft threshold power and adjacency relationship. Then, the adjacency matrix is converted into topological overlap matrix, and the corresponding dissimilarity is calculated to carry out hierarchical clustering analysis. The dynamic tree cutting method with a minimum module size of 20 were used to identify the co-expressed gene modules. The module associated with IVD degeneration was identified by correlation analysis. IVD degeneration genes were acquired by module membership within the modules.
Identification of IVD degeneration-mitochondria related genes
The IVD degeneration-mitochondria related genes (IVDD-Mito_Genes) were obtained by intersecting the IVD degeneration genes and Mito-Genes. We then used R package clusterProfiler to analyze the GO and KEGG signaling pathways of IVDD-Mito_Genes and output the top 10 GO and 15 KEGG signaling pathways.
Signature genes identification
Subsequently, we exerted machine learning algorithms called random forest to screen IVDD-Mito_Genes. We used R packet “randomforest” to classify IVDD-Mito_Genes. The random forest model calculated the average error rate of IVDD-Mito_Genes to determine the optimal number of variables (Xie et al., 2022; Yu et al., 2022). We then calculate the error rate for each tree and determine the optimal number of trees based on the lowest error rate. After the above parameters are determined, the random forest tree model is established. Finally, the characteristic importance score of each IVDD-Mito_Genes was performed, and the top10 genes were selected. We used the area under the curve (AUC) of the receiver operating characteristic curve (ROC) to evaluate the diagnostic efficiency of these signature genes. An AUC greater than 0.8 indicates a good diagnostic effect.
Gene set enrichment analysis
We grouped the IVD degeneration cohort according to the median expression of signature genes, and performed gene set enrichment analysis (GSEA) according to different subgroups to determine the relationship between signature genes and signaling pathways (Subramanian et al., 2005).
Immune cell infiltration analysis
We used the principle of linear support vector regression to deconvolution the expression matrix of 22 human immune cell subtypes by CIBERSORT method to investigate the differences of immune cell between IVD degeneration patients and normal subjects. Subsequently, immune cells with significant differences in infiltration between patients with IVD degeneration and normal subjects were screened, and spearman correlation analysis was used to identify their correlation with signature genes.
Signal cell RNA-seq data analysis
Raw RNA-seq data (GSE165722) were processed according to bioinformatics analysis principles. Cells with a unique feature count of <200 and a mitochondrial count of >20% in the sample were filtered out. t-distributed stochastic neighbor embedding (t-SNE) analysis, K-mean clustering and hierarchical clustering methods were used to analyze the data. The following R packages were used: limma, Seurat, dplyr, magrittr, celldex, SingleR, monocle.
Cell communication analysis
Cell communication is determined by assessing the expression of ligand and receptor pairs within a cell population (CellChat R package). Interactions between different cell types were examined, and gene expression 0.2 was set as an effective cut-off point.
Human NP cell culture
The NP tissue samples were separated and cut into 1 m pieces, and the nucleus pulposum tissue was digested with 0.25% pronase and 0.2% collagenase type II at 37°C. The digested suspension was filtered through a 70 μm pore mesh and cultured in DMEM medium containing 10% fetal bovine serum and 1% penicillin-streptomycin at a CO2 concentration of 5% and a temperature of 37°C.
Western blotting
Cells were lysed with a buffer containing a protease inhibitor, 0.25 M Tris-HCl, 20% glycerol, 4% sodium dodecyl sulfate (SDS), and 10% mercaptoethanol (pH 6.8). The same amount of total protein (10 μg) was separated by 10%–12% SDS-polyacrylamide gel and electrotransferred to polyvinylidene fluoride membrane. Then, 5% non-fat milk in Tris-buffered saline containing 0.1% Tween-20 (TBST) was employed for blocking at ambient (1 h), followed by incubation with primary antibodies in TBST containing 5% non-fat milk overnight at 4°C. Secondary antibody was added at room temperature for 1 h and Western blotting was performed using an enhanced chemiluminescence system.
ELISA
The contents of human IL-1β in the cell culture supernatant were detected according to the manufacturer’s instructions.
Quantitative reverse-transcription PCR (RT-qPCR)
After chloroform extraction and precipitation, the DNA in the sample was removed by DNase I treatment, followed by reverse transcription of purified RNA using RevertAid reverse transcriptase. The reverse transcription primers containing CMV promoter sequences were designed to target specific genes. RT-qPCR was performed with AriaMx Real-time PCR system and QuantStudio Real-time PCR system. Supplementary Table S3 showed the primers used in this study.
Cell transfection
Lipofectamine 3,000 (Invitrogen) was used to transfect the plasmids or siRNA respectively with third-generation NP cells as recommended by the manufacturer and described previously (Guo et al., 2020). Cells were collected 48 h after transfection.
Statistical analysis
All experiments were repeated three times or more. Continuous data are mean ± standard deviation (SD). Bioinformatic analyses in this study was performed by R software (version 4.1.2). Prism 7.0 (GraphPad Software, United States) and SPSS 22.0 (SPSS Inc., United States) were used for statistical analyses. *p < 0.05, **p < 0.01, and ***p < 0.001 were significance levels.
Results
Identification of IVD degeneration genes by WGCNA analysis
We use dynamic shear tree algorithm to segment the modules, and set the minimum size of the modules as 20 to build a scale-free co-representation network. Figure 1A shows the results of the cluster tree. Finally, the data were clustered into 19 modules (Figure 1B). Then the correlation between each module and IVD degradation is calculated. Supplementary Table S1 shows the genes in each module. Correlation analysis showed that blue module was significantly associated with IVD degradation. Therefore, genes in blue module was considered as IVD degradation genes. The crossover between the IVD degeneration genes and mitochondrial related genes was shown in Figure 1C (Supplementary Table S2).
FIGURE 1. The gene modules significantly associated with IVD degeneration were identified by WGCNA. (A) WGCNA cluster dendrogram. (B) WGCNA cluster modules. (C) Venn diagram shows the interaction between mitochondria-associated genes and genes associated with IVD degradation.
Function enrichment analysis
By GO and KEGG enrichment analysis of 58 IVD degeneration-Mito_genes, 81 biological processes (BP), 50 cellular components (CCs), 32 molecular functions (mf) and 20 KEGG signaling pathways were identified. The enrichment analysis results of the top 10 GO and 10 KEGG signaling pathways are shown in Figures 2A, B, respectively. The results showed that these genes are enriched in biological processes such as protein targeting to mitochondrial, mitochondrial transport, mitochondrial protein localization, protein transmembrane transport, and KEGG signaling pathways such as those in Parkinson’s disease, and oxidative phosphorylation.
FIGURE 2. Functional enrichment analysis of IVD degeneration-Mito_Genes. (A) Top 10 BP, MF, CC of GO functional enrichment analysis. (B) The KEGG analysis of degeneration-Mito_Genes.
Identification of signature genes by random forest algorithms
Eight signature genes with relative importance greater than 0.3 were identified by random forest analysis (Figures 3A, B), containing MRPS17, ECHDC2, MFN2, EFHD1, PDF, PXMP2, TOMM7, PANK2.
FIGURE 3. The machine algorithms for signature genes. (A) Confidence intervals for error rates of random forest models. (B) The relative importance of genes in random forest models.
Diagnostic efficacy of signature genes in predicting IVD degeneration
Figures 4A–H showed that four of the signature genes screened through random forest were more highly expressed in patients with IVD degeneration than normal people, including MRPS17, MFN2, PDF and PANK2, suggesting that these genes may play a key role in IVD degeneration. Figures 4I–P shows the area under curve (AUC) of these signature genes. The results showed that three of the signature genes (MRPS17, MFN2, PDF) had significant diagnostic effect in predicting the degeneration of IVD. These results were further confirmed in 10 patients and 10 controls human NP samples using RT-qPCR analysis and Western blotting, MRPS17, MFN2, PDF were significantly increased in IVD degeneration NP tissues, MFN2 showed significant diagnostic effect in predicting IVD degeneration (Supplementary Figure S1).
FIGURE 4. The performance of the signature genes. (A–H) Signature gene expression levels in the nucleus pulposus tissue of IVD degeneration patient and healthy subjects. (I–P) The diagnostic performance of signature genes was evaluated by ROC.
Gene set enrichment analysis
We analysis signaling pathways and GO associated with signature genes which had diagnostic efficiency by GSEA analysis. Figures 5A–C demonstrated the top six signal pathways. The results showed that MRPS17 was markedly associated with adherens junction, ether lipid metabolism, biosynthesis keratan sulfate, hematopoietic cell lineage, protein export, RNA degradation (Figure 5A). The expression of MFN2 significantly correlated with ether lipid metabolism, biosynthesis keratan sulfate, hematopoietic cell lineage, olfactory transduction, protein export and RNA degradation (Figure 5B). The expression of PDF significantly correlated with allograft rejection, autoimmune thyroid disease, graft versus host disease, propanoate metabolism and protein export (Figure 5C). The top six GO were demonstrated in Figures 5D–F. The results showed that MRPS17 was significantly correlated with myoblast fusion, positive regulation of monooxygenase activity, regulation of B cell mediated immunity, regulation of immunoglobulin production, regulation of muscle organ development and cytokine activity (Figure 5D). The expression of MFN2 significantly correlated with myoblast fusion, positive regulation of monooxygenase activity, regulation of B cell mediated immunity, regulation of immunoglobulin production, regulation of muscle organ development and cytokine activity (Figure 5E). The expression of PDF significantly correlated with kinetochore organization, neutral amino acid transport, pigment biosynthetic process, positive regulation of endothelial cell proliferation, regulation of extrinsic apoptotic signaling pathway via death domain receptors (Figure 5F).
FIGURE 5. The GSEA of the signature genes in IVD degeneration. (A–C) KEGG enrichment results of signature genes (MRPS17, MFN2, PDF). (D–F) GO enrichment results of signature genes (MRPS17, MFN2, PDF).
Immune cell infiltration
Compared with normal subjects, IVD degeneration patients have higher CD8+ T cells, NK cell activated, neutrophils infiltration and lower dendritic cell resting (Figure 6A). MFN2 was associated with the infiltration of CD8+ T cells, NK cell activated and neutrophils positively. PDF was associated with the infiltration of CD8+ T cells positively. MRPS17 was positively associated with the infiltration of CD8+ T cells and neutrophils (Figure 6B). Notably, MFN2 was positively correlated with all the IVD degeneration related immune cell infiltration. These results demonstrated that MFN2 might be a potential immune-related target, and further study of MFN2 might provide a better understanding of immune cell infiltration in IVD degeneration.
FIGURE 6. Signature genes related to immune cell infiltration. (A) The infiltration of immune cells in the NP tissues of the IVD degeneration and healthy subjects. (B) The relationship between signature genes and immune cell infiltration.
Single-cell transcriptome profiling of the IVD tissue
We used scRNA-seq data (GSE165722) to analyze and identify different cell types in NP tissues during IVD degeneration, and visualized the results with dimensionality reduction algorithm (t-SNE). We identified 21 cluster and 11 major distinct cell types which were determined using singleR and cell markers (Figures 7A, B). Consistent with above results, we found a lot of immune cells in the tSNE map which may suggest the immune cell infiltration during the progress of IVD degeneration. Subsequently, we examined MFN2 expression in these cell types. As expected, MFN2 was expressed in NP cells and most immune cells (Figures 7C, D).
FIGURE 7. Single cell sequencing data analysis and identification of different cell types. (A, B) t-SNE projections and cell type annotation of IVD tissue. (C, D) t-SNE and violin diagrams demonstrated MFN2 expression pattern at the single-cell level.
We used CellChat analysis to identified potential interactions between different types of immune cells and NP cells. The results showed that NP cells interact with neutrophils, myelocyte, monocyte, macrophage, endothelial cells closely (Figures 8A–K).
FIGURE 8. CellChat analysis. (A–K) Circle plots demonstrated the interactions in the cell-cell communication network between IVD degeneration and healthy subjects, respectively.
Mfn2 overexpression activated NLRP3 inflammasome
IL-1β secretion following exposure to MSU (a prototypical danger signals) was significantly decreased in MFN2 knockdown NP cells as compared to wild-type NP cells (Figure 9A). As shown in Figures 9B, C, MSU increased MFN2 mRNA and protein expression level in comparison with control group. Pyrosis is a programmed inflammatory cell death that occurs when inflammasome is activated. In this study, we found that the absence of MFN2 markedly attenuated the release of LDH from NP cells, suggesting that MFN2 had a regulatory effect on pyroptosis (Figure 9D). To elucidate the effect of MFN2 in the activation of NLRP3 inflammasome, we examined IL-1β cleavage and caspase-1 activation in MFN2 overexpressed NP cells. The results showed that the activation of Caspase-1 and the cleavage of IL-1β were significantly increased in MFN2 overexpressed NP cells (Figure 9E). Notably, MFN2 overexpression significantly increased NLRP3 protein and mRNA expression level, but had no effect on pro-IL-1β and pro-caspase-1 expression (Figures 9E, F). And NLRP3 knockdown markedly inhibited IL-1β secretion in MFN2 overexpression NP cells (Figure 9G). These results indicated that MFN2 can promote the expression of NLPR3 and activate the NLRP3 inflammasome, resulting in the pyroptosis of NP cells and the release of inflammatory cytokines.
FIGURE 9. Overexpression of MFN2 activated NLRP3 inflammasome in NP cells. (A) NP cells were transfected as indicated. IL-1β released was assessed by ELISA. (B, C) NP cells were stimulated with MSU. MFN2 mRNA expression was assessed by RT-qPCR and Western blot. (D) Pyroptosis was determined by LDH assay. (E) The Western blot assay for NLRP3, IL-1β cleavage and caspase-1 activation. (F) NP cells were transfected as indicated. NLRP3, pro- IL-1β and pro-caspase-1 mRNA were assessed by RT-qPCR. (G) NP cells were transfected as indicated. IL-1β released was assessed by ELISA.
Discussion
MFN2 is a multifunctional protein involved in a variety of physiological and pathological processes, including mitochondrial fusion transport, cell metabolism, pyroptssis, and autophagy (Yu et al., 2018). Dysregulation of MFN2 can lead to neurodegeneration, sarcopenia, metabolic disease, heart disease and many other diseases (Sebastian et al., 2016; Kim et al., 2017; Chandhok et al., 2018; Rocha et al., 2018). Increasing evidence have demonstrated that immune cell infiltration acts as an important role in IVD degeneration (Wang et al., 2021). This study is the first time to reveal the relationship between MFN2 and immune cell infiltration during intervertebral disc degeneration.
Previous studies have shown that MFN2 can be phosphorylated by Jun N-terminal kinase under cellular stress, leading to MFN2 degradation through the ubiquitin-proteasome system (Leboucher et al., 2012). The ubiquitination of MFN2 promotes mitochondrial degradation and prevents fusion of damaged mitochondria (Tanaka et al., 2010). In this study, our results indicated that MFN2 has the ability to regulate immune response in addition to the function of regulating mitochondrial fusion. Previously research demonstrated that IVD degeneration is characterized by infiltration of CD68+ macrophages, T cells (CD4+, CD8+), and neutrophils (Wang et al., 2021). In combination with scRNA-seq, cell heterogeneity can be identified for in-depth study of biological structure and function. The scRNA-seq technique can generate the expression profile of individual cells for the analysis of heterogeneous cell populations and the identification of cell types. Understanding the phenotype of immune cells in the IVD microenvironment is critical to understanding the mechanisms of IVD degeneration progression. Immune infiltrating cells in the IVD degeneration have been shown to be very important in the IVD degeneration (Risbud and Shapiro, 2014; Wang et al., 2021). In this research, we found significantly increased infiltration of CD8+ T cells and neutrophils in degenerative NP tissues compared to controls. We also found a strong positive correlation between the MFN2 and the level of immune infiltration of three types of invasive immune cells (including CD8+ T cells, NK cell and neutrophils). This implied a strong association of MFN2 with immune cells in the IVD degeneration, which was subsequently validated using scRNA-seq datasets. CellChat can quantitatively infer and analyze intercellular communication networks from scRNA-seq data. Our results showed the number of interactions between NP cells and immune cells in the IVD degeneration including neutrophils, myelocyte, monocyte, macrophage. IVD is considered as an immunologically privileged organ because their unique structure insulates the NP tissue from the host immune system (Sun et al., 2020). When NP cells pyroptosis, NP releases inflammatory factors such as IL-1β, which induces immune cell infiltration and triggers an immune response. This process act as an important role in IVD degeneration and result in multiple pathological processes that eventually lead to fibrotic of NP tissue (Wang et al., 2020b). At the same time, immune cells recruited into the IVD microenvironment will release more inflammatory mediators, further damaging the IVD microenvironment and causing nucleus pulposus cell death (Silva et al., 2019).
NLRP3 is a member of the pattern recognition receptors that mediate the activation of the innate immune system. The NLRP3 inflammasome is important for the immune defense system. NLRP3 initiates inflammasome assembly, promotes recruitment of inflammasome complex by caspase-1 and activates caspase-1. Activated caspase-1 cleave the IL-1β precursor protein and converts it into a biologically active mature form, which promotes an immune response. Previous studies have shown that MFN2 acts as a major regulator of the immune response by binding to NLRP3 and promoting IL-1β secretion after infection with the virus (Ichinohe et al., 2013; Tur et al., 2020). In our research, we observed that MFN2 was highly expressed during IVD degeneration. MFN2 overexpression associated with the infiltration of CD8+ T cells, NK cell activated and neutrophils. MFN2 overexpression also activate NLRP3 inflammasome which results in IL-1β release and NP cell pyroptosis. This may be due to the overexpression of MFN2 in NP cells inducing the expression of NLRP3 and the activation of inflammasome. Combining the results of this study with the conclusions of previous studies, we hypothesized that a special type of immune microenvironment is generated during the process of IVD degeneration by release of IL-1β from NP cells that recruit CD8+ T cells, NK cell activated and neutrophils. This immune microenvironment, in turn, further promotes NP cell death and IVD pathological changes.
This study also has some limitations. First of all, the data in this study are mainly from public databases and involve a small number of samples, which may lead to certain bias. However, the reliability of our analysis was confirmed by the results of in vitro experiments. Second, this study still needs a larger clinical samples to further validate the results. Third, the specific mechanism of MFN2 cause NLRP3 induction and increase inflammasome activation remains to be explored.
Our research shows that immune cell infiltration, including T cells, NK cell activated and neutrophils, is participated in the pathological process of IVD degeneration. Via microarray data analysis, single-cell sequencing data analysis and in vitro experiments, we identified MFN2 as a signature gene, which present good diagnostic value and may serve as a molecular target for treatment of IVD degeneration.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving human participants were reviewed and approved by Ethics committees of Hebei Province Cangzhou Hospital of Integrated Traditional Chinese Medicine-Western Medicine. The patients/participants provided their written informed consent to participate in this study.
Author contributions
WG designed the study and analysis. Experiments were performed by WG and W-SL. WG and KM analyzed data and typeset figures. The study was supervised by L-FW, X-ML, and J-YZ. WG wrote the manuscript with contributions from KM, W-SL, L-FW, X-ML, and J-YZ. All of the authors subsequently reviewed and edited the manuscript.
Funding
This study was supported by the National Natural Science Foundation of China (82002365) and Natural Science Foundation of Hebei (H2019110028).
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.2023.1135767/full#supplementary-material
SUPPLEMENTARY FIGURE S1 | (A) The mRNA expression level of signature genes in NP tissues was measured in 10 patients and 10 controls using qRT-PCR (**P < 0.01; ***P < 0.001). (B–I) The diagnostic performance of signature genes was evaluated by ROC. (J–Q) The protein expression level of signature genes in NP tissue were measured by Western blotting.
SUPPLEMENTARY FIGURE S2 | The efficiency of MFN2 knock down in NP cells.
References
Bridgen, D. T., Fearing, B. V., Jing, L., Sanchez-Adams, J., Cohan, M. C., Guilak, F., et al. (2017). Regulation of human nucleus pulposus cells by peptide-coupled substrates. Acta Biomater. 55, 100–108. doi:10.1016/j.actbio.2017.04.019
Chandhok, G., Lazarou, M., and Neumann, B. (2018). Structure, function, and regulation of mitofusin-2 in health and disease. Biol. Rev. Camb Philos. Soc. 93 (2), 933–949. doi:10.1111/brv.12378
Colombier, P., Clouet, J., Hamel, O., Lescaudron, L., and Guicheux, J. (2014). The lumbar intervertebral disc: From embryonic development to degeneration. Jt. Bone Spine 81 (2), 125–129. doi:10.1016/j.jbspin.2013.07.012
Deng, H., Zhu, S., Zhu, L., Sun, J., Ding, Y., Li, F., et al. (2022). Mfn2 is responsible for inhibition of the RIG-I/IRF7 pathway and activation of NLRP3 inflammasome in Seneca Valley virus-infected PK-15 cells to promote viral replication. Front. Immunol. 13, 955671. doi:10.3389/fimmu.2022.955671
Giorgi, C., Marchi, S., Simoes, I. C. M., Ren, Z., Morciano, G., Perrone, M., et al. (2018). Mitochondria and reactive oxygen species in aging and age-related diseases. Int. Rev. Cell Mol. Biol. 340, 209–344. doi:10.1016/bs.ircmb.2018.05.006
Guo, W., Mu, K., Zhang, B., Sun, C., Zhao, L., Li, H. R., et al. (2020). The circular RNA circ-GRB10 participates in the molecular circuitry inhibiting human intervertebral disc degeneration. Cell Death Dis. 11 (8), 612. doi:10.1038/s41419-020-02882-3
Huang, D., Peng, Y., Li, Z., Chen, S., Deng, X., Shao, Z., et al. (2020). Compression-induced senescence of nucleus pulposus cells by promoting mitophagy activation via the PINK1/PARKIN pathway. J. Cell Mol. Med. 24 (10), 5850–5864. doi:10.1111/jcmm.15256
Ichinohe, T., Yamazaki, T., Koshiba, T., and Yanagi, Y. (2013). Mitochondrial protein mitofusin 2 is required for NLRP3 inflammasome activation after RNA virus infection. Proc. Natl. Acad. Sci. U. S. A. 110 (44), 17963–17968. doi:10.1073/pnas.1312571110
Kim, Y. J., Park, J. K., Kang, W. S., Kim, S. K., Han, C., Na, H. R., et al. (2017). Association between mitofusin 2 gene polymorphisms and late-onset alzheimer's disease in the Korean population. Psychiatry Investig. 14 (1), 81–85. doi:10.4306/pi.2017.14.1.81
Knezevic, N. N., Candido, K. D., Vlaeyen, J. W. S., Van Zundert, J., and Cohen, S. P. (2021). Low back pain. Lancet 398 (10294), 78–92. doi:10.1016/S0140-6736(21)00733-9
Lan, T., Shiyu, H., Shen, Z., Yan, B., and Chen, J. (2021). New insights into the interplay between miRNAs and autophagy in the aging of intervertebral discs. Ageing Res. Rev. 65, 101227. doi:10.1016/j.arr.2020.101227
Langfelder, P., and Horvath, S. (2008). WGCNA: an R package for weighted correlation network analysis. BMC Bioinforma. 9, 559. doi:10.1186/1471-2105-9-559
Leboucher, G. P., Tsai, Y. C., Yang, M., Shaw, K. C., Zhou, M., Veenstra, T. D., et al. (2012). Stress-induced phosphorylation and proteasomal degradation of mitofusin 2 facilitates mitochondrial fragmentation and apoptosis. Mol. Cell 47 (4), 547–557. doi:10.1016/j.molcel.2012.05.041
Livshits, G., Popham, M., Malkin, I., Sambrook, P. N., Macgregor, A. J., Spector, T., et al. (2011). Lumbar disc degeneration and genetic factors are the main risk factors for low back pain in women: The UK twin spine study. Ann. rheumatic Dis. 70 (10), 1740–1745. doi:10.1136/ard.2010.137836
Lyu, F. J., Cheung, K. M., Zheng, Z., Wang, H., Sakai, D., and Leung, V. Y. (2019). IVD progenitor cells: A new horizon for understanding disc homeostasis and repair. Nat. Rev. Rheumatol. 15 (2), 102–112. doi:10.1038/s41584-018-0154-x
Martin, B. I., Deyo, R. A., Mirza, S. K., Turner, J. A., Comstock, B. A., Hollingworth, W., et al. (2008). Expenditures and health status among adults with back and neck problems. JAMA 299 (6), 656–664. doi:10.1001/jama.299.6.656
Prasun, P. (2020). Mitochondrial dysfunction in metabolic syndrome. Biochim. Biophys. Acta Mol. Basis Dis. 1866 (10), 165838. doi:10.1016/j.bbadis.2020.165838
Rath, S., Sharma, R., Gupta, R., Ast, T., Chan, C., Durham, T. J., et al. (2021). MitoCarta3.0: An updated mitochondrial proteome now with sub-organelle localization and pathway annotations. Nucleic Acids Res. 49 (D1), D1541–D1547. doi:10.1093/nar/gkaa1011
Risbud, M. V., and Shapiro, I. M. (2014). Role of cytokines in intervertebral disc degeneration: Pain and disc content. Nat. Rev. Rheumatol. 10 (1), 44–56. doi:10.1038/nrrheum.2013.160
Rocha, A. G., Franco, A., Krezel, A. M., Rumsey, J. M., Alberti, J. M., Knight, W. C., et al. (2018). MFN2 agonists reverse mitochondrial defects in preclinical models of Charcot-Marie-Tooth disease type 2A. Science 360 (6386), 336–341. doi:10.1126/science.aao1785
Sebastian, D., Sorianello, E., Segales, J., Irazoki, A., Ruiz-Bonilla, V., Sala, D., et al. (2016). Mfn2 deficiency links age-related sarcopenia and impaired autophagy to activation of an adaptive mitophagy pathway. EMBO J. 35 (15), 1677–1693. doi:10.15252/embj.201593084
Silva, A. J., Ferreira, J. R., Cunha, C., Corte-Real, J. V., Bessa-Goncalves, M., Barbosa, M. A., et al. (2019). Macrophages down-regulate gene expression of intervertebral disc degenerative markers under a pro-inflammatory microenvironment. Front. Immunol. 10, 1508. doi:10.3389/fimmu.2019.01508
Song, Y., Lu, S., Geng, W., Feng, X., Luo, R., Li, G., et al. (2021). Mitochondrial quality control in intervertebral disc degeneration. Exp. Mol. Med. 53 (7), 1124–1133. doi:10.1038/s12276-021-00650-7
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. U. S. A. 102 (43), 15545–15550. doi:10.1073/pnas.0506580102
Sun, K., Jing, X., Guo, J., Yao, X., and Guo, F. (2021). Mitophagy in degenerative joint diseases. Autophagy 17 (9), 2082–2092. doi:10.1080/15548627.2020.1822097
Sun, Z., Liu, B., and Luo, Z. J. (2020). The immune privilege of the intervertebral disc: Implications for intervertebral disc degeneration treatment. Int. J. Med. Sci. 17 (5), 685–692. doi:10.7150/ijms.42238
Tanaka, A., Cleland, M. M., Xu, S., Narendra, D. P., Suen, D. F., Karbowski, M., et al. (2010). Proteasome and p97 mediate mitophagy and degradation of mitofusins induced by Parkin. J. Cell Biol. 191 (7), 1367–1380. doi:10.1083/jcb.201007013
Troyanskaya, O., Cantor, M., Sherlock, G., Brown, P., Hastie, T., Tibshirani, R., et al. (2001). Missing value estimation methods for DNA microarrays. Bioinforma. Oxf. Engl. 17 (6), 520–525. doi:10.1093/bioinformatics/17.6.520
Tur, J., Pereira-Lopes, S., Vico, T., Marin, E. A., Munoz, J. P., Hernandez-Alvarez, M., et al. (2020). Mitofusin 2 in macrophages links mitochondrial ROS production, cytokine release, phagocytosis, autophagy, and bactericidal activity. Cell Rep. 32 (8), 108079. doi:10.1016/j.celrep.2020.108079
Wang, L., He, T., Liu, J., Tai, J., Wang, B., Zhang, L., et al. (2021). Revealing the immune infiltration landscape and identifying diagnostic biomarkers for lumbar disc herniation. Front. Immunol. 12, 666355. doi:10.3389/fimmu.2021.666355
Wang, Y., Che, M., Xin, J., Zheng, Z., Li, J., and Zhang, S. (2020). The role of IL-1β and TNF-α in intervertebral disc degeneration. Biomed. Pharmacother. 131, 110660. doi:10.1016/j.biopha.2020.110660
Wang, Y., Wang, H., Zhuo, Y., Hu, Y., Zhang, Z., Ye, J., et al. (2020). SIRT1 alleviates high-magnitude compression-induced senescence in nucleus pulposus cells via PINK1-dependent mitophagy. Aging (Albany NY) 12 (16), 16126–16141. doi:10.18632/aging.103587
Xie, X., Zhang, Y., Yu, J., Jiang, F., and Wu, C. (2022). Significance of m(6)A regulatory factor in gene expression and immune function of osteoarthritis. Front. Physiol. 13, 918270. doi:10.3389/fphys.2022.918270
Xu, X., Wang, D., Zheng, C., Gao, B., Fan, J., Cheng, P., et al. (2019). Progerin accumulation in nucleus pulposus cells impairs mitochondrial function and induces intervertebral disc degeneration and therapeutic effects of sulforaphane. Theranostics 9 (8), 2252–2267. doi:10.7150/thno.30658
Yu, F., Xu, T., Wang, M., Chang, W., Li, P., and Wang, J. (2018). Function and regulation of mitofusin 2 in cardiovascular physiology and pathology. Eur. J. Cell Biol. 97 (7), 474–482. doi:10.1016/j.ejcb.2018.07.003
Yu, J., Xie, X., Zhang, Y., Jiang, F., and Wu, C. (2022). Construction and analysis of a joint diagnosis model of random forest and artificial neural network for obesity. Front. Med. (Lausanne). 9, 906001. doi:10.3389/fmed.2022.906001
Keywords: MFN2, NLRP3 inflammasome, immune cells infiltration, intervertebral disc degeneration, signature genes
Citation: Guo W, Mu K, Li W-S, Gao S-X, Wang L-F, Li X-M and Zhao J-Y (2023) Identification of mitochondria-related key gene and association with immune cells infiltration in intervertebral disc degeneration. Front. Genet. 14:1135767. doi: 10.3389/fgene.2023.1135767
Received: 20 January 2023; Accepted: 27 February 2023;
Published: 08 March 2023.
Edited by:
Maurice H. T. Ling, Temasek Polytechnic, SingaporeReviewed by:
Feng Jiang, Fudan University, ChinaWanjing Shang, National Institutes of Health (NIH), United States
Copyright © 2023 Guo, Mu, Li, Gao, Wang, Li and Zhao. 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: Lin-Feng Wang, d2FuZ2xpbmZlbmdsYW9zaGlAMTYzLmNvbQ==; Xiao-Ming Li, bGl4aWFvbWluZ2N6QDEyNi5jb20=