ORIGINAL RESEARCH article

Front. Immunol., 05 October 2023

Sec. Microbial Immunology

Volume 14 - 2023 | https://doi.org/10.3389/fimmu.2023.1260780

Genome-wide Mendelian randomization identifies putatively causal gut microbiota for multiple peptic ulcer diseases

  • 1. Department of General Surgery, Xinhua Hospital, Affiliated to Shanghai Jiao Tong University School of Medicine, Shanghai, China

  • 2. Shanghai Key Laboratory of Biliary Tract Disease Research, Shanghai, China

  • 3. Department of Cardiovascular Surgery, The First Affiliated Hospital of Soochow University, Suzhou, Jiangsu, China

  • 4. Laboratory of Medical Science, School of Medicine, Nantong University, Nantong, Jiangsu, China

  • 5. Department of General Surgery, The First Affiliated Hospital of Soochow University, Suzhou, Jiangsu, China

  • 6. Department of Gastroenterology, Wuzhong People’s Hospital of Suzhou, Suzhou, China

Abstract

Objective:

The pathogenesis of peptic ulcer diseases (PUDs) involves multiple factors, and the contribution of gut microbiota to this process remains unclear. While previous studies have associated gut microbiota with peptic ulcers, the precise nature of the relationship, whether causal or influenced by biases, requires further elucidation.

Design:

The largest meta-analysis of genome-wide association studies was conducted by the MiBioGen consortium, which provided the summary statistics of gut microbiota for implementation in the Mendelian randomization (MR) analysis. Summary statistics for five types of PUDs were compiled using the FinnGen Consortium R8 release data. Various statistical techniques, including inverse variance weighting (IVW), MR-Egger, weighted median (WM), weighted mode, and simple mode, were employed to assess the causal relationships between gut microbiota and these five PUDs.

Result:

In the intestinal microbiome of 119 known genera, we found a total of 14 causal associations with various locations of PUDs and reported the potential pathogenic bacteria of Bilophila et al. Among them, four had causal relationships with esophageal ulcer, one with gastric ulcer, three with gastroduodenal ulcer, four with duodenal ulcer, and two with gastrojejunal ulcer.

Conclusion:

In this study, the pathogenic bacterial genera in the gut microbiota that promote the occurrence of PUDs were found to be causally related. There are multiple correlations between intestinal flora and PUDs, overlapping PUDs have overlapping associated genera. The variance in ulcer-related bacterial genera across different locations underscores the potential influence of anatomical locations and physiological functions.

1 Introduction

Peptic ulcer diseases (PUDs) represent a prevalent clinical condition characterized by multifactorial etiology and extremely complex pathogenesis, primarily related to Helicobacter pylori infection (). The incidence of PUDs is common among individuals between the age of 25 and 64 years and increases with age. These ulcers are predominantly located near the stomach or duodenum but can also occur in the esophagus or Meckel’s diverticulum (). In the general population, the lifetime prevalence of PUDs is estimated to range between 5% and 10%, while the annual incidence rate ranges from 0.1% to 0.3% (). The continued prevalence of peptic ulcers within the stomach and duodenum poses a significant threat to global public health. The diagnosis and treatment of PUDs remain a major healthcare problem with a significant disease burden ().

The human gut harbors an intricate and diverse microbial community that plays a crucial role in both health and diseases (, ), for instance, digestion and absorption of substances, synthesis of essential vitamins such as B and K, catabolism of compounds in vivo, coordination of innate and cell-mediated immune responses, and maintenance of intestinal barrier function (). The symbiotic relationship between these microbes and the host is indispensable for maintaining overall homeostasis; disruptions in this ecological equilibrium can lead to adverse health outcomes (). The correlation between alterations in the gut microbiome and peptic ulcers has been studied for a long time. Several studies have demonstrated the mechanism underlying H. pylori-induced PUDs (). At the same time, histological techniques have been utilized to examine the microbiome and metabolome of gastric biopsy tissues, identifying a distinct correlation between gastrointestinal ulcers and gastrointestinal bacteria (). Moreover, PUDs were significantly associated with abnormal microbiota compositions in the oropharynx, esophagus, and gastrointestinal tract (). Therapies to protect, adapt, shape, or restore the balance of the microbiome are critical aspects of the current and prospective approaches to gastrointestinal ulcer management (). However, the causal relationship between PUDs at different anatomical sites and the gut microbiota remains unclear and requires further elucidation.

The genome-wide association study (GWAS) has gained widespread acceptance as a pivotal approach for exploring potential genetic variants linked to diverse and complex traits and diseases (, ). Mendelian randomization (MR) analysis introduces an innovative paradigm to explore the potential causal association between exposure and outcome independent of confounding factors and ethical considerations. Through MR, genetic variants are leveraged as instrumental variables (IVs) for exposure, enabling the estimation of causality between the exposure and the resultant outcome (, ). An MR study mimics a randomized controlled trial (RCT), as genetic variations are randomly assigned during fertilization (). Furthermore, genotype formation occurs prior to disease onset and is typically unaffected by disease progression, reducing the likelihood of confounding influences.

Here, we employed MR analysis to investigate the correlation between gut microbiota and PUDs. We further explored the potential therapeutic implications of selective support or disorder of the gut microbiota.

2 Methods

2.1 Study design

The basic logic and analysis flow of the entire procedure were briefly described in Figure 1. The causal effects of the gut microbiota on five PUDs, including esophageal ulcer (OESU), gastric ulcer (GU), gastroduodenal ulcer (GASTRODU), duodenal ulcer (DU), and gastrojejunal ulcer (GJU), were evaluated. To comprehensively investigate the role of the gut microbiota in PUDs, MR analysis was performed at the classification of genera. The population information involved in the MR was detailed in Table 1.

Figure 1

Table 1

Exposure/OutcomeEthnic originSample size (case/control)GenderRegistry filter (ICD-10)Public releaseData source
Gut microbiota
(genus)
European, Hispanic, Middle Eastern,
Asian and African
18,340Mixed sex—2018MiBioGen consortium
OESUEurope1,840/292,256Mixed sexK22.12022FinnGen R8
GUEurope5,277/292,256Mixed sexK252022FinnGen R8
GASTRODUEurope8,240/292,256Mixed sexK2[5-8]2022FinnGen R8
DUEurope3,164/292,256Mixed sexK262022FinnGen R8
GJUEurope288/292,256Mixed sexK282022FinnGen R8

The population information involved in this Mendelian randomization.

OESU, esophageal ulcer; GU, gastric ulcer; GASTRODU, gastroduodenal ulcer; DU, duodenal ulcer; GJU, gastrojejunal ulcer; IVW, inverse variance weighted; SNP, single nucleotide polymorphism.

2.2 Data sources

A two-sample MR study was undertaken to explore the potential relationship between genus-level gut microbiota and PUDs, utilizing GWAS summary data. Studies received prior approval from their respective institutional review boards (IRBs), and informed consent was obtained from all participants and/or their legal guardians.

To obtain GWAS summary statistics for the gut microbiota, data from the MiBioGen consortium, the largest GWAS dataset published to date, were utilized (). This dataset consisted of 18,340 individuals spanning 24 population-based cohorts of diverse ancestry, including European, Middle Eastern, East Asian, American Hispanic/Latino, and American African. Microbial composition profiling and taxonomic classification were performed using direct taxonomic binning, targeting variable regions V4, V3–V4, and V1–V2 of the 16S rRNA gene. A microbiota quantitative trait loci (mbQTL) mapping analysis was employed to identify host genetic variants associated with the abundance of bacterial taxa within the gut microbiota. The genus with the least classification in the GWAS data for gut microbiota was selected for preprocessing. Out of 131 identified genera with an average abundance surpassing 1%, 119 genera were included for analysis, while 12 genera remain uncharacterized.

For the GWAS summary statistics of the five peptic ulcer types, data were obtained from the FinnGen consortium R8 release (https://www.finngen.fi/fi). The FinnGen consortium is a large public–private partnership aiming to collect and analyze genomic and health data from 500,000 Finnish biobank participants (). The dataset available for analysis was up to December 2022. The peptic ulcer types included were OESU (NCase = 1,840, NControl = 292,256), GU (NCase = 5,277, NControl = 292,256), GASTRODU (NCase = 8,240, NControl = 292,256), DU (NCase = 3,164, NControl = 292,256), and GJU (NCase = 288, NControl = 292,256). The classification of OESU, GU, GASTRODU, DU, and GJU adhered strictly to the guidelines outlined by the International Classification of Diseases 10th Revision (ICD-10) code.

2.3 Instrumental variable selection

The process of IV selection from the GWAS summary statistics of the gut microbiome adhered to the following criteria: 1) Single-nucleotide polymorphisms (SNPs) associated with each genus at the locus-wide significance threshold (P < 1e–5) were considered as potential IVs (); 2) Using reference panel data from the 1000 Genomes project European samples, linkage disequilibrium (LD) between SNPs was calculated. Among the SNPs with R2 < 0.001 (clumping window size = 10,000 kb), only those with the lowest P-values were retained to minimize biased genetic variation arising from residual LD; 3) The intensity of each IV and exclusion of weak instruments were evaluated by calculating the F-statistic (F > 10); 4) In the presence of palindromic SNPs, the alleles on the forward strand were deduced by utilizing information on allele frequencies; 5) In cases where exposure-associated SNPs were absent in outcome data, suitable proxy SNPs (r2 > 0.8) were identified and included in subsequent analyses; and 6) To address confounding, SNPs related to H. pylori infection, bile reflux, obesity, alcoholism, smoking, and stress factors were systematically removed during the MR analysis.

2.4 Statistical analysis

Five popular MR methods were utilized to analyze valid IVs: inverse variance-weighted (IVW) test, MR-Egger regression, weighted median, weighted mode, and simple mode (Figure 2). Among these, IVW was predominantly used due to its slightly higher power under certain conditions (). The IVW method utilized the inverse of the outcome variance as weights for fitting, regardless of the presence of an intercept term in the regression. Complementary assessments were performed using the remaining four methods, each of which was based on different assumptions about potential pleiotropy. If the results obtained by these complementary methods are consistent with the IVW estimation results, the robustness of the effect estimation can be reinforced.

Figure 2

Multiple methods of sensitivity analyses were undertaken to ensure robustness. Initially, Cochran’s Q statistics was applied to assess heterogeneity across diverse studies (). Statistically significant Cochran’s Q-test would indicate significant heterogeneity in the analytical outcomes. Secondly, MR-Pleiotropy Residual Sum and Outlier (MR-PRESSO) was used to detect instances of horizontal pleiotropy, with SNPs demonstrating horizontal pleiotropy outliers being systematically excluded to minimize pleiotropy-induced effects (). In cases where significant horizontal pleiotropy was detected in the MR-PRESSO global test, outliers with P < 0.05 were removed, and the remaining SNPs were reanalyzed with the IVW analysis. Thirdly, the MR-Egger regression intercept was employed to estimate the potential pleiotropy of SNPs, with a P-value >0.05 indicating no horizontal pleiotropy (). Fourthly, a leave-one-out analysis was performed to assess the impact of each SNP on the causal signal. Finally, funnel and forest plots were constructed to visually examine the presence of horizontal pleiotropy in the MR analysis, with P < 0.05 indicating potential causal associations. The statistical analyses were carried out using the R packages: two-sample MR () and MR-PRESSO ().

3 Results

3.1 Genetic correlations between gut microbiota and PUDs

Initially, 1,698 SNPs were screened as possible IVs for 131 bacterial genera including 12 unidentified genera. The genetic variants were then eliminated based on specific criteria. All F-statistic exceeded 10, suggesting the absence of weak ins. Following validation through the PhenoScanner database, the remaining SNPs exhibited no discernible associations with H. pylori infection, bile reflux, obesity, alcoholism, smoking, and stress factors, indicating that IVs were not resolved by confounding factors. Concurrently, the elimination of palindromic SNPs was performed. GWAS data for patients with PUDs were derived from corresponding cohorts.

3.2 Bidirectional causal relationship of gut bacteria on PUD development

Our results demonstrated an association of four bacterial genera with OESU, one with GU, three with GASTRODU, four with DU, and two with GJU. Reverse MR analysis demonstrated that peptic ulcers did not change the abundance of the above bacteria. The leave-one-out sensitivity test highlighted some continuity around the midpoint. Crossing the zero line indicated that the result may be insignificant or unstable. The overall assessment indicated an absence of SNPs having a dominant impact. Furthermore, the selected SNPs exhibited no significant heterogeneity, as indicated by Cochran’s Q statistics. Leave-one-out analysis did not identify a single SNP driving the association (Supplementary Figure S1). The application of MR-PRESSO yielded no outliers. The findings from the MR-Egger regression intercept analysis further corroborated the absence of significant directional horizontal pleiotropy (Table 2). The three main statistical results of the MR analysis were shown in Supplementary Table S1. Given the absence of significant statistical difference in reverse causality, relevant results were presented in Supplementary Table S2.

Table 2

ExposureOutcomeSNP (n)IVWCochran’s Q(MR-Egger)Pleiotropy_testF-statistic (median)
OR
(95% CI)
P-valueQQ_pvalEgger interceptSeP-value
Eubacterium halliiOESU140.71 (0.53-0.95)0.0245.9350.920.0280.0260.3021.10
FlavonifractorOESU51.69 (1.08-2.64)0.0201.3490.720.0280.0740.7321.93
Ruminiclostridium 6OESU151.39 (1.03-1.88)0.0308.6770.80-0.0270.0330.4220.91
Ruminococcaceae UCG013OESU111.82 (1.27-2.61)0.0017.7700.56-0.0600.0400.1721.52
Lachnospiraceae UCG004GU121.34 (1.09-1.65)0.0065.9000.82-0.0040.0280.8921.26
Lachnospiraceae FCS020GASTRODU120.85 (0.73-0.99)0.0402.9790.98-0.0020.0150.9221.66
Lachnospiraceae UCG004GASTRODU121.19 (1.00-1.40)0.0488.2170.610.0100.0230.6821.26
Ruminiclostridium 9GASTRODU80.77 (0.61-0.96)0.0192.7620.840.0070.0350.8421.46
CatenibacteriumDU41.31 (1.05-1.63)0.0182.3950.30-0.0700.1990.7621.28
Clostridium sensu stricto 1DU60.65 (0.42-1.00)0.0489.8150.040.0310.0620.6520.32
CollinsellaDU90.69 (0.50-0.95)0.0243.5490.83-0.0110.0430.8120.78
Ruminiclostridium 9DU80.68 (0.48-0.97)0.0313.5550.740.0020.0550.9721.46
ParabacteroidesGJU50.22 (0.06-0.84)0.0271.5900.660.0250.3740.9521.55
BilophilaGJU133.45 (1.52-7.81)0.0037.2300.78-0.0640.1500.6821.02

MR estimates for the association between gut microbiota and PUDs.

OESU, esophageal ulcer; GU, gastric ulcer; GASTRODU, gastroduodenal ulcer; DU, duodenal ulcer; GJU, gastrojejunal ulcer; IVW, inverse variance weighted; SNP, single nucleotide polymorphism.

3.2.1 Causal relationship of gut bacteria on OESU

Notably, the two-sample MR analysis unveiled a causal linkage between Eubacterium hallii and OESU [IVW odds ratio (OR) = 0.71, 95% CI: 0.53–0.95, P = 0.024]. Furthermore, three bacterial features exhibited potential associations with an increased OESU risk: Flavonifractor (IVW OR = 1.69, 95% CI: 1.08–2.64, P = 0.020), Ruminiclostridium 6 (IVW OR = 1.39, 95% CI: 1.03–1.88, P = 0.030), and Ruminococcaceae UCG013 (IVW OR = 1.82, 95% CI: 1.27–2.61, P = 0.001).

3.2.2 Causal relationship of gut bacteria on GU

Our results underscored a robust causal relationship between Lachnospiraceae UCG004 and GU (IVW OR = 1.34, 95% CI: 1.09–1.65, P = 0.006).

3.2.3 Causal relationship of gut bacteria on GASTRODU

We found that Lachnospiraceae UCG004 (IVW OR = 1.19, 95% CI: 1.00–1.40, P = 0.048) was associated with an increased GASTRODU risk. Moreover, Lachnospiraceae FCS020 (IVW OR = 0.85, 95% CI: 0.73–0.99, P = 0.040) and Ruminiclostridium 9 (IVW OR = 0.77, 95% CI: 0.61–0.96, P = 0.019) were associated with a lower GASTRODU risk.

3.2.4 Causal relationship of gut bacteria on DU

We found potential associations between one bacterial feature, Catenibacterium (IVW OR = 1.31, 95% CI: 1.05–1.63, P = 0.018), and increased DU risk. Meanwhile, three bacterial features, Clostridium sensu stricto 1 (IVW OR = 0.65, 95% CI: 0.42–1.00, P = 0.048), Collinsella (IVW OR = 0.69, 95% CI: 0.50–0.95, P = 0.024), and Ruminiclostridium 9 (IVW OR = 0.68, 95% CI: 0.48–0.97, P = 0.031), were associated with a reduced DU risk.

3.2.5 Causal relationship of gut bacteria on GJU

Furthermore, we found an association between Parabacteroides (IVW OR = 0.22, 95% CI: 0.06–0.84, P = 0.027) and a lower GJU risk, while Bilophila (IVW OR = 3.45, 95% CI: 1.52–7.81, P = 0.003) showed a correlation with an elevated risk of GJU.

4 Discussion

Prior to the recognition of H. pylori infection and the extensive utilization of nonsteroidal anti-inflammatory drugs (NSAIDs) during the latter time frame of the 20th century, PUDs were primarily attributed to a hypersecretory acidic environment, coupled with dietary factors or stress (, ). However, there is growing recognition that the etiology of PUDs extends beyond H. pylori infection in the stomach (). Rather, the genesis and progression of PUDs emerge as a result of the interplay of multiple factors, encompassing the presence of different H. pylori virulence proteins, ensuing human immune reactions, and imbalances in the gastrointestinal microbiota (, ). The role of intrinsic gut bacteria in PUD development is also noteworthy. Our findings emphasized the causal involvement of specific bacterial characteristics’ abundance in modulating the susceptibility to diverse peptic ulcer types. Remarkably, this study represents the first MR analysis to illuminate the multiple connections between gut microbiota and PUDs. We regarded it as a longitudinal microbiome investigation conducted before the onset of PUDs in humans. This study effectively identified robustly gene variants through the largest gut microbiome GWAS.

One pivotal role played by gut microbiota involves the synthesis of short-chain fatty acids (SCFAs), which can directly regulate host health through energy regulation, intestinal mucosal barrier, immune regulation, and induction of tumor cell differentiation and apoptosis (, ). Dysregulation in the equilibrium of SCFAs within the body results in a cascade of disease manifestations (28). Moreover, SCFAs can promote the expression of tight junction proteins, such as claudin, occludin, and Zonula occludens (ZOs) within the intestinal tract, decrease intestinal permeability, promote the proliferation of intestinal mucosal cells, and improve the mechanical barrier function of the intestine in animal models (29, 30). In this MR study, SCFA-producing bacteria included E. hallii (31), Flavonifractor (32), Ruminiclostridium (33), Ruminococcaceae (33), Collinsella (34), and Parabacteroides (35).

A previous study revealed reduced levels of Collinsella in patients with inflammatory bowel disease (IBD) and gut microbiota dysbiosis (36). This finding was consistent with our results, implicating a potentially beneficial role for Collinsella in gut health and highlighting its association with a reduced DU risk. Researchers have found that Parabacteroides could produce a molecule named rhamnose in the mouse gut to facilitate the repair and maintenance of the intestinal mucosal barrier in mice (37). This suggested that Parabacteroides may benefit gut health, aligning with our findings that it exhibited a negative association with GJU. E. hallii is a high-yielding butyrate producer in the gut, contributing significantly to the maintenance of intestinal metabolic equilibrium (38). In the context of aging populations characterized by a decrease in microbiota diversity, a reduction in the abundance of E. hallii has been noted, accompanied by decreased production of SCFAs and increased intestinal inflammation (39). Consistent with these studies, our results demonstrated a negative link with OESU. The presence of Eubacterium in the gut is primarily associated with increased dietary fiber intake. As previously reported, it may improve the intestinal mucosal barrier and metabolic diseases, making it a potential candidate strain for a new generation of probiotics (40).

Previous studies have shown an increased level of Flavonifractor in patients with early-onset colon cancer, while Lachnospiraceae UCG004 was significantly increased in patients with postmenopausal osteoporosis (PMO). While the adverse effects of Flavonifractor remain relatively underexplored, several studies hint at its role in stabilizing gut intestinal flora and immune modulation (41). Our study uncovered a positive correlation between Flavonifractor and OESU, hinting at its potential as a risk factor. Ruminiclostridium, a common anaerobic intestinal bacterium, plays a pivotal role in polysaccharide degradation and SCFA production, thereby influencing intestinal peristalsis, intestinal health, and immune modulation (42). In addition, Ruminiclostridium 9 inhibits the growth of other harmful bacteria, which is crucial for preserving gut microbiota equilibrium (43). Recent studies have shown that reduced Ruminiclostridium 9 abundance is also associated with some intestinal diseases, such as IBD and obesity (). Our MR study demonstrated it with the negative causal relationship between GASTRODU or DU, thus highlighting the potential protective role of Ruminiclostridium 9 during the development of PUDs. Despite the positive correlation between Ruminiclostrium 6 and OESU, we did not observe any significant negative effect in the gut. Ruminococcaceae and Lachnospiraceae are typical intestinal flora that are important in maintaining intestinal health (44). However, some studies have shown the increased abundance of Ruminococcaceae or Lachnospiraceae in metabolic disorders such as obesity and diabetes (45). Additionally, certain members of these genera have been associated with the production of inflammatory mediators, enterotoxins, and other harmful substances related to the occurrence and development of intestinal diseases (46, 47). Our study found that Ruminococcaceae UCG013 was positively associated with OESU, while Lachnospiraceae UCG004 was positively correlated with GU and GASTRODU, suggesting their potential role as risk factors. Intriguingly, Lachnospiraceae FCS020 exhibited a negative causal link with OESU.

C. sensu stricto, a beneficial intestinal bacterium, has many vital physiological and metabolic functions, such as participating in the metabolism of glucose and lactose as well as promoting the synthesis of biotin and vitamin K (48). C. sensu stricto can also promote the integrity of the intestinal mucosal barrier, regulate the intestinal immune response, and reduce intestinal inflammation (48). Our MR analysis corroborates its protective effect, revealing a negative correlation between C. sensu stricto 1 and DU. Conversely, Catenibacterium shows a positive causal association with DU. Catenibacterium is a Gram-positive bacterium (49). Although many species of Catenibacterium are unknown, different genera play different roles in intestinal diseases. For instance, certain strains have been implicated in the occurrence and development of IBD, wherein harmful substances, such as enterotoxins, contribute to intestinal mucosal barrier disruption and aggravated inflammatory responses (50). Our findings implicate Bilophila as a potential risk factor for GJU, which used to be mainly associated with metabolic diseases (51). Bilophila’s role in the gut requires further understanding; however, studies suggest that some members of Bilophila may be involved in the occurrence and development of intestinal inflammation. Some of these strains can produce harmful substances, such as hydrogen sulfide, breaking the intestinal mucosal barrier and increasing inflammatory responses (52). Notably, PUDs exhibited no significant association with the aforementioned bacteria in reverse MR analysis.

The variance in ulcer-related bacterial genera across different locations underscores the potential influence of anatomical locations and physiological functions. Interestingly, ulcers with overlapping sites, such as gastroduodenal ulcer, gastric ulcer, and duodenal ulcer, exhibited similar bacterial flora and associations, confirming the reliability of the results. Importantly, these findings prompted an exploration into the pathological mechanism underlying overlapping bacterial flora in ulcer development. Furthermore, there was no evidence of reverse causality between PUDs and gut bacterial genera. Overall, gene-based analysis from 119 gut bacterial genera revealed specific genera associated with PUDs in different locations and explained the multiple correlations between them. These findings supported the influence of gut microbiota on the development of PUDs and highlighted the putative association between specific bacterial genera and site-specific PUDs. Ultimately, these findings extended valuable implications for the clinical management of patients afflicted with PUDs.

5 Article summary

This study presented a comprehensive analysis of the causal relationships between 121 known gut bacterial genera and PUDs, utilizing both forward and reverse MR analyses. This approach effectively mitigated the influence of confounding variables and causal inference’s challenge of reverse causation. Notably, the genetic variants associated with gut microbiota were derived from the most extensive GWAS meta-analysis, enhancing robust IVs for the MR analysis. Multiple statistical methods were used to test the sensitivity, pleiotropy, and heterogeneity of this study.

However, despite our efforts to minimize confounding influences, the complete elimination of horizontal pleiotropy remains a challenge, largely attributed to our limited understanding of the disease. As knowledge and awareness evolve over time, perceptions about confounding factors may change. In the future, extending MR investigations on the causal relationship between gut microbiota and peptic ulcers in diverse European and non-European populations would enhance the generalization of our findings. Additionally, while this article only explores the problem from the perspective of genetics, higher-level RCTs are necessary to validate the causal relationship while also into the intricate mechanisms underlying specific bacterial contributions.

In summation, through a systematic investigation, this study constituted a pioneering MR analysis focused on gut microbiota and PUDs. Our findings shed light on the multiple correlations between gut microbiota genera and five PUD types. Moreover, it has established definitive links between 14 specific bacterial genera and their corresponding ulcer manifestations, the pathogenic intestinal bacteria deserved more attention. These findings hold significant implications for understanding the role of gut flora in PUDs, offering valuable insights for the formulation of preventive strategies in patients with this condition. Novel treatment avenues targeting specific intestinal bacterial genera may represent new treatment options for PUDs.

Statements

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: MiBioGen Consortium, https://mibiogen.gcc.rug.nl/ and FinnGen Consortium R8, https://r8.finngen.fi/.

Ethics statement

All studies were previously approved by respective Institutional Review Boards (IRBs). No new IRB approval was required. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

JZ: Conceptualization, Data curation, Formal Analysis, Methodology, Software, Validation, Writing – original draft, Writing – review & editing. YH: Data curation, Software, Writing – original draft, Writing – review & editing. YiZ: Conceptualization, Writing – original draft, Writing – review & editing. TX: Conceptualization, Data curation, Software, Writing – original draft. XF: Data curation, Formal Analysis, Software, Writing – original draft. ZY: Conceptualization, Investigation, Supervision, Visualization, Writing – original draft. YaZ: Funding acquisition, Investigation, Resources, Software, Supervision, Writing – review & editing. WG: Funding acquisition, Resources, Supervision, Writing – review & editing.

Funding

The authors declare financial support was received for the research, authorship, and/or publication of this article. This research was funded by the National Natural Science Foundation of China (No. 82172628, 81974371); Xinhua Hospital Funded Clinical Research (21XHDB10) and Shanghai Key Laboratory of Biliary Tract Disease Research Foundation (No.17DZ2260200).

Acknowledgments

We want to acknowledge the participants and investigators of the FinnGen study and MiBioGen consortium. All authors are thanked for their contributions and Bullet Edits Limited for the linguistic editing and proofreading of the manuscript. Support from the National Natural Science Foundation of China and Xinhua Hospital is gratuitous.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1260780/full#supplementary-material

Supplementary Figure 1

Sensitivity test performed by leave-one-out analysis.

References

  • 1

    LanasAChanFKL. Peptic ulcer disease. Lancet (2017) 390(10094):613–24. doi: 10.1016/S0140-6736(16)32404-7

  • 2

    RenJJinXLiJLiRGaoYZhangJet al. The global burden of peptic ulcer disease in 204 countries and territories from 1990 to 2019: a systematic analysis for the Global Burden of Disease Study 2019. Int J Epidemiol (2022) 51(5):1666–76. doi: 10.1093/ije/dyac033

  • 3

    LauJYSungJHillCHendersonCHowdenCWMetzDC. Systematic review of the epidemiology of complicated peptic ulcer disease: incidence, recurrence, risk factors and mortality. Digestion (2011) 84(2):102–13. doi: 10.1159/000323958

  • 4

    LozuponeCAStombaughJIGordonJIJanssonJKKnightR. Diversity, stability and resilience of the human gut microbiota. Nature (2012) 489(7415):220–30. doi: 10.1038/nature11550

  • 5

    RoundJLMazmanianSK. The gut microbiota shapes intestinal immune responses during health and disease. Nat Rev Immunol (2009) 9(5):313–23. doi: 10.1038/nri2515

  • 6

    WangCYuXLinHQWangGQLiuJMGaoCCet al. Integrating microbiome and metabolome revealed microbe-metabolism interactions in the stomach of patients with different severity of peptic ulcer disease. Front Immunol (2023) 14. doi: 10.3389/fimmu.2023.1134369

  • 7

    HuangYWangZMaHJiSChenZCuiZet al. Dysbiosis and implication of the gut microbiota in diabetic retinopathy. Front Cell Infect Microbiol (2021) 11. doi: 10.3389/fcimb.2021.646348

  • 8

    NIH Consensus Conference. Helicobacter pylori in peptic ulcer disease. NIH consensus development panel on helicobacter pylori in peptic ulcer disease. Jama (1994) 272(1):65–9. doi: 10.1001/jama.1994.03520010077036

  • 9

    GrahamDY. Campylobacter pylori and peptic ulcer disease. Gastroenterology (1989) 96(2 Pt 2 Suppl):615–25. doi: 10.1016/S0016-5085(89)80057-5

  • 10

    BelizarioJENapolitanoM. Human microbiomes and their roles in dysbiosis, common diseases, and novel therapeutic approaches. Front Microbiol (2015) 6. doi: 10.3389/fmicb.2015.01050

  • 11

    HowieBNDonnellyPMarchiniJ. A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. PloS Genet (2009) 5(6):e1000529. doi: 10.1371/journal.pgen.1000529

  • 12

    KurkiMIKarjalainenJPaltaPSipilaTPKristianssonKDonnerKMet al. FinnGen provides genetic insights from a well-phenotyped isolated population. Nature (2023) 613(7944):508. doi: 10.1038/s41586-022-05473-8

  • 13

    EmdinCAKheraAVKathiresanS. Mendelian randomization. Jama J Am Med Assoc (2017) 318(19):1925–6. doi: 10.1001/jama.2017.17219

  • 14

    GaoKCaoLFMaWZGaoYJLuoMSZhuJet al. Association between sarcopenia and cardiovascular disease among middle-aged and older adults: Findings from the China health and retirement longitudinal study. Eclinicalmedicine (2022) 44:101264. doi: 10.1016/j.eclinm.2021.101264

  • 15

    GreenlandS. An introduction to instrumental variables for epidemiologists. Int J Epidemiol (2000) 29(4):722–9. doi: 10.1093/ije/29.4.722

  • 16

    KurilshikovAMedina-GomezCBacigalupeRRadjabzadehDWangJDemirkanAet al. Large-scale association analyses identify host factors influencing human gut microbiome composition. Nat Genet (2021) 53(2):156–65. doi: 10.1038/s41588-020-00763-1

  • 17

    SannaSvan ZuydamNRMahajanAKurilshikovAVilaAVVosaUet al. Causal relationships among the gut microbiome, short-chain fatty acids and metabolic diseases. Nat Genet (2019) 51(4):600. doi: 10.1038/s41588-019-0350-x

  • 18

    BoehmFJZhouX. Statistical methods for Mendelian randomization in genome-wide association studies: A review. Comput Struct Biotechnol J (2022) 20:2338–51. doi: 10.1016/j.csbj.2022.05.015

  • 19

    Del GrecoMFMinelliCSheehancNAThompsoncJR. Detecting pleiotropy in Mendelian randomisation studies with summary data and a continuous outcome. Stat Med (2015) 34(21):2926–40. doi: 10.1002/sim.6522

  • 20

    VerbanckMChenC-YNealeBDoR. Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet (2018) 50(5):693. doi: 10.1038/s41588-018-0099-7

  • 21

    BurgessSThompsonSG. Interpreting findings from Mendelian randomization using the MR-Egger method. Eur J Epidemiol (2017) 32(5):377–89. doi: 10.1007/s10654-017-0255-x

  • 22

    HemaniGZhengnJElsworthBWadeKHHaberlandVBairdDet al. The MR-Base platform supports systematic causal inference across the human phenome. Elife (2018) 7:e34408. doi: 10.7554/eLife.34408

  • 23

    SuerbaumSMichettiP. Helicobacter pylori infection. N Engl J Med (2002) 347(15):1175–86. doi: 10.1056/NEJMra020542

  • 24

    PetersonWL. Helicobacter pylori and peptic ulcer disease. N Engl J Med (1991) 324(15):1043–8.

  • 25

    HuWDengCMaZQWangDJFanCXLiTet al. Utilizing melatonin to combat bacterial infections and septic injury. Br J Pharmacol (2017) 174(9):754–68. doi: 10.1111/bph.13751

  • 26

    SmithPMHowittMRPanikovNMichaudMGalliniCABohlooly-YMet al. The microbial metabolites, short-chain fatty acids, regulate colonic T-reg cell homeostasis. Science (2013) 341(6145):569–73. doi: 10.1126/science.1241165

  • 27

    MorrisonDJPrestonT. Formation of short chain fatty acids by the gut microbiota and their impact on human metabolism. Gut Microbes (2016) 7(3):189–200. doi: 10.1080/19490976.2015.1134082

  • 28

    TanCWuQWangHGaoXXuRCuiZet al. Dysbiosis of gut microbiota and short-chain fatty acids in acute ischemic stroke and the subsequent risk for poor functional outcomes. J Parenteral Enteral Nutr (2021) 45(3):518–29. doi: 10.1002/jpen.1861

  • 29

    den BestenGvan EunenKGroenAKVenemaKReijngoudD-JBakkerBM. The role of short-chain fatty acids in the interplay between diet, gut microbiota, and host energy metabolism. J Lipid Res (2013) 54(9):2325–40. doi: 10.1194/jlr.R036012

  • 30

    HamerHMJonkersDMAEBastAVanhoutvinSALWFischerMAJGKoddeAet al. Butyrate modulates oxidative stress in the colonic mucosa of healthy humans. Clin Nutr (2009) 28(1):88–93. doi: 10.1016/j.clnu.2008.11.002

  • 31

    EngelsCRuscheweyhH-JBeerenwinkelNLacroixCSchwabC. The Common Gut Microbe Eubacterium hallii also Contributes to Intestinal Propionate Formation. Front Microbiol (2016) 7. doi: 10.3389/fmicb.2016.00713

  • 32

    RosesCCuevas-SierraAQuintanaSRiezu-BojJIMartinezJAMilagroFIet al. Gut microbiota bacterial species associated with mediterranean diet-related food groups in a Northern Spanish population. Nutrients (2021) 13(2):636. doi: 10.3390/nu13020636

  • 33

    LiuDLiTZhengHYinXChenMLiaoZet al. Study on alterations of physiological functions in aged constipation rats with fluid-deficiency based on metabonomic and microbiology analysis. Rsc Adv (2017) 7(76):48136–50. doi: 10.1039/C7RA07651G

  • 34

    ChangYChenYZhouQWangCChenLDiWet al. Short-chain fatty acids accompanying changes in the gut microbiome contribute to he development o hypertension in patients with preeclampsia. Clin Sci (2020) 134(2):289–302. doi: 10.1042/CS20191253

  • 35

    WangKLiaoMZhouNBaoLMaKZhengZet al. Parabacteroides distasonis Alleviates Obesity and Metabolic Dysfunctions via Production of Succinate and Secondary Bile Acids. Cell Rep (2019) 26(1):222. doi: 10.1016/j.celrep.2018.12.028

  • 36

    GeversDKugathasanSDensonLAVazquez-BaezaYVan TreurenWRenBYet al. The treatment-naive microbiome in new-onset Crohn's disease. Cell Host Microbe (2014) 15(3):382–92. doi: 10.1016/j.chom.2014.02.005

  • 37

    SonnenburgJLXuJLeipDDChenCHWestoverBPWeatherfordJet al. Glycan foraging in vivo by an intestine-adapted bacterial symbiont. Science (2005) 307(5717):1955–9. doi: 10.1126/science.1109051

  • 38

    RussellWRHoylesLFlintHJDumasM-E. Colonic bacterial metabolites and human health. Curr Opin Microbiol (2013) 16(3):246–54. doi: 10.1016/j.mib.2013.07.002

  • 39

    LiYMaZQJiangSAHuWLiTDiSYet al. A global perspective on FOXO1 in lipid metabolism and lipid-related diseases. Prog Lipid Res (2017) 66:42–9. doi: 10.1016/j.plipres.2017.04.002

  • 40

    MocoSMartinFPJRezziS. Metabolomics view on gut microbiome modulation by polyphenol-rich foods. J Proteome Res (2012) 11(10):4781–90. doi: 10.1021/pr300581s

  • 41

    SanchezB. Bile acid-microbiota crosstalk in gastrointestinal inflammation and carcinogenesis: a role for bifidobacteria and lactobacilli? Nat Rev Gastroenterol Hepatol (2018) 15(4):205. doi: 10.1038/nrgastro.2018.23

  • 42

    FlintHJScottKPDuncanSHLouisPForanoE. Microbial degradation of complex carbohydrates in the gut. Gut Microbes (2012) 3(4):289–306. doi: 10.4161/gmic.19897

  • 43

    Gutierrez-CalabresEOrtega-HernandezAModregoJGomez-GordoRCaro-VadilloARodriguez-BobadaCet al. Gut microbiota profile identifies transition from compensated cardiac hypertrophy to heart failure in hypertensive rats. Hypertension (2020) 76(5):1545–54. doi: 10.1161/HYPERTENSIONAHA.120.15123

  • 44

    LouisPYoungPHoltropGFlintHJ. Diversity of human colonic butyrate-producing bacteria revealed by analysis of the butyryl-CoA:acetate CoA-transferase gene. Environ Microbiol (2010) 12(2):304–14. doi: 10.1111/j.1462-2920.2009.02066.x

  • 45

    QinJJLiYRCaiZMLiSHZhuJFZhangFet al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature (2012) 490(7418):55–60. doi: 10.1038/nature11450

  • 46

    BajajJSRidlonJMHylemonPBThackerLRHeumanDMSmithSet al. Linkage of gut microbiome with cognition in hepatic encephalopathy. Am J Physiolol Gastrointest Liver Physiol (2012) 302(1):G168–G75. doi: 10.1152/ajpgi.00190.2011

  • 47

    LepagePHaeslerRSpehlmannMERehmanAZvirblieneABegunAet al. Twin study indicates loss of interaction between microbiota and mucosa of patients with ulcerative colitis. Gastroenterology (2011) 141(1):227–36. doi: 10.1053/j.gastro.2011.04.011

  • 48

    ZhaoFFengJLiJZhaoLLiuYChenHet al. Alterations of the gut microbiota in Hashimoto's thyroiditis patients. Thyroid (2018) 28(2):175–86. doi: 10.1089/thy.2017.0395

  • 49

    KageyamaABennoY. Catenibacterium mitsuokai gen. nov., sp, nov., a Gram-positive anaerobic bacterium isolated from human faeces. Int J Syst Evol Microbiol (2000) 50:1595–9. doi: 10.1099/00207713-50-4-1595

  • 50

    Al-AmrahHSaadahOMosliMAnneseVAl-HindiREdrisSet al. Composition of the gut microbiota in patients with inflammatory bowel disease in Saudi Arabia: A pilot study. Saudi J Gastroenterol (2023) 29(2):102–10. doi: 10.4103/sjg.sjg_368_22

  • 51

    NatividadJMLamasBHang PhuongPMichelM-LRainteauDBridonneauCet al. Bilophila wadsworthia aggravates high fat diet induced metabolic dysfunctions in mice. Nat Commun (2018) 9(1):2802. doi: 10.1038/s41467-018-05249-7

  • 52

    FengZLongWHaoBDingDMaXZhaoLet al. A human stool-derived Bilophila wadsworthia strain caused systemic inflammation in specific-pathogen-free mice. Gut Pathog (2017) 9:59. doi: 10.1186/s13099-017-0208-7

Summary

Keywords

Mendelian randomization, gut microbiota, peptic ulcer, causal relationship, genus

Citation

Zhao J, Hou Y, Xie T, Zhu Y, Feng X, Zhang Y, Yang Z and Gong W (2023) Genome-wide Mendelian randomization identifies putatively causal gut microbiota for multiple peptic ulcer diseases. Front. Immunol. 14:1260780. doi: 10.3389/fimmu.2023.1260780

Received

18 July 2023

Accepted

04 September 2023

Published

05 October 2023

Volume

14 - 2023

Edited by

Ding Shi, Zhejiang University, China

Reviewed by

Ziheng Wang, University of Macau, China; Tian Li, Independent Researcher, Xi’an, China

Updates

Copyright

*Correspondence: Yong Zhang, ; Ziyi Yang, ; Wei Gong,

†These authors have contributed equally to this work and share first authorship

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics