Body Size, Physical Activity and Risk of Colorectal Cancer with or without the CpG Island Methylator Phenotype (CIMP)

Abstract

Background

We investigated how body size and physical activity influence the risk of the CpG island methylator phenotype (CIMP) in colorectal cancer (CRC).

Methods

In the Netherlands Cohort Study (n = 120,852), risk factors were self-reported at baseline in 1986. After 7.3 years of follow-up, 603 cases and 4,631 sub-cohort members were available. CIMP status according to the Weisenberger markers was determined using methylation specific PCR on DNA from paraffin embedded tumor tissue. Hazard rate ratios (HR) and 95% confidence intervals for CIMP (27.7%) and non-CIMP (72.3%) tumors were calculated according to BMI, BMI at age 20, BMI change, trouser/skirt size, height, and physical activity.

Results

BMI modeled per 5 kg/m2 increase was associated with both CIMP and non-CIMP tumors, however, HRs were attenuated when additionally adjusted for trouser/skirt size. Trouser/skirt size, per 2 size increase, was associated with both tumor subtypes, even after adjustment for BMI (CIMP HR: 1.20, 95%CI: 1.01–1.43; non-CIMP HR: 1.14, 95%CI: 1.04–1.28). Height per 5 cm was associated with both tumor sub-types, but HRs were attenuated when adjusted for body weight. BMI at age 20 was positively associated with increased risk of CIMP tumors and the association was significantly less pronounced for non-CIMP tumors (P-heterogeneity = 0.01). Physical activity was inversely associated with both subtypes, but a dose-response association was observed only for non-CIMP tumors (P-trend = 0.02).

Conclusions

Body size, especially central adiposity, may increase the risk of both CIMP and non-CIMP tumors. Body fat at young age may differentially influence risk. Physical activity appears to decrease the risk of CRC regardless of these molecular subtypes.

About the Authors:

Laura A. E. Hughes

* E-mail: [email protected]

Affiliation: Department of Epidemiology, GROW School for Oncology and Developmental Biology, Maastricht University, Maastricht, The Netherlands

Colinda C. J. M. Simons

Affiliation: Department of Epidemiology, GROW School for Oncology and Developmental Biology, Maastricht University, Maastricht, The Netherlands

Piet A. van den Brandt

Affiliation: Department of Epidemiology, GROW School for Oncology and Developmental Biology, Maastricht University, Maastricht, The Netherlands

R. Alexandra Goldbohm

Affiliation: Department of Prevention and Health, TNO Quality of Life, Leiden, The Netherlands

Anton F. de Goeij

Affiliation: Department of Pathology, GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre, Maastricht, The Netherlands

Adriaan P. de Bruïne

Affiliation: Department of Pathology, GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre, Maastricht, The Netherlands

Manon van Engeland

Affiliation: Department of Pathology, GROW School for Oncology and Developmental Biology, Maastricht University Medical Centre, Maastricht, The Netherlands

Matty P. Weijenberg

Affiliation: Department of Epidemiology, GROW School for Oncology and Developmental Biology, Maastricht University, Maastricht, The Netherlands

Introduction

It is well accepted that indicators of energy balance influence the risk of colorectal cancer (CRC). A high body mass index (BMI), waist circumference, and adult attained height are clear risk factors for CRC, while physical activity has been shown to be protective [1]. Although CRC is one of the best described cancers in terms of genetic and epigenetic events involved [2]–[5], little is known about how measures of anthropometry and physical activity are associated with different molecular subsets of this disease. Elucidating potential differences in risk between molecular subtypes of CRC may lead to a better understanding of CRC, treatment, and prevention. This is especially important as the global prevalence of overweight and obesity continues to rise.

A distinct characteristic of epigenetic instability in CRC is the CpG island methylator phenotype (CIMP), characterized by numerous promoter CpG island hypermethylated tumor suppressor- and DNA repair genes [6]–[9]. This in turn is associated with transcriptional silencing of gene expression [10]. Few studies have investigated associations between indicators of energy balance and CIMP status, and those that have, only considered BMI as an indicator of body size. In a case-control setting, Slattery et al. reported an association between a high BMI and CIMP low but not CIMP high colon tumors [11], and no association between BMI and CIMP status in rectal tumors [12]. Vigorous physical activity was associated with both CIMP high and CIMP low colon tumors, but not rectal tumors [11], [12]. It has been hypothesized that DNA methylation is a consequence of inflammation [13], [14]. Central adiposity is also associated with chronic inflammation [15]. Therefore, considering waist circumference as a risk factor for CIMP in addition to BMI is important. Additionally, methylation is thought to be an early event in CRC [16], so considering height and BMI at a young age may also be informative as these variables are indicative of in utero and early life exposures [17].

Using the prospective Netherlands Cohort Study on diet and cancer (NLCS), we investigated the association between BMI, clothing size (as a proxy for waist circumference) and physical activity and risk of developing a tumor characterized by CIMP. Furthermore, in attempt to elucidate if timing of exposure is important for modulating this risk, we also investigated associations according to BMI at age 20, BMI change, and adult-attained height.

Results

Baseline and molecular characteristics are described in Table 1. Proportionally, CIMP CRC cases were more likely to be female, have a tumor in the proximal colon, and be older than non-CIMP cases.

[Figure omitted. See PDF.]

Table 1. Characteristics of NLCS study participants according to colorectal cancer status after 7.3 years of follow-up.

https://doi.org/10.1371/journal.pone.0018571.t001

Associations between measures of anthropometry, physical activity and CRC risk according to CIMP status are shown in Table 2. After adjustment for age and sex, BMI modeled per 5 kg/m2 increase was statistically significantly associated with CIMP tumors (HR: 1.29, 95% CI: 1.01–1.66), and a statistically significant dose-response relationship was observed when modeling quartiles of BMI (P-trend = 0.02). However, when models were mutually adjusted for trouser/skirt size, these associations were greatly attenuated. BMI was also positively, though not statistically significantly, associated with non-CIMP tumors. These were also attenuated when mutually adjusted for trouser/skirt size. The HRs for CIMP and non-CIMP tumors did not differ significantly from one another (P-heterogeneity = 0.78). BMI at age 20 modeled per 5 kg/m2 increase was associated with increased risk of both subtypes, but HRs did not reach statistical significance. When BMI at age 20 was modeled in quartiles, a statistically significant trend was observed for CIMP tumors (P = 0.03). This trend became borderline significant when the model was adjusted for trouser/skirt size (P-trend  = 0.07). The test for heterogeneity between tumor subtypes was statistically significant with respect to BMI at age 20, even after adjustment for trouser/skirt size (P-heterogeneity = 0.01). No statistically significant observations were observed with respect to BMI change.

[Figure omitted. See PDF.]

Table 2. Hazard Ratios (HR) and 95% Confidence Intervals (CI) for colorectal cancer characterized by CIMP, according to measures of anthropometry in the NLCS after 7.3 years of follow-up.

https://doi.org/10.1371/journal.pone.0018571.t002

Trouser/skirt size, modeled per 2 unit size increase, was associated with both CIMP (HR: 1.20, 95%CI: 1.05–1.37) and non-CIMP tumors (HR: 1.14, 95% CI: 1.05–1.24) after multivariate adjustment, and these associations remained even when models were mutually adjusted for BMI. When trouser/skirt size was considered in approximate quartiles, the test for trend was significant for CIMP (p = 0.02) and borderline significant for non-CIMP (P = 0.06) tumors, although these were attenuated when models were mutually adjusted for BMI. Associations observed for CIMP and non-CIMP tumors were not statistically significantly different from each other (P-heterogeneity = 0.61).

Height, modeled per 5 cm increase was associated with both CIMP (HR: 1.16, 95% CI: 1.01–1.33) and non-CIMP (HR: 1.10, 95% CI: 1.00–1.21) tumors after multivariate adjustment. A borderline statistically significant trend was observed for risk of non-CIMP tumors when height was considered in quartiles (highest versus lowest quartile HR: 1.45, 95% CI: 1.03–2.06; P = 0.06) (p-heterogeneity = 0.98). These associations were attenuated when the model was mutually adjusted for body weight.

With low physical activity as the reference category, there was no inverse dose-response association between physical activity and CIMP tumors, although the HR for intermediate physical activity was statistically significant (HR: 0.50, 95%CI: 0.30–0.82) (Table 3). An inverse association was observed in a dose-response fashion for non-CIMP tumors (intermediate physical activity HR: 0.80, 95%CI: 0.61–1.04; high physical activity HR: 0.69, 95%CI: 0.47–0.96; P-trend = 0.02), however, associations with CIMP and non-CIMP tumors did not differ significantly from each other (P-heterogeneity = 0.33).

[Figure omitted. See PDF.]

Table 3. Hazard Ratios (HR) and 95% Confidence Intervals (CI) for colorectal cancer characterized by CIMP, according to physical activity status in the NLCS after 7.3 years of follow-up.

https://doi.org/10.1371/journal.pone.0018571.t003

Associations for a tumor methylation index in relation to anthropometric risk factors and physical activity are shown in Table 4. There was no clear pattern with respect to the degree of methylation, however, when modeled per 2 unit size increase, trouser/skirt size was associated in a dose-response manner with tumors displaying the highest level of methylation (4–7 genes methylated HR: 1.29, 95%CI: 1.06–1.58; P-trend = 0.08).

[Figure omitted. See PDF.]

Table 4. Hazard ratios (HR) * and 95% confidence intervals (CI) for a methylation index † according to anthropometric risk factors and physical activity in the NLCS 7.3 years of follow-up.

https://doi.org/10.1371/journal.pone.0018571.t004

Discussion

These data suggest that adult body fatness and height may increase the risk of CRC, but are not differentially associated with CIMP and non-CIMP tumors. Contrarily, BMI at age 20 may be a stronger risk factor for CIMP+ tumors. Baseline physical activity appears to decrease the risk of CRC regardless of CIMP status.

A major strength of this study is that we investigated associations in a prospective cohort setting. The NLCS has almost complete ascertainment of colorectal cancer and very little loss to follow-up. Although the number of total CRC after 7.3 years of follow-up in the NLCS was substantial, the number of cases with the CIMP phenotype was small. With limited power to detect associations, it is possible that some findings arose by chance. Another potential limitation of this study is that anthropometric variables were obtained by self-report. However, there are many examples in the literature showing that this method is a valid and reliable tool for assessing body weight and height in cohort studies [18]–[21].

To our knowledge, associations between indicators of energy balance and CIMP status of colorectal tumors have been reported only by Slattery et al. in a case-control setting [11], [12]. In addition to the study design, there are differences between our studies which should be taken into account when comparing results.

The NLCS utilized the Weisenberger panel of genes to define CIMP (CACNA1G, IGF2, NEUROG1, RUNX3 and SOCS1), whereas the study of Slattery et al., utilized the Classic panel (MINT1, MINT2, MINT31, p16 and hMLH1), as well as different cut-offs to define CIMP in tumors [22], [23]. The ‘right’ definition of CIMP is still a highly debated topic, as is the ideal gene panel and the appropriate method of methylation detection [12], [24]. While the Weisenberger panel has been validated, different markers of methylation may be more or less informative with respect to studying different exposures. The prevalence of CIMP in the NLCS population is higher than in the study by Slattery et al. (27% vs. 11%) [12]. However, a difference in primer designs and PCR conditions may substantially change sensitivity and specificity of a particular marker for the detection of CIMP in CRC [25]. Therefore, it is likely that differences in prevalence are not due to the different methods per se, but rather a difference in choice of primers. The MSP analyses that have been conducted in the NLCS have a high detection signal, and subsequently, a higher prevalence of CIMP has been observed. In the present study, we attempted to clarify our observations by constructing a methylation index with different cut-off points that included seven genes commonly methylated in CRC.

In the present study we considered colon and rectal tumors together to increase statistical power. A sensitivity analysis revealed that this did not bias our findings (data not shown). Furthermore, we suggest that idea of combining sub-localizations of tumors may be acceptable when studying molecular endpoints, because this may help explain differences in etiology according to sub-localization.

We observed that BMI was associated with both tumor CIMP and non-CIMP tumors; however, after adjustment for clothing size, these associations disappeared. This is in contrast to case-cohort data suggesting that BMI is associated only with CIMP negative colon tumors and not with rectal tumors [11], [12]. In our study, trouser/skirt size appears to be a strong, independent predictor of both tumor subtypes, even after adjustment for BMI. This is logical, because waist circumference, an indicator of central adiposity, is a stronger predictor of CRC than BMI [1], [26]. When we considered associations according to a methylation index, we did not observe clear associations with respect to BMI and degree of promoter methylation, however, we did observed that trouser/skirt size was associated with the highest level of methylation. That we observed associations between trouser/skirt size and both CIMP and non-CIMP tumors suggests that central adiposity may influence CRC risk through both a methylation and a non-methylation pathway.

Although very few studies have considered associations between BMI and CIMP, a number have considered endpoints in the same pathway as CIMP. Colorectal cancer has distinct molecular subsets, which evolve through different pathways [16]. The path to a serrated adenocarcinoma appears to take one of two major routes: the traditional serrated pathway or the sessile serrated pathway [27]. The sessile serrated pathway is characterized by a high degree of CIMP, BRAF V600E mutations, and ultimately develops into microsatellite instability (MSI) [27]. MSI may serve as a marker for CIMP or other molecular events in CRC [28], therefore, it may be informative to consider the findings of the present study in the context of that research. Two case-control studies have reported that BMI appears associated with microsatellite stable (MSS) tumors, and less with MSI tumors [29], [30]. Neither study reported associations according to waist circumference. In a recent pooled analysis of NLCS and data from the Melbourne Collaborative Cohort Study, we observed similar associations, although the test for heterogeneity between the tumor subtypes was not statistically significant (Hughes et al., submitted).

There is evidence to suggest that early life exposures influence epigenetic mechanisms associated with adult disease risk [31], [32]. Therefore, we also investigated associations between BMI at age 20, BMI change, height and CIMP status. Height is a marker of an aggregated fetal and childhood experience, and can be considered a proxy measure for important nutritional exposures, which affect several hormonal and metabolic axes [17]. In the NLCS, we have observed that childhood and adolescent energy restriction is associated with a decreased risk of CRC later in life [33], [34], which is supported by other population based studies [35]–[37]. We also recently reported that exposure to severe energy restriction during childhood and adolescence was associated with a low risk of developing a CIMP tumor [32]. Furthermore, pooled data suggest that taller individuals are at greater risk of developing a MSI tumor (Hughes et al. submitted). According to the present study, height is not differentially associated with the risk of tumors, however, we did observe significant heterogeneity with respect to BMI at age 20 and tumor subtypes. Taken together, our findings suggest that body size may differentially influence CIMP status during different periods of life, potentially affecting later CRC risk. The association between BMI at age 20 and CIMP tumors was stronger than with non-CIMP tumors, which is in line with previous findings for severe energy restriction during childhood and adolescence. Although our bootstrapping method is quite conservative, we did not observe a clear association with respect to BMI at age 20 and the methylation index and therefore we can not rule out that the differential association with CIMP status is a chance finding. The hypothesis that timing of exposure may influence epigenetic mechanisms requires further investigation.

That we did not observe any clear associations between BMI change and risk of tumors was surprising. This may indicate that metabolic changes in fat may be more important for modulating risk over time, rather than BMI. Alternatively, considering men and women together may have attenuated these observations. Campbell et al. report that adult weight gain was associated with CRC in men, but not in women, and only with respect to individuals who gained more than 21 kg since age 20 [29]. Finally, only considering two time points may not be indicative of true BMI change.

Our findings with respect to physical activity support those of Slattery et al. [11], and suggest that high levels of daily exercise are associated with a decreased risk of both CIMP and non-CIMP tumors. Observations with respect to our methylation index suggest that physical activity may be more protective of tumors with increasing methylation, but more research is required before firm conclusions can be drawn.

Preliminary evidence suggests that molecular markers can be used to classify colorectal cancers into distinct subtypes, which have implications for both etiology and prevention [28]. Fewer tumors arise from the sessile serrated pathway compared to the traditional adenoma pathway [27], [38]. As overweight and obesity are such strong risk factors for CRC, there is a possibility that these conditions may differentially influence risk through pathways and molecular mechanisms other than what we investigated here. More research is needed to clarify the association between indicators of energy balance and epigenetic mechanisms leading to CRC; preferably in a prospective cohort setting, with many cases [39]. Furthermore, as the field of molecular pathological epidemiology [40] continues to evolve, standardizing methods and definitions of molecular endpoints should be addressed. This will become especially critical as more opportunities for pooling data arise.

In conclusion, our findings suggest that measures of anthropometry reflecting a large body size increase the risk of both CIMP and non-CIMP tumors, and that body fat at young age may differentially influence risk. Physical activity appears to decrease the risk of CRC regardless of these molecular subtypes. Our observations reiterate the importance of a healthy body weight with respect to general CRC prevention.

Materials and Methods

Study populations and design

The NLCS is a prospective cohort study that was initiated in 1986 to investigate the association between diet and the development of cancer. It includes 58,279 men and 62,573 women between the ages of 55–69 years at baseline who completed a self-administered food frequency questionnaire involving 150 food items as well as questions on dietary habits, lifestyle, health, and demographics. Municipal registries from throughout the Netherlands were used to constitute an efficient sampling frame. The NLCS uses a case – cohort approach for data processing and analysis; case subjects were derived from the entire cohort, and the number of person-years at risk for the entire cohort was estimated from a sub cohort of 5000 men and women who were randomly sampled from the full cohort at baseline. All sub cohort members who reported prevalent cancer (excluding skin cancer) at baseline were excluded from analyses, leaving 4654. Further details of the NLCS design have been described [41]–[43].

Incident CRC cases were identified by annual record linkage to nine regional cancer registries and a national pathology database (PALGA) [44]. The completeness of cancer follow-up is almost 100% [45]. Paraffin embedded tumor material from CRC patients was retrieved, as described previously [46]. In total, 734 incident CRC patients were identified from a follow-up period of 7.3 years after baseline, excluding the first 2 years of follow-up, of whom a PALGA report of the lesion as well as sufficient DNA was available [46].

The study protocol was approved by the Medical Ethics Committees of the University Hospital Maastricht and TNO Nutrition. On recruitment, participants were informed in writing of the details of the study and its objectives. In accordance with the regulations at that time, written informed consent was obtained when participants returned the completed baseline questionnaire. Tumor material was collected after approval by the ethical review boards of Maastricht University, the National Cancer Registry, and PALGA.

Ascertainment of risk factors

Anthropometric variables.

Height (cm), body weight (kg), and body weight at age 20 (kg) were self-reported on the baseline questionnaire. From these variables, BMI and BMI at age 20, and BMI change were subsequently calculated. At baseline, individuals were also asked to report their lower body (trouser or skirt) clothing size from their clothing label (Dutch sizes). Trouser/skirt size has been shown to be an adequate proxy measure for waist circumference when predicting cancer risk in the NLCS, and details of how clothing size corresponds to waist measurements in men and women in this Dutch population has been published [47]. BMI, BMI at age 20, skirt/trouser size, and height were categorized into approximate sex-specific quartiles. As in previous NLCS analyses, BMI change was categorized as: <0 kg/m2, 0–4 kg/m2, 4–8 kg/m2, and >8 kg/m2 [48].

Physical activity.

With respect to physical activity and CRC risk in the NLCS population, occupational physical activity appears to be more important for men and recreational physical activity for women for predicting risk (Simons et al., submitted). Therefore, we used these two variables to create sex-specific categories of ‘low,’ ‘intermediate’ and ‘high’ physical activity.

Occupational physical activity was derived from data on participants' occupational history. Using information on the type of job and the duration, energy expenditure and sitting time was calculated for the longest and last held job. Energy expenditure was based on a rating system developed by Hettinger et al. [49] and distinguishes between work of low, moderate and high activity which corresponds to an energy expenditure of <8, 8–12 and >12 kJ/min. Men were categorized into the ‘low’ category if their occupational physical activity was <8 kj/minute, ‘intermediate’ if they fell into the 8–12 kJ/minute category, and high if their occupational physical activity was >12 kJ/minute.

Baseline non-occupational physical activity was calculated based on two questions. The first (open-ended) question was ‘How many minutes do you spend on average per day walking or cycling to your work, to go shopping or to take out your dog?’ The reported time spent per day was categorized into ≤10, >10–30, >30–60 and >60 minutes per day. The second question was ‘How many hours of your leisure time do you spend on average per week on 1) recreational cycling, walking, 2) gardening/doing odd jobs and 3) sports, gymnastics?’ Answering possibilities were never, <1 hour per week, 1–2 hours per week and >2 hours per week. The time spent on these activities and the minutes spent per day on walking or cycling to work, to go shopping or to take out the dog were summed to obtain an overall measure of baseline non-occupational physical activity, with categories <30, >30–60, >60–90 and >90 minutes per day. Low physical activity was defined as <30 minutes/day, intermediate as 30–90 minutes/day, and high as >90 minutes/day.

Promoter Methylation Analyses

CIMP in tumor tissue of CRC cases was defined by CpG island promoter hypermethylation of at least 3 out of 5 methylation markers (CACNA1G, IGF2, NEUROG1, RUNX3 and SOCS1), as proposed by Weisenberger et al. [9] were determined by bisulfite modification of 500 ng genomic DNA using a commercially available kit (Zymo Research), and subsequent methylation specific PCR (MSP) [50], [51]. We chose to use MSP as a method because it is effective, specific and does not require specific equipment. It has been shown that results from MSP are in accordance with other technologies, such as MethyLight [52]. Additionally, the methylation status of two other genes, APC and MLH1, were determined and we added them to the CIMP markers to create a methylation index of seven genes.

To facilitate MSP analysis on DNA retrieved from formalin-fixed, paraffin-embedded tissue, DNA was first amplified with flanking PCR primers that amplify bisulfite-modified DNA but do not preferentially amplify methylated or unmethylated DNA. The resulting fragment was used as a template for the MSP reaction. All PCRs were carried out with controls for unmethylated alleles (DNA from normal lymphocytes), methylated alleles [normal lymphocyte DNA treated in vitro with SssI methyltransferase (New England Biolabs, Ipswich, MA)] and a control without DNA. Ten microliters of each MSP reaction was directly loaded on to nondenaturing 6% polyacrylamide gels stained with ethidium bromide and visualised under UV illumination. The MSP analyses were successful for 81%, 79%, 79%, 90%, 83%, 93%, and 93% out of the 734 cases for CACNA1G, IGF2, NEUROG1, RUNX3, SOCS1, MLH1, and APC respectively.

Statistical analyses

Data were analyzed with Stata (version 10, Statacorp, College Station, TX, USA). Cox proportional hazards analysis using the case-cohort approach was used to obtain hazard ratios (HR) and 95% confidence intervals (CI) for the association between measures of anthropometry and physical activity and CRC characterized by CIMP status. To improve statistical power, we considered men and women together. Tests for effect modification by sex were not statistically significant. The proportional hazards assumption was tested using the scaled Schoenfeld residuals and visual inspection of the hazard curves. To account for the additional variance introduced by sampling the subcohort from the entire cohort, standard errors were estimated using the robust option. Statistical significance was tested at the 0.05 level.

For all anthropometric variables in question, three models were considered. The first was adjusted only for age and sex. The second was additionally adjusted for variables identified as being associated with both CRC and energy balance from previous literature. These included family history of CRC (yes/no), smoking status (never smoker, ex-smoker, current smoker), socioeconomic status (level of education: primary school, junior high school, senior high school, higher vocational school, or university), total energy intake (kcal/day), alcohol intake (0, 0.1–4, 5–14, 15–29, ≥30 g/day), physical activity (low, medium, high as previously described), and consumption of red meat, fruit, vegetables, and grains (g/day). Finally, models were mutually adjusted for other anthropometric variables. Models including BMI and BMI at age 20 were mutually adjusted for skirt/trouser size, BMI change was adjusted for BMI at age 20, trouser/skirt size was adjusted for BMI, and height was adjusted for body weight.

We modeled physical activity adjusted for age and sex, and additionally adjusted for trouser/skirt size, family history of CRC (yes/no), smoking status (never smoker, ex-smoker, current smoker), socioeconomic status (level of education: primary school, junior high school, senior high school, higher vocational school, or university), total energy intake (kcal/day), alcohol intake (0, 0.1–4, 5–14, 15–29, ≥30 g/day), clothing size, and consumption of red meat, fruit, vegetables, and fiber (g/day).

To assess how measures of anthropometry and physical activity were associated with the extent of promoter methylation in the CRC tumors, we used the aforementioned methylation index to categorize cases into one of three groups: ‘0–1 genes methylated’, ‘2–3 genes methylated’, or ‘4–7 genes methylated’. Of the 734 cases, 556 had sufficient information to be classified into one of the three categories. Models including anthropometric variables were adjusted for age, sex, and mutually adjusted for other anthropometric variables as previously described, and the model for physical activity was adjusted for age and sex.

Tests for heterogeneity were done to evaluate differences between subtypes of tumors (e.g., CIMP vs. non–CIMP) using the competing risks procedure in STATA. However, the SE for the difference of the log–hazard ratios from this procedure assumes independence of both estimated hazard ratios, which would overestimate that SE and thus overestimate the P values for their difference. Therefore, these P values and the associated confidence intervals were estimated based on a bootstrapping method that was developed for the case-cohort design, as described previously [53]. Each bootstrap analysis was based on 1000 replications.

Acknowledgments

We are indebted to the participants of this study and thank the cancer registries (IKA, IKL, IKMN, IKN, IKO, IKR, IKST, IKW, IKZ, and VIKC) and the Netherlands Nationwide Registry Pathology (PALGA). We also thank S van de Crommert and J Nelissen for their assistance with data entry and data management; A Volovics and A Kester for statistical advice; L Schouten, C de Zwart, M Moll, W van Dijk, M Jansen, and A Pisters for assistance; and H van Montfort, T van Moergastel, L van den Bosch, and R Schmeitz for programming assistance; K Wouters, S Hulsmans, and A Spiertz for laboratory assistance.

Author Contributions

Conceived and designed the experiments: PAvdB RAG APdB AFdG MPW MvE. Performed the experiments: MvE. Analyzed the data: LAEH CCJMS MPW MvE. Contributed reagents/materials/analysis tools: MvE. Wrote the paper: LAEH. Critically reviewed the manuscript: CCJMS PAvdB RAG AFdG APdB MvE MPW.

Citation: Hughes LAE, Simons CCJM, van den Brandt PA, Goldbohm RA, de Goeij AF, de Bruïne AP, et al. (2011) Body Size, Physical Activity and Risk of Colorectal Cancer with or without the CpG Island Methylator Phenotype (CIMP). PLoS ONE6(4): e18571. https://doi.org/10.1371/journal.pone.0018571

References

1. World Cancer Research Fund/American Institute for Cancer Research (2007) (2007) Food, Nutrition, Physical Activity, and the Prevention of Cancer: a Global Perspective. Washington DC: AIRC. World Cancer Research Fund/American Institute for Cancer Research2007Food, Nutrition, Physical Activity, and the Prevention of Cancer: a Global Perspective.Washington DCAIRC2007

2. Derks S, Postma C, Carvalho B, van den Bosch SM, Moerkerk PT, et al. (2008) Integrated analysis of chromosomal, microsatellite and epigenetic instability in colorectal cancer identifies specific associations between promoter methylation of pivotal tumour suppressor and DNA repair genes and specific chromosomal alterations. Carcinogenesis 29: 434–439.S. DerksC. PostmaB. CarvalhoSM van den BoschPT Moerkerk2008Integrated analysis of chromosomal, microsatellite and epigenetic instability in colorectal cancer identifies specific associations between promoter methylation of pivotal tumour suppressor and DNA repair genes and specific chromosomal alterations.Carcinogenesis29434439

3. Issa JP (2004) CpG island methylator phenotype in cancer. Nat Rev Cancer 4: 988–993.JP Issa2004CpG island methylator phenotype in cancer.Nat Rev Cancer4988993

4. Sieber OM, Heinimann K, Tomlinson IP (2003) Genomic instability–the engine of tumorigenesis? Nat Rev Cancer 3: 701–708.OM SieberK. HeinimannIP Tomlinson2003Genomic instability–the engine of tumorigenesis?Nat Rev Cancer3701708

5. Wong JJ, Hawkins NJ, Ward RL (2007) Colorectal cancer: a model for epigenetic tumorigenesis. Gut 56: 140–148.JJ WongNJ HawkinsRL Ward2007Colorectal cancer: a model for epigenetic tumorigenesis.Gut56140148

6. Samowitz WS, Albertsen H, Herrick J, Levin TR, Sweeney C, et al. (2005) Evaluation of a large, population-based sample supports a CpG island methylator phenotype in colon cancer. Gastroenterology 129: 837–845.WS SamowitzH. AlbertsenJ. HerrickTR LevinC. Sweeney2005Evaluation of a large, population-based sample supports a CpG island methylator phenotype in colon cancer.Gastroenterology129837845

7. Shen L, Toyota M, Kondo Y, Lin E, Zhang L, et al. (2007) Integrated genetic and epigenetic analysis identifies three different subclasses of colon cancer. Proc Natl Acad Sci U S A 104: 18654–18659.L. ShenM. ToyotaY. KondoE. LinL. Zhang2007Integrated genetic and epigenetic analysis identifies three different subclasses of colon cancer.Proc Natl Acad Sci U S A1041865418659

8. Toyota M, Ahuja N, Ohe-Toyota M, Herman JG, Baylin SB, et al. (1999) CpG island methylator phenotype in colorectal cancer. Proc Natl Acad Sci U S A 96: 8681–8686.M. ToyotaN. AhujaM. Ohe-ToyotaJG HermanSB Baylin1999CpG island methylator phenotype in colorectal cancer.Proc Natl Acad Sci U S A9686818686

9. Weisenberger DJ, Siegmund KD, Campan M, Young J, Long TI, et al. (2006) CpG island methylator phenotype underlies sporadic microsatellite instability and is tightly associated with BRAF mutation in colorectal cancer. Nat Genet 38: 787–793.DJ WeisenbergerKD SiegmundM. CampanJ. YoungTI Long2006CpG island methylator phenotype underlies sporadic microsatellite instability and is tightly associated with BRAF mutation in colorectal cancer.Nat Genet38787793

10. Ferguson LR, Karunasinghe N, Philpott M (2004) Epigenetic events and protection from colon cancer in New Zealand. Environ Mol Mutagen 44: 36–43.LR FergusonN. KarunasingheM. Philpott2004Epigenetic events and protection from colon cancer in New Zealand.Environ Mol Mutagen443643

11. Slattery ML, Curtin K, Sweeney C, Levin TR, Potter J, et al. (2007) Diet and lifestyle factor associations with CpG island methylator phenotype and BRAF mutations in colon cancer. Int J Cancer 120: 656–663.ML SlatteryK. CurtinC. SweeneyTR LevinJ. Potter2007Diet and lifestyle factor associations with CpG island methylator phenotype and BRAF mutations in colon cancer.Int J Cancer120656663

12. Slattery ML, Curtin K, Wolff RK, Herrick JS, Caan BJ, et al. (2010) Diet, physical activity, and body size associations with rectal tumor mutations and epigenetic changes. Cancer Causes Control 21: 1237–1245.ML SlatteryK. CurtinRK WolffJS HerrickBJ Caan2010Diet, physical activity, and body size associations with rectal tumor mutations and epigenetic changes.Cancer Causes Control2112371245

13. Hodge DR, Peng B, Cherry JC, Hurt EM, Fox SD, et al. (2005) Interleukin 6 supports the maintenance of p53 tumor suppressor gene promoter methylation. Cancer Res 65: 4673–4682.DR HodgeB. PengJC CherryEM HurtSD Fox2005Interleukin 6 supports the maintenance of p53 tumor suppressor gene promoter methylation.Cancer Res6546734682

14. Kang GH, Lee HJ, Hwang KS, Lee S, Kim JH, et al. (2003) Aberrant CpG island hypermethylation of chronic gastritis, in relation to aging, gender, intestinal metaplasia, and chronic inflammation. Am J Pathol 163: 1551–1556.GH KangHJ LeeKS HwangS. LeeJH Kim2003Aberrant CpG island hypermethylation of chronic gastritis, in relation to aging, gender, intestinal metaplasia, and chronic inflammation.Am J Pathol16315511556

15. Monteiro R, Azevedo I (2010) Chronic inflammation in obesity and the metabolic syndrome. Mediators Inflamm 2010. R. MonteiroI. Azevedo2010Chronic inflammation in obesity and the metabolic syndrome.Mediators Inflamm 2010

16. Jass JR (2007) Classification of colorectal cancer based on correlation of clinical, morphological and molecular features. Histopathology 50: 113–130.JR Jass2007Classification of colorectal cancer based on correlation of clinical, morphological and molecular features.Histopathology50113130

17. Okasha M, Gunnell D, Holly J, Davey Smith G (2002) Childhood growth and adult cancer. Best Pract Res Clin Endocrinol Metab 16: 225–241.M. OkashaD. GunnellJ. HollyG. Davey Smith2002Childhood growth and adult cancer.Best Pract Res Clin Endocrinol Metab16225241

18. Bolton-Smith C, Woodward M, Tunstall-Pedoe H, Morrison C (2000) Accuracy of the estimated prevalence of obesity from self reported height and weight in an adult Scottish population. J Epidemiol Community Health 54: 143–148.C. Bolton-SmithM. WoodwardH. Tunstall-PedoeC. Morrison2000Accuracy of the estimated prevalence of obesity from self reported height and weight in an adult Scottish population.J Epidemiol Community Health54143148

19. Nyholm M, Gullberg B, Merlo J, Lundqvist-Persson C, Rastam L, et al. (2007) The validity of obesity based on self-reported weight and height: Implications for population studies. Obesity (Silver Spring) 15: 197–208.M. NyholmB. GullbergJ. MerloC. Lundqvist-PerssonL. Rastam2007The validity of obesity based on self-reported weight and height: Implications for population studies.Obesity (Silver Spring)15197208

20. Spencer EA, Appleby PN, Davey GK, Key TJ (2002) Validity of self-reported height and weight in 4808 EPIC-Oxford participants. Public Health Nutr 5: 561–565.EA SpencerPN ApplebyGK DaveyTJ Key2002Validity of self-reported height and weight in 4808 EPIC-Oxford participants.Public Health Nutr5561565

21. Wada K, Tamakoshi K, Tsunekawa T, Otsuka R, Zhang H, et al. (2005) Validity of self-reported height and weight in a Japanese workplace population. Int J Obes (Lond) 29: 1093–1099.K. WadaK. TamakoshiT. TsunekawaR. OtsukaH. Zhang2005Validity of self-reported height and weight in a Japanese workplace population.Int J Obes (Lond)2910931099

22. Frazier ML, Xi L, Zong J, Viscofsky N, Rashid A, et al. (2003) Association of the CpG island methylator phenotype with family history of cancer in patients with colorectal cancer. Cancer Res 63: 4805–4808.ML FrazierL. XiJ. ZongN. ViscofskyA. Rashid2003Association of the CpG island methylator phenotype with family history of cancer in patients with colorectal cancer.Cancer Res6348054808

23. Park SJ, Rashid A, Lee JH, Kim SG, Hamilton SR, et al. (2003) Frequent CpG island methylation in serrated adenomas of the colorectum. Am J Pathol 162: 815–822.SJ ParkA. RashidJH LeeSG KimSR Hamilton2003Frequent CpG island methylation in serrated adenomas of the colorectum.Am J Pathol162815822

24. Samowitz WS (2007) The CpG island methylator phenotype in colorectal cancer. J Mol Diagn 9: 281–283.WS Samowitz2007The CpG island methylator phenotype in colorectal cancer.J Mol Diagn9281283

25. Ogino S, Kawasaki T, Kirkner GJ, Kraft P, Loda M, et al. (2007) Evaluation of markers for CpG island methylator phenotype (CIMP) in colorectal cancer by a large population-based sample. J Mol Diagn 9: 305–314.S. OginoT. KawasakiGJ KirknerP. KraftM. Loda2007Evaluation of markers for CpG island methylator phenotype (CIMP) in colorectal cancer by a large population-based sample.J Mol Diagn9305314

26. Pischon T, Lahmann PH, Boeing H, Friedenreich C, Norat T, et al. (2006) Body size and risk of colon and rectal cancer in the European Prospective Investigation Into Cancer and Nutrition (EPIC). J Natl Cancer Inst 98: 920–931.T. PischonPH LahmannH. BoeingC. FriedenreichT. Norat2006Body size and risk of colon and rectal cancer in the European Prospective Investigation Into Cancer and Nutrition (EPIC).J Natl Cancer Inst98920931

27. Noffsinger AE (2009) Serrated polyps and colorectal cancer: new pathway to malignancy. Annu Rev Pathol 4: 343–364.AE Noffsinger2009Serrated polyps and colorectal cancer: new pathway to malignancy.Annu Rev Pathol4343364

28. Ogino S, Stampfer M (2010) Lifestyle factors and microsatellite instability in colorectal cancer: the evolving field of molecular pathological epidemiology. J Natl Cancer Inst 102: 365–367.S. OginoM. Stampfer2010Lifestyle factors and microsatellite instability in colorectal cancer: the evolving field of molecular pathological epidemiology.J Natl Cancer Inst102365367

29. Campbell PT, Jacobs ET, Ulrich CM, Figueiredo JC, Poynter JN, et al. (2010) Case-control study of overweight, obesity, and colorectal cancer risk, overall and by tumor microsatellite instability status. JNCI 102: 391–400.PT CampbellET JacobsCM UlrichJC FigueiredoJN Poynter2010Case-control study of overweight, obesity, and colorectal cancer risk, overall and by tumor microsatellite instability status.JNCI102391400

30. Slattery ML, Curtin K, Anderson K, Ma KN, Ballard L, et al. (2000) Associations between cigarette smoking, lifestyle factors, and microsatellite instability in colon tumors. J Natl Cancer Inst 92: 1831–1836.ML SlatteryK. CurtinK. AndersonKN MaL. Ballard2000Associations between cigarette smoking, lifestyle factors, and microsatellite instability in colon tumors.J Natl Cancer Inst9218311836

31. Heijmans BT, Tobi EW, Stein AD, Putter H, Blauw GJ, et al. (2008) Persistent epigenetic differences associated with prenatal exposure to famine in humans. Proc Natl Acad Sci U S A 105: 17046–17049.BT HeijmansEW TobiAD SteinH. PutterGJ Blauw2008Persistent epigenetic differences associated with prenatal exposure to famine in humans.Proc Natl Acad Sci U S A1051704617049

32. Hughes LA, van den Brandt PA, de Bruine AP, Wouters KA, Hulsmans S, et al. (2009) Early life exposure to famine and colorectal cancer risk: a role for epigenetic mechanisms. PLoS One 4: e7951.LA HughesPA van den BrandtAP de BruineKA WoutersS. Hulsmans2009Early life exposure to famine and colorectal cancer risk: a role for epigenetic mechanisms.PLoS One4e7951

33. Dirx MJ, van den Brandt PA, Goldbohm RA, Lumey LH (2003) Energy restriction early in life and colon carcinoma risk: results of The Netherlands Cohort Study after 7.3 years of follow-up. Cancer 97: 46–55.MJ DirxPA van den BrandtRA GoldbohmLH Lumey2003Energy restriction early in life and colon carcinoma risk: results of The Netherlands Cohort Study after 7.3 years of follow-up.Cancer974655

34. Hughes LA, van den Brandt PA, Goldbohm RA, de Goeij AF, de Bruine AP, et al. (2010) Childhood and adolescent energy restriction and subsequent colorectal cancer risk: results from the Netherlands Cohort Study. Int J Epidemiol 39: 1333–1344.LA HughesPA van den BrandtRA GoldbohmAF de GoeijAP de Bruine2010Childhood and adolescent energy restriction and subsequent colorectal cancer risk: results from the Netherlands Cohort Study.Int J Epidemiol3913331344

35. Frankel S, Gunnell DJ, Peters TJ, Maynard M, Davey Smith G (1998) Childhood energy intake and adult mortality from cancer: the Boyd Orr Cohort Study. Bmj 316: 499–504.S. FrankelDJ GunnellTJ PetersM. MaynardG. Davey Smith1998Childhood energy intake and adult mortality from cancer: the Boyd Orr Cohort Study.Bmj316499504

36. Svensson E, Grotmol T, Hoff G, Langmark F, Norstein J, et al. (2002) Trends in colorectal cancer incidence in Norway by gender and anatomic site: an age-period-cohort analysis. Eur J Cancer Prev 11: 489–495.E. SvenssonT. GrotmolG. HoffF. LangmarkJ. Norstein2002Trends in colorectal cancer incidence in Norway by gender and anatomic site: an age-period-cohort analysis.Eur J Cancer Prev11489495

37. Svensson E, Moller B, Tretli S, Barlow L, Engholm G, et al. (2005) Early life events and later risk of colorectal cancer: age-period-cohort modelling in the Nordic countries and Estonia. Cancer Causes Control 16: 215–223.E. SvenssonB. MollerS. TretliL. BarlowG. Engholm2005Early life events and later risk of colorectal cancer: age-period-cohort modelling in the Nordic countries and Estonia.Cancer Causes Control16215223

38. Snover DC (2011) Update on the serrated pathway to colorectal carcinoma. Hum Pathol 42: 1–10.DC Snover2011Update on the serrated pathway to colorectal carcinoma.Hum Pathol42110

39. Ogino S, Cantor M, Kawasaki T, Brahmandam M, Kirkner GJ, et al. (2006) CpG island methylator phenotype (CIMP) of colorectal cancer is best characterised by quantitative DNA methylation analysis and prospective cohort studies. Gut 55: 1000–1006.S. OginoM. CantorT. KawasakiM. BrahmandamGJ Kirkner2006CpG island methylator phenotype (CIMP) of colorectal cancer is best characterised by quantitative DNA methylation analysis and prospective cohort studies.Gut5510001006

40. Ogino S, Chan AT, Fuchs CS, Giovannucci E (2010) Molecular pathological epidemiology of colorectal neoplasia: an emerging transdisciplinary and interdisciplinary field. Gut. S. OginoAT ChanCS FuchsE. Giovannucci2010Molecular pathological epidemiology of colorectal neoplasia: an emerging transdisciplinary and interdisciplinary field.Gut

41. Goldbohm RA, van den Brandt PA, Brants HA, van't Veer P, Al M, et al. (1994) Validation of a dietary questionnaire used in a large-scale prospective cohort study on diet and cancer. Eur J Clin Nutr 48: 253–265.RA GoldbohmPA van den BrandtHA BrantsP. van't VeerM. Al1994Validation of a dietary questionnaire used in a large-scale prospective cohort study on diet and cancer.Eur J Clin Nutr48253265

42. Goldbohm RA, van 't Veer P, van den Brandt PA, van 't Hof MA, Brants HA, et al. (1995) Reproducibility of a food frequency questionnaire and stability of dietary habits determined from five annually repeated measurements. Eur J Clin Nutr 49: 420–429.RA GoldbohmP. van 't VeerPA van den BrandtMA van 't HofHA Brants1995Reproducibility of a food frequency questionnaire and stability of dietary habits determined from five annually repeated measurements.Eur J Clin Nutr49420429

43. van den Brandt PA, Goldbohm RA, van 't Veer P, Volovics A, Hermus RJ, et al. (1990) A large-scale prospective cohort study on diet and cancer in The Netherlands. J Clin Epidemiol 43: 285–295.PA van den BrandtRA GoldbohmP. van 't VeerA. VolovicsRJ Hermus1990A large-scale prospective cohort study on diet and cancer in The Netherlands.J Clin Epidemiol43285295

44. Van den Brandt PA, Schouten LJ, Goldbohm RA, Dorant E, Hunen PM (1990) Development of a record linkage protocol for use in the Dutch Cancer Registry for Epidemiological Research. Int J Epidemiol 19: 553–558.PA Van den BrandtLJ SchoutenRA GoldbohmE. DorantPM Hunen1990Development of a record linkage protocol for use in the Dutch Cancer Registry for Epidemiological Research.Int J Epidemiol19553558

45. Goldbohm RA, van den Brandt PA, Dorant E (1994) Estimation of the coverage of municipalities by cancer registries and PALGA using hospital discharge data. Tijdschr Soc Gezondheidsz 72: 80–84.RA GoldbohmPA van den BrandtE. Dorant1994Estimation of the coverage of municipalities by cancer registries and PALGA using hospital discharge data.Tijdschr Soc Gezondheidsz728084

46. Brink M, de Goeij AF, Weijenberg MP, Roemen GM, Lentjes MH, et al. (2003) K-ras oncogene mutations in sporadic colorectal cancer in The Netherlands Cohort Study. Carcinogenesis 24: 703–710.M. BrinkAF de GoeijMP WeijenbergGM RoemenMH Lentjes2003K-ras oncogene mutations in sporadic colorectal cancer in The Netherlands Cohort Study.Carcinogenesis24703710

47. Hughes LA, Schouten LJ, Goldbohm RA, van den Brandt PA, Weijenberg MP (2009) Self-reported clothing size as a proxy measure for body size. Epidemiology 20: 673–676.LA HughesLJ SchoutenRA GoldbohmPA van den BrandtMP Weijenberg2009Self-reported clothing size as a proxy measure for body size.Epidemiology20673676

48. Pylypchuk RD, Schouten LJ, Goldbohm RA, Schouten HC, van den Brandt PA (2009) Body mass index, height, and risk of lymphatic malignancies: a prospective cohort study. Am J Epidemiol 170: 297–307.RD PylypchukLJ SchoutenRA GoldbohmHC SchoutenPA van den Brandt2009Body mass index, height, and risk of lymphatic malignancies: a prospective cohort study.Am J Epidemiol170297307

49. Hettinger TH, Mueller BH, Gebhard H (1989) Ermittlung des Arbeitsenergieumsatzes bei Dynamisch Muskulaerer Arbeit. Dortmund. TH HettingerBH MuellerH. Gebhard1989Ermittlung des Arbeitsenergieumsatzes bei Dynamisch Muskulaerer Arbeit.Dortmund

50. Derks S, Lentjes MH, Hellebrekers DM, de Bruine AP, Herman JG, et al. (2004) Methylation-specific PCR unraveled. Cell Oncol 26: 291–299.S. DerksMH LentjesDM HellebrekersAP de BruineJG Herman2004Methylation-specific PCR unraveled.Cell Oncol26291299

51. Herman JG, Graff JR, Myohanen S, Nelkin BD, Baylin SB (1996) Methylation-specific PCR: a novel PCR assay for methylation status of CpG islands. Proc Natl Acad Sci U S A 93: 9821–9826.JG HermanJR GraffS. MyohanenBD NelkinSB Baylin1996Methylation-specific PCR: a novel PCR assay for methylation status of CpG islands.Proc Natl Acad Sci U S A9398219826

52. Barault L, Charon-Barra C, Jooste V, de la Vega MF, Martin L, et al. (2008) Hypermethylator phenotype in sporadic colon cancer: study on a population-based series of 582 cases. Cancer Res 68: 8541–8546.L. BaraultC. Charon-BarraV. JoosteMF de la VegaL. Martin2008Hypermethylator phenotype in sporadic colon cancer: study on a population-based series of 582 cases.Cancer Res6885418546

53. de Vogel S, Bongaerts BW, Wouters KA, Kester AD, Schouten LJ, et al. (2008) Associations of dietary methyl donor intake with MLH1 promoter hypermethylation and related molecular phenotypes in sporadic colorectal cancer. Carcinogenesis 29: 1765–1773.S. de VogelBW BongaertsKA WoutersAD KesterLJ Schouten2008Associations of dietary methyl donor intake with MLH1 promoter hypermethylation and related molecular phenotypes in sporadic colorectal cancer.Carcinogenesis2917651773