{"paper_id":"5663eaaa-4cee-4032-af83-dfd3b295d741","body_text":"Endometriosis (EM) is a prevalent condition among women of reproductive age, affecting approximately 190 million women worldwide [ 1 ]. It is characterized by the presence of endometrial tissue (comprising glands and stroma) outside the uterine cavity, leading to symptoms such as dysmenorrhea, abnormal menstruation, dyspareunia, infertility, and more. Affecting approximately 10% of women [ 2 ], EM significantly impacts quality of life. These symptoms disrupt daily activities and reduce work efficiency, with pelvic pain and disease severity accounting for up to 60% of productivity loss.\nEM is known for its estrogen dependence, but its exact underlying causes remain incompletely understood. Currently, laparoscopic surgery is considered the gold standard for diagnosing EM. However, due to its invasive nature, it may not be suitable for early diagnosis or monitoring the progression of the disease. To date, no peripheral blood or endometrial biomarkers have suggested sufficient accuracy in diagnosing EM. Given the chronic nature of EM, long-term pharmacological management is essential. While surgical intervention can provide symptomatic relief and remove existing lesions, it does not offer a permanent solution and carries significant costs. Ahmed M Soliman et al. found that the average annual cost of medical care for patients with endometriosis was more than three times that of non-endometriosis controls ($16,573 vs $4,733;  p <  0.005) [ 3 ]. At the same time, the surgical recurrence rate was high. One systematic evaluation showed a persistence or recurrence rate of 22% at 2 years and 40%–50% at 5 years after surgery [ 4 ]. Therefore, the development and investigation of effective medications for long-term management and symptom control are important aspects of managing EM [ 5 ]. These medications can help alleviate pain and improve overall quality of life for individuals with EM. In summary, further exploration is needed to better understand the pathogenesis of EM, develop improved diagnostic methods, and establish effective disease control strategies. Given the unmet need for accurate diagnosis and effective treatment options for EM, exploring the potential of existing drugs to improve the condition of patients with EM may serve as a viable alternative. Considering the large population of individuals affected by EM, utilizing currently available medications may provide relief and help manage the symptoms associated with the condition.\nRecent clinical research has shown a significant elevation in triglyceride levels among patients with EM compared to healthy women [ 6 ]. Meanwhile, the findings of Hongling Zhang et al. provided evidence consistent with a potential causal relationship between TG and increased risk of EM (OR = 1.112, 95% CI: 1.033–1.198,  p  = 5.03 × 10 ⁻ ³). They concluded that higher serum TG levels were associated with a greater risk of EM [ 7 ]. Additionally, several studies have established a positive correlation between triglyceride levels and the severity of EM [ 8 ]. Triglycerides (TG) play a crucial role in the structural composition of living organisms and actively engage in diverse signaling pathways, exerting varying degrees of influence on cellular functions. TG is implicated in processes such as inflammation, immunity, and steroid hormone metabolism, and their close association with the onset and progression of certain diseases has been well-documented. Nevertheless, the precise nature of the association between TG and EM remains elusive. Consequently, it is imperative to investigate the genetic association between TG and EM, thus fostering a deeper comprehension of the pathogenesis of this condition and subsequently formulating a theoretical basis for its diagnosis and treatment.\nCurrently, the main triglyceride inhibitors encompass fibrates, ω-3 PUFA, and niacin medications.\nFibrates primarily target peroxisome proliferator-activated receptor-α (PPARA), leading to a reduction in triglyceride (TG) levels by 20% to 50% [ 9 ]. Fibrates function as agonists of PPARA, binding to and activating this nuclear receptor. PPARA exhibits predominant expression in tissues actively engaged in lipid metabolism, including the liver, adipose tissue, and skeletal muscle. Upon binding to PPARA, fibrates modulate the expression of genes associated with lipid and lipoprotein metabolism, leading to diverse physiological effects, such as heightened fatty acid oxidation, diminished triglyceride synthesis, and enhanced high-density lipoprotein cholesterol (HDL-C) production [ 10 ]. While fibrates primarily exert their action on PPARA, it is worth noting that PPARA itself possesses the ability to regulate the expression of numerous genes and influence diverse metabolic pathways involved in lipid and glucose metabolism. Therefore, despite fibrates having a singular primary target, their effects can transcend PPARA and indirectly engage multiple targets by means of downstream effects resulting from PPARA activation. For instance, fibrate-induced activation of PPARA upregulates the expression of genes involved in fatty acid transport, oxidation, and lipoprotein metabolism [ 11 ]. Consequently, this can instigate alterations in the activity of enzymes, transporters, and receptors implicated in these processes. Moreover, PPARA activation can modulate inflammatory and oxidative stress pathways, thereby additionally influencing gene expression and cellular processes [ 12 ]. In summary, while fibrates primarily target PPARA, they indirectly affect multiple targets by modulating gene expression and metabolic pathways associated with lipid and lipoprotein metabolism. Therefore, investigating the role of fibrates in EM can enhance our understanding of the relationship between PPARA and multiple metabolic targets, thereby providing a more scientific and rational basis for the use of fibrates.\nOmega-3 polyunsaturated fatty acids (ω-3 PUFA) can lower triglyceride levels in the body through various mechanisms and by targeting multiple pathways. They can reduce TG levels by approximately 25% to 30% [ 13 ]. Regarding triglyceride synthesis, omega-3 fatty acids inhibit diacylglycerol acyltransferase (DGAT), an enzyme involved in the terminal stage of triglyceride synthesis. Through the reduction of DGAT activity, omega-3 fatty acids decrease triglyceride production in the liver [ 14 ]. Simultaneously, omega-3 fatty acids upregulate the activity of enzymes implicated in fatty acid oxidation, such as acetyl-CoA carboxylase 1 (ACC1) and fatty acid synthase (FAS), facilitating the breakdown (oxidation) of fatty acids in the liver. This reduces the availability of glycerol and fatty acids for triglyceride synthesis, consequently lowering triglyceride levels [ 15 ]. Additionally, omega-3 fatty acids can enhance the activity of lipoprotein lipase (LPL), an enzyme that facilitates the breakdown of TG into fatty acids and glycerol for utilization by tissues. This, in turn, increases the rate at which TG are cleared from the bloodstream [ 16 ]. Furthermore, omega-3 fatty acids can promote a favorable distribution of lipids, which includes reducing triglyceride levels, by modulating the expression of genes involved in lipid metabolism. Additionally, omega-3 fatty acids possess anti-inflammatory properties and contribute to cardiovascular health and neuroprotection [ 17 ]. Exploring the multi-target and multi-mechanism characteristics of omega-3 fatty acids may yield novel insights and treatment strategies for addressing the chronic inflammatory state, pain, and hormonal imbalance associated with EM.\nNiacin is currently utilized in cases where fibrates and prescription-grade omega-3 fatty acids have been ineffective in managing TG, albeit its notable adverse effects. While the mechanism of action of niacin remains unclear, niacin appears to mitigate triglyceride and low-density lipoprotein cholesterol levels by modulating multiple targets involved in triglyceride metabolism within adipose tissue and the liver. Research indicates that niacin activates GPR109A, leading to the inhibition of lipolysis (the breakdown of stored fat) in adipose tissue. This, in turn, decreases the release of free fatty acids into the bloodstream, thereby aiding in the reduction of triglyceride synthesis in the liver [ 18 ]. Furthermore, niacin exerts a direct and non-competitive inhibitory effect on diacylglycerol acyltransferase 2 (DGAT2), a crucial enzyme involved in hepatic triglyceride synthesis.\nConsequently, this inhibition leads to decreased triglyceride synthesis and reduced secretion of atherosclerotic lipoproteins from the liver [ 19 ]. Moreover, niacin medications can activate the PI3K/PKB pathway, leading to the inhibition of FAS activity in the liver. As a result, triglyceride synthesis is reduced [ 20 ]. Simultaneously, niacin plays a role in energy metabolism within the body, influencing cellular redox status and ATP production [ 21 ]. Abnormalities in cell metabolism are associated with the development of EM [ 22 ]. Investigating the mechanism of action of niacin in this context can offer valuable insights into the progression of the disease. Furthermore, niacin regulates immune cell function and modulates inflammatory responses [ 23 ], while promoting tissue repair [ 24 ] and exerting antioxidant effects [ 25 ]. EM frequently involves inflammatory responses to tissue damage and heightened oxidative stress. Investigating the mechanisms by which niacin promotes tissue repair and mitigates oxidative stress may contribute to the development of novel treatment strategies for EM.\nIn conclusion, it is evident that a variety of triglyceride inhibitor drugs are currently available. However, selecting a suitable triglyceride inhibitor for EM patients remains a challenge. To provide genetic insights into the potential relevance of triglyceride-related drug targets in EM, our study employed the Mendelian randomization (MR) method to analyze the genetic associations between triglyceride inhibitor-related targets and EM risk. Our analysis aimed to evaluate the genetic evidence supporting the potential relevance of triglyceride-related targets in EM.  Fig 1  provides a detailed overview of the study rationale.\n\nFirstly, we utilized a two-sample univariate MR analysis to estimate the genetically predicted association of TG on EM. This analysis helped evaluate the genetic association between TG and EM and assess whether the findings were compatible with a potential causal effect.\nSecondly, we conducted a search in the drug-gene interaction database to identify target genes associated with triglyceride drugs. By leveraging this database and referring to previous studies on druggable genes, we aimed to characterize specific genes that are targeted by triglyceride drugs and potentially play a role in the development or treatment of EM.\nThirdly, we employed the genetic variants that are mediated by the target genes and associated with TG to serve as proxies for TG-lowering drug targets. Using a two-sample MR approach, we investigated the association between genetic proxies for TG-lowering drug targets and EM risk by analyzing data from two Genome-Wide Association Studies (GWAS).\nFourthly, we conducted colocalization analysis to examine whether the association between TG and EM was driven by the same genomic locus. This analysis could help determine if there was shared genetic regulation underlying both traits.\nFifthly, we explored the protein interaction relationships among the selected genes and examined their enrichment functions. Additionally, we evaluated the significance and importance of the identified genes in the context of EM.\nFinally, we searched for transcription factors (TFs) that regulated these genes, aiming to identify the regulatory mechanisms that influenced the expression and function of the identified genes in relation to EM. The overall study design and analysis procedure are displayed in  Fig 2 . The graphical summary is created by BioGDP [ 26 ].\nBased on the information provided, the data used in the study were obtained from public databases and specific GWAS consortia. The sources of the GWAS data used in the analysis are as follows:\nTG and high-density lipoprotein (HDL) data: The data were obtained from The Global Lipids Genetics Consortium (GLGC). This consortium conducts large-scale genetic studies to unravel the genetic architecture of lipid traits. The GLGC study included 188,577 participants, predominantly of European ancestry, consistent with the majority of its constituent studies. Plasma lipid concentrations were measured after a minimum fasting period of 8 hours. Estimates were adjusted for age, age squared, sex, and population stratification, with exclusion of participants known to be taking lipid-lowering medications. The selected genetic variants collectively explained 10–14% of the total trait variance. For individual variant association estimates, additive genetic models were fitted using inverse-normal transformed trait linear regression. For pooled effect estimates, weighted meta-analysis was performed using Stouffer’s method. In the European-specific analysis (approximately 1.3 million individuals, accounting for 80% of the total sample), the identified lipid-associated signals constituted 76% of all signals. The large sample size ensures that the effect estimates in the European population are relatively stable, and the overall results exhibit high consistency. HDL cholesterol was standardized, with the genetically predicted association reflecting the effect of a 1 standard deviation increase in HDL on the outcome. Triglycerides are log-transformed due to skewed distribution, with the genetically predicted association interpreted as the effect of doubling the triglyceride level [ 27 ].\nVery low-density lipoprotein (VLDL) data: The data were sourced from the IEU OpenGWAS project. The IEU OpenGWAS project is a platform that provides access to a diverse array of GWAS datasets for different traits. VLDL levels were analyzed as a continuous variable, typically reported in mmol/L (standard unit in European GWAS consortia). Sample size: 115,082 participants.\nEM and coronary heart disease (CHD) data: The data for these conditions were obtained from the FinnGen study. FinnGen is a research project in genomics and personalized medicine that aims to understand the genetic basis of diseases in the Finnish population. The FinnGen database includes 377,277 participants. In the FinnGen database, disease diagnoses are primarily derived from the Finnish National Health Registries and electronic health records (EHRs), encompassing inpatient and outpatient diagnoses, medication prescriptions, laboratory test results, and medical procedure documentation. All disease information is standardized using International Classification of Diseases (ICD) codes and cross-validated across multiple data sources to ensure accuracy. For specific conditions such as cancer, diagnoses are further validated through disease-specific registries (e.g., cancer registries) and pathological confirmation. Rare diseases are additionally corroborated using genomic data to enhance diagnostic precision [ 28 ].\nIt is important to note that the use of these data was approved by the respective institutional review boards, ensuring ethical considerations and participant privacy. Informed consent was obtained from all participants involved in the original GWAS studies.\nSee  Table 1  for detailed information on data sources.\nSingle nucleotide polymorphisms (SNPs) that meet the genome-wide significance level  (p  < 5 × 10 −8 ) were chosen as instrumental variables (IVs). Independent genetic variants were identified using a cutoff of linkage disequilibrium (LD) values (threshold set to r2 < 0.001, kb = 10,000) to ensure the independence of the IVs [ 29 ]. Proxy SNPs were not utilized in this study. SNPs lacking necessary statistical information were directly removed.\nThe inverse variance weighting (IVW) method was utilized as the primary statistical approach. In the case of a single SNP, the genetically predicted association was estimated using the Wald ratio test. Weighted median, Simple mode, and Weighted mode methods were employed to assess model robustness. MR-Egger can identify horizontal pleiotropy in models, although it has limited statistical power [ 30 ]. Therefore, in cases where MR-Egger yielded significant results, we investigated the included SNPs using PhenoScanner, a comprehensive genotype-to-phenotype database [ 31 ]. Additionally, to exclude SNPs that may influence EM through other phenotypes, we employed the leave-one-out method to assess whether a single SNP exerted a disproportionate impact on the overall MR estimate [ 32 ]. To mitigate the possibility of reverse causation, we conducted an additional analysis to assess whether EM might influence triglyceride levels.\nThe primary triglyceride-lowering drugs comprise fibrates, ω-3 PUFA, and niacin. Genes encoding the target proteins of these triglyceride drugs were acquired from DGIDB and ChEMBL. All targets gathered from these two sources were verified using NCBI Gene to determine if they were human genes. Non-human genes were excluded. We excluded genes that were not predicted with 90% confidence to be active according to ChEMBL’s Target Predictions. We selected SNPs that exhibited a significant association with TG levels ( p  < 5 × 10 −8 ) within a 100 Mb window surrounding the genomic region of the drug target [ 33 ]. To enhance the instrumental variable’s strength for each drug target gene, a less stringent threshold for independent clustering SNPs was employed (r 2  < 0.3, kb = 100).\nRestricting instrument selection to variants located within or near pharmacologically annotated target-gene regions was intended as a biologically informed strategy to strengthen the link between the selected variants and the predefined target-related pathways. Consistent with current methodological guidance, these regional instruments were treated as proxies for target-related genetic variation rather than as direct equivalents of pharmacological target inhibition or activation [ 34 , 35 ].\nPositive-control MR analyses were used to assess the biological plausibility of the pharmacogenetic instrument framework in a lipid-related disease context. Because triglyceride-lowering medications are used in the management of CHD [ 36 ], and genetic studies have provided evidence consistent with a potential causal relationship between TG and CHD risk [ 37 ], CHD was selected as a positive-control outcome [ 38 ]. Results from this analysis were interpreted only as a supportive consistency check for the target-region instrument framework, rather than as validation of gene-specific pharmacological effects or evidence of therapeutic efficacy for EM [ 34 , 35 ].\nIn clinical studies, fibrates and niacins have been shown to reduce the risk of atherosclerosis by increasing concentrations of HDL-C and decreasing concentrations of TG in plasma [ 39 ,  40 ]. Consequently, HDL-C is utilized as a downstream biomarker of the effects of fibrates and niacins. ω-3 PUFA reduces the synthesis and secretion of VLDL particles and enhances the removal of TG from VLDL and chylomicron particles through the upregulation of LPL [ 41 ]. As a result, VLDL is employed as a downstream biomarker of the effects of ω-3 PUFA.\nWe employed colocalization analysis to investigate whether TG and EM share the same causal variant within the drug target gene region. This analysis aimed to provide additional validation regarding the robustness of including IV. In this study, we assessed the possible probabilities using the following hypotheses: H0 represents the scenario where the gene is not associated with TG or EM; H1 suggests that the gene is associated with TG but not with EM; H2 indicates that the gene is associated with EM but not with TG; H3 suggests that the gene is associated with both TG and EM, but the SNPs are independent of each other; and H4 proposes that the gene is associated with both TG and EM and that they share common SNPs.\nAll analyses were conducted using R 4.2.2 statistical software. The ‘TwoSampleMR’ package was utilized for MR and drug target studies, while the ‘coloc’ package was employed for colocalization analysis.\nThe Gene Expression Omnibus (GEO) database is a publicly accessible repository that collects high-throughput gene expression data from research institutions worldwide. It is maintained by the National Center for Biotechnology Information (NCBI) [ 42 ]. The extensive collection of experimental data in the GEO database is widely utilized in various research fields, including gene regulation, disease mechanisms, and drug discovery.\nOur objective was to investigate the expression patterns and potential regulatory relationships among genes involved in EM using the GEO database. To accomplish this, we downloaded six transcriptome datasets related to EM from the GEO database. The specific datasets we obtained for analysis are  GSE23339  [ 43 ],  GSE25628  [ 44 ],  GSE58178  [ 45 ],  GSE11691  [ 46 ],  GSE7846  [ 47 ], and  GSE7305  [ 48 ]. The details of each data set are shown in  Table 2 .\nIn the statistical analysis conducted using R software, we performed preprocessing steps on the gene expression profiles obtained from the downloaded datasets. The following steps were carried out:\nID conversion: We performed ID conversion on the probes present in the expression profiles, converting them into corresponding gene names. This conversion enabled a more interpretable representation of the data.\nLog2 transformation: To standardize part of the unstandardized dataset, we applied the log2 function. This transformation helped normalize the expression values, ensuring a more symmetric distribution and facilitating subsequent statistical analyses.\nExtraction and merging of expression profiles: Next, we extracted the expression profiles of genes from each dataset separately and merged them. This consolidation allowed for a comprehensive analysis using a combined dataset.\nBatch effect removal: Batch effects were defined as technical differences resulting from the processing and measurement of samples in different batches, unrelated to any biological variation in the experiment. To enhance the accuracy of statistical inference, we used the ‘sva’ package to remove batch effects and other unwanted noise from the merged data. This step helped eliminate potential confounding factors and reduced the impact of batch effects on the results [ 49 ].\nPrincipal Component Analysis (PCA): Prior to and after removing the batch effect, we generated a PCA graph. This analysis helped assess the reduction of batch effects by visually comparing the distribution of samples. The aim was to minimize the influence of batch effects on the results and improve the robustness of subsequent analyses.\nAfter completing these preprocessing steps, we obtained a total gene expression profile comprising 101 samples, including data from 56 EM patients and 45 healthy controls. This processed dataset serves as the basis for further research and analysis.\nTo explore the expression patterns and mutual regulatory relationships of genes in the EM and control groups, we employed various R packages for visualization and analysis. The specific steps are as follows:\nHistograms of gene expression: We utilized the ‘reshape2’ package and the ‘ggpubr’ package to create histograms comparing the gene expression levels between EM patients and healthy controls. This visualization helped identify any differences in gene expression patterns between the two groups.\nSpearman correlation test: To assess the correlation of gene expression levels within the EM patient and healthy control groups separately, we employed the ‘corrplot’ package and performed Spearman correlation tests [ 50 ]. This analysis allowed us to evaluate the strength and direction of the relationships between gene expression levels. The results of these tests were visualized using the ‘ggplot2’ package.\nVisualization of gene positions and genome structure: We utilized the ‘RCircos’ package to visualize the relationship between the positions of the eight genes on the chromosome and the genome structure. This visualization approach provides insights into the spatial organization and arrangement of these genes within the genome [ 51 ].\nBy implementing these steps, we can gain a better understanding of the expression patterns, correlations, and genomic relationships of the genes involved in EM and control groups. The visualizations generated through these analyses facilitate the interpretation and communication of the results.\nTo explore the direct correlation between genes, the drug targets identified through MR analysis were input into the STRING database to establish a protein-protein interaction (PPI) network. This network helped uncover the interactions among the genes of interest.\nIn addition, a PPI network was also constructed in GeneMANIA to further examine the genes with significant indirect correlations and the functional roles they play. GeneMANIA provided insights into gene functions and interactions. The complementary nature of STRING and GeneMANIA has been suggested in several studies, with the former providing direct interactions and the latter extending indirect associations to improve the comprehensiveness and reliability of network-based functional annotation. STRING was used to initially construct the PPI network and identify core genes. GeneMANIA was further validated through enrichment analysis of the functional modules, which was cross-checked with the STRING results to reduce false positives [ 52 ].\nTo gain more insights from the constructed PPI networks, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment pathway analyses were performed. The thresholds for significance were set as false discovery rate (FDR) < 0.05 and strength > 0.01. These analyses helped identify the biological processes and pathways associated with the genes in the network. The threshold of > 0.01 for interaction strength in STRING was selected based on the database’s recommendations and established practices in network analysis [ 53 ]. In addition, a threshold of strength >0.01 preserves weak but potentially biologically significant associations while avoiding excessive stringency leading to the omission of potential pathways. In our study, this threshold allowed us to broadly capture potential interactions among lipid metabolism and inflammation-related genes, which were further validated through functional enrichment and sensitivity analyses [ 54 ]. To evaluate the importance of the target genes in the network, maximal clique centrality (MCC) analysis was conducted using the CytoHubba plugin in Cytoscape software. This analysis assessed the centrality of the target genes, indicating their significance within the network.\nFurthermore, the iRegulon plugin v1.3 [ 55 ] in Cytoscape software was employed to predict potential TFs associated with the key genes in the network. The enrichment score threshold was set to 3.0, area under the curve (AUC) threshold for the receiver operating characteristic (ROC) analysis was set to 0.03, ranking threshold was set to 5000, and the maximum FDR on motif similarity was set to 0.001. This analysis helped identify TFs that may regulate the expression of the key genes.\nBy performing these analyses and utilizing the mentioned tools and thresholds, this study aimed to gain insights into the gene interactions, functions, regulatory factors, and pathways associated with the identified drug targets.\n\nBased on the analysis conducted, a total of 350 genome-wide significant SNPs associated with TG levels and 26 SNPs associated with EM were included as instrumental variables (IVs). The inverse-variance weighted (IVW) estimation was employed to assess the association between genetically predicted TG levels and EM risk.\nThe results of the IVW estimation indicated that genetically predicted TG levels were significantly associated with EM risk (OR = 1.1853, 95% CI: 1.0704–1.3125,  p  = 0.0011). This suggests that higher genetically predicted TG levels are associated with an increased risk of EM ( Fig 3A ).\n(A) Genetically predicted associations of TG on EM; (B) Genetically predicted associations of EM on TG. The Fig shows the results of two data sets, A and B, involving different exposures and outcomes, demonstrating a comparison of several statistical methods. For group A, the exposure was TG and the outcome was EM. OR and  p  values show the effect estimates and statistical significance, respectively, for the different methods. ORs and their 95% confidence intervals (CIs) are represented as dots and lines in the Fig. For group B, the exposure was EM and the outcome was TG. Beta and  p  values show effect estimates and statistical significance, respectively. Nsnps represents the number of SNPs used for each method.\nOn the other hand, the analysis did not find a significant association between EM and TG levels (β = 0.0145, 95% CI: −0.0070 to 0.0361,  p  = 0.1862). This indicates that EM was not significantly associated with TG levels ( Fig 3B ).\nSince the MR-Egger analysis yielded significant results ( p <  0.05) in the TG-EM analysis, the leave-one-out method was utilized to test the significance of individual SNPs on the outcome. The results suggested that by eliminating SNPs related to TG-associated SNPs one by one and recalculating the effects of the remaining SNPs, the overall results remained consistent, suggesting that the pooled effects were not affected ( S1 Fig ).\nThese findings support the hypothesis that genetically predicted higher TG levels are associated with an increased risk of EM, while EM itself does not significantly impact TG levels.\nWe analyzed the target genes of 35 fibrates, 12 niacins, and 10 ω-3 PUFA drugs included in DGIDB and ChEMBL.\nOut of these, 21 genes showed genetic associations with triglyceride concentration, and among these TG-associated genes, 7 genes showed genetic associations with EM in the present analyses.\nCholine Fenofibrate and Ciprofibrate target TNF and SLC9A1 significantly among fibrates, while Fenofibrate targets  APOA1 ,  VEGFA , and  SCN3A . Among niacin drugs, Probucol significantly targets  ITGAV  and  SCN3A .  SCN3A  is a significant target for all three ω-3 PUFA drugs.  Table 3  displays the information about the included drug targets. Detailed explanations for gene exclusion or inclusion are provided in  S1 Table .\nAll significant targets associated with TG showed positive associations with CHD risk (OR > 1,  p  < 0.05) ( Fig 4 ).\nExposure: Indicates genes; nsnp: Number of single nucleotide polymorphisms (SNPs) associated with the gene in the sample;  p  value: Significance level, a value less than 0.05 is considered significant;  F  statistic: The  F  statistic, used to assess instrument strength. OR (95% CI): Odds Ratio and its 95% confidence interval. An OR greater than 1 indicates increased risk, while less than 1 indicates reduced risk.\nIn the TG-based drug-target MR analysis, genetically proxied effects related to  SCN3A, VEGFA, APOA1, TNF,  and  BACE1  were associated with increased EM risk (OR > 1,  p < 0.05 ), whereas genetically proxied effects related to  SLC9A1  and  ITGAV  were associated with reduced EM risk (OR < 1 , p < 0.01).  In sensitivity analysis,  SCN3A  and  BACE1 , which utilized very low density lipoprotein (VLDL) as downstream markers, yielded similar results in TG (OR < 1,  p  < 0.01). However, all included genes exhibited the opposite effect in HDL-C compared to TG and VLDL (OR > 1,  p  < 0.01) ( Fig 5 ).\nTarget genes: The genes being analyzed for their association with triglyceride levels. Exposure: The type of lipid measurements includes TG for triglycerides and HDL-C as high-density lipoprotein cholesterol. nsnp: Number of single nucleotide polymorphisms (SNPs) associated with the respective gene.  p  value: Significance level; values below 0.05 were considered statistically significant.  F  statistic: The  F  statistic, used to assess instrument strength. OR (95% CI): Odds Ratio and its 95% confidence interval. An OR greater than 1 indicates increased risk, while less than 1 indicates decreased risk. If the OR is less than 1 and the confidence interval does not include 1, this suggests an association with decreased risk. Narrower confidence interval indicates a more precise estimate.\nThe PP.H4 values for all seven genes were below 0.8, indicating limited evidence that TG and EM share a common causal variant within the examined target-gene regions. In contrast, the PP.H3 value for TNF exceeded 0.8, suggesting that TNF may be associated with both traits, but through independent genetic variants. For  APOA1 ,  VEGFA ,  SLC9A1 , and  BACE1 , PP.H1 values were higher than the other hypotheses, suggesting stronger evidence for association with TG than with EM in these regions. Overall, the colocalization results did not provide strong support for shared causal variants between TG and EM at the examined drug-target loci ( Table 4 ).\nGiven that fibrates serve as ligands for PPARA, we also examined the association between PPARA and the targets under investigation.\nThe batch effects were removed from the data. For details, see  S2 Fig  A and B. Compared with the healthy control group, the EM group exhibited upregulated expression of  SCN3A  ( p  < 0.001),  BACE1 , and  APOA1  ( p  < 0.0001), whereas no statistically significant differences were observed in the expression levels of  SLC9A1 ,  TNF ,  PPARA ,  VEGFA , and  ITGAV  ( p  > 0.05) ( Fig 6C ).\n(A) Gene-expression correlation matrix of the eight selected genes in EM patients. (B) Gene-expression correlation matrix of the eight selected genes in healthy controls. (C) Box plots comparing the expression levels of the eight selected genes between EM patients and healthy controls. Red indicates a positive correlation, whereas blue indicates a negative correlation. “ns” indicates no statistical significance. Significance:  p  < 0.05 (*),  p <  0.01 (**),  p <  0.001 (***), and  p  < 0.0001 (****).\nGene-expression correlation analysis revealed significant negative correlations between  SCN3A  and  PPARA  (correlation coefficient = −0.31,  p  < 0.05), as well as between  SCN3A  and  VEGFA  (correlation coefficient = −0.46,  p  < 0.01), within the healthy control group ( Fig 6B ).\nIn the EM patient group, positive correlations were observed between  SLC9A1  and  VEGFA  (correlation coefficient = 0.42,  p  < 0.01), as well as between  SCN3A  and  ITGAV  (correlation coefficient = 0.27,  p  < 0.05). Conversely, negative correlations were found between  SLC9A1  and  ITGAV  (correlation coefficient = −0.36,  p  < 0.01),  BACE1  and  TNF  (correlation coefficient = −0.33,  p  < 0.05),  APOA1  and  TNF  (correlation coefficient = −0.27,  p  < 0.05),  APOA1  and  VEGFA  (correlation coefficient = −0.43,  p  < 0.01), and  PPARA  and  VEGFA  (correlation coefficient = −0.34,  p  < 0.05) ( Fig 6A ). The chromosomal locations of the eight genes are shown in  S2 Fig  C.\nAfter including  PPARA , interactions were observed among seven out of the eight genes.  VEGFA, TNF , and  PPARA  were co-mentioned textual associations (i.e., they were jointly mentioned in the PubMed abstract) as well as co-expression interactions between pairs.  PPARA  and  APOA1  were textually associated, co-expressed, and experimentally confirmed to interact. A textual interaction relationship was indicated between  APOA1  and  TN F. While  TNF  and  BACE1  exhibited both textual and co-expression interaction relationships.  BACE1  and  SLC9A1  showed both textual and co-expression interaction relationships, as shown in  Fig 7A .\n(A) PPI network created with STRING. (B) MCC node-importance analysis. In  Fig 7A , background colors denote specific biological processes, and different edge styles represent different types of interactions. In  Fig 7B , the color coding indicates node importance.\nThe GO enrichment analysis results for the eight genes identified the top five enriched pathways with the smallest FDR values, listed in ascending order, as follows: Regulation of localization, Regulation of transport, Regulation of vesicle-mediated transport, Regulation of biological quality, and Response to external stimulus. All genes were associated with the Regulation of localization pathway. While seven genes ( Tnf, Slc9a1, Apoa1, Ppara, Vegfa, Itgav, Scn3a ) were enriched in the regulation of transport pathway, as depicted in  Fig 7A . Gene symbols in the STRING-derived network are shown in the original output format from the STRING database.\nSince STRING-based PPI network analysis showed no direct interaction between SCN3A and the other six genes, we utilized the GeneMANIA website to explore indirect relationships between  SCN3A  and the other six genes. Additionally, we employed GeneMANIA to predict potential associated genes that exhibit significant associations with these seven genes.\nIn the regional network,  PPARA ,  TNF , and  APOA1  exhibited clear physical interaction and pathway relationships. These three genes regulated lipid transport, while  PPARA, APOA1 , and  ITGAV  were involved in lipid transport.  APOA1, VEGFA, ITGAV , and  SLC9A1  collectively mediated cell-substrate adhesion,  TNF  and  VEGFA  collectively participated in the regulation of vasculature development.  BACE1  exhibited Pathway connections with  KLC1  to  KLC4 , whereas  SCN3A  suggested Genetic Interactions connections with  KLC4  and  SLC9A1  ( Fig 8 ).\nEach node represents a gene. Lines connecting the nodes signify different types of interactions among the genes. The colors and styles of the edges indicate the nature of these interactions. Larger nodes indicate greater importance and more interactions.\nThe main regulatory factors for the analyzed genes were identified using the Cytoscape plug-in iRegulon, showing that PPARA and estrogen receptor 1 (ESR1) exhibited the highest normalized enrichment score (NES) values among the transcription factors (TFs). PPARA functions as a transcription factor for 5 genes, including itself, while ESR1 acted as a transcription factor for the 7 genes (excluding  TNF) . The number of motifs associated with PPARA and ESR1 significantly exceeded that of the transcription factor ranked third in terms of NES ( Fig 9 ).\n(A) Transcription factor regulatory network. Purple circles indicate target genes regulated by TFs, and green octagons indicate TFs. (B) Transcription factor normalized enrichment score (NES).\nFig 9A  shows the transcription factor regulatory network. Purple circles show target genes regulated by transcription factors. Green octagons indicate transcription factors. Connecting lines indicate regulatory relationships between transcription factors and target genes. The direction of the line indicates the direction of regulation.  Fig 9B  shows the normalized enrichment score (NES) of transcription factors. The table lists the NES values, the number of target genes and the number of motifs associated with different transcription factors.\n\nThis MR study investigated genetic associations between lipid-lowering drug targets and EM risk. The findings suggest that several TG-related drug targets, including targets related to fibrates, niacin, and ω-3 PUFA, may show genetic associations with EM risk. Among these, fibrate-related targets showed relatively broader target coverage in the present analyses. In particular, fenofibrate-related targets were linked to multiple EM-associated signals and showed directionally similar associations in CHD and EM analyses. These findings suggest that fibrate-related pathways may warrant further investigation in EM research.\nBoth reverse MR and colocalization analysis did not support a bidirectional causal or shared genetic basis between TG levels, TG-lowering drug target genes, and EM. These findings suggest that selected TG-related drug targets may have potential relevance to EM-related biology through TG-associated pathways, although the evidence does not establish a definitive TG-mediated mechanism. This further supports the possibility that TG-related pathways may be linked to EM through multiple targets. The seven identified genes were not selected based on prior assumptions but were derived through a stepwise genetic screening framework integrating TG-related drug targets, genetic associations with TG, and EM-related MR analyses.\nImportantly, the colocalization analysis did not provide strong evidence that TG and EM share the same causal variants within the examined drug-target gene regions, as all PP.H4 values were below the commonly used threshold of 0.8. This finding weakens the certainty that the observed drug-target MR associations operate exclusively through shared TG-mediated genetic mechanisms. Therefore, our results should be interpreted as hypothesis-generating genetic evidence suggesting that TG-related pathways and selected lipid-lowering drug targets may be relevant to EM biology, rather than definitive proof of a direct TG-mediated mechanism [ 56 ]. This cautious interpretation is also consistent with the broader caution required when interpreting findings from recent endometriosis-related genetic research [ 57 ]. Further experimental and clinical studies are needed to validate these candidate targets and clarify their mechanistic roles in EM.\nThe transcriptome analysis of eight genes revealed that SCN3A, BACE1, and APOA1 exhibited significant up-regulation in EM samples. Moreover, considering the MR results, these three genes showed positive MR-estimated associations with TG, indicated by an OR greater than one. Therefore, in the context of EM, the upregulation of  SCN3A ,  BACE1 , and  APOA1  may reflect TG-related biological dysregulation relevant to EM, although the direction and mechanism of these relationships require experimental validation. These findings suggest that  SCN3A, BACE1,  and  APOA1  may be involved in TG-related biological pathways relevant to EM. However, whether modulation of these genes can alter TG levels or EM-related phenotypes requires further experimental validation. The following is a separate explanation of each target gene.\nOur study revealed that TNF exhibits the highest importance in the protein interaction network concerning the impact of TG on EM targets. As a cytokine primarily synthesized by activated macrophages and T cells, TNF exerts a pivotal role in the immune response, inflammatory response, and cellular apoptosis. The previous studies observed elevated TNF levels in the peritoneal fluid of women diagnosed with EM, and a positive correlation between TNF levels and the size of endometriotic lesions was identified [ 58 ], demonstrating the involvement of TNF in the pathogenesis and progression of EM. As an inflammatory factor, TNF initiated a cascade of inflammatory pathways that can influence the synthesis, breakdown, and transport of lipids [ 59 ]. This notion was further supported by the enrichment pathway analysis conducted in our study.\nIn the context of the inflammatory response, TNF, acting as a pro-inflammatory factor, initiated and intensified the inflammatory response to infections and trauma by attracting additional immune cells and pro-inflammatory mediator signals to the site of injury, thereby exacerbating the development of EM. Furthermore, TNF and interleukin-1 (IL-1) stimulated the expression of cyclooxygenase-2 (COX-2), which, in turn, regulates the synthesis of prostaglandin E2 (PGE2). COX-2 was typically absent under normal conditions and was only upregulated during instances of inflammation and injury [ 60 ]. The upregulation of COX-2 expression augments the synthesis of PGE2, which, in turn, can induce the expression of COX-2 and generate surplus PGE2 through a positive feedback loop. The surplus of PGE2 exacerbated inflammatory responses and pain, thereby fostering the progression of EM [ 61 ].\nAdditionally, the enrichment pathway analysis conducted in this study revealed the significant involvement of TNF in angiogenesis. TNF activated the vascular switch in vascular endothelial cells (ECs), thereby facilitating the process of angiogenesis [ 62 ]. The formation of new blood vessels facilitated the establishment of endometriotic implants in ectopic locations [ 63 ], thereby promoting the occurrence of EM.\nThe enriched pathway analysis of the GeneMANIA gene interaction network conducted in this study identified vascular endothelial growth factor A (VEGFA) as another significant contributor to angiogenesis, alongside  TNF .\nVEGFA is recognized as one of the most potent and indispensable angiogenic factors. It stimulates angiogenesis by stimulating the proliferation and migration of vascular endothelial cells (ECs), as well as enhancing vascular permeability [ 64 ]. The establishment of an effective blood supply is crucial for the formation, survival, and growth of ectopic endometrial implants, and angiogenesis plays a vital role in this process [ 65 ]. Thus, VEGFA may contribute to the development of EM through angiogenesis-related mechanisms. Furthermore, the peritoneal fluid of patients with EM exhibits an elevated presence of activated macrophages [ 66 ], which serve as the primary source of VEGFA [ 67 ]. This mechanism increases circulating VEGFA concentrations, amplifying its angiogenic effects. Additionally, several studies demonstrate that VEGFA can recruit macrophages [ 68 ], thereby establishing a positive feedback loop that results in the substantial accumulation of macrophages. These macrophages can secrete pro-inflammatory cytokines (TNF, IL-1β, and IL-6) [ 69 ], further exacerbating the inflammatory environment within the uterine membrane and intensifying the severity of EM [ 70 ].\nThe apolipoprotein A1 gene ( APOA1 ) encodes the major component of HDL in plasma. HDL, in turn, delivers APOC2, which activates low-density lipoprotein (LDL) to facilitate TG breakdown [ 71 ]. Consequently,  APOA1  is commonly regarded as a protective factor against CHD. Our research indicates that the SNPs within  APOA1  exhibit variations between TG and HDL, resulting in divergent outcomes. Consequently, the protective influence of the  APOA1  gene on CHD primarily stems from its role within HDL. Conversely, within TG, the  APOA1  genetic proxy showed a risk-increasing association with CHD, a pattern that is also observed in EM. In fact, APOA1 serves as a constituent not only of HDL but also of chylomicrons (CM). Given that 85% to 92% of CM consist of TG by mass [ 72 ], it is incorrect to employ APOA1 as a sole representation of HDL.\nTranscriptome analysis revealed a significant upregulation of  APOA1  in patients with EM, which aligns with the findings reported by JRA Sherwin. According to JRA Sherwin’s research,  APOA1, SFRP4,  and  PAPPA  are typically downregulated by human chorionic gonadotropin (hCG), but in animals with EM, they exhibit upregulation upon hCG stimulation [ 73 ]. Similarly, our study found that estrogen receptor 1 (ESR1) functions as a transcription factor for six target genes, including  APOA1 . Given that hCG can suppress the expression of ESR1 [ 74 ], the upregulation of  APOA1  in EM is associated with the atypical hormonal effects observed in EM.\nThe gene  ITGAV , which plays a role in the APOA1-mediated signaling pathway [ 75 ], is also a target of our study.  ITGAV  encodes the integrin αV subunit, which is a class of membrane protein receptors involved in cell adhesion to the extracellular matrix. Additionally,  ITGAV  has shown to promote adipocyte proliferation [ 76 ]. Notably, while high body mass index (BMI) represents a risk factor for coronary heart disease (CHD) [ 77 ], it is negatively correlated with the risk of EM [ 78 ]. This observation aligns with the findings regarding the performance of ITGAV in TG in CHD and EM, as suggested in our study. EM patients often experience lipodysfunction and fat loss [ 79 ]. The dysfunction and reduction of adipocytes, which are the primary cells responsible for storing TG, may lead to an excessive release of TG into the bloodstream. Therefore, the role of ITGAV in both TG and EM can partially explain the phenomenon of elevated TG levels and low BMI in EM patients. This finding also serves as a reminder that in the treatment of EM, it is crucial to reduce triglyceride levels in the blood without compromising the natural formation of adipocytes.\nThe  SCN3A  gene encodes Nav1.3, a voltage-gated sodium ion channel protein with selective permeability to sodium ions [ 80 ]. Previous studies have suggested that following spinal cord injury, there is abnormal expression of Nav1.3 Na (+) channels in second-order spinal dorsal horn neurons and third-order thalamic neurons along the pain pathway. This change leads to the overexcitability of these neurons, functioning as amplifiers and generators of pain [ 81 ]. In the context of EM, it is possible that Nav1.3 channels also contribute to the amplification of pain signals through a similar mechanism. The current study reveals that the expression of  SCN3A  in EM patients was significantly elevated compared to healthy individuals. These findings suggest that SCN3A-related pathways may be involved in pain-related mechanisms in EM; however, whether modulation of SCN3A can influence EM-associated pain requires further experimental investigation.\nThe sodium-hydrogen exchanger-1 (NHE-1), encoded by the SLC9A1 gene, is a plasma membrane protein involved in maintaining intracellular pH and volume [ 82 , 83 ]. NHE-1, along with the sodium-bicarbonate cotransporter NBCn1 (encoded by the  SLC4A7  gene), facilitates the efflux of acid from vascular smooth muscle cells (VSMCs) and ECs in arteries. This acid efflux is essential for maintaining intracellular pH homeostasis [ 84 ]. Damage to NHE-1 function disrupts acid excretion, leading to the accumulation of intracellular acid. This accumulation can directly affect the receptors of nociceptive afferent fibers, increasing their sensitivity and resulting in heightened pain perception [ 85 ]. Therefore, NHE-1 may represent a potential biological pathway related to EM-associated pain, although its therapeutic relevance requires further validation.\nIt is important to note that while this study supports the potential of NHE-1 in alleviating EM, activating NHE-1 can promote the proliferation of VSMCs, potentially leading to myocardial hypertrophy [ 86 ]. This is particularly concerning for patients with comorbid CHD. Therefore, caution should be exercised when considering drugs that act on NHE-1 in patients with both EM and CHD.\nβ-site amyloid precursor protein cleaving enzyme 1 (BACE1), also known as β-secretase, is a type 1 membrane protein involved in the cleavage of β-amyloid precursor protein in conjunction with γ-secretase. Its primary function is to produce beta-amyloid. BACE1 has been extensively studied in the context of Alzheimer’s disease (AD). While its expression is highest in the brain, it is also found in various other tissues, including endocrine tissue, pancreas, muscle tissue, respiratory tissue, bone marrow, and lymphoid tissue [ 87 ]. In this study, a gene interaction analysis using the GeneMANIA protein interaction network highlighted a connection between TNF and BACE1 through the gene RFFL. RFFL suppresses TNF-induced IKK/NF-κB activation by degrading internalized RIP in endocytic vesicles [ 88 ]. NF-κB, in turn, can increase the transcriptional activity of BACE1. Additionally, research has shown that TNF upregulates  BACE1  expression by activating PKC signaling and sequentially cleaving ST6 β-galactoside α-2,6-sialyltransferase (ST6Gal-I) [ 89 ], providing an explanation for the observed interconnection between BACE1 and TNF in the PPI network analysis.\nRegarding angiogenesis, existing studies have suggested that BACE1 can promote blood vessel formation [ 90 ]. In the case of EM, the growth and invasion of endometriotic implants into ectopic sites within the body require the development of new blood vessels [ 63 ]. This suggests that BACE1 may be involved in angiogenesis-related processes relevant to EM, which is consistent with the genetic association observed in this study.\nBoth AD and EM have a higher prevalence in women [ 91 ]. In the case of AD,  BACE1  has been a major focus of research as a potential therapeutic target. Another important target associated with AD is  ESR1  [ 92 , 93 ], which is the gene encoding the estrogen receptor alpha. This study suggests that ESR1 serves as a transcription factor for  BACE1 , indicating a regulatory relationship between the two genes. Dysfunction of ESR1 can lead to neuroinflammation, which in turn can result in abnormal gene expression, including the upregulation of  BACE1.  It is important to note that the damaging effects of ESR1 dysfunction are not limited to the brain but can also have an impact on other tissues, including the endometrium, which is geographically distant from the brain. The study implies that the dysregulation of ESR1 in endometrial tissue may contribute to neuroinflammatory processes [ 93 ] and abnormal  BACE1  expression, potentially linking the pathogenesis of AD and EM. This suggests that there may be shared molecular mechanisms and pathways between these two conditions, highlighting the importance of considering systemic effects and potential connections between diseases affecting different parts of the body.\nPeroxisome proliferator-activated receptors (PPARA, PPAR-β/δ, and PPAR-γ) belong to the nuclear receptor superfamily and are widely recognized for their significant involvement in lipid metabolism. Animal studies have suggested the presence of protein expression of the  PPAR  family, including  PPARA , in pig endometrial luminal epithelial cells and stromal cells [ 94 ]. Despite being an important target of fibrate drugs, this study found no direct evidence linking it to endometriosis (EM) pathogenesis. However,  APOA1, TNF , and  VEGFA , which directly interact with  PPARA , are targets of EM and are affected by fibrate drugs. Simultaneously, PPARA can function as a transcription factor for  SLC9A1, VEGFA, ITGAV , and  APOA1 . Studies have confirmed that PPARA activation can concomitantly decrease the expression of TNF and VEGFA [ 95 ]. As a transcription factor for APOA1, Gervois et al.‘s study suggested that PPARA activation induces the expression of APOA1 and APOA2 in the human liver, leading to an increase in plasma HDL-C and consequent reduction in plasma triglyceride levels [ 96 ], thereby potentially influencing EM-related biological processes through modulation of multiple downstream targets.\nPrevious studies have primarily identified the association between TG and EM through clinical investigations. For example, clinical studies have shown that patients with endometriosis (EM) have significantly abnormal lipid profiles, as evidenced by a 38% increase in low-density lipoprotein (LDL), a 30% increase in total cholesterol (TC), and a 27% increase in triglycerides (TG) compared to controls [ 97 ]. Patients with endometriosis have an increased risk of hyperlipidemia (OR = 1.25, 95% CI: 1.03–1.52) and a significantly higher risk of hypertension (OR = 1.29, 95% CI:1.07–1.56) [ 98 , 99 ]. The risk of microvascular dysfunction and atherosclerosis was significantly increased in patients with EM due to chronic inflammation and oxidative stress, which together with early menopause (OR = 1.32, 95% CI:1.07–1.64) exacerbated cardiovascular risk [ 100 ]. To further comprehend the genetic pathway and the role of TG in EM, the exploration was extended to include the interaction network between TG inhibitor targets. By examining the interactions and relationships among these targets, this study aims to elucidate the molecular mechanisms underlying the influence of TG on EM.\nThis study represents the first MR investigation to explore the potential association between lipid-lowering drug-targeted genes and EM risk, offering new insights into the relationship between lipid metabolism and EM. The findings suggest that several TG-related drug targets, including targets related to fibrates, niacin, and ω-3 PUFA, may show genetic associations with EM risk. Among these, fibrate-related targets showed relatively broader target coverage in the present analyses. Furthermore, transcriptome analyses identified  SCN3A ,  BACE1 , and  APOA1  as upregulated genes in EM, providing additional evidence that lipid-related pathways may be relevant to EM biology. However, these findings should be interpreted as hypothesis-generating and require further experimental and clinical validation.\nHowever, this study had several limitations. First, despite the elimination of numerous single nucleotide polymorphisms (SNPs) that might influence the results via phenotypes unrelated to triglyceride (TG) levels, the possibility of horizontal pleiotropy could still introduce bias to the MR findings. Nevertheless, our leave-one-out sensitivity analysis suggested that the association between genetically predicted TG levels and EM risk was not driven by any single SNP. Second, our MR analysis relies on existing genome-wide association study (GWAS) data, which only reflected the lifelong genetic exposure effect. In practical terms, the available data failed to distinguish TG levels at different life stages or within specific time windows. As a result, we were unable to perform sensitivity analyses or subgroup comparisons “across other exposure periods/timings.” This data limitation may have prevented us from uncovering potential variations in effect size related to exposure timing. Third, the study’s focus on the European population limited the generalizability of the findings to other ethnic groups, necessitating further investigation. In terms of genetic association, due to the absence of the VEGFA gene in STRING’s Homo sapiens database, we utilized STRING’s Mus musculus database to construct the PPI network, highlighting the need for further investigation into the applicability of this PPI network in humans. To compensate for this limitation, the study additionally conducted gene association predictions using GeneMANIA for Homo sapiens, providing readers with a point of reference. Fourth, this study was based solely on data from existing databases and lacked experimental validation; therefore, the findings should be interpreted with caution until confirmed by further experimental studies. Nevertheless, our results may offer valuable insights and directions for future experimental research.\nAlthough MR analysis can reduce confounding and reverse causation compared with conventional observational studies, its interpretation relies on several key assumptions, including the validity, independence, and exclusion restriction of genetic instruments [ 30 ]. Although multiple sensitivity analyses were performed in this study, residual horizontal pleiotropy or other unmeasured genetic effects cannot be completely excluded [ 30 ]. Therefore, the observed genetic associations should be interpreted cautiously and require further experimental and clinical validation. Additional caution is required when interpreting genetic variants located within or near pharmacologically annotated target genes. Such regional variants represent lifelong genetically influenced differences and cannot be regarded as direct equivalents of drug administration, target inhibition or activation, treatment dose, treatment duration, tissue-specific drug action, or off-target effects. Moreover, variants within a target-gene region may influence neighboring genes or regulatory elements, limiting gene-specific pharmacological interpretation. Accordingly, the observed target-region associations should be regarded as exploratory and hypothesis-generating evidence rather than direct estimates of the therapeutic effects of the corresponding drugs [ 34 , 35 , 101 ].\nThe transcriptome analyses in this study were based on publicly available datasets generated from different studies, which may involve variations in sample characteristics, experimental platforms, study designs, and experimental conditions. Although batch-effect correction was performed during data processing, residual heterogeneity cannot be completely excluded. Therefore, these transcriptomic findings should be further validated in independent cohorts and experimental studies.\nOur study provides MR-based evidence suggesting that TG-related pathways and lipid-lowering drug targets may be relevant to EM biology. These findings do not establish triglyceride inhibitors as therapeutic agents for EM, but they highlight fibrate-related pathways as candidates for further mechanistic and translational investigation. Future experimental and clinical studies are needed to determine whether modulation of TG-related pathways can influence EM-related phenotypes.\nHowever, the results also highlight several avenues for future research. Notably, the GWAS data employ in our analysis only provide an estimate of the lifelong exposure effect, without differentiation of TG levels across specific life stages or time windows. Consequently, it was not possible to assess how the associations might differ “across other exposure periods/timings.” Future studies capturing detailed, time-specific exposure information may further elucidate whether the relationship between TG levels and EM risk varies with different exposure periods. Enhancing the granularity of exposure data could contribute to developing more targeted and timely clinical interventions.\nIn our study, we used MR + multi-omics + drug repurposing strategy to investigate the role of lipid-lowering drugs in endometriosis. MR utilizes genetic variation as an instrumental variable, which can effectively avoid confounding bias and reverse causality in traditional observational studies. Transcriptome analysis, protein interaction networks and pathway enrichment can be used to reveal the “gene-pathway-phenotype” networks of other diseases. Target analysis lays the foundation for repurposing existing drugs and greatly shortens the development cycle of new drugs. Through the integration of genetic association analysis, multi-omics evidence, and drug-target annotation, this study provides a hypothesis-generating framework for investigating lipid-related pathways in complex diseases. Our study may provide a framework for integrating genetic association analysis, multi-omics evidence, and drug-target annotation to generate hypotheses for future translational research. What is more, in the future, we will incorporate electronic prediction of drug-gene interactions (e.g., molecular docking, pathway modeling) into prediction systems. This may improve the predictive value and biological interpretability of future drug-target prioritization studies.\n\nThis study presented MR evidence suggesting that TG-related pathways and triglyceride-lowering drug targets may be relevant to EM biology. Three classes of compounds were identified: fibrates, nicotinic acid, and ω-3 PUFA. Seven candidate drug-target-related genes were implicated:  TNF ,  SLC9A1, APOA1, VEGFA, ITGAV, SCN3A,  and  BACE1 . These selected targets formed an interaction network associated with biological processes potentially relevant to EM, including lipid transport, localization, cell adhesion, and blood vessel development. PPARA and ESR1 were identified as the two TFs with the highest enrichment scores in the regulatory analysis. Considering the genetic evidence and the observed involvement of fibrate-related targets, fenofibrate-related pathways may warrant further mechanistic and translational investigation in EM.\n\n(TIF)\n(A) Analysis of the PCA plot prior to batch effect removal revealed distinct clustering of samples from the same dataset and clear stratification between batches. These findings indicate the presence of substantial batch effects across datasets, which could potentially confound genuine biological variations. (B) Subsequent application of the sva package for batch effect correction resulted in a more uniform distribution of samples across the two-dimensional plane, thereby substantially mitigating the batch effect and enhancing result accuracy. (C) The chromosomal locations of the eight genes.\n(TIF)\n(DOCX)\nThe compressed archive contains Fibrate target gene information.xlsx, Niacin target gene information.xlsx, and ω-3 PUFA target gene information.\n(ZIP)","source_license":"public-domain-us","license_restricted":false}