- Department of Orthopedics, Tianjin Hospital, Tianjin, China
Background: Intervertebral disk degeneration (IVDD) is a major cause of low back pain and one of the most common health problems all over the world. However, the early diagnosis of IVDD is still restricted. The purpose of this study is to identify and validate the key characteristic gene of IVDD and analyze its correlation with immune cell infiltration.
Methods: 3 IVDD-related gene expression profiles were downloaded from the Gene Expression Omnibus database to screen for differentially expressed genes (DEGs). Gene Ontology (GO) and gene set enrichment analysis (GSEA) were conducted to explore the biological functions. Two machine learning algorithms were used to identify characteristic genes, which were tested to further find the key characteristic gene. The receiver operating characteristic curve was performed to estimate the clinical diagnostic value of the key characteristic gene. The excised human intervertebral disks were obtained, and the normal nucleus pulposus (NP) and degenerative NP were carefully separated and cultured in vitro. The expression of the key characteristic gene was validated by real-time quantitative PCR (qRT-PCR). The related protein expression in NP cells was detected by Western blot. Finally, the correlation was investigated between the key characteristic gene and immune cell infiltration.
Results: A total of 5 DEGs, including 3 upregulated genes and 2 downregulated genes, were screened between IVDD and control samples. GO enrichment analysis showed that DEGs were enriched to 4 items in BP, 6 items in CC, and 13 items in MF. They mainly included the regulation of ion transmembrane transport, transporter complex, and channel activity. GSEA suggested that the cell cycle, DNA replication, graft versus host disease, and nucleotide excision repair were enriched in control samples, while complement and coagulation cascades, Fc γ R–mediated phagocytosis, neuroactive ligand–receptor interaction, the NOD-like receptor signaling pathway, gap junctions, etc., were enriched in IVDD samples. Furthermore, ZNF542P was identified and tested as key characteristic gene in IVDD samples through machine learning algorithms and showed a good diagnostic value. The results of qRT-PCR showed that compared with normal NP cells, the expression of ZNF542P gene was decreased in degenerated NP cells. The results of Western blot suggested that compared with normal NP cells, the expression of NLRP3 and pro Caspase-1 was increased in degenerated NP cells. Finally, we found that the expression of ZNF542P was positively related to the proportions of T cells gamma delta (γδT cells).
Conclusion: ZNF542P is a potential biomarker in the early diagnosis of IVDD and may be associated with the NOD-like receptor signaling pathway and the infiltration of γδT cells.
1 Introduction
Clinically, low back pain is one of the most medical problems in spine surgery affecting 70%–85% of the population worldwide (1). Many causes could lead to low back pain, while intervertebral disk degeneration (IVDD) is considered as the most common reason (2). Currently, the main clinical therapeutic approaches regarding IVDD include conservative treatment and surgery. Symptomatic recurrence often occurs after conservative treatment, while surgery would destabilize the spine. Thus, IVDD has tremendous potential for improvement in early diagnostic aspects (3).
Bioinformatics is a discipline that studies the complex relationship of a large number of biological data, which is characterized by multidisciplinary interdisciplinary, using the internet as the medium and the database as the carrier. It takes advantage of various mathematical knowledge to establish different models and then uses the computer as a tool to store, retrieve, process, and analyze a large amount of biological data obtained from experiments and explains the results with biological knowledge (4). Machine learning is a multifield interdisciplinary subject that has risen in the past 20 years, involving the probability theory, statistics, approximation theory, convex analysis, algorithm complexity theory, and other disciplines. The machine learning algorithm is a kind of method that automatically analyzes and obtains rules from data and uses the rules to predict unknown data (5).
The purpose of our study is to use the RNA expression profiles obtained from the Gene Expression Omnibus (GEO) database to identify the key characteristic gene, validate the gene via an in vitro experiment, and analyze its correlation with immune cell infiltration in IVDD. In this study, we systematically analyzed gene expression profiles based on bioinformatics and machine learning algorithms between IVDD and healthy patients. In addition, we validated a key characteristic gene as an early diagnostic biomarker in IVDD.
2 Materials and methods
2.1 Dataset acquisition
3 IVDD-related gene expression profiles were downloaded from the GEO database, (https://www.ncbi.nlm.nih.gov/geo/). The GSE124272 (GPL21185) dataset included 8 patients with lumbar disk prolapse and 8 healthy volunteers. The GSE153761 (GPL22120) dataset included 3 patients with cervical spondylosis myelopathy and 3 healthy subjects. The GSE150408 (GPL21185) dataset including 17 patients with sciatica and 17 healthy volunteers was selected as a validation set. The GSE124272 and GSE153761 datasets were selected as a training set and to screen for differentially expressed genes (DEGs) between IVDD and control samples shown in Table 1.
2.2 Gene expression profile preprocessing
In this study, we downloaded 2 types of file for analysis: “platform” data and “series_matrix_file” data, among which “platform” data included GPL21185 and GPL22120, while “series_matrix_file” data contained GSE124272, GSE153761, and GSE150408. Platform data were used to probe the”series_matrix _file” data changing probe sets into gene symbols. Probe sets without corresponding gene symbols and repeated gene symbols were removed, and the median expression value was retained when multiple probe sets were mapped to the same gene.
2.3 Screening for differentially expressed genes
After the expression matrix of GSE124272 and GSE153761 datasets were merged, the data were further corrected and normalized. Genes between IVDD and control samples were screened by R software with the threshold of the adjusted p-value < 0.05 and ∣log2 fold change (FC)∣>0.5 were regarded as DEGs. Moreover, the heat map and volcano plot were plotted to visualize the date.
2.4 Gene Ontology and gene set enrichment analysis
Gene Ontology (GO) enrichment analysis was significantly related to the disease and could classify GO terms in 3 categories: biological processes (BP), cellular component (CC), and molecular function (MF). For gene set enrichment analysis (GSEA), it was performed to get more important biological functions or signaling pathways based on all genes and phenotype, which may be neglected by differential analysis. For both analyses, a p-value < 0.05 was considered as statistical significance.
2.5 Machine learning algorithms to identify characteristic genes
To further screen characteristic genes related to IVDD, least absolute shrinkage and selection operator (LASSO) logistic regression and support vector machine–recursive feature elimination (SVM-RFE) algorithms, two machine learning approaches, were combined to perform feature selection calculating the point with the smallest error in the training set (GSE124272 and GSE153761) by cross-validation. Next, the overlapping genes in the LASSO logistic regression and the SVM-RFE algorithm were regarded as characteristic genes for further analysis. Furthermore, to screen the key characteristic gene, we further tested the expression levels of characteristic genes in the validation set (GSE150408).
2.6 Estimation of clinical diagnostic value of key characteristic genes
To verify whether key characteristic genes could differentiate IVDD samples from control samples in the training set and validation set, the receiver operating characteristic (ROC) curves were plotted to estimate the clinical diagnostic value of key characteristic genes by calculating the area under the ROC curves (AUC).
2.7 Isolation and culture of human nucleus pulposus cells
The experimental protocol received approval from the Tianjin Hospital Ethics Committee. All patients participating in this study signed written informed consent allowing the researchers to use intervertebral disk tissues obtained from surgery. The intervertebral disk tissue was obtained from two patients (1 spinal trauma patient and 1 IVDD patient).
The NP tissue was cut into approximately 1 mm3 fragment, subsequently digested by 0.25% trypsin and 0.2% type II collagenase at 37°C, and filtered through a 100 μm nylon mesh. After isolation, NP cells were cultured in Dulbecco's Modified Eagle Medium (DMEM)/F12 medium containing 10% fetal bovine serum and penicillin at 37°C in a 5% CO2 incubator. The passage 3 (P3) cells were adopted for the following experiments.
2.8 Quantitative real-time polymerase chain reaction analysis
Total mRNA from NP cells was extracted with the TRIzol reagent (Invitrogen, USA), and cDNA was synthesized by the Rever Tra Ace qPCR RT Kit (TOYOBO, Japan). Real-time quantitative PCR (qRT-PCR) was performed with SYBR® Green Real-Time PCR Master Mix (TOYOBO, Japan). The GAPDH was used as a housekeeper gene. The ZNF542P and GAPDH PCR primer sequences were shown as follows: ZNF542P Forward (5’-3’) CTGCGAGTCCTCGCCTACT, Backward (5’-3’) CTGCCCAAGAAAGGAGAAACG. GAPDH Forward (5’-3’) GGTGGTCTCCTCTGACTTCAACAG, Backward (5’-3’) GTTGCTGTAGCCAAATTCGTTGTC.
2.9 Western blot
Total protein in NP cells was extracted with Radio Immunoprecipitation Assay (RIPA) lysis buffer, and the protein concentration was measured using the BCA Quantitation Kit (Boster, China). Proteins were separated by SDS-PAGE and transferred to polyvinylidene fluoride (PVDF) membranes. The membranes were blocked with 5% skim milk and then incubated with primary antibodies against NLRP3 (1:1,000; Abcam, UK), pro Caspase-1 (1:1,000; Abcam, UK), respectively, overnight at 4°C. The membranes were incubated with horseradish peroxidase–conjugated (HRP) secondary antibodies (1:1,000) for 1 h. Blots were developed using enhanced chemiluminescence (Santa Cruz Biotechnology, USA).
2.10 Correlation between key characteristic gene and immune cell infiltration
To further explore the correlation between the key characteristic gene and immune cell infiltration, we assessed the proportion and composition of infiltrating immune cells in IVDD and control samples in the training set by CIBERSORT analysis, a deconvolution algorithm based on a validated leukocyte signature matrix including 547 genes and 22 human immune cell subtypes (6). The difference between IVDD and control samples was shown using a boxplot and a violin plot. The correlation between the key characteristic gene and immune cell infiltration was shown using a scatter plot and a lollipop plot.
2.11 Statistical analysis
The statistical analysis and visualization of all data were processed using R language (version 4.1.2, based on corresponding R package). p-value < 0.05 was considered as statistically significant difference. Gene expression data were represented as mean ± standard deviation (±s). Statistical analysis was processed by the t-test using software SPSS 21.0.
3 Results
3.1 Data preprocessing and screening for differentially expressed genes
The GSE124272 dataset contains 16 samples (25,964 genes), and the GSE153761 dataset contains 6 samples (38,156 genes). After merging the data, the training set contains 22 samples (25,964 genes). After the standardization of the data, the training set contains 22 samples (22,833 genes). The validation set contains 34 samples (25,964 genes). A total of 5 DEGs were screened, including 3 upregulated genes (KCNH4, LGALSL, and ASMTL-AS1) and 2 downregulated genes (ZNF542P and TBC1D31). The heat map and volcano plot of DEGs expression were shown in Figures 1A, B.
Figure 1 The differentially expressed genes (DEGs) analysis of the training set. (A) The heat map of DEGs expression with high expression in red and low expression in blue. (B) The volcano plot of DEGs expression showed upregulated genes in red and downregulated genes in green.
3.2 Gene Ontology and gene set enrichment analyses
GO enrichment analysis includes 3 diverse categories: BP, CC, and MF. The top terms were selected according to the q-value from small to large, which could be observed in Figure 2A. With a total of 4 items in BP, DEGs were related to potassium ion transport, the regulation of membrane potential, and the regulation of ion transmembrane transport. Enriched to 6 items in CC, which mainly included the voltage-gated potassium channel complex, potassium channel complex, cation channel complex, ion channel complex, and transporter complex, and so on. Enriching 13 items in MF, it mainly contained potassium channel activity, potassium ion transmembrane transporter activity, voltage-gated channel activity, carbohydrate binding, channel activity, etc.
Figure 2 Gene Ontology (GO) enrichment analysis of DEGs and the gene set enrichment analysis (GSEA) of whole gene expressions. (A) Bubble plot of GO enrichment analysis, including BP, CC, and MF. (B, C) Line graph of GSEA in control and intervertebral disk degeneration (IVDD) samples.
As for GSEA, whole gene expressions were all included in analysis, the results was shown in Figures 2B, C; Table 2, which suggested that the cell cycle, DNA replication, graft versus host disease, nucleotide excision repair, and so on were enriched in control samples, while complement and coagulation cascades, Fc γ R–mediated phagocytosis, neuroactive ligand–receptor interaction, the NOD-like receptor (NLR) signaling pathway, gap junctions, etc., were enriched in IVDD samples.
3.3 Identification and validation of characteristic genes
The LASSO logistic regression and SVM-RFE algorithms were used to identify characteristic genes from the 5 DEGs. Firstly, 4 characteristic genes, namely, ZNF542P, KCNH4, ASMTL-AS1, and TBC1D31, were obtained via the LASSO logistic regression (Figure 3A). Moreover, 5 genes, namely, ASMTL-AS1, LGALSL, ZNF542P, TBC1D31 and KCNH4, were identified by the SVM-RFE algorithm (Figure 3B). Finally, ZNF542P, KCNH4, ASMTL-AS1, and TBC1D31 were screened as characteristic genes through overlapping the genes acquired by the two algorithms (Figure 3C). In order to further test the expression levels of these 4 characteristic genes, we performed the differential expression analysis in the validation set and found that the expression of ZNF542P was downregulated (p = 0.00017) in IVDD samples compared to control samples (Figure 3D). The expression of KCNH4, ASMTL-AS1, and TBC1D31 was not statistically significant (p = 0.92, 0.23, 0.33). Thus, ZNF542P was considered as a key characteristic gene.
Figure 3 Identification and validation of characteristic genes for IVDD. (A) The characteristic genes were obtained via the LASSO logistic regression. (B) The characteristic genes were identified by the SVM-RFE algorithm. (C) Venn diagram showed the overlapping genes acquired by the LASSO logistic regression and the SVM-RFE algorithm. (D) the expression of ZNF542P in the validation set.
3.4 Estimation of diagnostic value of key characteristic gene
Based on the key characteristic gene identified in this study, we performed ROC curve analysis to estimate the clinical diagnostic value. As shown in Figure 4, the results showed that the AUC of ZNF542P was 0.983 (95%CI: 0.926–1.000) in the training set (Figure 4A), while AUC was 0.879 (95%CI: 0.754–0.972) in the validation set (Figure 4B), respectively. Both of the AUCs were greater than 0.7, which indicated that ZNF542P could be used as the diagnostic biomarker.
Figure 4 Estimation of the diagnostic value of the key characteristic gene. (A) The ROC curve analysis of ZNF542P in the training set. (B) The ROC curve analysis of ZNF542P in the validation set.
3.5 Nucleus pulposus cell culture and staining identification
The primary NP cells began to adhere to the wall after 24 h (Figure 5A). The cells basically adhered to the wall 7 days later (Figure 5B). The cells were spindle shaped, with large and round nuclei and abundant cytoplasm. Compared with normal NP cells, the primary adherent time of degenerated NP cells was slightly longer. The cell morphology was poor and easy-to-appear aging cells after passage. Safranine O staining showed that the nucleus was red and cytoplasm staining was slight, indicating that cells could synthesize proteoglycans (Figure 5C). Toluidine blue staining showed that the nucleus was blue and cytoplasm staining was slight, demonstrating that cells could secrete glycosaminoglycans (Figure 5D).
Figure 5 Cell culture and staining identification. (A) Cells 7 days after inoculation. (B) Cells 21 days after inoculation. (C) Safranine O staining of NP cells (×200). (D) Toluidine blue staining of NP cells (×200).
3.6 The expression of the ZNF542P gene in nucleus pulposus cells
Compared with normal NP cells, the expression of the ZNF542P gene was decreased in degenerated NP cells (p < 0.05) shown in Figure 6.
Figure 6 The expression of the ZNF542P gene. Data presented as mean ± SD, *P < 0.05, vs. control group; as evaluated by the t-test.
3.7 The expression of NLRP3 and pro Caspase-1 in nucleus pulposus cells
Compared with normal NP cells, the expression of NLRP3 and pro Caspase-1 were increased in degenerated NP cells shown in Figure 7A, B.
Figure 7 The expression of NLRP3 and pro Caspase-1. (A) The protein band. (B) The protein quantification. Data presented as mean ± SD, *P < 0.05, vs. control group, as evaluated by the t-test. CON, control;DNP, degenerated NP.
3.8 Correlation between key characteristic genes and immune cell infiltration
The training set, including 11 IVDD samples and 11 control samples, was selected for the analysis of immune cell infiltration using the CIBERSORT algorithm, which was applied to calculate the relative percentage of 22 types of immune cells in each sample (Figure 8A). Compared with control samples, the results suggested that T cells gamma delta (γδT cells) were abnormally downregulated while neutrophils were upregulated in IVDD samples (Figure 8B). Moreover, the correlation analysis between the key characteristic gene and various infiltrating immune cells showed that the expression of ZNF542P was positively related to the proportions of γδT cells (R = 0.54, P = 0.0094), as shown in Figures 8C, D.
Figure 8 Correlation between key characteristic gene and immune cell infiltration. (A) Relative percentage of 22 immune cells in each samples. (B) Violin plot of the difference in each type of immune cell between IVDD and control samples. (C, D) The correlation analysis between key characteristic gene and various infiltrating immune cells.
4 Discussion
IVDD is a common spine disease associated with many risk factors, such as genetic susceptibility (7), mechanical stress (8), trauma (9), heavy workload (10), smoking (11), obesity (12), and aging (13) that may result in the occurrence of IVDD. Previous studies have suggested that a lot of pathological processes including ECM degradation (14), systemic inflammation (15), oxidative stress (16), mitochondrial dysfunction (17), telomere shortening (18) and DNA damage (19), nutritional deprivation (20), abnormal mechanical load (21), and epigenetic changes (22) may lead to IVDD. Currently, the early diagnosis of IVDD is restricted. To identify key characteristic gene for IVDD and explore the correlations of key characteristic gene and immune cell infiltration, we acquired 3 IVDD-related gene expression profiles from the GEO database and screened 5 DEGs including 3 upregulated genes (KCNH4, LGALSL, and ASMTL-AS1) and 2 downregulated genes (ZNF542P and TBC1D31).
Subsequently, GO enrichment analysis was conducted on the DEGs. As a result, these DEGs were enriched in BP, CC, and MF. BP mainly included potassium ion transport, the regulation of membrane potential, and the regulation of ion transmembrane transport. CC mainly contained the voltage-gated potassium channel complex, potassium channel complex, cation channel complex, ion channel complex, and transporter complex. Simultaneously, MF mainly included potassium channel activity, potassium ion transmembrane transporter activity, voltage-gated channel activity, carbohydrate binding, and channel activity.
In addition, we found that 42 pathways were related to whole gene expressions. As revealed by the GSEA, those genes were mostly associated with complement and coagulation cascades, Fc γ R–mediated phagocytosis, neuroactive ligand–receptor interaction, the NOD-like receptor signaling pathway, gap junctions, and so on. A study showed that there are a total of 1,785 proteins in normal lumbar disks (NDs) and 1,775 proteins in degenerate disks (DDs). The ND primarily participated in ECM maintenance and basic metabolic pathways, whereas the unique proteins of the DD group were involved in inflammatory pathways (complement and coagulation cascades)and infective channels that support the recent theories of inflammaging and subclinical infection as pathogens of IVDD (23). Fc γ R IIIa (FcγR3A)–mediated phagocytosis, the regulation of Toll-like receptors (TLRs) by endogenous ligand, neutrophil, and platelet degranulation, was a significantly changed pathway in patients with modic changes (24). The neuroactive ligand–receptor interaction pathway was a significant pathway response to osmotic stimuli in the intervertebral disk (25). The NOD-like receptor (NLR) signaling pathway is closely related to immune response. A study indicated that thioredoxin-interacting protein (TXNIP) could aggravate IVDD through the NLRP3/Caspase-1/IL-1β signaling pathway (26). The activation of the TLR4/NLRP3 signaling axis could cause apoptosis, oxidative stress, and inflammatory responses in NP cells. It has been reported that the activation of the NF-κB/NLRP3 signaling pathway could decrease cellular functions of NP cells and promoted cell pyroptosis (27). NF-κB signaling activation improved NLRP3 inflammasome activity, and NLRP3 activation increased NF-κB signaling activation induced by TNF-α. NLRP3 inhibition or NLRP3 knockdown induced autophagic flux in the presence of TNF-α (28). In our study, the results of Western blot suggested that the NOD-like receptor signaling pathway was activated. Another study illustrated that high mechanical loads could affect the intercellular communication via gap junctions in healthy annulus fibrosus tissue, while gap junctions additionally control chondrocyte physiology (29).
Next, the overlapping genes, ZNF542P, KCNH4, ASMTL-AS1, and TBC1D31, were identified as characteristic genes via LASSO logistic regression and SVM-RFE algorithms. Moreover, we tested the expression levels of these 4 characteristic genes, and found that the expression of ZNF542P was downregulated in IVDD samples compared to control samples in the validation set, whereas the expression of KCNH4, ASMTL-AS1, and TBC1D31, there were no changes. Thus, ZNF542P was considered as a key characteristic gene and had a good clinical diagnostic value via ROC curve analysis.
So far, there are few studies on the function of ZNF542P. It may be involved in transcriptional regulation. As a marker of statin response in vitro, its change is mostly correlated with statin-induced alternation in cellular cholesterol ester. Intracellular cholesterol ester levels upon simvastatin treatment were raised after ZNF542P knocking-down in a human hepatoma cell line (30). In our study, the results of qRT-PCR showed that compared with normal NP cells, the expression of ZNF542P gene was decreased in degenerated NP cells.
In view of the importance of immune response in IVDD, we assessed the proportion and composition of infiltrating immune cells in IVDD samples and control samples and found that neutrophils and γδT cells were significantly different. In particular, neutrophils were significantly increased while γδT cells was significantly decreased in IVDD samples compared to control samples. Neutrophils, a type of polymorphonuclear leukocyte existing in the human immune system, are well known as a major participant during acute inflammation. However, recent evidence has shown that they also participate in chronic inflammatory and adaptive immune responses (31). An epidural inflammation ranging from acute invasion of neutrophils to the formation of chronic granulation tissue existed in most dogs with IVDD (32).
γδT cells, as a part of the innate immune system, are a subset of T cells that play a vital role in various diseases, such as COVID-19 virus infections (33), cancer (34), hematologic malignancies (35), autoimmune diseases (36), and central nervous system diseases (37). γδT cells could respond to many molecular signals and have the ability to induce multiple cytokines, for example, GM-CSF, IL-4, IL-17, IL-21, IL-22, and IFN-γ (38). Recent data suggested the γδT cells were also involved in IVDD. Immune infiltration analysis showed that the imbalance of neutrophils and γδT cells were significantly related to IVDD (39). Finally, we found that the expression of ZNF542P was positively correlated with the proportions of γδT cells.
5 Conclusion
In conclusion, 5 DEGs were screened between IVDD samples and control samples. Moreover, ZNF542P was identified as a key characteristic gene in IVDD based on bioinformatics and machine learning algorithms, which were both downregulated and tested in the validation set. Furthermore, we found that the expression of ZNF542P may be associated with the NOD-like receptor signaling pathway and the infiltration of γδT cells. Therefore, our study may contribute to the understanding of IVDD and may help in improving the early diagnosis of IVDD. However, further studies are needed to study the roles of ZNF542P in IVDD.
Data availability statement
The original contributions presented in the study are included in the article/supplementary materials. Further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving human participants were reviewed and approved by Tianjin Hospital Ethics Committee. The patients/participants provided their written informed consent to participate in this study.
Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
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.
References
1. Traeger AC, Qaseem A, McAuley JH. Low back pain. JAMA (2021) 326:286. doi: 10.1001/jama.2020.19715
2. Li B, Yang Y, Wang L, Liu G. Stem cell therapy and exercise for treatment of intervertebral disc degeneration. Stem Cells Int (2021) 2021:7982333. doi: 10.1155/2021/7982333
3. Xin J, Wang Y, Zheng Z, Wang S, Na S, Zhang S. Treatment of intervertebral disc degeneration. Orthop Surg (2022) 14:1271-80. doi: 10.1111/os.13254
4. Wang Y, Zhao Y, Bollas A, Wang Y, Au KF. Nanopore sequencing technology, bioinformatics and applications. Nat Biotechnol (2021) 39:1348–65. doi: 10.1038/s41587-021-01108-x
5. Greener JG, Kandathil SM, Moffat L, Jones DTA. Guide to machine learning for biologists. Nat Rev Mol Cell Biol (2022) 23:40–55. doi: 10.1038/s41580-021-00407-0
6. Le T, Aronow RA, Kirshtein A, Shahriyari LA. Review of digital cytometry methods: estimating the relative abundance of cell types in a bulk of cells. Brief Bioinform (2021) 22:bbaa219. doi: 10.1093/bib/bbaa219
7. Wu F, Huang X, Zhang Z, Shao ZA. Meta-analysis assessing the association between COL11A1 and GDF5 genetic variants and intervertebral disc degeneration susceptibility. Spine (Phila Pa 1976) (2020) 45:E616–23. doi: 10.1097/BRS.0000000000003371
8. Zhao R, Liu W, Xia T, Yang L. Disordered mechanical stress and tissue engineering therapies in intervertebral disc degeneration. Polymers (Basel) (2019) 11:1151. doi: 10.3390/polym11071151
9. Zhang H-J, Liao H-Y, Bai D-Y, Wang Z-Q, Xie X-W. MAPK /ERK signaling pathway: a potential target for the treatment of intervertebral disc degeneration. Biomed Pharmacother (2021) 143:112170. doi: 10.1016/j.biopha.2021.112170
10. Han C, Kuang M-J, Ma J-X, Ma X-L. Prevalence of modic changes in the lumbar vertebrae and their associations with workload, smoking and weight in northern China. Sci Rep (2017) 7:46341. doi: 10.1038/srep46341
11. Kiraz M, Demir E. Relationship of lumbar disc degeneration with hemoglobin value and smoking. Neurochirurgie. (2020) 66:373–7. doi: 10.1016/j.neuchi.2020.06.133
12. Cannata F, Vadalà G, Ambrosio L, Fallucca S, Napoli N, Papalia R, et al. Intervertebral disc degeneration: a focus on obesity and type 2 diabetes. Diabetes. Metab Res Rev (2020) 36:e3224. doi: 10.1002/dmrr.3224
13. Francisco V, Pino J, González-Gay M.Á., Lago F, Karppinen J, Tervonen O, et al. A new immunometabolic perspective of intervertebral disc degeneration. Nat Rev Rheumatol (2022) 18:47–60. doi: 10.1038/s41584-021-00713-z
14. Lin Z, Ni L, Teng C, Zhang Z, Wu L, Jin Y, et al. Eicosapentaenoic acid-induced autophagy attenuates intervertebral disc degeneration by suppressing endoplasmic reticulum stress, extracellular matrix degradation, and apoptosis. Front Cell Dev Biol (2021) 9:745621. doi: 10.3389/fcell.2021.745621
15. Sharma A. The role of adipokines in intervertebral disc degeneration. Med Sci (Basel) (2018) 6:34. doi: 10.3390/medsci6020034
16. Lu S, Song Y, Luo R, Li S, Li G, Wang K, et al. Ferroportin-dependent iron homeostasis protects against oxidative stress-induced nucleus pulposus cell ferroptosis and ameliorates intervertebral disc degeneration. In Vivo Oxid Med Cell Longev (2021) 2021:6670497. doi: 10.1155/2021/6670497
17. Saberi M, Zhang X, Mobasheri A. Targeting mitochondrial dysfunction with small molecules in intervertebral disc aging and degeneration. Geroscience (2021) 43:517–37. doi: 10.1007/s11357-021-00341-1
18. Chen H-W, Zhang G-Z, Liu M-Q, Zhang L-J, Kang J-H, Wang Z-H, et al. Natural products of pharmacology and mechanisms in nucleus pulposus cells and intervertebral disc degeneration. Evid Based Complement Alternat Med (2021) 2021:9963677. doi: 10.1155/2021/9963677
19. Pitcher LE, Wong AK, Robbins PD. Targeting ATM to mitigate intervertebral disc degeneration. Aging (Albany NY) (2021) 13:10814–5. doi: 10.18632/aging.203015
20. Yurube T, Ito M, Kakiuchi Y, Kuroda R, Kakutani K. Autophagy and mTOR signaling during intervertebral disc aging and degeneration. JOR Spine (2020) 3:e1082. doi: 10.1002/jsp2.1082
21. Fu F, Bao R, Yao S, Zhou C, Luo H, Zhang Z, et al. Aberrant spinal mechanical loading stress triggers intervertebral disc degeneration by inducing pyroptosis and nerve ingrowth. Sci Rep (2021) 11:772. doi: 10.1038/s41598-020-80756-6
22. Kawarai Y, Jang SH, Lee S, Millecamps M, Kang H, Gregoire S, et al. Exercise attenuates low back pain and alters epigenetic regulation in intervertebral discs in a mouse model. Spine J (2021) 21:1938–49. doi: 10.1016/j.spinee.2021.06.002
23. Rajasekaran S, Soundararajan DCR, Tangavel C, Nayagam SM, K S SV, R S, et al. Uncovering molecular targets for regenerative therapy in degenerative disc disease: do small leucine-rich proteoglycans hold the key? Spine J (2021) 21:5–19. doi: 10.1016/j.spinee.2020.04.011
24. Rajasekaran S, Soundararajan DCR, Nayagam SM, Tangavel C, Raveendran M, Thippeswamy PB, et al. Modic changes are associated with activation of intense inflammatory and host defense response pathways - molecular insights from proteomic analysis of human intervertebral discs. Spine J (2022) 22:19–38. doi: 10.1016/j.spinee.2021.07.003
25. Wang L, Liu C, Tian JW, Li GQ, Zhao QH. Significant pathways response to osmotic stimuli in the intervertebral disc. Eur Rev Med Pharmacol Sci (2013) 17:1620–5.
26. Zhou Y, Chen Z, Yang X, Cao X, Liang Z, Ma H, et al. Morin attenuates pyroptosis of nucleus pulposus cells and ameliorates intervertebral disc degeneration via inhibition of the TXNIP/NLRP3/Caspase-1/IL-1β signaling pathway. Biochem Biophys Res Commun (2021) 559:106–12. doi: 10.1016/j.bbrc.2021.04.090
27. Wang D, Cai X, Xu F, Kang H, Li Y, Feng R. Ganoderic acid a alleviates the degeneration of intervertebral disc via suppressing the activation of TLR4/NLRP3 signaling pathway. Bioengineered (2022) 13:11684–93. doi: 10.1080/21655979.2022.2070996
28. Chen J, Yan J, Li S, Zhu J, Zhou J, Li J, et al. Atorvastatin inhibited TNF-α induced matrix degradation in rat nucleus pulposus cells by suppressing NLRP3 inflammasome activity and inducing autophagy through NF-κB signaling. Cell Cycle (2021) 20:2160–73. doi: 10.1080/15384101.2021.1973707
29. Desrochers J, Duncan NA. Intercellular communication via gap junctions affected by mechanical load in the bovine annulus fibrosus. Comput Methods Biomech Biomed Engin (2014) 17:64–71. doi: 10.1080/10255842.2012.717268
30. Kim K, Theusch E, Kuang Y-L, Dose A, Mitchel K, Cubitt C, et al. ZNF542P is a pseudogene associated with LDL response to simvastatin treatment. Sci Rep (2018) 8:12443. doi: 10.1038/s41598-018-30859-y
31. Kolaczkowska E, Kubes P. Neutrophil recruitment and function in health and inflammation. Nat Rev Immunol (2013) 13:159–75. doi: 10.1038/nri3399
32. Fadda A, Oevermann A, Vandevelde M, Doherr MG, Forterre F, Henke D. Clinical and pathological analysis of epidural inflammation in intervertebral disk extrusion in dogs. J Vet Intern Med (2013) 27:924–34. doi: 10.1111/jvim.12095
33. Massow G. Von; oh, s.; lam, a.; gustafsson, k. gamma delta T cells and their involvement in COVID-19 virus infections. Front Immunol (2021) 12:741218. doi: 10.3389/fimmu.2021.741218
34. Deng J, Yin H. Gamma delta (γδ) T cells in cancer immunotherapy; where it comes from, where it will go? Eur J Pharmacol (2022) 919:174803. doi: 10.1016/j.ejphar.2022.174803
35. Goyal A, Nardi V. Gamma/Delta (γδ) T cells: the role of the T-cell receptor in diagnosis and prognosis of hematologic malignancies. Am J Dermatopathol (2022) 44:237–48. doi: 10.1097/DAD.0000000000002041
36. Bank I. The role of gamma delta T cells in autoimmune rheumatic diseases. Cells (2020) 9:462. doi: 10.3390/cells9020462
37. Park JH, Kang In, Lee HK. γδ T cells in brain homeostasis and diseases. Front Immunol (2022) 13:886397. doi: 10.3389/fimmu.2022.886397
38. Wo J, Zhang F, Li Z, Sun C, Zhang W, Sun G. The role of gamma-delta T cells in diseases of the central nervous system. Front Immunol (2020) 11:580304. doi: 10.3389/fimmu.2020.580304
Keywords: intervertebral disc degeneration, biomarker, nucleus pulposus cell, immune cell infiltration, bioinformatics, machine learning algorithms
Citation: Yang Z, Yuan Z-Z and Ma X-L (2023) Identification of a potential novel biomarker in intervertebral disk degeneration by bioinformatics analysis and experimental validation. Front. Immunol. 14:1136727. doi: 10.3389/fimmu.2023.1136727
Received: 04 January 2023; Accepted: 10 May 2023;
Published: 30 May 2023.
Edited by:
Wei Xu, Texas A&M University Corpus Christi, United StatesReviewed by:
He Zhang, Shengjing Hospital of China Medical University, ChinaJie Zhao, Tianjin University, China
Copyright © 2023 Yang, Yuan and Ma. 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: Zhao Yang, YW5ueS5hbGxhbkAxMjYuY29t; Xin-Long Ma, bWF4aW5sb25nODZAMTI2LmNvbQ==
†These authors have contributed equally to this work and share first authorship