Background: Primary ovarian insufficiency (POI) is a common clinical endocrine disorder with a high heterogeneity in both endocrine hormones and etiological phenotypes. However, the etiology of POI remains unclear. Herein, we unraveled the causality of genetically determined metabolites (GDMs) on POI through Mendelian randomization (MR) study with the overarching goal of disclosing underlying mechanisms.

Methods: Genetic links with 486 metabolites were retrieved from GWAS data of 7824 European participants as exposures, while GWAS data concerning POI were utilized as the outcome. Via MR analysis, we selected inverse-variance weighted (IVW) method for primary analysis and several additional MR methods (MR-Egger, weighted median, and MR-PRESSO) for sensitivity analyses. MR-Egger intercept and Cochran's Q statistical analysis were conducted to assess potential heterogeneity and pleiotropy. In addition, genetic variations in the key target metabolite were scrutinized further. We conducted replication, meta-analysis, and linkage disequilibrium score regression (LDSC) to reinforce our findings. The MR Steiger test and reverse MR analysis were utilized to assess the robustness of genetic directionality. Furthermore, to deeply explore causality, we performed colocalization analysis and metabolic pathway analysis.

Results: Via IVW methods, our study identified 33 metabolites that might exert a causal effect on POI development. X-11437 showed a robustly significant relationship with POI in four MR analysis methods (P _IVW=0.0119; P _{weighted-median} =0.0145; P_MR-Egger =0.0499; P_MR-PRESSO =0.0248). Among the identified metabolites, N-acetylalanine emerged as the most significant in the primary MR analysis using IVW method, reinforcing its pivotal status as a serum biomarker indicative of an elevated POI risk with the most notable P-value (P _IVW=0.0007; P_MR-PRESSO =0.0022). Multiple analyses were implemented to further demonstrate the reliability and stability of our deduction of causality. Reverse MR analysis did not provide evidence for the causal effects of POI on 33 metabolites. Colocalization analysis revealed that some causal associations between metabolites and POI might be driven by shared genetic variants.

Conclusion: By incorporating genomics with metabolomics, this study sought to offer a comprehensive analysis in causal impact of serum metabolome phenotypes on risks of POI with implications for underlying mechanisms, disease screening and prevention.