- Research article
- Open Access
Adiposity, metabolites, and colorectal cancer risk: Mendelian randomization study
BMC Medicine volume 18, Article number: 396 (2020)
Higher adiposity increases the risk of colorectal cancer (CRC), but whether this relationship varies by anatomical sub-site or by sex is unclear. Further, the metabolic alterations mediating the effects of adiposity on CRC are not fully understood.
We examined sex- and site-specific associations of adiposity with CRC risk and whether adiposity-associated metabolites explain the associations of adiposity with CRC. Genetic variants from genome-wide association studies of body mass index (BMI) and waist-to-hip ratio (WHR, unadjusted for BMI; N = 806,810), and 123 metabolites from targeted nuclear magnetic resonance metabolomics (N = 24,925), were used as instruments. Sex-combined and sex-specific Mendelian randomization (MR) was conducted for BMI and WHR with CRC risk (58,221 cases and 67,694 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium, Colorectal Cancer Transdisciplinary Study, and Colon Cancer Family Registry). Sex-combined MR was conducted for BMI and WHR with metabolites, for metabolites with CRC, and for BMI and WHR with CRC adjusted for metabolite classes in multivariable models.
In sex-specific MR analyses, higher BMI (per 4.2 kg/m2) was associated with 1.23 (95% confidence interval (CI) = 1.08, 1.38) times higher CRC odds among men (inverse-variance-weighted (IVW) model); among women, higher BMI (per 5.2 kg/m2) was associated with 1.09 (95% CI = 0.97, 1.22) times higher CRC odds. WHR (per 0.07 higher) was more strongly associated with CRC risk among women (IVW OR = 1.25, 95% CI = 1.08, 1.43) than men (IVW OR = 1.05, 95% CI = 0.81, 1.36). BMI or WHR was associated with 104/123 metabolites at false discovery rate-corrected P ≤ 0.05; several metabolites were associated with CRC, but not in directions that were consistent with the mediation of positive adiposity-CRC relations. In multivariable MR analyses, associations of BMI and WHR with CRC were not attenuated following adjustment for representative metabolite classes, e.g., the univariable IVW OR for BMI with CRC was 1.12 (95% CI = 1.00, 1.26), and this became 1.11 (95% CI = 0.99, 1.26) when adjusting for cholesterol in low-density lipoprotein particles.
Our results suggest that higher BMI more greatly raises CRC risk among men, whereas higher WHR more greatly raises CRC risk among women. Adiposity was associated with numerous metabolic alterations, but none of these explained associations between adiposity and CRC. More detailed metabolomic measures are likely needed to clarify the mechanistic pathways.
Colorectal cancer (CRC) is one of the most commonly diagnosed cancers among adults globally [1,2,3]. Obesity is viewed as a likely cause of CRC by the International Agency for Research on Cancer (IARC), the American Institute for Cancer Research (AICR), and the World Cancer Research Fund (WCRF) [3, 4], based largely on positive associations between adiposity and CRC risk from observational epidemiology. Further, the limited data available from observational studies suggest that intentional weight loss lowers the risk of CRC in postmenopausal women . Mendelian randomization (MR) studies, which use genetic variants as instruments (proxies) for adiposity given their randomly allocated and fixed nature , further support causality [7,8,9]. Despite this growing consensus, it remains unclear whether the effect of adiposity on CRC risk differs among men and women, whether the relationship varies by CRC sub-site, and what the underlying biological mechanisms are. These are important to clarify given the ongoing obesity epidemic and difficulties in reducing adiposity itself [10, 11].
Observationally, body mass index (BMI) relates more strongly to CRC risk among men and waist-to-hip ratio (WHR) relates similarly to CRC risk among men and women . However, recent MR studies suggest that higher BMI more greatly raises CRC risk among women, while higher WHR more greatly raises CRC risk among men [7, 8]. Whether these MR estimates are robust is unclear because they were based on relatively small sample sizes, genetic instruments that were not sex-specific, and genetic instruments for WHR that were conditioned on BMI—all potential sources of bias [13,14,15,16,17].
Adiposity alters the systemic metabolism [18,19,20], but evidence for the effects of adiposity-altered metabolites on CRC is scarce. One MR study suggested that total cholesterol raises CRC risk , while others suggested no effect of blood glucose  and mixed support for fatty acids . Overall, the scope of metabolic traits examined has also been narrow. Targeted metabolomics allows deeper phenotyping at a large scale , and its recent integration with genotype data  enables us to examine the associations of metabolites with CRC using MR. Expanded genotype data for CRC is also available , affording a sample size six times larger than used in previous MR studies (58,221 cases, 67,694 controls).
This study has two aims. First, we aimed to better estimate sex-specific effects of adiposity on CRC risk using two-sample MR. We examined associations of BMI and WHR with CRC risk using expanded GWAS data and genetic instruments for exposures that were sex-specific and were not mutually conditioned, to reduce bias [13,14,15,16,17]. Second, we aimed to identify potential metabolic mediators of effects of adiposity on CRC risk using two-step MR (by examining associations of BMI and WHR with metabolites, and of BMI- or WHR-related metabolites with CRC risk) and multivariable MR (by adjusting associations of BMI and WHR with CRC for representative metabolites).
We used two-sample MR to examine the associations (pertaining to estimates of the effect predicted from genetic variants used as instruments) of adiposity with CRC risk, of adiposity with metabolites, of adiposity-associated metabolites with CRC risk, and finally of adiposity with CRC risk adjusted for representative metabolites. In two-sample MR, SNP-exposure and SNP-outcome associations are obtained from different study sources and combined as a ratio to estimate the effects of exposures on outcomes [13, 27]. Our study aims and assumptions are shown in Fig. 1.
We identified SNPs that were independently associated (low linkage disequilibrium (LD), R2 < 0.001) with BMI and WHR (unadjusted for BMI) at P < 5 × 10−8 from a recent large-scale genome-wide association study (GWAS) meta-analysis of 221,863 to 806,810 male and female adults of European ancestry from the Genetic Investigation of ANthropometric Traits (GIANT) consortium and the UK Biobank  (Additional file 1: Table S1). BMI and WHR are expressed in standard deviation (SD) units. For sex-combined analyses of BMI and WHR, 312 and 209 SNPs were used, respectively. For sex-specific analyses of BMI, 185 and 152 SNPs were used for women and men, respectively. For sex-specific analyses of WHR, 153 and 64 SNPs were used for women and men, respectively. The proportion of variance explained in adiposity traits by instruments ranged from 0.3 to 5.04% (these were based on approximations for BMI using the equation described by Shim et al. ), and F-statistics (a formal test of whether variance explained is sufficiently high to avoid weak instrument bias) for adiposity instruments ranged from 75.81 to 124.49 (Additional file 1: Table S2) which indicated instrument strength above the recommended minimum levels .
We identified SNPs that were independently associated (R2 < 0.001 and P < 5 × 10−8) with metabolites from a GWAS of 123 traits from targeted nuclear magnetic resonance (NMR) metabolomics (Additional file 1: Table S1); these included lipoprotein subclass-specific lipids, amino acids, fatty acids, inflammatory glycoproteins, and others . Between 13,476 and 24,925 adults (men and women combined) of European ancestry were included. Metabolic traits are expressed in SD units. The proportion of variance explained in metabolites by instruments ranged from 0.44 to 12.49%, and F-statistics for metabolite instruments ranged from 30.2 to 220.8 (Additional file 1: Table S2) which indicated sufficient instrument strength for univariable analyses.
Colorectal cancer GWAS data
We obtained SNP estimates from the most comprehensive GWAS of CRC to date , including 58,221 cases and 67,694 controls (sexes combined) from 45 studies within 3 consortia: Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO), Colorectal Cancer Transdisciplinary Study (CORECT), and Colon Cancer Family Registry (CCFR). Across these studies, there were 28,207 CRC cases and 22,204 controls among men, and 24,568 CRC cases and 23,736 controls among women. Cases were diagnosed by a physician and recorded overall and by site (colon, proximal colon, distal colon, rectum). Approximately 92% of the participants were White-European (~ 8% were East Asian). Case distributions are outlined in Additional file 1: Table S3; other study characteristics are detailed elsewhere . Ethics were approved by respective institutional review boards.
First, we examined the associations of BMI and WHR with overall and site-specific CRC using SNP estimates from sex-combined GWAS of exposures as well as outcomes. We then examined the associations of BMI and WHR with overall CRC based on SNP estimates from sex-specific GWAS of exposure as well as outcome (sex-specific GWAS were not available and thus not used for site-specific CRC). Summary statistics were harmonized using the harmonise_data function within the TwoSampleMR R package . All GWAS were assumed to be coded on the forward strand, and harmonization was confirmed as consistent using option 2 of the “action” argument. As sensitivity analyses, up to four MR methods were used to generate effect estimates using the TwoSampleMR R package  which make differing pleiotropy assumptions (detailed in Fig. 1 legend) [29, 36, 37]. When only a single SNP was available, the Wald ratio was used . When ≥ 2 SNPs were available, random-effects inverse-variance-weighted (IVW) , MR-Egger , weighted median (WM) , and weighted mode  models were used. Cochrane’s Q-statistic was used to assess the heterogeneity of SNP effects (smaller P values indicating higher heterogeneity and higher potential for directional pleiotropy ). Scatter plots were used to compare MR models, and “leave-one-SNP-out” analyses were used to detect SNP outliers .
Second, we examined associations of BMI and WHR with metabolites using results from sex-combined GWAS for exposures as well as outcomes (sex-specific GWAS were not available for metabolites, and so sex-specific analyses were not conducted) and the MR models described above. Each metabolite (analyzed as an outcome) that was associated with either BMI or WHR based on an IVW model P value ≤ 0.05 following a false discovery rate (FDR) correction (Benjamini-Hochberg method ) was taken forward and examined for association with CRC risk using the IVW model (if ≥ 2 SNPs) or the Wald ratio (if 1 SNP). Multivariable MR  was also used to examine the associations of BMI and WHR with CRC risk, adjusting for single metabolites that were representative of various metabolite classes based on previous network analyses  and that had the highest instrument strength based on the F-statistic (Additional file 1: Table S2). As a positive control, we adjusted BMI for WHR as a covariate (which is expected to attenuate the association of BMI with CRC risk), and likewise, we adjusted WHR for BMI as a covariate with the same expectation. A smaller set of SNPs for BMI and WHR based on earlier GWAS [44, 45] was used for these multivariable models to avoid a relative dilution of metabolite instrument strength given that the number of SNPs for BMI and WHR from expanded GWAS far outnumbered those for metabolites. Conditional F-statistics were calculated for exposures in multivariable models .
In each instance, MR estimates are interpreted as the change in outcome per SD unit change in the exposure. Estimates for metabolite outcomes reflect SD unit change, and estimates for CRC outcomes reflect odds ratios (OR). Statistical analyses were performed using R (version 3.5.2).
Associations of BMI and WHR with CRC risk
In sex-combined analyses (Fig. 2; Additional file 1: Table S4), higher BMI (per 4.8 kg/m2) was associated with a higher risk of overall CRC (IVW OR = 1.16, 95% CI = 1.07, 1.26). The WM estimate was similar, but the MR-Egger and weighted mode estimates were both reduced (e.g., MR-Egger OR = 1.02, 95% CI = 0.84, 1.25). BMI associations were consistent across CRC sites. Associations were directionally consistent for WHR as for BMI but were marginally stronger—e.g., higher WHR (per 0.09 ratio) was associated with 1.28 (95% CI = 1.16, 1.42) times higher odds of CRC in an IVW model (MR-Egger and weighted mode estimates were each positive but of a smaller magnitude with wide intervals spanning the null). WHR associations were more consistent for colon rather than rectal sub-sites. SNP heterogeneity was similarly high for BMI and WHR (P value range across models = 9.54 × 10−10 to 1.97 × 10−8).
In sex-specific IVW models (Fig. 2; Additional file 1: Table S4), higher BMI (per 4.2 kg/m2) was associated with 1.23 (95% CI = 1.08, 1.38) times higher odds of CRC among men and 1.09 (95% CI = 0.97, 1.22) times higher odds of CRC (per 5.2 kg/m2) among women. In a WM model, this BMI estimate was robust among men (OR = 1.22, 95% CI = 1.02, 1.46) but reduced among women (OR = 1.04, 95% CI = 0.86, 1.26). MR-Egger and weighted mode estimates were similarly imprecise among men and women, and SNP heterogeneity was similar for both. In IVW models, higher WHR (per 0.07 ratio) was associated with 1.25 (95% CI = 1.08, 1.43) times higher odds of CRC among women; this estimate was 1.05 (95% CI = 0.81, 1.36) among men (per 0.07 ratio). This pattern was also supported by WM estimates (OR = 1.14, 95% CI = 0.91, 1.42 among women and OR = 0.95, 95% CI = 0.90, 1.29 among men), and by MR-Egger and weighted mode estimates. SNP heterogeneity was similarly high among men and women.
Scatter plots comparing different MR models and results of the “leave-one-SNP-out” analyses are presented in Additional file 2: Figures S1-42.
Associations of BMI and WHR with metabolites
In sex-combined analyses, higher BMI (per 4.8 kg/m2) or WHR (per 0.09 ratio) was associated with 104 metabolites based on FDR-corrected P value ≤ 0.05 in IVW models (Additional file 2: Figures S43-47; Additional file 1: Table S5). Evidence was strong in relation to lipids including total cholesterol and triglycerides in very low-density lipoproteins (VLDL), low-density lipoproteins (LDL), and high-density lipoproteins (HDL)—e.g., 0.23 SD (95% CI = 0.15, 0.31) higher triglycerides in large VLDL from higher BMI. Associations of higher BMI were also strong with lactate, pyruvate, and branched-chain amino acids—e.g., 0.19 SD (95% CI = 0.13, 0.25) higher isoleucine—and with inflammatory glycoproteins (0.28 SD, 95% CI = 0.20, 0.36 higher). Similar patterns were seen for WHR.
Associations of BMI- or WHR-related metabolites with CRC
Of 104 metabolites associated (as outcomes) with BMI or WHR in sex-combined analyses, 100 had SNPs for use in Wald or IVW models. As shown in Additional file 1: Table S1, 321 unique SNPs were used to instrument 100 metabolites (3 metabolites had 1 SNP, 13 metabolites had < 5 SNPs, and 51 metabolites had < 10 SNPs; SNP counts across metabolites ranged from 1 to 26). Lipid traits showed generally weak associations with CRC which were also in directions inconsistent with the mediation of the adiposity-CRC relationship—e.g., lipids in medium HDL were positively associated with CRC, but these had been negatively associated with BMI or WHR (Fig. 3; Additional file 1: Table S6). In contrast, there was more consistent evidence of a positive association of lipids in intermediate-density lipoprotein (IDL), VLDL, and LDL with a risk of distal colon cancer, and these lipids had been positively associated with higher BMI or WHR. For example, higher total lipids in IDL (per SD) were associated with 1.09 (95% CI = 1.02, 1.15) times higher odds of distal colon cancer. Lipids were unassociated with the risk of proximal colon cancer. Fatty acids were unassociated with CRC risk except for higher monounsaturated fatty acid levels which were associated with a lower risk of rectal cancer (IVW OR = 0.85, 95% CI = 0.75, 0.95; Fig. 4). Lactate and pyruvate were inversely associated with CRC at 0.66 (95% CI = 0.42, 1.03) times lower odds and 0.64 (95% CI = 0.52, 0.80) times lower odds, respectively. However, these metabolites were positively associated with BMI, and so directions were inconsistent with the mediation of the adiposity-CRC relationship. Amino acids and glycoprotein acetyls were unassociated with CRC risk.
Associations of BMI and WHR with CRC risk independent of metabolites
The association of BMI with overall CRC was not attenuated following adjustment for various metabolite classes (Fig. 5; Additional file 1: Table S7). The univariable IVW OR for BMI (per 4.77 kg/m2 higher, based on 67 SNPs) in relation to CRC was 1.12 (95% CI = 1.00, 1.26), whereas this IVW OR was 1.14 (95% CI = 1.01, 1.29) adjusting for VLDL lipids and 1.11 (95% CI = 0.99, 1.26) adjusting for IDL and LDL lipids. Attenuation was greater when adjusting the BMI-CRC association for WHR (positive control), at IVW OR = 0.93 (95% CI = 0.78, 1.11). Results for WHR in relation to CRC were directionally consistent as seen for BMI, with a lack of attenuation upon adjustment for metabolite classes.
We aimed to better estimate sex-specific effects of adiposity on CRC risk and to identify potential metabolic mediators of the effects of adiposity on CRC, using two-sample MR methods and expanded sample sizes. Our results, based on genetic instruments for adiposity that were sex-specific and were not mutually conditioned, suggest that higher BMI more greatly raises CRC risk among men, whereas higher WHR more greatly raises CRC risk among women. In sex-combined mediation analyses, adiposity was associated with numerous metabolic alterations, but none of these alterations explained the associations between adiposity and CRC. More detailed metabolomic measures are likely needed to clarify the mechanistic pathways.
Observational [3, 47] and MR [7,8,9] studies have suggested adverse effects of adiposity on CRC risk, but causal evidence has been lacking regarding sex specificity. Previous MR studies suggested stronger effects of BMI on CRC risk among women [7,8,9], which contradicts observational suggestions of stronger effects among men . The genetic regulation of BMI and WHR shows strong sexual dimorphism, thought attributable to the influence of sex hormones, namely estrogen, and it is important to capture these differences in MR estimates [48, 49]. Our new results are based on instruments for BMI and WHR that were sex-specific and a sample size for CRC that was six times larger than used previously which enabled higher power relative to two previous MR studies of BMI, WHR, and CRC risk [7, 8] (Additional file 2: Figure S48). These new results suggest that BMI more greatly raises CRC risk among men—a reversal of previous MR estimates. This new pattern for BMI and CRC (22% higher risk among men per 4.2 kg/m2 and 9% higher risk among women per 5.2 kg/m2) is highly consistent with observational estimates reviewed by IARC (22% higher risk in men and 9% higher risk in women per 5 kg/m2 ). Our results also support a reversal of previous MR estimates for WHR, with risk now appearing higher among women than among men. This is unexpected since BMI and abdominal fat measures correlate highly [50, 51]; however, given that fat storage is more peripheral in women [18, 19], WHR (unadjusted for BMI) may be a better proxy for the extremeness of fat volume among women since fat may be stored more abdominally only when peripheral fat stores are overwhelmed. As a post hoc comparison, we repeated analyses of the main effects of adiposity on CRC using the sex-combined adiposity instruments in relation to split samples of men and women (Additional file 1: Table S8, A) to examine the potential for biased results. These suggest that use of sex-combined instruments for BMI and WHR would lead to the conclusion that both are associated with higher CRC risk in males as well as females, but with still higher risk with BMI among males and with WHR among females, in contrast to previous MR studies [7, 8]. This suggests that discrepancies in the result patterns are most likely due to the differences in the power of the main adiposity-CRC relationship (Additional file 2: Figure S48).
SNP heterogeneity was high for BMI and WHR with CRC, although this was similar between sexes and directions of effect from sensitivity models were consistent, suggesting balanced SNP heterogeneity. One cause of heterogeneity may be pleiotropy in violation of the exclusion restriction criteria (assumption 3, Fig. 1). This is not unexpected due to the large number of SNPs included in the adipose trait instruments and the many underlying biological pathways that explain variation in adiposity. A future approach to minimizing heterogeneity in instrument selection could be to analyze the association between subsets of genetic variants related to specific pathways of BMI and WHR in relation to CRC; this requires more biological knowledge of these genetic variants than currently exists.
Given the difficulty of weight loss  and the ongoing obesity epidemic, it is increasingly important to identify the biological pathways which explain the effect of adiposity on the risk of chronic diseases including CRC . Adipose tissue is highly metabolically active and secretes pro-inflammatory cytokines such as interleukin (IL)-6 and tumor necrosis factor (TNF)-alpha which may promote tumor initiation . Adipose tissue-derived inflammation also promotes insulin resistance in glucose storage tissues that can lead to hyperinsulinemia , and insulin and insulin-like growth factors (IGF) such as IGF-1 have pro-mitogenic and anti-apoptotic effects that are cancer promotive [47, 54,55,56,57,58]. Our current results suggest effects of BMI or WHR on numerous lipids and pre-glycemic traits; however, few of these traits had any strong association with CRC risk, and the few that did were in a direction that was inconsistent with a mediating role in the adiposity-CRC relationship. Results of a series of multivariable MR models, which adjusted for various metabolites considered representative of broader metabolite classes , suggested that associations of BMI and WHR with CRC risk were highly independent of these metabolites. However, this analysis may be limited by weak instrument bias  given that F-statistics for metabolite instruments included in each multivariable MR model were relatively low. Nevertheless, the results of two complementary approaches to mediation (two-step MR and multivariable MR) provide little evidence that the effects of adiposity on CRC risk are mediated by adiposity-related metabolites that are detectable by NMR metabolomics. Future studies could examine metabolites, proteins, hormones, and inflammatory factors that are detectable by other metabolomic and proteomic platforms.
The few traits that did show consistent directions of effect included total lipids in IDL, LDL, and VLDL particles which were raised by BMI and which in turn raised the risk of distal colon cancer specifically (not proximal colon or rectal cancer). If robust, this pattern may reflect differential sensitivity of colon regions to lipid exposure owing to divergent functions (the distal colon functions primarily in the storage of resultant fecal matter whereas the proximal colon functions primarily in water absorption and fecal solidification ), or it may reflect differential detectability through screening (proximal colon tumors tend to be detected in older ages and at more advanced stages ). Colorectal anatomical regions may also have distinct molecular features , e.g., the distal colon may be more susceptible to p53 mutations and chromosomal instability , whereas the proximal colon may be more mucinous and susceptible to microsatellite instability and B-Raf proto-oncogene expression [63, 64]. Several meta-analyses of long-term follow-ups of randomized controlled trials of LDL cholesterol-lowering statin use suggested no strong evidence of a protective effect of statin used on CRC risk [65,66,67]; CRC sub-sites were largely unexamined. One previous MR study suggested an adverse effect of higher LDL cholesterol, and a protective effect of genetically proxied statin use, on overall CRC risk ; again, CRC sub-sites were not examined. Prospective observational evidence for LDL cholesterol and CRC risk is less consistent than for total cholesterol or triglycerides; heterogeneity in meta-analyzed effect estimates is much higher for LDL cholesterol (82.7% based on an I2 statistic) compared with total cholesterol and triglycerides (46.7% and 47.8%, respectively) . Prospective estimates of lipoprotein subclass measures from metabolomic platforms are lacking as these are only recently available at scale.
The limitations of this study include the non-specificity of genetic variants used as instruments for some metabolites which stems from their expectedly correlated nature (e.g., rs1260326, a SNP in GCKR, was included in genetic instruments for 54 metabolites). A total of 321 unique SNPs was used to instrument 100 metabolites, but the number of instruments available for a given metabolite was typically small. This limits causal inference for individual traits but should not prevent the identification of relevant classes of traits (e.g., lipid, amino acid). It should also be stressed that genetic variants used for metabolites may alter the enzyme expression and so serve as instruments for the metabolizing enzyme itself, not factors influenced downstream of that enzyme. Since inference in MR applies to the most proximal trait that the genetic variant relates to , directing inference to specific glycolytic traits as distinct from their downstream consequences like insulin resistance  (a key result of higher fatness and trigger of tumorigenesis ) is difficult and requires stronger genetic instruments alongside mechanistic insights from preclinical studies . Adiposity was measured indirectly using BMI and WHR because these correlate highly with more objectively measured fat indexes [50, 51] and allow much larger GWAS sample sizes than otherwise possible (comparably strong GWAS were unavailable for waist circumference). UK Biobank data are included within GWAS for both the exposure and outcome used for MR estimates of adiposity for CRC risk. Sample overlap in a two-sample MR setting is reported to contribute to weak instrument bias and inflated type one error rates, resulting in MR estimates that are biased towards confounding-prone observational estimates . However, given that the proportion of sample overlap is presently low (< 5%) and estimated F-statistics are relatively high (each > 70 for adiposity traits), we do not expect considerable bias here. As a post hoc comparison, we obtained CRC summary GWAS statistics with UK Biobank excluded and repeated MR analyses of adiposity for CRC risk. Estimates were largely consistent with or without the inclusion of UK Biobank data (Additional file 1: Table S8, B). Our sex-specific MR investigations were confined to effects of adiposity on overall CRC because sex-specific GWAS were unavailable for site-specific CRC and metabolite outcomes. Sex-stratified GWAS of such outcomes would enable these in the future.
Our results based on sex-specific MR instruments and expanded sample sizes suggest that higher BMI more greatly raises CRC risk among men, whereas higher WHR more greatly raises CRC risk among women. In sex-combined mediation analyses, adiposity was associated with numerous metabolic alterations, but none of these alterations explained the associations between adiposity and CRC. More detailed metabolomic measures are likely needed to clarify the mechanistic pathways.
Availability of data and materials
The summary-level GWAS data on outcomes used in this study are available following an application to the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO): https://www.fredhutch.org/en/research/divisions/public-health-sciences-division/research/cancer-prevention/genetics-epidemiology-colorectal-cancer-consortium-gecco.html.
A copy of the code used in this analysis is available at https://github.com/cb12104/adiposity_metabolites_crc.
Sung H, Siegel RL, Rosenberg PS, Jemal A. Emerging cancer trends among young adults in the USA: analysis of a population-based cancer registry. Lancet Public Health. 2019;4(3):E137-E147. https://doi.org/10.1016/S2468-2667(18)30267-6.
Mauri G, Sartore-Bianchi A, Russo AG, Marsoni S, Bardelli A, Siena S. Early-onset colorectal cancer in young individuals. Mol Oncol. 2019;13(2):109–31.
World Cancer Research Fund/American Institute for Cancer Research. Continuous Update Project Expert Report. Diet, nutrition, physical activity and colorectal cancer. 2018.
Lauby-Secretan B, Scoccianti C, Loomis D, Grosse Y, Bianchini F, Straif K. Body fatness and cancer—viewpoint of the IARC Working Group. N Engl J Med. 2016;375(8):794–8.
Luo J, Hendryx M, Manson JE, Figueiredo JC, LeBlanc ES, Barrington W, et al. Intentional weight loss and obesity-related cancer risk. JNCI Cancer Spectrum. 2019;3(4):pkz054.
Davey Smith G, Ebrahim S. ‘Mendelian randomization’: can genetic epidemiology contribute to understanding environmental determinants of disease? Int J Epidemiol. 2003;32(1):1–22.
Thrift AP, Gong J, Peters U, Chang-Claude J, Rudolph A, Slattery ML, et al. Mendelian randomization study of body mass index and colorectal cancer risk. Cancer Epidemiol Biomark Prev. 2015;24(7):1024–31.
Jarvis D, Mitchell JS, Law PJ, Palin K, Tuupanen S, Gylfe A, et al. Mendelian randomisation analysis strongly implicates adiposity with risk of developing colorectal cancer. Br J Cancer. 2016;115(2):266–72.
Gao C, Patel CJ, Michailidou K, Peters U, Gong J, Schildkraut J, et al. Mendelian randomization study of adiposity-related traits and risk of breast, ovarian, prostate, lung and colorectal cancer. Int J Epidemiol. 2016;45(3):896–908.
Gunter MJ, Riboli E. Obesity and gastrointestinal cancers—where do we go from here? Nature Rev Gastroenterol Hepatol. 2018;15(11):651.
Dombrowski SU, Knittle K, Avenell A, Araujo-Soares V, Sniehotta FF. Long term maintenance of weight loss with non-surgical interventions in obese adults: systematic review and meta-analyses of randomised controlled trials. BMJ. 2014;348:g2646.
World Cancer Research Fund/American Institute for Cancer Research. Diet, nutrition, physical activity and colorectal cancer: continuous update project. 2017.
Lawlor DA. Commentary: two-sample Mendelian randomization: opportunities and challenges. Int J Epidemiol. 2016;45(3):908.
Aschard H, Vilhjálmsson BJ, Joshi AD, Price AL, Kraft P. Adjusting for heritable covariates can bias effect estimates in genome-wide association studies. Am J Hum Genet. 2015;96(2):329–39.
Holmes MV, Ala-Korpela M, Davey SG. Mendelian randomization in cardiometabolic disease: challenges in evaluating causality. Nat Rev Cardiol. 2017;14:577–90.
Hartwig FP, Tilling K, Davey-Smith G, Lawlor DA, Borges M-CJB. Bias in two-sample Mendelian randomization by using covariable-adjusted summary associations. bioRxiv. 2019. p. 816363.
Holmes MV, Davey SG. Problems in interpreting and using GWAS of conditional phenotypes illustrated by ‘alcohol GWAS’. Mol Psych. 2019;24(2):167.
Rosen ED, Spiegelman BM. What we talk about when we talk about fat. Cell. 2014;156(1):20–44.
Kahn SE, Hull RL, Utzschneider KM. Mechanisms linking obesity to insulin resistance and type 2 diabetes. Nature. 2006;444(7121):840–6.
Würtz P, Wang Q, Kangas AJ, Richmond RC, Skarp J, Tiainen M, et al. Metabolic signatures of adiposity in young adults: Mendelian randomization analysis and effects of weight change. PLoS Med. 2014;11(12):e1001765.
Rodriguez-Broadbent H, Law PJ, Sud A, Palin K, Tuupanen S, Gylfe A, et al. Mendelian randomisation implicates hyperlipidaemia as a risk factor for colorectal cancer. Int J Cancer. 2017;140(12):2701–8.
Song M, Lu Y, Gunter M, Murphy N, Banbury BL, Ma W, et al. Type 2 diabetes and glycemic traits in relation to colorectal cancer risk: a Mendelian randomization study. AACR; 2018.
May-Wilson S, Sud A, Law PJ, Palin K, Tuupanen S, Gylfe A, et al. Pro-inflammatory fatty acid profile and colorectal cancer risk: a Mendelian randomisation analysis. Eur J Cancer. 2017;84:228–38.
Würtz P, Kangas AJ, Soininen P, Lawlor DA, Davey Smith G, Ala-Korpela M. Quantitative serum NMR metabolomics in large-scale epidemiology: a primer on-omic technology. Am J Epidemiol. 2017:kwx016.
Kettunen J, Demirkan A, Würtz P, Draisma HH, Haller T, Rawal R, et al. Genome-wide study for circulating metabolites identifies 62 loci and reveals novel systemic effects of LPA. Nat Commun. 2016;7:11122.
Huyghe JR, Bien SA, Harrison TA, Kang HM, Chen S, Schmit SL, et al. Discovery of common and rare genetic risk variants for colorectal cancer. Nature Genet. 2019;51(1):76.
Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Gen. 2014;23(R1):R89–98.
Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40(4):304–14.
Hartwig FP, Davey Smith G, Bowden J. Robust inference in summary data Mendelian randomization via the zero modal pleiotropy assumption. Int J Epidemiol. 2017;46(6):1985–98.
Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44(2):512–25.
Sterne JA, Davey SG. Sifting the evidence—what’s wrong with significance tests? BMJ. 2001;322(7280):226–31.
Wasserstein RL, Lazar NA. The ASA’s statement on p-values: context, process, and purpose. Am Statistician. 2016;70(2):129–33.
Pulit SL, Stoneman C, Morris AP, Wood AR, Glastonbury CA, Tyrrell J, et al. Metaanalysis of genome-wide association studies for body fat distribution in 694 649 individuals of European ancestry. Hum Mol Gen. 2019;28(1):166–74.
Shim H, Chasman DI, Smith JD, Mora S, Ridker PM, Nickerson DA, et al. A multivariate genome-wide association analysis of 10 LDL subfractions, and their response to statin treatment, in 1868 Caucasians. PLoS One. 2015;10(4):e0120758.
Haycock PC, Burgess S, Wade KH, Bowden J, Relton C, Davey SG. Best (but oft-forgotten) practices: the design, analysis, and interpretation of Mendelian randomization studies. Am J Clin Nutr. 2016;103(4):965–78.
Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-base platform supports systematic causal inference across the human phenome. eLife. 2018;7:e34408.
Burgess S, Bowden J, Fall T, Ingelsson E, Thompson SG. Sensitivity analyses for robust causal inference from Mendelian randomization analyses with multiple genetic variants. Epidemiology. 2017;28(1):30.
Wald A. The fitting of straight lines if both variables are subject to error. Ann Mathematical Statistics. 1940;11(3):284–300.
Bowden J, Hemani G, Davey Smith GJAjoe. Invited commentary: Detecting individual and global horizontal pleiotropy in Mendelian randomization—a job for the humble heterogeneity statistic? 2018;187(12):2681–5.
Zheng J, Baird D, Borges M-C, Bowden J, Hemani G, Haycock P, et al. Recent developments in Mendelian randomization studies. Curr Epidemiol Rep. 2017;4(4):330–45.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Statistic Soc: Series B (Methodological). 1995;57(1):289–300.
Sanderson E, Davey Smith G, Windmeijer F, Bowden J. An examination of multivariable Mendelian randomization in the single-sample and two-sample summary data settings. Int J Epidemiol. 2018;dyy262:1–15.
Kujala UM, Mäkinen V-P, Heinonen I, Soininen P, Kangas AJ, Leskinen TH, et al. Long-term leisure-time physical activity and serum metabolome. Circulation. 2012:CIRCULATIONAHA. 112.105551.
Locke AE, Kahali B, Berndt SI, Justice AE, Pers TH, Day FR, et al. Genetic studies of body mass index yield new insights for obesity biology. Nature. 2015;518(7538):197–206.
Shungin D, Winkler TW, Croteau-Chonka DC, Ferreira T, Locke AE, Mägi R, et al. New genetic loci link adipose and insulin biology to body fat distribution. Nature. 2015;518(7538):187–96.
Sanderson E, Spiller W, Bowden J. Testing and correcting for weak and pleiotropic instruments in two-sample multivariable Mendelian randomisation. BioRxiv. 2020. https://doi.org/10.1101/2020.04.02.021980.
Murphy N, Jenab M, Gunter MJ. Adiposity and gastrointestinal cancers: epidemiology, mechanisms and future directions. Nat Rev Gastroenterol Hepatol. 2018;15:659–70.
Wells JC. Sexual dimorphism of body composition. Best Pract Res Clin Endocrinol Metab. 2007;21(3):415–30.
Lovejoy J, Champagne C, De Jonge L, Xie H, Smith S. Increased visceral fat and decreased energy expenditure during the menopausal transition. Int J Obes. 2008;32(6):949–58.
Bell JA, Carslake D, O’Keeffe LM, Frysz M, Howe LD, Hamer M, et al. Associations of body mass and fat indexes with cardiometabolic traits. J Am Coll Cardiol. 2018;72(24):3142–54.
Flegal KM, Shepherd JA, Looker AC, Graubard BI, Borrud LG, Ogden CL, et al. Comparisons of percentage body fat, body mass index, waist circumference, and waist-stature ratio in adults. Am J Clin Nutr. 2009;89(2):500–8.
Wei H-J, Zeng R, Lu J-H, Lai W-FT, Chen W-H, Liu H-Y, et al. Adipose-derived stem cells promote tumor initiation and accelerate tumor growth by interleukin-6 production. Oncotarget. 2015;6:7713–26.
Hotamisligil G. Inflammation and metabolic disorders. Nature. 2006;444(7121):860–7.
Rinaldi S, Cleveland R, Norat T, Biessy C, Rohrmann S, Linseisen J, et al. Serum levels of IGF-I, IGFBP-3 and colorectal cancer risk: results from the EPIC cohort, plus a meta-analysis of prospective studies. Int J Cancer. 2010;126:NA-NA.
Tran TT, Naigamwalla D, Oprescu AI, Lam L, McKeown-Eyssen G, Bruce WR, et al. Hyperinsulinemia, but not other factors associated with insulin resistance, acutely enhances colorectal epithelial proliferation in vivo. Endocrinol. 2006;147:1830–7.
Kiunga GA, Raju J, Sabljic N, Bajaj G, Good CK, Bird RP. Elevated insulin receptor protein expression in experimentally induced colonic tumors. Cancer Lett. 2004;211:145–53.
Kaaks R, Toniolo P, Akhmedkhanov A, Lukanova A, Biessy C, Dechaud H, et al. Serum C-peptide, insulin-like growth factor (IGF)-I, IGF-binding proteins, and colorectal cancer risk in women. J Natl Cancer Inst. 2000;92:1592–600.
Murphy N, Carreras-Torres R, Song M, Chan AT, Martin RM, Papadimitriou N, et al. Circulating levels of insulin-like growth factor 1 and insulin-like growth factor binding protein 3 associate with risk of colorectal cancer based on serologic and Mendelian randomization analyses. Gastroenterology. 2019;158(5):1300–312.e20. https://doi.org/10.1053/j.gastro.2019.12.020.
Burgess S, Thompson SG, CRP CHD Genetics Collaboration. Avoiding bias from weak instruments in Mendelian randomization studies. International Journal of Epidemiology 2011;40(3):755–64.
Gonzalez EC, Roetzheim RG, Ferrante JM, Campbell R. Predictors of proximal vs. distal colorectal cancers. Dis Colon Rectum. 2001;44:251–8.
Jacobs ET, Thompson PA, Martínez MaE. Diet, gender, and colorectal neoplasia. J Clin Gastroenterol 2007;41:731–746.
Okubo R, Masuda H, Nemoto N. p53 mutation found to be a significant prognostic indicator in distal colorectal cancer. Oncol Rep. 2001;8(3):509-14.
Pekow J, Meckel K, Dougherty U, Butun F, Mustafi R, Lim J, et al. Tumor suppressors miR-143 and miR-145 and predicted target proteins API5, ERK5, K-RAS, and IRS-1 are differentially expressed in proximal and distal colon. Am J Physiol-Gastrointestinal Liver Physiol. 2015;308:G179–G87.
Missiaglia E, Jacobs B, D’Ario G, Di Narzo AF, Soneson C, Budinska E, et al. Distal and proximal colon cancers differ in terms of molecular, pathological, and clinical features. Ann Oncol. 2014;25:1995–2001.
Dale KM, Coleman CI, Henyan NN, Kluger J, White CM. Statins and cancer risk: a meta-analysis. JAMA. 2006;295(1):74–80.
Liu Y, Tang W, Wang J, Xie L, Li T, He Y, et al. Association between statin use and colorectal cancer risk: a meta-analysis of 42 studies. Cancer Causes Control. 2014;25(2):237–49.
Lytras T, Nikolopoulos G, Bonovas S. Statins and the risk of colorectal cancer: An updated systematic review and meta-analysis of 40 studies. World J Gastroenterol. 2014;20(7):1858–70.
Yao X, Tian Z. Dyslipidemia and colorectal cancer risk: a meta-analysis of prospective studies. Cancer Causes Control. 2015;26(2):257–68.
Lee S, Zhang C, Kilicarslan M, Piening BD, Bjornson E, Hallström BM, et al. Integrated network analysis reveals an association between plasma mannose levels and insulin resistance. Cell Metab. 2016;24(1):172–84.
Lawlor DA, Tilling K, Davey SG. Triangulation in aetiological epidemiology. Int J Epidemiol. 2016;45(6):1866–86.
Burgess S, Davies NM, Thompson SG. Bias due to participant overlap in two-sample Mendelian randomization. Genet Epidemiol. 2016;40(7):597–608.
ASTERISK: We are very grateful to Dr. Bruno Buecher without whom this project would not have existed. We also thank all those who agreed to participate in this study, including the patients and the healthy control persons, as well as all the physicians, technicians, and students.
COLON and NQplus: The authors would like to thank the COLON and NQplus investigators at Wageningen University & Research and the involved clinicians in the participating hospitals.
CCFR: The Colon CFR graciously thanks the generous contributions of their 42,505 study participants, dedication of study staff, and the financial support from the US National Cancer Institute, without which this important registry would not exist.
CORSA: We kindly thank all those who contributed to the screening project Burgenland against CRC. Furthermore, we are grateful to Doris Mejri and Monika Hunjadi for the laboratory assistance.
CPS-II: The authors thank the CPS-II participants and Study Management Group for their invaluable contributions to this research. The authors would also like to acknowledge the contribution to this study from central cancer registries supported through the Centers for Disease Control and Prevention National Program of Cancer Registries, and cancer registries supported by the National Cancer Institute Surveillance Epidemiology and End Results program.
Czech Republic CCS: We are thankful to all clinicians in major hospitals in the Czech Republic, without whom the study would not be practicable. We are also sincerely grateful to all patients participating in this study.
DACHS: We thank all participants and cooperating clinicians, and Ute Handte-Daub, Utz Benscheid, Muhabbet Celik, and Ursula Eilber for the excellent technical assistance.
EDRN: We acknowledge all the following contributors to the development of the resource: University of Pittsburgh School of Medicine, Department of Gastroenterology, Hepatology and Nutrition: Lynda Dzubinski; University of Pittsburgh School of Medicine, Department of Pathology: Michelle Bisceglia; and University of Pittsburgh School of Medicine, Department of Biomedical Informatics.
Harvard cohorts (HPFS, NHS, PHS): The study protocol was approved by the institutional review boards of the Brigham and Women’s Hospital and Harvard T.H. Chan School of Public Health, and those of participating registries as required. We would like to thank the participants and staff of the HPFS, NHS, and PHS for their valuable contributions as well as the following state cancer registries for their help: AL, AZ, AR, CA, CO, CT, DE, FL, GA, ID, IL, IN, IA, KY, LA, ME, MD, MA, MI, NE, NH, NJ, NY, NC, ND, OH, OK, OR, PA, RI, SC, TN, TX, VA, WA, and WY. The authors assume full responsibility for analyses and interpretation of these data.
Kentucky: We would like to acknowledge the staff at the Kentucky Cancer Registry.
LCCS: We acknowledge the contributions of Jennifer Barrett, Robin Waxman, Gillian Smith, and Emma Northwood in conducting this study.
NCCCS I & II: We would like to thank the study participants and the NC Colorectal Cancer Study staff.
PLCO: The authors thank the PLCO Cancer Screening Trial screening center investigators and the staff from Information Management Services Inc. and Westat Inc. Most importantly, we thank the study participants for their contributions that made this study possible.
The SCCFR and the PMH study graciously thanks the generous contributions of their study participants, dedication of study staff, and the financial support from the US National Cancer Institute, without which this important research was not possible. The content of this manuscript does not necessarily reflect the views or policies of the NIH or any of the collaborating centers in the CCFR, nor does mention of trade names, commercial products, or organizations imply endorsement by the US government, any cancer registry, or the CCFR.
SEARCH: We thank the SEARCH team.
SELECT: We thank the research and clinical staff at the sites that participated in the SELECT study, without whom the trial would not have been successful. We are also grateful to the 35,533 dedicated men who participated in SELECT.
Women’s Health Initiative: The authors thank the WHI investigators and staff for their dedication and the study participants for making the program possible. A full listing of WHI investigators can be found at http://www.whi.org/researchers/Documents%20%20Write%20a%20Paper/WHI%20Investigator%20Short%20List.pdf.
This publication is the work of the authors who are guarantors for its contents. CJB is supported by Diabetes UK (17/0005587), the Wellcome Trust (202802/Z/16/Z), and the World Cancer Research Fund (WCRF UK), as part of the World Cancer Research Fund International grant programme (IIG_2019_2009). JAB is supported by the Cancer Research UK (C18281/A19169) and the Elizabeth Blackwell Institute for Health Research, University of Bristol and the Wellcome Trust Institutional Strategic Support Fund (204813/Z/16/Z). NJT is a Wellcome Trust Investigator (202802/Z/16/Z); is the PI of the Avon Longitudinal Study of Parents and Children (MRC & WT 217065/Z/19/Z); is supported by the University of Bristol NIHR Biomedical Research Centre (BRC-1215-20011) and the MRC Integrative Epidemiology Unit (MC_UU_12013/3); and works within the CRUK Integrative Cancer Epidemiology Programme (C18281/A19169). EEV is supported by Diabetes UK (17/0005587) and the World Cancer Research Fund (WCRF UK), as part of the World Cancer Research Fund International grant programme (IIG_2019_2009). GDS works in a unit funded by the UK Medical Research Council (MC_UU_00011/1) and the University of Bristol. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Where authors are identified as personnel of the International Agency for Research on Cancer/World Health Organization, the authors alone are responsible for the views expressed in this article and they do not necessarily represent the decisions, policy, or views of the International Agency for Research on Cancer/World Health Organization.
Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO): National Cancer Institute, National Institutes of Health, US Department of Health and Human Services (U01 CA164930, U01 CA137088, R01 CA059045, R21 CA191312).
ASTERISK: A Hospital Clinical Research Program (PHRC-BRD09/C) from the University Hospital Center of Nantes (CHU de Nantes) and supported by the Regional Council of Pays de la Loire, the Groupement des Entreprises Françaises dans la Lutte contre le Cancer (GEFLUC), and the Association Anne de Bretagne Génétique and the Ligue Régionale Contre le Cancer (LRCC).
The ATBC Study is supported by the Intramural Research Program of the US National Cancer Institute, National Institutes of Health, and by US Public Health Service contract HHSN261201500005C from the National Cancer Institute, Department of Health and Human Services.
CLUE II: This research was funded by the American Institute for Cancer Research and the Maryland Cigarette Restitution Fund at Johns Hopkins, the NCI (U01 CA86308, P30 CA006973 to W.G. Nelson), and the National Institute on Aging (U01 AG18033).
COLO2&3: National Institutes of Health (R01 CA60987).
ColoCare: This work was supported by the National Institutes of Health (grant numbers R01 CA189184 (Li/Ulrich), U01 CA206110 (Ulrich/Li/Siegel/Figueireido/Colditz, 2P30CA015704-40 (Gilliland), R01 CA207371 (Ulrich/Li)), the Matthias Lackas Foundation, the German Consortium for Translational Cancer Research, and the EU TRANSCAN initiative.
The Colon Cancer Family Registry (CCFR, www.coloncfr.org) was supported in part by funding from the National Cancer Institute (NCI), National Institutes of Health (NIH) (award U01 CA167551), and through U01/U24 cooperative agreements from NCI with the following CCFR centers: Australasian (CA074778 and CA097735), Ontario (OFCCR) (CA074783), Seattle (SFCCR) (CA074794 (and R01 CA076366 to PAN)), USC Consortium (CA074799), Mayo Clinic (CA074800), and Hawaii (CA074806). Support for case ascertainment was provided in part from the Surveillance, Epidemiology, and End Results (SEER) Program and the following US state cancer registries: AZ, CO, MN, NC, and NH; and by the Victoria Cancer Registry (Australia) and Ontario Cancer Registry (Canada). Additional funding for the OFCCR/ARCTIC was through award GL201-043 from the Ontario Research Fund (to BWZ), award 112746 from the Canadian Institutes of Health Research (to TJH), through a Cancer Risk Evaluation (CaRE) Program grant from the Canadian Cancer Society (to SG), and through generous support from the Ontario Ministry of Research and Innovation. The SCCFR Illumina HumanCytoSNP array was supported through NCI award R01 CA076366 (to PAN). The CCFR Set-1 (Illumina 1 M/1 M-Duo) and Set-2 (Illumina Omni1-Quad) scans were supported by NIH awards U01 CA122839 and R01 CA143247 (to GC). The CCFR Set-3 (Affymetrix Axiom CORECT Set array) was supported by NIH award U19 CA148107 and R01 CA81488 (to SBG). The CCFR Set-4 (Illumina OncoArray 600 K SNP array) was supported by NIH award U19 CA148107 (to SBG) and by the Center for Inherited Disease Research (CIDR), which is funded by the NIH to the Johns Hopkins University, contract number HHSN268201200008I. Colon Cancer Family Registry (CCFR): The content of this manuscript does not necessarily reflect the views or policies of the NIH or any of the collaborating centers in the CCFR, nor does mention of trade names, commercial products, or organizations imply endorsement by the US government, any cancer registry, or the CCFR.
COLON: The COLON study is sponsored by Wereld Kanker Onderzoek Fonds, including funds from grant 2014/1179 as part of the World Cancer Research Fund International Regular Grant Programme, by Alpe d’Huzes and the Dutch Cancer Society (UM 2012–5653, UW 2013-5927, UW2015-7946), and by TRANSCAN (JTC2012-MetaboCCC, JTC2013-FOCUS). The Nqplus study is sponsored by a ZonMW investment grant (98-10030); by PREVIEW, the project PREVention of diabetes through lifestyle intervention and population studies in Europe and around the World (PREVIEW) project which received funding from the European Union Seventh Framework Programme (FP7/2007–2013) under grant no. 312057; by funds from TI Food and Nutrition (cardiovascular health theme), a public-private partnership on precompetitive research in food and nutrition; and by FOODBALL, the Food Biomarker Alliance, a project from JPI Healthy Diet for a Healthy Life.
Colorectal Cancer Transdisciplinary (CORECT) Study: The CORECT Study was supported by the National Cancer Institute, National Institutes of Health (NCI/NIH), US Department of Health and Human Services (grant numbers U19 CA148107, R01 CA81488, P30 CA014089, R01 CA197350, P01 CA196569, R01 CA201407), and National Institutes of Environmental Health Sciences, National Institutes of Health (grant number T32 ES013678).
CORSA: “Österreichische Nationalbank Jubiläumsfondsprojekt” (12511) and Austrian Research Funding Agency (FFG) grant 829675.
CPS-II: The American Cancer Society funds the creation, maintenance, and updating of the Cancer Prevention Study-II (CPS-II) cohort. This study was conducted with Institutional Review Board approval.
CRCGEN: Colorectal Cancer Genetics & Genomics, Spanish study was supported by Instituto de Salud Carlos III, co-funded by FEDER funds—a way to build Europe—(grants PI14-613 and PI09-1286), Agency for Management of University and Research Grants (AGAUR) of the Catalan Government (grant 2017SGR723), and Junta de Castilla y León (grant LE22A10-2). The sample collection of this work was supported by the Xarxa de Bancs de Tumors de Catalunya sponsored by Pla Director d’Oncología de Catalunya (XBTC), Plataforma Biobancos PT13/0010/0013, and ICOBIOBANC, sponsored by the Catalan Institute of Oncology.
Czech Republic CCS: This work was supported by the Grant Agency of the Czech Republic (grants CZ GA CR: GAP304/10/1286 and 1585) and by the Grant Agency of the Ministry of Health of the Czech Republic (grants AZV 15-27580A, AZV 17-30920A, and NV18/03/00199).
DACHS: This work was supported by the German Research Council (BR 1704/6-1, BR 1704/6-3, BR 1704/6-4, CH 117/1-1, HO 5117/2-1, HE 5998/2-1, KL 2354/3-1, RO 2270/8-1, and BR 1704/17-1); the Interdisciplinary Research Program of the National Center for Tumor Diseases (NCT), Germany; and the German Federal Ministry of Education and Research (01KH0404, 01ER0814, 01ER0815, 01ER1505A, and 01ER1505B).
DALS: National Institutes of Health (R01 CA48998 to M. L. Slattery).
EDRN: This work is funded and supported by the NCI, EDRN Grant (U01 CA 84968-06).
EPIC: The coordination of EPIC is financially supported by the European Commission (DGSANCO) and the International Agency for Research on Cancer. The national cohorts are supported by Danish Cancer Society (Denmark); Ligue Contre le Cancer, Institut Gustave Roussy, Mutuelle Générale de l’Education Nationale, Institut National de la Santé et de la Recherche Médicale (INSERM) (France); German Cancer Aid, German Cancer Research Center (DKFZ), Federal Ministry of Education and Research (BMBF), Deutsche Krebshilfe, Deutsches Krebsforschungszentrum and Federal Ministry of Education and Research (Germany); the Hellenic Health Foundation (Greece); Associazione Italiana per la Ricerca sul Cancro-AIRCItaly and National Research Council (Italy); Dutch Ministry of Public Health, Welfare and Sports (VWS), Netherlands Cancer Registry (NKR), LK Research Funds, Dutch Prevention Funds, Dutch ZON (Zorg Onderzoek Nederland), World Cancer Research Fund (WCRF), Statistics Netherlands (The Netherlands); ERC-2009-AdG 232997 and Nordforsk, Nordic Centre of Excellence programme on Food, Nutrition and Health (Norway); Health Research Fund (FIS), PI13/00061 to Granada, PI13/01162 to EPIC-Murcia, Regional Governments of Andalucía, Asturias, Basque Country, Murcia and Navarra, ISCIII RETIC (RD06/0020) (Spain); Swedish Cancer Society, Swedish Research Council and County Councils of Skåne and Västerbotten (Sweden); Cancer Research UK (14136 to EPIC-Norfolk; C570/A16491 and C8221/A19170 to EPIC-Oxford), Medical Research Council (1000143 to EPIC-Norfolk, MR/M012190/1 to EPICOxford) (UK).
EPICOLON: This work was supported by grants from Fondo de Investigación Sanitaria/FEDER (PI08/0024, PI08/1276, PS09/02368, P111/00219, PI11/00681, PI14/00173, PI14/00230, PI17/00509, 17/00878, Acción Transversal de Cáncer), Xunta de Galicia (PGIDIT07PXIB9101209PR), Ministerio de Economia y Competitividad (SAF07-64873, SAF 2010-19273, SAF2014-54453R), Fundación Científica de la Asociación Española contra el Cáncer (GCB13131592CAST), Beca Grupo de Trabajo “Oncología” AEG (Asociación Española de Gastroenterología), Fundación Privada Olga Torres, FP7 CHIBCHA Consortium, Agència de Gestió d’Ajuts Universitaris i de Recerca (AGAUR, Generalitat de Catalunya, 2014SGR135, 2014SGR255, 2017SGR21, 2017SGR653), Catalan Tumour Bank Network (Pla Director d’Oncologia, Generalitat de Catalunya), PERIS (SLT002/16/00398, Generalitat de Catalunya), CERCA Programme (Generalitat de Catalunya), and COST Action BM1206 and CA17118. CIBERehd is funded by the Instituto de Salud Carlos III.
ESTHER/VERDI. This work was supported by grants from the Baden-Württemberg Ministry of Science, Research and Arts and the German Cancer Aid.
Harvard cohorts (HPFS, NHS, PHS): HPFS is supported by the National Institutes of Health (P01 CA055075, UM1 CA167552, U01 CA167552, R01 CA137178, R01 CA151993, R35 CA197735, K07 CA190673, and P50 CA127003), NHS by the National Institutes of Health (R01 CA137178, P01 CA087969, UM1 CA186107, R01 CA151993, R35 CA197735, K07CA190673, and P50 CA127003), and PHS by the National Institutes of Health (R01 CA042182).
Hawaii Adenoma Study: NCI grants R01 CA72520.
HCES-CRC: The Hwasun Cancer Epidemiology Study–Colon and Rectum Cancer (HCES-CRC; grants from Chonnam National University Hwasun Hospital, HCRI15011-1).
Kentucky: This work was supported by the following grant support: Clinical Investigator Award from Damon Runyon Cancer Research Foundation (CI-8); NCI R01CA136726.
LCCS: The Leeds Colorectal Cancer Study was funded by the Food Standards Agency and Cancer Research UK Programme Award (C588/A19167).
MCCS cohort recruitment was funded by VicHealth and Cancer Council Victoria. The MCCS was further supported by Australian NHMRC grants 509348, 209057, 251553, and 504711 and by infrastructure provided by Cancer Council Victoria. Cases and their vital status were ascertained through the Victorian Cancer Registry (VCR) and the Australian Institute of Health and Welfare (AIHW), including the National Death Index and the Australian Cancer Database.
MEC: National Institutes of Health (R37 CA54281, P01 CA033619, and R01 CA063464).
MECC: This work was supported by the National Institutes of Health, US Department of Health and Human Services (R01 CA81488 to SBG and GR).
MSKCC: The work at Sloan Kettering in New York was supported by the Robert and Kate Niehaus Center for Inherited Cancer Genomics and the Romeo Milio Foundation. Moffitt: This work was supported by funding from the National Institutes of Health (grant numbers R01 CA189184, P30 CA076292), Florida Department of Health Bankhead-Coley Grant 09BN-13, and the University of South Florida Oehler Foundation. Moffitt contributions were supported in part by the Total Cancer Care Initiative, Collaborative Data Services Core, and Tissue Core at the H. Lee Moffitt Cancer Center & Research Institute, a National Cancer Institute-designated Comprehensive Cancer Center (grant number P30 CA076292).
The Multiethnic Cohort Study is supported by the National Cancer Institute grant (CA164973).
NCCCS I & II: We acknowledge the funding support for this project from the National Institutes of Health, R01 CA66635 and P30 DK034987.
NFCCR: This work was supported by an Interdisciplinary Health Research Team award from the Canadian Institutes of Health Research (CRT 43821), the National Institutes of Health, US Department of Health and Human Services (U01 CA74783), and National Cancer Institute of Canada grants (18223 and 18226). The authors wish to acknowledge the contribution of Alexandre Belisle and the genotyping team of the McGill University and Génome Québec Innovation Centre, Montréal, Canada, for genotyping the Sequenom panel in the NFCCR samples. Funding was provided to Michael O. Woods by the Canadian Cancer Society Research Institute.
NSHDS: Swedish Cancer Society; Cancer Research Foundation in Northern Sweden; Swedish Research Council; J C Kempe Memorial Fund; Faculty of Medicine, Umeå University, Umeå, Sweden; and Cutting-Edge Research Grant from the County Council of Västerbotten, Sweden.
OFCCR: The Ontario Familial Colorectal Cancer Registry was supported in part by the National Cancer Institute (NCI) of the National Institutes of Health (NIH) under award U01 CA167551 and award U01/U24 CA074783 (to SG). Additional funding for the OFCCR and ARCTIC testing and genetic analysis was through a Canadian Cancer Society CaRE (Cancer Risk Evaluation) program grant and Ontario Research Fund award GL201-043 (to BWZ), through the Canadian Institutes of Health Research award 112746 (to TJH), and through generous support from the Ontario Ministry of Research and Innovation.
OSUMC: OCCPI funding was provided by Pelotonia, and HNPCC funding was provided by the NCI (CA16058 and CA67941).
PLCO: Intramural Research Program of the Division of Cancer Epidemiology and Genetics and supported by contracts from the Division of Cancer Prevention, National Cancer Institute, NIH, DHHS. Funding was provided by the National Institutes of Health (NIH), Genes, Environment and Health Initiative (GEI) Z01 CP 010200, NIH U01 HG004446, and NIH GEI U01 HG 004438.
SCCFR: The Seattle Colon Cancer Family Registry was supported in part by the National Cancer Institute (NCI) of the National Institutes of Health (NIH) under awards U01 CA167551. Additional support for the SFCCR, Postmenopausal Hormones and Colon Cancer (PMH) study and the SCCFR Illumina HumanCytoSNP array were through NCI/NIH awards U01/U24 CA074794 and R01 CA076366 (to PAN).
SEARCH: The University of Cambridge has received salary support in respect of PDPP from the NHS in the East of England through the Clinical Academic Reserve. Cancer Research UK (C490/A16561); the UK National Institute for Health Research Biomedical Research Centres at the University of Cambridge.
SELECT: Research reported in this publication was supported in part by the National Cancer Institute of the National Institutes of Health under Award Numbers U10 CA37429 (CD Blanke) and UM1 CA182883 (CM Tangen/IM Thompson). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
SMS: This work was supported by the National Cancer Institute (grant P01 CA074184 to J.D.P. and P.A.N.; grants R01 CA097325, R03 CA153323, and K05 CA152715 to P.A.N.) and the National Center for Advancing Translational Sciences at the National Institutes of Health (grant KL2 TR000421 to A.N.B.-H.)
The Swedish Low-risk Colorectal Cancer Study: The study was supported by grants from the Swedish research council: K2015-55X-22674-01-4, K2008-55X-20157-03-3, K2006-72X-20157-01-2, and the Stockholm County Council (ALF project).
Swedish Mammography Cohort and Cohort of Swedish Men: This work is supported by the Swedish Research Council /Infrastructure grant, the Swedish Cancer Foundation, and the Karolinska Institute’s Distinguished Professor Award to Alicja Wolk.
UK Biobank: This research has been conducted using the UK Biobank Resource under Application Number 8614.
VITAL: National Institutes of Health (K05 CA154337).
Women’s Health Initiative: The WHI program is funded by the National Heart, Lung, and Blood Institute; National Institutes of Health; and US Department of Health and Human Services through contracts HHSN268201600018C, HHSN268201600001C, HHSN268201600002C, HHSN268201600003C, and HHSN268201600004C.
Ethics approval and consent to participate
All participants provided written informed consent, and each study was approved by the relevant research ethics committee or institutional review board.
Consent for publication
None to declare.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Genetic variants used to instrument BMI, WHR and metabolites. Table S2. Assesment of instrument strength. Table S3. Colorectal cancer case distributions by study, sex and site. Table S4. LogOR colorectal cancer per SD higher BMI or WHR. Table S5. Beta change in NMR-detected metabolite per SD higher BMI or WHR. Table S6. LogOR colorectal cancer per SD higher BMI or WHR-driven NMR-detected metabolite. Table S7. Risk of overall colorectal cancer per SD higher adipose or metabolite trait, estimated using multivariable Mendelian randomization. Table S8. Posthoc investigations.
Figure S1. Scatter plot of SNP-BMI and SNP-CRC associations. Figure S2. Scatter plot of SNP-BMI and SNP-CRC associations (female specific). Figure S3. Scatter plot of SNP-BMI and SNP-CRC associations (male specific). Figure S4. Scatter plot of SNP-BMI and SNP-colon cancer associations. Figure S5. Scatter plot of SNP-BMI and SNP-proximal colon cancer associations. Figure S6. Scatter plot of SNP-BMI and SNP-distal colon cancer associations. Figure S7. Scatter plot of SNP-BMI and SNP-rectal cancer associations. Figure S8. Forest plot showing individual SNP (black) and combined MR estimates (red; Egger and IVW) for the effect of BMI on CRC. Figure S9. Forest plot showing individual SNP (black) and combined MR estimates (red; Egger and IVW) for the effect of BMI on CRC (female specific). Figure S10. Forest plot showing individual SNP (black) and combined MR estimates (red; Egger and IVW) for the effect of BMI on CRC (male specific). Figure S11. Forest plot showing individual SNP (black) and combined MR estimates (red; Egger and IVW) for the effect of BMI on colon cancer. Figure S12. Forest plot showing individual SNP (black) and combined MR estimates (red; Egger and IVW) for the effect of BMI on proximal colon cancer. Figure S13. Forest plot showing individual SNP (black) and combined MR estimates (red; Egger and IVW) for the effect of BMI on distal colon cancer. Figure S14. Forest plot showing individual SNP (black) and combined MR estimates (red; Egger and IVW) for the effect of BMI on rectal cancer. Figure S15. Leave-one-out plot showing the association between BMI and CRC, following SNP-by-SNP removal from the model. Figure S16. Leave-one-out plot showing the association between BMI and CRC (femalespecific), following SNP-by-SNP removal from the model. Figure S17. Leave-one-out plot showing the association between BMI and CRC (malespecific), following SNP-by-SNP removal from the model. Figure S18. Leave-one-out plot showing the association between BMI and colon cancer, following SNP-by-SNP removal from the model. Figure S19. Leave-one-out plot showing the association between BMI and proximal colon cancer, following SNP-by-SNP removal from the model. Figure S20. Leave-one-out plot showing the association between BMI and distal colon cancer, following SNP-by-SNP removal from the model. Figure S21. Leave-one-out plot showing the association between BMI and rectal cancer, following SNP-by-SNP removal from the model. Figure S22. Scatter plot of SNP-WHR and SNP-CRC associations. Figure S23. Scatter plot of SNP-WHR and SNP-CRC associations (female specific). Figure S24. Scatter plot of SNP-WHR and SNP-CRC associations (male specific). Figure S25. Scatter plot of SNP-WHR and SNP-colon cancer associations. Figure S26. Scatter plot of SNP-WHR and SNP-proximal colon cancer associations. Figure S27. Scatter plot of SNP-WHR and SNP-distal colon cancer associations. Figure S28. Scatter plot of SNP-WHR and SNP-rectal cancer associations. Figure S29. Forest plot showing individual SNP (black) and combined MR estimates (red; Egger and IVW) for the effect of WHR on CRC. Figure S30. Forest plot showing individual SNP (black) and combined MR estimates (red; Egger and IVW) for the effect of WHR on CRC (female specific). Figure S31. Forest plot showing individual SNP (black) and combined MR estimates (red; Egger and IVW) for the effect of WHR on CRC (male specific). Figure S32. Forest plot showing individual SNP (black) and combined MR estimates (red; Egger and IVW) for the effect of WHR on colon cancer. Figure S33. Forest plot showing individual SNP (black) and combined MR estimates (red; Egger and IVW) for the effect of WHR on proximal colon cancer. Figure S34. Forest plot showing individual SNP (black) and combined MR estimates (red; Egger and IVW) for the effect of WHR on distal colon cancer. Figure S35. Forest plot showing individual SNP (black) and combined MR estimates (red; Egger and IVW) for the effect of WHR on rectal cancer. Figure S36. Leave-one-out plot showing the association between WHR and CRC, following SNP-by-SNP removal from the model. Figure S37. Leave-one-out plot showing the association between WHR and CRC, following SNP-by-SNP removal from the model (female specific). Figure S38. Leave-one-out plot showing the association between WHR and CRC, following SNP-by-SNP removal from the model (male specific). Figure S39. Leave-one-out plot showing the association between WHR and colon cancer, following SNP-by-SNP removal from the model. Figure S40. Leave-one-out plot showing the association between WHR and proximal colon cancer, following SNP-by-SNP removal from the model. Figure S41. Leave-one-out plot showing the association between WHR and distal colon cancer, following SNP-by-SNP removal from the model. Figure S42. Leave-one-out plot showing the association between WHR and rectal cancer, following SNP-by-SNP removal from the model. Figure S43. Effects of BMI and WHR on circulating metabolite levels (NMR-detected metabolites, 1 of 5), based on two-sample MR (IVW models) in summary GWAS consortia data. Figure S44. Effects of BMI and WHR on circulating metabolite levels (NMR-detected metabolites, 2 of 5), based on two-sample MR (IVW models) in summary GWAS consortia data. Figure S45. Effects of BMI and WHR on circulating metabolite levels (NMR-detected metabolites, 3 of 5), based on two-sample MR (IVW models) in summary GWAS consortia data. Figure S46. Effects of BMI and WHR on circulating metabolite levels (NMR-detected metabolites, 4 of 5), based on two-sample MR (IVW models) in summary GWAS consortia data. Figure S47. Effects of BMI and WHR on circulating metabolite levels (NMR-detected metabolites, 5 of 5), based on two-sample MR (IVW models) in summary GWAS consortia data. Figure S48. Power curves for MR analyses, based on samples sizes for colorectal cancer in the present study (black), Thrift et al., 2015 (blue) and Jarvis et al., 2016 (purple). Upper and lower power curves describe genetic instruments explaining 5% and 0.3% of variance respectively for each study.
About this article
Cite this article
Bull, C.J., Bell, J.A., Murphy, N. et al. Adiposity, metabolites, and colorectal cancer risk: Mendelian randomization study. BMC Med 18, 396 (2020). https://doi.org/10.1186/s12916-020-01855-9
- Body mass index
- Waist-to-hip ratio
- Colorectal cancer
- Mendelian randomization