Introduction

Hypermethylation of CpG islands in gene promoters is an important epigenetic alteration involved in carcinogenesis [1]. In colorectal cancer (CRC), the CpG island methylator phenotype (CIMP) is characterized by frequent promoter CpG island hypermethylation [2]. However, little is known about potential determinants of this type of aberrant DNA methylation in CRC.

Folate and methionine are dietary methyl group donors that may be hypothesized to influence DNA methylation, whereas vitamins B2 and B6 potentially modulate the bioavailability of methyl groups [3, 4]. Low folate status or intake was suggested to decrease genomic methylation [58], while folate supplementation resulted in increased global DNA methylation in the colonic mucosa [9]. Adequate methyl donor intake possibly also prevents aberrant CpG island promoter hypermethylation. In this respect, a weak inverse association with gene promoter hypermethylation was suggested [10], although methyl donor intake and alcohol consumption, which may reduce the bioavailability of folate, were not associated with CIMP in CRC [11]. Conversely, folate supplementation was suggested to increase promoter hypermethylation of multiple genes in colorectal mucosa [12], and circulating folate concentration was associated with increased gene promoter hypermethylation in colorectal tumors [13]. In addition, high vitamin B6 intake may be associated with increased MutL homologue 1 (MLH1) promoter methylation in CRC [14]. Apparently, the precise effect of methyl group bioavailability on gene promoter hypermethylation is still unclear and should be investigated further.

A potential effect of methyl donor intake on DNA methylation may be modified by polymorphisms in folate metabolizing enzymes. For example, the catalytic activity of the methylene tetrahydrofolate reductase (MTHFR) enzyme may be reduced in individuals carrying rare variants of the MTHFR c.665C > T (rs1801133) and c.1286A > C (rs1801131) polymorphisms [15, 16], which were also associated with the CIMP phenotype in colorectal cancer [1719]. We previously observed inverse associations between methionine synthase (MTR) c.2756A > G (rs1805087) and CIMP, and between methionine synthase reductase (MTRR) c.66A > G (rs1801394) with MLH1 hypermethylation [20]. Other enzymes involved in epigenetic regulation of gene expression are DNA methyltransferases (DNMTs) and histone methyltransferases (HMTs). However, whether an influence of methyl donor intake is modified by polymorphisms in such epigenetic regulators has not previously been studied in relation to CRC.

Here, we aimed to investigate associations between dietary folate, methionine, vitamins B2 and B6 with overall CRC, and risk of CRCs harboring CIMP, accounting for the occurrence of any, or combinations of rare variants of folate metabolizing enzymes MTHFR, MTR and MTRR, the DNA methyltransferase DNMT3B, and histone methyltransferases Euchromatin histone methyltransferase 1 (EHMT1), Euchromatin histone methyltransferase 2 (EHMT2) and PR domain zinc finger protein 2 (PRDM2) in the Netherlands Cohort Study on diet and cancer.

Methods

Study population

The participants of this study were incident CRC patients from the Netherlands Cohort Study on diet and cancer (NLCS), which has been described in detail elsewhere [21]. Briefly, this prospective cohort study was initiated in September 1986 and includes 58,279 men and 62,573 women aged 55–69 years and free of disease at baseline. The cohort is followed for cancer occurrence by annual record linkage to the Netherlands Cancer Registry (NCR) and to the Pathologisch Anatomisch Landelijk Geautomatiseerd Archief (PALGA), a nationwide network and registry of histopathology and cytopathology reports [22, 23]. A subcohort of 5,000 subjects was randomly selected after baseline exposure measurement, to estimate accumulation of person-time in the cohort through biennial follow-up of vital status. Cases with prevalent cancer other than non-melanoma skin cancer were excluded from this subcohort, which left 4,774 men and women eligible for analysis.

Food frequency questionnaire

At baseline, participants filled out a self-administered, 150-item semi-quantitative food frequency questionnaire (FFQ), which concentrated on habitual consumption of food and beverages during the year preceding the start of the study, and also contained questions about age, sex, body weight and length, smoking status and family history of CRC. Daily mean nutrient intakes were calculated as the cumulated product of the frequencies and portion sizes of all food items and their tabulated nutrient contents from the Dutch Food Composition Table (NEVO table, 1986) [24]. The questionnaire was validated through comparison with a 9-day diet record [25]. Reproducibility and stability of dietary habits were determined by five annually repeated measurements [26]. In order to minimize observer bias in coding and interpretation of the data, questionnaire data were key-entered twice for all incident cases in the cohort and for all subcohort members in a blinded manner with respect to case/subcohort status.

Folate data were derived from a validated liquid chromatography trienzyme method [27] used to analyze the 125 most important Dutch foods contributing to folate intake [28]. Dietary supplement data were also obtained via the food frequency questionnaire. However, the use of B-vitamin supplements was low (7%) and folic acid was generally not included in these supplements in the Netherlands in the late 1980s. Therefore, folic acid supplement use most likely plays a very minor role in our study population, and supplement use was not further accounted for in the analyses.

Sample collection

Subcohort members still alive in December 2000 (n = 3,579) were contacted and asked to collect mouth swabs, of whom 1,929 (54%) responded and returned the mouth swab with informed consent. In total, DNA could successfully be isolated of 1,829 subcohort members who also had complete follow-up information [20].

Tumor material of the CRC patients was collected after approval by the ethical review boards of Maastricht University, the NCR and PALGA. During a follow-up period of 7.3 years after baseline, 734 incident CRC patients were identified who had an available PALGA report of the lesion as well as a sufficient amount of isolated DNA needed for molecular analyses.

Genotyping analyses

MTHFR (rs1801133 and rs1801131), MTR (rs1805087), MTRR (rs1801394), DNMT3B (rs2424913 and rs406193), EHMT1 (rs4634736), EHMT2 (rs535586) and PRDM2 (rs2235515) genotypes were determined using multiplex polymerase chain reaction (PCR) amplification and single base extension (SBE) reactions as described previously [20, 29]. Genotype data were validated by sequencing of fragments containing specific SNPs, which were similar to the main results for all but one (99.6%) of the 9 SNPs within a subset of 30 samples [20]. Reproducibility of the analysis was established among 93 samples, and we observed that the analyses could be reproduced in 99.5% of these cases [20]. In total, genotyping analyses were successful from 1,736 subcohort members and 659 CRC patients.

Promoter methylation analyses

The CpG island methylator phenotype (CIMP) was defined by promoter hypermethylation of at least 3 out of 5 methylation markers (CACNA1G, IGF2, NEUROG1, RUNX3, and SOCS1), as suggested by Weisenberger et al. [2]. Hypermethylation of the CpG islands of these five CIMP markers and of the MLH1 gene was determined by Methylation Specific PCR (MSP) [30] and described in detail by de Vogel et al [20]. The MSP analyses were successful of 81, 79, 79, 90, 83, and 93% out of the 734 patients for CACNA1G, IGF2, NEUROG1, RUNX3, SOCS1, and MLH1, respectively.

Microsatellite instability

MSI was determined by a pentaplex PCR, using the MSI markers BAT-26, BAT-25, NR-21, NR-22 and NR-24, as described in detail by Suraweera et al. [31]. MSI analyses were successful on 662 (90%) out of the 734 available samples.

Statistical analyses

Cox proportional hazards regression models were used to estimate multivariate-adjusted incidence rate ratios (RR) and corresponding 95% confidence intervals (CI) over tertiles of dietary folate, methionine, vitamins B2 and B6, using the lowest tertiles as reference. Tests for dose response trends over the tertiles of intake were estimated by fitting the ordinal exposure variables as continuous variables and evaluated using the Wald test. Standard errors of the RR were estimated using the robust Huber-White sandwich estimator to account for additional variance introduced by sampling from the cohort [32]. The proportional hazards assumption was tested using the scaled Schoenfeld residuals [33] and by fitting the main determinants as time-dependent variables. The dietary variables were adjusted for total energy intake by calculating nutrient residuals from the regression of nutrient intake on total energy intake, as described by Willett et al. [34]. The analyses were stratified according to genetic status of individuals, i.e. among those homozygous to common genetic variants and among subjects carrying rare alleles. Interactions were tested between dietary folate, methionine, vitamins B2 and B6, and each of the genetic variants. Associations between dietary factors and CRC were also estimated for combinations of genotypes per functional group (i.e. based on the number of rare alleles in any of the folate metabolizing enzymes MTHFR, MTR and MTRR, in the DNA methyltransferase DNMT3B, or in any of the histone methyltransferases EHMT1, EHMT2 and PRDM2.

To investigate whether dietary methyl donors have an effect on promoter hypermethylation in CRC, associations of folate, methionine, vitamins B2, and B6 with the CIMP phenotype were estimated. The associations with MLH1 hypermethylation and MSI were reported previously [14]. Furthermore, it was investigated whether the associations with MLH1 hypermethylation, MSI, or CIMP would be modified by genetic status, by estimating the associations with methylation endpoints within genotypes of folate metabolizing enzymes, DNMT3B and histone methyltransferases.

All models included the co-variates dietary folate, methionine, vitamin B2 and B6 and were additionally adjusted for age, sex, family history of CRC, smoking status, body mass index (BMI), alcohol consumption, and energy intake. After excluding subjects with missing information on these covariates or subjects who did not completely filled out the questionnaire, 1,663 subcohort members and 609 CRC cases remained for statistical analyses. All analyses were performed with the Stata statistical software package (version 10).

Results

CRC risk was estimated over tertiles of folate intake, methionine, vitamins B2 and B6, among subjects homozygous for common alleles and among carriers of rare alleles. Folate or methionine intakes were not associated with CRC within either common homozygotes or within heterozygotes and rare homozygotes of any of the genotypes (Table 1). However, we observed a non-significant inverse association between vitamin B2 intake and CRC risk among subjects with the MTHFR c.665CC (rs1801133) common genotype (RR for the highest versus the lowest tertile of intake = 0.66, p trend = 0.08, Table 2), and an inverse association among subjects with the common GG genotype of PRDM2 G > A (rs2235515, RR = 0.67, p trend = 0.05). In addition, vitamin B2 was associated with reduced CRC risk in individuals carrying the variant allele of DNMT3B C > T (rs2424913, RR = 0.69, p trend = 0.05). Conversely, subjects in the third tertile of vitamin B6 intake were at increased CRC risk when they carried the rare allele of DNMT3B C > T (rs406193, RR = 1.90, p trend = 0.04), or the common allele of PRDM2 G > A (rs2235515, RR = 1.49, p trend = 0.03). However, interactions between these dietary factors and genotypes were not statistically significant.

Table 1 Intake of folate and methionine and CRC risk stratified by genetic status
Table 2 Intake of vitamins B2 and B6 and CRC risk stratified by genetic status
Table 3 Dietary folate, methionine, vitamins B2 and B6 and CRC risk for combinations of genotypes DNA methyltransferase 3B

We also investigated the associations between methyl donor intake and CRC risk according to the number of rare alleles within each functional group (i.e. folate metabolizing enzymes, DNMT3B and histone methyltransferases). It appeared that methionine was inversely associated with CRC if subjects were homozygous to both of the common variants of the DNMT3B rs2424913 and rs406193 C > T SNPs (RR = 0.44, p trend = 0.05, P interaction = 0.07, Table 3). Moreover, relatively high vitamin B2 intake was associated with reduced CRC risk in subjects carrying less than one rare variant of folate metabolizing enzymes (RR = 0.30, p trend = 0.005, P interaction = 0.36, Table 4). No dietary associations were observed according to the number of rare alleles in the studied histone methyltransferases (Table 5).

Table 4 Dietary folate, methionine, vitamins B2 and B6 and CRC risk for combinations of genotypes in folate metabolizing enzymes
Table 5 Dietary folate, methionine, vitamins B2 and B6 and CRC risk for combinations of genotypes in histone methyltransferases

With respect to CpG island promoter hypermethylation, we observed no overall associations between folate, methionine, vitamins B2 or B6 with CIMP (Table 6). Moreover, there were no clear associations between methyl donor intake and CIMP, MLH1 hypermethylation or MSI when accounting for genetic status of individuals (data not shown).

Table 6 Associations of folate, methionine, vitamins B2 and B6 with CIMP in CRC

Discussion

In the current prospective case-cohort study, we observed no clear associations between dietary folate and vitamin B6 with CRC risk when accounting for genetic variants of folate metabolizing enzymes, DNA methyltransferases, or histone methyltransferases. However, relatively high methionine intake may protect against CRC if enzymatic activity of DNMT3B is not affected by two C > T SNPs in its encoding gene. In addition, subjects with high vitamin B2 intake may be at reduced CRC risk in combination with optimal MTHFR activity in individuals homozygous for the common c.665CC (rs1801133) variant, with common PRDM2 GG (rs2235515) genotype and among those with the variant allele of DNMT3B C > T (rs2424913). We observed a strong inverse association between vitamin B2 intake and CRC risk among individuals carrying ≤ 1 rare allele in the combination of any of the folate metabolizing enzymes MTHFR, MTR, or MTRR. There were no associations with the CIMP phenotype overall, or within strata of the studied genotypes.

The MTHFR c.665C > T (rs1801133) polymorphism reduces binding of the MTHFR enzyme to its cofactor flavin adenine dinucleotide (FAD), a metabolite of vitamin B2, resulting in loss of enzymatic activity [15]. The potentially resulting reduced bioavailability of methyl groups may induce DNA hypomethylation in for example blood cells [35, 36] or CpG island promoter hypermethylation in CRC [19, 37]. We observed an inverse association between vitamin B2 and CRC risk, predominantly among subjects homozygous for the MTHFR c.665CC (rs1801133) variant, suggesting that vitamin B2 may maximize the catalytic activity of MTHFR when binding to FAD is optimal. Similarly, it was recently observed that high vitamin B2 plasma concentrations, in combination with MTHFR c.665CC or CT genotypes, may reduce risk of CRA recurrence, whereas such an inverse association was not observed among individuals with the MTHFR c.665TT variant [38].

The rare variant of another MTHFR polymorphism, MTHFR c.1286A > C (rs1801131), may also reduce enzymatic MTHFR activity [16], and was associated with CIMP in colorectal cancer [18], possibly in combination with low folate and methionine intakes and high alcohol consumption [17]. However, we previously observed that this polymorphism was neither associated with overall CRC or with the CIMP phenotype [39], nor when methyl donor intake was accounted for in the current study. Possibly, the use of different panels to identify CIMP-high (a “classic” panel [17] or a new panel [18, 39] which may be more robust [2]) may have contributed to this inconsistency. Moreover, different essays to measure DNA methylation were used, i.e. MSP [17, 39] or a quantitative method [18]. However, in addition to this variety of approaches, it is also important to realize that the one-carbon metabolism is involved in both DNA synthesis as well as DNA methylation, both of which may have an effect on colorectal carcinogenesis [40]. The relative contribution of each of these biological processes in carcinogenesis remains to be established and may not have been similar in the investigated study populations. Furthermore, global DNA hypomethylation and CIMP are possibly inversely associated in CRC [41], and methyl group donors may have an effect on both of these potentially distinct methylation-associated pathways in colorectal carcinogenesis. In this respect, low folate and high alcohol intakes were associated with LINE-1 hypomethylation as an indicator for global DNA hypomethylation [8], which is in agreement with in vivo experimental data [7].

We did not observe associations between methyl donor intake and the CIMP phenotype in CRC, either overall or after stratifying the analyses for the genetic variants of folate metabolizing enzymes or methyltransferases. Moreover, overall associations between methyl donor intake and CIMP in CRC were not observed in another population-based study [11], while in the same cohort, a diet-gene association with CIMP in CRC was observed for only one out of thirteen one-carbon metabolism genes [17]. However, these studies, as well as our study, may have lacked adequate power to demonstrate such associations. The effect of methyl donor intake on gene promoter hypermethylation may indeed be weak though, and to demonstrate whether such an effect is modified by genetic variability of the methyl metabolism requires studies with large numbers of cases. Nonetheless, we did observe an inverse association between vitamin B2 and CRC risk in individuals carrying ≤1 variant allele out of the four studied SNPs of folate metabolizing enzymes, suggesting that the combination of common wild-type genotypes, which possibly results in higher bioavailability of methyl groups, protects against CRC in these people.

The findings of our study may indicate that relatively high methionine intake protects against CRC if enzymatic DNMT3B activity is not affected by two polymorphisms. DNMT3B activity, which may be increased by the DNMT3B C > T (rs2424913) polymorphism [42], was associated with CIMP-high in CRC [43], and with increased risk of various other types of cancer [42, 44, 45]. In addition, experimental research suggested that DNMT3B overexpression induced formation of tumors with promoter hypermethylation [46]. DNMT3B C > T (rs2424913) was also associated with increased colorectal adenoma risk in individuals with low folate and methionine intakes [47], suggesting a nutrient–gene interaction in colorectal carcinogenesis. In view of the function of the DNMT3B enzyme of incorporating methyl groups into DNA, an interaction between methionine intake, DNMT3B polymorphisms and CpG island hypermethylation may be expected, but we did not observe clear associations between methyl donor intake and CIMP, MLH1 hypermethylation or MSI when accounting for DNMT3B genotypes.

The potential protective effects of vitamin B2 or methionine may only be present among individuals with ≤1 polymorphism in folate metabolizing enzymes or among those with common wild-type genotypes of DNMT3B, respectively. This suggests that the occurrence of only one rare variant may be compensated for, but that the combination of several polymorphic genes may lead to disruption of a particular metabolic or regulatory function and to the abolishment of beneficial effects of nutrients. However, we should be careful in drawing definite conclusions because the sample size of our study may have been insufficient to conduct stratified analyses with adequate precision. Moreover, the P-values for interaction were not statistically significant, suggesting the absence of heterogeneity of diet associations with CRC across genotypes. In addition, we conducted several stratified analyses, and these multiple comparisons do not exclude the possibility of reporting chance findings. Nonetheless, although these observations are based on subgroup analyses, and thus have to be interpreted with some caution, this study may indicate that combining genotypes is important to reveal associations of dietary factors with cancer risk. Such an approach has not been followed in previous studies investigating associations between genetic factors and cancer risk, and we recommend that combinations of genotypes should be considered in addition to overall analyses in future studies.

Subgroup analyses in the present study indicated that vitamin B2 and methionine may protect against CRC among individuals who do not carry rare variants of folate metabolizing enzymes and a DNA methyltransferase. However, larger studies are needed to investigate a potential interaction between dietary methyl donor intake, genetic variation of folate metabolizing enzymes and epigenetic regulators, and methylation endpoints in CRC with more precision. Because multiple genes may collectively affect the folate metabolism, combining genotypes of related genes is a useful approach of investigating associations of dietary methyl donors and CRC.