- Research article
- Open Access
Identification of key ferroptosis-related biomarkers in steroid-induced osteonecrosis of the femoral head based on machine learning
Journal of Orthopaedic Surgery and Research volume 18, Article number: 327 (2023)
This study was aimed to identify key ferroptosis-related biomarkers in steroid-induced osteonecrosis of the femoral head (SONFH) based on machine learning algorithm.
The SONFH dataset GSE123568 (including 30 SONFH patients and 10 controls) was used in this study. The differentially expressed genes (DEGs) were selected between SONFH and control groups, which were subjected to WGCNA. Ferroptosis-related genes were downloaded from FerrDb V2, which were then compared with DEGs and module genes. Two machine learning algorithms were utilized to identify key ferroptosis-related genes, and the underlying mechanisms were analyzed by GSEA. Correlation analysis between key ferroptosis-related genes and immune cells was analyzed by Spearman method. The drug–gene relationships were predicted in CTD.
Total 2030 DEGs were obtained. WGCNA identified two key modules and obtained 1561 module genes. Finally, 43 intersection genes were identified as disease-related ferroptosis-related genes. After LASSO regression and RFE-SVM algorithms, 4 intersection genes (AKT1S1, BACH1, MGST1 and SETD1B) were considered as key ferroptosis-related gene. The 4 genes were correlated with osteoclast differentiation pathway. Twenty immune cells with significant differences were obtained between the groups, and the 4 key ferroptosis-related genes were correlated with most immune cells. In CTD, 41 drug–gene relationship pairs were finally obtained.
The 4 key ferroptosis-related genes, AKT1S1, BACH1, MGST1 and SETD1B, were identified to play a critical role in SONFH progression through osteoclast differentiation and immunologic mechanisms. Additionally, all the 4 genes had good disease prediction effect and could act as biomarkers for the diagnosis and treatment of SONFH.
Steroid-induced osteonecrosis of the femoral head (SONFH) is a devastating disease that usually progresses to osteoarthritis of the hip joint and femoral head collapse . The most common symptoms of SONFH are movement limitation and severe pain . The male gender and longer symptom duration are the risk factor for poor prognosis . The current treatment measures for osteonecrosis of the femoral head (ONFH) mainly included operative management  and conservative treatment, such as osteotomy , core decompression  and bone marrow-derived cell therapies . Given that the outcomes of ONFH were heterogeneous, there is no most approved therapy for ONFH patients. In China, there are around 1.5 × 105–2 × 105 new SONFH cases annually . Thus, it is imperative to discover the novel biomarkers for the diagnosis and treatment of SONFH. Although the etiology and pathogenesis of SONFH have been extensively studied , there is no clear consensus on its exact origin.
Ferroptosis is a recently discovered form of iron-mediated cell death, which causes much attention as new regulated necrosis . Ferroptosis is characterized by an increased level of lipid peroxidation products and reactive oxygen species [11, 12]. The morphological characteristics of ferroptosis include the obvious shrinkage of cell mitochondria, reduction or disappearance of mitochondrial crest, but the cell membrane is intact with normal nucleus size . Dysregulation of ferroptosis is correlated with a lot of pathological processes, such as inflammation-related diseases, neurodegenerative diseases, and cancers [14–16]. Recently, Sun et al.  reported that dexamethasone could induce ferroptosis through the pathway of P53/SLC7A11/GPX4 in glucocorticoid-induced osteonecrosis of the femoral head. Whereas, the study on the roles of ferroptosis in SONFH is still scarce.
Presently, this study was aimed to analyze key ferroptosis-related biomarkers in SONFH based on machine learning algorithm. The current findings may offer new insights on the SONFH pathogenesis, thereby providing new strategies for its diagnosis and treatment.
The SONFH dataset GSE123568 was downloaded from NCBI Gene Expression Omnibus database, which included 30 SONFH patients and 10 non-SONFH patients as control. The GPL15207 [PrimeView] Affymetrix Human Gene Expression Array platform was used for gene-chip assays. This dataset was used as analytical dataset.
For the above gene-chip dataset, the preprocessed, standardized and log2 transformed probe expression matrix were downloaded, and then, the annotation files were downloaded. Through one-to-one matching of probe and gene symbol, the probe that did not match gene symbol was deleted. For different probes that mapped to the same gene, the mean value was taken as the expression value of this gene for subsequent analysis.
Differentially expressed gene analysis
Based on the analysis dataset, the classical Bayesian method in limma 3.10.3  was used for differentially expressed gene (DEGs) analysis of SONFH vs. control. The p value was corrected by Benjamini & Hochberg algorithm. The differential expression thresholds were set as adj.p.value < 0.05 and |logFC (fold change)|< 0.5. After obtaining the DEGs, the volcano map and heatmap were drawn, respectively, for visual display.
Screening of SONFH correlated module genes
Weighted gene co-expression network analysis (WGCNA) is an analysis method to cluster the genes with similar expression patterns and then distinguish modules by gene expression similarity. The correlations between modules and modules, as well as modules and sample traits were calculated, so as to screen modules with highly correlated traits. Additionally, the genes in the significant module could be analyzed to find the target genes related to the research.
In order to find out the module genes highly related to diseases in the dataset GSE123568, all the genes were ranked in the dataset according to the variance from largest to smallest. Then, the genes with the top 25% variance (total number:12050) were selected and the disease status of the samples was taken as traits for analysis using the R package WGCNA1.61 .
SONFH-related ferroptosis-related gene (FRG) analysis
Firstly, FRGs were retrieved from FerrDb V2 . After that, the intersection of DEGs, FRGs, and disease-highly associated module genes obtained from WGCNA analysis was used as SONFH-related FRGs for subsequent analysis.
Functional enrichment analyses
Based on these disease-related FRGs, the GO (Gene Ontology) function  and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways  were analyzed using DAVID . The number of enriched genes in each term was set as at least 2. A p < 0.05 was considered as the threshold vale. The top 10 GO entries and top 20 KEGG pathways were selected for display.
Protein–protein interaction (PPI) network analysis
In order to understand the protein interactions between disease-related FRGs, the online database STRING 11.0  was used to predict and analyze whether there was any interaction between gene-encoded proteins. PPI score was set as 0.4. Cytoscape 3.4.0  was used to construct PPI networks. In addition, CytoNCA  plug-in 2.1.6 was used to analyze the topological properties, including “degree,” “betweenness” and “closeness of nodes” in the network.
Screening of key FRGs
SONFH-related FRGs (number: 43) screened above were further filtered by two machine learning algorithms as previously described [27, 28].
The expression values of disease-related FRGs in various samples were extracted. Then, LASSO algorithm was used to screen the feature genes combined with sample grouping. The glmnet package 4.2–2  of R 3.6.1 was used for regression analysis, and the parameter was set as nfold = 20.
The support vector machine (SVM) algorithm in R package "e1071" (version1.7–9)  was used to sort the disease-related FRGs. The recursive feature elimination (RFE) method was used to select the best gene from GSE123568 cohort to avoid the overfitting. Briefly, the importance ranking of each gene were achieved, and the error rate and accuracy of each iteration combination were obtained. The lowest error rate was selected as the best combination, and the corresponding gene was taken as the feature gene.
Finally, the overlapped feature genes screened from LASSO regression algorithm and RFE-SVM algorithm were obtained as the key FRGs.
Evaluation of diagnostic efficacy of key FRGs
For the key FRGs obtained from the above analyses, the expression values of the key genes in the data set were extracted, and then, the diagnostic ROC curves of the key genes were drawn by combining the sample grouping.
Gene set enrichment analysis (GSEA)
In dataset GSE123568, Pearson correlation coefficients between each key FRG and all the other genes were calculated, respectively, and then, the correlation coefficients were sorted from largest to smallest. After obtaining the related gene set for each FRG, GSEA based of KEGG pathway enrichment was conducted by using R clusterProfiler 3.8.1 . Meanwhile, Benjamini & Hochberg method was adopted to conduct multiple inspection correction, and adj.p < 0.05 was regarded as threshold. The top 5 KEGG pathways according to the ranking of significance were displayed.
Correlation analysis of key FRGs and immunity
Firstly, 28 kinds of immune cells and their corresponding marker genes were obtained from a previous literature . Further, based on the gene expression matrix of marker genes in all samples, ssGSEA algorithm was adopted and R package GSVA 1.36.2  was used to calculate the enrichment fraction of each immune cell. Then, the immune cells with differential infiltration levels between SONFH and control group were analyzed by Wilcoxon test. The differentially infiltrating immune cells with p < 0.05 were considered to be closely related to SONFH.
Then, Spearman method  was applied to analyze the correlation between key FRGs and infiltrating immune cells, and the corresponding p value and correlation coefficient cor were obtained.
Targeted drug prediction of key FRGs
The online database CTD  was used to search chemical interactions of the key genes. The drug–gene relationships supported by at least two references were selected. Cytoscape 3.4.0 was used for network building.
According to the cutoff value of adj.p.val < 0.05 and |logFC|> 0.5, 1380 up-regulated and 650 down-regulated DEGs between SONFH and control groups were selected. The differential volcano and heatmaps are shown in Fig. 1A, B.
Disease-related gene analysis
In order to observe the overall correlation of all samples in the dataset, the samples were clustered and the outlier samples were eliminated to ensure the accuracy of the analysis. The sample clustering based on gene expression is shown in Fig. 2A. The SONFH and control samples were clustered separately, so the samples were not eliminated. Then, in order to ensure that the interaction between genes met the scale-free distribution, we first determined the soft threshold of the data, as shown in Fig. 2B. We selected the optimal power value recommended by the WGCNA package, that was power = 25.
The adjacency among genes was analyzed, and the similarity was calculated based on the adjacency. Then, the differentiation coefficient among genes was derived, and the systematic clustering tree among genes was obtained. Then, according to the dynamic tree cutting standard, the fewest genes were set to 70 for each gene module, obtaining 22 modules. After that, MEDissThres was set to 0.2 to combine similar modules analyzed by dynamic tree algorithm. Cluster dendrogram of genes in SONFH is presented in Fig. 2C. Further, by calculating the correlation between the feature vector gene of each module and the clinical traits, 9 modules closely correlated with SONFH were obtained, of which cyan module (1425 genes) and grey60 module (136 genes) showed the strongest positive and negative correlation with the disease traits, respectively (Fig. 2D). Therefore, these two modules were regarded as the SONFH-closely related modules. The genes in the modules were considered as SONFH-related genes.
Disease-related FRGs analysis
Based on FerrDb V2 database, total 621 FRGs were retrieved, including 264 driver, 9 markers, 238 suppressor and 110 unclassified genes. After removing the duplicate genes, 564 FRGs were obtained.
The intersection of DEGs, FRGs, and module genes obtained from WGCNA analysis was taken, and 43 genes were obtained (Fig. 3A), which were considered as disease-related FRGs for subsequent analysis.
GO and KEGG pathway analysis showed that the 43 disease-related FRGs were enriched in 84 GO-BP (biological process), 8 GO-CC (cellular component), 10 GO-MF (molecular function) terms and 52 pathways. The top 10 GO and top 20 pathways are displayed in Fig. 3B, C. The disease-related FRGs were closely related to inflammation and apoptosis-related functions, as well as TNF signaling pathway, HIF-1 signaling pathway, and osteoclast differentiation.
Protein interaction network and correlation analysis
The protein interaction pairs of disease-related FRGs were analyzed, and 99 interaction pairs were identified, involving 30 proteins (Fig. 3D). The top ten nodes (hub nodes) were ILIB (degree, 15), HMOX1 (degree, 15), STAT3 (degree, 14), TLR4 (degree, 14), MAPK3 (degree, 14), PTGS2 (degree, 13), CREB1 (degree, 12), PTEN (degree, 12), CYBB (degree, 11), DUSP1 (degree, 10), and MAPK5 (degree, 10).
Screening of key FRGs
LASSO regression analysis was conducted based on the 43 disease-related FRGs, and 6 key feature genes (AKT1S1, ARNTL, BACH1, MGST1, SETD1B and NNMT) were obtained (Fig. 4A). Additionally, RFE-SVM algorithm was also used to screen key genes, and 17 feature genes were obtained, including BACH1, AKT1S1, SETD1B, KLF2, RICTOR, MGST1, ACSL1, SAT1, LRRFIP1, PTEN, FTH1, RBMS1, PTGS2, TXNIP, TMBIM4, NCF2, and DPEP1 (Fig. 4B). Finally, the intersection of feature genes obtained by LASSO regression algorithm and RFE-SVM algorithm was taken, and 4 intersection genes (AKT1S1, BACH1, MGST1 and SETD1B) were considered as key FRGs (Fig. 4C).
According to the method, diagnostic ROC curves of the 4 key FRGs were drawn, as shown in Fig. 4D. AUC areas were all above 0.85, indicating that the 4 genes had good disease prediction effect.
Potential mechanism of key genes explored by GSEA
The KEGG pathways of 4 key FRGs were explored by GSEA. Total 47 KEGG pathways were positively associated with AKT1S1 and 118 were negatively correlated with AKT1S1; 145 positively correlated and 34 negatively correlated KEGG pathways were obtained by BACH1 enrichment; MGST1 enrichment resulted in 118 positive and 16 negative KEGG pathways; 127 positive and 28 negative KEGG pathways were obtained by SETD1B enrichment. Top 5 KEGG pathways of each gene were selected for visual display, as shown in Fig. 5. The enrichment results are detailed in Additional file 1: Tables S1–S4. The osteoclast differentiation pathway was significantly enriched, indicating that this pathway may play a key role in SONFH.
Correlation analysis of key FRGs and immunity
The enrichment scores of 28 types of immune cells were calculated to compare the differences in immune cell infiltration levels between the disease group and the normal group. As shown in Fig. 6A, there were 20 immune cells with differential infiltration levels between the groups, such as activated B cell, activated CD8 T cell, central memory CD4 T cell, eosinophil, immature dendritic cell, and neutrophil. The correlation between differential immune cells is shown in Fig. 6B. Type 17 T helper cell and CD56 bright natural killer cell were positively correlated with most types of immune cells.
In order to further understand the correlation between the screened biomarkers (4 key FRGs) and differential immune cells, the Spearman method was applied, and the corresponding p value and correlation coefficient cor were obtained. The 4 key FRGs were correlated with most immune cells (Fig. 6C).
Targeted drug prediction
According to the threshold value, 41 drug–gene relationship pairs were finally obtained, including 34 drug molecules, such as phenobarbital, resveratrol, and clofibrate, and 4 key FRGs. Specially, carbon tetrachloride, bisphenol A, acetaminophen, benzo(a)pyrene, tetrachlorodibenzodioxin, and titanium dioxide had interactions with two genes (Fig. 7).
SONFH as a debilitating disease has become a public health burden. SONFH is usually asymptomatic; thus, the early diagnosis and prevention are imperative for SONFH patients. Emerging effects have been made for the identification of diagnostic marker for SONFH by bioinformatics approach. The study of Zhao et al. has applied WGCNA and CIBERSORT to screen the key diagnostic markers underlying the immune-related mechanism, such as TYROBP, TLR2 and P2RY13 . Ferroptosis is a newly discovered type of cell death, which has been found to be implicated in various pathophysiological processes, such as infection, inflammation and immunity . However, the diagnostic value of FRGs in SONFH has not clarified. Thus, in the present study, we attempted to identify the key FRGs valuable for the diagnosis and treatment for SONFH.
Our results showed that there were 4 key FRGs: AKT1S1, BACH1, MGST1 and SETD1B through machine learning algorithm. They were positively correlated with the pathway of osteoclast differentiation. Additionally, 20 immune cells showed differential infiltration levels between SONFH and control groups, and the 4 key FRGs were correlated with most immune cells. Findings of this study may help to understand the pathogenesis of SONFH as well as develop novel diagnosis and treatment targets.
AKT1S1 is a substrate of Akt and a component of the mTOR complex 1 . It has been reported that AKT1S1 can suppress the activity of mTOR complex 1 . mTOR is a member of mTOR complex 1, and the mTOR signaling pathway can regulate many stem cell processes, including cell survival, proliferation, and differentiation . Recent study revealed that the mTOR pathway may mediate the bone homeostasis . Activation of the mTOR signaling pathway can impair the osteogenic/adipogenic lineage differentiation of bone marrow mesenchymal stem cells (BMSCs) . Importantly, involvement of BMSCs in osteoblast lineage and bone formation is thought to be a major mediator of SONFH . In this study, AKT1S1 was down-regulated. We speculated that the down-regulation of AKT1S1 may increase the activity of mTOR complex 1, thereby resulting in SONFH.
All of BACH1, MGST1 and SETD1B were positively correlated with osteoclast differentiation. Bone remodeling is the result of a balance between osteoblast and osteoclast differentiation . Bone remodeling has played a key role in the development and severity of osteonecrosis of the femoral head . Nevertheless, little is known about the underlying molecular pathophysiology. Interestingly, the three key FRGs have not been investigated in SONFH to our best knowledge. Thus, we speculated that they may regulate osteoclast differentiation in SONFH.
Recently, more and more studies have suggested that infiltration of immune cell is associated with the occurrence of SONFH. For example, the activated B cells were found to have a higher frequency in SONFH patients peripheral blood compared with healthy controls . Inhibitory T cells are related to the pathogenesis of nontraumatic osteonecrosis of the femoral head by regulating the bone mass balance of the femoral head . In the present study, among the 28 immune cells, 20 were differential between SONFH and control groups, further suggesting the critical roles of immune cell infiltration in SONFH progression. Importantly, among these 20 immune cells, activated B cell was included and also had a higher infiltration level in SONFH group, being consistent with the report above. Mover, correlation analysis revealed that 4 key FRGs were correlated with most immune cells, suggesting that the identified key FRGs may play a role in SONFH pathogenesis through immunologic mechanisms.
Finally, we predicted the targeted drugs of the 4 FRGs. For instance, titanium dioxide had interactions with BACH1 and AKT1S1. Titanium dioxide is a common component of orthopedic prosthesis . Acetaminophen was predicted to have interactions with BACH1 and MGST1. Acetaminophen has a central analgesic effect and has been used in osteoarthritis . Investigating the interaction between key gene and targeted drug may help the development of medicine for SONFH.
Four key FRGs, AKT1S1, BACH1, MGST1 and SETD1B, were identified to play a critical role in the progression of SONFH through osteoclast differentiation and immunologic mechanisms. Additionally, all the 4 genes had good disease prediction effect and may act as biomarkers of the diagnosis and treatment of SONFH.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Differentially expressed genes
Steroid-induced osteonecrosis of the femoral head
Weighted gene co-expression network analysis
Kyoto Encyclopedia of Genes and Genomes
Support vector machine
Recursive feature elimination
van der Jagt D, Mokete L, Pietrzak J, Zalavras CG, Lieberman JR. Osteonecrosis of the femoral head: evaluation and treatment. J Am Acad Orthop Surg. 2015;23:69–70.
Zhang Q-Y, Li Z-R, Gao F-Q, Sun W. Pericollapse stage of osteonecrosis of the femoral head: a last chance for joint preservation. Chin Med J. 2018;131:2589–98.
Migliorini F, Maffulli N, Baroncini A, Eschweiler J, Tingart M, Betsch M. Prognostic factors in the management of osteonecrosis of the femoral head: a systematic review. Surgeon. 2022;21:85–98.
Migliorini F, La Padula G, Oliva F, Torsiello E, Hildebrand F, Maffulli N. Operative management of avascular necrosis of the femoral head in skeletally immature patients: a systematic review. Life. 2022;12:179.
Quaranta M, Miranda L, Oliva F, Aletto C, Maffulli N. Osteotomies for avascular necrosis of the femoral head. Br Med Bull. 2021;137:98–111.
Sadile F, Bernasconi A, Russo S, Maffulli N. Core decompression versus other joint preserving treatments for osteonecrosis of the femoral head: a meta-analysis. Br Med Bull. 2016;118:33–49.
Migliorini F, Maffulli N, Eschweiler J, Tingart M, Baroncini A. Core decompression isolated or combined with bone marrow-derived cell therapies for femoral head osteonecrosis. Expert Opin Biol Ther. 2021;21:423–30.
Song Y, Du Z, Ren M, Yang Q, Wang Q, Chen G, et al. Association of gene variants of transcription factors PPARγ, RUNX2, Osterix genes and COL2A1, IGFBP3 genes with the development of osteonecrosis of the femoral head in Chinese population. Bone. 2017;101:104–12.
Liu F, Wang W, Yang L, Wang B, Wang J, Chai W, et al. An epidemiological study of etiology and clinical characteristics in patients with nontraumatic osteonecrosis of the femoral head. J Res Med Sci. 2017;22:85.
Bernhard ME, Barnes CL, DeFeo BM, Kaste SC, Wang X, Lu Z, et al. Total hip arthroplasty in adolescents and young adults for management of advanced corticosteroid-induced osteonecrosis secondary to treatment for hematologic malignancies. J Arthroplasty. 2021;36:1352–60.
Li J, Cao F, Yin H-L, Huang Z-J, Lin Z-T, Mao N, et al. Ferroptosis: past, present and future. Cell Death Dis. 2020;11:1–13.
Hirschhorn T, Stockwell BR. The development of the concept of ferroptosis. Free Radical Biol Med. 2019;133:130–43.
Yoon B-H, Jones LC, Chen C-H, Cheng EY, Cui Q, Drescher W, et al. Etiologic classification criteria of ARCO on femoral head osteonecrosis part 1: glucocorticoid-associated osteonecrosis. J Arthroplasty. 2019;34(163–8): e1.
Mou Y, Wang J, Wu J, He D, Zhang C, Duan C, et al. Ferroptosis, a new form of cell death: opportunities and challenges in cancer. J Hematol Oncol. 2019;12:1–16.
Reichert CO, de Freitas FA, Sampaio-Silva J, Rokita-Rosa L, Barros PDL, Levy D, et al. Ferroptosis mechanisms involved in neurodegenerative diseases. Int J Mol Sci. 2020;21:8765.
Sun Y, Chen P, Zhai B, Zhang M, Xiang Y, Fang J, et al. The emerging role of ferroptosis in inflammation. Biomed Pharmacother. 2020;127: 110108.
Sun F, Lin Zhou J, Lin Liu Z, Wen Jiang Z, Peng H. Dexamethasone induces ferroptosis via P53/SLC7A11/GPX4 pathway in glucocorticoid-induced osteonecrosis of the femoral head. Biochem Biophys Res Commun. 2022;602:149–55.
Smyth GK, Ritchie M, Thorne N, Wettenhall J. LIMMA: linear models for microarray data. In: Bioinformatics and computational biology solutions using R and bioconductor. Statistics for biology and health. 2005.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinform. 2008;9:1–13.
Zhou N, Yuan X, Du Q, Zhang Z, Shi X, Bao J, et al. FerrDb V2: update of the manually curated database of ferroptosis regulators and ferroptosis-disease associations. Nucleic Acids Res. 2022;51:571–82.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25:25–9.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.
Sherman BT, Hao M, Qiu J, Jiao X, Baseler MW, Lane HC, et al. DAVID: a web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. 2022;50:216–21.
Mering CV, Huynen M, Jaeggi D, Schmidt S, Bork P, Snel B. STRING: a database of predicted functional associations between proteins. Nucleic Acids Res. 2003;31:258–61.
Kohl M, Wiese S, Warscheid B. Cytoscape: Software for Visualization and Analysis of Biological Networks. Data Mining in Proteomics. Springer; 2011. p. 291–303.
Tang Y, Li M, Wang J, Pan Y, Wu F-X. CytoNCA: a cytoscape plugin for centrality analysis and evaluation of protein interaction networks. Biosystems. 2015;127:67–72.
Li S, Zhou C, Xu Y, Wang Y, Li L, Pelekos G, et al. Similarity and potential relation between periimplantitis and rheumatoid arthritis on transcriptomic level: results of a bioinformatics study. Front Immunol. 2021;12: 702661.
Liu C, Liu R, Cao Z, Guo Q, Huang H, Liu L, et al. Identification of MMP9 as a novel biomarker to mitochondrial metabolism disorder and oxidative stress in calcific aortic valve stenosis. Oxidat. Med. Cell. Longev. 2022.
Jiang B, Sun P, Tang J, Luo B. Glmnet: graph learning-matching networks for feature matching. arXiv preprint arXiv:191107681. 2019.
Meyer D, Dimitriadou E, Hornik K, Weingessel A, Leisch F, Chang C-C et al. Package ‘e1071’. R J. 2019.
Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics J Integr Biol. 2012;16:284–7.
Charoentong P, Finotello F, Angelova M, Mayer C, Efremova M, Rieder D, et al. Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep. 2017;18:248–62.
Hänzelmann S, Castelo R, Guinney J. GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinform. 2013;14:1–15.
Zhao J, Zhang X, Guan J, Su Y, Jiang J. Identification of key biomarkers in steroid-induced osteonecrosis of the femoral head and their correlation with immune infiltration by bioinformatics analysis. BMC Musculoskelet Disord. 2022;23:67.
Davis AP, Grondin CJ, Johnson RJ, Sciaky D, Wiegers J, Wiegers TC, et al. Comparative toxicogenomics database (CTD): update 2021. Nucleic Acids Res. 2021;49:D1138–43.
Chen X, Kang R. Ferroptosis in infection, inflammation, and immunity. J Exp Med. 2021;218:e20210518.
Völkers M, Sussman M. mTOR/PRAS40 interaction: hypertrophy or proliferation. Taylor & Francis; 2013. p. 3579–80.
Laplante M, Sabatini DM. mTOR signaling at a glance. J Cell Sci. 2009;122:3589–94.
Xian L, Wu X, Pang L, Lou M, Rosen CJ, Qiu T, et al. Matrix IGF-1 maintains bone mass by activation of mTOR in mesenchymal stem cells. Nat Med. 2012;18:1095–101.
Chen C, Akiyama K, Wang D, Xu X, Li B, Moshaverinia A, et al. mTOR inhibition rescues osteopenia in mice with systemic sclerosis. J Exp Med. 2015;212:73–91.
Liu Y, Kou X, Chen C, Yu W, Su Y, Kim Y, et al. Chronic high dose alcohol induces osteopenia via activation of mTOR signaling in bone marrow mesenchymal stem cells. Stem Cells. 2016;34:2157–68.
Yu Z, Fan L, Li J, Ge Z, Dang X, Wang K. Lithium chloride attenuates the abnormal osteogenic/adipogenic differentiation of bone marrow-derived mesenchymal stem cells obtained from rats with steroid-related osteonecrosis by activating the β-catenin pathway. Int J Mol Med. 2015;36:1264–72.
Manolagas SC. Role of cytokines in bone resorption. Bone. 1995;17:S63–7.
Radke S, Battmann A, Jatzke S, Eulert J, Jakob F, Schütze N. Expression of the angiomatrix and angiogenic proteins CYR61 CTGF, and VEGF in osteonecrosis of the femoral head. J Orthop Res. 2006;24:945–52.
Zhang H, Xiao F, Liu Y, Zhao D, Shan Y, Jiang Y. A higher frequency of peripheral blood activated B cells in patients with non-traumatic osteonecrosis of the femoral head. Int Immunopharmacol. 2014;20:95–100.
Ma J, Ge J, Gao F, Wang B, Yue D, Sun W, et al. The role of immune regulatory cells in nontraumatic osteonecrosis of the femoral head: a retrospective clinical study, BioMed Res Int. 2019.
Borghi SM, Mizokami SS, Pinho-Ribeiro FA, Fattori V, Crespigio J, Clemente-Napimoga JT, et al. The flavonoid quercetin inhibits titanium dioxide (TiO2)-induced chronic arthritis in mice. J Nutr Biochem. 2018;53:81–95.
Towheed T, Maxwell L, Judd M, Catton M, Hochberg MC, Wells GA. Acetaminophen for osteoarthritis. Cochrane Database Syst Rev. 2006.
Ethics approval and consent to participate
This study was approved by the ethics committee of Sunshine Union Hospital, and informed consent has been signed by all participants. All methods were performed in accordance with relevant guidelines and regulations.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Liu, J., Han, X., Qu, L. et al. Identification of key ferroptosis-related biomarkers in steroid-induced osteonecrosis of the femoral head based on machine learning. J Orthop Surg Res 18, 327 (2023). https://doi.org/10.1186/s13018-023-03800-x
- Steroid-induced osteonecrosis of the femoral head
- Machine learning