Identification and analysis of mitochondria-related central genes in steroid-induced osteonecrosis of the femoral head, along with drug prediction
Authors: Ma Zheru , Sun Jing , Jiang Qi , Zhao Yao , Jiang Haozhuo , Sun Peng , Feng Wei 
Published: February 2024

Purpose

Steroid-induced osteonecrosis of the femoral head (SONFH) is a refractory orthopedic hip joint disease that primarily affects middle-aged and young individuals. SONFH may be caused by ischemia and hypoxia of the femoral head, where mitochondria play a crucial role in oxidative reactions. Currently, there is limited literature on whether mitochondria are involved in the progression of SONFH. Here, we aim to identify and validate key potential mitochondrial-related genes in SONFH through bioinformatics analysis. This study aims to provide initial evidence that mitochondria play a role in the progression of SONFH and further elucidate the mechanisms of mitochondria in SONFH.

Methods

The GSE123568 mRNA expression profile dataset includes 10 non-SONFH (non-steroid-induced osteonecrosis of the femoral head) samples and 30 SONFH samples. The GSE74089 mRNA expression profile dataset includes 4 healthy samples and 4 samples with ischemic necrosis of the femoral head. Both datasets were downloaded from the Gene Expression Omnibus (GEO) database. The mitochondrial-related genes are derived from MitoCarta3.0, which includes data for all 1136 human genes with high confidence in mitochondrial localization based on integrated proteomics, computational, and microscopy approaches. By intersecting the GSE123568 and GSE74089 datasets with a set of mitochondrial-related genes, we screened for mitochondrial-related genes involved in SONFH. Subsequently, we used the good Samples Genes method in R language to remove outlier genes and samples in the GSE123568 dataset. We further used WGCNA to construct a scale-free co-expression network and selected the hub gene set with the highest connectivity. We then intersected this gene set with the previously identified mitochondrial-related genes to select the genes with the highest correlation. A total of 7 mitochondrial-related genes were selected. Next, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis on the selected mitochondrial-related genes using R software. Furthermore, we performed protein network analysis on the differentially expressed proteins encoded by the mitochondrial genes using STRING. We used the GSEA software to group the genes within the gene set in the GSE123568 dataset based on their coordinated changes and evaluate their impact on phenotype changes. Subsequently, we grouped the samples based on the 7 selected mitochondrial-related genes using R software and observed the differences in immune cell infiltration between the groups. Finally, we evaluated the prognostic significance of these features in the two datasets, consisting of a total of 48 samples, by integrating disease status and the 7 gene features using the cox method in the survival R package. We performed ROC analysis using the roc function in the pROC package and evaluated the AUC and confidence intervals using the ci function to obtain the final AUC results.

Results

Identification and analysis of 7 intersecting DEGs (differentially expressed genes) were obtained among peripheral blood, cartilage samples, hub genes, and mitochondrial-related genes. These 7 DEGs include FTH1, LACTB, PDK3, RAB5IF, SOD2, and SQOR, all of which are upregulated genes with no intersection in the downregulated gene set. Subsequently, GO and KEGG pathway enrichment analysis revealed that the upregulated DEGs are primarily involved in processes such as oxidative stress, release of cytochrome C from mitochondria, negative regulation of intrinsic apoptotic signaling pathway, cell apoptosis, mitochondrial metabolism, p53 signaling pathway, and NK cell-mediated cytotoxicity. GSEA also revealed enriched pathways associated with hub genes. Finally, the diagnostic value of these key genes for hormone-related ischemic necrosis of the femoral head (SONFH) was confirmed using ROC curves.

Conclusion

BID, FTH1, LACTB, PDK3, RAB5IF, SOD2, and SQOR may serve as potential diagnostic mitochondrial-related biomarkers for SONFH. Additionally, they hold research value in investigating the involvement of mitochondria in the pathogenesis of ischemic necrosis of the femoral head.