Skip to main content

Pyroptosis -related potential diagnostic biomarkers in steroid-induced osteonecrosis of the femoral head

Abstract

Purpose

Steroid-induced necrosis of the femoral head (SONFH) is a refractory orthopedic hip disease occurring in young and middle-aged people, with glucocorticoids being the most common cause. Previous experimental studies have shown that cell pyroptosis may be involved in the pathological process of SONFH, but its pathogenesis in SONFH is still unclear. This study aims to screen and validate potential pyroptosis-related genes in SONFH diagnosis by bioinformatics analysis to further elucidate the mechanism of pyroptosis in SONFH.

Methods

There were 33 pyroptosis-related genes obtained from the prior reviews. The mRNA expression was downloaded from GSE123568 dataset in the Gene Expression Omnibus (GEO) database, including 10 non-SONFH (following steroid administration) samples and 30 SONFH samples. The pyroptosis-related differentially expressed genes involved in SONFH were identified with “affy” and “limma” R package by intersecting the GSE123568 dataset with pyroptosis genes. In addition, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the pyroptosis-related differentially expressed genes involved in SONFH were conducted by “clusterProfiler” R package and visualized by “GOplot” R package. Then, the correlations between the expression levels of the pyroptosis-related differentially expressed genes involved in SONFH were confirmed with “corrplot” R package. Moreover, the protein–protein interaction (PPI) network was analysed by using GeneMANIA database. Next, The ROC curve of pyroptosis-related differentially expressed genes were analyzed by “pROC” R package.

Results

A total of 10 pyroptosis-related differentially expressed genes were identified between the peripheral blood samples of SONFH patients and non-SONFH patients based on the defined criteria, including 20 upregulated genes and 10 downregulated genes. The GO and KEGG pathway enrichment analyses revealed that these 10 pyroptosis-related differentially expressed genes involved in SONFH were particularly enriched in cysteine-type endopeptidase activity involved in apoptotic process, positive regulation of interleukin-1 beta secretion and NOD-like receptor signaling pathway. Correlation analysis revealed significant correlations among the 10 differentially expressed pyroptosis-related genes involved in SONFH. The PPI results demonstrated that the 10 pyroptosis-related differentially expressed genes interacted with each other. Compared to non-SONFH samples, these pyroptosis-related differentially expressed genes had good predictive diagnostic efficacy (AUC = 1.000, CI = 1.000–1.000) in the SONFH samples, and NLRP1 had the highest diagnostic value (AUC: 0.953) in the SONFH samples.

Conclusions

There were 10 potential pyroptosis-related differentially expressed genes involved in SONFH were identified via bioinformatics analysis, which might serve as potential diagnostic biomarkers because they regulated pyroptosis. These results expand the understanding of SONFH associated with pyroptosis and provide new insights to further explore the mechanism of action and diagnosis of pyroptosis associated in SONFH.

Peer Review reports

Introduction

Steroid induced osteonecrosis of the femoral head (SONFH) is a chronic refractory orthopaedic disease, which has extremely high disability rate [1, 2]. SONFH seriously threaten human health, especially in young and middle-aged men aged 30 to 50. If not treated promptly and effectively, about 80% of SONFH patients would evolve into femoral head collapse within 1 to 4 years, causing a heavy social and economic burden [3, 4]. Due to the lack of effective screening tools and the difficulties in early diagnosis, many patients are already in the middle to late stage the onset of the corresponding symptoms and at the time of X-ray diagnosis. The main current treatments for SONFH are medication, hip-conserving surgery, and hip replacement surgery. Patients may need to face multiple revision total hip surgeries due to the presence of dislocation, aseptic loosening, infection, and limitations on the service life of the prosthesis. Therefore, reliable new predictive and prognostic models of SONFH are urgently needed for early diagnosis to make targeted therapy more feasible and achieve early treatment.

Pyroptosis, a novel pattern of non-apoptotic programmed cell death due to cellular inflammatory necrosis [5,6,7], which the morphological features are swelling of cells, forming pyroptotic bodies, cell membrane rupture and releasing inflammatory cytokines. Pyroptosis was widely found in eukaryotic cell organisms and was a self-protective mode of cells against external damage, but excessive activation could lead to damage. A large number of experimental studies have shown that pyroptosis was closely related to the development of diabetes mellitus, rheumatoid arthritis, ankylosing spondylitis, osteomyelitis, infection, Parkinson's disease, tumor and other diseases [8,9,10,11,12,13,14]. The extensive use of glucocorticoids could activate bone marrow mesenchymal stem cells (BMSC) pyroptosis, which was closely related to the development of SONFH [15]. The cellular pyroptosis pathway was a multifactorial and multi-level and complex regulatory process [16] Wang et al. showed that the expression of NLRP3 and caspase-1 were increased in human BMSCs treated with lipopolysaccharide (LPS), partly suggesting a pyroptosis in BMSCs [17]. Chen Y et al. reported that the NLRP3 inflammasome could mediate cell pyroptosis in BMSCs and triggered the imbalance in osteoblasts and osteoclast differentiation, which might play a critical role in the development of SONFH [12]. As demonstrated by Shahzad K et al., the caspase-1-dependent inflammasome activation occupies an important role in the development and development of SONFH [16]. Although many scholars at home and abroad have revealed the role of different types of pyroptosis in SONFH, the specific regulatory mechanism of pyroptosis in SONFH is still unclear, and the specific relationship between pyroptosis and SONFH needs further study. To date, there have been no bioinformatics studies on the genes and pathogenesis of SONFH associated with pyroptosis. A large sample of patient data helped to determine the adapted disease stage and prognostic value of therapeutic interventions, and Zhang Y et al. performed a genome-wide analysis of human peripheral serum from SONFH patients and healthy individuals and uploaded the data to the dataset GSE123568 in the National Center for Biotechnology Information Gene Expression Omnibus database (NCBI GEO). Therefore, we proposed to conduct systematic studies using bioinformatics and experimental validation means to determine the expression levels of cellular pyroptosis-related genes in normal and SONFH tissues, to explore the prognostic value of these genes,and understand the pathogenesis of pyroptosis in SONFH, and provide potential biomarkers for the clinical diagnosis and treatment of SONFH.

Materials and methods

Pyroptosis -related gene datasets and microarray data

The 33 pyroptosis-related genes were extracted from prior reviews [18]. We obtained the mRNA dataset of 10 non-SONFH (following steroid administration) samples and 30 SONFH samples from NCBI GEO (https://www.ncbi.nlm.nih.gov/geo/, ID: GSE123568).

Differentially expressed analysis of pyroptosis-related genes in SONFH

The “affy” and “limma” R package were used to identify the pyroptosis-related differentially expressed genes between non-SONFH samples and SONFH samples in GSE123568. The screening criterion as an adjusted P-value < 0.01 and |log2 fold change (FC)|> log21.5. Additionally, the results were computed and visualized in the heatmaps and volcano plots by “pheatmap” and “ggplot2” R package.

Functional enrichment analysis of pyroptosis-related differentially expressed genes in SONFH

Based on these pyroptosis-related differentially expressed genes in SONFH, Gene Ontology (GO) [19] and Kyoto Encyclopedia of Genes and Genomes (KEGG) [20] pathway enrichment analysis were carried out by “clusterProfiler” R package and visualized by “GOplot” R package. The GO enrichment analysis mainly included biological process (BP), cellular composition (CC), and molecular function (MF). The statistically significant difference of enriched terms was filtered according to specific criteria of an adjusted P-value < 0.05.

Correlation analysis of the pyroptosis-related differentially expressed genes in SONFH

To explore the intrinsic relationship among these pyroptosis-related differentially expressed genes in SONFH, the Spearman correlation coefficient among each gene were calculated by the “corrplot” R package.

Construction of PPI network of the pyroptosis-related differentially expressed genes in SONFH

To gain more insight into the in-depth relationships of pyroptosis-related differentially expressed genes in SONFH, a PPI network was predicted and constructed by the GeneMANIA database (http://www.genemania.org), [21] and which was constructed and visualized by Cytoscape software (version: 3.7.2, http://cyto-scape.org/). The minimum required interaction score for the PPI analysis was set at in a combined score >  = 0.4. The nodes represented genes, and the edges represented the interactive relations between the genes. and the isolated nodes were discarded.

ROC curve analysis of the pyroptosis-related differentially expressed genes in SONFH

The “pROC” R package were employed to perform a ROC curve analysis of these pyroptosis-related differentially expressed genes in SONFH. The area under the curve (AUC) was calculated to analyse the intrinsic effectiveness of diagnostic tests, which was an indicator combining sensitivity and specificity between 0.5 and 1. The closer the AUC was to 1, the better the diagnosis was. The AUC had low accuracy at 0.5–0.7, and moderate accuracy at 0.7–0.9, and high accuracy above 0.9.

Statistical analysis

All statistical analysis was accomplished with R software (version 4.0.2) in this study. Gene expression levels of samples in GSE123568 were employed by the Kruskal–Wallis test. The data are shown as means ± standard deviation. Independent-sample t tests were used to compare means between two different groups. Fisher's exact test was used to compare the incidence of disease between two different groups. Statistical analysis was conducted using GraphPad Prism (Version 9.0.0). A threshold for statistical significance was set as P < 0.05.

Results

Pyroptosis-related differentially expressed genes in SONFH

A total of 33 non-duplicated genes were excavated referring to the article of Ye et al., and a total of 31 expression of pyroptosis-related genes were intersected them with GSE123568. There were 10 pyroptosis-related differentially expressed genes were identified based on the defined criteria of an adjusted P-value < 0.05 and |logFC|> log21.5, including 9 up-regulated genes and 1 down-regulated genes in SONFH samples compared to non-SONFH Samples (following steroid administration) in Table 1.

Table 1 The 10 pyroptosis-related differentially expressed genes in SONFH samples compared to non-SONFH samples

The 10 pyroptosis-related differentially expressed genes of GSE123568 were presented as the volcano plots and the circo heatmap by the ggplot and pheatmap package R software in Fig. 1A and B. Immediately after, the expression patterns of 10 pyroptosis-related genes differentially expressed between SONFH and non-SONFH samples were shown in bean plot in Fig. 2.

Fig. 1
figure 1

Pyroptosis-related differentially expressed genes in SONFH and non-SONFH samples. A Volcano plot of the pyroptosis-related genes intersected with GSE123568. The red dots represent the significantly up-regulated genes and the blue dots indicate the significantly down-regulated genes. C Circo heatmap of the 10 pyroptosis-related differentially expressed genes in SONFH and non-SONFH samples

Fig. 2
figure 2

Bean plot of the 10 differentially expressed pyroptosis -related genes in SONFH and non-SONFH samples

GO and KEGG enrichment analysis of pyroptosis-related differentially expressed genes

We utilized the “clusterProfiler” R package to conducte GO and KEGG enrichment analyses of these pyroptosis-related differentially expressed genes. There were 481 significantly enriched BP terms, which were involved in activation of cysteine-type endopeptidase activity involved in apoptotic process, positive regulation of interleukin-1 beta secretion, interleukin-1 beta production, positive regulation of interleukin-1 secretion, and interleukin-1 production, as shown in Fig. 3A, B, D, and Table 2. A total of 9 significantly enriched CC terms were mainly correlated with inflammasome complex, cytosolic part, azurophil granule lumen, primary lysosome and azurophil granule, as shown in Fig. 3A. There were 22 significantly enriched MF terms including cysteine-type endopeptidase activity involved in apoptotic process, cysteine-type endopeptidase activator activity involved in apoptotic process, CARD domain binding, peptidase activator activity involved in apoptotic process and cysteine-type endopeptidase activity, as shown in Fig. 3A. In the KEGG enrichment analysis, these pyroptosis-related differentially expressed genes were significantly enriched in 13 KEGG pathway terms, such as NOD-like receptor signaling pathway, Salmonella infection, Pertussis, Pathogenic Escherichia coli infection and Shigellosis, as shown in Fig. 3A, C, E, and Table 2.

Fig. 3
figure 3

GO analysis and KEGG pathway functional enrichment analysis of 10 pyroptosis-related differentially expressed genes in SONFH and non-SONFH samples. A Bar plot of the top 10 enriched BP terms, CC terms MF terms and KEGG pathway terms of these 30 genes. B Circle plot of the top 10 enriched BP terms of these 10 genes. C Circle plot of top 10 enriched KEGG terms of these 10 genes. D Chord plot of relationship among the top 10 enriched BP terms and targets. E Chord plot of relationship among the top 10 enriched KEGG pathway terms and targets. The colours of the nodes range from red to blue in descending order of logFC values. The genes are ordered according to logFC values

Table 2 The GO and KEGG Enrichment Analysis of 10 pyroptosis-related differentially expressed genes in SONFH samples compared to non-SONFH samples

Correlation analysis of pyroptosis-related differentially expressed genes

To explore the expression correlation among these pyroptosis-related differentially expressed genes in the GSE123568 dataset, we carried out the correlation analysis. The results suggested that some of these pyroptosis-related genes showed weak-to-moderate correlations. Among these 10 genes, CASP1 and CASP4 were most positively correlated (Cor = 0.92), while NLRP1 and ELANE (Cor =  − 0.24) were most negatively correlated as shown in Fig. 4.

Fig. 4
figure 4

Spearman correlation analysis of the 10 differentially expressed pyroptosis-related genes in SONFH and non-SONFH samples

PPI network of pyroptosis-related differentially expressed genes

We conducted a PPI analysis by using the GeneMANIA database, to determine the in-depth interactions of these pyroptosis-related differentially expressed genes. As shown in Fig. 5, the results of GeneMANIA also revealed that the functions of these pyroptosis-related genes interacted with each other and the number of interactions of each gene, which was found that 57.91% had co-expression, 16.47% physical interactions, 9.02% pathway, 8.49% predicted, 7.74% shared protein domians, and 0.37% co-localization.

Fig. 5
figure 5

PPI network of 10 pyroptosis-related differentially expressed genes in SONFH and non-SONFH samples. Black nodes represent target proteins, and connecting colors indicate different correlations. Functional associations between targets were investigated using GeneMANIA. Genes in black circles were query terms while these in gray circle indicate genes associated with query genes

ROC curve analysis of pyroptosis-related differentially expressed genes

The ROC curve analysis of these pyroptosis-related differentially expressed genes of SONFH samples were analyzed, established and drawn by “pROC” and “ggplot2” R package. Area under the curve (AUC) values indicated an indicator combining sensitivity and specificity metrics, which could describe the intrinsic effectiveness of diagnostic tests. Compared to non-SONFH samples, these pyroptosis-related differentially expressed genes had good predictive diagnostic efficacy (AUC = 1.000, CI = 1.000–1.000) in the SONFH samples as shown in Fig. 6A, which might be potential pyroptosis-related biomarkers for SONFH based on our present samples. Among them, NLRP1 had the highest diagnostic value (AUC: 0.953) in the SONFH samples. The diagnostic values of other genes were also performed as follows in SONFH samples as shown in Fig. 6B-K: NLRP3 (AUC: 0.897), CASP4 (AUC:0.830), CASP1 (AUC:0.827), CASP5 (AUC:0.800), PYCARD (AUC: 0.793), SCAF11 (AUC: 0.787), NOD2 (AUC: 0.763), ELANE (AUC: 0.730) and IL1B (AUC:0.723).

Fig. 6
figure 6

ROC curve of the 10 pyroptosis-related differentially expressed genes in SONFH samples

Discussion

The most common cause of non-traumatic osteonecrosis of the femoral head is alcohol and glucocorticoids [22]. Since Pietrogrande et al. reported the first case of SONFH [23], glucocorticoids as one of the pathogenic factors have been shown to be directly associated with the onset of SONFH [24]. With the wide application of glucocorticoids in clinical practice, the incidence of SONFH has increased year by year, which has surpassed alcoholic osteonecrosis to become the first major cause of non-traumatic osteonecrosis [25]. Through a follow-up study, Li Zirong et al. found that most SONFH occurred within 2 years after medication, mainly within 6 months. Moreover, high-dose glucocorticoids shock therapy and long-term application of glucocorticoids were closely related to SONFH, while the short-term application of low-dose glucocorticoids had no obvious correlation with SONFH [25,26,27,28]. Patients in the early stage of SONFH did not have any clear symptoms, and the symptoms cannot be reversed, and the incidence rate increased year by year [29]. But some studies have shown that under the action of glucocorticoid, normal connection mechanism, signal transmission, osteocytes, osteoblasts cannot according to the external stimulus orderly repair, reshape the normal structure, the mechanical properties of the femoral head cannot effectively reconstruct, and cannot restore the normal femoral head mechanical support structure, and exhibit local sclerosis cystic changes, and even collapse [30,31,32,33,34]. Therefore, the early and accurate timing of hip-preserving treatment is the key to the success treatment of SONFH [25].

Although glucocorticoids have been recognized as the pathological cause of SONFH, the exact mechanisms underlying their specific pathogenesis and progression are not fully understood, which was a dynamic pathological process jointly regulated by multi-target, multi-pathway, multi-link, and multi-stage cooperation, involving the interaction of various molecular and biological processes [35]. With the further deepening of basic and clinical research of SONFH, the increasingly highly recognized pathogenic mechanisms were proposed to explain the pathogenic mechanisms, including intravascular coagulation, uncoupled bone remodeling, lipid metabolism, abnormal differentiation of MSCs, oxidative stress injury, programmed cell death, and so on [36]. Among them, changes in osteocyte cell activity due to programmed cell death were the intermediate link in various theories. Programmed cell death was an active process controlled by relevant genes that regulates body development and maintains a stable internal environment [37]. Many scholars believe that programmed cell death may play a key role in the process of SONFH. Pyroptosis was a programmed cell necrosis different from apoptosis and autophagic inflammation and proposed by Boise in 2001, which is a natural immune barrier for the human body [7]. The morphological aspect of pyroptosis is characterized by swelling of cells, formation of pyroptotic bodies, membrane rupture and releasing inflammatory factors [38]. At present, it has been confirmed that the extensive use of GC could activate the pyroptosis of BMSCs, which is closely related to the occurrence of SONFH [15]. It was found that high Nlrp3D301N expression causes bone mass loss, mainly due to the role of inflammatory factors such as IL-1β, facilitated by NLRP3 in mouse models specifically expressing Nlrp3D301N in osteoclasts, to further activate osteoclasts to cause bone resorption [39]. Liu et al. demonstrated the presence of pyroptosis from the perspective of LPS inducing oxidative stress in osteoblasts, and the inhibition of pyroptosis enhances osteoblast differentiation and activity [40]. Glucocorticoids on the signal regulation of programmed cell death, especially cell pyroptosis, was a complex process, cells in different states of glucocorticoid effect, under a certain concentration of glucocorticoids could effectively promote bone tissue related cell proliferation activity, rebuild normal signaling network. Tao et al. indicated that pyroptosis was proposed as a novel hypothesis of bone-related diseases pathogenesis for the first time, thus providing a new direction for the treatment of bone-related diseases in the future [41]. Therefore, this study focuses on the expression level changes of key pyroptosis-related genes in SONFH, and the achievement of tone results will have important theoretical and practical significance for the pathological mechanism, early diagnosis and even prognosis of SONFH.

To the best of our knowledge, no relevant articles have yet identified pyroptosis-related genes in SONFH based on bioinformatics analysis and animal experimental validation so far. There were 10 pyroptosis-related differentially expressed genes associated with SONFH based on the intersection of the 33 non-duplicated pyroptosis-related genes referring to the article of Ye et al. and differentially expressed genes from the GSE123568 dataset according to the defined criteria of an adjusted P-value < 0.05 and |logFC|> log21.5, including 9 up-regulated genes and 1 down-regulated genes. The 10 pyroptosis-related differentially expressed genes in SONFH samples compared to non-SONFH samples included CASP1, CASP4, CASP5, ELANE, IL1B, NLRP1, NLRP3, NOD2, PYCARD and SCAF11.

To further elucidate the underlying mechanism of the 10 pyroptosis-related differentially expressed genes in SONFH, we evaluated their biological functions and signaling pathways GO and KEGG enrichment analysis by the “clusterProfiler” R package. The GO and KEGG enrichment analyses showed that the 481 BP terms, the 9 CC terms, the 22 MF terms and 13 KEGG terms were significantly enriched, such as cysteine-type endopeptidase activity involved in apoptotic process, positive regulation of interleukin-1 beta secretion, and NOD-like receptor signaling pathway, and so on. The 10 genes involved in pyroptosis were further carried out by the correlation analysis, and the result showed CASP1 and CASP4 were most positively correlated, while NLRP1 and ELANE were most negatively correlated. At the same time, we constructed the PPI network on the above differentially expressed genes using the GeneMANIA database. The PPI networks might consist of several functional modules, which had co-expression (57.91%), physical interactions (16.47%), pathway(9.02%), predicted(8.49%), shared protein domians(7.74%) and co-localization(0.37%). Besides, compared to non-SONFH samples, the diagnostic values of these above 10 genes were performed as follows: NLRP1, NLRP3, CASP4, CASP1, CASP5, PYCARD, SCAF11, NOD2, ELANE and IL1B; and these 10 pyroptosis-related differentially expressed genes had higher diagnostic value (AUC = 1.000, CI = 1.000–1.000) in the SONFH samples, which might be potential pyroptosis-related biomarkers for SONFH.

Some of these SONFH and pyroptosis-related genes have been previously studied. For example, Yu and his colleagues showed that IL1B gene influences the genetic susceptibility of SONFH [42] Li et al. identified that NOD2 might be potential inflammation-related biomarkers of SONFH by weighted gene co-expression network analysis (WGCNA) and inflammation-related genes (IRGs) from the Molecular Signatures Database (MSigDB) [43]. Further studies are needed to reveal the roles and mechanisms of potential pyroptosis-related differentially expressed genes in SONFH.

Although the pyroptosis-related genes potentially involved in SONFH have been mined via bioinformatics and validation of animal experiments, there were still some limitations in this study. First, although there was a large amount of high-throughput data stored in the GEO database, we still did not find more homo spapiens datasets about SONFH, which led to a particularly limited number of samples. Second, further experiments of human blood specimens are lacking to pyroptosis-related validate the differentially expressed genes and their accuracy and specificity as diagnostic biomarkers. Third, the specific their underlying molecular mechanisms of differential genes in the pathogenesis and progression of SONFH are still unclear, and more basic experiments were needed to elaborate. Therefore, further research is necessary to address the possible limitations in terms of biased results and conclusions.

Conclusion

In conclusion, there were 10 potential pyroptosis-related differentially expressed genes involved in SONFH identified based on bioinformatic analysis in the present study, which might serve as potential diagnostic biomarkers of pyroptosis after SONFH. These results would help to further explore the development of SONFH by regulating pyroptosis, and to provide data support and new research ideas for the diagnosis of SONFH.

Availability of data and materials

Publicly available datasets were analysed in this study. These datasets can be found at the following URLs: GSE123568 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE123568).

Abbreviations

AUC:

Area Under the Curve

BP:

Biological Process

BMSC:

Bone Marrow Mesenchymal Stem Cells

CC:

Cellular Component

FC:

Fold Change

GC:

Glucocorticoids

GEO:

Gene Expression Omnibus

GO:

Gene Ontology

KEGG:

Kyoto Encyclopedia of Genes and Genomes

LPSL:

Lipopolysaccharide

MF:

Molecular Function

NCBI:

National Center for Biotechnology Information

PPI:

Protein–Protein Interaction

SD:

Sprague–Dawley

SONFH:

Steroid-Induced Osteonecrosis of the Femoral Head

References

  1. Hauzeur JP, Malaise M, de Maertelaer V. A prospective cohort study of the clinical presentation of non-traumatic osteonecrosis of the femoral head: spine and knee symptoms as clinical presentation of hip osteonecrosis [J]. Int Orthop. 2016;40(7):1347–51.

    Article  PubMed  Google Scholar 

  2. Zhang QY, Li ZR, Gao FQ, et al. Pericollapse Stage of Osteonecrosis of the Femoral Head: A Last Chance for Joint Preservation [J]. Chin Med J. 2018;131(21):2589–98.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Mont MA, Zywiel MG, Marker DR, et al. The natural history of untreated asymptomatic osteonecrosis of the femoral head: a systematic literature review [J]. J Bone Joint Surg Am. 2010;92(12):2165–70.

    Article  PubMed  Google Scholar 

  4. Gou WL, Lu Q, Wang X, et al. Key pathway to prevent the collapse of femoral head in osteonecrosis [J]. Eur Rev Med Pharmacol Sci. 2015;19(15):2766–74.

    PubMed  Google Scholar 

  5. Kovacs SB, Miao EA. Gasdermins: Effectors of Pyroptosis [J]. Trends Cell Biol. 2017;27(9):673–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Kumar S. Caspase function in programmed cell death [J]. Cell Death Differ. 2007;14(1):32–43.

    Article  CAS  PubMed  Google Scholar 

  7. Boise LH, Collins CM. Salmonella-induced cell death: apoptosis, necrosis or programmed cell death? [J]. Trends Microbiol. 2001;9(2):64–7.

    Article  CAS  PubMed  Google Scholar 

  8. Vande Walle L, Van Opdenbosch N, Jacques P, et al. Negative regulation of the NLRP3 inflammasome by A20 protects against arthritis [J]. Nature, 2014, 512(7512): 69–73.

  9. Tian ZG, Yao M, Chen J. Micheliolide alleviates ankylosing spondylitis (AS) by suppressing the activation of the NLRP3 inflammasome and maintaining the balance of Th1/Th2 via regulating the NF-κB signaling pathway [J]. Ann Transl Med. 2020;8(16):991.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Zhu X, Zhang K, Lu K, et al. Inhibition of pyroptosis attenuates Staphylococcus aureus-induced bone injury in traumatic osteomyelitis [J]. Ann Transl Med. 2019;7(8):170.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. van Opdenbosch N, Lamkanfi M. Caspases in Cell Death, Inflammation, and Disease [J]. Immunity. 2019;50(6):1352–64.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Chen Y, Qin X, An Q, et al. Mesenchymal Stromal Cells Directly Promote Inflammation by Canonical NLRP3 and Non-canonical Caspase-11 Inflammasomes [J]. EBioMedicine. 2018;32:31–42.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Li YF, Huang X, Li X, et al. Caspase-1 mediates hyperlipidemia-weakened progenitor cell vessel repair [J]. Front Bioscience (Landmark edition). 2016;21(1):178–91.

    Article  CAS  Google Scholar 

  14. Yu ZW, Zhang J, Li X, et al. A new research hot spot: The role of NLRP3 inflammasome activation, a key step in pyroptosis, in diabetes and diabetic complications [J]. Life Sci. 2020;240: 117138.

    Article  CAS  PubMed  Google Scholar 

  15. Wang Duoxian, Cao Linzhong, Jiang Wei, et al. Inflammasome activation and pyrolysis in hormone-induced femoral head necrosis. Chem Life. 2021;41(08): 1733–9.( in Chinese)

  16. Shahzad K, Bock F, Al-Dabet MM, et al. Caspase-1, but Not Caspase-3, promotes diabetic nephropathy [J]. J Am Soc Nephrol. 2016;27(8):2270–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Wang L, Chen K, Wan X, et al. NLRP3 inflammasome activation in mesenchymal stem cells inhibits osteogenic differentiation and enhances adipogenic differentiation [J]. Biochem Biophys Res Commun. 2017;484(4):871–7.

    Article  CAS  PubMed  Google Scholar 

  18. Ye Y, Dai Q, Qi H. A novel defined pyroptosis-related gene signature for predicting the prognosis of ovarian cancer [J]. Cell death discovery. 2021;7(1):71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Ashburner M, Ball C A, Blake J A, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium [J]. Nat Genet, 2000;25(1): 25–9.

  20. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes [J]. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Warde-Farley D, Donaldson S L, Comes O, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function [J]. Nucleic acids research, 2010, 38(Web Server issue): W214–20.

  22. Mont MA, Jones LC, Hungerford DS. Nontraumatic osteonecrosis of the femoral head: ten years later [J]. J Bone Joint Surg Am. 2006;88(5):1117–32.

    PubMed  Google Scholar 

  23. Pietrogrande V, Mastromarino R. Osteopatia da prolunga trattamento cortisonico. Ortop Traumatol. 1957;25:791.

    Google Scholar 

  24. Zizic TM, Marcoux C, Hungerford DS, et al. Corticosteroid therapy associated with ischemic necrosis of bone in systemic lupus erythematosus [J]. Am J Med. 1985;79(5):596–604.

    Article  CAS  PubMed  Google Scholar 

  25. Fu Weimin, Liu Baoyi, Wang Benjie, et al. Preserving hip treatment of steroid⁃induced osteonecrosis of the femoral head. Chin J Orthop. 2019, (23): 1424–5–6–7–8–9–30–31. ( in Chinese)

  26. Li Zirong. Neconecrosis of femoral head: early diagnosis and individualized treatment .Orthop J China. 2013;21(19):1909–11. ( in Chinese)

  27. Chinese Guideline for the Diagnosis and Treatment of Osteonecrosis of the Femoral Head in Adults [J]. Orthop Surg, 2017;9(1): 3–12.

  28. Murakami J, Ishii M, Suehiro F, et al. Vascular endothelial growth factor-C induces osteogenic differentiation of human mesenchymal stem cells through the ERK and RUNX2 pathway [J]. Biochem Biophys Res Commun. 2017;484(3):710–8.

    Article  CAS  PubMed  Google Scholar 

  29. Zhao D, Cui D, Wang B, et al. Treatment of early stage osteonecrosis of the femoral head with autologous implantation of bone marrow-derived and cultured mesenchymal stem cells [J]. Bone. 2012;50(1):325–30.

    Article  PubMed  Google Scholar 

  30. Zhang SY, Wang F, Zeng XJ, et al. Astragalus polysaccharide ameliorates steroid-induced osteonecrosis of femoral head through miR-206/HIF-1α/BNIP3 axis [J]. Kaohsiung J Med Sci. 2021;37(12):1089–100.

    Article  CAS  PubMed  Google Scholar 

  31. Xu K, Lu C, Ren X, et al. Overexpression of HIF-1α enhances the protective effect of mitophagy on steroid-induced osteocytes apoptosis [J]. Environ Toxicol. 2021;36(11):2123–37.

    Article  CAS  PubMed  Google Scholar 

  32. Huang Z, Wang Q, Zhang T, et al. Hyper-activated platelet lysates prevent glucocorticoid-associated femoral head necrosis by regulating autophagy [J]. Biomed Pharmacother. 2021;139.

    Article  CAS  PubMed  Google Scholar 

  33. Petek D, Hannouche D, Suva D. Osteonecrosis of the femoral head: pathophysiology and current concepts of treatment [J]. EFORT Open Rev. 2019;4(3):85–97.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Liu GB, Li R, Lu Q, et al. Three-dimensional distribution of cystic lesions in osteonecrosis of the femoral head [J]. J Orthop Transl. 2020;22:109–15.

    Google Scholar 

  35. Li JC, Liang XZ, Luo D, et al. Study on the molecular mechanism of BuShenHuoXue capsule in treatment of steroid-induced osteonecrosis of the femoral head [J]. Ann Transl Med. 2020;8(24):1680.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Luo P, Gao F, Han J, et al. The role of autophagy in steroid necrosis of the femoral head: a comprehensive research review [J]. Int Orthop. 2018;42(7):1747–53.

    Article  PubMed  Google Scholar 

  37. Yoon-Seok Y, Soo-Youn L. SOO-HO L J C O S. Apoptosis in the Osteonecrosis of the Femoral Head [J]. 2010;2(4):250–5.

    Google Scholar 

  38. Ma Zijian, Ma mingling, Dong Hui, et al. the effect and prospect of pyroptosis on bone homeostasis .Chin J Osteoporos. 2022, 28(06): 922–6. ( in Chinese)

  39. Qu C, Bonar SL, Hickman-Brecks CL, et al. NLRP3 mediates osteolysis through inflammation-dependent and -independent mechanisms [J]. FASEB J. 2015;29(4):1269–79.

    Article  CAS  PubMed  Google Scholar 

  40. Liu S, Du J, Li D, et al. Oxidative stress induced pyroptosis leads to osteogenic dysfunction of MG63 cells [J]. J Mol Histol. 2020;51(3):221–32.

    Article  CAS  PubMed  Google Scholar 

  41. Tao Z, Wang J, Wen K, et al. Pyroptosis in Osteoblasts: A Novel Hypothesis Underlying the Pathogenesis of Osteoporosis [J]. Front Endocrinol. 2020;11: 548812.

    Article  Google Scholar 

  42. Yu Y, Zhang Y, Wu J, et al. Genetic polymorphisms in IL1B predict susceptibility to steroid-induced osteonecrosis of the femoral head in Chinese Han population [J]. Osteoporosis Int. 2019;30(4):871–7.

    Article  CAS  Google Scholar 

  43. Li T, Huang C, Ma J, et al. Identification of Inflammation-Related Genes and Exploration of Regulatory Mechanisms in Patients with Osteonecrosis of the Femoral Head [J]. Biomed Res Int. 2022;2022:4501186.

    PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We would like to acknowledge the following financial support: the National Natural Science Foundation of China (No. 82205154 and No. 82074453); the National Natural Science Foundation of Shandong Province (No. ZR2021QH004 and No. ZR2021LZY002); the Development Plan of Shandong Medical and Health Technology (No. 2019WS577 and No.202104070600); the Development Plan of Shandong Traditional Chinese Medicine Science and Technology (No. 2020Q009 and No. 2020Q130); and the Innovation Program of Jinan Clinical Medicine Science and Technology (No.202019056). We thank the NCBI GEO database for sharing the data.

Funding

This work was supported by grants from the National Natural Science Foundation of China (No. 82205154 and No. 82074453); the National Natural Science Foundation of Shandong Province (No. ZR2021QH004 and No. ZR2021LZY002); the Development Plan of Shandong Medical and Health Technology (No. 2019WS577 and No.202104070600); the Development Plan of Shandong Traditional Chinese Medicine Science and Technology (No. 2020Q009 and No. 2020Q130); and the Innovation Program of Jinan Clinical Medicine Science and Technology (No.202019056).

Author information

Authors and Affiliations

Authors

Contributions

All authors made a significant contribution to the work reported and agreed to be accountable for all aspects of the work. C.J.L., L.X.Z., and W.P. designed the experiments. C.J.L., L.B.W., D.H.T., W.M.T. and L.X.Z performed the experiments. C.J.L. and L.X.Z. prepared the initial draft of the manuscript. L.X.Z and W.P. gave critical feedback during the study or during the submission of the manuscript. All authors provided final approval of the version to be submitted and agreed on the journal for publication.

Corresponding authors

Correspondence to Xue-Zhen Liang or Ping Wang.

Ethics declarations

Ethics approval and consent to participate

The study was conducted in agreement with the Declaration of Helsinki and its later amendments or comparable ethical standards.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Chai, JL., Lu, BW., Du, HT. et al. Pyroptosis -related potential diagnostic biomarkers in steroid-induced osteonecrosis of the femoral head. BMC Musculoskelet Disord 24, 609 (2023). https://doi.org/10.1186/s12891-023-06729-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12891-023-06729-8

Keywords