Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Microbiol., 23 January 2026

Sec. Microorganisms in Vertebrate Digestive Systems

Volume 16 - 2025 | https://doi.org/10.3389/fmicb.2025.1730313

This article is part of the Research TopicMicrobiome and its Roles in Disease Diagnosis and Treatment: Pathogen Resistance Spectrum, Metabolism, Risk Model, and Vaccine DesignView all 14 articles

Gut microbiota, circulating metabolites, and pancreatic cancer risk: a multi-method causal inference study with cross-population validation

Shicheng LinShicheng Lin1Enze ShiEnze Shi1Yuxin ZhangYuxin Zhang1Xiaofan WangXiaofan Wang1Zhen TianZhen Tian2Jing Han
Jing Han3*Quanwang Li
Quanwang Li2*
  • 1Graduate School, Beijing University of Chinese Medicine, Beijing, China
  • 2Department of Oncology, Dongfang Hospital, Beijing University of Chinese Medicine, Beijing, China
  • 3Department of Laboratory Medicine, Dongfang Hospital, Beijing University of Chinese Medicine, Beijing, China

Pancreatic cancer (PC) is a lethal malignancy with limited early detection strategies and poor therapeutic response. Emerging evidence implicates the gut microbiota in carcinogenesis, yet whether microbial alterations are causal or secondary remains uncertain. In this study, we integrated cross-sectional 16S rDNA sequencing, two-sample Mendelian randomization (MR), and mediation analysis to investigate the causal role of gut microbiota in PC risk. We profiled fecal microbiota in a Beijing-based cohort of 26 newly diagnosed PC patients and 9 healthy controls, revealing significant dysbiosis characterized by reduced microbial diversity, depletion of butyrate-producing genera (e.g., Faecalibacterium), and enrichment of pro-inflammatory taxa such as Olsenella. Using European GWAS summary data, MR analysis identified 17 gut microbial taxa causally associated with PC risk, including Olsenella and Pauljensenia sp000411415. Notably, higher abundance of Pauljensenia sp000411415 was associated with increased PC risk, an effect partially mediated by reduced circulating levels of octanoylcarnitine (C8) and glutarylcarnitine (C5-DC)—metabolites independently linked to lower PC risk. Population-matched MR in East Asian cohorts validated several causal associations, enhancing ancestral relevance. Our findings support a causal role for specific gut microbes in pancreatic carcinogenesis and highlight a Pauljensenia–acylcarnitine axis whereby microbial suppression of protective metabolites may contribute to disease development. This integrative approach bridges microbial dysbiosis with functional mechanisms, offering novel insights for microbiome-informed strategies in PC prevention and early detection.

Introduction

Pancreatic cancer (PC) is one of the most lethal malignancies worldwide, with a mortality rate nearly equal to its incidence. It ranks as the sixth leading cause of cancer-related death globally (Bray et al., 2024). Due to its insidious onset, the majority of patients are diagnosed at an advanced, surgically unresectable stage. Furthermore, PC exhibits limited responsiveness to conventional therapies—including chemotherapy, radiotherapy, targeted therapy, and immunotherapy—resulting in a dismal prognosis and a 5-year survival rate of only 10% (Wu et al., 2022; Wang et al., 2023). Early detection offers the best opportunity for curative intervention and improved survival, highlighting the urgent need to identify modifiable risk factors that could inform prevention, early diagnosis, and therapeutic strategies.

Accumulating evidence has revealed the presence of a low-biomass intratumoral microbiome in pancreatic tissue—previously considered sterile—implicating microbial communities in tumor progression, diagnosis, and treatment response (Luo et al., 2025). The gut, the largest reservoir of commensal microorganisms in the human body, is believed to be a primary source of these intratumoral microbes, potentially seeding the pancreas via the portal circulation or enterohepatic pathways (Guan et al., 2023; Thomas et al., 2018). Although tumor colonization may originate in the upper gastrointestinal tract, microbial signals from proximal sites can propagate downstream and are detectable in stool (Vilà-Quintana et al., 2024; Tavanaeian et al., 2025; Chung et al., 2021). For instance, Gemella morbillorum and Fusobacterium nucleatum subsp. vincentii—taxa linked to pancreatic tumors (Sidiropoulos et al., 2024)—have also been reported in fecal samples from patients with pancreatic cancer (Chung et al., 2021). While fecal profiles may not fully recapitulate the spatial origin of tumor-associated microbes, their non-invasive accessibility supports their use as a practical surrogate for risk stratification and early detection.

The gut microbiota interacts dynamically with the host, modulating immune responses, regulating inflammation, and influencing carcinogenesis through the production of bioactive metabolites and modulation of signaling pathways (Tilg et al., 2020; Koh et al., 2016). However, microbial composition is highly susceptible to environmental influences such as diet, antibiotic use, chronic inflammation, and nutritional status, raising uncertainty about whether observed microbial alterations in PC are causal drivers or secondary consequences of disease. Moreover, bidirectional crosstalk between the gut microbiome and host metabolism may contribute to tumor metabolic reprogramming. Gut microbes participate in diverse metabolic processes, influence cancer-related signaling cascades, and shape host immune function, thereby exerting either pro- or anti-tumorigenic effects (Meng et al., 2023). Yet, it remains unclear whether circulating metabolites might mediate the influence of gut microbiota on PC susceptibility.

Despite growing interest in the gut–pancreas axis, critical questions regarding causality and mechanism remain unresolved. Large-scale human studies are often limited by confounding and reverse causation, particularly in observational designs. Recently, Mendelian randomization (MR) has emerged as a powerful tool for causal inference (Smith, 2010; Emdin et al., 2017), leveraging germline genetic variants as instrumental variables to mimic randomized controlled trials (Burgess et al., 2023). By minimizing confounding and avoiding reverse causality, MR enables robust assessment of whether gut microbiota alterations are causally linked to PC risk—not merely associated.

However, current research on the causal relationship between gut microbiota and pancreatic cancer remains limited, largely confined to single-cohort studies or observational analyses. Critical gaps include the lack of multi-population validation and systematic exploration of underlying biological mechanisms. For instance, no study to date has systematically integrated microbiome profiling with serum metabolomics to dissect causal pathways along the “gut–pancreas axis” (Chen Z. et al., 2024; Li et al., 2023; Li and Liang, 2024). Moreover, existing Mendelian randomization (MR) evidence is predominantly derived from European populations, with a striking paucity of data from Chinese or East Asian cohorts (Jiang et al., 2023; Zhu and Ding, 2024). This geographical bias severely limits the generalizability of findings to non-European ancestries.

To address these limitations, we integrated 16S rDNA sequencing data from a Beijing-based cohort with international genomic resources, employing two-sample Mendelian randomization (MR) and mediation analysis to assess causal links between gut microbiota and PC. We identified causally associated taxa in European-ancestry populations—such as the pro-inflammatory genus Olsenella and the risk-associated Pauljensenia sp000411415—with effects partially mediated by serum acylcarnitines, whose depletion may contribute to pancreatic carcinogenesis. In East Asian MR analyses, while the specific set of significantly associated taxa differed from that in Europeans, several microbial signals showed directionally consistent effects, suggesting potential shared biological pathways alongside population-specific associations. These cross-ancestry comparisons help contextualize our findings within diverse genetic and environmental backgrounds and highlight the need for inclusive, multi-ethnic studies in microbiome–disease research. This work provides mechanistic insights into the “gut–pancreas axis” and underscores microbial–metabolite pathways as promising targets for early PC detection.

Methods

Study approval, registration, and informed consent

An overview of the study design is presented in Figure 1. This study was approved by the IRB of Dongfang Hospital Beijing University of Chinese Medicine (approval number: JDF-IRB-2024055002) and registered at the Chinese Clinical Trial Registry (ChiCTR, ID: ITMCTR2025001710). Written informed consent was obtained from all participants in accordance with the Declaration of Helsinki.

Figure 1
Diagram showing the relationship between gut microbiota and pancreatic cancer through two studies. The cross-sectional study includes

Figure 1. Study design flowchart. This study integrates a cross-sectional 16S rDNA sequencing analysis with two-sample MR to investigate the association between gut microbiota and pancreatic cancer, and to explore causal mechanisms and mediating roles of serum metabolites. PC, pancreatic cancer; HC, healthy controls.

Study design and sample collection

We recruited patients with newly diagnosed PC from Dongfang Hospital, Beijing University of Chinese Medicine and Cancer Hospital of the Chinese Academy of Medical Sciences, between October 2024 and May 2025. Healthy controls (HC) were enrolled from community residents and health check-up centers in southern Beijing, matched to PC patients by age, sex, and body mass index (BMI).

Inclusion criteria for PC patients were: (1) histologically confirmed diagnosis of pancreatic cancer; (2) no prior chemotherapy, radiotherapy, targeted therapy, immunotherapy, or other anti-tumor treatments; (3) age 18–75 years; and (4) provision of written informed consent. Controls were free of chronic diseases and gastrointestinal disorders. Exclusion criteria for both groups included: (1) history of malignancy treated with anti-cancer therapies; (2) prior gastrointestinal conditions (e.g., inflammatory bowel disease, peptic ulcer, gastroenteritis, irritable bowel syndrome, or gastrointestinal tumors); (3) use of antibiotics or corticosteroids within the past 2 months; or (4) intake of probiotics or prebiotics within the past week. Demographic data—including age, sex, and BMI—were collected during recruitment. Pancreatic cancer diagnoses were independently verified by two pathologists at Beijing University of Chinese Medicine.

DNA extraction and 16S rDNA gene sequencing

Total genomic DNA was extracted from fecal samples using the Magnetic Soil and Stool DNA Kit (BioTeke, DP712) following the manufacturer’s protocol. DNA concentration was measured using Qubit (Invitrogen, United States). The V3–V4 hypervariable region of the 16S rDNA gene was amplified via PCR using primers 341F (5′-CCTACGGGNGGCWGCAG-3′) and 805R (5′-GACTACHVGGGTATCTAATCC-3′), under the following cycling conditions: initial denaturation at 98 °C for 30 s; 32 cycles of 98 °C for 10 s, 54 °C for 30 s, and 72 °C for 45 s; followed by final extension at 72 °C for 10 min. PCR products were purified using AMPure XT Beads (Beckman Coulter Genomics), quantified by Qubit, and assessed on an Agilent 2100 Bioanalyzer. Libraries were pooled and sequenced on an Illumina NovaSeq 6000 platform (PE250 mode) at LC-Bio Technology (Hangzhou, China).

16S rDNA gene sequence analysis

Raw sequencing reads were demultiplexed and processed using cutadapt (v1.9) and FLASH (v1.2.8) for primer trimming and paired-end read merging. Low-quality sequences (Phred score <20, length <100 bp, or containing >5% N bases) were filtered using fqtrim (v0.94). Chimeric sequences were removed with VSEARCH (v2.3.4). Amplicon sequence variants (ASVs) were inferred using DADA2 for high-resolution taxonomic profiling. Taxonomic annotation was performed using the QIIME2 plugin feature-classifier against the SILVA and NT-16S databases (16S rDNA gene, 99% identity). Alpha and beta diversity metrics were calculated in QIIME2. Differential abundance analysis at the genus level was performed using ANCOM-BC2 to account for the compositional nature of microbiome data and control for false discovery rate (FDR). Taxa with an FDR-corrected q-value <0.1 were considered statistically significant. All visualizations were generated in R (v3.4.4).

Gut microbiota and pancreatic cancer GWAS summary data

Summary-level genetic data for gut microbiota were obtained from a genome-wide association study (GWAS) of the Finnish FINRISK 2002 cohort (Qin et al., 2022), a population-based study comprising 5,959 individuals. Genome-wide association testing of 7,979,834 human genetic variants was conducted to identify host genetic determinants of gut microbial composition, as profiled by metagenomic sequencing. This analysis revealed 567 independent SNP–taxon associations, implicating host genetics in the regulation of gut microbiome structure. The study covered six geographic regions in Finland and reported results for 2,801 microbial taxa, classified within 473 GTDB taxa (including 59 phyla, 95 classes, 187 orders, 415 families, 922 genera, and 1,123 species). Comprehensive summary statistics are publicly available in the NHGRI-EBI GWAS Catalog (accession numbers GCST90032172–GCST90032644). The original study was approved by the Coordinating Ethics Committee of the Helsinki and Uusimaa Hospital District (reference number 558/E3/2001). Pancreatic cancer data were obtained from FinnGen release R12 (pancreatic malignant neoplasms, MNP; cases = 3,139, controls = 378,749).

Gut microbiome and metabolic pathway data for Asian populations were obtained from the China Nucleotide Sequence Archive (CNSA) (Liu et al., 2022). The pancreatic cancer outcome dataset for East Asian populations was sourced from Biobank Japan, with identifier GCST90018673 (Sakaue et al., 2021).

Mendelian randomization analysis

We employed the TwoSampleMR R package1 to assess causal relationships between gut microbiota and pancreatic cancer risk. Forward MR analyses were conducted with microbial taxa as exposures and PC as the outcome; reverse MR analyses tested the inverse direction. Instrumental variables (IVs) were defined as genome-wide significant (p < 1 × 10−5) SNPs with minor allele frequency (MAF) >0.01. To minimize weak instrument bias, only SNPs with F-statistics >10 were retained. Clumping was performed to ensure independence (r2 < 0.001, window size = 10 Mb).

The inverse-variance weighted (IVW) method was used as the primary approach to estimate causal effects. Continuous outcomes are reported as beta coefficients (β) with standard errors; binary outcomes are presented as odds ratios (ORs) with 95% confidence intervals (CIs). Significance was set at p < 0.05. Heterogeneity was evaluated using Cochran’s Q statistic (p < 0.05 in IVW or MR-Egger indicating heterogeneity). Horizontal pleiotropy was assessed via the MR-Egger intercept test; a significant deviation from zero (p < 0.05) indicated presence of pleiotropy. Effect estimates and 95% CIs were derived using fixed-effect IVW models.

Mediation analysis of gut microbiota, serum metabolites, and pancreatic cancer

To investigate whether circulating metabolites mediate the gut microbiota–PC association, we performed a two-step Mendelian randomization mediation analysis. Genetic instruments for serum metabolites were obtained from a large-scale GWAS of blood metabolomics in European populations (Sollis et al., 2023) (GWAS Catalog: GCST90199621–GCST90201020), encompassing 1,091 individual metabolites and 309 metabolite ratios after stringent quality control (Chen et al., 2023).

The total effect (βEO) of gut microbiota on PC risk was decomposed into: (1) the indirect effect mediated through a given serum metabolite, estimated as βEM × βMO, where βEM is the causal effect of the microbial exposure on the metabolite, and βMO is the causal effect of the metabolite on PC; and (2) the direct effect. The proportion mediated was calculated as (βEM × βMO)/βEO. All mediation tests were corrected for multiple testing using the false discovery rate (FDR) method.

Results

Demographic and clinical characteristics

A total of 26 patients with PC and 9 HC were included in this study. All participants were of Han Chinese ethnicity and resided in northern China, where the typical diet includes wheat, rice, meat, vegetables, and legumes. Although dietary patterns varied among individuals, none of the participants were vegetarians. Demographic and clinical characteristics are summarized in Table 1. No significant differences were observed between the two groups in terms of sex, age, BMI, smoking status, alcohol consumption, or history of diabetes (p > 0.05).

Table 1
www.frontiersin.org

Table 1. Demographic and clinical data.

Gut microbiota dysbiosis in pancreatic cancer patients

To investigate gut microbiota differences between PC patients and HC, we collected fecal samples from individuals at Dongfang Hospital, Beijing University of Chinese Medicine, and the Cancer Hospital of the Chinese Academy of Medical Sciences. 16S rDNA sequencing was performed to profile microbial composition. Alpha diversity analysis revealed that, compared with the HC group, the PC group exhibited significantly lower Shannon and Simpson indices (p = 0.028 and p = 0.018, respectively; Figures 2A,B), indicating reduced microbial diversity, while no significant differences in community richness were observed between the two groups (Supplementary Figures S1A–C). Principal coordinates analysis (PCoA) was used to assess overall microbial community structure, and analysis of similarities (ANOSIM) revealed a significant separation between the PC and HC groups (p = 0.001, Figure 2C), suggesting marked gut microbiota dysbiosis in PC patients.

Figure 2
Graphs and visualizations comparing microbiome data between two groups labeled HC and PC. A and B show box plots for Shannon and Simpson diversity indices. C is a Principal Coordinates Analysis plot indicating differences in microbial composition. D and E are bar charts showing taxonomy composition at the phylum and genus levels. F is a bar graph showing mean proportions with 95% confidence intervals for various genera. G is a cladogram visualizing phylogenetic relationships and abundance.

Figure 2. Bioinformatics analysis of 16S rDNA gene sequencing data from a cross-sectional cohort of pancreatic cancer patients (PC, n = 26) and healthy controls (HC, n = 9). (A) Shannon index and (B) Simpson index were significantly lower in the PC group (p < 0.05), indicating reduced alpha diversity at the ASV level. (C) PCoA based on Bray–Curtis dissimilarity revealed significant separation of microbial community structure between groups (PERMANOVA, R2 = 0.07, p = 0.001), indicating altered beta diversity in pancreatic cancer. (D) Relative abundance of gut microbiota at the phylum level. (E) Relative abundance of gut microbiota at the genus level. Taxonomic composition was visualized using stacked bar plots, showing distinct microbial profiles between PC and HC groups. (F) Differential abundance plot showing mean relative abundances and between-group differences of genera identified as differentially abundant by ANCOM-BC2 (q < 0.05) in PC patients versus HC. The PC group exhibited significantly lower mean abundance of g__Faecalibacterium, whereas g__Olsenella, g__Escherichia–Shigella, and g__Veillonella were markedly enriched. (G) Cladograms depicting the phylogenetic relationships among taxa identified as differentially abundant by ANCOM-BC2, with branches colored according to the group (PC or HC) in which each taxon is enriched.

Further taxonomic analysis revealed significant alterations in gut microbial composition between PC patients and HC. At the phylum level, the PC group exhibited decreased relative abundances of Firmicutes, Actinobacteriota, and Bacteroidota, accompanied by a notable increase in Proteobacteria (Figure 2D). At the genus level, PC patients showed a marked depletion of multiple butyrate-producing taxa, including Faecalibacterium, Roseburia, Lachnospira, Fusicatenibacter, Butyricicoccus, and Ruminococcus, while genera such as Veillonella, Megasphaera, Fusobacterium, and Olsenella were significantly enriched (Figure 2E), indicating a shift toward a pro-inflammatory microbial profile and intestinal dysbiosis. To rigorously identify differentially abundant taxa while accounting for the compositional nature of microbiome data, we performed ANCOM-BC2 analysis (q < 0.05) (Figures 2F,G). This analysis revealed distinct microbial signatures: in the PC group, Escherichia-Shigella and Olsenella were significantly enriched, taxa often associated with mucosal inflammation and pathogenic potential. In contrast, Faecalibacterium—a well-characterized anti-inflammatory commensal known for its role in maintaining gut barrier integrity and immune homeostasis—was significantly depleted in PC patients. These findings collectively indicate a functionally distinct gut microbiota composition in PC patients compared to healthy individuals.

Causal relationships between gut microbiota and pancreatic cancer

To further investigate the causal relationships between gut microbiota and pancreatic cancer, we first performed forward MR analysis, which identified 17 bacterial taxa significantly associated with pancreatic cancer risk. A circular heatmap was generated to visually summarize all findings (Figure 3), with the IVW method used as the primary analytical approach (Figure 4). Bacterial taxa are reported using GTDB release 89 nomenclature. Specifically, 11 taxa were associated with an increased risk of pancreatic cancer: CAG-83 sp000435555 (OR = 1.16, 95% CI: 1.02–1.33; p = 0.029), Clostridium M sp001304855 (OR = 1.58, 95% CI: 1.05–2.38; p = 0.029), GCA-900066495 (OR = 1.30, 95% CI: 1.00–1.69; p = 0.047), Olsenella C (OR = 1.51, 95% CI: 1.06–2.14; p = 0.021), Pauljensenia sp000411415 (OR = 1.30, 95% CI: 1.03–1.64; p = 0.026), TMED109 (OR = 1.63, 95% CI: 1.09–2.43; p = 0.017), Treponema D (OR = 1.25, 95% CI: 1.01–1.55; p = 0.042), UBA1033 sp001695555 (OR = 2.29, 95% CI: 1.28–4.11; p = 0.006), UBA1777 sp002320035 (OR = 1.63, 95% CI: 1.03–2.58; p = 0.036), and UBA8517 (OR = 1.94, 95% CI: 1.21–3.13; p = 0.006).

Figure 3
Circular visualization illustrating odds ratios (OR) and inverse-variance weight (IVW) p-values. Concentric rings compare different methods, including MR Egger, weighted median, IVW, simple mode, and weighted mode, with color gradients from green to red indicating statistical significance.

Figure 3. Circular heatmap summarizing the causal associations between gut microbial taxa and pancreatic cancer risk from Mendelian randomization analysis. Each radial segment represents a microbial taxon at genus or higher level, ordered by phylogenetic relatedness. The inner ring displays the odds ratio (OR) for causal effect on pancreatic cancer, with red indicating increased risk (OR >1) and green indicating protective effects (OR <1). The outer ring shows the statistical significance (p-value) of each association, with red segments denoting p < 0.05. Taxa with significant causal effects are highlighted in red along the outer edge.

Figure 4
Forest plot showing odds ratios (OR) with 95% confidence intervals (CI) for various traits analyzed using the inverse variance weighted method. Traits like MF0015:fructose_degradation and s_Oribacterium_sinus show ORs greater than one, indicating increased risk, while others like f_Veillonellaceae show ORs less than one, indicating protection. The reference line is at one. Red dashed line aids comparison.

Figure 4. Forest plot of causal effects of gut microbial taxa on pancreatic cancer risk estimated by IVW Mendelian randomization analysis. Each row represents a microbial taxon at genus or higher level. The OR and 95% CI are shown for each association. Horizontal lines indicate the 95% CI, with black squares representing the OR estimate. Dashed vertical line marks OR = 1 (no effect). Taxa with statistically significant associations (p < 0.05) are indicated by red dots, and FDR-adjusted p-values are reported. Positive ORs indicate increased risk; negative ORs would indicate protective effects.

Conversely, six gut microbial taxa showed protective associations against pancreatic cancer, of which five remained significant after FDR correction. These included Poseidoniaceae (OR = 0.41, 95% CI: 0.20–0.83; p = 0.034), Provencibacterium (OR = 0.76, 95% CI: 0.58–0.99; p = 0.047), Tepidanaerobacteraceae (OR = 0.50, 95% CI: 0.31–0.79; p = 0.002), UBA1066 (OR = 0.49, 95% CI: 0.28–0.85; p = 0.011), and UBA1409 (OR = 0.73, 95% CI: 0.57–0.94; p = 0.013). Sensitivity analyses indicated that Treponema D exhibited significant horizontal pleiotropy, as revealed by Cochran’s Q test and MR-Egger regression intercept analysis. Despite this, scatter plots and leave-one-out analysis further supported the robustness of the remaining associations (Supplementary Figures S2, S3). Additionally, MR-PRESSO analysis detected no outliers.

To assess reverse causation, we conducted reverse MR analyses with pancreatic cancer as the exposure. Significant reverse causal associations were observed for Treponema D and UBA1409 (p < 0.05), indicating potential bidirectional effects. To minimize bias from reverse causality, these two taxa were excluded from further interpretation in the forward MR analysis.

Causal relationships between serum metabolites and pancreatic cancer

To explore the potential mechanisms underlying the relationship between gut microbiota and pancreatic cancer, we conducted a two-step mediation MR analysis on serum metabolites. Initially, we performed a two-sample MR analysis to investigate the causal relationships between serum metabolites and pancreatic cancer. Subsequently, we assessed the causal impact of gut microbiota on significant serum metabolites.

MR analysis of 1,400 serum metabolites, with correction for multiple testing using the FDR, identified 10 metabolites as risk factors for pancreatic cancer (Figure 5): taurine to cysteine ratio, alpha-ketoglutarate to pyruvate ratio, X-11850, pentadecanoate (15:0), palmitoyl-sphingosine-phosphoethanolamine (d18:1/16:0), X-12714, 2-linoleoylglycerol (18:2), 1-methylhistidine, ceramide (d18:1/24:1), and 4-ethylcatechol sulfate. Conversely, 22 metabolites were identified as protective factors against pancreatic cancer: hexadecenedioate (C16:1-DC), methyl glucopyranoside (alpha + beta), 5-methyluridine (ribothymidine), glutarylcarnitine (C5-DC), 1-oleoyl-GPG (18:1), octanoylcarnitine (c8), hexanoylglutamine, adenosine 5′-monophosphate (AMP) to threonine ratio, ethyl beta-glucopyranoside, X-21733, X-12707, X-12411, glycosyl-N-tricosanoyl-sphingadienine (d18:2/23:0), 1-oleoyl-2-docosahexaenoyl-GPC (18:1/22:6), glycochenodeoxycholate, ornithine, S-methylcysteine sulfoxide, X-16580, 1-lignoceroyl-GPC (24:0), phenol sulfate, hydroxyasparagine, and citrulline to phosphate ratio.

Figure 5
Table with various traits and their respective statistical analysis results. Traits are listed alongside the method used, sample size (nsnp), p-value (pval), false discovery rate (FDR), and odds ratio with confidence interval (OR[95%CI]). A horizontal forest plot visually represents the OR values, showing traits' associations with risk or protection.

Figure 5. Forest plot of causal effects of serum metabolites on pancreatic cancer risk estimated by IVW Mendelian randomization analysis. Each row represents a metabolite or metabolic ratio. Odds ratios (OR) and 95% confidence intervals (CI) are shown, with red dots indicating statistically significant associations (p < 0.05). The dashed vertical line marks OR = 1 (no effect).

Sensitivity analyses using Cochran’s Q test and MR-Egger regression intercept showed no significant heterogeneity or horizontal pleiotropy. However, reverse MR analysis using the IVW method indicated a significant reverse causal association between phenol sulfate and pancreatic cancer (p < 0.05), suggesting potential bidirectional effects. These results highlight a complex metabolic profile associated with pancreatic cancer risk.

Gut microbiota-mediated effects on pancreatic cancer via serum metabolites

We further identified significant causal relationships between specific gut microbial taxa and serum metabolites. Notably, a decrease in Pauljensenia sp000411415 was associated with elevated levels of octanoylcarnitine (c8) and glutarylcarnitine (C5-DC) (β = −0.175 and −0.155, FDR = 0.032 and 0.034, respectively) (Supplementary Table S5). Integrating these findings with the earlier MR results for metabolites and pancreatic cancer, we observed that higher levels of both octanoylcarnitine (c8) and glutarylcarnitine (C5-DC) are causally linked to a reduced risk of pancreatic cancer. To formally test for mediation, we applied the product-of-coefficients method, which confirmed that octanoylcarnitine (c8) and glutarylcarnitine (C5-DC) partially mediate the effect of Pauljensenia sp000411415 on pancreatic cancer risk. Specifically, 6.5% via octanoylcarnitine (c8) (indirect effect = 0.017) and 5.5% via glutarylcarnitine (C5-DC) (indirect effect = 0.015) (Figure 6 and Table 2). Given that acylcarnitines like octanoylcarnitine (c8) and glutarylcarnitine (C5-DC) are critical for mitochondrial fatty acid oxidation and their depletion is linked to metabolic stress and inflammation—a known driver of pancreatic carcinogenesis—our findings provide quantitative evidence that Pauljensenia sp000411415 may promote pancreatic cancer risk, in part, by suppressing systemic acylcarnitine levels.

Figure 6
Two mediation analysis diagrams. (A) Pauljensenia sp000411415 as exposure affects pancreatic cancer with Octanoylcarnitine (c8) as mediator. Path coefficients: βEM=-0.17, βEO=0.263, βMO=-0.10. 6.5% effect mediation. (B) Similar setup with Glutarylcarnitine (c5-dc) as mediator. Path coefficients: βEM=-0.15, βEO=0.263, βMO=-0.09. 5.5% effect mediation.

Figure 6. Mediation analysis showing the role of serum metabolites as mediators in the causal pathway from gut microbiota to pancreatic cancer risk. (A) Octanoylcarnitine (C8) mediates 6.5% of the effect of Paulijennia sp000411415 on pancreatic cancer. (B) Glutarylcarnitine (C5-DC) mediates 5.5% of the effect. Solid arrows indicate significant associations (p < 0.05); β coefficients are shown for each path. βEM, effects of exposure on mediator; βMO, effects of mediator on outcome; βEO, effects of exposure on outcome.

Table 2
www.frontiersin.org

Table 2. The mediation effect of gut microbiota on pancreatic cancer via affecting serum metabolites.

Validation of causal associations in an East Asian context using population-matched data

Given the potential limitations of applying European-derived genetic instruments to infer causal effects in a Chinese discovery cohort, we sought to validate our findings using population-matched data. We performed a two-sample Mendelian randomization analysis to assess the causal effects of the differentially abundant gut microbial taxa—identified in our Chinese cross-sectional cohort—on pancreatic cancer risk in East Asian populations.

Using genetic instruments derived from East Asian GWAS, we identified 12 microbial taxa and 2 metabolic pathways significantly associated with pancreatic cancer risk (Figure 7). Among them, 8 taxa and one pathway were inversely associated with disease risk and thus considered potential protective factors. These included s_Clostridium_saccharolyticum, s_Acinetobacter_baumannii, g_Eremococcus, f_Veillonellaceae, f_Brachyspiraceae, g_Escherichia, s_Escherichia_coli, and the threonine degradation I pathway. In contrast, 4 microbial taxa and one metabolic pathway were positively associated with pancreatic cancer risk, indicating potential risk factors. These were s_Bacteroides_caccae, g_Kingella, s_Pseudomonas_stutzeri, g_Acidovorax, s_Oribacterium_sinus, and the fructose degradation pathway. Sensitivity analysis using MR-Egger intercept and MR-PRESSO detected evidence of horizontal pleiotropy for s_Oribacterium_sinus (p < 0.05), suggesting that the observed association may be influenced by pleiotropic effects.

Figure 7
Forest plot displaying odds ratios and confidence intervals for various traits analyzed using the inverse variance weighted method. Traits include UBA1033, Sorangium, UBA8517, and others. The odds ratios range from protecting to risk, marked with a red dashed line at 1. Key values: UBA1033 has an odds ratio of 2.29 (1.28-4.11), Thermoprotei has an odds ratio of 0.32 (0.10-0.96). The plot highlights significant associations determined by p-values and false discovery rates.

Figure 7. Forest plot of causal effects of gut microbial taxa on pancreatic cancer risk in East Asian populations, estimated by inverse variance-weighted Mendelian randomization analysis. Each row represents a microbial taxon at genus or higher level. Odds ratios (OR) and 95% confidence intervals (CI) are shown for the association between each taxon and pancreatic cancer. Horizontal lines indicate 95% CIs; results are based on genetic instrumental variables derived from Chinese cohort data and East Asian GWAS of pancreatic cancer.

Discussion

The human gut harbors one of the most complex and densely populated microbial ecosystems, collectively known as the gut microbiota, which has co-evolved in a symbiotic relationship with its host (Athanasopoulou et al., 2023). Beyond its role in digestion and nutrient metabolism, the gut microbiome functions as a dynamic metabolic and immunomodulatory organ, producing a vast array of bioactive molecules—including neurotransmitters, organic acids, and immune mediators—that can influence host physiology both locally and at distant sites (Kim et al., 2024). Disruption of this finely tuned microbial balance, a state termed dysbiosis, has been increasingly implicated as a potential trigger or amplifier in the development of chronic diseases, including cancer (Asseri et al., 2023). In the context of pancreatic cancer, emerging evidence suggests that alterations in gut microbial composition and function may play a significant role in both tumor initiation and progression (Guo and Shao, 2025).

PC remains one of the most lethal malignancies worldwide, characterized by intricate pathogenesis and a paucity of early-detection biomarkers and curative therapies (Siegel et al., 2021). Here, we integrated microbiome profiling with rigorous causal-inference algorithms and mediation analysis to delineate the contribution of the gut microbiota to pancreatic carcinogenesis. We provide convergent evidence that microbial dysbiosis is not merely an epiphenomenon but functionally implicated in tumor initiation and progression through defined microbial taxa and metabolite-signaling axes, thereby unveiling tractable targets for mechanistic dissection and therapeutic intervention. Our findings not only confirm the presence of profound microbial dysbiosis in patients but also move beyond association to identify specific microbes and metabolite pathways with causal relevance, offering new insights into potential mechanisms and intervention strategies.

The gut microbial ecosystem in PC patients exhibits a consistent pattern of structural and functional dysbiosis, marked by depletion of butyrate-producing Firmicutes—such as Faecalibacterium, Roseburia, and Lachnospira—and expansion of pro-inflammatory genera including Escherichia–Shigella, Veillonella, and Fusobacterium. This pattern is corroborated by a recent European study linking reduced Lachnospira abundance to intraductal papillary mucinous neoplasm (IPMN), a key PC precursor (Sidiropoulos et al., 2024). Strikingly, many of these taxa are also enriched in other microbial niches. Oral microbiome studies report elevated Prevotella, Rothia, Veillonella, Porphyromonas, Actinomyces, and Fusobacterium in PC patients (Meng et al., 2025; Uguz et al., 2025), while tumor tissue analyses have detected increased abundances of Escherichia coli and Shigella—members of the Escherichia–Shigella complex identified in our cohort—particularly in pancreatic head cancers (Tajpara et al., 2025). Together, these multi-compartment observations suggest systemic selection of specific bacterial lineages across anatomically distinct sites under shared immunometabolic pressures of malignancy, underscoring a broader paradigm in gastrointestinal carcinogenesis: a systemic shift from a symbiotic to an inflammogenic microbiota (Mori et al., 2018; Pötgens et al., 2018). Butyrate and other short-chain fatty acids (SCFAs) are critical for maintaining intestinal barrier function, regulating T-regulatory cell differentiation, and suppressing NF-κB signaling (Jin et al., 2025; Zheng and Xu, 2025; Xu et al., 2023). Their reduction may facilitate systemic inflammation and immune dysregulation, creating a permissive environment for tumor initiation and progression. The pronounced reduction of Faecalibacterium—a dominant anti-inflammatory commensal typically enriched in healthy individuals—aligns with findings in pancreatic and colorectal cancers (Nagata et al., 2022; Li J. et al., 2024); notably, reduced SNP density in its predominant species F. prausnitzii has been observed in specific pancreatic cancer subgroups (Han et al., 2025), hinting at altered strain-level dynamics beyond simple abundance loss.

While 16S rDNA profiling captures microbial community shifts associated with PC, it cannot distinguish whether these changes are drivers or consequences of disease. To address this, we applied two-sample MR to investigate causal relationships between gut microbiota and PC risk, identifying in European populations 17 bacterial taxa with significant causal associations—11 conferring increased risk and 6 exerting protective effects. Notably, among the taxa showing differential abundance in our 16S analysis, Olsenella—which was enriched in PC patients—also demonstrated a significant causal effect on elevated PC risk in MR analysis. This concordance between observational and causal evidence strengthens the plausibility of Olsenella as a potential contributor to pancreatic carcinogenesis. Olsenella, a genus within the Actinobacteriota phylum, has been previously associated with inflammatory conditions such as periodontitis, osteoporosis, and colorectal cancer (Vogtmann et al., 2025; Sheng et al., 2019; Qin et al., 2021; Dewhirst et al., 2001). A key characteristic of this genus is its capacity to produce lactate as a major fermentation product and to degrade mucin (Kraatz et al., 2011; Barandouzi et al., 2020). The overproduction of lactate may acidify the local gut environment, potentially favoring the expansion of other pathobionts and disrupting the symbiotic microbial community structure. More critically, mucin degradation by Olsenella is likely to compromise the integrity of the gastrointestinal mucus barrier—a crucial first line of host defense (Li S. et al., 2024). This barrier impairment could facilitate the translocation of microbial metabolites or even live bacteria into systemic circulation, potentially reaching the pancreas via the portal vein or enterohepatic circulation (Wang et al., 2019). This process may, in turn, trigger chronic systemic inflammation and activate immune cells, thereby fostering a pro-tumorigenic microenvironment that could promote pancreatic carcinogenesis (Mager et al., 2020).

In contrast to abundance-altered, causal taxa like Olsenella, our MR findings unveil a subtler mechanism: microbial influence on PC can operate through functional shifts independent of population size, a paradigm exemplified by Pauljensenia sp000411415. Mediation analysis further revealed that increased abundance of Pauljensenia sp000411415 was positively associated with PC risk, and this effect was partially mediated by its association with reduced levels of circulating octanoylcarnitine (C8) and glutarylcarnitine (C5-DC). This identified Pauljensenia-acylcarnitine axis reveals a novel microbiome-host metabolic interaction mechanism influencing PC risk. While the precise biological functions of Pauljensenia sp000411415 remain incompletely characterized, its genetically predicted link to decreased octanoylcarnitine (C8) and glutarylcarnitine (C5-DC) levels suggests a potential role in disrupting mitochondrial fatty acid β-oxidation and energy metabolism. Medium-chain acylcarnitines, such as octanoylcarnitine (C8), participate in medium-chain fatty acid metabolism, which may originate from dietary sources or microbial fermentation (Solo et al., 2025). Glutarylcarnitine (C5-DC) is associated with lysine and tryptophan degradation pathways (Li Q. et al., 2021). Beyond serving as metabolic intermediates, these specific acylcarnitines may exert direct protective effects—for instance, by supplying alternative energy substrates to damaged tissues, alleviating oxidative stress through enhanced mitochondrial efficiency, or modulating immune cell function and differentiation (Brosnan and Brosnan, 2009; Waagsbø et al., 2016). We propose that Pauljensenia sp000411415 may promote PC development by persistently depleting these functional metabolites, thereby compromising host metabolic homeostasis and tissue repair capacity. These findings highlight a gut microbiota-acylcarnitine-host metabolic axis as a potential pathogenic mechanism in pancreatic carcinogenesis.

However, we also observed substantial population-specific associations. Several microbial taxa showed significant causal effects exclusively in European populations (e.g., Pauljensenia sp000411415, UBA1033 sp001695555) or East Asian populations (e.g., s_Escherichia_coli, s_Acinetobacter_baumannii). Notably, while the specific taxa differed between ancestries, several microbial signals exhibited directionally concordant effects—such as enrichment of pro-inflammatory or mucus-degrading genera and depletion of butyrate producers—suggesting convergent dysbiosis patterns along shared pathogenic pathways. This heterogeneity likely reflects differences in genetic background, dietary patterns, and environmental exposures between populations. For example, the higher intake of plant-based fibers and fermented foods in traditional East Asian diets may promote the colonization or metabolic activity of beneficial microbial groups—particularly butyrate-producing taxa such as Faecalibacterium and Roseburia—which are known to support intestinal barrier integrity and immune homeostasis (Chen H. et al., 2024; Li G. et al., 2021). In our Beijing cohort, these taxa were markedly depleted in pancreatic cancer patients, consistent with a loss of protective functions that may be shaped by long-term dietary habits. Concurrently, population-specific genetic variants in host genes involved in mucosal immunity or bile acid metabolism can shape distinct host–microbiota interaction landscapes, thereby giving rise to ancestry-dependent causal associations (Graves and Milovanova, 2019). Importantly, the absence of exact taxonomic overlap does not undermine causality; rather, it underscores the necessity of accounting for population-specific factors in microbiota-targeted strategies for cancer prevention and treatment. Precision medicine initiatives should incorporate inter-population differences in gut microbiota composition and host–microbe interactions to develop tailored interventions that maximize therapeutic efficacy while minimizing unintended risks.

It is also instructive to contrast our MR findings with the 16S results. For instance, while MR suggested protective effects of f_Veillonellaceae and g_Escherichia (s_Escherichia_coli), their unclassified counterparts—g__Veillonellaceae_unclassified and g__Escherichia-Shigella—were enriched in pancreatic cancer patients. Conversely, s_Bacteroides_caccae was identified as a risk factor in MR, yet other Bacteroides species showed reduced abundance in disease. These inconsistencies likely reflect the distinction between causal drivers and disease-associated ecological shifts. MR captures lifelong genetic predisposition to microbial abundance, reflecting upstream influences on carcinogenesis. In contrast, cross-sectional 16S data may capture downstream consequences of pancreatic dysfunction—such as biliary obstruction or immune suppression—that favor opportunistic colonization regardless of initial causality. Moreover, intra-genus functional heterogeneity means that not all members of a taxon share the same biological role. Thus, microbial enrichment in disease does not imply pathogenicity, reinforcing the need for causal inference to disentangle true contributors from bystanders.

Notably, our 16S rDNA gene-based approach is inherently limited to bacterial profiling and cannot assess changes in the virome or mycobiome. Nevertheless, fungi and viruses—through cross-kingdom interactions—are increasingly recognized as key modulators of bacterial community structure and host immunity, potentially influencing pancreatic carcinogenesis (Meng et al., 2025). Our team is currently conducting large-scale shotgun metagenomic sequencing studies to enable a more comprehensive, multi-kingdom characterization of the gut microbiome in pancreatic cancer.

Several limitations should be acknowledged. First, although we employed multiple sensitivity methods (e.g., MR-Egger, weighted median, MR-PRESSO) to mitigate horizontal pleiotropy, residual confounding through alternative pathways cannot be fully excluded. Second, our two-sample MR relies on summary-level GWAS data, which precludes adjustment for individual covariates, assessment of non-linear effects, and direct inference of temporality—highlighting the need for longitudinal validation (Visscher et al., 2017). Third, while we leveraged both European and East Asian datasets, population-specific genetic and environmental factors may limit generalizability to other ancestries. Fourth, our observational 16S cohort is modest in size and enriched for late-stage disease (only 6 early-stage cases), limiting power to detect stage-related microbial differences; findings require replication in larger, stage-stratified cohorts. Finally, although mediation analysis implicates a Pauljensenia–acylcarnitine axis in pancreatic cancer risk, its causal mechanism awaits functional validation in experimental models.

In conclusion, our integrative analysis reveals that gut microbiota influences pancreatic cancer risk through both compositional shifts and causal microbial effects, exemplified by the dual-evidence role of Olsenella, the functionally mediated protection of Pauljensenia, and the complex interplay between causal effects and disease-associated dysbiosis. These findings advance our understanding of the gut–pancreas axis and highlight the potential of combining observational and causal inference approaches to identify novel targets for early intervention and microbiome-based therapeutics.

Data availability statement

The data presented in this study are publicly available. The data can be found here: https://ngdc.cncb.ac.cn/gsa/browse/CRA031736. GSA Accession Number: CRA031736.

Ethics statement

The studies involving humans were approved by IRB of Dongfang Hospital Beijing University of Chinese Medicine. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

SL: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft. ES: Conceptualization, Data curation, Formal analysis, Software, Writing – original draft. YZ: Conceptualization, Data curation, Investigation, Writing – original draft. XW: Conceptualization, Data curation, Visualization, Writing – original draft. ZT: Conceptualization, Data curation, Investigation, Writing – original draft. JH: Investigation, Methodology, Software, Supervision, Validation, Writing – review & editing. QL: Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Validation, Visualization, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This project was supported by the National High Level Chinese Medicine Hospital Clinical Research Funding (DFRCZY-2024ARC-003).

Acknowledgments

The DNA libraries were sequenced on the illumina NovaseqTM 6000 platform by LC Bio Technology Co., Ltd. (Hangzhou, China). We are grateful to LC Bio Technology Co., Ltd. for assisting in sequencing and bioinformatics analysis.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that Generative AI was used in the creation of this manuscript. Generative AI was used to assist with language editing, refinement of scientific phrasing, and formatting of tables and figures to improve clarity and readability. All content was critically reviewed, validated, and approved by the authors. The use of AI did not extend to data analysis, interpretation of results, or writing of original scientific content, which were conducted entirely by the human authors.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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/fmicb.2025.1730313/full#supplementary-material

Footnotes

References

Asseri, A. H., Bakhsh, T., Abuzahrah, S. S., Ali, S., and Rather, I. A. (2023). The gut dysbiosis-cancer axis: illuminating novel insights and implications for clinical practice. Front. Pharmacol. 14:1208044. doi: 10.3389/fphar.2023.1208044,

PubMed Abstract | Crossref Full Text | Google Scholar

Athanasopoulou, K., Adamopoulos, P. G., and Scorilas, A. (2023). Unveiling the human gastrointestinal tract microbiome: the past, present, and future of metagenomics. Biomedicine 11:827. doi: 10.3390/biomedicines11030827,

PubMed Abstract | Crossref Full Text | Google Scholar

Barandouzi, Z. A., Starkweather, A. R., Henderson, W. A., Gyamfi, A., and Cong, X. S. (2020). Altered composition of gut microbiota in depression: a systematic review. Front. Psychiatry 11:541. doi: 10.3389/fpsyt.2020.00541,

PubMed Abstract | Crossref Full Text | Google Scholar

Bray, F., Laversanne, M., Sung, H., Ferlay, J., Siegel, R. L., Soerjomataram, I., et al. (2024). Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 74, 229–263. doi: 10.3322/caac.21834,

PubMed Abstract | Crossref Full Text | Google Scholar

Brosnan, M. E., and Brosnan, J. T. (2009). Hepatic glutamate metabolism: a tale of 2 hepatocytes. Am. J. Clin. Nutr. 90, 857s–861s. doi: 10.3945/ajcn.2009.27462Z,

PubMed Abstract | Crossref Full Text | Google Scholar

Burgess, S., Davey Smith, G., Davies, N., Dudbridge, F., Gill, D., Glymour, M. M., et al. (2023). Guidelines for performing Mendelian randomization investigations: update for summer 2023. Wellcome Open Res. 4:186. doi: 10.12688/wellcomeopenres.15555.3

Crossref Full Text | Google Scholar

Chen, Y., Lu, T., Pettersson-Kymmer, U., Stewart, I. D., Butler-Laporte, G., Nakanishi, T., et al. (2023). Genomic atlas of the plasma metabolome prioritizes metabolites implicated in human diseases. Nat. Genet. 55, 44–53. doi: 10.1038/s41588-022-01270-1,

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, H., Qian, Y., Jiang, C., Tang, L., Yu, J., Zhang, L., et al. (2024). Butyrate ameliorated ferroptosis in ulcerative colitis through modulating Nrf2/GPX4 signal pathway and improving intestinal barrier. Biochim. Biophys. Acta Mol. basis Dis. 1870:166984. doi: 10.1016/j.bbadis.2023.166984,

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, Z., Wang, Z., Bao, H., and Ma, S. (2024). Gut microbiota and pancreatic cancer risk, and the mediating role of immune cells and inflammatory cytokines: a Mendelian randomization study. Front. Immunol. 15:1408770. doi: 10.3389/fimmu.2024.1408770,

PubMed Abstract | Crossref Full Text | Google Scholar

Chung, M., Zhao, N., Meier, R., Koestler, D. C., Wu, G., del Castillo, E., et al. (2021). Comparisons of oral, intestinal, and pancreatic bacterial microbiomes in patients with pancreatic cancer and other gastrointestinal diseases. J. Oral Microbiol. 13:1887680. doi: 10.1080/20002297.2021.1887680,

PubMed Abstract | Crossref Full Text | Google Scholar

Dewhirst, F. E., Paster, B. J., Tzellas, N., Coleman, B., Downes, J., Spratt, D. A., et al. (2001). Characterization of novel human oral isolates and cloned 16S rDNA sequences that fall in the family Coriobacteriaceae: description of olsenella gen. nov., reclassification of Lactobacillus uli as Olsenella uli comb. nov. and description of Olsenella profusa sp. nov. Int. J. Syst. Evol. Microbiol. 51, 1797–1804. doi: 10.1099/00207713-51-5-1797,

PubMed Abstract | Crossref Full Text | Google Scholar

Emdin, C. A., Khera, A. V., and Kathiresan, S. (2017). Mendelian randomization. JAMA 318, 1925–1926. doi: 10.1001/jama.2017.17219,

PubMed Abstract | Crossref Full Text | Google Scholar

Graves, D. T., and Milovanova, T. N. (2019). Mucosal immunity and the FOXO1 transcription factors. Front. Immunol. 10:2530. doi: 10.3389/fimmu.2019.02530,

PubMed Abstract | Crossref Full Text | Google Scholar

Guan, S. W., Lin, Q., and Yu, H. B. (2023). Intratumour microbiome of pancreatic cancer. World J. Gastrointest. Oncol. 15, 713–730. doi: 10.4251/wjgo.v15.i5.713,

PubMed Abstract | Crossref Full Text | Google Scholar

Guo, X., and Shao, Y. (2025). Role of the oral-gut microbiota axis in pancreatic cancer: a new perspective on tumor pathophysiology, diagnosis, and treatment. Mol. Med. 31:103. doi: 10.1186/s10020-025-01166-w,

PubMed Abstract | Crossref Full Text | Google Scholar

Han, Y., Cao, B., Tang, J., and Wang, J. (2025). A comprehensive multi-omics analysis uncovers the associations between gut microbiota and pancreatic cancer. Front. Microbiol. 16:1592549. doi: 10.3389/fmicb.2025.1592549,

PubMed Abstract | Crossref Full Text | Google Scholar

Jiang, Z., Mou, Y., Wang, H., Li, L., Jin, T., Wang, H., et al. (2023). Causal effect between gut microbiota and pancreatic cancer: a two-sample Mendelian randomization study. BMC Cancer 23:1091. doi: 10.1186/s12885-023-11493-y,

PubMed Abstract | Crossref Full Text | Google Scholar

Jin, Y., Lu, L., Hua, H., Chen, B., Fang, W., Lin, K., et al. (2025). Qingxie Fuzheng granules attenuate cancer cachexia by restoring gut microbiota homeostasis and suppressing IL-6/NF-κB signaling in colorectal adenocarcinoma. Hereditas 162:178. doi: 10.1186/s41065-025-00541-1,

PubMed Abstract | Crossref Full Text | Google Scholar

Kim, S., Seo, S. U., and Kweon, M. N. (2024). Gut microbiota-derived metabolites tune host homeostasis fate. Semin. Immunopathol. 46:2. doi: 10.1007/s00281-024-01012-x,

PubMed Abstract | Crossref Full Text | Google Scholar

Koh, A., De Vadder, F., Kovatcheva-Datchary, P., and Bäckhed, F. (2016). From dietary Fiber to host physiology: short-chain fatty acids as key bacterial metabolites. Cell 165, 1332–1345. doi: 10.1016/j.cell.2016.05.041,

PubMed Abstract | Crossref Full Text | Google Scholar

Kraatz, M., Wallace, R. J., and Svensson, L. (2011). Olsenella umbonata sp. nov., a microaerotolerant anaerobic lactic acid bacterium from the sheep rumen and pig jejunum, and emended descriptions of Olsenella, Olsenella uli and Olsenella profusa. Int. J. Syst. Evol. Microbiol. 61, 795–803. doi: 10.1099/ijs.0.022954-0,

PubMed Abstract | Crossref Full Text | Google Scholar

Li, S., Chen, M., Wang, Z., Abudourexiti, W., Zhang, L., Ding, C., et al. (2024). Ant may well destroy a whole dam: glycans of colonic mucus barrier disintegrated by gut bacteria. Microbiol. Res. 281:127599. doi: 10.1016/j.micres.2023.127599,

PubMed Abstract | Crossref Full Text | Google Scholar

Li, J., Chen, Z., Wang, Q., du, L., Yang, Y., Guo, F., et al. (2024). Microbial and metabolic profiles unveil mutualistic microbe-microbe interaction in obesity-related colorectal cancer. Cell Rep. Med. 5:101429. doi: 10.1016/j.xcrm.2024.101429,

PubMed Abstract | Crossref Full Text | Google Scholar

Li, X. J., Gao, M. G., Chen, X. X., Rong, Y. M., Huang, L. L., and Huang, J. S. (2023). Genetically predicted causal effects of gut microbiota and gut metabolites on digestive tract cancer: a two-sample Mendelian randomization analysis. World. J. Oncol. 14, 558–569. doi: 10.14740/wjon1737,

PubMed Abstract | Crossref Full Text | Google Scholar

Li, X., and Liang, Z. (2024). Causal effect of gut microbiota on pancreatic cancer: a Mendelian randomization and colocalization study. J. Cell. Mol. Med. 28:e18255. doi: 10.1111/jcmm.18255,

PubMed Abstract | Crossref Full Text | Google Scholar

Li, G., Lin, J., Zhang, C., Gao, H., Lu, H., Gao, X., et al. (2021). Microbiota metabolite butyrate constrains neutrophil functions and ameliorates mucosal inflammation in inflammatory bowel disease. Gut Microbes 13:1968257. doi: 10.1080/19490976.2021.1968257,

PubMed Abstract | Crossref Full Text | Google Scholar

Li, Q., Yang, C., Feng, L., Zhao, Y., Su, Y., Liu, H., et al. (2021). Glutaric acidemia, pathogenesis and nutritional therapy. Front. Nutr. 8:704984. doi: 10.3389/fnut.2021.704984,

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, X., Tong, X., Zou, Y., Lin, X., Zhao, H., Tian, L., et al. (2022). Mendelian randomization analyses support causal relationships between blood metabolites and the gut microbiome. Nat. Genet. 54, 52–61. doi: 10.1038/s41588-021-00968-y,

PubMed Abstract | Crossref Full Text | Google Scholar

Luo, D., Chen, Q., Li, Y., Yang, J., Tao, Y., Ji, L., et al. (2025). Microbiome-metabolome interplay in pancreatic cancer progression: insights from multi-omics analysis. Mol. Cancer 24:240. doi: 10.1186/s12943-025-02458-9,

PubMed Abstract | Crossref Full Text | Google Scholar

Mager, L. F., Burkhard, R., Pett, N., Cooke, N. C. A., Brown, K., Ramay, H., et al. (2020). Microbiome-derived inosine modulates response to checkpoint inhibitor immunotherapy. Science 369, 1481–1489. doi: 10.1126/science.abc3421,

PubMed Abstract | Crossref Full Text | Google Scholar

Meng, Y. F., Fan, Z. Y., Zhou, B., and Zhan, H. X. (2023). Role of the intratumoral microbiome in tumor progression and therapeutics implications. Biochim. Biophys. Acta Rev. Cancer 1878:189014. doi: 10.1016/j.bbcan.2023.189014,

PubMed Abstract | Crossref Full Text | Google Scholar

Meng, Y., Wu, F., Kwak, S., Wang, C., Usyk, M., Freedman, N. D., et al. (2025). Oral bacterial and fungal microbiome and subsequent risk for pancreatic cancer. JAMA Oncol. 11, 1331–1340. doi: 10.1001/jamaoncol.2025.3377,

PubMed Abstract | Crossref Full Text | Google Scholar

Mori, G., Rampelli, S., Orena, B. S., Rengucci, C., de Maio, G., Barbieri, G., et al. (2018). Shifts of faecal microbiota during sporadic colorectal carcinogenesis. Sci. Rep. 8:10329. doi: 10.1038/s41598-018-28671-9,

PubMed Abstract | Crossref Full Text | Google Scholar

Nagata, N., Nishijima, S., Kojima, Y., Hisada, Y., Imbe, K., Miyoshi-Akiyama, T., et al. (2022). Metagenomic identification of microbial signatures predicting pancreatic cancer from a multinational study. Gastroenterology 163, 222–238. doi: 10.1053/j.gastro.2022.03.054,

PubMed Abstract | Crossref Full Text | Google Scholar

Pötgens, S. A., Brossel, H., Sboarina, M., Catry, E., Cani, P. D., Neyrinck, A. M., et al. (2018). Klebsiella oxytoca expands in cancer cachexia and acts as a gut pathobiont contributing to intestinal dysfunction. Sci. Rep. 8:12321. doi: 10.1038/s41598-018-30569-5,

PubMed Abstract | Crossref Full Text | Google Scholar

Qin, Y., Havulinna, A. S., Liu, Y., Jousilahti, P., Ritchie, S. C., Tokolyi, A., et al. (2022). Combined effects of host genetics and diet on human gut microbiota and incident disease in a single population cohort. Nat. Genet. 54, 134–142. doi: 10.1038/s41588-021-00991-z,

PubMed Abstract | Crossref Full Text | Google Scholar

Qin, Q., Yan, S., Yang, Y., Chen, J., Yan, H., Li, T., et al. (2021). The relationship between osteoporosis and intestinal microbes in the Henan Province of China. Front. Cell Dev. Biol. 9:752990. doi: 10.3389/fcell.2021.752990,

PubMed Abstract | Crossref Full Text | Google Scholar

Sakaue, S., Kanai, M., Tanigawa, Y., Karjalainen, J., Kurki, M., Koshiba, S., et al. (2021). A cross-population atlas of genetic associations for 220 human phenotypes. Nat. Genet. 53, 1415–1424. doi: 10.1038/s41588-021-00931-x,

PubMed Abstract | Crossref Full Text | Google Scholar

Sheng, Q., Du, H., Cheng, X., Tang, Y., Pan, L., Wang, Q., et al. (2019). Characteristics of fecal gut microbiota in patients with colorectal cancer at different stages and different sites. Oncol. Lett. 18, 4834–4844. doi: 10.3892/ol.2019.10841

Crossref Full Text | Google Scholar

Sidiropoulos, T., Dovrolis, N., Katifelis, H., Michalopoulos, N. V., Kokoropoulos, P., Arkadopoulos, N., et al. (2024). Dysbiosis signature of fecal microbiota in patients with pancreatic adenocarcinoma and pancreatic intraductal papillary mucinous neoplasms. Biomedicine 12:1040. doi: 10.3390/biomedicines12051040,

PubMed Abstract | Crossref Full Text | Google Scholar

Siegel, R. L., Miller, K. D., Fuchs, H. E., and Jemal, A. (2021). Cancer statistics, 2021. CA Cancer J. Clin. 71, 7–33. doi: 10.3322/caac.21654,

PubMed Abstract | Crossref Full Text | Google Scholar

Smith, G. D. (2010). Mendelian randomization for strengthening causal inference in observational studies: application to gene × environment interactions. Perspect. Psychol. Sci. 5, 527–545. doi: 10.1177/1745691610383505,

PubMed Abstract | Crossref Full Text | Google Scholar

Sollis, E., Mosaku, A., Abid, A., Buniello, A., Cerezo, M., Gil, L., et al. (2023). The NHGRI-EBI GWAS catalog: knowledgebase and deposition resource. Nucleic Acids Res. 51, D977–D985. doi: 10.1093/nar/gkac1010,

PubMed Abstract | Crossref Full Text | Google Scholar

Solo, K. J., Zhang, Y., Bharathi, S. S., Zhang, B. B., Richert, A. C., Schmidt, A. V., et al. (2025). Oral octanoylcarnitine alleviates exercise intolerance in mouse models of long-chain fatty acid oxidation disorders. JCI Insight 10:e199443. doi: 10.1172/jci.insight.199443,

PubMed Abstract | Crossref Full Text | Google Scholar

Tajpara, P., Sobkowiak, M. J., Healy, K., Naud, S., Gündel, B., Halimi, A., et al. (2025). Patient-derived pancreatic tumor bacteria exhibit oncogenic properties and are recognized by MAIT cells in tumor spheroids. Front. Immunol. 16:1553034. doi: 10.3389/fimmu.2025.1553034,

PubMed Abstract | Crossref Full Text | Google Scholar

Tavanaeian, S., Feizabadi, M. M., Falsafi, S., Aghdaei, H. A., and Houri, H. (2025). Oral and fecal microbiome alterations in pancreatic cancer: insights into potential diagnostic biomarkers. BMC Microbiol. 25:624. doi: 10.1186/s12866-025-04344-2,

PubMed Abstract | Crossref Full Text | Google Scholar

Thomas, R. M., Gharaibeh, R. Z., Gauthier, J., Beveridge, M., Pope, J. L., Guijarro, M. V., et al. (2018). Intestinal microbiota enhances pancreatic carcinogenesis in preclinical models. Carcinogenesis 39, 1068–1078. doi: 10.1093/carcin/bgy073,

PubMed Abstract | Crossref Full Text | Google Scholar

Tilg, H., Zmora, N., Adolph, T. E., and Elinav, E. (2020). The intestinal microbiota fuelling metabolic inflammation. Nat. Rev. Immunol. 20, 40–54. doi: 10.1038/s41577-019-0198-4,

PubMed Abstract | Crossref Full Text | Google Scholar

Uguz, A., Muftuoglu, C., Mert, U., Gumus, T., Ece, D., Asadi, M., et al. (2025). Unveiling microbiota profiles in saliva and pancreatic tissues of patients with pancreatic cancer. Microorganisms 13:119. doi: 10.3390/microorganisms13010119,

PubMed Abstract | Crossref Full Text | Google Scholar

Vilà-Quintana, L., Fort, E., Pardo, L., Albiol-Quer, M. T., Ortiz, M. R., Capdevila, M., et al. (2024). Metagenomic study reveals phage-bacterial interactome dynamics in gut and oral microbiota in pancreatic diseases. Int. J. Mol. Sci. 25:10988. doi: 10.3390/ijms252010988,

PubMed Abstract | Crossref Full Text | Google Scholar

Visscher, P. M., Wray, N. R., Zhang, Q., Sklar, P., McCarthy, M. I., Brown, M. A., et al. (2017). 10 years of GWAS discovery: biology, function, and translation. Am. J. Hum. Genet. 101, 5–22. doi: 10.1016/j.ajhg.2017.06.005,

PubMed Abstract | Crossref Full Text | Google Scholar

Vogtmann, E., Yano, Y., Zouiouich, S., Hua, X., Wan, Y., Purandare, V., et al. (2025). The human oral microbiome and risk of colorectal cancer within three prospective cohort studies in the United States. Cancer 131:e35802. doi: 10.1002/cncr.35802,

PubMed Abstract | Crossref Full Text | Google Scholar

Waagsbø, B., Svardal, A., Ueland, T., Landrø, L., Øktedalen, O., Berge, R. K., et al. (2016). Low levels of short- and medium-chain acylcarnitines in HIV-infected patients. Eur. J. Clin. Investig. 46, 408–417. doi: 10.1111/eci.12609,

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, S., Dong, W., Liu, L., Xu, M., Wang, Y., Liu, T., et al. (2019). Interplay between bile acids and the gut microbiota promotes intestinal carcinogenesis. Mol. Carcinog. 58, 1155–1167. doi: 10.1002/mc.22999,

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, Y., Yan, Q., Fan, C., Mo, Y., Wang, Y., Li, X., et al. (2023). Overview and countermeasures of cancer burden in China. Sci. China Life Sci. 66, 2515–2526. doi: 10.1007/s11427-022-2240-6,

PubMed Abstract | Crossref Full Text | Google Scholar

Wu, S. L., Zhang, Y., Fu, Y., Li, J., and Wang, J. S. (2022). Gastric cancer incidence, mortality and burden in adolescents and young adults: a time-trend analysis and comparison among China, South Korea, Japan and the USA. BMJ Open 12:e061038. doi: 10.1136/bmjopen-2022-061038,

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, C., Gu, L., Hu, L., Jiang, C., Li, Q., Sun, L., et al. (2023). FADS1-arachidonic acid axis enhances arachidonic acid metabolism by altering intestinal microecology in colorectal cancer. Nat. Commun. 14:2042. doi: 10.1038/s41467-023-37590-x,

PubMed Abstract | Crossref Full Text | Google Scholar

Zheng, Z., and Xu, X. (2025). Integrative multi-omics identification and functional validation of potential targets linking metabolism-immune-colorectal cancer causal pathway. Front. Immunol. 16:1649788. doi: 10.3389/fimmu.2025.1649788,

PubMed Abstract | Crossref Full Text | Google Scholar

Zhu, S. J., and Ding, Z. (2024). Association between gut microbiota and seven gastrointestinal diseases: a Mendelian randomized study. J. Gene Med. 26:e3623. doi: 10.1002/jgm.3623,

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: 16S rDNA sequencing, gut microbiota, mendelian randomization, metabolites, pancreatic cancer

Citation: Lin S, Shi E, Zhang Y, Wang X, Tian Z, Han J and Li Q (2026) Gut microbiota, circulating metabolites, and pancreatic cancer risk: a multi-method causal inference study with cross-population validation. Front. Microbiol. 16:1730313. doi: 10.3389/fmicb.2025.1730313

Received: 22 October 2025; Revised: 27 November 2025; Accepted: 09 December 2025;
Published: 23 January 2026.

Edited by:

Lucinda Janete Bessa, Egas Moniz Center for Interdisciplinary Research (CiiEM), Portugal

Reviewed by:

Maria Gazouli, National and Kapodistrian University of Athens, Greece
George Grant, Independent Researcher, Aberdeen, United Kingdom

Copyright © 2026 Lin, Shi, Zhang, Wang, Tian, Han and Li. 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: Quanwang Li, cXVhbndhbmdsaUAxMjYuY29t; Jing Han, aGpoZXhmMDQxMUAxNjMuY29t

Disclaimer: 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.