Abstract
Background:
Intervertebral Disc Degeneration (IDD) is a major cause of lower back pain and a significant global health issue. However, the specific mechanisms of IDD remain unclear. This study aims to identify key genes and pathways associated with IDD using bioinformatics and machine learning algorithms.
Methods:
Gene expression profiles, including those from 35 LDH patients and 43 healthy volunteers, were downloaded from the GEO database (GSE124272, GSE150408, GSE23130, GSE153761). After merging four microarray datasets, differentially expressed genes (DEGs) were selected for GO and KEGG pathway enrichment analysis. Weighted Gene Co-expression Network Analysis (WGCNA) was then applied to the merged dataset to identify relevant modules and intersect with DEGs to discover candidate genes with diagnostic value. A LASSO model was established to select appropriate genes, and ROC curves were drawn to elucidate the diagnostic value of genetic markers. A Protein-Protein Interaction (PPI) network was constructed and visualized to determine central genes, followed by external validation using qRT-PCR.
Results:
Differential analysis of the preprocessed dataset identified 244 genes, including 183 upregulated and 61 downregulated genes. WGCNA analysis revealed the most relevant module intersecting with DEGs, yielding 9 candidate genes. The lasso-cox method was used for regression analysis, ultimately identifying 6 genes: ASPH, CDC42EP3, FOSL2, IL1R1, NFKBIZ, TCF7L2. A Protein-Protein Interaction (PPI) network created with GENEMANIA identified IL1R1 and TCF7L2 as central genes.
Conclusion:
Our study shows that IL1R1 and TCF7L2 are the core genes of IDD, offering new insights into the pathogenesis and therapeutic development of IDD.
1 Introduction
Lower back pain is one of the most common medical issues worldwide, affecting 70% to 85% of the population globally (). Intervertebral Disc Degeneration (IDD) is the primary cause of lower back pain (). Currently, the treatment strategies for IDD mainly include conservative and surgical treatments, but these methods often only alleviate pain without halting the progression of IDD (). Surgical treatment can have long-term adverse effects on biomechanics, while conservative treatment is usually temporary (). IDD is a complex, multifactorial process involving changes in phenotype and genotype, hence, gene-targeted therapy may be a promising direction for future IDD treatments (). Recent studies have found that immune cell infiltration and the increase in inflammatory mediators and oxidative stress play a crucial role in the development of IDD, accelerating the inflammatory cascade and ultimately leading to discogenic pain (). Therefore, a deeper understanding of the pathological changes and mechanisms of IDD is essential for diagnosis and the development of new treatment methods.
Bioinformatics, an emerging biotechnology that has grown with the Human Genome Project, utilizes computer technology to store, retrieve, analyze, and interpret large volumes of biological data (). In the last two decades, with the advent of the big data era and the continuous growth and inherent complexity of biological data, machine learning, as an interdisciplinary field, has been widely applied to construct informational models and predictions of biological processes, significantly impacting the field of biology.
The aim of this study is to identify core genes associated with IDD by analyzing RNA expression profiles from the Gene Expression Omnibus (GEO) database and to validate these genes through in vitro experiments, with the goal of establishing these core genes as potential biomarkers for IDD.
2 Methods
2.1 Collection and preprocessing of microarray data
This study downloaded four gene expression profile datasets related to Intervertebral Disc Degeneration (IDD) from the Gene Expression Omnibus (GEO) database. The datasets include: GSE124272 (GPL21185), GSE153761 (GPL22120), GSE150408 (GPL21185), and GSE23130 (GPL1352), covering gene expression data from patients with lumbar disc herniation, cervical spondylotic myelopathy, and lumbar disc degeneration, as well as healthy volunteers. Data grouping was based on the Pfirrmann grading system (), with grades I-III as the control group and grades IV-V as the experimental group. Initially, the R package inSilicoMerging was used to merge these datasets, followed by the application of a method proposed by Johnson et al (). to remove batch effects, resulting in a corrected data matrix containing samples from 35 experimental and 43 control group participants.
2.2 Selection of differentially expressed genes
After merging the four datasets, further correction and normalization of the expression matrix were performed. Differential analysis was conducted using the Limma R package, with |log2FC| > 1.3 and p < 0.05 set as the criteria for selecting differentially expressed genes (DEGs). Subsequently, heatmaps and volcano plots were utilized to visually display the expression patterns of DEGs.
2.3 GO and KEGG analysis
Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were conducted to systematically analyze gene functions and key biological pathways closely related to Intervertebral Disc Degeneration (). GO analysis results were categorized into Biological Processes (BP), Cellular Components (CC), and Molecular Functions (MF) (). Enrichment analysis was performed using the DAVID database to identify statistically significant GO terms and KEGG pathways, with a selection threshold set at p < 0.05. All results were visualized using the ggplot2 R package.
2.4 WGCNA analysis
Weighted Gene Co-expression Network Analysis (WGCNA) was performed using the WGCNA package in R. Initially, outlier data were filtered to enhance model stability, and an appropriate soft threshold β was selected to construct the Topological Overlap Matrix (TOM). Hierarchical clustering methods were then used to generate a dendrogram of genes, and gene significance (GS) and module membership (MM) were calculated to assess the correlation between genes and clinical information, as well as to analyze the significant association between modules and disease characteristics.
2.5 Identification of characteristic genes and application of machine learning algorithms
Characteristic module genes were identified through cross-analysis with differentially expressed genes (DEGs) from the WGCNA analysis. Subsequently, the glmnet package in R was used to implement the LASSO machine learning algorithm, selecting significant genes with the highest predictive value for IDD. LASSO regression aims to simplify high-dimensional data by selecting characteristic genes with non-zero coefficients, followed by multivariate logistic regression analysis to identify all significant risk DEGs.
2.6 Construction of protein-protein interaction networks
The GeneMANIA network platform was utilized to construct Protein-Protein Interaction (PPI) networks to explore interactions and functions among core genes. GeneMANIA aids in the analysis and understanding of PPI networks by generating gene function predictions and identifying genes with similar functions ().
2.7 Clinical diagnostic value assessment of core genes
Receiver Operating Characteristic (ROC) curves, drawn using the “pROC” package in R, were employed to evaluate the efficacy of core genes in the diagnosis of IDD. ROC curves depict the relationship between sensitivity and specificity, estimating the clinical diagnostic value of core genes by calculating the area under the curve (AUC).
2.8 Sample collection
This study protocol was approved by the Ethics Committee of Zhongshan Hospital affiliated with Dalian University. All participating patients signed informed consent forms. The severity of IDD in patients was assessed using magnetic resonance imaging (MRI) according to the Pfirrmann grading system (). Intervertebral disc tissue samples were obtained from patients undergoing surgery for lumbar disc herniation.Finally, two patients with Pfirrmann grades III and V were selected. After the interdisc tissue was obtained, it was washed with phosphate buffer (PBS) to remove the impurities completely. The nucleus pulposum tissue was cut into 1mm fragments, digested with 0.2% type II collagenase, and incubated in a 37° incubator for digestion and incubation for 2 hours. The digestion was stopped with fetal bovine serum and phosphate buffer (PBS) containing more than 20% by volume, and then centrifuged. The superserum was discarded, transferred to dmem medium containing 10% fetal bovine serum and 1% penicillin streptomycin, and cultured in a co2 incubator. After 3 days of inoculation, the cell morphology of nucleus pulposum cells was observed with an inverted microscope. If the fusion degree reached 80%-90%, the subculture was performed again. Repeat this step to pass to the third generation.
2.9 Real-time quantitative polymerase chain reaction
Total RNA was extracted from nucleus pulposum tissues using StradyPure Universal RNA extraction kit, then transcribed into cDNA using the Cofitt Reverse Transcription Kit. The expression levels of target genes were determined using a high-specificity qRT-PCR detection kit, with specific PCR primer sequences including IL1R1 and TCF7L2.
2.10 Statistical analysis
Data processing and chart generation were performed using R 4.2.1 and GraphPad Prism 9 software. Statistical evaluation of the data was conducted using one-way ANOVA in SPSS 25.0 software, with differences considered statistically significant at P < 0.05.
3 Results
3.1 Preprocessing and selection of differentially expressed genes
After merging all datasets and removing batch effects, the data were normalized. The results showed that, after processing, the average expression levels of genes across the datasets were consistent (see Figures 1A, B). Differential analysis revealed a total of 244 differentially expressed genes (DEGs), including 183 upregulated and 61 downregulated genes. These results were visually displayed using a volcano plot (see Figure 1C). Additionally, the expression patterns of the identified DEGs were further presented using a heatmap (see Figure 1D).
Figure 1
3.2 Gene ontology enrichment analysis
The Gene Ontology (GO) enrichment analysis revealed three main categories: Biological Processes (BP), Cellular Components (CC), and Molecular Functions (MF), selecting GO terms with a p-value of less than 0.05 and ranking them based on the number of involved genes (see Figures 2A–C). In the category of Biological Processes, a total of 10 GO terms were enriched, mainly involving immune system processes, developmental processes, response to stimulus, and response to stress. In the category of Cellular Components, 10 main GO terms were identified, including endomembrane systems, vesicles, extracellular regions, and intrinsic components of membranes. In the category of Molecular Functions, 6 items were enriched, covering glycosaminoglycan binding, carbohydrate binding, protease binding, extracellular matrix structural constituents, hyaluronic acid binding, and cell-cell adhesion mediator activity. KEGG enrichment analysis showed that differentially expressed genes (DEGs) are highly related to complement and coagulation cascades, phagosome pathways, and osteoclast differentiation pathways, with specific results shown in Figure 2D).
Figure 2
3.3 WGCNA analysis
In this study, Weighted Gene Co-expression Network Analysis (WGCNA) was utilized to identify gene modules potentially related to Intervertebral Disc Degeneration (IDD). Initially, all candidate genes and samples were screened to confirm the absence of abnormally high proportions of missing values. Subsequently, through sample clustering analysis, evident outlier samples were excluded with a height cutoff value set at 140 (see Figure 3A). The “sft$powerEstimate” was used to determine the optimal soft threshold, selecting 10 as the soft threshold (scale independence reached 0.85) to better differentiate between IDD and normal tissues (see Figures 3B, C). With a soft threshold of 10, a minimum module size of 30, and a deepSplit parameter of 3, 17 gene co-expression modules were ultimately identified (see Figure 3D). By drawing the inter-module relationship map (see Figure 3E) and analyzing the module-IDD feature association using Spearman correlation coefficients, the darkgrey module was found to have the most significant correlation with IDD (see Figure 3F). Additionally, a significant correlation was observed between the module membership (MM) and gene significance (GS) within the darkgrey module (see Figure 3G). Within the darkgrey module, 56 genes important for module function were identified. Through Venn diagram analysis, co-expressed genes between the differentially expressed genes (DEGs) and the WGCNA-derived feature module genes were identified (see Figure 3H).
Figure 3
3.4 Identification of characteristic genes and application of machine learning algorithms
Using the R package glmnet, combined with survival time, status, and gene expression data, characteristic genes were identified through lasso-cox regression analysis. A 5-fold cross-validation was utilized to determine the optimal model, selecting a Lambda value of 0.0462682649698763, ultimately identifying 6 genes: ASPH, CDC42EP3, FOSL2, IL1R1, NFKBIZ, TCF7L2 (see Figures 4A, B). To visually present the expression levels of these characteristic genes, box plots were drawn, showing an upward trend for these 6 genes in IDD samples compared to control groups (see Figure 4C).
Figure 4
3.5 Construction of PPI networks to identify core genes
Using the GeneMANIA database, we analyzed the Protein-Protein Interaction (PPI) network composed of six genes, ASPH, CDC42EP3, FOSL2, IL1R1, NFKBIZ, TCF7L2, and their interacting genes (see Figure 5). In this network, core genes are represented by circles at the center. Around these core genes are 20 circles, representing genes closely related to them, determined based on physical interactions, predictions, biological pathways, and genetic interactions. Through a comprehensive analysis of the PPI network, IL1R1 and TCF7L2 were identified as the final core genes of this study.
Figure 5
3.6 Evaluation of the clinical diagnostic value of core genes
Receiver Operating Characteristic (ROC) curve analysis was conducted on two core genes, IL1R1 and TCF7L2, to estimate their clinical diagnostic value. As shown in Figures 6A, B, the Area Under the Curve (AUC) values for both genes were close to 0.7, It is suggested that IL1R1 and TCF7L2 have certain clinical diagnostic value.
Figure 6
3.7 Expression validation of core genes
Using real-time quantitative polymerase chain reaction (qRT-PCR) analysis, we validated the expression levels of two core genes, IL1R1 and TCF7L2, in the Intervertebral Disc Degeneration (IDD) group compared to the normal intervertebral disc group. The results indicated that the expression levels of IL1R1 and TCF7L2 were significantly higher in the IDD group than in the normal group (see Figure 7).
Figure 7
4 Discussion
Intervertebral Disc Degeneration (IDD) is a significant cause of lower back pain, imposing a substantial burden on society and the economy. Globally, approximately 700 million patients are affected, with lower back pain not only impacting the physical and mental health of patients but also causing economic losses by reducing social productivity (). Therefore, studying the pathological changes and molecular mechanisms of IDD is crucial for improving clinical diagnoses and treatment methods. This study delves into the molecular mechanisms of IDD through bioinformatics analysis and machine learning algorithms.
Through differential analysis, we identified 183 upregulated and 61 downregulated genes from four expression profiles obtained from the GEO database, and these differentially expressed genes (DEGs) underwent GO and KEGG enrichment analyses. These genes were enriched in biological processes (BP), cellular components (CC), and molecular functions (MF), unveiling key biological pathways involved in the progression of IDD. WGCNA analysis further produced 17 co-expressed modules, among which the dark grey module showed the highest relevance to IDD, a finding confirmed by the correlation analysis between module membership (MM) and gene significance (GS). Through LASSO regression analysis and PPI network analysis, IL1R1 and TCF7L2 were determined as the core genes of this study. ROC curve analysis indicated that these two genes have good clinical diagnostic value.
Moreover, in vitro experiments through RT-qPCR validated the expression differences of IL1R1 and TCF7L2 in IDD patients, consistent with the bioinformatics analysis results.
IL1R1, also known as interleukin-1 receptor type 1, is a gene belonging to the cytokine receptor family of the interleukin-1 receptor family. The protein encoded is a receptor for interleukin 1α, interleukin 1β, and the interleukin 1 receptor antagonist, playing a pivotal role as a mediator in many cytokine-induced immune and inflammatory responses.The pain caused by intervertebral disc degeneration (IDD) is associated with the interaction among nucleus pulposus cells, annulus fibrosus cells, and immune cells such as macrophages, T cells, and neutrophils, leading to the release of pro-inflammatory cytokines (). The primary cytokines include TNF-α, IL-1, IL-6, IL-17, IL-2, IL-4, and IFN-γ (, ), among which TNF-α and IL-1 are the most extensively studied.Research has shown an increase in the gene and protein production of the interleukin-1 receptor type I (IL1R1), rather than the TNF receptor type I, in degenerative intervertebral discs (IVD). Thus, although both cytokines may be involved in the pathogenesis of IVD degeneration, IL-1 is likely to play a more significant role than TNF-α (). The IL-1 cytokines, IL-1α and IL-1β, act as soluble factors that interact with IL-1R1 and, together with IL-1R3, form a trimeric signaling IL-1 receptor complex. The biological effects of IL-1 are primarily inflammatory, activating target cells at the onset of defensive inflammatory responses. These effects include the induction of inflammatory cytokines, the production of toxic reactive oxygen and nitrogen species, prostaglandins, proteolytic enzymes, etc. The expression of IL-1R1 can be positively and negatively regulated by inflammation-related factors and other stress signals (). Studies have shown that IL-1 is produced in degenerating intervertebral discs. In degenerative discs, an imbalance between IL-1 and its inhibitor IL-1Ra disrupts the maintenance of matrix homeostasis, leading to an increase in degradative enzymes and a decrease in matrix protein gene expression, characteristic of disc degeneration. IL-1 regulates inflammatory responses through IL-1R1, thus playing a significant role in the development of IDD. Therefore, further research into IL-1R1 as a potential biomarker for IDD is warranted ().
TCF7L2 is a member of the T-cell factor/lymphoid enhancer-binding factor family (TCF/LEF), a protein-coding gene known as one of the most potent risk loci for type 2 diabetes, involved in several fundamental processes including adipogenesis, autophagy, and alternative splicing (). Chen S et al. discovered that TCF7L2 is overexpressed in degenerative intervertebral disc tissue samples (). Silencing TCF7L2 improved lipopolysaccharide (LPS)-induced nucleus pulposus (NP) cell proliferation, extracellular matrix (ECM) synthesis, and alleviated LPS-induced NP cell senescence, marking the first demonstration of the miR-1260b/TCF7L2 axis’s role in maintaining the chondrocyte-like phenotype of NP cells under LPS stimulation and in the synthesis of the extracellular matrix by NP cells.Sun J et al.’s findings indicate that the inhibition of TCF7L2 regulated by miR-155, through the p65/NF-κB signaling pathway, suppresses matrix degradation (). Inhibiting TCF7L2 may hold therapeutic potential for intervertebral disc degeneration (IDD), as related research has shown a close association between TCF7L2 and the development of IDD.
Nonetheless, this study has limitations, such as the in vitro experiments only validating the mRNA levels of the core genes without confirming protein levels through western blot. Therefore, more in vivo and in vitro experiments are needed to elucidate the roles and potential mechanisms of these two core genes in IDD.
In summary, this study applied bioinformatics and machine learning methods to analyze IDD’s differential genes, including functional enrichment, WGCNA, LASSO regression, PPI network, and ROC curve analyses, and validated the mRNA expression of IL1R1 and TCF7L2 through RT-qPCR experiments, consistent with bioinformatics results. This provides new insights into the pathogenesis of IDD and aids in identifying therapeutic targets.
5 Conclusion
This study comprehensively analyzed the pathogenesis of Intervertebral Disc Degeneration (IDD) and identified potential target genes, IL1R1 and TCF7L2, closely related to the progression of IDD. These findings not only enhance our understanding of the mechanisms underlying IDD but also provide potential therapeutic targets for the early diagnosis and treatment of IDD, warranting further in-depth research.
Statements
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 humans were approved by Ethics Committee of Affiliated Zhongshan Hospital of Dalian University. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
XY: Writing – original draft, Writing – review & editing. HZ: Writing – original draft, Writing – review & editing. SS: Writing – original draft, Writing – review & editing. XH: Writing – original draft. CG: Writing – original draft. ZJZ: Writing – original draft. ZTZ: Writing – review & editing. JG: Writing – review & editing. MZ: Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. The authors disclose receipt of the following financial or material support for the research, authorship, and/or publication of this article: This study was supported by Dalian Science and Technology Innovation Fund (2021JJ13SN68).
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
KnezevicNNCandidoKDVlaeyenJWSVan ZundertJCohenSP. Low back pain. Lancet. (2021) 398:78–92. doi: 10.1016/S0140-6736(21)00733-9
2
SmoldersLAMeijBPOnisDRiemersFMBergknutNWubboltsRet al. Gene expression profiling of early intervertebral disc degeneration reveals a down-regulation of canonical Wnt signaling and caveolin-1 expression: implications for development of regenerative strategies. Arthritis Res Ther. (2013) 15:R23. doi: 10.1186/ar4157
3
StichSJagielskiMFleischmannAMeierCBussmannPKohlBet al. Degeneration of lumbar intervertebral discs: characterization of anulus fibrosus tissue and cells of different degeneration grades. Int J Mol Sci. (2020) 21(6). doi: 10.3390/ijms21062165
4
HermantinFUPetersTQuartararoLKambinP. A prospective, randomized study comparing the results of open discectomy with those of video-assisted arthroscopic microdiscectomy. J Bone Joint Surg Am. (1999) 81:958–65. doi: 10.2106/00004623-199907000-00008
5
SamparaPBanalaRRVemuriSKAvGRGpvS. Understanding the molecular biology of intervertebral disc degeneration and potential gene therapy strategies for regeneration: a review. Gene Ther. (2018) 25:67–82. doi: 10.1038/s41434-018-0004-0
6
WangLHeTLiuJTaiJWangBZhangLet al. Revealing the immune infiltration landscape and identifying diagnostic biomarkers for lumbar disc herniation. Front Immunol. (2021) 12:666355. doi: 10.3389/fimmu.2021.666355
7
WangYZhaoYBollasAWangYAuKF. Nanopore sequencing technology, bioinformatics and applications. Nat Biotechnol. (2021) 39:1348–65. doi: 10.1038/s41587-021-01108-x
8
WangYXJ. Several concerns on grading lumbar disc degeneration on MR image with Pfirrmann criteria. J Orthop Translat. (2022) 32:101–2. doi: 10.1016/j.jot.2021.12.003
9
JohnsonWELiCRabinovicA. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. (2007) 8:118–27. doi: 10.1093/biostatistics/kxj037
10
KanehisaMGotoS. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. (2000) 28:27–30. doi: 10.1093/nar/28.1.27
11
AshburnerMBallCABlakeJABotsteinDButlerHCherryJMet al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. (2000) 25:25–9. doi: 10.1038/75556
12
FranzMRodriguezHLopesCZuberiKMontojoJBaderGDet al. GeneMANIA update 2018. Nucleic Acids Res. (2018) 46:W60–w4. doi: 10.1093/nar/gky311
13
PfirrmannCWMetzdorfAZanettiMHodlerJBoosN. Magnetic resonance classification of lumbar intervertebral disc degeneration. Spine. (2001) 26:1873–8. doi: 10.1097/00007632-200109010-00011
14
FranciscoVPinoJGonzález-GayMLagoFKarppinenJTervonenOet al. A new immunometabolic perspective of intervertebral disc degeneration. Nat Rev Rheumatol. (2022) 18:47–60. doi: 10.1038/s41584-021-00713-z
15
YamamotoJMaenoKTakadaTKakutaniKYurubeTZhangZet al. Fas ligand plays an important role for the production of pro-inflammatory cytokines in intervertebral disc nucleus pulposus cells. J Orthop Res. (2013) 31:608–15. doi: 10.1002/jor.22274
16
ShamjiMFSettonLAJarvisWSoSChenJJingLet al. Proinflammatory cytokine expression profile in degenerated and herniated human intervertebral disc tissues. Arthritis Rheum. (2010) 62:1974–82. doi: 10.1002/art.27444
17
SéguinCAPilliarRMRoughleyPJKandelRA. Tumor necrosis factor-alpha modulates matrix production and catabolism in nucleus pulposus tissue. Spine. (2005) 30:1940–8. doi: 10.1097/01.brs.0000176188.40263.f9
18
Le MaitreCLHoylandJAFreemontAJ. Catabolic cytokine expression in degenerate and herniated human intervertebral discs: IL-1beta and TNFalpha expression profile. Arthritis Res Ther. (2007) 9:R77. doi: 10.1186/ar2275
19
BoraschiDItalianiPWeilSMartinMU. The family of the interleukin-1 receptors. Immunol Rev. (2018) 281:197–232. doi: 10.1111/imr.12606
20
Le MaitreCLFreemontAJHoylandJA. The role of interleukin-1 in the pathogenesis of human intervertebral disc degeneration. Arthritis Res Ther. (2005) 7:R732–45. doi: 10.1186/ar1732
21
Del Bosque-PlataLMartínez-MartínezEEspinoza-CamachoMGragnoliC. The role of TCF7L2 in type 2 diabetes. Diabetes. (2021) 70:1220–8. doi: 10.2337/db20-0573
22
ChenSShiGZengJLiPHPengYDingZet al. MiR-1260b protects against LPS-induced degenerative changes in nucleus pulposus cells through targeting TCF7L2. Hum Cell. (2022) 35:779–91. doi: 10.1007/s13577-021-00655-4
23
SunJHongJSunSWangXPengYZhouJet al. Transcription factor 7-like 2 controls matrix degradation through nuclear factor κB signaling and is repressed by microRNA-155 in nucleus pulposus cells. BioMed Pharmacother. (2018) 108:646–55. doi: 10.1016/j.biopha.2018.09.076
Summary
Keywords
intervertebral disc degeneratio1, bioinformatics, machine learning, genes, IL1R15, TCF7L2
Citation
Zhang H, Shi S, Huang X, Gong C, Zhang Z, Zhao Z, Gao J, Zhang M and Yu X (2024) Identification of core genes in intervertebral disc degeneration using bioinformatics and machine learning algorithms. Front. Immunol. 15:1401957. doi: 10.3389/fimmu.2024.1401957
Received
16 March 2024
Accepted
27 June 2024
Published
10 July 2024
Volume
15 - 2024
Edited by
Bo Gao, Air Force Military Medical University, China
Reviewed by
Fuxin Wei, Sun Yat-sen University, China
Chen Xu, Shanghai Changzheng Hospital, China
Updates
Copyright
© 2024 Zhang, Shi, Huang, Gong, Zhang, Zhao, Gao, Zhang and Yu.
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: Xiaobing Yu, yuxiaobing1976@hotmail.com
†These authors have contributed equally to this work and share first authorship
Disclaimer
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.