Campbell, Peter T.

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 <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 <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 <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 <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<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 <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 < 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 (>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 > 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.

AbstractBackgroundThe incidence of colorectal cancer (CRC) among individuals aged younger than 50 years has been increasing. As screening guidelines lower the recommended age of screening initiation, concerns including the burden on screening capacity and costs have been recognized, suggesting that an individualized approach may be warranted. We developed risk prediction models for early-onset CRC that incorporate an environmental risk score (ERS), including 16 lifestyle and environmental factors, and a polygenic risk score (PRS) of 141 variants.MethodsRelying on risk score weights for ERS and PRS derived from studies of CRC at all ages, we evaluated risks for early-onset CRC in 3486 cases and 3890 controls aged younger than 50 years. Relative and absolute risks for early-onset CRC were assessed according to values of the ERS and PRS. The discriminatory performance of these scores was estimated using the covariate-adjusted area under the receiver operating characteristic curve.ResultsIncreasing values of ERS and PRS were associated with increasing relative risks for early-onset CRC (odds ratio per SD of ERS = 1.14, 95% confidence interval [CI] = 1.08 to 1.20; odds ratio per SD of PRS = 1.59, 95% CI = 1.51 to 1.68), both contributing to case-control discrimination (area under the curve = 0.631, 95% CI = 0.615 to 0.647). Based on absolute risks, we can expect 26 excess cases per 10 000 men and 21 per 10 000 women among those scoring at the 90th percentile for both risk scores.ConclusionsPersonal risk scores have the potential to identify individuals at differential relative and absolute risk for early-onset CRC. Improved discrimination may aid in targeted CRC screening of younger, high-risk individuals, potentially improving outcomes.

AbstractElevated blood levels of C‐reactive protein (CRP) have been linked to colorectal cancer (CRC) survival. We evaluated genetic variants associated with CRP levels and their interactions with sex and lifestyle factors in association with CRC‐specific mortality. Our study included 16 142 CRC cases from the International Survival Analysis in Colorectal Cancer Consortium. We identified 618 common single nucleotide polymorphisms (SNPs) associated with CRP levels from the NHGRI‐EBI GWAS Catalog. Cox proportional hazards regression was used to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for associations between SNPs and CRC‐specific mortality adjusting for age, sex, genotyping platform/study and principal components. We investigated their interactions with sex and lifestyle factors using likelihood ratio tests. Of 5472 (33.9%) deaths accrued over up to 10 years of follow‐up, 3547 (64.8%) were due to CRC. No variants were associated with CRC‐specific mortality after multiple comparison correction. We observed strong evidence of interaction between variant rs1933736 at FRK gene and sex in relation to CRC‐specific mortality (corrected Pinteraction = .0004); women had higher CRC‐specific mortality associated with the minor allele (HR = 1.11, 95% CI = 1.04‐1.19) whereas an inverse association was observed for men (HR = 0.88, 95% CI = 0.82‐0.94). There was no evidence of interactions between CRP‐associated SNPs and alcohol, obesity or smoking. Our study observed a significant interaction between sex and a CRP‐associated variant in relation to CRC‐specific mortality. Future replication of this association and functional annotation of the variant are needed.

Abstract Background Observational studies have shown higher folate consumption to be associated with lower risk of colorectal cancer (CRC). Understanding whether and how genetic risk factors interact with folate could further elucidate the underlying mechanism. Aggregating functionally relevant genetic variants in set-based variant testing has higher power to detect gene-environment (GxE) interactions and may provide information on the underlying biological pathway. Objective We investigated interactions between folate consumption and predicted gene expression on colorectal cancer risk across the genome. Methods We used variant weights from the PrediXcan models of colon tissue-specific gene expression as a priori variant information for a set-based GxE approach. We harmonized total folate intake (mcg/day) based on dietary intake and supplemental use across cohort and case-control studies and calculated sex and study specific quantiles. Analyses were performed using a mixed effects score tests for interactions between folate and genetically predicted expression of 4,839 genes with available genetically predicted expression. We meta-analyzed results across 23 studies for a total of 13,498 cases with colorectal tumors and 13,918 controls of European ancestry. Results We found suggestive evidence of interaction with folate intake for genes including glutathione S-Transferase Alpha 1 (GSTA1; p=4.3E-4), Tonsuko Like, DNA Repair Protein (TONSL; p=4.3E-4), and Aspartylglucosaminidase (AGA: p=4.5E-4). Glutathione is an antioxidant, preventing cellular damage and is a downstream metabolite of homocysteine and metabolized by GSTA1. TONSL is part of a complex that functions in the recovery of double strand breaks and AGA plays a role in lysosomal breakdown of glycoprotein. Conclusion We identified three genes involved in preventing or repairing DNA damage that may interact with folate consumption to alter CRC risk.

Salicylic acid (SA) has observationally been shown to decrease colorectal cancer (CRC) risk. Aspirin (acetylsalicylic acid, that rapidly deacetylates to SA) is an effective primary and secondary chemopreventive agent. Through a Mendelian randomization (MR) approach, we aimed to address whether levels of SA affected CRC risk, stratifying by aspirin use. A two-sample MR analysis was performed using GWAS summary statistics of SA (INTERVAL and EPIC-Norfolk, N = 14,149) and CRC (CCFR, CORECT, GECCO and UK Biobank, 55,168 cases and 65,160 controls). The DACHS study (4410 cases and 3441 controls) was used for replication and stratification of aspirin-use. SNPs proxying SA were selected via three methods: (1) functional SNPs that influence the activity of aspirin-metabolising enzymes; (2) pathway SNPs present in enzymes’ coding regions; and (3) genome-wide significant SNPs. We found no association between functional SNPs and SA levels. The pathway and genome-wide SNPs showed no association between SA and CRC risk (OR: 1.03, 95% CI: 0.84–1.27 and OR: 1.08, 95% CI: 0.86–1.34, respectively). Results remained unchanged upon aspirin use stratification. We found little evidence to suggest that an SD increase in genetically predicted SA protects against CRC risk in the general population and upon stratification by aspirin use.