NMF typing and machine learning algorithm-based exploration of preeclampsia-related mechanisms on ferroptosis signature genes.

Background: Globally, pre-eclampsia (PE) poses a major threat to the health and survival of pregnant women and fetuses, contributing significantly to morbidity and mortality. Recent studies suggest a pathological link between PE and ferroptosis. We aim to utilize non-negative matrix factorization (NMF) clustering and machine learning algorithms to pinpoint disease-specific genes related to the process of ferroptosis in PE and investigate likely underlying biochemistry mechanisms.

Methods: The acquisition of four microarray datasets from the Gene Expression Omnibus (GEO) repository, the integration of these datasets, and the elimination of batch effects formed the core procedure. Genes related to ferroptosis in PE (DE-FRG) were identified. NMF clustering was performed on DE-FRG for unsupervised analysis, generating a heatmap for clustering validation via principal component analysis. Immunocyte infiltration differences between different subtypes were compared to elucidate the impact of ferroptosis on immune infiltration in the placental tissue of PE patients. The application of weighted gene co-expression network analysis (WGCNA) revealed important module genes linked to sample subtypes and disease status. The screening of PE feature genes involved employing SVM, RF, GLM, and XGB machine learning algorithms, and their predictive performance was validated using various analyses and an external dataset. The iRegulon tool was utilized to predict upstream transcription factors associated with ferroptosis feature genes, from which differentially expressed transcription factors were screened to construct a "Transcription Factor-FRG-ferroptosis" regulatory network. Finally, in vitro (cultured cells) and in vivo (rat) models were utilized to evaluate the regulatory mechanisms of ferroptosis in normal and PE placental tissues.

Results: Differential analysis of the four merged GEO datasets identified 41 DE-FRGs. NMF clustering based on DE-FRGs revealed two PE subtypes. Immunocyte infiltration analysis indicated significant differences in immune levels between these subtypes. Further WGCNA analysis identified module genes associated with PE and these two subtypes. Subsequently, we developed an integrated machine learning model incorporating five FRGs and validated its predictive efficacy using various analyses and an external validation dataset. Finally, based on the transcription factor ARID3A and ferroptosis feature genes EPHB3 and PAPPA2, we constructed a "Transcription Factor-FRG-ferroptosis" regulatory network, with in vitro and in vivo experiments confirming that ARID3A promotes the progression of PE and ferroptosis by activating the expression of EPHB3 and PAPPA2.

Conclusion: This analytical journey illuminated a critical regulatory nexus in PE, underscoring the central influence of ARID3A on PE through ferroptosis-mediated pathways.