- Jinling Hospital, Affiliated Hospital of Medical School, Nanjing University, Nanjing, China
Aim: Pancreatic cancer (PC) is a devastating malignancy characterized by its aggressive nature and poor prognosis. However, the relationship of PC with peripheral metabolites remains not fully investigated. The study aimed to explore the causal linkage between PC and peripheral metabolite profiles.
Methods: Employing publicly accessible genome-wide association studies (GWAS) data, we conducted a bidirectional two-sample Mendelian randomization (MR) analysis. The primary analysis employed the inverse-variance weighted (IVW) method. To address potential concerns about horizontal pleiotropy, we also employed supplementary methods such as maximum likelihood, weighted median, MR-Egger regression, and MR pleiotropy residual sum and outlier (MR-PRESSO).
Results: We ascertained 20 genetically determined peripheral metabolites with causal linkages to PC while high-density lipoprotein (HDL) and very low-density lipoprotein (VLDL) particles accounted for the vast majority. Specifically, HDL particles exhibited an elevated PC risk while VLDL particles displayed an opposing pattern. The converse MR analysis underscored a notable alteration in 17 peripheral metabolites due to PC, including branch chain amino acids and derivatives of glycerophospholipid. Cross-referencing the bidirectional MR results revealed a reciprocal causation of PC and X-02269 which might form a self-perpetuating loop in PC development. Additionally, 1-arachidonoylglycerophosphocholine indicated a reduced PC risk and an increase under PC influence, possibly serving as a negative feedback regulator.
Conclusion: Our findings suggest a complex interplay between pancreatic cancer and peripheral metabolites, with potential implications for understanding the etiology of pancreatic cancer and identifying novel early diagnosis and therapeutic targets. Moreover, X-02269 may hold a pivotal role in PC onset and progression.
1 Introduction
Pancreatic cancer (PC) is one of the most lethal forms of solid malignancy with a 5-year survival rate lower than 9%, which poses a formidable challenge to global health (Qin et al., 2020; Casolino et al., 2021). The annual incidence and mortality rates of PC are on a steady rise, showing minimal advancements in overall survival outcomes (Khalaf et al., 2021). The asymptomatic nature and rapid metastatic spread of PC contribute to the diagnosis at advanced stages, significantly limiting treatment options and leading to dismal survival rates (Huang et al., 2019; Chen et al., 2021). Another reason of PC’s unfavorable prognosis is that surgery stands as the sole curative approach. Thus, early diagnosis and novel treatment strategies for PC are eagerly wanted.
Metabolic reprogramming, recognized as an emerging hallmark of cancer (Hanahan and Weinberg, 2011), has garnered renewed attention over the last decades (Leal-Esteban and Fajas, 2020; Infantino et al., 2021). Metabolites, whether present within the local tumor or circulating in the peripheral blood, play a multifaceted role. They not only modulated energy dynamics but also orchestrated signal transduction among malignant and nonmalignant cells in tumor microenvironment (Garcia-Bermudez et al., 2018; Schworer et al., 2019; Pavlova et al., 2022). Beyond the local tumor site, the intricate reprogramming of tumor metabolism can exert far-reaching consequences on the host’s systemic metabolic landscape, and further leads to cachexia and antitumor immunity deficiency (Porporato, 2016; Drijvers et al., 2020; DePeaux and Delgoffe, 2021). Indeed, observational studies have reported associations between specific peripheral metabolites in PC patients, including branched-chain amino acids (BCAAs), glycerophospholipids, lipids and lipoprotein (Mayers et al., 2014; Shu et al., 2018; Stolzenberg-Solomon et al., 2020; Elebo et al., 2021; Encarnacion-Rosado and Kimmelman, 2021; Wang et al., 2022). Specifically, elevated plasma levels of BCAAs, including leucine, isoleucine and valine, are associated with an increased risk of future PC diagnosis (Mayers et al., 2014). The glycerophospholipids were preferred to be associated with decreased risk of PC (Shu et al., 2018). Lipids exhibited various effects on PC. Steroid hormones, such as androgenic steroids, pregnenolone steroids and corticosteroids, showed an inverse association with incident PC. Conversely, fatty acids and glycoursodeoxycholic acid sulfate displayed a positive association with PC (Wang et al., 2022). Regarding lipoproteins, various lipoprotein parameters, such as size, particle number, and lipid concentration (cholesterol and triglycerides) in the primary classes of lipoproteins, including very-low-density lipoprotein (VLDL), intermediate-density lipoprotein (IDL), low-density lipoprotein (LDL), and high-density lipoprotein (HDL), showed associations with PC. The size of HDL particles exhibited a positive correlation with PC stage, whereas HDL cholesterol levels and HDL particles showed opposite associations (Elebo et al., 2021). However, these findings originate from observational evidences, which are constrained by sample size limitations and beset by inherent biases, notably including reverse causation and confounding factors. Addressing this causal uncertainty of peripheral metabolites and PC is essential for understanding the potential therapeutic implications and identifying novel targets for interventions.
Recognized for its robustness, Mendelian randomization (MR) leverages genetic variants as instrumental variables (IVs) to unravel causal associations within observational datasets. In this study, we present a comprehensive bidirectional two-sample MR analysis to disentangle the reciprocal effects between peripheral metabolites and pancreatic cancer development. We identified 20 peripheral metabolites as potent contributors to PC while 17 were regulated by PC. Notably, our results revealed mutually positive causal relationship of X-02269 and PC, highlighting that they might form a vicious circle in PC onset and development. Through this investigation, we endeavor to unravel novel insights into the dynamic interplay between peripheral metabolites and PC, contributing to the broader understanding of PC etiology. The identified metabolites hold promise for the development of early diagnosis and targeted interventions for PC, thus mitigating the burden imposed by this formidable malignancy.
2 Materials and methods
2.1 Study design
We employed a bidirectional two-sample MR analysis (Figure 1) to investigate the causal relationship between PC and 842 peripheral metabolites in European-ancestry datasets. The genetic instruments utilized in our study met the three essential assumptions of MR analysis. Firstly, the selected genetic variations, serving as instrumental variables (IVs), exhibited significant associations with the exposure of interest. Secondly, these genetic variations were independent of both known and unknown confounding factors. Lastly, the genetic instruments were specifically associated with the outcome, establishing a plausible causal pathway.
FIGURE 1. Study design of exploring the causal relationship between pancreatic cancer and peripheral metabolites by bidirectional two-sample Mendelian randomization. COPD, chronic obstructive pulmonary disease; SNPs, single-nucleotide polymorphisms.
2.2 Data source
Genetic associations pertaining to human peripheral metabolites were extracted from three distinct research studies (Shin et al., 2014; Kettunen et al., 2016; Borges et al., 2020), which were accessible through MRC Integrative Epidemiology Unit (IEU) OpenGWAS project (https://gwas.mrcieu.ac.uk/). Correspondingly, genetic associations concerning pancreatic cancer were acquired from the FinnGen database (R8 level), encompassing a sample size of 1249 cases and 259583 controls. It is noteworthy that all conducted GWAS involved participants of European ancestry. The comprehensive specifics regarding the sourced GWAS datasets are meticulously outlined in Table 1.
2.3 Mendelian randomization analysis
In preparation for the MR analysis, a meticulous selection process was meticulously undertaken to identify genetic IVs in adherence to specific criteria. Initially, the selection criterion mandated that single nucleotide polymorphisms (SNPs) demonstrate a statistically significant association with the exposure at the genome-wide significance threshold. Subsequently, an extensive range of tests spanning from 5 × 10−8 to 1 × 10−6 were meticulously performed to discern the optimal p-value threshold for peripheral metabolites and pancreatic cancer, respectively. The chosen thresholds, 5 × 10−8 for peripheral metabolites and 5 × 10−6 for pancreatic cancer, were primarily grounded in the abundance of informative SNPs. Furthermore, SNPs with a limited association strength with telomere length (F-statistic ≤10) were judiciously excluded from consideration. The F statistic was calculated employing the equation:
Upon the harmonization of exposure and outcome summary data, the central MR analysis was conducted employing the inverse variance-weighted (IVW) MR to elucidate the potential causal link between peripheral metabolites and pancreatic cancer (Burgess et al., 2013). Additional methods were applied to assessed potential violations of MR assumptions, including maximum likelihood (Burgess et al., 2016), weighted median (Bowden et al., 2016), and MR-Egger (Bowden et al., 2015). The primary method employed was the IVW approach, which assumed the validity of all SNPs used as instrumental variables. In cases where a minimum of 50% of the weight was derived from valid instrument variants, the weighted median (WM) method was employed to ensure consistent estimations. For the evaluation of causal effects while accounting for potential horizontal pleiotropy, the MR-Egger regression method was utilized, albeit with a potential trade-off in precision. Simultaneously, the simple mode method was also utilized, positing a direct causal link between genetic variants associated with the exposure and the outcome, without accommodating the intricacies or biases intrinsic to more intricate MR methodologies. Comprehensive explanations of these methodologies can be obtained from previous scholarly publications (Burgess et al., 2013; Burgess et al., 2016; Bowden et al., 2015; Bowden et al., 2016). To explore the causal relationship between peripheral metabolites and pancreatic cancer, the Steiger test was employed (Hemani et al., 2017).
To verify the robustness of our findings, several sensitivity analyses were executed. Initially, the Cochran’s Q test was applied to evaluate the presence of heterogeneity within the data (Bowden et al., 2019). Subsequently, MR-Pleiotropy Residual Sum and Outlier (MR-PRESSO) methodology were employed to detect and eliminate SNPs with potential horizontal pleiotropic outliers, thus mitigating the potential impact of pleiotropy on the causal estimations (Verbanck et al., 2018). In cases where significant horizontal pleiotropy was detected in the MR-PRESSO global test, SNPs identified as outliers (with a significance level of p < 0.05) were excluded, and the remaining SNPs underwent re-analysis. Furthermore, the MR-Egger regression intercept was employed to assess potential pleiotropy in the SNPs, with a p-value ≥0.05 indicating the absence of horizontal pleiotropy. Additionally, we conducted a two-sample MR analysis on each individual SNP. Lastly, a leave-one-out analysis was carried out to detect any pleiotropy stemming from each individual SNP. All these analytical procedures were performed utilizing the TwoSampleMR packages. The outcomes were presented in terms of odds ratios (OR) alongside their corresponding 95% confidence intervals (CI). All p-values reported were two-sided, and statistical significance was defined at the 5% threshold.
3 Results
3.1 Causal effects of the peripheral metabolites on PC
The GWAS data of peripheral metabolites was obtained from three previously published study, which contained 7824, 24925, and 115078 samples of European ancestry, respectively. Supplementary Table S1 contains comprehensive information on the peripheral-metabolite-associated SNPs, including beta values, standard errors, effect alleles, other alleles, etc. The F-statistic values of selected SNPs ranged from 22.81 to 4176.17, supporting the absence of weak instrument bias.
We first performed MR analysis to investigate the causal effects of peripheral metabolites on PC. The IVW-MR results demonstrated 5 and 15 metabolites which were positively and negatively causal associated with PC, respectively (Figure 2). Leucine (OR 114.74, 95% CI 1.42–9266.7, p = 0.034), X-11469 (OR 5.22, 95% CI 1.78–15.29, p = 0.003) and X-02269 (OR 4.61, 95% CI 1.69–12.57, p = 0.003) were the top three metabolites indicating an increasing PC risk. Remarkably, the HDL and VLDL particles occupied 45% of all significant 20 metabolites. Specifically, cholesterol esters in large HDL (OR 1.29, 95% CI 1.05–1.59, p = 0.017), total cholesterol in large HDL (OR 1.26, 95% CI 1.02–1.56, p = 0.031) and free cholesterol in large HDL (OR 1.25, 95% CI 1.02–1.55, p = 0.036) were associated with an elevated risk of PC whilst free cholesterol in large VLDL (OR 0.74, 95% CI 0.55–0.99, p = 0.043), triglycerides in very large VLDL (OR 0.71, 95% CI 0.51–0.99, p = 0.044), phospholipids in very large VLDL (OR 0.71, 95% CI 0.51–0.99, p = 0.042), total lipids in very large VLDL (OR 0.7, 95% CI 0.5–0.99, p = 0.046), total lipids in chylomicrons and largest VLDL particles (OR 0.66, 95% CI 0.46–0.95, p = 0.023) and concentration of chylomicrons and largest VLDL particles (OR 0.64, 95% CI 0.44–0.94, p = 0.023) demonstrated the opposite trend. In addition, omega-3 fatty acids (OR 0.66, 95% CI 0.47–0.91, p = 0.013), docosahexaenoic acid (22:6) (OR 0.64, 95% CI 0.45–0.92, p = 0.017), X-12063 (OR 0.53, 95% CI 0.32–0.87, p = 0.012), 3-Hydroxybutyrate (OR 0.47, 95% CI 0.26–0.86, p = 0.014), acetoacetate (OR 0.38, 95% CI 0.16–0.92, p = 0.031), 1-arachidonoylglycerophosphocholine (OR 0.25, 95% CI 0.07–0.94, p = 0.041), 1-arachidonoylglycerophosphoinositol (OR 0.18, 95% CI 0.04–0.86, p = 0.032) and uridine (OR 0.01, 95% CI 0–0.76, p = 0.038) were related to reduced PC risk.
FIGURE 2. The effect of genetically determined peripheral metabolites on the risk of pancreatic cancer in the IVW results. OR, odds ratios; CI, confidence interval; IVW, inverse variance-weighted.
Next, potential violations of MR assumptions were assessed by performing additional MR analyses using three alternative methods: maximum likelihood, weighted median and MR Egger (Supplementary Table S2). The overall trends of the effects of peripheral metabolites on PC analyzed by the additional methods mentioned above remained consistent, although some of them did not reach statistical significance (Supplementary Tables S2, S3). In contrast, MR-Egger regression yielded that leucine (OR 1.45e-05, 95% CI 1.73e-17-1.22e+7), X-02269 (OR 0.19, 95% CI 6.30e-5-564.38) and uridine (OR 1.41, 95% CI 1.20e-9-1.65e+11) was in the opposite direction of the IVW estimation, but it also did not reach statistical significance.
To evaluate heterogeneity and identify potential outliers, we conducted Cochran’s Q test and MR-PRESSO test, respectively. The results (Supplementary Table S3) indicated no significant heterogeneity (PQ > 0.05) and no outliers were detected. Additionally, all p-values from the MR-Egger intercept test were greater than 0.05, suggesting the absence of horizontal pleiotropy. Steiger filtering revealed that all genetic instrumental variables used for PC explained more variance in peripheral metabolites than in PC (Supplementary Table S2). To further assess the robustness of our findings, we employed several additional analyses. Firstly, we visually examined the forest plot (Supplementary Figure S1). Secondly, we utilized four methods to assess the results of the MR analysis, and a scatter plot specifically for PC was generated (Supplementary Figure S2). Finally, leave-one-out analysis further confirmed the robustness of our main results (Supplementary Figure S3).
3.2 Causal effects of the PC on peripheral metabolites
We incorporated a total of 11 independent SNPs as IVs for PC from the FinnGen database (Infantino et al., 2021). The detailed information of selected SNP was presented in Supplementary Table S4. No bias of weak instrument was confirmed by F statistics, which ranged from 20.87 to 32.06.
The IVW-MR analysis demonstrated a significant causal relationship of PC with 17 peripheral metabolites (Figure 3). Specifically, PC was associated with increasing of X-12007 (OR 1.05, 95% CI 1.01–1.1, p = 0.019), valine (OR 1.04, 95% CI 1.01–1.08, p = 0.009), X-02269 (OR 1.04, 95% CI 1–1.08, p = 0.044), guanosine (OR 1.04, 95% CI 1–1.08, p = 0.035), X-12544 (OR 1.04, 95% CI 1–1.07, p < 0.05), isoleucine (OR 1.03, 95% CI 1–1.07, p = 0.038), 1-arachidonoylglycerophosphocholine (OR 1.02, 95% CI 1.01–1.04, p = 0.005), 2-linoleoylglycerophosphocholine (OR 1.02, 95% CI 1–1.04, p = 0.033), pyruvate (OR 1.02, 95% CI 1–1.04, p = 0.026) and 1-docosahexaenoylglycerophosphocholine (OR 1.02, 95% CI 1–1.03, p = 0.044); while decreasing of propionylcarnitine (OR 0.99, 95% CI 0.98–1, p = 0.038), isovalerate (OR 0.98, 95% CI 0.97–1, p = 0.036), X-13553 (OR 0.98, 95% CI 0.97–1, p = 0.025), beta-hydroxyisovalerate (OR 0.98, 95% CI 0.96–0.99, p = 0.003), X-12847 (OR 0.97, 95% CI 0.94–1, p = 0.043), aspartylphenylalanine (OR 0.96, 95% CI 0.94–0.99, p = 0.004) and ursodeoxycholate (OR 0.95, 95% CI 0.92–0.98, p = 0.001) (Figure 3A).
FIGURE 3. The effect of genetically determined pancreatic cancer on the alternation of peripheral metabolites in the IVW results. OR, odds ratios; CI, confidence interval; IVW, inverse variance-weighted.
The additional analyses using the maximum likelihood, weighted-median and MR Egger methods consistently yielded a comparable consistent trend of PC on peripheral metabolites to the IVW method (Supplementary Table S5). The MR-Egger results demonstrated that PC-related metabolites, including X-02269 (OR 0.39), guanosine (OR 0.83), X-12544 (OR 0.48), 2-linoleoylglycerophosphocholine (OR 0.57), 1-docosahexaenoylglycerophosphocholine (OR 0.99), propionylcarnitine (OR 1.58), X-13553 (OR 1.05), X-12847 (OR 1.27), aspartylphenylalanine (OR 3.08) and ursodeoxycholate (OR 3.98), were contrast to IVW results although neither of them reached significant. The MR-Egger intercept test suggested no evidence of pleiotropy and MR-PRESSO global test indicated the absence of horizontal pleiotropic (Supplementary Table S6). No significant heterogeneity was observed based on Cochran’s Q statistics when considering non-outlier single IVs within the PC cohort. Scatter plots present the individual SNP effect and combined effect from each method for each outcome dataset (Supplementary Figure S4). The leave-one-out analyses confirmed the stability of the results even after removing one SNP at a time (Supplementary Figure S5). Steiger filtering indicated that the selected genetic IVs for PC exhibited a higher degree of variance explanation in PC compared to peripheral metabolites (Supplementary Table S5).
We next compared the bi-causal relationships of PC and peripheral metabolites to find out the existence of positive feedback in PC development. Interestingly, X-02269 and PC formed a positive causal relationship pair (Figure 4), highlighting that X-02269 and PC might form a vicious circle in PC progression. We also found that 1-arachidonoylglycerophosphocholine indicated a lower risk of PC while PC demonstrated a higher level of it (Figure 2 and Figure 3). The result implied that 1-arachidonoylglycerophosphocholine might serve as a negative feedback node in PC development.
FIGURE 4. The Venn plot of positive (left) and negative (right) bidirectional two-sample Mendelian randomization IVW results. Met2PC represents the effect of peripheral metabolites on pancreatic cancer; PC2Met represents the effect of pancreatic cancer on peripheral metabolites.
4 Discussion
PC stands as an ongoing challenge to global health due to its aggressive behavior and bleak prognosis (Qin et al., 2020; Casolino et al., 2021). The systemic metabolic status might play a crucial role in both PC development and non-neoplastic organ function. Despite advancements in cancer research, the intricate interplay between PC and peripheral metabolites remains largely unexplored. To illuminate the causal relationships between PC and peripheral metabolite profiles, we applied a comprehensive bidirectional tow-sample MR analysis. In our investigation, we discerned a spectrum of 20 peripheral metabolites intricately linked with the risk of PC. Furthermore, our study unveiled dynamic alterations in the levels of 17 metabolites in response to PC, thereby offering promising avenues for both early diagnostic strategies and therapeutic interventions. Notably, leucine, exhibiting the highest OR for PC causation in our analysis, has been previously reported to stimulate the growth of pancreatic cancer (Otsuru et al., 2019). Conversely, studies have reported that the restriction of leucine intake can suppress PC cells growth through the induction of apoptosis (Kamphorst et al., 2015). The detailed mechanisms underlying leucine’s promotion effect on PC may involve the activation of mTOR signaling and the promotion of tumor-related inflammation (Liu et al., 2014; Otsuru et al., 2019).
Interestingly, our findings underscore the role of lipid metabolic dysregulation as a significant regulator in PC progression. HDL particles exhibited an augmented PC risk, while VLDL particles demonstrated the converse trend. Cholesterol and cholesterol esters (CEs) within large HDL emerged as indicators of heightened PC risk in our analysis. This observation is corroborated by studies suggesting that HDL-mediated cholesterol depletion from PC cells restrains tumor development both in vivo and in vitro (Oberle et al., 2022). Moreover, the overexpression of SR-B1, a physiologically pertinent HDL receptor, in numerous tumors, including PC, has been implicated in promoting tumor growth (Rajora and Zheng, 2016; Oberle et al., 2022). This occurred through its binding to lipoproteins on the surface of HDL and facilitating lipid delivery, as well as the uptake of HDL-CE into the cell via a pore formed by SR-B1 (Valacchi et al., 2011). In contrast, the causal relationship of VLDL particles with PC was reversed compared to HDL particles in our results. VLDL, originating in the liver, facilitates the transfer of triglycerides from the liver to peripheral tissues. During its transit in the bloodstream, triglycerides are progressively removed, leading to the conversion of VLDL to intermediate-density lipoprotein (IDL) and subsequently to LDL. Despite this, scant research has focused on VLDL as opposed to LDL in PC. Some studies have suggested that genetically elevated LDL-cholesterol levels are associated with pancreatic cancer, and inhibiting LDLR in PC cells considerably reduces cholesterol uptake and alters its distribution, leading to decreased cancer cell proliferation (Jung et al., 2021; Yin et al., 2022). These observations may underscore the potential pathological significance of the transformation from VLDL to LDL in PC development. In conclusion, these results highlighted the pronounced heterogeneity of the HDL and VLDL particle pool and their interaction with PC cells or non-malignant cells in the context of disease progression, warranting further in-depth exploration.
The interplay between PC and peripheral metabolites was further substantiated through reverse MR analysis. Seventeen metabolites displayed altered levels in response to PC. Notably, two BCAAs, valine and isoleucine, displayed heightened levels under PC influence, contrary to prior studies reporting BCAA inhibition in PC, even though serum BCAA levels were found to be elevated in PC rats compared to those with chronic pancreatitis (Fang et al., 2007; Li et al., 2021). Moreover, enhanced BCAA catabolism by BCAT2 plays a pivotal role in the development of PC, and this role is further exacerbated by KRAS mutations (Li et al., 2020). Another significant metabolite type elevated by PC was derivatives of glycerophospholipids, including 1-arachidonoylglycerophosphocholine, 2-linoleoylglycerophosphocholine, and 1-docosahexaenoylglycerophosphocholine. Targeting glycerophospholipid synthesis has been suggested as a potential strategy to enhance cytotoxic drug sensitivity (Kaoutari et al., 2021). Furthermore, HIF-1 has the capacity to promote lipid droplet accumulation and enhance cancer cell viability under hypoxic conditions by directly targeting AGPAT2, an enzyme crucially involved in the glycerophospholipid/triacylglycerol biosynthesis pathway (Triantafyllou et al., 2018). The notable activity of pyruvate kinase type M2, culminating in pyruvate synthesis in PC, aligns with our findings of increased peripheral pyruvate levels (Cameron et al., 2018).
By cross-referencing the results of the bidirectional two-sample MR analysis, we identified two shared metabolites—X-02269 and 1-arachidonoylglycerophosphocholine. The mutual causation observed between X-02269 and PC implies their potential involvement in a vicious circle in PC onset and progression. Unfortunately, X-02269 remains inadequately investigated as a novel metabolite. Additionally, 1-arachidonoylglycerophosphocholine indicated a reduced PC risk and an increase under PC influence, possibly serving as a negative feedback regulator. However, there is a lack of research on 1-arachidonoylglycerophosphocholine and its relationship with PC, as far as we know, indicating the need for further studies in this area.
Despite the comprehensive nature of our MR analysis in elucidating the causal links between peripheral metabolites and PC, several limitations warrant consideration. Firstly, our conclusions are solely derived from MR analysis, necessitating experimental validation. Secondly, while we focused on the relationship between PC and peripheral metabolites, the role of local metabolites within the PC tissue remains unexplored. Thirdly, although we established causal connections between PC and metabolites, the underlying mechanisms necessitate further investigation, particularly regarding specific cell types and signaling pathways. Lastly, our study’s confinement to a single cancer type and dataset prompts us to acknowledge the potential for intriguing revelations through systematic analyses across diverse cancer types and cohorts.
5 Conclusion
The bidirectional MR analysis uncovers the intricate interplay between peripheral metabolites and PC, offering novel insights into PC etiology. Our findings contribute to the broader understanding of PC biology and may guide the development of innovative strategies for early diagnosis and intervention. Further research is warranted to elucidate the precise molecular mechanisms underlying the observed causal relationships and to explore the clinical implications of the identified metabolites.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics statement
Ethical approval was not required for the study involving humans in accordance with the local legislation and institutional requirements. Written informed consent to participate in this study was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and the institutional requirements.
Author contributions
RS: Writing–review and editing, Formal Analysis, Investigation, Writing–original draft. HX: Formal Analysis, Investigation, Methodology, Writing–original draft. FL: Conceptualization, Data curation, Resources, Writing–original draft. BZ: Conceptualization, Data curation, Writing–original draft. ML: Conceptualization, Data curation, Resources, Writing–review and editing. XS: Conceptualization, 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 appreciation to all the researchers who generously shared the data utilized in this study. We also express our gratitude to MRC IEU OpenGWAS project and FinnGen database for their invaluable contribution in providing the data to the public.
Conflict of interest
The reviewer SS declared a shared parent affiliation with the authors to the handling editor at the time of review.
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/fmolb.2023.1279157/full#supplementary-material
References
Borges, M. C., Schmidt, A. F., Jefferis, B., Wannamethee, S. G., Lawlor, D. A., Kivimaki, M., et al. (2020). Circulating fatty acids and risk of coronary heart disease and stroke: individual participant data meta-analysis in up to 16 126 participants. J. Am. Heart Assoc. 9, e013131. doi:10.1161/JAHA.119.013131
Bowden, J., Davey, S. G., and Burgess, S. (2015). Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int. J. Epidemiol. 44, 512–525. doi:10.1093/ije/dyv080
Bowden, J., Davey, S. G., Haycock, P. C., and Burgess, S. (2016). Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet. Epidemiol. 40, 304–314. doi:10.1002/gepi.21965
Bowden, J., Del, G. M. F., Minelli, C., Zhao, Q., Lawlor, D. A., Sheehan, N. A., et al. (2019). Improving the accuracy of two-sample summary-data Mendelian randomization: moving beyond the NOME assumption. Int. J. Epidemiol. 48, 728–742. doi:10.1093/ije/dyy258
Burgess, S., Butterworth, A., and Thompson, S. G. (2013). Mendelian randomization analysis with multiple genetic variants using summarized data. Genet. Epidemiol. 37, 658–665. doi:10.1002/gepi.21758
Burgess, S., Dudbridge, F., and Thompson, S. G. (2016). Combining information on multiple instrumental variables in Mendelian randomization: comparison of allele score and summarized data methods. Stat. Med. 35, 1880–1906. doi:10.1002/sim.6835
Cameron, M. E., Yakovenko, A., and Trevino, J. G. (2018). Glucose and lactate transport in pancreatic cancer: glycolytic metabolism revisited. J. Oncol. 2018, 6214838. doi:10.1155/2018/6214838
Casolino, R., Paiella, S., Azzolina, D., Beer, P. A., Corbo, V., Lorenzoni, G., et al. (2021). Homologous recombination deficiency in pancreatic cancer: a systematic review and prevalence meta-analysis. J. Clin. Oncol. 39, 2617–2631. doi:10.1200/JCO.20.03238
Chen, X., Liu, F., Xue, Q., Weng, X., and Xu, F. (2021). Metastatic pancreatic cancer: mechanisms and detection (Review). Oncol. Rep. 46, 231. doi:10.3892/or.2021.8182
DePeaux, K., and Delgoffe, G. M. (2021). Metabolic barriers to cancer immunotherapy. Nat. Rev. Immunol. 21, 785–797. doi:10.1038/s41577-021-00541-y
Drijvers, J. M., Sharpe, A. H., and Haigis, M. C. (2020). The effects of age and systemic metabolism on anti-tumor T cell responses. Elife 9, e62420. doi:10.7554/eLife.62420
Elebo, N., Omoshoro-Jones, J., Fru, P. N., Devar, J., De Wet, V. Z. C., Vorster, B. C., et al. (2021). Serum metabolomic and lipoprotein profiling of pancreatic ductal adenocarcinoma patients of african ancestry. Metabolites 11, 663. doi:10.3390/metabo11100663
Encarnacion-Rosado, J., and Kimmelman, A. C. (2021). Harnessing metabolic dependencies in pancreatic cancers. Nat. Rev. Gastroenterol. Hepatol. 18, 482–492. doi:10.1038/s41575-021-00431-7
Fang, F., He, X., Deng, H., Chen, Q., Lu, J., Spraul, M., et al. (2007). Discrimination of metabolic profiles of pancreatic cancer from chronic pancreatitis by high-resolution magic angle spinning 1H nuclear magnetic resonance and principal components analysis. Cancer Sci. 98, 1678–1682. doi:10.1111/j.1349-7006.2007.00589.x
Garcia-Bermudez, J., Baudrier, L., La, K., Zhu, X. G., Fidelin, J., Sviderskiy, V. O., et al. (2018). Aspartate is a limiting metabolite for cancer cell proliferation under hypoxia and in tumours. Nat. Cell Biol. 20, 775–781. doi:10.1038/s41556-018-0118-z
Hanahan, D., and Weinberg, R. A. (2011). Hallmarks of cancer: the next generation. Cell 144, 646–674. doi:10.1016/j.cell.2011.02.013
Hemani, G., Tilling, K., and Davey, S. G. (2017). Orienting the causal relationship between imprecisely measured traits using GWAS summary data. Plos Genet. 13, e1007081. doi:10.1371/journal.pgen.1007081
Huang, L., Jansen, L., Balavarca, Y., Molina-Montes, E., Babaei, M., van der Geest, L., et al. (2019). Resection of pancreatic cancer in Europe and USA: an international large-scale study highlighting large variations. Gut 68, 130–139. doi:10.1136/gutjnl-2017-314828
Infantino, V., Santarsiero, A., Convertini, P., Todisco, S., and Iacobazzi, V. (2021). Cancer cell metabolism in hypoxia: role of HIF-1 as key regulator and therapeutic target. Int. J. Mol. Sci. 22, 5703. doi:10.3390/ijms22115703
Jung, Y. Y., Ko, J. H., Um, J. Y., Chinnathambi, A., Alharbi, S. A., Sethi, G., et al. (2021). LDL cholesterol promotes the proliferation of prostate and pancreatic cancer cells by activating the STAT3 pathway. J. Cell Physiol. 236, 5253–5264. doi:10.1002/jcp.30229
Kamphorst, J. J., Nofal, M., Commisso, C., Hackett, S. R., Lu, W., Grabocka, E., et al. (2015). Human pancreatic cancer tumors are nutrient poor and tumor cells actively scavenge extracellular protein. Cancer Res. 75, 544–553. doi:10.1158/0008-5472.CAN-14-2211
Kaoutari, A. E., Fraunhoffer, N. A., Hoare, O., Teyssedou, C., Soubeyran, P., Gayet, O., et al. (2021). Metabolomic profiling of pancreatic adenocarcinoma reveals key features driving clinical outcome and drug resistance. Ebiomedicine 66, 103332. doi:10.1016/j.ebiom.2021.103332
Kettunen, J., Demirkan, A., Wurtz, P., Draisma, H. H., Haller, T., Rawal, R., et al. (2016). Genome-wide study for circulating metabolites identifies 62 loci and reveals novel systemic effects of LPA. Nat. Commun. 7, 11122. doi:10.1038/ncomms11122
Khalaf, N., El-Serag, H. B., Abrams, H. R., and Thrift, A. P. (2021). Burden of pancreatic cancer: from Epidemiology to practice. Clin. Gastroenterol. Hepatol. 19, 876–884. doi:10.1016/j.cgh.2020.02.054
Leal-Esteban, L. C., and Fajas, L. (2020). Cell cycle regulators in cancer cell metabolism. Biochim. Biophys. Acta Mol. Basis Dis. 1866, 165715. doi:10.1016/j.bbadis.2020.165715
Li, J. T., Yin, M., Wang, D., Wang, J., Lei, M. Z., Zhang, Y., et al. (2020). BCAT2-mediated BCAA catabolism is critical for development of pancreatic ductal adenocarcinoma. Nat. Cell Biol. 22, 167–174. doi:10.1038/s41556-019-0455-6
Li, J. Y., Sun, F., Yang, C. L., Zhou, H. F., Gao, M., Zhang, Q., et al. (2021). GEO data mining and TCGA analysis reveal altered branched chain amino acid metabolism in pancreatic cancer patients. Aging (Albany NY) 13, 11907–11918. doi:10.18632/aging.202892
Liu, K. A., Lashinger, L. M., Rasmussen, A. J., and Hursting, S. D. (2014). Leucine supplementation differentially enhances pancreatic cancer growth in lean and overweight mice. Cancer Metab. 2, 6. doi:10.1186/2049-3002-2-6
Mayers, J. R., Wu, C., Clish, C. B., Kraft, P., Torrence, M. E., Fiske, B. P., et al. (2014). Elevation of circulating branched-chain amino acids is an early event in human pancreatic adenocarcinoma development. Nat. Med. 20, 1193–1198. doi:10.1038/nm.3686
Oberle, R., Kührer, K., Österreicher, T., Weber, F., Steinbauer, S., Udonta, F., et al. (2022). The HDL particle composition determines its antitumor activity in pancreatic cancer. Life Sci. Alliance 5, e202101317. doi:10.26508/lsa.202101317
Otsuru, T., Kobayashi, S., Wada, H., Takahashi, T., Gotoh, K., Iwagami, Y., et al. (2019). Epithelial-mesenchymal transition via transforming growth factor beta in pancreatic cancer is potentiated by the inflammatory glycoprotein leucine-rich alpha-2 glycoprotein. Cancer Sci. 110, 985–996. doi:10.1111/cas.13918
Pavlova, N. N., Zhu, J., and Thompson, C. B. (2022). The hallmarks of cancer metabolism: still emerging. Cell Metab. 34, 355–377. doi:10.1016/j.cmet.2022.01.007
Porporato, P. E. (2016). Understanding cachexia as a cancer metabolism syndrome. Oncogenesis 5, e200. doi:10.1038/oncsis.2016.3
Qin, C., Yang, G., Yang, J., Ren, B., Wang, H., Chen, G., et al. (2020). Metabolism of pancreatic cancer: paving the way to better anticancer strategies. Mol. Cancer 19, 50. doi:10.1186/s12943-020-01169-7
Rajora, M. A., and Zheng, G. (2016). Targeting SR-BI for cancer diagnostics, imaging and therapy. Front. Pharmacol. 7, 326. doi:10.3389/fphar.2016.00326
Schworer, S., Vardhana, S. A., and Thompson, C. B. (2019). Cancer metabolism drives a stromal regenerative response. Cell Metab. 29, 576–591. doi:10.1016/j.cmet.2019.01.015
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, 543–550. doi:10.1038/ng.2982
Shu, X., Zheng, W., Yu, D., Li, H. L., Lan, Q., Yang, G., et al. (2018). Prospective metabolomics study identifies potential novel blood metabolites associated with pancreatic cancer risk. Int. J. Cancer 143, 2161–2167. doi:10.1002/ijc.31574
Stolzenberg-Solomon, R., Derkach, A., Moore, S., Weinstein, S. J., Albanes, D., and Sampson, J. (2020). Associations between metabolites and pancreatic cancer risk in a large prospective epidemiological study. Gut 69, 2008–2015. doi:10.1136/gutjnl-2019-319811
Triantafyllou, E. A., Georgatsou, E., Mylonis, I., Simos, G., and Paraskeva, E. (2018). Expression of AGPAT2, an enzyme involved in the glycerophospholipid/triacylglycerol biosynthesis pathway, is directly regulated by HIF-1 and promotes survival and etoposide resistance of cancer cells under hypoxia. Biochim. Biophys. Acta Mol. Cell Biol. Lipids 1863, 1142–1152. doi:10.1016/j.bbalip.2018.06.015
Valacchi, G., Sticozzi, C., Lim, Y., and Pecorelli, A. (2011). Scavenger receptor class B type I: a multifunctional receptor. Ann. N. Y. Acad. Sci. 1229, E1–E7. doi:10.1111/j.1749-6632.2011.06205.x
Verbanck, M., Chen, C. Y., Neale, B., and Do, R. (2018). Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat. Genet. 50, 693–698. doi:10.1038/s41588-018-0099-7
Wang, S., Li, M., Yan, L., He, M., Lin, H., Xu, Y., et al. (2022). Metabolomics study reveals systematic metabolic dysregulation and early detection markers associated with incident pancreatic cancer. Int. J. Cancer 150, 1091–1100. doi:10.1002/ijc.33877
Keywords: pancreatic cancer, peripheral metabolites, two-sample Mendelian randomization, high-density lipoprotein, very low-density lipoprotein
Citation: Sun R, Xu H, Liu F, Zhou B, Li M and Sun X (2023) Unveiling the intricate causal nexus between pancreatic cancer and peripheral metabolites through a comprehensive bidirectional two-sample Mendelian randomization analysis. Front. Mol. Biosci. 10:1279157. doi: 10.3389/fmolb.2023.1279157
Received: 17 August 2023; Accepted: 16 October 2023;
Published: 25 October 2023.
Edited by:
Ioannis Ntanasis-Stathopoulos, National and Kapodistrian University of Athens Medical School, GreeceReviewed by:
Shiyu Song, Nanjing University, ChinaEleni Zografos, National and Kapodistrian University of Athens, Greece
Copyright © 2023 Sun, Xu, Liu, Zhou, Li and Sun. 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: Minli Li, bGltaW5saV94aEAxNjMuY29t; Xiangdong Sun, c3VueGRfODFAMTI2LmNvbQ==
†These authors have contributed equally to this work