Purpose: Steroid-induced osteonecrosis of the femoral head (SONFH) is a refractory orthopaedic hip joint disease that occurs in young- and middle-aged people. Previous experimental studies have shown that autophagy might be involved in the pathological process of SONFH, but the pathogenesis of autophagy in SONFH remains unclear. We aimed to identify and validate the key potential autophagy-related genes involved in SONFH to further illustrate the mechanism of autophagy in SONFH through bioinformatics analysis. Methods: The GSE123568 mRNA expression profile dataset, including 10 non-SONFH (following steroid administration) samples and 30 SONFH samples, was downloaded from the Gene Expression Omnibus (GEO) database. Autophagy-related genes were obtained from the Human Autophagy Database (HADb). The autophagy-related genes involved in SONFH were screened by intersecting the GSE123568 dataset with the set of autophagy genes. The differentially expressed autophagy-related genes involved in SONFH were identified with R software. In addition, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of the differentially expressed autophagy-related genes involved in SONFH were conducted by using R software. Then, the correlations between the expression levels of the differentially expressed autophagy-related genes involved in SONFH were confirmed with R software. Moreover, the protein–protein interaction (PPI) network was analysed by using the Search Tool for the Retrieval of Interacting Genes (STRING), significant gene cluster modules were identified with the MCODE Cytoscape plugin, and hub genes among the differentially expressed autophagy-related genes involved in SONFH were screened by using the CytoHubba Cytoscape plugin. Finally, the expression levels of the hub genes of the differentially expressed autophagy-related genes involved in SONFH were validated in hip articular cartilage specimens from necrotic femur heads (NFHs) by using the GSE74089 dataset and further verification by qRT-PCR. Results: A total of 34 differentially expressed autophagy-related genes were identified between the peripheral blood samples of SONFH patients and non-SONFH patients based on the defined criteria, including 25 upregulated genes and 9 downregulated genes. The GO and KEGG pathway enrichment analyses revealed that these 34 differentially expressed autophagy-related genes involved in SONFH were particularly enriched in death domain receptors, the FOXO signalling pathway and apoptosis. Correlation analysis revealed significant correlations among the 34 differentially expressed autophagy-related genes involved in SONFH. The PPI results demonstrated that the 34 differentially expressed autophagy-related genes interacted with each other. Ten hub genes were identified by using the MCC algorithms of CytoHubba. The GSE74089 dataset showed that TNFSF10, PTEN and CFLAR were significantly upregulated while BCL2L1 was significantly downregulated in the hip cartilage specimens, which was consistent with the GSE123568 dataset. TNFSF10, PTEN and BCL2L1 were detected with consistent expression by qRT-PCR. Conclusions: Thirty-four potential autophagy-related genes involved in SONFH were identified via bioinformatics analysis. TNFSF10, PTEN and BCL2L1 might serve as potential drug targets and biomarkers because they regulate autophagy. These results expand the autophagy-related understanding of SONFH and might be useful in the diagnosis and prognosis of SONFH.

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.