• KAFM
  • Contact us
  • E-Submission
ABOUT
ARTICLE CATEGORY
BROWSE ARTICLES
AUTHOR INFORMATION

Articles

Original Article

Assessing the impact of metabolomic markers on gastric cancer risk: a two-sample Mendelian randomization study

Published online: January 14, 2026

1Faculty of Pharmacy, University of Health Sciences, Vietnam National University, Ho Chi Minh City, Vietnam

2Faculty of Odonto-Stomatology, University of Health Sciences, Vietnam National University, Ho Chi Minh City, Vietnam

3Department of Gastroenterology and Hepatology, Nghe An Oncology Hospital, Nghe An, Vietnam

4Department of Oncology, Hanoi Medical University, Ha Noi, Vietnam

*Corresponding Author: Tung Hoang Tel: +84-368-730-588, Fax: +84-368-730-588, E-mail: htung@uhsvnu.edu.vn
• Received: July 29, 2025   • Revised: August 23, 2025   • Accepted: September 3, 2025

© 2026 The Korean Academy of Family Medicine

This is an open-access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted noncommercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

  • 412 Views
  • 15 Download
  • Background
    This study aimed to examine the relationship between genetically predicted metabolite levels and gastric cancer (GC) risk using Mendelian randomization (MR), and to identify the metabolic pathways potentially involved.
  • Methods
    We selected genetic instruments for metabolites from 64 genome-wide association studies covering 362,750 participants. A two-sample MR design was applied to evaluate the associations with GC using summary-level data from a combined analysis of the UK Biobank and FinnGen. The primary analysis relied on the inverse-variance weighted method, while the median-weighted and MR-Egger methods were used to account for potential violations of instrumental variable assumptions and provide the estimate even when a subset of instruments was invalid. The MR-Egger intercept test was performed to detect directional pleiotropy. Metabolites showing significant associations with GC were further examined using pathway enrichment analysis to identify relevant metabolic and lipid processes.
  • Results
    MR analyses identified 25 and 17 metabolites that were positively and inversely associated with GC risk, respectively. Notably, hexanoylcarnitine and cis-4-decenoylcarnitine were strongly associated with increased risk, whereas pregnanediol disulfate, acetylcarnitine, prolyl-hydroxyproline, and X-18914 were associated with reduced risk, with no evidence of heterogeneity or directional pleiotropy. Enrichment analyses highlighted key metabolic pathways, including cysteine and methionine catabolism, beta-oxidation of pristanoyl-CoA (coenzyme A), oxidation of branched-chain fatty acids, and peroxisomal lipid metabolism.
  • Conclusion
    This study identified a set of genetically predicted metabolites associated with GC risk, highlighting the potential utility of metabolite panels and lipid-based biomarkers for risk stratification and early detection. However, further standardization and extensive validation are necessary prior to clinical application.
According to GLOBOCAN 2022, gastric cancer (GC) is the fifth most commonly diagnosed cancer worldwide, with 968,350 new cases reported by 2022 [1]. Globally, there is a high incidence of GC in East Asian countries, especially where national cancer screening programs have been implemented [2]. In high-risk areas such as Japan and Korea, endoscopy has been included in the GC screening recommendations and has been shown to reduce the risk of GC by 61% and 47%, respectively [3]. However, endoscopy remains expensive and infeasible in many areas because of its invasiveness and complexity. Additionally, serious complications such as bleeding, infection, and patient discomfort are important factors to consider when implementing endoscopy for the early detection of GC [4]. Thus, non-invasive modalities such as Helicobacter pylori eradication, serum pepsinogen, and stool antigen tests have been introduced as potential screening strategies in low- and intermediate-risk areas [3]. However, these modalities have relatively low sensitivity and specificity for GC detection [3]. Therefore, identifying reliable biomarkers for more accurate screening and diagnostic tests is of great interest for the early detection of GC.
Metabolites, byproducts, or intermediates of metabolic processes are closely linked to the phenotype of a biological system and play key roles in regulating their functional expression [5]. Evidence from the European Prospective Investigation into Cancer and Nutrition study, which included data from 10 cohorts (238 cases and 626 matched controls), found positive associations between certain fatty acids, including oleic acid, di-homo-γ-linolenic acid, and α-linolenic acid with GC. In a prospective analysis of 400 participants, cross-validated findings reported inverse associations between α-linolenic acid, linoleic acid, palmitic acid, arachidonic acid, sn-1 lysophosphatidylcholine (LysoPC)(18:3), and sn-2 LysoPC(20:3) with GC [6]. Incorporating these metabolites into predictive models enhances the accuracy of gastric lesion progression and early GC risk assessment beyond models that only include age, sex, H. pylori infection, and histopathological findings [6]. A recent study of 1,800 participants also identified a metabolic signature of 26 metabolites related to H. pylori, which was linked to an increased risk of GC [7]. However, large-scale human studies remain limited and observational studies cannot fully control for confounding factors and metabolic changes, thereby restricting causal inference.
Mendelian randomization (MR) uses genetic variants as instrumental variables to infer causal relationships between exposure and outcome, thereby minimizing confounding and avoiding reverse causation, which can affect observational studies [8]. By leveraging the random allocation of alleles at conception, MR can simulate certain aspects of a randomized controlled trial and reduce bias from unmeasured or residual confounders as well as underlying disease effects [9]. Three recent MR studies addressed the causal effects of circulating metabolites on GC risk using genome-wide association study (GWAS) data for 486 blood metabolites from the same source and GC outcome data from the UK Biobank [1012]. In these studies, three metabolic pathways related to GC were identified: linoleic acid, valine, leucine, isoleucine, and histidine metabolism. More recently, large-scale summary statistics have become available, comprising 1,400 metabolites and metabolite ratios from the UK Biobank and GC outcome data from the FinnGen study, with histidine metabolism emerging as the most significant pathway associated with GC risk [13].
Building on this evidence and the hypothesis that additional metabolites may be involved in the metabolic pathways linked to GC risk, we expanded the data sources for both metabolites and GC by incorporating multiple GWAS datasets. Accordingly, we employed an MR approach to systematically examine metabolites associated with GC and explore the underlying metabolic pathways involved.
Data sources
Exposure data (metabolite GWAS summary statistics) were sourced from the MetaboAnalyst database (https://www.mgwas.ca/mGWAS/upload/MGBrowseView.xhtml). The initial dataset contained 337,881 independent variants associated with over 4,000 metabolite levels across 64 GWAS (r2<0.01 within 10,000 kb) excluding the UK Biobank [14]. After filtering for missing allele information and zero standard errors, 191,606 variants remained. A meta-analysis then yielded 175,130 significant and strong instruments (F-statistics>10 and P<5×10−8) across 2,286 metabolites [14]. Harmonizing these metabolite summary statistics with GC summary statistics (32,342,963 variants from FinnGen and the UK Biobank) resulted in 139,891 instrumental variables for the final MR association between 2,123 metabolites and GC (Figure 1).
As this study analyzed secondary data from previous published studies with anonymized data, ethical review and informed consent were not required.
Statistical analysis
The inverse-variance weighted (IVW) method assuming balanced pleiotropy was applied for the primary analysis [15]. In this method, we utilized single-nucleotide polymorphism (SNP)-specific Wald estimates derived by dividing the SNP-outcome association by the SNP-exposure association. These estimates were then pooled using a random-effects method, and heterogeneity was assessed using I2 (%) and Cochran’s Q test. Sensitivity analyses were conducted using the weighted median method [16] and MR-Egger regression [17]. Whereas the IVW method provides the lowest variance estimate under the assumption that all instrumental variables are valid, the weighted median method offers a balance between efficiency and robustness when some instrumental variables may be invalid. In addition, we applied MR-Egger regression, which can yield consistent estimates even in the presence of directional pleiotropy; however, this method is less efficient than IWV [1517]. A series of complementary analytical strategies were applied to evaluate the potential influence of horizontal pleiotropy. The MR-Egger regression intercept was first estimated, with the statistical significance of the intercept differing from zero interpreted as evidence of directional pleiotropy [17]. All analyses were performed using the R statistical ver. 4.2.0 (The R Foundation) [18]. A false discovery rate correction was applied to control for multiple comparisons, with a significance threshold of Q<0.05.
To further explore the mechanisms linking metabolites and GC, we conducted enrichment analyses of GC-associated metabolic and lipid pathways using the web-based tool MetaboAnalyst ver. 6.0 (Wishart Research Group, University of Alberta) [19]. This platform integrates curated annotations from multiple databases (including Human Metabolome Database and Small Molecule Pathway Database) via the Relational Metabolic Pathway Database, ensuring broad coverage of metabolites and lipids. Its principle relies on over-representation analysis (hypergeometric testing) to determine whether significant metabolites are enriched in specific pathways more than expected by chance [19]. Compared to conventional single-database tools (e.g., Kyoto Encyclopedia of Genes and Genomes-only approaches), MetaboAnalyst provides a more comprehensive, cross-referenced metabolite library, greater flexibility in pathway definitions, and robust statistical frameworks, thereby increasing the reliability and interpretability of pathway-level insights into GC.
Metabolites positively associated with GC risk and underlying pathways
MR analysis revealed 183 genetically predicted metabolites associated with an increased risk of GC under the IVW method. Of these, 25 metabolites remained significant in the sensitivity analyses using median-weighted and Egger regression methods (Table 1, Supplement 1). Among them, hexanoylcarnitine showed robust associations (odds ratio [OR], 1.06; 95% confidence interval [CI], 1.04–1.07 for the IVW method; OR, 1.09; 95% CI, 1.07–1.11 for the weighted median method; OR, 1.04; 95% CI, 1.01–1.07 for the MR-Egger method) (Figures 2A, 3A). Similarly, cis-4-decenoylcarnitine was strongly associated (OR, 1.07; 95% CI, 1.05–1.09 for the IVW method; OR, 1.14; 95% CI, 1.11–1.18 for the weighted median method; OR, 1.09; 95% CI, 1.04–1.14 for the MR-Egger method) (Figures 2B, 3B). There was no evidence of heterogeneity (Cochrane’s Q test, all P>0.999) or directional pleiotropy (MR-Egger intercept, P=0.055, P=0.077, and P=0.427, respectively). Enrichment analysis emphasized the role of metabolic and lipid processes in GC (Figure 4), with cysteine and methionine catabolism being the most significant pathways (P=6.91×10−6).
Metabolites negatively associated with GC risk and underlying pathways
Using the IVW method, MR analysis identified 175 genetically predicted metabolites that were inversely associated with the risk under the IVW method. Of these, 17 metabolites remained significant in the sensitivity analyses using the median-weighted and Egger regression methods (Table 2, Supplement 2). Among them, pregnanediol disulfate (OR, 0.95; 95% CI, 0.92–0.98 for the IVW method; OR, 0.93; 95% CI, 0.89–0.97 for the weighted median method; OR, 0.89; 95% CI, 0.83–0.96 for the MR-Egger method), acetylcarnitine (OR, 0.48; 95% CI, 0.42–0.53 for the IVW method; OR, 0.46; 95% CI, 0.39–0.53 for the weighted median method; OR, 0.31; 95% CI, 0.17–0.58 for the MR-Egger method), prolyl-hydroxyproline (OR, 0.67; 95% CI, 0.53–0.84 for the IVW method; OR, 0.66; 95% CI, 0.51–0.87 for the weighted median method; OR, 0.66; 95% CI, 0.52–0.84 for the MR-Egger method), and X-18914 (OR, 0.92; 95% CI, 0.87–0.98 for the IVW method; OR, 0.92; 95% CI, 0.87–0.97 for the weighted median method; OR, 0.82; 95% CI, 0.72–0.94 for the MR-Egger method) consistently showed robust associations (Figure 2C–F, Figure 3D–F). Furthermore, there was no evidence of heterogeneity (Cochran’s Q test, all P>0.90) or directional pleiotropy (MR-Egger intercept, P=0.096, P=0.184, P=0.790, and P=0.064, respectively). Enrichment analysis also emphasized the role of metabolic and lipid processes in GC (Figure 5), with beta-oxidation of pristanoyl-CoA (coenzyme A) (P=6.2×10−5), oxidation of branched-chain fatty acids (9.01×10−5), and peroxisomal lipid metabolism (P=0.001) emerging as the top three significant pathways.
This study found that specific genetically predicted metabolites were associated with GC risk. Among these, 25 metabolites, especially hexanoylcarnitine and cis-4-decenoylcarnitine, were associated with a higher risk, whereas 17 metabolites, especially pregnanediol disulfate, acetylcarnitine, prolyl-hydroxyproline, and X-18914, were associated with a lower risk. The key potentially relevant metabolic pathways include cysteine and methionine catabolism, beta-oxidation of pristanoyl-CoA, oxidation of branched-chain fatty acids, and peroxisomal lipid metabolism, suggesting that alterations in these processes may be associated with GC pathogenesis.
To date, several omics studies have focused on the potential of DNA, RNA, and protein biomarkers to explore the characteristics of GC [20]. In our current review, the predictive performance of metabolites in comparisons with protein biomarkers was elucidated by Jung et al. [21] and Matsumoto et al. [22] only. These findings are in line with existing knowledge regarding the low sensitivity of carcinoembryonic antigen (9.5% vs. 95.9% for urine metabolic profiling and 10.5% vs. 69.0% for plasma metabolic profiling) and carbohydrate antigen 19-9 (12.7% vs. 95.9% for urine metabolic profiling and 10.5% and 2.86% vs. 69.0% for plasma metabolic profiling) biomarkers [21,22]. Therefore, it is suggested that the metabolomics approach has a better predictive performance than the protein biomarker approach.
Metabolic reprogramming, particularly through the Warburg effect, plays a critical role in GC by promoting glucose uptake, glycolysis, and lactate production, even under aerobic conditions [23]. This shift not only supports tumor growth, invasion, and drug resistance but also disrupts mitochondrial function and alters lipid and amino acid metabolism, contributing to poor prognosis [23]. Given that prompt prognostic surveillance may contribute to better clinical outcomes in patients with GC, and that prognosis prediction is affected by physicians’ judgment based on various clinical evaluations, such as TNM stages and histopathology, a recent systematic review of seven studies reported the predictive value of metabolites in identifying GC metastasis [24]. Accordingly, metabolites are classified into the main pathogenic pathways, including energy production, structural protein synthesis, pro-apoptotic effects, cell cycle arrest, antioxidant capacity, tumor angiogenesis, enhanced degradation of collagen extracellular matrix, mitochondrial enzyme impairment, T-cell dysfunction/inactivation, and cell viability, migration, and invasion [24].
In the present study, we determined the oxidation processes related to GC risk, notably beta-oxidation of pristanoyl-CoA, oxidation of branched-chain fatty acids, and the pathways associated with cysteine and methionine catabolism and peroxisomal lipid metabolism. These findings highlight the metabolic flexibility of GC cells, which reprogram the oxidative pathways to meet their increased energy and biosynthetic demands. Betaoxidation of pristanoyl-CoA, a peroxisomal process that breaks down branched-chain fatty acids like pristanic acid, provides an alternative energy source that supports tumor cell survival and growth, particularly under metabolic stress [25]. In addition, aberrant branched-chain fatty acid oxidation can alter cellular energy balance, redox homeostasis, and lipid signaling, thereby influencing tumor growth and progression [26,27]. Disrupted branched-chain fatty acid catabolism may also modulate oncogenic pathways such as PI3K/AKT/mTOR and AMPK/mTOR [26,27], both of which are critical in GC biology. This finding aligns with the inverse associations between total branched-chain amino acid levels and GC from the UK Biobank [28]. Besides, alterations in amino acid and lipid metabolism are central to GC biology, with methionine and cysteine catabolism influencing one-carbon metabolism, redox balance, and sulfur metabolite production, while peroxisomal lipid metabolism sustains tumor growth through fatty acid oxidation and ether lipid biosynthesis [25]. Methionine restriction suppresses proliferation, enhances antitumor immunity, and sensitizes tumors to therapy, whereas peroxisomal activity promotes survival by balancing oxidative stress, energy demands, and mitochondrial function primarily via PPAR signaling [25].
Previous MR analyses explored the associations between 486 genetically predicted circulating metabolites and GC using summary statistics from the UK Biobank (1,029 cases and 474,841 controls) or the FinnGen cohort (1,307 GC cases and 287,137 controls) [1012]. These analyses highlighted three metabolic pathways, linoleic acid metabolism, valine/leucine/isoleucine metabolism, and histidine metabolism, and suggested potential protective roles for specific metabolites, such as 3-methyl-2-oxovalerate, piperine, and the dipeptide Phe-Phe, although the findings were inconsistent across datasets, with additional inverse associations reported for tryptophan, non-adecanoate, and erythritol. To overcome the limitations of earlier analyses that relied on single GWAS sources [29], our study integrated summary statistics from 64 metabolite GWAS and employed GC data obtained via meta-analysis of the UK Biobank and FinnGen cohorts, thereby increasing statistical power and generalizability. However, the significant metabolites and pathways identified in our study differed from those reported in previous studies, thereby offering novel insights into the metabolic basis of GC risk.
Although the increased sample size enhanced the statistical power, heterogeneity in metabolite measurements across different metabolomic platforms may introduce variability and complicate the harmonization of exposure data. In addition, discrepancies in summary statistics between FinnGen and UK Biobank, arising from differences in population structure, case definitions, and analytical pipelines, could contribute to inconsistencies in MR estimates. Another potential limitation of this study was the absence of bidirectional MR analyses. While metabolites can indeed be influenced by disease states, such as GC, our analyses were restricted to testing the causal direction from metabolites to GC risk. This was primarily due to data availability; we obtained only SNP-level associations for significant variants from the metabolite GWAS meta-analysis (64 studies) rather than complete summary statistics for each metabolite, which are required to examine the effect of GC on metabolite levels. Therefore, we cannot completely exclude the possibility of reverse causality. Future studies with publicly available full GWAS summary statistics for metabolites are essential to validate the causal relationships between metabolites and GC. Furthermore, an important limitation of our study is that both the exposure and outcome datasets were primarily derived from individuals of European ancestry. This may restrict the generalizability of our findings to populations with different genetic backgrounds and etiologies of GC. In particular, gastric adenocarcinoma is more prevalent in East Asian populations, such as Koreans and Japanese, where H. pylori infection, dietary patterns, and genetic susceptibility factors contribute more prominently to disease risk than in European populations [30]. Given the epidemiological and etiological differences in GC between European and Asian populations, replication in Asian cohorts is a critical future step to validate the causal metabolites identified in this study and better understand population-specific disease mechanisms.
In summary, this study found that genetically predicted metabolites are associated with GC risk and may serve as potential biomarkers for early detection, with 25 metabolites being linked to an increased risk and 17 to a reduced risk. These findings underscore the possible role of the disruption of amino acid metabolism and other key metabolic pathways in GC pathogenesis.

Conflict of interest

No potential conflict of interest relevant to this article was reported.

Funding

Tho Thi Anh Tran was funded by the Master, PhD Scholarship Programme of Vingroup Innovation Foundation (VINIF), code VINIF.2023.TS.119. This research was also funded by Vietnam National University Ho Chi Minh City (VNU-HCM) under the grant number C2024-44-34.

Data availability

All data for this study are available upon reasonable request to the corresponding author.

Author contribution

Conceptualization: TH, VMT, TTAT. Data curation: TH. Formal analysis: TH. Methodology: TH, VMT, TTAT. Writing–original draft: TH. Writing–review & editing: VMT, TTAT. Final approval of the manuscript: all authors.

Supplementary materials can be found via https://doi.org/10.4082/kjfm.25.0229.
Supplement 1.
Sensitivity analysis identified metabolites positively associated with gastric cancer risk.
kjfm-25-0229-Supplement-1.pdf
Supplement 2.
Sensitivity analysis identified metabolites negatively associated with gastric cancer risk.
kjfm-25-0229-Supplement-2.pdf
Figure 1
Flowchart for identification of instrumental variables.
kjfm-25-0229f1.jpg
Figure 2
Funnel plots for Mendelian randomization (MR) analysis of hexanoylcarnitine (A), cis-4-decenoylcarnitine (B), pregnanediol disulfate (C), acetylcarnitine (D), prolyl-hydroxyproline (E), and X-18914 metabolites and gastric cancer (F). The horizontal axis represents the estimated effect of each single-nucleotide polymorphism on the exposure variable and the vertical axis represents the precision or uncertainty of these estimates. IVW, inverse-variance weighted; SE, standard error.
kjfm-25-0229f2.jpg
Figure 3
Scatter plots for Mendelian randomization (MR) analysis of hexanoylcarnitine (A), cis-4-decenoylcarnitine (B), pregnanediol disulfate (C), acetylcarnitine (D), prolyl-hydroxyproline (E), and X-18914 metabolites and gastric cancer (F). The horizontal axis denotes the impact of instrumental variables on metabolites, and the vertical axis represents the effect of instrumental variables on gastric cancer. Each black dot represents an individual single-nucleotide polymorphism (SNP), and the vertical and horizontal lines represent the corresponding 95% confidence intervals. The slope of the line represents the estimated causal effect of the various MR methods. IVW, inverse-variance weighted.
kjfm-25-0229f3.jpg
Figure 4
Enriched pathways of metabolites positively associated with gastric cancer. IMD, inherited metabolic disorder.
kjfm-25-0229f4.jpg
Figure 5
Enriched pathways of metabolites negatively associated with gastric cancer. CoA, coenzyme A; ESR, estrogen receptor.
kjfm-25-0229f5.jpg
kjfm-25-0229f6.jpg
Table 1
Genetically predicted metabolites positively associated with gastric cancer risk
Metabolite No. of instruments Median F-statistics (range) Inverse-variance weighted Heterogeneity
OR (95% CI) P-value Q-value I2 (%) P-value
Urate/histidine 322 52.04 (30.11–67.40) 2.44 (2.00–2.97) 9.87E-19 3.74E-17 0 >0.999
Decanoylcarnitine 296 39.53 (29.73–116.51) 1.04 (1.02–1.07) 8.13E-05 0.001 0 >0.999
Hexanoylcarnitine (fumarylcarnitine) 423 46.85 (29.73–67.40) 1.06 (1.04–1.07) 1.5E-15 4.62E-14 0 >0.999
cis-4-decenoyl carnitine 284 45.12 (29.91–67.40) 1.07 (1.05–1.09) 2.44E-10 4.79E-09 0 >0.999
Methyl glucopyranoside (α+β) 147 44.23 (29.76–67.40) 1.02 (1.01–1.04) 0.001 0.007 0 0.999
X-23314 98 32.49 (29.85–67.40) 1.04 (1.02–1.06) 0.001 0.004 0 0.999
Apolipoprotein B 1,099 35.06 (29.72–67.40) 1.27 (1.24–1.30) 3.38E-88 1.44E-85 0 >0.999
LDL-cholesterol 1,207 41.11 (29.73–219.46) 1.22 (1.20–1.25) 8.92E-78 1.89E-75 0 >0.999
Large LDL-cholesterol 1,166 39.37 (29.72–111.42) 1.21 (1.19–1.24) 6.38E-66 9.04E-64 0 >0.999
Large LDL cholesteryl esters 1,070 40.87 (29.72–108.62) 1.23 (1.21–1.26) 6.83E-78 1.61E-75 0 >0.999
Medium LDL-cholesterol 1,090 38.00 (29.72–108.96) 1.26 (1.23–1.29) 1.45E-88 7.71E-86 0 >0.999
Glycine 1,540 40.23 (29.72–67.40) 1.05 (1.03–1.07) 1.06E-09 1.99E-08 0 0.999
VLDL-D 604 46.98 (29.78–110.53) 1.12 (1.09–1.15) 1.26E-18 4.68E-17 0 >0.999
1-Ribosyl-imidazoleacetate 65 57.89 (42.74–67.40) 1.23 (1.16–1.31) 2.14E-11 4.78E-10 0 0.919
X-08402/Cholesterol 84 42.63 (29.72–67.40) 2.43 (1.61–3.66) 2.25E-05 0.0002 0 >0.999
Lysine 302 43.11 (30.09–77.15) 1.31 (1.26–1.37) 5.73E-42 4.68E-40 7.7 0.153
Pyruvate 69 33.16 (29.88–46.96) 1.26 (1.13–1.41) 2.6E-05 0.0003 0 >0.999
N1-methyl-3-pyridone-4-carboxamide 23 43.57 (31.79–66.03) 2.55 (1.68–3.86) 9.6E-06 0.0001 0 0.949
trans-urocanate 77 48.92 (29.98–66.03) 1.06 (1.03–1.09) 7.71E-05 0.001 0 0.772
Butyrylcarnitine 443 39.96 (29.76–80.31) 1.03 (1.02–1.05) 0.0002 0.002 0 0.512
Sphingomyelin 154 38.02 (29.81–67.40) 1.15 (1.09–1.21) 8.48E-08 1.34E-06 0.028 0.385
X-11905 103 35.22 (29.76–66.03) 1.46 (1.23–1.74) 2.11E-05 0.0002 0 >0.999
Dimethylglycine 116 46.72 (30.02–124.98) 1.07 (1.05–1.09) 3.68E-10 7.11E-09 0 >0.999
X-21849 245 30.29 (29.78–43.89) 1.07 (1.05–1.08) 3.99E-15 1.18E-13 0 >0.999
Acetylcarnitine/hexanoylcarnitine 80 49.16 (31.14–67.40) 1.51 (1.14–2.01) 0.004267 0.027 0 >0.999

OR, odds ratio; CI, confidence interval; LDL, low-density lipoprotein; VLDL, very low-density lipoprotein.

Table 2
Genetically metabolites negatively associated with gastric cancer risk
Metabolite No. of instruments Median F-statistics (range) Inverse-variance weighted Heterogeneity
OR (95% CI) P-value Q-value I2 (%) P-value
Glycylphenylalanine 92 44.00 (29.73–67.40) 0.97 (0.94–0.99) 0.002 0.014 0 >0.999
Alanine/tyrosine 290 47.64 (41.40–66.03) 0.88 (0.86–0.91) 3.1E-19 1.24E-17 0 >0.999
1-Arachidonoyl-glycosylphosphatidylinositol (20:4) 447 33.75 (30.12–65.23) 0.89 (0.86–0.92) 3.24E-12 7.74E-11 0 >0.999
X-11444 195 36.04 (29.80–67.40) 0.85 (0.83–0.88) 5.8E-31 3.42E-29 0.011 0.443
X-17178 102 43.31 (31.02–63.87) 0.88 (0.86–0.91) 3.83E-16 1.23E-14 0.352 0.0004
Pregnanediol disulfate 62 41.21 (30.16–64.22) 0.95 (0.92–0.98) 0.001 0.008 0 >0.999
X-11470 35 35.03 (29.72–64.22) 0.85 (0.80–0.90) 1.31E-08 2.31E-07 0 >0.999
Acetylcarnitine 59 34.87 (30.94–59.69) 0.48 (0.42–0.53) 3.54E-37 2.51E-35 0 >0.999
Propionylcarnitine 58 36.63 (29.95–63.30) 0.55 (0.50–0.61) 8.2E-33 5.12E-31 0 0.976
Imidazole lactate 138 38.45 (29.90–51.98) 0.98 (0.97–0.99) 0.001 0.008 0 >0.999
Glycocholenate sulfate 150 52.19 (30.96–67.40) 0.80 (0.75–0.84) 4.6E-16 1.46E-14 0.232 0.008
Phenylalanylserine 57 33.80 (29.81–67.40) 0.92 (0.89–0.94) 1.97E-09 3.61E-08 0 0.779
Glycolithocholate sulfate 65 53.70 (30.64–55.53) 0.74 (0.68–0.80) 3.63E-12 8.43E-11 0.524 6.28E-07
Estriol 238 43.86 (30.15–62.86) 0.98 (0.97–0.99) 3.77E-06 4.57E-05 0 >0.999
X-18914 47 42.91 (29.80–64.66) 0.92 (0.89–0.96) 3.05E-05 0.0003 0 0.974
X-12822 167 56.20 (29.88–66.03) 0.96 (0.94–0.98) 2.47E-06 3.15E-05 0 0.983
Prolyl-hydroxyproline 4 37.54 (29.91–38.77) 0.67 (0.53–0.84) 0.001 0.004 0 0.904

OR, odds ratio; CI, confidence interval.

  • 1. Bray F, Laversanne M, Sung H, Ferlay J, Siegel RL, Soerjomataram I, et al. Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin 2024;74:229-63.
  • 2. Wong MC, Huang J, Chan PS, Choi P, Lao XQ, Chan SM, et al. Global incidence and mortality of gastric cancer, 1980–2018. JAMA Netw Open 2021;4:e2118457.
  • 3. Januszewicz W, Turkot MH, Malfertheiner P, Regula J. A global perspective on gastric cancer screening: which concepts are feasible, and when? Cancers (Basel) 2023;15:664.
  • 4. Waddingham W, Kamran U, Kumar B, Trudgill NJ, Tsiamoulos ZP, Banks M. Complications of diagnostic upper Gastrointestinal endoscopy: common and rare: recognition, assessment and management. BMJ Open Gastroenterol 2022;9:e000688.
  • 5. Qiu S, Cai Y, Yao H, Lin C, Xie Y, Tang S, et al. Small molecule metabolites: discovery of biomarkers and therapeutic targets. Signal Transduct Target Ther 2023;8:132.
  • 6. Huang S, Guo Y, Li ZW, Shui G, Tian H, Li BW, et al. Identification and validation of plasma metabolomic signatures in precancerous gastric lesions that progress to cancer. JAMA Netw Open 2021;4:e2114186.
  • 7. Yu Y, Zheng Z, Gao X, Gu Y, Zhang M, Hu B, et al. Plasma metabolomic signatures of H. pylori infection, alcohol drinking, smoking, and risk of gastric cancer. Mol Carcinog 2025;64:463-74.
  • 8. Hu X, Cai M, Xiao J, Wan X, Wang Z, Zhao H, et al. Benchmarking Mendelian randomization methods for causal inference using genome-wide association study summary statistics. Am J Hum Genet 2024;111:1717-35.
  • 9. Sobczyk MK, Zheng J, Davey Smith G, Gaunt TR. Systematic comparison of Mendelian randomisation studies and randomised controlled trials using electronic databases. BMJ Open 2023;13:e072087.
  • 10. Chen S, Zhang B, Wang S, Yang M, Shen Q, Zhang R, et al. Blood metabolites and gastric cancer risk: a bidirectional 2-sample Mendelian randomization study. Medicine (Baltimore) 2025;104:e43339.
  • 11. Zhang C, Huang DL, Zhou K, Cai JT, Liu D, Tan MH, et al. Human blood metabolites and gastric cancer: a Mendelian randomization analysis. BMC Gastroenterol 2024;24:478.
  • 12. Qian L, Song J, Zhang X, Qiao Y, Tan Z, Li S, et al. Elucidating the causal relationship between 486 genetically predicted blood metabolites and the risk of gastric cancer: a comprehensive Mendelian randomization analysis. Front Oncol 2024;14:1418283.
  • 13. Yu H, Liu Y, Xie P, Feng W, Xu D. Genetic evidence of circulating metabolites linking to gastric cancers: a two-sample Mendelian randomization study. Ann Med Surg 2025;87:571-6.
  • 14. Burgess S, Davey Smith G, Davies NM, Dudbridge F, Gill D, Glymour MM, et al. Guidelines for performing Mendelian randomization investigations: update for summer 2023. Wellcome Open Res 2019;4:186.
  • 15. Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan N, Thompson J. A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization. Stat Med 2017;36:1783-802.
  • 16. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol 2016;40:304-14.
  • 17. Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol 2015;44:512-25.
  • 18. R Foundation for Statistical Computing. A language and environment for statistical computing [Internet] R Foundation for Statistical Computing; 2024 [cited 2025 Nov 10]. Available from: https://www.R-project.org/
  • 19. Pang Z, Lu Y, Zhou G, Hui F, Xu L, Viau C, et al. MetaboAnalyst 6.0: towards a unified platform for metabolomics data processing, analysis and interpretation. Nucleic Acids Res 2024;52(W1):W398-406.
  • 20. Chen Y, Wang B, Zhao Y, Shao X, Wang M, Ma F, et al. Metabolomic machine learning predictor for diagnosis and prognosis of gastric cancer. Nat Commun 2024;15:1657.
  • 21. Jung J, Jung Y, Bang EJ, Cho SI, Jang YJ, Kwak JM, et al. Noninvasive diagnosis and evaluation of curative surgery for gastric cancer by using NMR-based metabolomic profiling. Ann Surg Oncol 2014;21(Suppl 4):S736-42.
  • 22. Matsumoto T, Fukuzawa M, Itoi T, Sugimoto M, Aizawa Y, Sunamura M, et al. Targeted metabolomic profiling of plasma samples in gastric cancer by liquid chromatography-mass spectrometry. Digestion 2023;104:97-108.
  • 23. Liu Y, Zhang Z, Wang J, Chen C, Tang X, Zhu J, et al. Metabolic reprogramming results in abnormal glycolysis in gastric cancer: a review. Onco Targets Ther 2019;12:1195-204.
  • 24. Ursu S, Ciocan A, Ursu CP, Gherman CD, Ciocan RA, Pop RS, et al. Role of metabolomics in pathogenesis and prompt diagnosis of gastric cancer metastasis: a systematic review. Diagnostics (Basel) 2023;13:3401.
  • 25. Pratama AM, Sharma M, Naidu S, Bommel H, Prabhuswamimath SC, Madhusudhan T, et al. Peroxisomes and PPARs: emerging role as master regulators of cancer metabolism. Mol Metab 2024;90:102044.
  • 26. Xiong H, Liu R, Xu K, Chen X, Huang L, Shou Y, et al. Branched-chain amino acid and cancer: metabolism, immune microenvironment and therapeutic targets. J Transl Med 2025;23:636.
  • 27. Xu E, Ji B, Jin K, Chen Y. Branched-chain amino acids catabolism and cancer progression: focus on therapeutic interventions. Front Oncol 2023;13:1220638.
  • 28. Yu L, Bao S, Zhu F, Xu Y, Liu Y, Jiang R, et al. Association between branched-chain amino acid levels and gastric cancer risk: large-scale prospective cohort study. Front Nutr 2024;11:1479800.
  • 29. Shin SY, Fauman EB, Petersen AK, Krumsiek J, Santos R, Huang J, et al. An atlas of genetic influences on human blood metabolites. Nat Genet 2014;46:543-50.
  • 30. Shin WS, Xie F, Chen B, Yu P, Yu J, To KF, et al. Updated epidemiology of gastric cancer in Asia: decreased incidence but still a big challenge. Cancers (Basel) 2023;15:2639.

Figure & Data

References

    Citations

    Citations to this article as recorded by  

      Download Citation

      Download a citation file in RIS format that can be imported by all major citation management software, including EndNote, ProCite, RefWorks, and Reference Manager.

      Format:

      Include:

      Assessing the impact of metabolomic markers on gastric cancer risk: a two-sample Mendelian randomization study
      Download Citation
      Download a citation file in RIS format that can be imported by all major citation management software, including EndNote, ProCite, RefWorks, and Reference Manager.

      Format:
      • RIS — For EndNote, ProCite, RefWorks, and most other reference management software
      • BibTeX — For JabRef, BibDesk, and other BibTeX-specific software
      Include:
      • Citation for the content below
      Assessing the impact of metabolomic markers on gastric cancer risk: a two-sample Mendelian randomization study
      Close

      Figure

      • 0
      • 1
      • 2
      • 3
      • 4
      • 5
      Assessing the impact of metabolomic markers on gastric cancer risk: a two-sample Mendelian randomization study
      Image Image Image Image Image Image
      Figure 1 Flowchart for identification of instrumental variables.
      Figure 2 Funnel plots for Mendelian randomization (MR) analysis of hexanoylcarnitine (A), cis-4-decenoylcarnitine (B), pregnanediol disulfate (C), acetylcarnitine (D), prolyl-hydroxyproline (E), and X-18914 metabolites and gastric cancer (F). The horizontal axis represents the estimated effect of each single-nucleotide polymorphism on the exposure variable and the vertical axis represents the precision or uncertainty of these estimates. IVW, inverse-variance weighted; SE, standard error.
      Figure 3 Scatter plots for Mendelian randomization (MR) analysis of hexanoylcarnitine (A), cis-4-decenoylcarnitine (B), pregnanediol disulfate (C), acetylcarnitine (D), prolyl-hydroxyproline (E), and X-18914 metabolites and gastric cancer (F). The horizontal axis denotes the impact of instrumental variables on metabolites, and the vertical axis represents the effect of instrumental variables on gastric cancer. Each black dot represents an individual single-nucleotide polymorphism (SNP), and the vertical and horizontal lines represent the corresponding 95% confidence intervals. The slope of the line represents the estimated causal effect of the various MR methods. IVW, inverse-variance weighted.
      Figure 4 Enriched pathways of metabolites positively associated with gastric cancer. IMD, inherited metabolic disorder.
      Figure 5 Enriched pathways of metabolites negatively associated with gastric cancer. CoA, coenzyme A; ESR, estrogen receptor.
      Graphical abstract
      Assessing the impact of metabolomic markers on gastric cancer risk: a two-sample Mendelian randomization study

      Genetically predicted metabolites positively associated with gastric cancer risk

      Metabolite No. of instruments Median F-statistics (range) Inverse-variance weighted Heterogeneity
      OR (95% CI) P-value Q-value I2 (%) P-value
      Urate/histidine 322 52.04 (30.11–67.40) 2.44 (2.00–2.97) 9.87E-19 3.74E-17 0 >0.999
      Decanoylcarnitine 296 39.53 (29.73–116.51) 1.04 (1.02–1.07) 8.13E-05 0.001 0 >0.999
      Hexanoylcarnitine (fumarylcarnitine) 423 46.85 (29.73–67.40) 1.06 (1.04–1.07) 1.5E-15 4.62E-14 0 >0.999
      cis-4-decenoyl carnitine 284 45.12 (29.91–67.40) 1.07 (1.05–1.09) 2.44E-10 4.79E-09 0 >0.999
      Methyl glucopyranoside (α+β) 147 44.23 (29.76–67.40) 1.02 (1.01–1.04) 0.001 0.007 0 0.999
      X-23314 98 32.49 (29.85–67.40) 1.04 (1.02–1.06) 0.001 0.004 0 0.999
      Apolipoprotein B 1,099 35.06 (29.72–67.40) 1.27 (1.24–1.30) 3.38E-88 1.44E-85 0 >0.999
      LDL-cholesterol 1,207 41.11 (29.73–219.46) 1.22 (1.20–1.25) 8.92E-78 1.89E-75 0 >0.999
      Large LDL-cholesterol 1,166 39.37 (29.72–111.42) 1.21 (1.19–1.24) 6.38E-66 9.04E-64 0 >0.999
      Large LDL cholesteryl esters 1,070 40.87 (29.72–108.62) 1.23 (1.21–1.26) 6.83E-78 1.61E-75 0 >0.999
      Medium LDL-cholesterol 1,090 38.00 (29.72–108.96) 1.26 (1.23–1.29) 1.45E-88 7.71E-86 0 >0.999
      Glycine 1,540 40.23 (29.72–67.40) 1.05 (1.03–1.07) 1.06E-09 1.99E-08 0 0.999
      VLDL-D 604 46.98 (29.78–110.53) 1.12 (1.09–1.15) 1.26E-18 4.68E-17 0 >0.999
      1-Ribosyl-imidazoleacetate 65 57.89 (42.74–67.40) 1.23 (1.16–1.31) 2.14E-11 4.78E-10 0 0.919
      X-08402/Cholesterol 84 42.63 (29.72–67.40) 2.43 (1.61–3.66) 2.25E-05 0.0002 0 >0.999
      Lysine 302 43.11 (30.09–77.15) 1.31 (1.26–1.37) 5.73E-42 4.68E-40 7.7 0.153
      Pyruvate 69 33.16 (29.88–46.96) 1.26 (1.13–1.41) 2.6E-05 0.0003 0 >0.999
      N1-methyl-3-pyridone-4-carboxamide 23 43.57 (31.79–66.03) 2.55 (1.68–3.86) 9.6E-06 0.0001 0 0.949
      trans-urocanate 77 48.92 (29.98–66.03) 1.06 (1.03–1.09) 7.71E-05 0.001 0 0.772
      Butyrylcarnitine 443 39.96 (29.76–80.31) 1.03 (1.02–1.05) 0.0002 0.002 0 0.512
      Sphingomyelin 154 38.02 (29.81–67.40) 1.15 (1.09–1.21) 8.48E-08 1.34E-06 0.028 0.385
      X-11905 103 35.22 (29.76–66.03) 1.46 (1.23–1.74) 2.11E-05 0.0002 0 >0.999
      Dimethylglycine 116 46.72 (30.02–124.98) 1.07 (1.05–1.09) 3.68E-10 7.11E-09 0 >0.999
      X-21849 245 30.29 (29.78–43.89) 1.07 (1.05–1.08) 3.99E-15 1.18E-13 0 >0.999
      Acetylcarnitine/hexanoylcarnitine 80 49.16 (31.14–67.40) 1.51 (1.14–2.01) 0.004267 0.027 0 >0.999

      OR, odds ratio; CI, confidence interval; LDL, low-density lipoprotein; VLDL, very low-density lipoprotein.

      Genetically metabolites negatively associated with gastric cancer risk

      Metabolite No. of instruments Median F-statistics (range) Inverse-variance weighted Heterogeneity
      OR (95% CI) P-value Q-value I2 (%) P-value
      Glycylphenylalanine 92 44.00 (29.73–67.40) 0.97 (0.94–0.99) 0.002 0.014 0 >0.999
      Alanine/tyrosine 290 47.64 (41.40–66.03) 0.88 (0.86–0.91) 3.1E-19 1.24E-17 0 >0.999
      1-Arachidonoyl-glycosylphosphatidylinositol (20:4) 447 33.75 (30.12–65.23) 0.89 (0.86–0.92) 3.24E-12 7.74E-11 0 >0.999
      X-11444 195 36.04 (29.80–67.40) 0.85 (0.83–0.88) 5.8E-31 3.42E-29 0.011 0.443
      X-17178 102 43.31 (31.02–63.87) 0.88 (0.86–0.91) 3.83E-16 1.23E-14 0.352 0.0004
      Pregnanediol disulfate 62 41.21 (30.16–64.22) 0.95 (0.92–0.98) 0.001 0.008 0 >0.999
      X-11470 35 35.03 (29.72–64.22) 0.85 (0.80–0.90) 1.31E-08 2.31E-07 0 >0.999
      Acetylcarnitine 59 34.87 (30.94–59.69) 0.48 (0.42–0.53) 3.54E-37 2.51E-35 0 >0.999
      Propionylcarnitine 58 36.63 (29.95–63.30) 0.55 (0.50–0.61) 8.2E-33 5.12E-31 0 0.976
      Imidazole lactate 138 38.45 (29.90–51.98) 0.98 (0.97–0.99) 0.001 0.008 0 >0.999
      Glycocholenate sulfate 150 52.19 (30.96–67.40) 0.80 (0.75–0.84) 4.6E-16 1.46E-14 0.232 0.008
      Phenylalanylserine 57 33.80 (29.81–67.40) 0.92 (0.89–0.94) 1.97E-09 3.61E-08 0 0.779
      Glycolithocholate sulfate 65 53.70 (30.64–55.53) 0.74 (0.68–0.80) 3.63E-12 8.43E-11 0.524 6.28E-07
      Estriol 238 43.86 (30.15–62.86) 0.98 (0.97–0.99) 3.77E-06 4.57E-05 0 >0.999
      X-18914 47 42.91 (29.80–64.66) 0.92 (0.89–0.96) 3.05E-05 0.0003 0 0.974
      X-12822 167 56.20 (29.88–66.03) 0.96 (0.94–0.98) 2.47E-06 3.15E-05 0 0.983
      Prolyl-hydroxyproline 4 37.54 (29.91–38.77) 0.67 (0.53–0.84) 0.001 0.004 0 0.904

      OR, odds ratio; CI, confidence interval.

      Table 1 Genetically predicted metabolites positively associated with gastric cancer risk

      OR, odds ratio; CI, confidence interval; LDL, low-density lipoprotein; VLDL, very low-density lipoprotein.

      Table 2 Genetically metabolites negatively associated with gastric cancer risk

      OR, odds ratio; CI, confidence interval.

      TOP