Woods, Michael O.

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 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 Colorectal cancer risk can be impacted by genetic, environmental, and lifestyle factors, including diet and obesity. Gene-environment interactions (G × E) can provide biological insights into the effects of obesity on colorectal cancer risk. Here, we assessed potential genome-wide G × E interactions between body mass index (BMI) and common SNPs for colorectal cancer risk using data from 36,415 colorectal cancer cases and 48,451 controls from three international colorectal cancer consortia (CCFR, CORECT, and GECCO). The G × E tests included the conventional logistic regression using multiplicative terms (one degree of freedom, 1DF test), the two-step EDGE method, and the joint 3DF test, each of which is powerful for detecting G × E interactions under specific conditions. BMI was associated with higher colorectal cancer risk. The two-step approach revealed a statistically significant G×BMI interaction located within the Formin 1/Gremlin 1 (FMN1/GREM1) gene region (rs58349661). This SNP was also identified by the 3DF test, with a suggestive statistical significance in the 1DF test. Among participants with the CC genotype of rs58349661, overweight and obesity categories were associated with higher colorectal cancer risk, whereas null associations were observed across BMI categories in those with the TT genotype. Using data from three large international consortia, this study discovered a locus in the FMN1/GREM1 gene region that interacts with BMI on the association with colorectal cancer risk. Further studies should examine the potential mechanisms through which this locus modifies the etiologic link between obesity and colorectal cancer. Significance: This gene-environment interaction analysis revealed a genetic locus in FMN1/GREM1 that interacts with body mass index in colorectal cancer risk, suggesting potential implications for precision prevention strategies.

INTRODUCTION: Early-onset colorectal cancer diagnosed before the age of 50 years has been increasing. Likely reflecting the pathogenic role of the intestinal microbiome, which gradually changes across the entire colorectal length, the prevalence of certain tumor molecular characteristics gradually changes along colorectal subsites. Understanding how colorectal tumor molecular features differ by age and tumor location is important in personalized patient management. METHODS: Using 14,004 cases with colorectal cancer including 3,089 early-onset cases, we examined microsatellite instability (MSI), CpG island methylator phenotype (CIMP), and KRAS and BRAF mutations in carcinomas of the cecum, ascending colon, transverse colon, descending colon, sigmoid colon, and rectum and compared early-onset cases with later-onset cases. RESULTS: The proportions of MSI-high, CIMP-high, and BRAF -mutated early-onset tumors were lowest in the rectum (8.8%, 3.4%, and 3.5%, respectively) and highest in the ascending colon (46% MSI-high; 15% CIMP-high) or transverse colon (8.6% BRAF -mutated) (all P trend <0.001 across the rectum to ascending colon). Compared with later-onset tumors, early-onset tumors showed a higher prevalence of MSI-high status and a lower prevalence of CIMP-high status and BRAF mutations in most subsites. KRAS mutation prevalence was higher in the cecum compared with that in the other subsites in both early-onset and later-onset tumors ( P < 0.001). Notably, later-onset MSI-high tumors showed a continuous decrease in KRAS mutation prevalence from the rectum (36%) to ascending colon (9%; P trend <0.001), followed by an increase in the cecum (14%), while early-onset MSI-high cancers showed no such trend. DISCUSSION: Our findings support biogeographical and pathogenic heterogeneity of colorectal carcinomas in different colorectal subsites and age groups.

Abstract Background: Tobacco smoking is an established risk factor for colorectal cancer. However, genetically defined population subgroups may have increased susceptibility to smoking-related effects on colorectal cancer. Methods: A genome-wide interaction scan was performed including 33,756 colorectal cancer cases and 44,346 controls from three genetic consortia. Results: Evidence of an interaction was observed between smoking status (ever vs. never smokers) and a locus on 3p12.1 (rs9880919, P = 4.58 × 10−8), with higher associated risk in subjects carrying the GG genotype [OR, 1.25; 95% confidence interval (CI), 1.20–1.30] compared with the other genotypes (OR &lt;1.17 for GA and AA). Among ever smokers, we observed interactions between smoking intensity (increase in 10 cigarettes smoked per day) and two loci on 6p21.33 (rs4151657, P = 1.72 × 10−8) and 8q24.23 (rs7005722, P = 2.88 × 10−8). Subjects carrying the rs4151657 TT genotype showed higher risk (OR, 1.12; 95% CI, 1.09–1.16) compared with the other genotypes (OR &lt;1.06 for TC and CC). Similarly, higher risk was observed among subjects carrying the rs7005722 AA genotype (OR, 1.17; 95% CI, 1.07–1.28) compared with the other genotypes (OR &lt;1.13 for AC and CC). Functional annotation revealed that SNPs in 3p12.1 and 6p21.33 loci were located in regulatory regions, and were associated with expression levels of nearby genes. Genetic models predicting gene expression revealed that smoking parameters were associated with lower colorectal cancer risk with higher expression levels of CADM2 (3p12.1) and ATF6B (6p21.33). Conclusions: Our study identified novel genetic loci that may modulate the risk for colorectal cancer of smoking status and intensity, linked to tumor suppression and immune response. Impact: These findings can guide potential prevention treatments.

Abstract Background: Evidence suggests that, compared to later-onset colorectal cancer (CRC), early-onset CRC (diagnosed before 50 years of age) occur more frequently in rectal location and less frequently in proximal colon location. Studies also indicate heterogeneity of molecular characteristics within early-onset CRCs and between early-onset and later-onset CRCs. Understanding how tumor molecular features and biological aggressiveness of early-onset CRC differ by detailed tumor location is important in personalized patient management. Methods: Using 14,004 CRC cases including 3,089 early-onset cases from The Cancer Genome Atlas and the Genetics and Epidemiology of CRC Consortium (GECCO), we examined microsatellite instability (MSI), CpG island methylator phenotype (CIMP), and KRAS and BRAF mutations in cancers of the cecum, ascending colon, transverse colon, descending colon, sigmoid colon, and rectum. We further conducted CRC-specific survival analyses of detailed tumor location stratified by molecular characteristics using multivariable Cox proportional hazards models. Results: The proportions of MSI-high, CIMP-high, and BRAF-mutated early-onset tumors gradually increased from the rectum (8.8%, 3.4%, and 3.5%, respectively) to the ascending colon (46% MSI-high; 15% CIMP-high) or transverse colon (8.6% BRAF-mutated), followed by declines in the cecum for MSI-high (36%) and BRAF-mutated (4.6%) tumors. Notably, later-onset MSI-high tumors showed a continuous decrease in KRAS mutation prevalence from the rectum (36%) to ascending colon (9%; Ptrend &lt;0.001) followed by an increase in the cecum (14%), while early-onset MSI-high cancer showed no such trend. In survival analyses, there was a trend for better CRC-specific survival from the cecum to sigmoid colon (Ptrend&lt;0.001) (but not to rectum), which did not significantly differ by age of diagnosis (Pinteraction=0.61). Multivariable-adjusted hazard ratios (HRs) for sigmoid vs. cecum were 0.86 (95% CI, 0.64-1.14) for early-onset cases and 0.81 (95% CI, 0.70-0.93) for later-onset cases. Notably, the prognostic association of tumor location differed by MSI status (Pinteraction=0.001). Non-MSI-high tumors exhibited the cecum-to-sigmoid trend for better CRC-specific survival [Ptrend &lt;0.001; multivariable-adjusted HR for the sigmoid colon (vs. cecum), 0.80; 95% CI, 0.70-0.92], whereas MSI-high tumors demonstrated an opposite cecum-to-sigmoid trend for worse survival (Ptrend = 0.020; the corresponding HR, 2.13; 95% CI, 1.15-3.92). Conclusions: Our findings support biogeographical and pathogenic heterogeneity of CRC in different age groups and colorectal subsites. Both detailed colorectal location and tumor molecular features need to be accounted for prognostication of early-onset and later-onset CRCs to advance precision medicine. Citation Format: Tomotaka Ugai, Yin Cao, Koichiro Haruki, Tabitha A. Harrison, Daniel D. Buchanan, Peter T. Campbell, Andrew T. Chan, Steven Gallinger, Marc J. Gunter, Michael Hoffmeister, Mark A. Jenkins, Roger L. Milne, Polly A. Newcomb, Rish K. Pai, Andrew J. Pellatt, Robert E. Schoen, Bethany Van Guelpen, Michael O. Woods, Amanda I. Phipps, Ulrike Peters, Shuji Ogino. Molecular characteristics and prognosis of early-onset colorectal cancer according to detailed anatomical locations: Comparison to later-onset cases [abstract]. In: Proceedings of the AACR Special Conference on Colorectal Cancer; 2022 Oct 1-4; Portland, OR. Philadelphia (PA): AACR; Cancer Res 2022;82(23 Suppl_1):Abstract nr A010.

Abstract Background Obesity is an established risk factor for colorectal cancer (CRC), but the evidence for the association is inconsistent across molecular subtypes of the disease. Methods We pooled data on body mass index (BMI), tumor microsatellite instability status, CpG island methylator phenotype status, BRAF and KRAS mutations, and Jass classification types for 11 872 CRC cases and 11 013 controls from 11 observational studies. We used multinomial logistic regression to estimate odds ratios (OR) and 95% confidence intervals (CI) adjusted for covariables. Results Higher BMI was associated with increased CRC risk (OR per 5 kg/m2 = 1.18, 95% CI = 1.15 to 1.22). The positive association was stronger for men than women but similar across tumor subtypes defined by individual molecular markers. In analyses by Jass type, higher BMI was associated with elevated CRC risk for types 1-4 cases but not for type 5 CRC cases (considered familial-like/Lynch syndrome microsatellite instability-H, CpG island methylator phenotype-low or negative, BRAF-wild type, KRAS-wild type, OR = 1.04, 95% CI = 0.90 to 1.20). This pattern of associations for BMI and Jass types was consistent by sex and design of contributing studies (cohort or case-control). Conclusions In contrast to previous reports with fewer study participants, we found limited evidence of heterogeneity for the association between BMI and CRC risk according to molecular subtype, suggesting that obesity influences nearly all major pathways involved in colorectal carcinogenesis. The null association observed for the Jass type 5 suggests that BMI is not a risk factor for the development of CRC for individuals with Lynch syndrome.

Abstract Background The use of menopausal hormone therapy (MHT) may interact with genetic variants to influence colorectal cancer (CRC) risk. Methods We conducted a genome-wide, gene-environment interaction between single nucleotide polymorphisms and the use of any MHT, estrogen only, and combined estrogen-progestogen therapy with CRC risk, among 28 486 postmenopausal women (11 519 CRC patients and 16 967 participants without CRC) from 38 studies, using logistic regression, 2-step method, and 2– or 3–degree-of-freedom joint test. A set-based score test was applied for rare genetic variants. Results The use of any MHT, estrogen only and estrogen-progestogen were associated with a reduced CRC risk (odds ratio [OR] = 0.71, 95% confidence interval [CI] = 0.64 to 0.78; OR = 0.65, 95% CI = 0.53 to 0.79; and OR = 0.73, 95% CI = 0.59 to 0.90, respectively). The 2-step method identified a statistically significant interaction between a GRIN2B variant rs117868593 and MHT use, whereby MHT-associated CRC risk was statistically significantly reduced in women with the GG genotype (OR = 0.68, 95% CI = 0.64 to 0.72) but not within strata of GC or CC genotypes. A statistically significant interaction between a DCBLD1 intronic variant at 6q22.1 (rs10782186) and MHT use was identified by the 2–degree-of-freedom joint test. The MHT-associated CRC risk was reduced with increasing number of rs10782186-C alleles, showing odds ratios of 0.78 (95% CI = 0.70 to 0.87) for TT, 0.68 (95% CI = 0.63 to 0.73) for TC, and 0.66 (95% CI = 0.60 to 0.74) for CC genotypes. In addition, 5 genes in rare variant analysis showed suggestive interactions with MHT (2-sided P &lt; 1.2 × 10−4). Conclusion Genetic variants that modify the association between MHT and CRC risk were identified, offering new insights into pathways of CRC carcinogenesis and potential mechanisms involved.

Abstract Background: Currently known associations between common genetic variants and colorectal cancer explain less than half of its heritability of 25%. As alcohol consumption has a J-shape association with colorectal cancer risk, nondrinking and heavy drinking are both risk factors for colorectal cancer. Methods: Individual-level data was pooled from the Colon Cancer Family Registry, Colorectal Transdisciplinary Study, and Genetics and Epidemiology of Colorectal Cancer Consortium to compare nondrinkers (≤1 g/day) and heavy drinkers (&gt;28 g/day) with light-to-moderate drinkers (1–28 g/day) in GxE analyses. To improve power, we implemented joint 2df and 3df tests and a novel two-step method that modifies the weighted hypothesis testing framework. We prioritized putative causal variants by predicting allelic effects using support vector machine models. Results: For nondrinking as compared with light-to-moderate drinking, the hybrid two-step approach identified 13 significant SNPs with pairwise r2 &gt; 0.9 in the 10q24.2/COX15 region. When stratified by alcohol intake, the A allele of lead SNP rs2300985 has a dose–response increase in risk of colorectal cancer as compared with the G allele in light-to-moderate drinkers [OR for GA genotype = 1.11; 95% confidence interval (CI), 1.06–1.17; OR for AA genotype = 1.22; 95% CI, 1.14–1.31], but not in nondrinkers or heavy drinkers. Among the correlated candidate SNPs in the 10q24.2/COX15 region, rs1318920 was predicted to disrupt an HNF4 transcription factor binding motif. Conclusions: Our study suggests that the association with colorectal cancer in 10q24.2/COX15 observed in genome-wide association study is strongest in nondrinkers. We also identified rs1318920 as the putative causal regulatory variant for the region. Impact: The study identifies multifaceted evidence of a possible functional effect for rs1318920.

AbstractDiabetes is an established risk factor for colorectal cancer. However, colorectal cancer is a heterogeneous disease and it is not well understood whether diabetes is more strongly associated with some tumor molecular subtypes than others. A better understanding of the association between diabetes and colorectal cancer according to molecular subtypes could provide important insights into the biology of this association. We used data on lifestyle and clinical characteristics from the Colorectal Cancer Family Registry (CCFR) and the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO), including 9756 colorectal cancer cases (with tumor marker data) and 9985 controls, to evaluate associations between reported diabetes and risk of colorectal cancer according to molecular subtypes. Tumor markers included BRAF and KRAS mutations, microsatellite instability and CpG island methylator phenotype. In the multinomial logistic regression model, comparing colorectal cancer cases to cancer‐free controls, diabetes was positively associated with colorectal cancer regardless of subtype. The highest OR estimate was found for BRAF‐mutated colorectal cancer, n = 1086 (ORfully adj: 1.67, 95% confidence intervals [CI]: 1.36‐2.05), with an attenuated association observed between diabetes and colorectal cancer without BRAF‐mutations, n = 7959 (ORfully adj: 1.33, 95% CI: 1.19‐1.48). In the case only analysis, BRAF‐mutation was differentially associated with diabetes (Pdifference = .03). For the other markers, associations with diabetes were similar across tumor subtypes. In conclusion, our study confirms the established association between diabetes and colorectal cancer risk, and suggests that it particularly increases the risk of BRAF‐mutated tumors.