Yi Liu
Zhihan Xiao
Kun Ye1,2†- 1Department of Digestive System, Anqing Municipal Hospital, Anqing, China
- 2Department of Digestive System, Wannan Medical College, Wuhu, China
- 3Department of Thoracic Surgery, The First Affiliated Hospital of Nanjing Medical University, Nanjing, China
Background: Observational evidence has shown that smoking, alcohol consumption, type 2 diabetes, and body mass index (BMI) are risk factors for peptic ulcer disease (PUD), including gastric ulcer (GU) and duodenal ulcer (DU). However, the observed associations may be confounding factors. Herein, we use Mendelian randomization (MR) to examine causal associations such as smoking, alcohol, type 2 diabetes, BMI, and risks of PUD.
Methods: We used 8,17,41,325,82, 231, and 616 identified genetic variants as proxies for age of smoking initiation (AgeSmk), smoking cessation (SmkCes, current/former), number of cigarettes smoked per day (CigDay), smoking status (SmkIni, ever/never), alcohol consumption, type 2 diabetes, and BMI to obtain unconfounded effect estimates on the GU and DU levels among 452,264 participants from the Gene ATLAS. The causal relationship was estimated by using inverse-variance weighted (IVW) as the main method. Sensitivity analysis includes Cochran’s Q test, the MR-Egger test, MR pleiotropy residual sum and outlier (MR-PRESSO), and MR-robust adjusted profile score (MR-RAPS). In addition, secondary MR analysis was conducted within summary data using genetic risk scores (GRSs) as instrumental variables (IVs).
Results: In our two-sample MR analyses, genetic predisposition to smoking (SmkInit) and BMI were associated with an increased risk of GU. The beta values were 0.0035 (95% CI, 0.0021, 0.0049, p = 1.56E-06) for smoking (SmkInit) and 0.0021 (95% CI, 0.0009, 0.0033, p = 0.0008) for BMI. Genetic predisposition to smoking (SmkInit) and higher genetically predicted BMI were associated with an increased risk of DU. The beta values of DU were 0.0029 (95% CI, 0.0017, 0.0041, p = 2.43E-06) for smoking (SmkInit) and 0.0018 (95% CI, 0.0007, 0.0029, p = 0.001) for BMI. No other causal association between smoking (AgeSmk, CigDay, and SmkCes), alcohol consumption, type 2 diabetes, and GU or DU was observed. Consistent results were obtained in sensitivity analyses. Furthermore, the GRS approach showed similar results in the several MR methods.
Conclusion: These findings do not support a causal role of AgeSmk, CigDay, SmkCes, alcohol consumption, and type 2 diabetes in the development of GU and DU. However, it is confirmed that SmkInit and BMI have a causal part in the development of GU and DU.
1 Introduction
Symptoms of the common gastrointestinal condition known as peptic ulcer disease (PUD) include nausea, vomiting, abdominal pain, and even bleeding or perforation in the digestive tract. PUD often develops in the proximal duodenum and stomach, but it can also occur in the esophagus, distal duodenum, or jejunum (Lanas and Chan 2017). PUD is generally referred to in clinical medicine primarily as gastric ulcer (GU) or duodenal ulcer (DU), and therefore, GU and DU are collectively referred to as PUD in this study. According to estimates, the overall population’s lifetime prevalence of PUD is between 5 and 10 percent, with an annual incidence of between 0.1 and 0.3 percent (Kurata et al., 1992; Rosenstock and Jørgensen 1995). In the previous 30 years, PUD-related morbidity and mortality have dramatically decreased, whereas in the last 15 years, there has been a progressive rising trend (Xie et al., 2022).
The pathogenesis of PUD is multifactorial and may be related to genetic and environmental factors. In previous studies, it was proposed that smoking (Levenstein et al., 1997; Shigemi et al., 1999; Andersen et al., 2000; Rosenstock et al., 2003), alcohol consumption (Levenstein et al., 1997; Andersen et al., 2000; Rosenstock et al., 2003), type 2 diabetes (Tseng et al., 2012), and BMI (Holtmann et al., 1992; Medalie et al., 1992; Boylan et al., 2014) may be possible elements for the development of PUD. The majority of the data on potential risk factors for PUD are based on observational research, which can be subject to problems from residual confounding. We now evaluate the causal association between smoking, alcohol consumption, type 2 diabetes, BMI, and the likelihood of developing GU and DU using Mendelian randomization (MR).
The drawbacks of observational studies, such as reverse causality, residual confounding, and recall bias, have been overcome by MR analysis, in which genetic variants are used as a surrogate for lifestyle/environmental exposures. Consequently, it is a sophisticated tool for enhancing causal inference (Vandebergh et al., 2022). The random distribution of genetic variants at conception and the fixed status of genetic variants that cannot be changed by the onset or progression of the disease, respectively, allow for the reduction of residual confounding and the diminishing of reverse causality, which are two main benefits of this method (Yuan et al., 2021). Here, we did a two-sample MR study to examine the relationships between the risk of GU and DU and genetically proxied smoking, alcohol consumption, type 2 diabetes, and BMI.
2 Materials and methods
2.1 Genetic instruments for smoking, alcohol consumption, type 2 diabetes, and BMI
The most recent meta-analysis on tobacco and alcohol consumption based on over 30 GWASs identified 566 and 99 genetic variants associated with four smoking phenotypes and alcohol consumption phenotypes in over 1.2 million people of European ancestry. These genetic instruments for smoking phenotypes and alcohol consumption were obtained from this study (Liu et al., 2019). There have been several descriptions of smoking and alcohol consumption phenotypes in detail (Liu et al., 2019). In summary, the smoking status (SmkIni, Ever/Never), age at smoking initiation (AgeSmk), smoking cessation (SmkCes, Current/Former), and number of cigarettes smoked per day (CigDay, both for current and former smokers, quantitative measures were binned into five bins or collected with predefined bins as follows: 1 = 1–5, 2 = 6–15, 3 = 16–25, 4 = 26–35, and 5 = 36 + cigarettes per day) were included (Zhou et al., 2021). For all phenotypes in family studies and quasi-continuous phenotypes (age of smoking initiation and cigarettes per day) in unrelated individuals, the association data were produced using a linear mixed model (Zhou et al., 2021). In research of unrelated individuals, a logistic model was used to calculate the additive genetic effects for the binary phenotypes of the smoking status and smoking cessation. Meta-analysis was carried out using rareGWAMA (McKay et al., 2017). Supplementary Tables S1–S5 provide comprehensive data regarding instrumental variables (IVs) related to alcohol intake and smoking phenotypes.
Based on a meta-analysis of 32 GWASs including 74,124 type 2 diabetes cases and 824,006 control people of European ancestry, the IVs selection for type 2 diabetes was conducted (known as the DIAGRAM consortium) (Mahajan et al., 2018). Single-nucleotide polymorphisms (SNPs) were identified as IVs for type 2 diabetes (n = 403) when they reached the genome-wide statistical significance criteria (p < 5 * 10–8). According to earlier research, BMI-mediation effects were entirely for the impacts of type 2 diabetes-related genetic variations in the FTO, MC4R, TMEM18, SEC16B, and GNPDA2 genes (Mahajan et al., 2018). The FTO, MC4R, TMEM18, SEC16B, and GNPDA2 gene regions were thus removed, leaving 394 SNPs as type 2 diabetes instrumental factors. In the sensitivity analysis for type 2 diabetes, 289 SNPs that reached the genome-wide significance threshold were included (FTO, MC4R, TMEM18, SEC16B, and GNPDA2 variations were removed). In Supplementary Table S6, comprehensive details for type 2 diabetes IVs are provided.
The exposure variable data for the genetic variants associated with BMI was derived from a GWAS meta-analysis in the Genetic Investigation of Anthropometric Traits (https://portals.broadinstitute.org/collaboration/giant/index.php/GIANT_consortium_data_files, accessed on 28 April 2021) consortium (n = 681,275 individuals of European ancestry) (Yengo et al., 2018). Detailed information for IVs of BMI is presented in Supplementary Table S7.
2.2 Genetic associations of SNPs with gastric ulcer and duodenal ulcer risk
Summary statistics on GU and DU risks, including standard errors (SE) and odds ratio (OR) estimates for instrumental SNPs, were available from Gene ATLAS (http://geneatlas.roslin.ed.ac.uk/) in 452,264 participants from European ancestry. Participants with more derived European ancestry, as indicated by these studies, were only allowed to participate.
2.3 Instrumental variable selection
Based on the following criteria, the IVs for MR studies were chosen: 1) a 500-kb window with r2 measure of LD among instruments <0.001; 2) p value lower than the genome-wide significant level (5 * 10–8 for smoking, alcohol consumption, type 2 diabetes, and BMI) found in the relevant study; 3) minor allele frequency (MAF) > 0.01; 4) palindromic SNPs (A/T and G/C polymorphisms were considered for elimination (MAF is between 0.4 and 0.6), but there was no ambiguity in the effects).
2.4 Statistical analyses
2.4.1 Mendelian randomization analysis
In our two-sample MR study, we used the inverse-variance weighted (IVW) approach as the main method to assess the causal relationship of smoking, alcohol consumption, type 2 diabetes, and BMI with GU and DU. IVW calculates the exposure-outcome effect corresponding to each SNP using the Wald ratio method, then performs a weighted linear regression with a forced intercept of zero. It achieved higher estimate accuracy and test power when IVs satisfied the following three underlying assumptions (Burgess et al., 2013): 1) the SNPs are significantly associated with the exposure (the “relevance” assumption); 2) the SNPs are not correlated with potential confounders of the exposure-outcome relation (the “exchangeability” assumption); and 3) the SNPs are not directly associated with the outcome (the “exclusion” restriction) (Figure 1). If there was significant heterogeneity, a random-effects IVW model would be implemented, which was less prone to the bias of weaker SNP-exposure associations (Bowden et al., 2017). To avoid the interference of unknown and unmeasurable confounders, we also performed the MR-Egger regression (MR-Egger) (Bowden et al., 2015), weighted median estimator (Verbanck et al., 2018b), simple median estimator (Verbanck et al., 2018b), MR pleiotropy residual sum and outlier (MR-PRESSO) test (Verbanck et al., 2018a), and the robust adjusted profile score (RAPS) (Zhao et al., 2020a) to test the robustness of our results. The study design overview is shown in Figure 1.
FIGURE 1. Study design overview and assumptions of the Mendelian randomization framework. IVW, inverse-variance weighted; LD, linkage disequilibrium; SNPs, single-nucleotide polymorphisms. MR-PRESSO, Mendelian randomization pleiotropy residual sum and outlier; MR-RAPS, Mendelian randomization robust-adjusted profile score. It is shown that (A) the genetic variants proposed as instrumental variables should be robustly associated with the risk factor of interest; (B) the used genetic variants should not be associated with potential confounders; and (C) the selected genetic variants should affect the risk of the outcome merely through the risk factor, not via alternative pathways.
The Bonferroni method was used to correct for multiple testing in our two-sample MR models. The association with two-sided p values <0.007 (where α = 0.05/7) were deemed statistically significant. Associations with p values between 0.05 and 0.007 were regarded as suggestive associations, requiring confirmation. We conducted all analyses using R (version 4.1.2) the R packages “TwoSampleMR” (Hemani et al., 2018) and “MendelianRandomization” (Yavorska and Burgess 2017). We used only freely accessible summarized data in this study; therefore, this work did not require ethical approval.
2.4.2 Associations between exposure genetic risk score (GRS) and outcome
MR analysis was carried out using weighted GRS as IVs by using the same summary data to produce the combined estimate of the connection between exposure-influencing alleles and the outcome. We used R (version 4.1.2) and the “gtx” R package (Windows version 0.0.8), whose grs. summary module contains the GRS function, to carry out the analyses. Similar to an approach where an outcome is regressed onto an additive GRS, the grs. summary module only used single SNP association summarized data from the outcomes of the GWAS research (Luo et al., 2020). This method is explained in depth elsewhere (Dastani et al., 2012). Additionally, prior research claimed that this MR approach using meta-GWAS summary data was just as effective as that using individual-level data (Voight et al., 2012).
2.4.3 Sensitivity analyses
2.4.3.1 Weighted median estimates
We provide weighted median estimates, which are accurate estimates when at least 50% of the data comes from reliable SNPs (Verbanck et al., 2018b).
2.4.3.2 MR-Egger
If the instrument strength independent of direct effect (InSIDE) condition is met, MR-Egger regression yields accurate estimates even when all the SNPs are incorrect (Bowden et al., 2015). The MR-Egger method relaxes the requirement that there is no pleiotropy between genetic variants in the IVW method. It assumes that the IVs–exposure and IVs–outcome associations are independent. This is referred to as the InSIDE assumption and is relatively weak compared to the strict exclusion restriction criteria. MR-Egger intercepts with p < 0.05, indicating that horizontal pleiotropy are present. Horizontal pleiotropy refers to the effect of genetic variation on a trait through multiple different paths. When using pleiotropy genetic variation as IVs, it is equivalent to constructing other pathways other than the “IVs–exposure–outcome” path, thus making the IVs fail due to the breach of core assumptions. The application of pleomorphic IVs can lead to bias in the estimation of the effect of the studied causal pathway. We concentrated more on the consistency of the estimate direction between MR-Egger and IVW because MR-Egger has less statistical power than IVW (Yeung et al., 2019).
2.4.3.3 MR pleiotropy residual sum and outlier (MR-PRESSO)
By examining outliers among the included SNPs that contribute to the MR estimate, MR-PRESSO can identify horizontal pleiotropy. This approach makes the assumptions that the InSIDE assumption is true, balanced pleiotropy, and that at least 50% of the SNPs are valid SNPs. Additionally, MR-PRESSO offers adjusted estimates and identifies the outlier SNPs. This approach was used by us as a statistical tool to detect and eliminate potential pleiotropic SNPs (Bowden et al., 2016). In summary, the MR-PRESSO method has several core functions: 1) to detect the presence of horizontal pleiotropy using the “MR-PRESSO global test,” MR-PRESSO global test p < 0.05 indicates the presence of horizontal pleiotropy; 2) to remove abnormal SNPs (outliers) and estimate the corrected results (which remove horizontal pleiotropy) using the “MR-PRESSO outlier test.” The MR-PRESSO (raw) method is the result before the removal of the abnormal SNPs, and MR-PRESSO (outlier-corrected) is the result after the removal of the abnormal SNPs. When the results of the MR-PRESSO (outlier-corrected) method show “NA,” this indicates that there are no outlier SNPs in the IVs we selected.
2.4.3.4 MR-robust adjusted profile score (MR-RAPS)
The MR-RAPS method corrected for horizontal pleiotropy in the IVW analysis by using robust-adjusted profile scores. In addition, MR-RAPS with Huber loss function which can model random-effects distribution of the pleiotropic effects of genetic variants is discussed (Zhao et al., 2020b).
2.4.3.5 Cochran’s Q test
To assess the extent of heterogeneity among the individual impact estimates produced from each genetic mutation, the Cochran’s Q test was calculated (Haycock et al., 2016). Cochran’s Q statistic was used to determine heterogeneity statistics. Heterogeneity was defined as a p value <0.05 from a Cochran’s Q calculation. A random-effects IVW model would be used if there was high heterogeneity since it was less susceptible to the bias of weaker SNPs-exposure relationships (Bowden et al., 2017).
2.4.3.6 F statistics
In order to detect the strength of the IVs at a threshold of F > 10, which is normally advocated in MR analysis, we calculated the F statistics of the identified all SNPs. The formula for the F statistics is
3 Results
3.1 Instrumental variables and their validity
A total of 231 independent SNPs for type 2 diabetes, 82 independent SNPs for alcohol consumption, 616 independent SNPs for BMI, and 8, 41, 17, and 325 SNPs for the AgeSmk, CigDay, SmkCes, and SmkInit of smoking phenotypes, respectively, were included as IVs (Supplementary Tables S1–S7).
3.2 Mendelian randomization
3.2.1 MR results of smoking to GU and DU
Genetic predisposition to SmkInit was associated with an increased risk of GU and DU. These associations are consistent across multiple MR methods. The beta and corresponding 95% CIs of GU were 0.0035 (95% CI, 0.0021, 0.0049) for one SD increase for SmkInit and DU were 0.0029 (95% CI, 0.0017, 0.0041) for one SD increase for SmkInit (Table 1, Supplementary Figures S1–S2). Genetic predisposition to AgeSmk, CigDay, and SmkCes was not associated with GU and DU susceptibilities (Table 1, Supplementary Figures S1–S2).
TABLE 1. Summarized results of the Mendelian randomization study on smoking, alcohol consumption, diabetes, and BMI and GU and DU risks.
3.2.2 MR results of alcohol consumption to GU and DU
Increased alcohol consumption had no discernible influence on the risk of GU and DU (beta = -0.0030, 95% CI = -0.0080, 0.0021, p = 0.25; beta = 0.0027, 95% CI = -0.0016, 0.0070, p = 0.21, respectively) (Table 1). The MR-Egger regression and Cochran’s Q test did not reveal any evidence that supported horizontal pleiotropy and heterogeneity. MR-PRESSO shows an SNP with horizontal pleiotropic, with one SNP identified as outliers (rs281379). When we removed rs281379, our results were stable, and neither horizontal pleiotropic nor heterogeneity was detected (Table 3). The results of all methods were consistent, and MR-PRESSO did not show the existence of pleiotropic SNPs (Table 1, Supplementary Figures S1–S2).
3.2.3 MR results of type 2 diabetes to GU and DU
Genetic predisposition to type 2 diabetes was not associated with an increased risk of GU and DU. These associations are consistent across multiple MR methods. The beta and corresponding 95% CIs of GU were 0.0002 (95% CI, -0.0003, 0.0007) for type 2 diabetes and DU were 0.0002 (95% CI, -0.0002, 0.0007) for type 2 diabetes (Table 1, Supplementary Figures S1–S2). The MR-Egger regression intercept provided no support for directional pleiotropy (Table 3). The Cochran’s Q test did show some heterogeneity among individual SNPs effect estimates in type 2 diabetes on DU; therefore, a random-effects IVW model would be implemented. Random-effects IVW model results were consistent with other methods (fixed-effects IVW model weighted median estimator, simple median estimator, MR-RAPS, and MR-PRESSO) (Table 1, Supplementary Figures S1–S2).
3.2.4 MR results of BMI to GU and DU
There was a link between a higher risk of GU and DU and a genetic tendency to BMI. The beta and corresponding 95% CIs of GU were 0.0021 (95% CI, 0.0009, 0.0033) for one SD increase for BMI and DU were 0.0018 (95% CI, 0.0007, 0.0029) for one SD increase for BMI (Table 1, Supplementary Figures S1–S2). The Cochran’s Q test and MR-Egger regression intercept for the relationship between BMI and GU did not reveal any evidence of directional pleiotropy or heterogeneity, and the Cochran’s Q test and MR-Egger regression intercept did show heterogeneity and directional pleiotropy among individual SNPs effect estimates in BMI on DU (Table 3). However, no SNPs were identified as an outlier in the MR-PRESSO method in the BMI to DU MR analysis. In addition, random-effects IVW model results were consistent with other methods (fixed-effects IVW model, weighted median estimator, simple median estimator, MR-RAPS, and MR-PRESSO) (Table 1, Supplementary Figures S1–S2).
3.3 GRSsmoking, GRSalcohol, GRStype 2 diabetes, and GRSBMI with GU and DU
Consistent with the MR results of smoking (AgeSmk, CigDay, SmkCes, and Smklnit), alcohol consumption, type 2 diabetes, and BMI to GU and DU, the GRS revealed no causal effect of smoking (AgeSmk, CigDay, and SmkCes), alcohol consumption, and type 2 diabetes on GU and DU risks (Table 2). The GRS revealed significant causal effect of smoking (Smklnit) and BMI on GU and DU (Table 2). The beta and corresponding 95% CIs of GU were 0.0035 (95% CI, 0.0022, 0.0048) for one SD increase for Smklnit and DU were 0.0029 (95% CI, 0.0018, 0.0041) for one SD increase for Smklnit (Table 2, Supplementary Figure S3). The beta and corresponding 95% CIs of GU were 0.0021 (95% CI, 0.0009, 0.0033) for one SD increase for BMI and DU were 0.0018 (95% CI, 0.0008, 0.0028) for one SD increase for BMI (Table 2, Supplementary Figure S3).
TABLE 2. Effect of the GRS instrument of smoking, alcohol consumption, diabetes, and BMI on GU, DU.
3.4 Heterogeneity and sensitivity analysis
Cochran’s Q tests revealed that there was no heterogeneity between the smoking, alcohol consumption, type 2 diabetes, and BMI IVs in the study of smoking, alcohol consumption, type 2 diabetes, BMI, and GU MR (Table 3). In smoking and alcohol consumption to DU MR analysis, the results of Cochran’s Q tests showed that there was no heterogeneity among the smoking and alcohol consumption-related IVs (Table 3). In type 2 diabetes and BMI to DU MR analysis, the Cochran’s Q tests revealed some heterogeneity between both the IVs for type 2 diabetes (Q = 284.7, p = 0.01) and BMI (Q = 689.6, p = 0.02) (Table 3); however, random-effects IVW model results were consistent with other methods (fixed-effects IVW model, weighted median estimator, simple median estimator, MR-RAPS, and MR-PRESSO). Therefore, our results are still very reliable and this little heterogeneity has little effect on our results.
TABLE 3. Potential pleiotropy evaluation using different methods.
The MR-Egger regression analysis revealed that alcohol consumption, type 2 diabetes, smoking, and BMI to GU MR did not indicate horizontal pleiotropy of the IVs (Table 3), and the horizontal pleiotropy of the IVs was not present for alcohol consumption, smoking (CigDay, SmkCes, and Smklnit), and type 2 diabetes to DU MR and it was present for smoking (AgeSmk, MR-Egger p = 0.01) and BMI (MR-Egger p = 0.049) (Table 3). The MR-PRESSO global test p was >0.05 in smoking (AgeSmk) and global test p was <0.05 in BMI to DU MR, with no SNPs identified as outliers (Table 3).
3.5 Excluding associations between individual exposure factors
To exclude associations between individual exposure factors, we searched the phenoscanner website for SNPs of each exposure factor (Kamat et al., 2019). SNPs that were strongly correlated with other exposure factors were eliminated. Specific screening results of SNPs for each exposure factor are presented in Supplementary Table S8. After removing these SNPs, which were strongly correlated with other exposure factors, we performed MR analysis again, and we got the same results as before (Supplementary Table S9). Therefore, our results can be considered very reliable.
4 Discussion
Increasing evidence indicates the association between smoking, alcohol consumption, type 2 diabetes, and BMI with GU and DU risks, and several plausible mechanisms underlying are proposed. According to several studies, smoking constricts the mucosa’s blood vessels while ischemia lowers the mucosa’s resistance (Andersen et al., 2000). Given that smokers’ levels of carboxyhemoglobin are higher, carbon monoxide may have a further role in the mucosal ischemia (Nolan et al., 1985; Smith et al., 1998). The substantial link between smoking and PUD may be explained by these latter mechanisms. When exposed to high levels of ethanol quickly, the result has significant damage to the gastric epithelium, necrosis of deeper layers of the mucosa, and microvascular damage that causes engorgement, increased permeability, and intramucosal bleeding (Tarnawski et al., 1987). Apart from its irritant properties locally, ethanol has been demonstrated to slow down stomach emptying at moderate to high dosages (Jian et al., 1986). Type 2 diabetes-related angiopathy could compromise mucosal integrity and result in more serious ulcers (Weil et al., 2000). Past researches looked into the relationship between PUD and BMI, but the results were still debatable. Past researches have reported obesity as a separate risk factor for PUD (Aro et al., 2006; Garrow and Delegge 2010; Boylan et al., 2014). Due to the observational study’s design, these researches are still unable to control the impact of potential biases. This has led to growing interest in proving a causal relationship between smoking, alcohol consumption, type 2 diabetes, and BMI with GU and DU risks. Our research is the first largest MR analysis the causal association between smoking, alcohol consumption, type 2 diabetes, and BMI with GU and DU to the best of our knowledge.
According to the results of the current MR study, genetic predisposition to BMI and smoking (Smklnit) both enhance the likelihood of developing GU and DU. However, there is no evidence linking a genetic propensity to smoking (AgeSmk, CigDay, and SmkCes), alcohol consumption, or type 2 diabetes to an increased risk of GU and DU.
The following are this study’s strengths and limitations. First, the MR method, which can reduce reverse causality and diminish residual confounding, is the study’s main strength. Second, to assure an accurate outcome, we chose the largest GWAS databases on smoking, alcohol consumption, type 2 diabetes, and BMI. Third, for the first time, we used MR to draw the conclusion that genetically predicted BMI and smoking (Smklnit) increase the risk of GU and DU, which also offers some reasons for certain earlier epidemiological studies. At last, secondary MR analyses were conducted within summary data using GRSs as IVs, and we got the same result as the previous analysis. However, we also recognize that there are some drawbacks. First, since only people of European ancestry were included in the GWASs from which we chose our IVs, it is difficult to extrapolate our results to other ethnicities. Second, potential pleiotropy is a drawback of our work. Third, we did not conduct further multivariable MR analysis to make our research results more convincing. However, in order to exclude associations between exposure individual factors, we searched the phenoscanner website for the SNPs of each exposure factor, excluded SNPs that were strongly correlated with other exposure factors and product the MR analysis again, which was found to be consistent with the previous results, thus, to some extent, compensating for the lack of multivariable MR analysis.
Summary of findings: In a large, well-powered study, we did not find strong evidence for a possible association for genetically predicted smoking (AgeSmk, CigDay, and SmkCes), alcohol consumption, or type 2 diabetes in the risk of GU and DU. This suggests that previous study may be confounded by potential biases or due to reverse causation. However, we found a causal relationship between smoking (Smklnit) and BMI on the one hand and GU and DU on the other. Nevertheless, our results require further study of large sample data in the future.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
Author contributions
The study’s conception and design were influenced by YL and ZX. Statistics were analyzed and the text was written by YL and KY. The tables and figures were made by ZX and LX, and YZ provided writing guidance, revised the book draft, and oversaw the entire endeavor. The final draft of the work was approved by all authors.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.992080/full#supplementary-material
References
Andersen, I. B., Jørgensen, T., Bonnevie, O., Grønbaek, M., and Sørensen, T. I. (2000). Smoking and alcohol intake as risk factors for bleeding and perforated peptic ulcers: A population-based cohort study. Epidemiology 11, 434–439. doi:10.1097/00001648-200007000-00012
Aro, P., Storskrubb, T., Ronkainen, J., Bolling-Sternevald, E., Engstrand, L., Vieth, M., et al. (2006). Peptic ulcer disease in a general adult population: The kalixanda study: A random population-based study. Am. J. Epidemiol. 163, 1025–1034. doi:10.1093/aje/kwj129
Bowden, J., Davey Smith, G., and Burgess, S. (2015). Mendelian randomization with invalid instruments: Effect estimation and bias detection through egger regression. Int. J. Epidemiol. 44, 512–525. doi:10.1093/ije/dyv080
Bowden, J., Davey Smith, G., Haycock, P. C., and Burgess, S. (2016). Consistent estimation in mendelian randomization with some invalid instruments using a weighted median estimator. Genet. Epidemiol. 40, 304–314. doi:10.1002/gepi.21965
Bowden, J., Del Greco, M. F., Minelli, C., Davey Smith, G., Sheehan, N., and Thompson, J. (2017). A framework for the investigation of pleiotropy in two-sample summary data Mendelian randomization. Stat. Med. 36, 1783–1802. doi:10.1002/sim.7221
Boylan, M. R., Khalili, H., Huang, E. S., and Chan, A. T. (2014). Measures of adiposity are associated with increased risk of peptic ulcer. Clin. Gastroenterol. Hepatol. 12, 1688–1694. doi:10.1016/j.cgh.2014.03.021
Burgess, S., Butterworth, A., and Thompson, S. G. (2013). Mendelian randomization analysis with multiple genetic variants using summarized data. Genet. Epidemiol. 37, 658–665. doi:10.1002/gepi.21758
Dastani, Z., Hivert, M. F., Timpson, N., Perry, J. R., Yuan, X., Scott, R. A., et al. (2012). Novel loci for adiponectin levels and their influence on type 2 diabetes and metabolic traits: A multi-ethnic meta-analysis of 45, 891 individuals. PLoS Genet. 8, e1002607. doi:10.1371/journal.pgen.1002607
Garrow, D., and Delegge, M. H. (2010). Risk factors for gastrointestinal ulcer disease in the US population. Dig. Dis. Sci. 55, 66–72. doi:10.1007/s10620-008-0708-x
Haycock, P. C., Burgess, S., Wade, K. H., Bowden, J., Relton, C., and Davey Smith, G. (2016). Best (but oft-forgotten) practices: The design, analysis, and interpretation of mendelian randomization studies. Am. J. Clin. Nutr. 103, 965–978. doi:10.3945/ajcn.115.118216
Hemani, G., Zheng, J., Elsworth, B., Wade, K. H., Haberland, V., Baird, D., et al. (2018). The MR-Base platform supports systematic causal inference across the human phenome. Elife 7, e34408. doi:10.7554/eLife.34408
Holtmann, G., Armstrong, D., Pöppel, E., Bauerfeind, A., Goebell, H., Arnold, R., et al. (1992). Influence of stress on the healing and relapse of duodenal ulcers. A prospective, multicenter trial of 2109 patients with recurrent duodenal ulceration treated with ranitidine. RUDER Study Group. Scand. J. Gastroenterol. 27 (11), 917–923. doi:10.3109/00365529209000163
Jian, R., Cortot, A., Ducrot, F., Jobin, G., Chayvialle, J. A., Modigliani, R., et al. (1986). Effect of ethanol ingestion on postprandial gastric emptying and secretion, biliopancreatic secretions, and duodenal absorption in man. Dig. Dis. Sci. 31, 604–614. doi:10.1007/BF01318691
Kamat, M. A., Blackshaw, J. A., Young, R., Surendran, P., Burgess, S., Danesh, J., et al. (2019). PhenoScanner V2: An expanded tool for searching human genotype-phenotype associations. Bioinformatics 35, 4851–4853. doi:10.1093/bioinformatics/btz469
Kurata, J. H., Nogawa, A. N., Abbey, D. E., and Petersen, F. (1992). A prospective study of risk for peptic ulcer disease in Seventh-Day Adventists. Gastroenterology 102, 902–909. doi:10.1016/0016-5085(92)90176-y
Lanas, A., and Chan, F. K. L. (2017). Peptic ulcer disease. Lancet 390, 613–624. doi:10.1016/S0140-6736(16)32404-7
Levenstein, S., Kaplan, G. A., and Smith, M. W. (1997). Psychological predictors of peptic ulcer incidence in the Alameda County Study. J. Clin. Gastroenterol. 24, 140–146. doi:10.1097/00004836-199704000-00004
Liu, M., Jiang, Y., Wedow, R., Li, Y., Brazel, D. M., Chen, F., et al. (2019). Association studies of up to 1.2 million individuals yield new insights into the genetic etiology of tobacco and alcohol use. Nat. Genet. 51, 237–244. doi:10.1038/s41588-018-0307-5
Luo, Q., Wen, Z., Li, Y., Chen, Z., Long, X., Bai, Y., et al. (2020). Assessment causality in associations between serum uric acid and risk of schizophrenia: A two-sample bidirectional mendelian randomization study. Clin. Epidemiol. 12, 223–233. doi:10.2147/CLEP.S236885
Mahajan, A., Taliun, D., Thurner, M., Robertson, N. R., Torres, J. M., Rayner, N. W., et al. (2018). Fine-mapping type 2 diabetes loci to single-variant resolution using high-density imputation and islet-specific epigenome maps. Nat. Genet. 50, 1505–1513. doi:10.1038/s41588-018-0241-6
McKay, J. D., Hung, R. J., Han, Y., Zong, X., Carreras-Torres, R., Christiani, D. C., et al. (2017). Large-scale association analysis identifies new lung cancer susceptibility loci and heterogeneity in genetic susceptibility across histological subtypes. Nat. Genet. 49, 1126–1132. doi:10.1038/ng.3892
Medalie, J. H., Stange, K. C., Zyzanski, S. J., and Goldbourt, U. (1992). The importance of biopsychosocial factors in the development of duodenal ulcer in a cohort of middle-aged men. Am. J. Epidemiol. 136, 1280–1287. doi:10.1093/oxfordjournals.aje.a116436
Nolan, J., Jenkins, R. A., Kurihara, K., and Schultz, R. C. (1985). The acute effects of cigarette smoke exposure on experimental skin flaps. Plast. Reconstr. Surg. 75, 544–551. doi:10.1097/00006534-198504000-00018
Pierce, B. L., and Burgess, S. (2013). Efficient design for mendelian randomization studies: Subsample and 2-sample instrumental variable estimators. Am. J. Epidemiol. 178, 1177–1184. doi:10.1093/aje/kwt084
Rosenstock, S. J., and Jørgensen, T. (1995). Prevalence and incidence of peptic ulcer disease in a Danish County--a prospective cohort study. Gut 36, 819–824. doi:10.1136/gut.36.6.819
Rosenstock, S., Jørgensen, T., Bonnevie, O., and Andersen, L. (2003). Risk factors for peptic ulcer disease: A population based prospective cohort study comprising 2416 Danish adults. Gut 52, 186–193. doi:10.1136/gut.52.2.186
Shigemi, J., Mino, Y., and Tsuda, T. (1999). The role of perceived job stress in the relationship between smoking and the development of peptic ulcers. J. Epidemiol. 9, 320–326. doi:10.2188/jea.9.320
Smith, C. J., Guy, T. D., Stiles, M. F., Morton, M. J., Collie, B. B., Ingebrethsen, B. J., et al. (1998). A repeatable method for determination of carboxyhemoglobin levels in smokers. Hum. Exp. Toxicol. 17, 29–34. doi:10.1177/096032719801700105
Tarnawski, A., Hollander, D., Stachura, J., Klimczyk, B., Mach, T., and Bogdal, J. (1987). Alcohol injury to the normal human gastric mucosa: Endoscopic, histologic and functional assessment. Clin. Invest. Med. 10, 259–263.
Tseng, P. H., Lee, Y. C., Chiu, H. M., Chen, C. C., Liao, W. C., Tu, C. H., et al. (2012). Association of diabetes and HbA1c levels with gastrointestinal manifestations. Diabetes Care 35, 1053–1060. doi:10.2337/dc11-1596
Vandebergh, M., Becelaere, S., Dubois, B., and Goris, A. (2022). Body mass index, interleukin-6 signaling and multiple sclerosis: A mendelian randomization study. Front. Immunol. 13, 834644. doi:10.3389/fimmu.2022.834644
Verbanck, M., Chen, C. Y., Neale, B., and Do, R. (2018a). Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat. Genet. 50, 693–698. doi:10.1038/s41588-018-0099-7
Verbanck, M., Chen, C. Y., Neale, B., and Do, R. (2018b). Publisher Correction: Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat. Genet. 50, 1196. doi:10.1038/s41588-018-0164-2
Voight, B. F., Peloso, G. M., Orho-Melander, M., Frikke-Schmidt, R., Barbalic, M., Jensen, M. K., et al. (2012). Plasma HDL cholesterol and risk of myocardial infarction: A mendelian randomisation study. Lancet 380, 572–580. doi:10.1016/S0140-6736(12)60312-2
Weil, J., Langman, M. J., Wainwright, P., Lawson, D. H., Rawlins, M., Logan, R. F., et al. (2000). Peptic ulcer bleeding: Accessory risk factors and interactions with non-steroidal anti-inflammatory drugs. Gut 46, 27–31. doi:10.1136/gut.46.1.27
Xie, X., Ren, K., Zhou, Z., Dang, C., and Zhang, H. (2022). The global, regional and national burden of peptic ulcer disease from 1990 to 2019: A population-based study. BMC Gastroenterol. 22, 58. doi:10.1186/s12876-022-02130-2
Yavorska, O. O., and Burgess, S. (2017). MendelianRandomization: an R package for performing Mendelian randomization analyses using summarized data. Int. J. Epidemiol. 46, 1734–1739. doi:10.1093/ije/dyx034
Yengo, L., Sidorenko, J., Kemper, K. E., Zheng, Z., Wood, A. R., Weedon, M. N., et al. (2018). Meta-analysis of genome-wide association studies for height and body mass index in ∼700000 individuals of European ancestry. Hum. Mol. Genet. 27, 3641–3649. doi:10.1093/hmg/ddy271
Yeung, C. H. C., Au Yeung, S. L., Fong, S. S. M., and Schooling, C. M. (2019). Lean mass, grip strength and risk of type 2 diabetes: A bi-directional mendelian randomisation study. Diabetologia 62, 789–799. doi:10.1007/s00125-019-4826-0
Yuan, S., Giovannucci, E. L., and Larsson, S. C. (2021). Gallstone disease, diabetes, calcium, triglycerides, smoking and alcohol consumption and pancreatitis risk: Mendelian randomization study. NPJ Genom. Med. 6, 27. doi:10.1038/s41525-021-00189-6
Zhao, Q., Wang, J., Hemani, G., Bowden, J., and Small, D. S. (2020a). Statistical inference in two-sample summary-data Mendelian randomization using robust adjusted profile score. Ann. Stat. 48. doi:10.1214/19-aos1866
Zhao, Q., Wang, J., Hemani, G., Bowden, J., and Small, D. S. (2020b). Statistical inference in two-sample summary-data Mendelian randomization using robust adjusted profile score. Inst. Math. Statistics 48. doi:10.1214/19-AOS1866
Keywords: smoking, alcohol consumption, diabetes, body mass index, peptic ulcer, Mendelian randomization
Citation: Liu Y, Xiao Z, Ye K, Xu L and Zhang Y (2023) Smoking, alcohol consumption, diabetes, body mass index, and peptic ulcer risk: A two-sample Mendelian randomization study. Front. Genet. 13:992080. doi: 10.3389/fgene.2022.992080
Received: 12 July 2022; Accepted: 28 November 2022;
Published: 06 January 2023.
Edited by:
Pingzhao Hu, Western University, CanadaReviewed by:
Osvaldo Espin-Garcia, Western University, CanadaQian Liu, University of Manitoba, Canada
Copyright © 2023 Liu, Xiao, Ye, Xu and Zhang. 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: Yanping Zhang, zyping001@126.com
†These authors have contributed equally to this work and share first authorship