Huyghe, Jeroen R.

Background: Considering the pathophysiological importance of the intestine-passing microbiota, we hypothesized that the mutational and microbial profiles of colorectal cancer might gradually change along detailed anatomical locations and that these trends might differ by the age of diagnosis and the microsatellite instability (MSI) status. Methods: We assessed somatic mutations and selected bacterial species in colorectal carcinomas along the cecum, ascending colon, transverse colon, descending colon, sigmoid colon, and rectum, utilizing targeted sequencing data ( N = 5685) from the Genetics and Epidemiology of Colorectal Cancer Consortium and tissue microbial data ( N = 1389) from the Health Professionals Follow-up Study and Nurses’ Health Study. Results: There were decreasing trends from the cecum to rectum in the prevalence of gene mutations in the IGF/PI3K, receptor tyrosine kinase (RTK)/RAS, TGFB, and WNT signaling pathways (all P trend < 0.0001), as well as high levels of Bacteroides fragilis ( P trend < 0.0001) and high levels of Fusobacterium nucleatum ( P trend = 0.0029) (but not enterotoxigenic Bacteroides fragilis or pks+ Escherichia coli ). Conversely, TP53 pathway mutations showed an increasing cecum-to-rectum trend ( P trend < 0.0001). Similar to later-onset cases diagnosed at ages 50 years and older, in early-onset cases diagnosed under the age of 50 years, we found generally decreasing cecum-to-rectum trends of IGF/PI3K, RTK/RAS, TGFB, and WNT pathway mutations and an increasing cecum-to-rectum trend in TP53 pathway mutations. By contrast, the cecum-to-rectum prevalence-changing trends in IGF/PI3K, RTK/RAS, TGFB, and TP53 pathway mutations were present in non–MSI-high tumors but not in MSI-high tumors (all P interaction < 0.0052). Conclusion: Non–MSI-high colorectal cancers exhibited decreasing cecum-to-rectum trends of IGF/PI3K, RTK/RAS, and TGFB pathway mutations (and an increasing trend of TP53 pathway mutations), in contrast to MSI-high tumors that demonstrated no such trends. Our study highlights the importance of detailed colorectal subsite information combined with tumor molecular and microbial characteristics in colorectal cancer research.

Abstract Background Physical activity is associated with lower colorectal cancer (CRC) risk, but its association with molecular subtypes defined by genetic and epigenetic alterations of the disease is unclear. Such information may enhance the understanding of the mechanisms related to the benefits of physical activity. Methods Pooled observational (cases: n = 5386; controls: n = 6798; studies n = 5) and genome-wide association data (cases: n = 8178; controls: n = 10 472; studies n = 5) were used. We used multivariable logistic regression models and Mendelian randomization to assess the association between physical activity and the risk of CRC subtypes defined by individual tumor markers (and marker combinations), namely microsatellite instability status, CpG island methylator phenotype status, and BRAF and KRAS mutations. We used case-only analysis to test for differences between molecular subtypes. We applied Bonferroni correction to account for multiple tests. Results In the pooled observational analysis, higher levels of physical activity were associated with lower CRC risk (Obs-per 1SD, odds ratio [OR] = 0.94, 95% confidence interval [CI] = 0.90 to 0.97), with an association that was stronger in males (Obs-per 1SD, OR = 0.91, 95% CI = 0.87 to 0.96) than in females (Obs-per 1SD, OR = 0.97, 95% CI = 0.91 to 1.03; Pinteraction = .04). Higher physical activity was associated with a lower risk of CRC across all molecular subtypes, especially in males. There was no difference in the associations by subtypes by pooled observational or Mendelian randomization analyses. The findings did not differ by study design, anatomical site, and early or late age onset of CRC. Conclusions Our findings suggest that physical activity is not differentially associated with the 4 major molecular subtypes involved in colorectal carcinogenesis, indicating that its benefits extend broadly across colorectal cancer pathogenesis.

Abstract Recent research suggests higher circulating lymphocyte counts may protect against colorectal cancer (CRC). However, the role of specific lymphocyte subtypes and activation states remain unclear. CD4+ T cells—a highly dynamic lymphocyte subtype—undergo gene expression changes upon activation that are critical to their effector function. Previous studies using bulk tissue have limited our understanding of their role in CRC risk to static associations. We applied Mendelian randomization (MR) and genetic colocalisation to evaluate causal relationships of gene expression on CRC risk across multiple CD4+ T cell subtypes and activation states. Genetic proxies were obtained from single-cell transcriptomic data, allowing us to investigate the causal effect of expression of 1,805 genes across CD4+ T cell activation states on CRC risk (78,473 cases; 107,143 controls). Analyses were stratified by CRC anatomical subsites and sex, with sensitivity analyses assessing whether the observed effect estimates were likely to be CD4+ T cell-specific. We identified 6 genes—FADS2, FHL3, HLA-DRB1, HLA-DRB5, RPL28, and TMEM258—with strong evidence for a causal role in CRC development (FDR-P &lt; 0.05; colocalisation H4 &gt; 0.8). Causal estimates varied by CD4+ T cell subtype, activation state, CRC subsite and sex. However, many of genetic proxies used to instrument gene expression in CD4+ T cells also act as eQTLs in other tissues, highlighting the challenges of using genetic proxies to instrument tissue-specific expression changes. We demonstrate the importance of capturing the dynamic nature of CD4+ T cells in understanding CRC risk, and prioritize genes for further investigation in cancer prevention.

Abstract Background Physical activity (PA) is an established protective factor for colorectal cancer (CRC), but it is unclear if genetic variants modify this effect. To investigate this possibility, we conducted a genome-wide gene–PA interaction analysis. Methods Using logistic regression and two-step and joint tests, we analyzed interactions between common genetic variants across the genome and PA in relation to CRC risk. Self-reported PA levels were categorized as active (≥ 8.75 MET-h/wk) vs. inactive (< 8.75 MET-h/wk) and as study- and sex-specific quartiles of activity. Results PA had an overall protective effect on CRC (OR [active vs. inactive] = 0.85; 95%CI = 0.81–0.90). The two-step GxE method identified an interaction between rs4779584, an intergenic variant near the GREM1 and SCG5 genes, and PA for CRC risk (p-interaction = 2.6×10 − 8 ). Stratification by genotype at this locus showed a significant reduction in CRC risk by 20% in active vs. inactive participants with the CC genotype (OR = 0.80; 95%CI = 0.75–0.85), but no significant PA–CRC association among CT or TT carriers. When PA was modeled as quartiles, the 1-d.f. GxE test identified that rs56906466, an intergenic variant near the KCNG1 gene, modified the association between PA and CRC (p-interaction = 3.5×10 − 8 ). Stratification at this locus showed that increase in PA (highest vs. lowest quartile) was associated with a lower CRC risk solely among TT carriers (OR = 0.77; 95%CI = 0.72–0.82). Conclusions In summary, we identified two genetic variants that modified the association between PA and CRC risk. One of them, related to GREM1 and SCG5 , suggests that the bone morphogenetic protein (BMP)-related, inflammatory, and/or insulin signaling pathways may be associated with the protective influence of PA on colorectal carcinogenesis.

Abstract Background: Waist circumference (WC) and its allometric counterpart, “a body shape index” (ABSI), are risk factors for colorectal cancer; however, it is uncertain whether associations with these body measurements are limited to specific molecular subtypes of the disease. Methods: Data from 2,772 colorectal cancer cases and 3,521 controls were pooled from four cohort studies within the Genetics and Epidemiology of Colorectal Cancer Consortium. Four molecular markers (BRAF mutation, KRAS mutation, CpG island methylator phenotype, and microsatellite instability) were analyzed individually and in combination (Jass types). Multivariable logistic and multinomial logistic models were used to assess the associations of WC and ABSI with overall colorectal cancer risk and, in case-only analyses, to evaluate heterogeneity by molecular subtype, respectively. Results: Higher WC (ORper 5 cm = 1.06, 95% confidence interval, 1.04–1.09) and ABSI (ORper 1-SD = 1.07, 95% confidence interval, 1.00–1.14) were associated with elevated colorectal cancer risk. There was no evidence of heterogeneity between the molecular subtypes. No difference was observed regarding the influence of WC and ABSI on the four major molecular markers in proximal colon, distal colon, and rectal cancers, as well as in early- and late-onset colorectal cancers. Associations did not differ in the Jass-type analysis. Conclusions: Higher WC and ABSI were associated with elevated colorectal cancer risk; however, they do not differentially influence all four major molecular mutations involved in colorectal carcinogenesis but underscore the importance of maintaining a healthy body weight in colorectal cancer prevention. Impact: The proposed results have potential utility in colorectal cancer prevention.

Regular, long-term aspirin use may act synergistically with genetic variants, particularly those in mechanistically relevant pathways, to confer a protective effect on colorectal cancer (CRC) risk. We leveraged pooled data from 52 clinical trial, cohort, and case-control studies that included 30,806 CRC cases and 41,861 controls of European ancestry to conduct a genome-wide interaction scan between regular aspirin/nonsteroidal anti-inflammatory drug (NSAID) use and imputed genetic variants. After adjusting for multiple comparisons, we identified statistically significant interactions between regular aspirin/NSAID use and variants in 6q24.1 (top hit rs72833769 ), which has evidence of influencing expression of TBC1D7 (a subunit of the TSC1-TSC2 complex, a key regulator of MTOR activity), and variants in 5p13.1 (top hit rs350047 ), which is associated with expression of PTGER4 (codes a cell surface receptor directly involved in the mode of action of aspirin). Genetic variants with functional impact may modulate the chemopreventive effect of regular aspirin use, and our study identifies putative previously unidentified targets for additional mechanistic interrogation.

Abstract Background: The genotoxin colibactin causes a tumor single-base substitution (SBS) mutational signature, SBS88. It is unknown whether epidemiologic factors’ association with colorectal cancer risk and survival differs by SBS88. Methods: Within the Genetic Epidemiology of Colorectal Cancer Consortium and Colon Cancer Family Registry, we measured SBS88 in 4,308 microsatellite stable/microsatellite instability low tumors. Associations of epidemiologic factors with colorectal cancer risk by SBS88 were assessed using multinomial regression (N = 4,308 cases, 14,192 controls; cohort-only cases N = 1,911), and with colorectal cancer–specific survival using Cox proportional hazards regression (N = 3,465 cases). Results: 392 (9%) tumors were SBS88 positive. Among all cases, the highest quartile of fruit intake was associated with lower risk of SBS88-positive colorectal cancer than SBS88-negative colorectal cancer [odds ratio (OR) = 0.53, 95% confidence interval (CI) 0.37-0.76; OR = 0.75, 95% CI 0.66–0.85, respectively, Pheterogeneity = 0.047]. Among cohort studies, associations of body mass index (BMI), alcohol, and fruit intake with colorectal cancer risk differed by SBS88. BMI ≥30 kg/m2 was associated with worse colorectal cancer–specific survival among those SBS88-positive [hazard ratio (HR) = 3.40, 95% CI 1.47–7.84], but not among those SBS88-negative (HR = 0.97, 95% CI 0.78–1.21, Pheterogeneity = 0.066). Conclusions: Most epidemiologic factors did not differ by SBS88 for colorectal cancer risk or survival. Higher BMI may be associated with worse colorectal cancer–specific survival among those SBS88-positive; however, validation is needed in samples with whole-genome or whole-exome sequencing available. Impact: This study highlights the importance of identification of tumor phenotypes related to colorectal cancer and understanding potential heterogeneity for risk and survival.

Background: Recognizing the early signs of cancer risk is vital for informing prevention, early detection, and survival. Methods: To investigate whether changes in circulating metabolites characterize the early stages of colorectal cancer (CRC) development, we examined the associations between a genetic risk score (GRS) associated with CRC liability (72 single-nucleotide polymorphisms) and 231 circulating metabolites measured by nuclear magnetic resonance spectroscopy in the Avon Longitudinal Study of Parents and Children (N = 6221). Linear regression models were applied to examine the associations between genetic liability to CRC and circulating metabolites measured in the same individuals at age 8 y, 16 y, 18 y, and 25 y. Results: The GRS for CRC was associated with up to 28% of the circulating metabolites at FDR-P < 0.05 across all time points, particularly with higher fatty acids and very-low- and low-density lipoprotein subclass lipids. Two-sample reverse Mendelian randomization (MR) analyses investigating CRC liability (52,775 cases, 45,940 controls) and metabolites measured in a random subset of UK Biobank participants (N = 118,466, median age 58 y) revealed broadly consistent effect estimates with the GRS analysis. In conventional (forward) MR analyses, genetically predicted polyunsaturated fatty acid concentrations were most strongly associated with higher CRC risk. Conclusions: These analyses suggest that higher genetic liability to CRC can cause early alterations in systemic metabolism and suggest that fatty acids may play an important role in CRC development. Funding: This work was supported by the Elizabeth Blackwell Institute for Health Research, University of Bristol, the Wellcome Trust, the Medical Research Council, Diabetes UK, the University of Bristol NIHR Biomedical Research Centre, and Cancer Research UK. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. This work used the computational facilities of the Advanced Computing Research Centre, University of Bristol - http://www.bristol.ac.uk/acrc/ .

Recognizing the early signs of cancer risk is vital for informing prevention, early detection, and survival. To investigate whether changes in circulating metabolites characterise the early stages of colorectal cancer (CRC) development, we examined associations between a genetic risk score (GRS) associated with CRC liability (72 single nucleotide polymorphisms) and 231 circulating metabolites measured by nuclear magnetic resonance spectroscopy in the Avon Longitudinal Study of Parents and Children (N=6,221). Linear regression models were applied to examine associations between genetic liability to colorectal cancer and circulating metabolites measured in the same individuals at age 8, 16, 18 and 25 years. The GRS for CRC was associated with up to 28% of the circulating metabolites at FDR-P<0.05 across all time points, particularly with higher fatty acids and very-low- and low-density lipoprotein subclass lipids. Two-sample reverse Mendelian randomization (MR) analyses investigating CRC liability (52,775 cases, 45,940 controls) and metabolites measured in a random subset of UK Biobank participants (N=118,466, median age 58y) revealed broadly consistent effect estimates with the GRS analysis. In conventional (forward) MR analyses, genetically predicted polyunsaturated fatty acid concentrations were most strongly associated with higher CRC risk. These analyses suggest that higher genetic liability to CRC can cause early alterations in systemic metabolism, and suggest that fatty acids may play an important role in CRC development.