The anti-osteoporosis effects of Vitamin K in postmenopausal women

Identifying potential ferroptosis key genes for diagnosis and treatment of postmenopausal osteoporosis through competitive endogenous RNA network analysis

Objective

Postmenopausal osteoporosis (PMOP) is a common systemic metabolic bone disorder that is owing to the reduced estrogen secretion and imbalance of bone absorption and bone formation in postmenopausal women. Ferroptosis has been identified as a novel modulatory mechanism of osteoporosis. Nevertheless, the particular modulatory mechanism between ferroptosis and PMOP is still unclear. The objective of the current investigation was to detect potential biomarkers connected to ferroptosis in PMOP and discover its probable mechanism through bioinformatics.

Methods

We downloaded PMOP-related microarray datasets from the database of Gene Expression Omnibus (GEO) and obtained the differentially expressed genes (DEGs). Utilizing bioinformatics analysis, the DEGs were intersected with the ferroptosis dataset to obtain ferroptosis-connected mRNAs. Enrichment analysis employing KOBAS 3.0 was conducted to comprehend the biological functions and enrichment pathways of the DEGs. The generation of the protein-protein interaction (PPI) network was conducted with the aim of identifying central genes. Lastly, the coexpression and competitive endogenous RNA (ceRNA) networks were built using Cytoscape. With the help of external datasets GSE56815 to verify the reliability of the hub genes by plotting ROC curves.

Results

We identified 178 DE microRNAs (miRNAs), 138 DE circular RNAs (circRNAs), and 86 ferroptosis-related mRNAs. Enrichment analysis exhibited that mRNAs were primarily connected with the signaling pathways of PI3K/Akt, metabolism, mTOR, FoxO, HIF-1, AMPK, MAPK, ferroptosis, VEGF, and NOD-like receptors. Generation of the PPI network detected eight hub genes. The circRNA/miR-23b-3p/PTEN axis may relieve PMOP by inhibiting ferroptosis through targeting the pathway of PI3K/Akt signaling, which is a vital modulatory pathway for PMOP progression. Moreover, the ROC curves ultimately indicates that the four hub genes have greater diagnostic importance in PMOP samples in contrast to the normal group samples, which may be possible markers for PMOP diagnosis.

Conclusions

Bioinformatics analysis identified four hub genes, namely, PTEN, SIRT1, VEGFA, and KRAS, as potential biomarkers for PMOP diagnosis and management. Moreover, the circRNA/miR-23b-3p/PTEN axis may relieve PMOP by suppressing ferroptosis through targeting the pathway of PI3K/Akt signaling, providing a new avenue to explore the pathogenesis of PMOP.

Identification of ferroptosis-associated biomarkers for the potential diagnosis and treatment of postmenopausal osteoporosis

OBJECTIVE

Postmenopausal osteoporosis (PMOP) is one of the most commonly occurring conditions worldwide and is characterized by estrogen deficiency as well as persistent calcium loss with age. The aim of our study was to identify significant ferroptosis-associated biomarkers for PMOP.

METHODS AND MATERIALS

We obtained our training dataset from the Gene Expression Omnibus (GEO) database using GSE56815 expression profiling data. Meanwhile, we extracted ferroptosis-associated genes for further analysis. Differentially expressed ferroptosis-associated genes (DEFAGs) between OP patients and normal controls were selected using the "limma" package. We established a ferroptosis-associated gene signature using training models, specifically, random forest (RF) and support vector machine (SVM) models. It was further validated in another dataset (GSE56814) which also showed a high AUC: 0.98, indicating high diagnostic value. Using consensus clustering, the OP patient subtypes were identified. A ferroptosis associated gene (FAG)-Scoring scheme was developed by PCA. The important candidate genes associated with OP were also compared between different ferrclusters and geneclusters.

RESULTS

There were significant DEFAGs acquired, of which five (HMOX1, HAMP, LPIN1, MAP3K5, FLT3) were selected for establishing a ferroptosis-associated gene signature. Analyzed from the ROC curve, our established RF model had a higher AUC value than the SVM model (RF model AUC:1.00). Considering these results, the established RF model was chosen to be the most appropriate training model. Later, based on the expression levels of the five DEFAGs, a clinical application nomogram was established. The OP patients were divided into two subtypes (ferrcluster A, B and genecluster A, B, respectively) according to the consensus clustering method based on DEFAGs and differentially expressed genes (DEGs). Ferrcluster B and genecluster B had higher ferroptosis score than ferrcluster A and genecluster A, respectively. The expression of COL1A1 gene was significantly higher in ferrcluster B and gencluster B compared with ferrcluster A and gencluster A, respectively, while there is no statistical difference in term of VDR gene, COL1A2 genes, and PTH gene expressions between ferrcluster A and B, together with gencluster A and B.

CONCLUSIONS

On the basis of five explanatory variables (HMOX1, HAMP, LPIN1, MAP3K5 and FLT3), we developed a diagnostic ferroptosis-associated gene signature and identified two differently categorized OP subtypes that may potentially be applied for the early diagnosis and individualized treatment of PMOP. The ER gene, VDR gene, IL-6 gene, COL1A1 and COL1A2 genes, and PTH gene are important candidate gene of OP, however, more studies are still anticipated to further elucidate the relationship between these genes and ferroptosis in OP.