ORIGINAL RESEARCH article

Front. Nutr., 01 May 2024

Sec. Nutritional Epidemiology

Volume 11 - 2024 | https://doi.org/10.3389/fnut.2024.1388732

Causal associations between dietary factors and colorectal cancer risk: a Mendelian randomization study

  • 1. Department of Gastroenterology, The First Affiliated Hospital of Xi’an Jiaotong University, Xi’an, China

  • 2. Shannxi Clinical Research Center of Digestive Disease (Cancer Division), Xi’an, China

  • 3. Health Science Center, Xi’an Jiaotong University, Xi’an, China

  • 4. Shanxi Institute of Science and Technology, Jincheng, China

  • 5. Department of Surgery, Dangtu Central Health Center, Ma’anshan, China

Article metrics

View details

7

Citations

4,2k

Views

1,4k

Downloads

Abstract

Background:

Previous epidemiological studies have found a link between colorectal cancer (CRC) and human dietary habits. However, the inherent limitations and inevitable confounding factors of the observational studies may lead to the inaccurate and doubtful results. The causality of dietary factors to CRC remains elusive.

Methods:

We conducted two-sample Mendelian randomization (MR) analyses utilizing the data sets from the IEU Open GWAS project. The exposure datasets included alcoholic drinks per week, processed meat intake, beef intake, poultry intake, oily fish intake, non-oily fish intake, lamb/mutton intake, pork intake, cheese intake, bread intake, tea intake, coffee intake, cooked vegetable intake, cereal intake, fresh fruit intake, salad/raw vegetable intake, and dried fruit intake. In our MR analyses, the inverse variance weighted (IVW) method was employed as the primary analytical approach. The weighted median, MR-Egger, weighted mode, and simple mode were also applied to quality control. Heterogeneity and pleiotropic analyses were implemented to replenish the accuracy of the results.

Results:

MR consequences revealed that alcoholic drinks per week [odds ratio (OR): 1.565, 95% confidence interval (CI): 1.068–2.293, p = 0.022], non-oily fish intake (OR: 0.286; 95% CI: 0.095–0.860; p = 0.026), fresh fruit intake (OR: 0.513; 95% CI: 0.273–0.964; p = 0.038), cereal intake (OR: 0.435; 95% CI: 0.253–0.476; p = 0.003) and dried fruit intake (OR: 0.522; 95% CI: 0.311–0.875; p = 0.014) was causally correlated with the risk of CRC. No other significant relationships were obtained. The sensitivity analyses proposed the absence of heterogeneity or pleiotropy, demonstrating the reliability of the MR results.

Conclusion:

This study indicated that alcoholic drinks were associated with an increased risk of CRC, while non-oily fish intake, fresh fruit intake, cereal intake, and dried fruit were associated with a decreased risk of CRC. This study also indicated that other dietary factors included in this research were not associated with CRC. The current study is the first to establish the link between comprehensive diet-related factors and CRC at the genetic level, offering novel clues for interpreting the genetic etiology of CRC and replenishing new perspectives for the clinical practice of gastrointestinal disease prevention.

1 Introduction

Colorectal carcinoma (CRC) is the third most commonly diagnosed cancer worldwide, accounting for 9.4% of cancer-related fatalities globally (1). CRC patients exhibit clinical manifestations, including bowel habits changes, occult or overt rectal bleeding, abdominal pain, and anemia. However, in the early phase, patients are primarily asymptomatic or exhibit minor symptoms like common bowel diseases. When their bodies present a series of perceptible abnormalities, the cancer has already progressed to an advanced stage, even metastasized. Localized CRC patients have a high 5-year survival rate, decreasing from approximately 90% for primary tumors to 14% for metastatic CRC (2). With the incidence increasing constantly worldwide (1), CRC poses a significant challenge to public health globally. Individuals affected by CRC, including the patients and their families, fall into physical as well as financial adversities that ensue (3). Furthermore, CRC patients face psychological distress, including anxiety and depression (4). Eventually, the prolonged physical and mental issues may worsen the quality of life of patients. In addition, this disease not only presents a severe threat to personal health but also consumes substantial social and medical resources and heavily burdens society and healthcare systems (5). Clarifying the pathogenesis and etiology, including potential risk and protective factors, has excellent significance for the clinical practice of disease prevention and management.

Although the cause of CRC is still unclear, several researches have revealed some risk factors functionally integrated in the progression of this gastrointestinal disease. The Global Burden of Diseases, Injuries, and Risk Factors Study (GBD) indicated that the incidence rates of CRC increased with age, particularly surging in individuals aged 50–54 years and older (6). Additionally, a genome-wide association study identified 155 high-confidence effector genes that were functionally related to CRC risk, such as ARHGEF4, GNA12, LRIG1, GAB1, CNIH2, etc. These genes have multiple functions and affect tumor biology through various biological processes, including proliferation, homeostasis, migration, cell adhesion, immunity, and microbial interactions (7). Previous studies also found that environmental risk factors, Sedentary behavior (RR: 1.30, 95% CI: 1.22–1.39) (8) and smoking (RR: 1.17; 95% CI: 1.15–1.20) (9), could potentially impact the risk of CRC. Notably, in the realm of diet and nutrition, many experimental and epidemiological studies have made significant findings. For instance, Diets low in milk or calcium have been identified as primary contributors to the CRC disability-adjusted life years (6). Moreover, it has been found that nutritional supplements, such as omega-3 and arginine supplementation, could also modify the risk of CRC development (10).

According to previous studies, alcohol intake (11), red meat intake, processed meat intake (12), vegetable intake, and fruit intake (13) were associated with the pathophysiology of CRC. The potential mechanism of these pathologies is complicated and may contain direct biological effects on epithelial cells, modifications in inflammation and immune reactions, and diet-induced regulation in the composition and abundance of human gut microbes (14).

Current observational and meta-analysis studies on dietary factors and CRC face inherent limitations. Firstly, the sample sizes are typically small, affecting the reliability of results. Additionally, potential confounders may interfere with the interpretation of findings. Due to these factors, it’s challenging for these studies to conclusively demonstrate the epidemiological link between dietary habits and CRC risk. Hence, more robust and high-quality evidence is necessary to bridge the existing research gap.

Since the relationship between dietary factors and CRC has not been explored by any genetic instruments, we hypothesized there was a causative association of CRC with dietary factors. Similar to randomized controlled trials, the Mendelian randomization (MR) study is a novel research method that uses single-nucleotide polymorphisms (SNPs) as instrumental variables (IVs) to infer causal relationships between risk factors and health outcomes (15). This research methodology draws upon Mendel’s second law of genetics. It involves categorizing the study cohort according to the presence of specific genetic variations and subsequently comparing the occurrence of outcomes across these categories (16). SNPs follow the principle of being randomly allocated during the process of meiosis. This helps to eliminate the influence of confounding factors and the possibility of reverse causation, as genetic variants exist before the onset of the disease (17, 18). Recent MR studies suggest that dietary habits have a significant effect on several cardiovascular diseases (16) and five major mental disorders (19). Through MR studies, more diet-related factors for various diseases could be investigated. Herein, we performed a two-sample MR design to investigate the possible association of CRC with dietary factors.

2 Materials and methods

2.1 Study design

A flowchart (Figure 1) presents our study design concisely, including the procedure of selecting IVs, conducting MR studies using five methods, and carrying out sensitivity analyses. To provide a better understanding of our study design, it’s important to detail the foundation of MR, which consists of three essential assumptions. The first assumption is that the SNPs employed as IVs are supposed to be closely related to exposure factors. The second assumption indicates that the screened IVs should not be associated with any confounding factors. The third assumption requires that the proposed genetic variants should impact the risk of the health outcome only via exposure we focused on Chen et al. (15). The three crucial assumptions guaranteed that the MR results would not be interfered with by other confounding factors, such as the population’s characteristics, environment, and socioeconomic status. Also, since the genetic variation explains the formation of the exposure part before the outcome, reverse causality can be eliminated, thus compensating for the limitations of traditional methods. The two-sample MR analysis was performed to identify the causal relationship between traits utilizing publicly available genetic datasets in several genome-wide association studies (GWAS).

Figure 1

2.2 Data sources

Dietary factors employed in our study covered drinks intake (alcoholic drinks per week, tea intake, and coffee intake), vegetable and fruit intake (salad/raw vegetable intake, cooked vegetable intake, fresh fruit intake and dried fruit intake), meat intake (pork intake, beef intake, lamb/mutton intake, poultry intake, oily fish intake, non-oily fish intake, and processed meat intake), staple diet intake (bread intake and cereal intake), and dairy product intake (cheese intake). These GWAS summary-level data were obtained from the UK Biobank by the IEU open GWAS project, supported by the MRC Integrative Epidemiology Unit (IEU) at the University of Bristol. The GWAS summary-level data of CRC was extracted from the European Bioinformatics Institute by the IEU open GWAS project. More relevant information about the original datasets is shown in Table 1 and Supplementary Table S1. All the data used in this work are publicly available and were obtained from studies with the consent and ethical approval of the relevant participants. As a result, this study did not require the ethical approval of an institutional review board.

Table 1

GWAS IDExposure/outcomeIdentified SNPsParticipants included in the analysisF-statistic
ieu-b-73Alcoholic drinks per week32335,394 European-descent individuals80.012
ukb-b-6324Processed meat intake23461,981 European-descent individuals38.536
ukb-b-8006Poultry intake7461,900 European-descent individuals32.539
ukb-b-2862Beef intake15461,053 European-descent individuals40.509
ukb-b-17627Non-oily fish intake11460,880 European-descent individuals44.802
ukb-b-2209Oily fish intake59460,443 European-descent individuals45.100
ukb-b-5460Pork intake14460,162 European-descent individuals37.686
ukb-b-14179Lamb/mutton intake31460,006 European-descent individuals39.797
ukb-b-11348Bread intake29452,236 European-descent individuals41.884
ukb-b-1489Cheese intake62451,486 European-descent individuals39.238
ukb-b-8089Cooked vegetable intake17448,651 European-descent individuals37.584
ukb-b-6066Tea intake40447,485 European-descent individuals61.576
ukb-b-3881Fresh fruit intake51446,462 European-descent individuals46.191
ukb-b-15926Cereal intake39441,640 European-descent individuals46.028
ukb-b-1996Salad/raw vegetable intake18435,435 European-descent individuals39.230
ukb-b-5237Coffee intake38428,860 European-descent individuals74.525
ukb-b-16576Dried fruit intake41421,764 European-descent individuals42.200
ebi-a-GCST012876Colorectal cancerNA11,895 European-descent cases and 14,695 European-descent controlsNA

Information of the exposures and outcome datasets.

GWAS, Genome-Wide Association Studies; SNPs, single-nucleotide polymorphisms; N/A, not applicable.

2.3 Genetic variants

In order to meet the three assumptions of the MR analysis, the quality control steps below were applied to screen the related SNPs. We selected SNPs that are closely associated with various dietary factors. This selection was based on a genome-wide significant level (p < 5 × 10−8). We also performed the clumping process [distance window of 10,000 kb, linkage disequilibrium (LD) coefficient r2 < 0.001] (20). This step was crucial to avoid LD between SNPs and to ensure the independence of genetic variants. We selected the SNPs closely associated with various dietary factors at the significant level of genome-wide (p < 5 × 10−8) and conducted the clumping process [distance window 10,000 kb, linkage disequilibrium (LD) coefficient r2 < 0.001] to avoid LD between SNPs and ensure the independence of genetic variants (20). If no SNP was intensely related to any dietary factors found in the CRC database, proxy SNPs were allowed with a minimum LD R2 = 0.8 (21). Palindrome SNPs were reserved based on the threshold that the minor allele frequency (MAF) <0.3 (22). Notably, if the allele frequency contained in the details of an SNP is close to 0.5, we could hardly pinpoint the minor allele exactly, as there is sampling variance around the allele frequency. To enhance the accuracy of our study, we excluded such SNPs at the outset of MR analyses. In addition, to measure the power of the screened IVs and ensure their close relationships with exposures, we calculated the F-statistics and the proportion of variance interpreted (R2) for each SNP. Genetic variants (F-statistics <10) were generally considered as weak instruments, which should be removed from our MR analysis (23). Finally, MR-PRESSO tests were also employed to recognize potential horizontal pleiotropy, and the identified outliers would be ruled out to prevent the impact of pleiotropy (24).

2.4 Statistical analysis

We first performed an inverse variance weighted (IVW) test. This test is recognized for its strongest ability to determine causation (25). We applied it as the primary method to identify the causal effect between diet-related factors and CRC. We performed the inverse variance weighted (IVW) test, which possessed the most substantial ability to determine causation (25), as the significant method to detect the causal effect between diet-related factors and CRC. The evidence from the IVW method was complemented with the MR-Egger, weighted mode, weighted median, and simple mode. The conclusion would be more credible, stable, and precise when the consequences of these methods were consistent (26). For the IVW test and MR-Egger model, Cochran’s Q test was conducted to assess heterogeneity (27). Cochran’s Q test p < 0.05 indicated the existence of heterogeneity. Besides the MR-PRESSO test, as stated earlier, we also used the MR-Egger intercept test to detect directional pleiotropy. The absence of non-zero intercepts (p > 0.05) indicated that IVs did not affect CRC through other confounders (28). Leave-one-out analysis was applied to judge whether the causal link was affected by eliminating a particular SNP (29).

Statistical analysis was carried out with R software using the “TwoSampleMR” (20) package and “MR-PRESSO” (24). The significant threshold of the existence of causation is p < 0.05.

3 Results

3.1 Selection of instrumental variables

The causal associations of dietary factors with CRC were analyzed with 17 different exposures. The number of SNPs employed in our study ranged from 7 to 62. The F-statistics were greater than 10 for all the IVs (range: 32.539 to 80.012), suggesting that the selected IVs fulfilled the requirements of intense association with exposure. The amounts of European participants included in the exposure datasets ranged from 335,394 to 461,981. The outcome dataset covered 11,895 European-descent CRC cases and 14,695 European-descent controls. It was sourced from the European Bioinformatics Institute. Compared with the exposure datasets, there was little potential deviation in population stratification. More detailed information is presented in Table 1. Due to the non-significant conclusions of the MR-PRESSO global test (p > 0.05), no outlier was eliminated through MR-PRESSO.

3.2 MR analysis of dietary factors for CRC

In our study, a total of 5 causal relationships were discovered (p < 0.05 by IVW method). We identified that alcoholic drinks per week (OR: 1.565; 95% CI: 1.068–2.293; p = 0.022) was relevant to a higher risk of CRC. Non-oily fish intake (OR: 0.286; 95% CI: 0.095–0.860; p = 0.026), fresh fruit intake (OR: 0.513; 95% CI: 0.273–0.964; p = 0.038), cereal intake (OR: 0.435; 95% CI: 0.253–0.476; p = 0.003) and dried fruit intake (OR: 0.522; 95% CI: 0.311–0.875; p = 0.014) were all recognized as significantly protective factors. In addition, we have also reached positive conclusions in the weighted median model of cereal intake (OR: 0.299; 95% CI: 0.147–0.607; p = 0.001), oily fish intake (OR: 0.572; 95% CI: 0.332–0.985; p = 0.044) and the MR Egger model of cheese intake (OR: 5.490; 95% CI: 1.325–22.751; p = 0.022), although the IVW results of oily fish intake and cheese intake were non-significant. This study also found that other dietary factors were not associated with CRC. More specific analysis results are in Figures 2, 3 and Table 2.

Figure 2

Figure 3

Table 2

ExposureMR methodnSNPOR (95% CI)p-val
Alcoholic drinks per weekMR Egger321.455 (0.809–2.618)0.220
Weighted median321.417 (0.886–2.265)0.146
Inverse variance weighted321.565 (1.068–2.293)0.022
Simple mode321.598 (0.433–5.89)0.487
Weighted mode321.472 (0.889–2.435)0.143
Processed meat intakeMR Egger230.347 (0.016–7.337)0.504
Weighted median231.300 (0.572–2.955)0.531
Inverse variance weighted231.732 (0.945–3.174)0.076
Simple mode231.181 (0.259–5.373)0.832
Weighted mode231.160 (0.224–6.011)0.861
Poultry intakeMR Egger73.263 × 1012 (1.000 × 10−3–1.450 × 1030)0.172
Weighted median71.540 (0.250–9.473)0.641
Inverse variance weighted71.951 (0.505–7.542)0.333
Simple mode70.941 (0.063–14.106)0.966
Weighted mode70.941 (0.059–14.976)0.967
Beef intakeMR Egger153.686 (0.038–361.881)0.587
Weighted median152.244 (0.736–6.839)0.155
Inverse variance weighted152.036 (0.872–4.757)0.100
Simple mode151.592 (0.174–14.579)0.687
Weighted mode151.667 (0.208–13.325)0.638
Non-oily fish intakeMR Egger110.066 (0.000–15.012)0.352
Weighted median110.295 (0.077–1.125)0.074
Inverse variance weighted110.286 (0.095–0.860)0.026
Simple mode110.263 (0.030–2.272)0.253
Weighted mode110.233 (0.036–1.498)0.156
Oily fish intakeMR Egger590.349 (0.071–1.712)0.200
Weighted median590.572 (0.332–0.985)0.044
Inverse variance weighted590.700 (0.479–1.022)0.064
Simple mode590.488 (0.152–1.572)0.234
Weighted mode590.445 (0.146–1.355)0.159
Pork intakeMR Egger141.123 (0.001–1538.854)0.975
Weighted median141.687 (0.407–6.987)0.471
Inverse variance weighted141.208 (0.419–3.483)0.727
Simple mode141.955 (0.182–21.001)0.589
Weighted mode141.626 (0.12–22.026)0.720
Lamb/mutton intakeMR Egger311.660 (0.093–29.637)0.733
Weighted median310.749 (0.293–1.919)0.547
Inverse variance weighted310.562 (0.280–1.129)0.105
Simple mode310.964 (0.158–5.879)0.968
Weighted mode310.854 (0.152–4.815)0.860
Bread intakeMR Egger290.428 (0.024–7.732)0.570
Weighted median291.006 (0.461–2.193)0.989
Inverse variance weighted290.842 (0.469–1.510)0.563
Simple mode290.893 (0.232–3.438)0.870
Weighted mode291.066 (0.303–3.752)0.922
Cheese intakeMR Egger625.490 (1.325–22.751)0.022
Weighted median620.956 (0.588–1.554)0.855
Inverse variance weighted620.876 (0.613–1.252)0.467
Simple mode620.872 (0.263–2.897)0.824
Weighted mode620.976 (0.364–2.619)0.962
Cooked vegetable intakeMR Egger170.079 (0.000–1635.243)0.624
Weighted median170.778 (0.217–2.788)0.700
Inverse variance weighted170.672 (0.271–1.667)0.391
Simple mode171.068 (0.083–13.724)0.960
Weighted mode170.879 (0.086–8.949)0.914
Tea intakeMR Egger400.960 (0.425–2.168)0.922
Weighted median400.793 (0.470–1.337)0.385
Inverse variance weighted400.777 (0.548–1.100)0.155
Simple mode400.955 (0.346–2.633)0.929
Weighted mode400.922 (0.503–1.691)0.795
Fresh fruit intakeMR Egger510.539 (0.054–5.399)0.602
Weighted median510.474 (0.182–1.236)0.127
Inverse variance weighted510.513 (0.273–0.964)0.038
Simple mode510.143 (0.019–1.099)0.067
Weighted mode510.226 (0.045–1.132)0.076
Cereal intakeMR Egger390.509 (0.047–5.530)0.583
Weighted median390.299 (0.147–0.607)0.001
Inverse variance weighted390.435 (0.253–0.746)0.003
Simple mode390.250 (0.054–1.156)0.084
Weighted mode390.232 (0.058–0.923)0.045
Salad/raw vegetable intakeMR Egger184.179 (0.013–1388.484)0.636
Weighted median180.677 (0.168–2.724)0.583
Inverse variance weighted180.395 (0.127–1.227)0.108
Simple mode181.814 (0.107–30.647)0.685
Weighted mode182.346 (0.173–31.775)0.530
Coffee intakeMR Egger381.099 (0.420–2.876)0.849
Weighted median381.080 (0.583–1.998)0.808
Inverse variance weighted381.029 (0.647–1.637)0.903
Simple mode381.202 (0.357–4.044)0.768
Weighted mode381.043 (0.522–2.082)0.906
Dried fruit intakeMR Egger410.939 (0.088–10.052)0.959
Weighted median410.540 (0.262–1.113)0.095
Inverse variance weighted410.522 (0.311–0.875)0.014
Simple mode410.544 (0.125–2.356)0.420
Weighted mode410.544 (0.130–2.272)0.409

The results of Mendelian randomization analyses.

p-val, p-value; SNP, single nucleotide polymorphism; OR, odds ratio; CI, confidence interval.

3.3 Sensitivity analysis

Meanwhile, no heterogeneity was discovered in Cochran’s Q tests (p > 0.05 for all the consequences). MR-Egger intercept test indicated that except for the causality calculation between cheese intake and CRC, no statistically significant horizontal pleiotropy was observed in other remaining research (Figure 3 and Supplementary Table S3). Leave-one-out results suggested that no particular SNP could independently affect the MR positive conclusions (Figure 4). All the sensitivity analyses ensured the reliability of our results.

Figure 4

4 Discussion

We executed two-sample MR analyses utilizing large-scale GWAS summary statistics. These analyses observed genetic evidence for a causal association of CRC risk with 17 genetically predicted diet-related factors. Specifically, we noticed suggestive evidence that weekly alcoholic drinks may elevate CRC risk while a higher intake of non-oily fish, cereals, fresh fruit, and dried fruit may reduce risk. Apart from these five exposures, there was no evidence that other dietary factors affected CRC risk significantly. Clarifying these relationships had a vital impact on developing nutritional recommendations for CRC management and prevention.

The relationship between dietary factors and CRC remains controversial. Previously, some observational studies indicated that alcohol intake was a risk factor for colorectal cancer. For instance, a nested case-control study in South Asia revealed that current or former drinkers had a higher risk of CRC (OR: 5.4; 95% CI: 1.1–27.8; p = 0.043) (30). Similar conclusions were reported from other methods and regions (31, 32). However, a previous European MR study found no evidence of a pronounced relationship (OR: 1.60; 95% CI: 0.85–3.04; p = 0.146) (33). Whereas a total of 3 IVs representing weekly alcohol consumption were utilized, and only 0.2% of the genetic variation could be explained, which might lead to a weak statistical power and the absence of robustness. Our study, using 32 SNPs in total, preliminarily demonstrated that alcohol drinks per week was causally associated with about a 56.5% increase in the risk of CRC in European individuals. Some experimental evidence indicated that alcohol might result in the development of CRC by disrupting the composition of gut microbacteria. The possible acetaldehyde accumulation in the Ruminococcus and Coriobacterium located in the colorectum would contribute to mutagenesis and the enablement of carcinogenesis (34). Simultaneously, alcohol metabolites might trigger DNA-adduct formation, lipid peroxidation, and oxidative stress, leading to the initiation of cancer-promoting cascades (35). Additionally, an epigenetic analysis and a gene-alcohol interaction analysis revealed that alcohol consumption could affect DNA methylation by regulating the expression of the COLCA1/COLCA2 gene, which would also increase CRC risk (36). Further investigations are necessary to identify the role of alcohol intake in the genetic and metabolic effects of CRC.

The consequences are also inconsistent between fruit intake and the CRC risk. A European prospective investigation covering 2,819 incident CRC cases has shown that fruit consumption was inversely linked with CRC. The CRC risk was compared between the highest and the lowest EPIC-wide quintile of consumption over an 8.8-year follow-up (HR: 0.86; 95% CI: 0.75–1.00; p trend = 0.04) (37). Similarly, a cohort study on Chinese males obtained the same result (HR: 0.67; 95% CI: 0.48–0.95; p trend = 0.03) (38). On the contrary, a meta-analysis containing 16 cohort studies indicated the absence of significant association (39). Notably, the aforementioned conclusions might not be reliable due to the inherent drawbacks of the observational study design. Removing the underlying confounding factors and focusing on the fresh and dried fruit separately, our MR analyses suggested both fresh fruit (OR: 0.513; 95% CI: 0.273–0.964; p = 0.038) and dried fruit intake (OR: 0.522; 95% CI: 0.311–0.875; p = 0.014) were genetically correlated with a lower risk of CRC. The casual relationship may be attributable to several physiological mechanisms. Specifically, apigenin, a flavonoid that widely existed in fruits, targeted the K433 site of PKM2, thus restricted the glycolysis of HCT-8 and LS-174T cells, thereby serving the crucial function of anti-CRC in vivo and in vitro and markedly attenuating tumor growth in the meantime (40). Moreover, anthocyanins are phenolic pigments that give red and purple fruits their vivid colors. It has been demonstrated to protect against CRC by suppressing the activity and expression of DNA methyltransferase enzymes (DNMT1 and DNMT3B) and demethylating WNT upstream regulators (CDKN2A, SFRP2, SFRP5, and WIF1) (41). Further explorations were necessary to confirm the existence of the causality and investigate the concrete mechanism.

To date, the role of cereal intake in CRC has been widely studied, and a certain amount of epidemiological studies have yielded similar conclusions. A meta-analysis containing 7 European studies suggested a 10% decreased risk of CRC for each 10 g/day intake of cereal and more obvious reductions with higher intake (42). A prospective study of the UK Biobank deduced that intake of fiber from breakfast cereals was a statistically protective factor to CRC (HR: 0.86, 95% CI: 0.76–0.98, p trend = 0.005) with the multivariable model (43). Our results further confirmed the significant causal effect of cereal consumption (OR: 0.435; 95% CI: 0.253–0.746; p = 0.003) against the development of CRC. Mechanism studies reported that cereal foods could increase stool bulk, dilute fecal carcinogens, and decrease transit time. These procedures could offer the lining of the colorectum protective effects against carcinogens (44), which supported our discovery. Specifically, cereal foods’ regulatory effects on CRC development were mediated by activating AHR and GPCRs and inhibiting STAT3 phosphorylation (45). Analogically, other cereal components, including vitamins, phytoestrogens, and trace minerals, have also been associated with a lower risk of CRC (46). More underlying anticarcinogenic mechanisms of high levels of cereal intake could be investigated in the future.

In contrast, there is only a limited number of clinical studies focusing on non-oily fish and CRC. A large European cohort investigation observed an inverse association with CRC incidence (HR: 0.91; 95% CI: 0.83–1.00; p trend = 0.016) (47), which was compatible with our present study (OR: 0.286; 95% CI: 0.095–0.860, p = 0.0026). Additionally, pathophysiological evidence proposed that the ω-3 polyunsaturated fatty acids (PUFAs) contained in the fish might regulate eicosanoid metabolism (48). It was revealed that eicosapentaenoic acid (EPA), which is a type of ω-3 PUFAs, could lead to a decrease in the number and size of colorectal tumors by inhibiting COX-2, reducing β-catenin nuclear translocation and increasing apoptosis (49). ω-3 PUFAs could also promote a higher gut microbial diversity, thus ameliorating the body’s metabolic and immune functions and eventually reducing the CRC risk (34, 50). Subsequent high-quality analyses are required to deduce potential causalities and biological mechanisms.

Notably, some food of animal origins, such as dairy products and eggs, are susceptible to contamination by persistent organic pollutants (POPs), including polychlorinated dibenzo-p-dioxins (PCDDs), polychlorinated dibenzofurans (PCDFs) (51, 52), and polychlorinated biphenyls (PCBs) (53, 54). Long-term exposure to those POPs could damage the immune system and interfere with endocrine functions, thus causing a range of adverse health effects, especially cancer (51–54). Given that dietary intake is the primary route of exposure for humans, contaminated food of animal origin poses a significant risk to public health. Possible interventions, including vigilant monitoring, improved agricultural practices, regulatory enforcement, and public education, should be taken to reduce the risks associated with these contaminants.

There are multiple critical advantages of this work as follows: for all we know, this is the first work to elucidate the causal associations between CRC and diet-related factors by the two-sample MR method. This method addressed the debate of the prior epidemiologic studies and avoided the inherent deficiencies of previous traditional observational research, such as reverse causality and confounders. It also provided novel insights and methods for assessing the health benefits associated with dietary configurations. Secondly, benefiting from the large-scale GWAS database, the massive sample size of our analyses and the solid statistical evaluation effect of each IV (F-statistic >10) guaranteed the statistical validity of the current study. Moreover, we restricted the participants of this study to European-descent individuals, which minimized the potential bias induced by population stratification. Eventually, 5 MR methods and diverse sensitivity analyses were applied to assess the consistency of causal effects and obtained similar results, ensuring the robustness and stability of our discovery.

Some possible limitations in this study should also be considered. First, Mendel’s second law is not universally applicable to all genetic variants because not all genes determining traits are isolated independently. The inherent presence of developmental compensation bias also contributes to the potential inaccuracy of Mendelian randomization studies. Second, all analyses conducted in the current study were merely based on the European participants. Thus, it remained to be seen whether our results could be extrapolated to non-European populations. Third, due to the lack of classified population GWAS data for different sexes and ages, we could not execute a sex- or age-stratified analysis. Specifically, owing to the limited details provided by the original research, it was difficult to predict the generalizability of the study results across different exposure periods and levels. Analogically, diet-related information obtained from surveys may be prone to recall bias, which could possibly render our results unreliable. Additionally, given the complexity of dietary habits, we were unable to distinguish the impacts of diverse dietary combinations. Hence, it was challenging to identify the specific role of these interested dietary factors in the etiology and pathogenesis of CRC. Further investigation will focus on conducting more comprehensive studies to gather high-quality evidence regarding the idiographic mechanisms through which dietary factors affect CRC risk. This involves expanding the scope of research to include a broader range of dietary factors, identifying potential biomarkers that could help in understanding the effect of diet on CRC development, exploring genetic predispositions that may modify the impact of dietary factors, and longitudinal studies to track dietary habits over time and their direct correlation with CRC incidence.

5 Conclusion

Based on the GWAS summary data of CRC and European dietary habits, this study was implemented to identify the potential associations of colorectal cancer with 17 dietary factors using genetic instruments. The causal relationship between alcoholic drinks per week and an increased risk of CRC and the inverse causality of non-oily fish intake, cereal intake, fresh fruit, and dried fruit intake with CRC were determined by performing the two-sample MR analyses. The current study is the first to build the link between comprehensive diet-related factors and CRC at the genetic level, offering novel clues for interpreting the genetic etiology of CRC and replenishing new perspectives for managing gastrointestinal diseases. The result also prompts future explorations, including longitudinal studies and nutritional interventions, highlights the importance of interdisciplinary collaboration for clinical diagnostics, comprehensive patient care, and genetic counseling and education, and helps develop public health recommendations and tailored nutrition and prevention strategies.

Statements

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Author contributions

XZ: Formal analysis, Project administration, Conceptualization, Data curation, Investigation, Methodology, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. ZW: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. XW: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. BX: Formal analysis, Investigation, Resources, Writing – review & editing, Data curation, Project administration, Validation, Writing – original draft. PH: Conceptualization, Investigation, Resources, Software, Writing – original draft, Formal analysis, Visualization. YY: Conceptualization, Supervision, Project administration, Writing – review & editing. SH: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Software, Supervision, Validation, Visualization, Writing – original draft. MR: Methodology, Resources, Software, Validation, Writing – review & editing, Conceptualization, Funding acquisition, Supervision.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The research was funded by the Key Research and Development Program of Shaanxi (No. 2021ZDLSF02-06) and the Shaanxi Provincial International Scientific and Technological Cooperation Projects (No. 2018 KW-014).

Acknowledgments

Special thanks to the IEU open GWAS project developed by the MRC Integrative Epidemiology Unit at the University of Bristol for extracting related GWAS summary-level data from published articles, the UK Biobank, and the European Bioinformatics Institute. The authors also acknowledge all the investigators and participants of the relevant study.

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/fnut.2024.1388732/full#supplementary-material

References

  • 1.

    SungHFerlayJSiegelRLLaversanneMSoerjomataramIJemalAet al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2021) 71:20949. doi: 10.3322/caac.21660

  • 2.

    SiegelRLMillerKDGoding SauerAFedewaSAButterlyLFAndersonJCet al. Colorectal cancer statistics, 2020. CA Cancer J Clin. (2020) 70:14564. doi: 10.3322/caac.21601

  • 3.

    McDougallJABlairCKWigginsCLGoodwinMBChiuVKRajputAet al. Socioeconomic disparities in health-related quality of life among colorectal cancer survivors. J Cancer Surviv. (2019) 13:45967. doi: 10.1007/s11764-019-00767-9

  • 4.

    HarringtonCBHansenJAMoskowitzMToddBLFeuersteinM. It’s not over when it’s over: long-term symptoms in cancer survivors—a systematic review. Int J Psychiatry Med. (2010) 40:16381. doi: 10.2190/PM.40.2.c

  • 5.

    RabeneckLChiuH-MSenoreC. International perspective on the burden of colorectal cancer and public health effects. Gastroenterology. (2020) 158:44752. doi: 10.1053/j.gastro.2019.10.007

  • 6.

    SharmaRAbbasi-KangevariMAbd-RabuRAbidiHAbu-GharbiehEAcunaJMet al. Global, regional, and national burden of colorectal cancer and its risk factors, 1990–2019: a systematic analysis for the global burden of disease study 2019. Lancet Gastroenterol Hepatol. (2022) 7:62747. doi: 10.1016/S2468-1253(22)00044-9

  • 7.

    Fernandez-RozadillaCTimofeevaMChenZLawPThomasMSchmitSet al. Deciphering colorectal cancer genetics through multi-omic analysis of 100, 204 cases and 154,587 controls of European and East Asian ancestries. Nat Genet. (2023) 55:8999. doi: 10.1038/s41588-022-01222-9

  • 8.

    AnSParkS. Association of physical activity and sedentary behavior with the risk of colorectal cancer. J Korean Med Sci. (2022) 37:e158. doi: 10.3346/jkms.2022.37.e158

  • 9.

    BotteriEBorroniESloanEKBagnardiVBosettiCPeveriGet al. Smoking and colorectal cancer risk, overall and by molecular subtypes: a meta-analysis. Rev Gastroenterol. (2020) 115:19409. doi: 10.14309/ajg.0000000000000803

  • 10.

    WangFUgaiTHarukiKWanYAkimotoNArimaKet al. Healthy and unhealthy plant-based diets in relation to the incidence of colorectal cancer overall and by molecular subtypes. Clin Transl Med. (2022) 12:e893. doi: 10.1002/ctm2.893

  • 11.

  • 12.

    ChanDSMLauRAuneDVieiraRGreenwoodDCKampmanEet al. Red and processed meat and colorectal cancer incidence: meta-analysis of prospective studies. PLoS One. (2011) 6:e20456. doi: 10.1371/journal.pone.0020456

  • 13.

    SongMGarrettWSChanAT. Nutrients, foods, and colorectal cancer prevention. Gastroenterology. (2015) 148:124460.e16. doi: 10.1053/j.gastro.2014.12.035

  • 14.

    TuanJChenY-X. Dietary and lifestyle factors associated with colorectal cancer risk and interactions with microbiota: fiber, red or processed meat and alcoholic drinks. Gastrointest Tumors. (2015) 3:1724. doi: 10.1159/000442831

  • 15.

    ChenZLiuCWuJKongF. Association of coffee and caffeine consumption with risk and prognosis of endometrial cancer and its subgroups: a Mendelian randomization. Front Nutr. (2023) 10:1291355. doi: 10.3389/fnut.2023.1291355

  • 16.

    ZengYCaoSYangH. Causal associations between dried fruit intake and cardiovascular disease: a Mendelian randomization study. Front Cardiovasc Med. (2023) 10:1080252. doi: 10.3389/fcvm.2023.1080252

  • 17.

    ZhengJBairdDBorgesM-CBowdenJHemaniGHaycockPet al. Recent developments in Mendelian randomization studies. Curr Epidemiol Rep. (2017) 4:33045. doi: 10.1007/s40471-017-0128-6

  • 18.

    ShengJLiuJChanKHK. Evaluating the causal effects of gestational diabetes mellitus, heart disease, and high body mass index on maternal Alzheimer’s disease and dementia: multivariable Mendelian randomization. Front Genet. (2022) 13:833734. doi: 10.3389/fgene.2022.833734

  • 19.

    DuZGuoSSunYZhouQJiangYShenYet al. Causal relationships between dietary habits and five major mental disorders: a two-sample Mendelian randomization study. J Affect Disord. (2023) 340:60715. doi: 10.1016/j.jad.2023.08.098

  • 20.

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

  • 21.

    WoottonRELawnRBMillardLACDaviesNMTaylorAEMunafòMRet al. Evaluation of the causal effects between subjective wellbeing and cardiometabolic health: Mendelian randomisation study. BMJ. (2018) 362:k3788. doi: 10.1136/bmj.k3788

  • 22.

    LiCNiuMGuoZLiuPZhengYLiuDet al. A mild causal relationship between tea consumption and obesity in general population: a two-sample Mendelian randomization study. Front Genet. (2022) 13:795049. doi: 10.3389/fgene.2022.795049

  • 23.

    WangXGlubbDMO’MaraTA. Dietary factors and endometrial cancer risk: a Mendelian randomization study. Nutrients. (2023) 15:603. doi: 10.3390/nu15030603

  • 24.

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

  • 25.

    YangJWuWAmierYLiXWanWYuX. Causal relationship of genetically predicted circulating micronutrients levels with the risk of kidney stone disease: a Mendelian randomization study. Front Nutr. (2023) 10:1132597. doi: 10.3389/fnut.2023.1132597

  • 26.

    FuSZhangLMaFXueSSunTXuZ. Effects of selenium on chronic kidney disease: a Mendelian randomization study. Nutrients. (2022) 14:4458. doi: 10.3390/nu14214458

  • 27.

    BowdenJDel GrecoMFMinelliCZhaoQLawlorDASheehanNAet al. Improving the accuracy of two-sample summary-data Mendelian randomization: moving beyond the NOME assumption. Int J Epidemiol. (2018) 48:72842. doi: 10.1093/ije/dyy258

  • 28.

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

  • 29.

    WangZChenMWeiYZZhuoCGXuHFLiWDet al. The causal relationship between sleep traits and the risk of schizophrenia: a two-sample bidirectional Mendelian randomization study. BMC Psychiatry. (2022) 22:399. doi: 10.1186/s12888-022-03946-8

  • 30.

    SamarakoonYMGunawardenaNSPathiranaA. Behavioral, familial and comorbid illness risk factors of colorectal cancer: a case control study. Ceylon Med J. (2018) 63:1138. doi: 10.4038/cmj.v63i3.8725

  • 31.

    WangJWangHChenYHaoPZhangY. Alcohol ingestion and colorectal neoplasia: a meta-analysis based on a Mendelian randomization approach. Color Dis. (2011) 13:e718. doi: 10.1111/j.1463-1318.2010.02530.x

  • 32.

    LiYYeDZhouWLiuBMaoYSunX. Alcohol consumption and colorectal cancer risk: a Mendelian randomization study. Front Genet. (2022) 13:967229. doi: 10.3389/fgene.2022.967229

  • 33.

    CornishAJLawPJTimofeevaMPalinKFarringtonSMPallesCet al. Modifiable pathways for colorectal cancer: a Mendelian randomisation analysis. Lancet Gastroenterol Hepatol. (2020) 5:5562. doi: 10.1016/S2468-1253(19)30294-8

  • 34.

    SongMChanAT. Environmental factors, gut microbiota, and colorectal cancer prevention. Clin Gastroenterol Hepatol. (2019) 17:27589. doi: 10.1016/j.cgh.2018.07.012

  • 35.

    RossiMJahanzaib AnwarMUsmanAKeshavarzianABishehsariF. Colorectal cancer and alcohol consumption—populations to molecules. Cancers. (2018) 10:38. doi: 10.3390/cancers10020038

  • 36.

    ZhouXWangLXiaoJSunJYuLZhangHet al. Alcohol consumption, DNA methylation and colorectal cancer risk: results from pooled cohort studies and Mendelian randomization analysis. Int J Cancer. (2022) 151:8394. doi: 10.1002/ijc.33945

  • 37.

    van DuijnhovenFJBBueno-De-MesquitaHBFerrariPJenabMBoshuizenHCRosMMet al. Fruit, vegetables, and colorectal cancer risk: the European prospective investigation into cancer and nutrition 123. Am J Clin Nutr. (2009) 89:144152. doi: 10.3945/ajcn.2008.27120

  • 38.

    VogtmannEXiangY-BLiH-LLevitanEBYangGWaterborJWet al. Fruit and vegetable intake and the risk of colorectal cancer: results from the Shanghai Men’s Health Study. Cancer Causes Control. (2013) 24:193545. doi: 10.1007/s10552-013-0268-z

  • 39.

    HuxleyRRAnsary-MoghaddamACliftonPCzernichowSParrCLWoodwardM. The impact of dietary and lifestyle risk factors on risk of colorectal cancer: a quantitative overview of the epidemiological evidence. Int J Cancer. (2009) 125:17180. doi: 10.1002/ijc.24343

  • 40.

    ShiJJiXShanSZhaoMBiCLiZ. The interaction between apigenin and PKM 2 restrains progression of colorectal cancer. J Nutr Biochem. (2023) 121:109430. doi: 10.1016/j.jnutbio.2023.109430

  • 41.

    de Paula do NascimentoRda Fonseca MachadoAP. The preventive and therapeutic effects of anthocyanins on colorectal cancer: a comprehensive review based on up-to-date experimental studies. Food Res Int. (2023) 170:113028. doi: 10.1016/j.foodres.2023.113028

  • 42.

    AuneDChanDSMLauRVieiraRGreenwoodDCKampmanEet al. Dietary fibre, whole grains, and risk of colorectal cancer: systematic review and dose-response meta-analysis of prospective studies. BMJ. (2011) 343:d6617. doi: 10.1136/bmj.d6617

  • 43.

    BradburyKEMurphyNKeyTJ. Diet and colorectal cancer in UK Biobank: a prospective study. Int J Epidemiol. (2019) 49:24658. doi: 10.1093/ije/dyz064

  • 44.

    SlavinJL. Mechanisms for the impact of whole grain foods on cancer risk. J Am Coll Nutr. (2000) 19:300S7S. doi: 10.1080/07315724.2000.10718964

  • 45.

    ZhangBXuYLiuSLvHHuYWangYet al. Dietary supplementation of foxtail millet ameliorates colitis-associated colorectal cancer in mice via activation of gut receptors and suppression of the stat 3 pathway. Nutrients. (2020) 12:2367. doi: 10.3390/nu12082367

  • 46.

    CotterchioMBoucherBAMannoMGallingerSOkeyAHarperP. Dietary phytoestrogen intake is associated with reduced colorectal cancer risk. J Nutr. (2006) 136:304653. doi: 10.1093/jn/136.12.3046

  • 47.

    AglagoEKHuybrechtsIMurphyNCasagrandeCNicolasGPischonTet al. Consumption of fish and long-chain N-3 polyunsaturated fatty acids is associated with reduced risk of colorectal cancer in a large European cohort. Clin Gastroenterol Hepatol. (2020) 18:65466.e6. doi: 10.1016/j.cgh.2019.06.031

  • 48.

    YangPJiangYFischerSM. Prostaglandin E3 metabolism and cancer. Cancer Lett. (2014) 348:111. doi: 10.1016/j.canlet.2014.03.010

  • 49.

    IrúnPLanasAPiazueloE. Omega-3 polyunsaturated fatty acids and their bioactive metabolites in gastrointestinal malignancies related to unresolved inflammation. A review. Front Pharmacol. (2019) 10:852. doi: 10.3389/fphar.2019.00852

  • 50.

    PiazziGProssomaritiABaldassarreMMontagnaCVitaglionePFoglianoVet al. A Mediterranean diet mix has chemopreventive effects in a murine model of colorectal cancer modulating apoptosis and the gut microbiota. Front Oncol. (2019) 9:140. doi: 10.3389/fonc.2019.00140

  • 51.

    GiannicoOVBaldacciSDesianteFBasileFCFrancoEFragnelliGRet al. PCDD/Fs and PCBs in Mytilus galloprovincialis from a contaminated area in Italy: the role of mussel size, temperature and meteorological factors. Food Addit Contam A. (2022) 39:112335. doi: 10.1080/19440049.2022.2059108

  • 52.

    GiannicoOVBaldacciSBasileFCPellegrinoADesianteFFrancoEet al. PCDD/Fs and PCBs in hen eggs from a contaminated area in Italy: a 9 years spatio-temporal monitoring study. Food Addit Contam A. (2023) 40:294304. doi: 10.1080/19440049.2022.2157051

  • 53.

    GiannicoOVDesianteFBasileFCFrancoEBaldacciSFragnelliGRet al. Dioxins and PCBs contamination in mussels from Taranto (Ionian Sea, Southern Italy): a seven years spatio-temporal monitoring study. Ann Ist Super Sanita. (2020) 56:45261. doi: 10.4415/ann_20_04_07

  • 54.

    GiannicoOVFragnelliGRBaldacciSDesianteFPellegrinoABasileFCet al. Dioxins and PCBs contamination in milk and dairy products from province of Taranto (Puglia Region, Southern Italy): a six years spatio-temporal monitoring study. Ann Ist Super Sanita. (2021) 57:2338. doi: 10.4415/ann_21_03_06

Summary

Keywords

dietary factors, colorectal cancer, Mendelian randomization, causality, GWAS

Citation

Zhang X, Wu Z, Wang X, Xin B, Hu P, Yin Y, He S and Ren M (2024) Causal associations between dietary factors and colorectal cancer risk: a Mendelian randomization study. Front. Nutr. 11:1388732. doi: 10.3389/fnut.2024.1388732

Received

20 February 2024

Accepted

19 April 2024

Published

01 May 2024

Volume

11 - 2024

Edited by

Roberta Zupo, University of Bari Aldo Moro, Italy

Reviewed by

Orazio Valerio Giannico, Local Health Authority of Taranto, Italy

Youjie Zeng, Central South University, China

Updates

Copyright

*Correspondence: Shuixiang He, ; Mudan Ren,

†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