- 1Department of Breast Surgery, General Surgery Center, First Hospital of Jilin University, Changchun, Jilin, China
- 2Department of Ophthalmology, The Second Hospital of Jilin University, Changchun, Jilin, China
Objective: Previous studies have established a causal relationship between metabolites and breast cancer (BC), but the underlying mechanisms remain unclear. Thus, we aimed to investigate the genetic relationship between metabolites and BC, including its subtypes, using Mendelian randomization (MR) analysis.
Methods: Utilizing the latest and most comprehensive summary statistics from genome-wide association studies we conducted an Mendelian randomization study. Data on 233 metabolites, used as exposure variables, were obtained from a study involving 136,016 participants. BC data, used as outcome variables, were sourced from a study comprising 122,977 cases and 105,974 controls. We used the inverse-variance weighted method as the primary approach, along with three supplementary methods, to assess the causal relationship. We also used Cochran’s Q test to detect heterogeneity and MR-Egger regression to examine the presence of horizontal pleiotropy.
Results: Upon analyzing 233 metabolites across 11 classes in relation to BC, we found six classes of metabolites (fatty acids glycerides and phospholipids, lipoprotein subclasses, lipids, apolipoproteins, and lipoprotein particle size) associated with overall BC. Five classes of metabolites (fatty acids glycerides and phospholipids, lipoprotein subclasses, lipids, and lipoprotein particle size) were related to estrogen receptor (ER) + BC, and eight classes of metabolites (fatty acids, amino acids, glycerides and phospholipids, lipoprotein subclasses, lipids, apolipoproteins, glycolysis-related metabolites, and lipoprotein particle size) were linked to ER- BC.
Conclusion: Our study demonstrates a genetic causal relationship between most metabolites and BC, confirming the link between these factors. This research provides a significant foundation for the prevention and treatment of BC.
1 Introduction
Breast cancer (BC) remains the most common cancer among women globally. In 2024, the United States is expected to see 313,510 new cases—2,790 in men and 310,720 in women—and about 42,780 deaths, with 530 involving men and 42,250 women (Siegel et al., 2024). Key risk factors for BC include gender, with women being at higher risk, advancing age, and genetic predispositions. Additional risk factors encompass obesity, alcohol consumption, tobacco use, and hormone replacement therapy (Britt et al., 2020). Symptoms of BC can vary but often include the presence of lumps, changes in breast shape or size, skin indentations, or redness, although early stages may be asymptomatic. The management of BC typically involves a combination of surgical intervention, radiation therapy, and pharmacotherapy, with the choice of treatment depending on the cancer subtype and stage (Li et al., 2022). BC can be classified into estrogen receptor (ER) positive and ER negative types based on the presence of ER in the cancer cells. Approximately 60%–75% of BC are ER + BC, while the remaining 25%–40% are ER- BC(Haibe-Kains et al., 2012). ER- BC tends to grow faster and have a worse short-term prognosis. They often recur within the first few years after treatment, making their management more challenging. Even with the emergence of many novel methods (Faris et al., 2021; Salam et al., 2023), early detection and treatment are critical for improving survival rates.
In recent years, circulating metabolic biomarkers have gained significant attention for their role in understanding metabolic processes and disease mechanisms (Gieger et al., 2008; Illig et al., 2010). These biomarkers include a diverse array of fatty acids, amino acids, lipoprotein subclasses, apolipoproteins, lipids, glycerides and phospholipids, glycolysis-related metabolites, inflammation and lipoprotein particle size. These biomarkers provide detailed insights into the metabolic state of the human body. They play crucial roles in fundamental metabolic activities and are closely linked with various diseases, such as diabetes, and metabolic syndrome. Metabolites are crucial components of metabolic pathways that can influence cancer development and progression. Changes in lipid metabolism have a significant role in BC by affecting cell membrane composition, signaling pathways, and energy production (Zipinotti dos Santos et al., 2023). Understanding the relationship between circulating metabolic biomarkers and BC is essential for uncovering potential metabolic pathways involved in cancer development. Identifying causal links can provide insights into how metabolic alterations influence BC risk. Additionally, this knowledge may inform personalized prevention strategies and therapeutic interventions targeting metabolic processes. Therefore, investigating these relationships is vital for advancing BC research and improving patient outcomes. Several observational studies have investigated the association between metabolic biomarkers and BC. For instance, an observational study by Jennifer C. Melvin et al., involving 1,824 Swedish women diagnosed with BC, reported a modest positive association between serum glucose and the ApoB/ApoA-1 ratio with BC severity (Melvin et al., 2017). This finding suggests that these factors do not significantly contribute to the association between obesity and BC severity. Additionally, a case-control study by Julia Debik et al., which included 1,199 case-control pairs from the Trøndelag Health Study (HUNT), found that several lipoprotein subfractions, particularly VLDL subfractions, were inversely linked with the BC in premenopausal women (Debik et al., 2022). Meanwhile, another study demonstrated that, compared to healthy women, serum TC, LDL-C, and TG levels were significantly elevated in BC patients, while HDL-C levels were significantly reduced, suggesting that lipids may play an important role in the development of BC (Faris et al., 2023). This indicates that alterations in lipid metabolism may occur well in advance of diagnosis of BC. However, there are few observational studies on other types of metabolites, highlighting the need for new approaches to investigate the relationship between metabolites and BC. However, observational studies have limitations, such as the inability to establish the direction of causality and the potential for reverse causation (Sheehan et al., 2008; Sobczyk et al., 2023).
Mendelian randomization (MR) is an epidemiological method that uses genetic variants as instrumental variables (IVs) to determine causal relationships between risk factors and health outcome (Shin et al., 2014; Murphy et al., 2020). By leveraging the random assortment of genes from parents to offspring, MR minimizes reverse causation and confounding which are common issues in observation study (Davies et al., 2018). The primary advantage of MR is its ability to provide more robust causal inferences, similar to randomized controlled trials, but using observational data. This method has been successfully used to investigate the causal effects of various exposures on diseases, offering valuable insights into potential mechanisms and therapeutic targets. By elucidating whether specific metabolic profiles influence BC risk, this study can inform targeted prevention strategies and potential metabolic interventions for BC.
2 Materials and methods
2.1 Study design and MR assumptions
We conducted this MR study to investigate the association between circulating metabolic biomarkers and BC. MR analysis must adhere to three core assumptions: the relevance assumption, the independence assumption, and the exclusion restriction assumption (Bowden and Holmes, 2019). Firstly, the relevance assumption requires that the chosen IVs are strongly correlated with the circulating metabolic biomarkers. Secondly, the independence assumption stipulates that the IVs must not influence BC through any confounding factors that could affect the development of BC. Lastly, the exclusion restriction assumption asserts that the IVs do not directly cause the development of BC. The illustration of the MR process is provided in Figure 1. In this study, we first adhered to the three core assumptions of MR analysis. To address weak IVs, we applied the F-statistic to exclude IVs with an F-value less than 10, ensuring the robustness of our analysis. Additionally, we mitigated potential bias due to linkage disequilibrium by retaining independent IVs. After obtaining the filtered IVs, we conducted analyses using four distinct MR methods: IVW, MR-Egger regression, weighted median, and weighted mode method, to ensure the robustness and consistency of the results. Following the preliminary analysis, we performed various sensitivity analyses to assess the reliability of the findings. First, Cochran’s Q test was used to assess heterogeneity and determine the consistency of effects across different IVs. Second, MR-Egger regression and MR-PRESSO were applied to detect horizontal pleiotropy, ruling out the possibility that the IVs directly influenced the outcome. Furthermore, we conducted reverse MR analysis to check for reverse causality, ensuring the validity of the causal direction. Finally, we performed leave-one-out analysis by removing each IV one at a time and reanalyzing the data to verify that the results were not driven by any single IV. To assess potential publication bias, we utilized a funnel plot, ensuring the accuracy and impartiality of the results. Through this comprehensive series of sensitivity analyses and robustness checks, we arrived at reliable causal inference outcomes.
2.2 Data sources for exposure and outcome data
In this MR study, circulating metabolic biomarkers were selected as exposures, and BC including ER + BC and ER- BC as outcomes. Data on 233 circulating metabolic biomarkers were obtained from the study conducted by Minna et al. (Karjalainen et al., 2024). A total of 136,016 participants from 33 cohorts participated in this study and the authors finally identified more than 400 independent loci and assigned probable causal genes at two-thirds of these using manual curation of plausible biological candidates. The GWAS IDs of the data we used in this study were GCST90301941–GCST90302173 from the GWAS Catalog (https://www.ebi.ac.uk/gwas/) (Supplementary Table S1). Most metabolites have the unit mmol/l, with specific units for each metabolite detailed in Supplementary Table S17. For the outcomes, BC was categorized into three groups: overall BC, ER + BC, and ER- BC. BC genetic association data were sourced from the Breast Cancer Association Consortium (BCAC). This study, conducted by Kyriaki et al., included 122,977 cases and 105,974 controls of European ancestry, combining genetic data from iCOGS, Oncoarray, and GWAS meta-analysis (Michailidou et al., 2017). Detailed information about the data sources is provided in Table 1.
2.3 Selection of instrumental variable
During this process, we selected single nucleotide polymorphisms (SNPs) as IVs. To satisfy the relevance assumption, we set the p-value threshold for IVs at genome-wide significance (p < 1 × 10⁻⁸). Additionally, we chose R2 < 0.001 and a distance of 10,000 kb to avoid linkage disequilibrium with other SNPs. We calculated the F-statistic to assess the validity of the association between the IVs and the exposure (Levin et al., 2020). The formula for the F-statistic is
3 Statistical analysis
In this MR analysis, four methods were used to determine the causal relationship between circulating metabolic biomarkers and BC, including the inverse-variance weighted (IVW) method, the weighted mode method, the weighted median method, and the MR-Egger method. The IVW method, a fixed-effect meta-analysis model, is the primary approach used to combine ratio estimates for each variant into a single causal estimate (Burgess et al., 2015). Weighted median method, weighted mode, and method MR-egger method are complements to the IVW method. The false discovery rate (FDR) method is used in MR studies to control for multiple comparisons, ensuring that findings are robust despite the testing of many interrelated phenotypes (Benjamini and Hochberg, 1995). By setting an FDR threshold of 5%, the method helps confirm the significance of observed effects, adjusting for the potential of false positives due to the correlation among phenotypes. Furthermore, Cochran’s Q test was employed to assess heterogeneity among the IVs. A p-value greater than 0.05 indicated the absence of heterogeneity, justifying the use of the fixed-effects IVW method. Conversely, a p-value less than 0.05 suggested the presence of heterogeneity, necessitating the application of the random-effects IVW method. In addition, odds Ratio (OR) and confidence interval (CI) were used to determine the effect of exposure on outcomes as protective or risk factors. MR-Egger regression is employed to detect horizontal pleiotropy. A p-value greater than 0.05 indicates the absence of horizontal pleiotropy. We use MR-PRESSO to detect outliers of SNPs and remove these. We use the Steiger test to detect whether there is an inverse causality, and p less than 0.05 indicates that there is no inverse causality. Sensitivity analysis included leave-one-out analysis, funnel plots, forest plots, and scatter plots in this MR analysis. Leave-one-out analysis is used to observe the stability of causal relationships.
In this study, R software 4.4.0 (https://www.R-project.org) with MRPRESSO (1.0), MendelianRandomization (version 0.10.0), and TwoSampleMR (version 0.5.8) packages were used to perform the above analysis.
4 Results
233 circulating metabolic biomarkers have been classified into 11 categories: lipoprotein subclasses, lipids, inflammation, fluid balance, ketone bodies, glycolysis-related metabolites, amino acids, fatty acids, glycerides and phospholipids, apolipoproteins, and lipoprotein particle size. Eight of these categories were causally linked to all BC types. Six categories—fatty acids glycerides and phospholipids, lipoprotein subclasses, lipids, apolipoproteins, and lipoprotein particle size—showed a causal association with overall BC, whereas five categories—fatty acids glycerides and phospholipids, lipoprotein subclasses, lipids, and lipoprotein particle size—were associated with ER + BC, and another eight categories (fatty acids, amino acids, glycerides and phospholipids, lipoprotein subclasses, lipids, apolipoproteins, glycolysis-related metabolites, and lipoprotein particle size) are related to ER- BC. This finding suggests that variations in the concentrations of different types of circulating metabolites may contribute to the promotion or inhibition of BC development.
4.1 Causal effects of 233 circulating metabolic biomarkers and overall BC
In this MR study, a comprehensive analysis of 223 metabolic biomarkers was conducted to investigate their potential causal associations with BC risk. Among these, 59 biomarkers demonstrated significant causal relationships with BC, with 35 acting as protective factors (OR < 1) and 24 as risk factors (OR > 1). The metabolites associated with BC were categorized into six major groups: fatty acids, lipids, glycerides and phospholipids, apolipoproteins, lipoprotein subclasses, and lipoprotein particle size. Notably, within the Lipoprotein subclasses category, 48 biomarkers were identified as having causal associations. Of these, 29 were protective factors, including ratio of phospholipids to total lipids in large HDL (P-fdr: 0.035, OR: 0.936, 95%CI: 0.891–0.982), total lipids in large VLDL (P-fdr: 0.026, OR: 0.923, 95%CI: 0.874–0.976), concentration of large VLDL particles (P-fdr: 0.014, OR: 0.919, 95%CI: 0.871–0.969), triglycerides in large VLDL (P-fdr: 0.028, OR: 0.924, 95%CI: 0.874–0.977), the ratio of triglycerides to total lipids in medium HDL (P-fdr: 0.022, OR: 0.923, 95%CI: 0.874–0.974), total lipids in medium VLDL (P-fdr: 0.035, OR: 0.933, 95%CI: 0.887–0.982), concentration of medium VLDL particles (P-fdr: 0.027, OR: 0.930, 95%CI: 0.884–0.978), phospholipids in medium VLDL (P-fdr: 0.040, OR: 0.935, 95%CI: 0.888–0.984), triglycerides in medium VLDL (P-fdr: 0.010, OR: 0.921, 95%CI: 0.879–0.965), triglycerides in small HDL (P-fdr: 0.006, OR: 0.920, 95%CI: 0.880–0.962), triglycerides ratio of small HDL to total lipids (P-fdr: 0.012, OR: 0.922, 95%CI: 0.877–0.968), concentration of small VLDL particles (P-fdr: 0.048, OR: 0.943, 95%CI: 0.900–0.987), phospholipids in small VLDL (P-fdr: 0.035, OR: 0.939, 95%CI: 0.896–0.983), triglycerides in small VLDL (P-fdr: 0.040, OR: 0.939, 95%CI: 0.895–0.985), triglyceride levels in VLDL (P-fdr: 0.015, OR: 0.927, 95%CI: 0.884–0.973), ratio of total cholesterol to total lipids in very large HDL (P-fdr: 0.006, OR: 0.900, 95%CI: 0.851–0.952), ratio of cholesteryl esters to total lipids in very large HDL (P-fdr: 0.004, OR: 0.899, 95%CI: 0.855–0.947), ratio of triglycerides to total lipids in very large HDL (P-fdr: 0.012, OR: 0.903, 95%CI: 0.850–0.960), total cholesterol in very large VLDL (P-fdr: 0.028, OR: 0.920, 95%CI: 0.867–0.976), free cholesterol in very large VLDL (P-fdr: 0.040, OR: 0.920, 95%CI: 0.864–0.980), total lipids in very large VLDL (P-fdr = 0.023, OR: 0.915, 95%CI: 0.862–0.971), concentration of very large VLDL particles (P-fdr: 0.017, OR: 0.912, 95%CI: 0.860–0.968), phospholipids in very large VLDL (P-fdr: 0.028, OR: 0.917, 95%CI: 0.863–0.975), triglycerides in very large VLDL (P-fdr: 0.013, OR: 0.908, 95%CI: 0.855–0.964), free cholesterol levels in chylomicrons and extremely large VLDL (P-fdr: 0.029, OR: 0.918, 95%CI: 0.863–0.976), total lipid levels in chylomicrons and extremely large VLDL (P-fdr: 0.027, OR: 0.921, 95%CI: 0.869–0.975), concentration of chylomicrons and extremely large VLDL particles (P-fdr: 0.029, OR: 0.922, 95%CI: 0.871–0.977), phospholipid levels in chylomicrons and extremely large VLDL (P-fdr: 0.026, OR: 0.925, 95%CI: 0.877–0.976) and triglyceride levels in chylomicrons and extremely large VLDL (P-fdr: 0.015, OR: 0.919, 95%CI: 0.870–0.970). Conversely, 19 biomarkers were identified as risk factors, including free cholesterol to total lipids ratio in IDL (P-fdr: 0.026, OR: 1.074, 95%CI: 1.023–1.129), total cholesterol in large HDL (P-fdr: 0.004, OR: 1.087, 95%CI: 1.043–1.132), ratio of total cholesterol to total lipids in large HDL (P-fdr: 0.012, OR: 1.100, 95%CI: 1.038–1.165), cholesterol esters in large HDL (P-fdr: 0.004, OR: 1.086, 95%CI: 1.043–1.131), free cholesterol in large HDL (P-fdr: 0.007, OR: 1.082, 95%CI: 1.036–1.130), ratio of free cholesterol to total lipids in large HDL (P-fdr: 0.012, OR: 1.107, 95%CI: 1.042–1.177), total lipids in large HDL (P-fdr: 0.005, OR: 1.087, 95%CI: 1.041–1.135), concentration of large HDL particles (P-fdr: 0.010, OR: 1.077, 95%CI: 1.032–1.125), phospholipids in large HDL (P-fdr: 0.009, OR: 1.087, 95%CI: 1.038–1.140), ratio of total cholesterol to total lipids in medium HDL (P-fdr: 0.035, OR: 1.077, 95%CI: 1.020–1.137), ratio of medium VLDL cholesterol to medium VLDL lipids (P-fdr: 0.040, OR: 1.067, 95%CI: 1.016–1.120), total cholesterol levels in very large HDL (P-fdr: 0.013, OR: 1.075, 95%CI: 1.028–1.125), cholesterol esters in very large HDL (P-fdr: 0.042, OR: 1.067, 95%CI: 1.015–1.122), free cholesterol in very large HDL (P-fdr: 0.005, OR: 1.084, 95%CI: 1.039–1.130), total lipids in very large HDL (P-fdr: 0.012, OR: 1.070, 95%CI: 1.028–1.115), concentration of very large HDL particles (P-fdr: 0.012, OR: 1.070, 95%CI: 1.029–1.113), phospholipids in very large HDL (P-fdr: 0.007, OR: 1.079, 95%CI: 1.035–1.125), ratio of phospholipids to total lipids in very large HDL (P-fdr: 0.004, OR: 1.123, 95%CI: 1.065–1.184) and ratio of phospholipids to total lipids in very small VLDL (P-fdr: 0.045, OR: 1.065, 95%CI: 1.014–1.118). In the lipid category, both total cholesterol levels in HDL (OR: 1.083, 95%CI: 1.032–1.136, P-fdr: 0.012) and total cholesterol in HDL2 (P-fdr: 0.010, OR: 1.085, 95%CI: 1.036–1.136) were identified as risk factors for BC. Within the glycerides and phospholipids category, the ratio of diacylglycerol to triglycerides (P-fdr: 0.013, OR: 0.829, 95%CI: 0.738–0.931) and the ratio of triglycerides to phosphoglycerides (P-fdr: 0.034, OR: 0.931, 95%CI: 0.884–0.981) were recognized as protective factors. Among the apolipoproteins, elevated apolipoprotein A-I levels (P-fdr: 0.013, OR: 1.076, 95%CI: 1.028-1.126) were linked to a higher risk of BC. In the fatty acids category, three biomarkers, namely, ratio of conjugated linoleic acid to total fatty acids (P-fdr: 0.013, OR: 0.805, 95%CI: 0.705–0.920), conjugated linoleic acid (P-fdr: 0.012, OR: 0.873, 95%CI: 0.805–0.946) and ratio of monounsaturated fatty acids to total fatty acids (P-fdr: 0.016, OR: 0.922, 95%CI: 0.876–0.972) were found to be protective factors. However, the ratio of omega-6 fatty acids to total fatty acids (P-fdr: 0.012, OR: 1.131, 95%CI: 1.051–1.218) was associated with increased risk. In the category of lipoprotein particle size, the mean diameter of VLDL particles (P-fdr: 0.026, OR: 0.927, 95%CI: 0.881–0.977) was found to be protective, whereas the mean diameter of HDL particles (P-fdr: 0.005, OR: 1.083, 95%CI: 1.040–1.129) posed an increased risk. Sensitivity analyses indicated significant heterogeneity in 56 of the metabolic biomarkers. Consequently, the IVW method was employed for analysis. Additionally, the Steiger test did not detect any reverse causality (Supplementary Table S14). No pleiotropy was detected in the study. Detailed findings on heterogeneity and pleiotropy are presented in Supplementary Tables S5, 8, and the causal effects of the circulating metabolic biomarkers on overall breast cancer risk are depicted in Figure 2. The results are shown in Table 2.
4.2 Causal effects of 233 circulating metabolic biomarkers and ER + BC
In our study on ER + BC, we identified a total of 60 circulating metabolic biomarkers with causal relationships to BC. Among these, 39 were found to be protective factors against BC (OR < 1), while 21 were identified as risk factors (OR > 1). Specifically, within the lipoprotein subclasses, 34 biomarkers including total cholesterol in large VLDL (P-fdr: 0.036, OR: 0.916, 95%CI: 0.859–0.977), cholesterol esters in large VLDL (P-fdr: 0.034, OR: 0.920, 95%CI: 0.867–0.977), total lipids in large VLDL (P-fdr: 0.022, OR: 0.912, 95%CI: 0.858–0.970), concentration of large VLDL particles (P-fdr: 0.014, OR: 0.908, 95%CI: 0.854–0.964), triglycerides in large VLDL (P-fdr: 0.034, OR: 0.916, 95%CI: 0.859–0.976), the ratio of triglycerides to total lipids in medium HDL (P-fdr: 0.014, OR: 0.907, 95%CI: 0.853–0.963), total lipids in medium VLDL (P-fdr: 0.030, OR: 0.919, 95%CI: 0.866–0.975), concentration of medium VLDL particles (P-fdr: 0.027, OR: 0.917, 95%CI: 0.864–0.974), phospholipids in medium VLDL (P-fdr: 0.034, OR: 0.921, 95%CI: 0.868–0.978), triglycerides in medium VLDL (P-fdr: 0.014, OR: 0.909, 95%CI: 0.858–0.962), triglycerides in small HDL (P-fdr: 0.012, OR: 0.907, 95%CI: 0.860–0.956), triglycerides ratio of small HDL to total lipids (P-fdr: 0.014, OR: 0.906, 95%CI: 0.856–0.960), total lipids in small VLDL (P-fdr: 0.042, OR: 0.928, 95%CI: 0.877–0.982), concentration of small VLDL particles (P-fdr: 0.035, OR: 0.927, 95%CI: 0.877–0.980), phospholipids in small VLDL (P-fdr: 0.030, OR: 0.921, 95%CI: 0.869–0.976), triglycerides in small VLDL (P-fdr: 0.042, OR: 0.928, 95%CI: 0.877–0.982), triglyceride levels in VLDL (P-fdr: 0.019, OR: 0.917, 95%CI: 0.867–0.971), ratio of total cholesterol to total lipids in very large HDL (P-fdr: 0.012, OR: 0.896, 95%CI: 0.844–0.952), ratio of cholesteryl esters to total lipids in very large HDL (P-fdr: 0.014, OR: 0.902, 95%CI: 0.852–0.956), ratio of triglycerides to total lipids in very large HDL (P-fdr: 0.014, OR: 0.882, 95%CI: 0.821–0.948), total cholesterol in very large VLDL (P-fdr: 0.014, OR: 0.905, 95%CI: 0.851–0.962), cholesterol esters in very large VLDL (P-fdr: 0.042, OR: 0.919, 95%CI: 0.861–0.980), free cholesterol in very large VLDL (P-fdr: 0.014, OR: 0.901, 95%CI: 0.845–0.962), total lipids in very large VLDL (P-fdr: 0.014, OR: 0.901, 95%CI: 0.845–0.960), concentration of very large VLDL particles (P-fdr: 0.014, OR: 0.898, 95%CI: 0.843–0.957), phospholipids in very large VLDL (P-fdr: 0.015, OR: 0.903, 95%CI: 0.846–0.963), triglycerides in very large VLDL (P-fdr: 0.014, OR: 0.895, 95%CI: 0.840–0.954), total cholesterol levels in chylomicrons and extremely large VLDL (P-fdr: 0.048, OR: 0.919, 95%CI: 0.860–0.982), cholesteryl ester levels in chylomicrons and extremely large VLDL (P-fdr: 0.043, OR: 0.918, 95%CI: 0.860–0.980), free cholesterol levels in chylomicrons and extremely large VLDL (P-fdr: 0.014, OR: 0.904, 95%CI: 0.849–0.963), total lipid levels in chylomicrons and extremely large VLDL (P-fdr: 0.014, OR: 0.905, 95%CI: 0.853–0.961), concentration of chylomicrons and extremely large VLDL particles (P-fdr: 0.014, OR: 0.907, 95%CI: 0.856–0.962), phospholipid levels in chylomicrons and extremely large VLDL (P-fdr: 0.014, OR: 0.913, 95%CI: 0.864–0.965) and triglyceride levels in chylomicrons and extremely large VLDL (P-fdr: 0.014, OR: 0.905, 95%CI: 0.854–0.958) exhibited a negative association with BC, and 17 biomarkers including free cholesterol to total lipids ratio in IDL (P-fdr: 0.048, OR: 1.074, 95%CI: 1.015–1.136), total cholesterol in large HDL (P-fdr: 0.011, OR: 1.095, 95%CI: 1.044–1.148), ratio of total cholesterol to total lipids in large HDL (P-fdr: 0.019, OR: 1.105, 95%CI: 1.035–1.180), cholesterol esters in large HDL (P-fdr: 0.011, OR: 1.094, 95%CI: 1.043–1.147), ratio of cholesteryl esters to total lipids in large HDL (P-fdr: 0.049, OR: 1.086, 95%CI: 1.018–1.159), free cholesterol in large HDL (P-fdr: 0.014, OR: 1.088, 95%CI: 1.036–1.144), ratio of free cholesterol to total lipids in large HDL (P-fdr: 0.014, OR: 1.117, 95%CI: 1.043–1.196), total lipids in large HDL (P-fdr: 0.011, OR: 1.100, 95%CI: 1.046–1.157), concentration of large HDL particles (P-fdr: 0.017, OR: 1.081, 95%CI: 1.028–1.136), phospholipids in large HDL (P-fdr: 0.014, OR: 1.097, 95%CI: 1.039–1.157), ratio of total cholesterol to total lipids in medium HDL (P-fdr: 0.014, OR: 1.102, 95%CI: 1.038–1.171), ratio of cholesteryl esters to total lipids in medium HDL (P-fdr: 0.030, OR: 1.085, 95%CI: 1.024–1.148), free cholesterol in very large HDL (P-fdr: 0.021, OR: 1.081, 95%CI: 1.027–1.138), total lipids in very large HDL (P-fdr: 0.038, OR: 1.068, 95%CI: 1.017–1.122), concentration of very large HDL particles (P-fdr: 0.037, OR: 1.068, 95%CI: 1.017–1.120), phospholipids in very large HDL (P-fdr: 0.014, OR: 1.079, 95%CI: 1.030–1.132) and ratio of phospholipids to total lipids in very large HDL (P-fdr: 0.002, OR: 1.141, 95%CI: 1.078–1.209) showed a positive association. Regarding lipids, total cholesterol in HDL2 (P-fdr: 0.019, OR: 1.084, 95%CI: 1.028–1.142) and total cholesterol levels in HDL (P-fdr: 0.034, OR: 1.079, 95%CI: 1.021–1.140) were determined to be risk factors for ER + BC. Within the category of glycerides and phospholipids, ratio of diacylglycerol to triglycerides (P-fdr: 0.033, OR: 0.833, 95%CI: 0.731–0.949), ratio of triglycerides to phosphoglycerides (P-fdr: 0.040, OR: 0.920, 95%CI: 0.864–0.979) were identified as protective factors against ER + BC. In the fatty acids category, conjugated linoleic acid (P-fdr: 0.022, OR: 0.859, 95%CI: 0.776–0.951), monounsaturated fatty acids (16:1, 18:1) levels (P-fdr: 0.045, OR: 0.926, 95%CI: 0.872–0.983) was found to be a protective factor for ER + BC, whereas the ratio of omega-6 fatty acids to total fatty acids (P-fdr: 0.014, OR: 1.146, 95%CI: 1.057–1.243) was identified as a risk factor. For lipoprotein particle size, the mean diameter of VLDL particles (P-fdr: 0.022, OR: 0.917, 95%CI: 0.865–0.972) was found to be a protective factor for ER + BC, while the mean diameter of HDL particles (P-fdr: 0.011, OR: 1.090, 95%CI: 1.041–1.141) was a risk factor. Our analysis identified 58 instances of heterogeneity but did not reveal any evidence of horizontal pleiotropy. Additionally, the Steiger test did not detect any reverse causality (Supplementary Table S15). Detailed findings on heterogeneity and pleiotropy are presented in Supplementary Tables 6, 10, and the causal effects of the circulating metabolic biomarkers on ER + BC risk are depicted in Figure 3. The results are shown in Table 3.
4.3 Causal effects of 233 circulating metabolic biomarkers and ER- BC
Using the IVW method as the primary analytical approach, we found 30 metabolites with causal relationships to ER- BC, although these did not pass the Bonferroni correction. Among these 30 metabolites, 21 were risk factors for ER- BC, and 9 were protective factors. Within the lipoprotein subclasses, there were 18 metabolites in total, with 5 metabolites including ratio of phospholipids to total lipids in large HDL (P: 0.005, OR: 0.901, 95%CI: 0.837–0.970), the ratio of triglycerides to total lipids in medium HDL (P: 0.041, OR: 0.924, 95%CI: 0.857–0.997), ratio of total cholesterol to total lipids in very large HDL (P: 0.009, OR: 0.884, 95%CI: 0.805–0.970), ratio of cholesteryl esters to total lipids in very large HDL (P: 0.0005, OR: 0.860, 95%CI: 0.791–0.936) and ratio of triglycerides to total lipids in very large HDL (P: 0.043, OR: 0.911, 95%CI: 0.832–0.997) identified as protective factors, and 13 metabolites including total cholesterol in large HDL (P: 0.022, OR: 1.089, 95%CI: 1.012–1.172), ratio of total cholesterol to total lipids in large HDL (P: 0.001, OR: 1.134, 95%CI: 1.052–1.222), cholesterol esters in large HDL (P: 0.023, OR: 1.089, 95%CI: 1.012–1.172), free cholesterol in large HDL (P: 0.028, OR: 1.087, 95%CI: 1.009–1.170), ratio of free cholesterol to total lipids in large HDL (P: 0.031, OR: 1.108, 95%CI: 1.009–1.216), concentration of large HDL particles (P: 0.047, OR: 1.074, 95%CI: 1.001–1.152), phospholipids in large HDL (P: 0.042, OR: 1.081, 95%CI: 1.003–1.165), total cholesterol levels in very large HDL (P: 0.030, OR: 1.082, 95%CI: 1.008–1.162), free cholesterol in very large HDL (P: 0.013, OR: 1.085, 95%CI: 1.017–1.158), total lipids in very large HDL (P: 0.010, OR: 1.080, 95%CI: 1.018–1.146), concentration of very large HDL particles (p = 0.012, OR: 1.078, 95%CI: 1.017–1.142), phospholipids in very large HDL (P: 0.009, OR: 1.085, 95%CI: 1.021–1.154) and ratio of phospholipids to total lipids in very large HDL (P: 0.038, OR: 1.109, 95%CI: 1.006–1.223) identified as risk factors. In the lipids category, total cholesterol in HDL2 (P: 0.019, OR: 1.092, 95%CI: 1.014–1.175) and total cholesterol levels in HDL (P: 0.005, OR: 1.105, 95%CI: 1.030–1.185) were protective factors for ER- BC. For glycerides and phospholipids, the ratio of diacylglycerol to triglycerides (P: 0.049, OR: 0.794, 95%CI: 0.631–0.999) was a protective factor for ER- BC. Within the apolipoproteins, the apolipoprotein A-I levels (P: 0.024, OR: 1.107, 95%CI: 1.013–1.209) were identified as a risk factor for ER- BC. In the fatty acids category, the ratio of conjugated linoleic acid to total fatty acids (P: 0.008, OR: 0.718, 95%CI: 0.563–0.916), ratio of monounsaturated fatty acids to total fatty acids (p = 0.002, OR: 0.876, 95%CI: 0.807–0.951) were protective factors for ER- BC, whereas the ratio of omega-6 fatty acids to total fatty acids (P: 0.020, OR: 1.130, 95%CI: 1.019–1.253) and the ratio of polyunsaturated fatty acids to total fatty acids (P: 0.036, OR: 1.115, 95%CI: 1.007–1.234) were risk factors. Among glycolysis-related metabolites, the citrate levels (P: 0.0003, OR: 1.260, 95%CI: 1.111–1.429) were found to be a risk factor for ER- BC. In the amino acids category, the isoleucine levels (P: 0.012, OR: 0.838, 95%CI: 0.730–0.962) were protective factors, while the glycine levels (P: 0.011, OR: 1.067, 95%CI: 1.015–1.121) were risk factors. For lipoprotein particle size, the mean diameter of HDL particles (P: 0.005, OR: 1.088, 95%CI: 1.025–1.155) was identified as a risk factor for ER- BC. Nineteen instances of heterogeneity were detected (Supplementary Table S7). Two instances of horizontal pleiotropy were identified using the MR-Egger method: the ratio of diacylglycerol to triglycerides (P: 0.049, OR: 0.794, 95%CI: 0.631–0.999) and the ratio of phospholipids to total lipids in very large HDL (P: 0.038, OR: 1.109, 95%CI: 1.006–1.223) (Supplementary Table S10). Consequently, these two metabolites were excluded from the analysis. No other horizontal pleiotropy was detected. Additionally, the Steiger test did not identify any reverse causality (Supplementary Table S16). The causal effects of the metabolites on ER- BC risk are depicted in Figure 4. The results are shown in Table 4.
5 Discussion
In our study, we utilized MR analysis to investigate the causal relationships between 233 circulating metabolic biomarkers and BC at a genetic level. We identified 59 circulating metabolic biomarkers that exhibit a causal relationship with overall BC following FDR method. Additionally, 60 circulating metabolic biomarkers were found to have a causal association with ER + BC, also confirmed through Bonferroni correction. Furthermore, 30 circulating metabolic biomarkers demonstrated a causal link with ER- BC via the IVW method, although they did not pass the Bonferroni correction. As far as we are aware, this is the pioneering study that employs MR analysis to examine the causal links between 233 circulating metabolic biomarkers and BC at the genetic level. In discussing the relationship between circulating metabolic biomarkers and BC, numerous studies have indicated that specific metabolic biomarkers are significantly associated with the onset, progression, and prognosis of BC. For example, research has shown that plasma levels of Transforming Growth Factor β1 are significantly elevated in BC individuals and correlate with treatment response and risk of recurrence (Veyssière et al., 2022). Additionally, serum amino acid metabolism, particularly levels of serine and glutamine, has been linked with drug resistance and overall survival in BC patients (Wang et al., 2023).
In our investigation, we established a causal relationship at the genetic level between 233 circulating metabolic biomarkers—including lipoproteins, glycerides and phospholipids, apolipoproteins, and lipoprotein particle size—and all types of BC. This evidence highlights the significant genetic interactions involving metabolic biomarkers in the pathogenesis of BC, providing a robust basis for the development of targeted therapies and preventive measures.
In our study, we observed a significant correlation between lipoprotein levels and the incidence of BC. Lipoproteins play a significant role in cancer through their influence on lipid metabolism and antioxidative activities. Similarly, evidence from previous observational studies aligns with our MR findings. A case-control study involving 1,199 pairs by Julia Debik et al., demonstrated that various lipoprotein subfractions, especially VLDL subfractions, were inversely linked with BC risk in premenopausal women. (Debik et al., 2022). Moreover, lipid metabolism is deeply intertwined with tumorigenesis. Changes in lipid metabolism may affect several cancer-related pathways including inflammation, oxidative stress, and angiogenesis (Mazzuferi et al., 2021; Allegra et al., 2024). The relationship between lipoprotein levels and BC risk has been extensively studied, with various findings suggesting a complex interaction. For instance, an MR study led by Johnson et al. found that elevated levels of both LDL and HDL cholesterol are linked with a higher risk of BC. This study underlines the potential causal relationship of lipoprotein levels on breast cancer, emphasizing the need for further investigation into the underlying mechanisms and potential therapeutic targets for modifying these lipid levels (Johnson et al., 2020). Another study conducted by Debik et al. reported that specific subfractions of VLDL particles, including VLDL-2, VLDL-3, and VLDL-4, are inversely associated with the risk of BC in premenopausal women. This study underscores significant metabolic changes in the endogenous lipid pathway that can appear long before a breast cancer diagnosis (Debik et al., 2021). These findings suggest that while some lipoproteins may increase risk, others could potentially have a protective effect. This complexity underscores the importance of detailed lipid profiling in understanding breast cancer risk and developing targeted interventions.
Glycerides and phospholipids, essential components of cell membranes, are significantly altered in cancerous tissues compared to normal tissues. In exploring the causal relationship between glycerides and phospholipids, and BC, it is essential to consider various biochemical and genetic factors that may contribute to disease progression. Research indicates that glycerides and phospholipids may play significant roles in the development of BC. A study by Lin et al. explored how alterations in glyceride and phospholipid metabolism might influence BC risk. Their findings suggested that disruptions in these lipid pathways could be linked to the development of BC, proposing a potential mechanism involving lipid metabolism in cancer progression (Au Yeung and Gill, 2023). Another significant piece of research by Kim et al. investigated the relationship between blood levels of certain phospholipids and BC risk. They observed that elevated levels of specific phospholipids were linked to an increased risk of BC (Kerber et al., 2008). These insights underline the complexity of lipid-related pathways in BC and highlight the potential for targeting these mechanisms in therapeutic strategies. As such, understanding the interplay of genetic factors and lipid metabolism could open new avenues for personalized medicine strategies in BC treatment and prevention.
Apolipoproteins are involved in the management of lipid transport and metabolism and have shown varying impacts on BC based on their structure and function. The relationship between apolipoproteins and BC is increasingly recognized in scientific literature, highlighting their roles in both promoting and inhibiting tumor growth depending on their specific types and interactions within the body. For instance, the distribution of apolipoproteins in different lipoprotein fractions appears associated with the severity of BC, suggesting their potential utility in understanding disease progression (Bobin-Dubigeon et al., 2022). Moreover, certain apolipoproteins like ApoA1 and ApoE have shown potential as biomarkers in male BC due to their upregulated levels in serum, which correlates with disease presence and progression (He et al., 2022). Studies have highlighted how the dysregulation of apolipoproteins can influence the onset of BC, with different apolipoproteins playing distinct roles in cancer cell dynamics. For example, ApoC1 is found to promote breast tumorigenesis by altering cell adhesion and migration processes, while ApoM has inhibitory effects on cancer cell invasion and proliferation, potentially by its interaction with vitamin D receptors (He et al., 2022). These studies underline the importance of apolipoproteins in BC, indicating that their modulation could offer new avenues for targeted therapies and diagnostics. As research continues to unravel the specific mechanisms by which apolipoproteins influence cancer biology, their potential as therapeutic targets and diagnostic markers will become increasingly important in managing BC.
Lipoprotein particle size refers to the dimensions of lipoproteins in the blood, which vary in size and composition. These particles, such as HDL and LDL, play roles in transporting cholesterol. Research indicates a connection between the sizes and levels of these particles, particularly HDL, and cancer risk (Pedersen et al., 2020). The relationship between lipoprotein particle size and BC risk appears to be nuanced, with different subclasses of lipoproteins demonstrating varied associations with cancer risk. Research using MR analysis and nuclear magnetic resonance (NMR) has indicated causal links between certain lipoprotein traits and BC risk. Notably, LDL receptor expression in BC cells, such as in certain ER-cell lines, suggests a potential role in lipid uptake and storage which may influence tumor characteristics and progression (Okekunle et al., 2022).
Consistent with our experimental findings, several studies have shown that elevated LDL levels are associated with increased malignancy in breast cancer, potentially through mechanisms related to lipid metabolism (Yuan et al., 2023). One plausible mechanism involves the LDL receptor, a common receptor for LDL. After binding to LDL, LDLR facilitates its internalization into the cell via endocytosis, thereby enhancing fatty acid metabolism in cancer cells. Additional studies have demonstrated that LDLR is highly expressed in BC, contributing to elevated blood cholesterol levels in patients and correlating with poor prognosis in those with breast cancer (de Gonzalo-Calvo et al., 2015; Pires et al., 2012). The overexpression of LDLR increases the uptake of circulating LDL by BC cells, promoting increased metabolic activity and further enhancing the malignancy of the cancer (Rodrigues dos Santos et al., 2023). Moreover, another mechanism that may promote breast cancer development is related to epithelial-mesenchymal transition (EMT). LDL uptake can lead to elevated cholesterol levels in cancer cells, increasing the production of 27-hydroxycholesterol, which has been shown to stimulate both cancer cell proliferation and EMT in BC cells (Cruz et al., 2010; Sierralta et al., 2011). These findings suggest that lipid metabolism plays a crucial role in BC progression.
Previous analyses have already examined the relationship between metabolites and BC. For example, the study by Xiaosheng Zhu et al. analyzed 249 plasma metabolites from the UK Biobank, finding that total cholesterol in HDL and LDL lipid levels were associated with both ER+ and ER- BC (Zhu et al., 2024). However, our study did not find a relationship between acetate and BC, which was confirmed in their experiments. Addressing this issue requires a larger database and more detailed subgroup analyses. Additionally, another MR study indicated that elevated HDL-C (OR-1.08) might increase the risk of BC (Zhou et al., 2023). However, this study did not further analyze the relationship between specific components of HDL and BC. These discrepancies highlight the need for further research to better understand the complex relationship between metabolites and BC.
In our MR study, we successfully leveraged genetic proxies to elucidate the causal relationships between 233 circulating metabolic biomarkers and BC, including its subtypes ER + BC and ER- BC. Our approach has several advantages as follows. Firstly, by utilizing genetic variants as instruments in MR analysis, we significantly reduced the impact of confounding factors that typically affect observational studies. This approach enables a more precise understanding of the causal effects of metabolites on BC risk. Secondly, we analyzed 233 circulating metabolic biomarkers, providing a broad and detailed panorama of circulating metabolic biomarkers’ role in BC. This comprehensive approach helps in identifying potential biomarkers and therapeutic targets.
Our experimental results have significant clinical implications. Firstly, these biomarkers can be used for the early diagnosis and risk prediction of BC. Since we have identified several metabolic biomarkers significantly associated with BC, future developments in blood-based testing tools could enable more precise risk assessment and personalized screening strategies. For instance, elevated levels of certain metabolic biomarkers may serve as warning signs for high-risk individuals, aiding clinicians in adopting more proactive monitoring and preventive measures. Secondly, these metabolic biomarkers may serve as potential therapeutic targets. By targeting the metabolic pathways causally linked to breast cancer, new treatment strategies could be developed. Particularly, the biomarkers identified in our study could help formulate tailored treatment plans for different subtypes of BC, such as ER+ and ER-subtypes. In summary, this research not only provides new insights into the metabolic mechanisms of breast cancer but also offers potential avenues for future clinical applications and translational research.
While our MR study offers valuable insights into the causal links between circulating metabolic biomarkers and BC, it is important to acknowledge a few limitations. The genetic instruments used might not always be strong enough, which could compromise the reliability of our causal inferences, especially if these variants explain only a minimal part of the variation in lipid levels. Additionally, despite employing MR-Egger to check for pleiotropy, the issue of horizontal pleiotropy remains, where genetic variants might influence multiple traits through unrelated pathways, potentially skewing the results. Unfortunately, in this experiment, the ancestry of the population for the metabolite data was diverse, yet for the outcome analysis, we only selected individuals of European descent. We searched through numerous databases, and the metabolite GWAS database we used was the most comprehensive one we could find. However, this may introduce bias into the results of the study. Moreover, the generalizability of our findings could be limited, as they may not apply to different populations due to variations in genetic background, lifestyle, or environmental factors that influence lipid metabolism and BC risk. This necessitates cautious interpretation and underscores the importance of validating our results across diverse groups. At the last, there are multiple classification methods for BC, including the commonly recognized subtypes: Luminal A, Luminal B, HER2-overexpressing, and TNBC. Additionally, BC can also be categorized based on ER expression levels into ER-negative BC and ER-low BC. These molecular subtypes differ significantly in terms of characteristics, prognosis, and aggressiveness. However, publicly available databases currently lack SNP data with such detailed classifications. Therefore, this study focuses on exploring the relationship between metabolites and these two BC subtypes based on ER expression levels. Access to more detailed classification data in the future would allow for further subtype-specific research, deepening our understanding of the relationship between metabolites and BC.
6 Conclusion
Our MR study offers new insights into the relationship between metabolites and BC. In our investigation, most metabolites were found to have a causal relationship with the occurrence of BC. Specifically, the six classes of metabolites—fatty acids glycerides and phospholipids, lipoprotein subclasses, lipids, apolipoproteins, and lipoprotein particle size—showed causal associations with overall BC. Five classes of metabolites—fatty acids glycerides and phospholipids, lipoprotein subclasses, lipids, and lipoprotein particle size—were associated with ER + BC, while eight classes of metabolites—fatty acids, amino acids, glycerides and phospholipids, lipoprotein subclasses, lipids, apolipoproteins, glycolysis-related metabolites, and lipoprotein particle size—were linked to ER- BC. This study enhances our insight into the association between metabolites and BC, offering a significant basis for future research and potential prevention strategies for BC.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
BW: Formal Analysis, Methodology, Writing–original draft, Writing–review and editing. YL: Data curation, Investigation, Software, Writing–review and editing. HZ: Conceptualization, Data curation, Visualization, Writing–review and editing. MY: Conceptualization, Methodology, Supervision, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
We extend our gratitude to the participants and investigators of the IEU Open GWAS project and the Breast Cancer Association Consortium for their efforts in making the summary statistics publicly accessible.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2024.1448748/full#supplementary-material
References
Allegra, A., Murdaca, G., Mirabile, G., and Gangemi, S. (2024). Protective effects of high-density lipoprotein on cancer risk: focus on multiple myeloma. Biomedicines 12 (3), 514. doi:10.3390/biomedicines12030514
Au Yeung, S. L., and Gill, D. (2023). Standardizing the reporting of Mendelian randomization studies. BMC Med. 21 (1), 187. doi:10.1186/s12916-023-02894-8
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B Methodol. 57 (1), 289–300. doi:10.1111/j.2517-6161.1995.tb02031.x
Bobin-Dubigeon, C., Nazih, H., Blanchard, V., Croyal, M., and Bard, J.-M. (2022). Circulating HDL and non-HDL associated apolipoproteins and breast cancer severity. J. Clin. Med. 11 (5), 1345. doi:10.3390/jcm11051345
Bowden, J., and Holmes, M. V. (2019). Meta-analysis and Mendelian randomization: a review. Res. synthesis methods 10 (4), 486–496. doi:10.1002/jrsm.1346
Britt, K. L., Cuzick, J., and Phillips, K. A. (2020). Key steps for effective breast cancer prevention. Nat. Rev. Cancer 20 (8), 417–436. doi:10.1038/s41568-020-0266-x
Burgess, S., Scott, R. A., Timpson, N. J., Davey Smith, G., and Thompson, S. G.EPIC- InterAct Consortium (2015). Using published data in Mendelian randomization: a blueprint for efficient identification of causal risk factors. Eur. J. Epidemiol. 30 (7), 543–552. doi:10.1007/s10654-015-0011-z
Cruz, P., Torres, C., RamÍRez, M. E., EpuÑÁN, M. J., Valladares, L. E., and Sierralta, W. D. (2010). Proliferation of human mammary cancer cells exposed to 27-hydroxychole sterol. Exp. Ther. Med. 1 (3), 531–536. doi:10.3892/etm_00000084
Davies, N. M., Holmes, M. V., and Smith, G. D. (2018). Reading Mendelian randomisation studies: a guide, glossary, and checklist for clinicians. bmj 362, k601. doi:10.1136/bmj.k601
Debik, J., Schaefer, H., Andreassen, T., Wang, F., Fang, F., Cannet, C., et al. (2021). Lipoprotein and metabolite associations to breast cancer risk in the HUNT2 study. medRxiv. doi:10.1101/2021.10.08.21264729
Debik, J., Schafer, H., Andreassen, T., Wang, F., Fang, F., Cannet, C., et al. (2022). Lipoprotein and metabolite associations to breast cancer risk in the HUNT2 study. Br. J. Cancer 127 (8), 1515–1524. doi:10.1038/s41416-022-01924-1
de Gonzalo-Calvo, D., López-Vilaró, L., Nasarre, L., Perez-Olabarria, M., Vázquez, T., Escuin, D., et al. (2015). Intratumor cholesteryl ester accumulation is associated with human bre ast cancer proliferation and aggressive potential: a molecular and cli nicopathological study. BMC Cancer 15 (1), 460. doi:10.1186/s12885-015-1469-5
Faris, R. A., Mahdi, Z. F., and Al-Layla, A. M. T. (2021). Preliminary study of the insulin growth factor binding protein-3 (IGFBP3) level in Iraqi women with breast cancer. AIP Conf. Proc. 2372 (1). doi:10.1063/5.0065368
Faris, R. A., Mahdi, Z. F., Al-Tu’ma, M. M. K., and Al-Tu’ma, F. J. (2023). Surface enhanced Raman spectroscopy based sensitive and specific detection of vitamin D3, glycated hemoglobin, and serum lipid profile of breast cancer patients. AIP Conf. Proc. 2830 (1). doi:10.1063/5.0157207
Gieger, C., Geistlinger, L., Altmaier, E., Hrabé de Angelis, M., Kronenberg, F., Meitinger, T., et al. (2008). Genetics meets metabolomics: a genome-wide association study of metabolite profiles in human serum. PLoS Genet. 4 (11), e1000282. doi:10.1371/journal.pgen.1000282
Haibe-Kains, B., Desmedt, C., Loi, S., Culhane, A. C., Bontempi, G., Quackenbush, J., et al. (2012). A three-gene model to robustly identify breast cancer molecular subtypes. J. Natl. Cancer Inst. 104 (4), 311–325. doi:10.1093/jnci/djr545
He, Y., Chen, J., Ma, Y., and Chen, H. (2022). Apolipoproteins: new players in cancers. Front. Pharmacol. 13, 1051280. doi:10.3389/fphar.2022.1051280
Illig, T., Gieger, C., Zhai, G., Römisch-Margl, W., Wang-Sattler, R., Prehn, C., et al. (2010). A genome-wide perspective of genetic variation in human metabolism. Nat. Genet. 42 (2), 137–141. doi:10.1038/ng.507
Johnson, K. E., Siewert, K. M., Klarin, D., Damrauer, S. M., the, V. A. M. V. P., Chang, K.-M., et al. (2020). The relationship between circulating lipids and breast cancer risk: a Mendelian randomization study. PLOS Med. 17 (9), e1003302. doi:10.1371/journal.pmed.1003302
Karjalainen, M. K., Karthikeyan, S., Oliver-Williams, C., Sliz, E., Allara, E., Fung, W. T., et al. (2024). Genome-wide characterization of circulating metabolic biomarkers. Nature 628 (8006), 130–138. doi:10.1038/s41586-024-07148-y
Kerber, R. A., Amos, C. I., Yeap, B. Y., Finkelstein, D. M., and Thomas, D. C. (2008). Design considerations in a sib-pair study of linkage for susceptibility loci in cancer. BMC Med. Genet. 9 (1), 64. doi:10.1186/1471-2350-9-64
Levin, M. G., Judy, R., Gill, D., Vujkovic, M., Verma, S. S., Bradford, Y., et al. (2020). Genetics of height and risk of atrial fibrillation: a Mendelian randomization study. PLoS Med. 17 (10), e1003288. doi:10.1371/journal.pmed.1003288
Li, Y., Zhang, H., Merkher, Y., Chen, L., Liu, N., Leonov, S., et al. (2022). Recent advances in therapeutic strategies for triple-negative breast cancer. J. Hematol. and Oncol. 15 (1), 121. doi:10.1186/s13045-022-01341-0
Mazzuferi, G., Bacchetti, T., Islam, M. O., and Ferretti, G. (2021). High density lipoproteins and oxidative stress in breast cancer. Lipids Health Dis. 20 (1), 143. doi:10.1186/s12944-021-01562-1
Melvin, J. C., Garmo, H., Holmberg, L., Hammar, N., Walldius, G., Jungner, I., et al. (2017). Glucose and lipoprotein biomarkers and breast cancer severity using data from the Swedish AMORIS cohort. BMC Cancer 17 (1), 246. doi:10.1186/s12885-017-3232-6
Michailidou, K., Lindström, S., Dennis, J., Beesley, J., Hui, S., Kar, S., et al. (2017). Association analysis identifies 65 new breast cancer risk loci. Nature 551 (7678), 92–94. doi:10.1038/nature24284
Murphy, N., Knuppel, A., Papadimitriou, N., Martin, R. M., Tsilidis, K., Smith-Byrne, K., et al. (2020). Insulin-like growth factor-1, insulin-like growth factor-binding protein-3, and breast cancer risk: observational and Mendelian randomization analyses with∼ 430 000 women. Ann. Oncol. 31 (5), 641–649. doi:10.1016/j.annonc.2020.01.066
Okekunle, A. P., Yie, G.-E., Song, S., Kim, Z., Youn, H. J., Cho, J., et al. (2022). Association of lipid profile with obesity among breast cancer survivors: a cross-sectional study. Lipids Health Dis. 21 (1), 66. doi:10.1186/s12944-022-01674-2
Pedersen, K. M., Çolak, Y., Bojesen, S. E., and Nordestgaard, B. G. (2020). Low high-density lipoprotein and increased risk of several cancers: 2 population-based cohort studies including 116,728 individuals. J. Hematol. and Oncol. 13 (1), 129. doi:10.1186/s13045-020-00963-6
Pires, L. A., Hegg, R., Freitas, F. R., Tavares, E. R., Almeida, C. P., Baracat, E. C., et al. (2012). Effect of neoadjuvant chemotherapy on low-density lipoprotein (LDL) re ceptor and LDL receptor-related protein 1 (LRP-1) receptor in locally advanced breast cancer. Braz. J. Med. Biol. Res. 45 (6), 557–564. doi:10.1590/s0100-879x2012007500068
Rodrigues dos Santos, C., Domingues, G., Matias, I., Matos, J., Fonseca, I., de Almeida, J. M., et al. (2014). LDL-cholesterol signaling induces breast cancer proliferation and inva sion. Lipids Health Dis. 13 (1), 16. doi:10.1186/1476-511x-13-16
Salam, A., Faris, R. A., and Al-janabi, A. (2023). Ultrafast lithium disilicate veneer debonding time assisted by a CO2 laser with temperature control. Opt. Contin. 2 (4), 825–837. doi:10.1364/OPTCON.485281
Sheehan, N. A., Didelez, V., Burton, P. R., and Tobin, M. D. (2008). Mendelian randomisation and causal inference in observational epidemiology. PLoS Med. 5 (8), e177. doi:10.1371/journal.pmed.0050177
Shin, S.-Y., Fauman, E. B., Petersen, A.-K., Krumsiek, J., Santos, R., Huang, J., et al. (2014). An atlas of genetic influences on human blood metabolites. Nat. Genet. 46 (6), 543–550. doi:10.1038/ng.2982
Siegel, R. L., Giaquinto, A. N., and Jemal, A. (2024). Cancer statistics, 2024. CA Cancer J. Clin. 74 (1), 12–49. doi:10.3322/caac.21820
Sierralta, W., Ramírez, M. E., Cruz, P., Epuñan, M. J., and Valladares, L. E. (2011). 27-Hydroxycholesterol induces the transition of MCF7 cells into a mese nchymal phenotype. Oncol. Rep. 26, 389–397. doi:10.3892/or.2011.1284
Sobczyk, M. K., Zheng, J., Smith, G. D., and Gaunt, T. R. (2023). Systematic comparison of Mendelian randomisation studies and randomised controlled trials using electronic databases. BMJ open 13 (9), e072087. doi:10.1136/bmjopen-2023-072087
Veyssière, H., Bidet, Y., Penault-Llorca, F., Radosevic-Robin, N., and Durando, X. (2022). Circulating proteins as predictive and prognostic biomarkers in breast cancer. Clin. Proteomics 19 (1), 25. doi:10.1186/s12014-022-09362-0
Wang, W., Rong, Z., Wang, G., Hou, Y., Yang, F., and Qiu, M. (2023). Cancer metabolites: promising biomarkers for cancer liquid biopsy. Biomark. Res. 11 (1), 66. doi:10.1186/s40364-023-00507-3
Yuan, Q., Lu, X., Guo, H., Sun, J., Yang, M., Liu, Q., et al. (2023). Low-density lipoprotein receptor promotes crosstalk between cell stemness and tumor immune microenvironment in breast cancer: a large data-based multi-omics study. J. Transl. Med. 21 (1), 871. doi:10.1186/s12967-023-04699-y
Zhou, M., Henricks, M., Loch, V., Zhang, G., Lu, Y., and Li, X. (2023). Mendelian randomization analysis revealed potential metabolic causal factors for breast cancer. Sci. Rep. 13 (1), 14290. doi:10.1038/s41598-023-41130-4
Zhu, X., Huang, H., Zou, M., Luo, H., Liu, T., Zhu, S., et al. (2024). Identification of circulating metabolites linked to the risk of breast cancer: a mendelian randomization study. Front. Pharmacol. 15, 1442723. doi:10.3389/fphar.2024.1442723
Zipinotti dos Santos, D., de Souza, J. C., Pimenta, T. M., da Silva Martins, B., Junior, R. S. R., Butzene, S. M. S., et al. (2023). The impact of lipid metabolism on breast cancer: a review about its role in tumorigenesis and immune escape. Cell Commun. Signal. 21 (1), 161. doi:10.1186/s12964-023-01178-1
Keywords: circulating metabolic biomarkers, metabolite, lipoprotein, lipids, mendelian randomization, breast cancer
Citation: Wang B, Ling Y, Zhang H and Yang M (2024) Assessing the causal relationships between circulating metabolic biomarkers and breast cancer by using mendelian randomization. Front. Genet. 15:1448748. doi: 10.3389/fgene.2024.1448748
Received: 18 June 2024; Accepted: 06 December 2024;
Published: 18 December 2024.
Edited by:
Ferdinand Molnár, Nazarbayev University, KazakhstanReviewed by:
Guochong Jia, Vanderbilt University Medical Center, United StatesRawaa A. Faris, University of Baghdad, Iraq
Copyright © 2024 Wang, Ling, Zhang and Yang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ming Yang, eWFuZ21pbmdAamx1LmVkdS5jbg==