Potter, John D.

Abstract Protective associations of fruits, vegetables, and fiber intake with colorectal cancer risk have been shown in many, but not all epidemiologic studies. One possible reason for study heterogeneity is that dietary factors may have distinct effects by colorectal cancer molecular subtypes. Here, we investigate the association of fruit, vegetables, and fiber intake with four well-established colorectal cancer molecular subtypes separately and in combination. Nine observational studies including 9,592 cases with molecular subtypes for microsatellite instability (MSI), CpG island methylator phenotype (CIMP), and somatic mutations in BRAF and KRAS genes, and 7,869 controls were analyzed. Both case-only logistic regression analyses and polytomous logistic regression analyses (with one control set and multiple case groups) were used. Higher fruit intake was associated with a trend toward decreased risk of BRAF-mutated tumors [OR 4th vs. 1st quartile = 0.82 (95% confidence interval, 0.65–1.04)] but not BRAF-wildtype tumors [1.09 (0.97–1.22); P difference as shown in case-only analysis = 0.02]. This difference was observed in case–control studies and not in cohort studies. Compared with controls, higher fiber intake showed negative association with colorectal cancer risk for cases with microsatellite stable/MSI-low, CIMP-negative, BRAF-wildtype, and KRAS-wildtype tumors (Ptrend range from 0.03 to 3.4e-03), which is consistent with the traditional adenoma-colorectal cancer pathway. These negative associations were stronger compared with MSI-high, CIMP-positive, BRAF-mutated, or KRAS-mutated tumors, but the differences were not statistically significant. These inverse associations for fruit and fiber intake may explain, in part, inconsistent findings between fruit or fiber intake and colorectal cancer risk that have previously been reported. Significance: These analyses by colorectal cancer molecular subtypes potentially explain the inconsistent findings between dietary fruit or fiber intake and overall colorectal cancer risk that have previously been reported.

Abstract Background: Regular use of nonsteroidal anti-inflammatory drugs (NSAID) is associated with lower risk of colorectal cancer. Genome-wide interaction analysis on single variants (G × E) has identified several SNPs that may interact with NSAIDs to confer colorectal cancer risk, but variations in gene expression levels may also modify the effect of NSAID use. Therefore, we tested interactions between NSAID use and predicted gene expression levels in relation to colorectal cancer risk. Methods: Genetically predicted gene expressions were tested for interaction with NSAID use on colorectal cancer risk among 19,258 colorectal cancer cases and 18,597 controls from 21 observational studies. A Mixed Score Test for Interactions (MiSTi) approach was used to jointly assess G × E effects which are modeled via fixed interaction effects of the weighted burden within each gene set (burden) and residual G × E effects (variance). A false discovery rate (FDR) at 0.2 was applied to correct for multiple testing. Results: Among the 4,840 genes tested, genetically predicted expression levels of four genes modified the effect of any NSAID use on colorectal cancer risk, including DPP10 (PG×E = 1.96 × 10−4), KRT16 (PG×E = 2.3 × 10−4), CD14 (PG×E = 9.38 × 10−4), and CYP27A1 (PG×E = 1.44 × 10−3). There was a significant interaction between expression level of RP11-89N17 and regular use of aspirin only on colorectal cancer risk (PG×E = 3.23 × 10−5). No interactions were observed between predicted gene expression and nonaspirin NSAID use at FDR < 0.2. Conclusions: By incorporating functional information, we discovered several novel genes that interacted with NSAID use. Impact: These findings provide preliminary support that could help understand the chemopreventive mechanisms of NSAIDs on colorectal cancer.

Interindividual differences in DNA repair systems may play a role in modulating the individual risk of developing colorectal cancer. To better ascertain the role of DNA repair gene polymorphisms on colon and rectal cancer risk individually, we evaluated 15,419 single nucleotide polymorphisms (SNPs) within 185 DNA repair genes using GWAS data from the Colon Cancer Family Registry (CCFR) and the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO), which included 8,178 colon cancer, 2,936 rectum cancer cases and 14,659 controls. Rs1800734 (in MLH1 gene) was associated with colon cancer risk (p‐value = 3.5 × 10−6) and rs2189517 (in RAD51B) with rectal cancer risk (p‐value = 5.7 × 10−6). The results had statistical significance close to the Bonferroni corrected p‐value of 5.8 × 10−6. Ninety‐four SNPs were significantly associated with colorectal cancer risk after Binomial Sequential Goodness of Fit (BSGoF) procedure and confirmed the relevance of DNA mismatch repair (MMR) and homologous recombination pathways for colon and rectum cancer, respectively. Defects in MMR genes are known to be crucial for familial form of colorectal cancer but our findings suggest that specific genetic variations in MLH1 are important also in the individual predisposition to sporadic colon cancer. Other SNPs associated with the risk of colon cancer (e.g., rs16906252 in MGMT) were found to affect mRNA expression levels in colon transverse and therefore working as possible cis‐eQTL suggesting possible mechanisms of carcinogenesis.

Abstract Background: Different molecules within coffee may be pro-carcinogenic, while others may have anti-oxidant and anti-neoplastic effects. Hence, it can be hypothesized that coffee may have differing risk associations with specific colorectal cancer (CRC) subtypes. The aim of this study was to evaluate associations between coffee consumption and incident CRC overall, by anatomic subsite, and by molecularly defined tumor subtypes, in the population-based Iowa Women's Health Study (IWHS). Methods: IWHS participants, ages 55-69 years, provided dietary intake estimates using self-reported food frequency questionnaires at baseline (1986). From 35,216 women included the final analysis, 1298 developed CRC (673 with proximal CRC, 597 with distal CRC) through December 31, 2002 (end of the IWHS molecular epidemiology project study period). Archived, paraffin-embedded tissue specimens for 732 representative CRC cases were collected and analyzed to determine if coffee consumption (none or <1 cup/month, <1 cup/week, 1 cup/day, 2-3 cups/day and ≥4 cups/day) was associated with microsatellite instability (MSI-L/MSS or MSI-H), CpG island methylator phenotype (CIMP-negative or CIMP-positive), BRAF mutation status (BRAF-wildtype or BRAF-mutated), KRAS mutation status (KRAS-wildtype and KRAS-mutated), p53 protein expression level (p53 absent, low or high) and estrogen receptor β protein expression level (ERβ absent, low or high). Multivariable Cox regression models were fit to estimate relative risks (RR) and 95% confidence intervals (CI) with adjustment for age, smoking, alcohol use, diabetes, body mass index and waist-hip ratio, physical activity, estrogen use, dietary intake of fat, calcium, vitamin E, sucrose, red meat and methionine. P-values were not adjusted for multiple testing. Results: There was no significant difference in the risk of overall (RR, 1.05; 95% CI, 0.85-1.29; Ptrend=0.50), proximal (RR, 0.91; 95% CI, 0.68-1.21; Ptrend=0.57) or distal CRC (RR, 1.18; 95% CI, 0.86-1.64; Ptrend=0.19) in women who consumed ≥4 cups of coffee/day (caffeinated or decaffeinated), as compared to non-consumers. By CRC molecular subtypes for these same two exposure groups, women with ≥4 cups of coffee/day were more likely to develop ERβ-low CRC (RR, 2.98; 95% CI, 1.04-8.51; Ptrend=0.06), though no significant associations were observed for caffeinated or decaffeinated coffee specifically. In contrast, women who consumed ≥4 cups of decaffeinated coffee/day were less likely to develop MSI-L/MSS (RR, 0.57; 95% CI, 0.36-0.92; Ptrend=0.04) and BRAF-WT (RR, 0.69; 95% CI, 0.44-1.06; Ptrend=0.06) CRCs than non-consumers. Women who consumed ≥4 cups of caffeinated coffee/day were also less likely to develop p53-low CRCs (RR, 0.53; 95% CI, 0.24-1.15; Ptrend=0.02). Null associations were observed between decaffeinated coffee consumption and risk of CIMP-, KRAS- and p53- defined CRC subtypes, as well as between caffeinated coffee consumption and risk of MSI-, BRAF-, CIMP- or KRAS-defined CRC subtypes. Conclusion: In this prospective cohort study of older women, though coffee consumption was not associated with risk of overall or site-specific CRC, unique associations were observed with molecularly defined CRC subtypes, supporting a need for further investigation of the potential pro- and anti-carcinogenic effects of coffee constituents on biologic targets involved in colorectal carcinogenesis. Citation Format: Siddharth Singh, Robert Vierkant, Lori Tillmans, Alice Wang, Daniel Weisenberger, Peter Laird, Maki Inoue-Choi, Andrea Mariani, Charles Lynch, Kristin Anderson, Amy French, Robert Haile, Lisa Harnack, John Potter, Susan Slager, Thomas Smyrk, Stephen Thibodeau, James Cerhan, Paul Limburg. Coffee consumption and colorectal cancer risk by molecularly defined subtypes in a prospective cohort study of older women. [abstract]. In: Proceedings of the Twelfth Annual AACR International Conference on Frontiers in Cancer Prevention Research; 2013 Oct 27-30; National Harbor, MD. Philadelphia (PA): AACR; Can Prev Res 2013;6(11 Suppl): Abstract nr A44.

Abstract Genome-wide association studies of colorectal cancer (CRC) have identified over 50 susceptibility loci. These variants represent only a small fraction of total heritability for CRC. Gene-environment (GxE) interaction studies may help identify novel loci and biological interactions that give insight to the pathogenesis of CRC. Previous genome-wide GxE studies with dietary factors have identified interactions between loci and processed meat consumption and alcohol; however, limited statistical power remains a primary concern. Set-based SNP testing has the potential to increase statistical power to detect GxE interactions by aggregating functionally relevant SNPs. In this large pooled analysis using 14 case-control and nested case-control studies, we incorporated functional information from the transcriptome prediction tool, PrediXcan, into a novel set-based approach for testing GxE interactions. We used variant weights from the PrediXcan models of tissue-specific gene expression in the transverse colon as a priori variant information for a set-based GxE approach. We restricted our analysis to variants in the PrediXcan transverse colon gene models (n = 4,842). This discovery phase included 10,360 CRC and advanced adenoma cases, and 11,183 controls of European ancestry from the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO) and the Colon Cancer Family Registry. All 14 studies were analyzed together in a pooled data set using the Mixed Effects Score Tests for interactions. We tested for gene interactions with sex- and study- specific quartiles of dietary intake of red meat (servings/day), processed meat (servings/day), vegetables (servings/day), fruits (servings/day), and fiber (g/day). We detected two genes with suggestive interactions (false discovery rate (FDR) < 0.2) with intake of red meat and risk of CRC: Superoxide Dismutase 2 (SOD2) and Ubiquitin Conjugating Enzyme E2 H (UBE2H). No interactions at FDR < 0.2 were observed for processed meat, vegetables, fruits, or fiber. The SOD2 gene, which encodes an enzyme important in apoptotic signaling and clearing of reactive oxygen species, may regulate response to colonic exposure to heme iron or increased bile acid from high-fat content in red meat. UBE2H is part of the ubiquitin-proteasome system that has a role in Wnt signaling, which can be mediated by heme iron in red meat and is commonly found dysregulated in cancer. These findings highlight the efficacy of integrating functional information and set-based testing for novel discovery of genes interacting with known dietary risk factors of CRC. We plan to replicate these findings in additional studies. Citation Format: Paneen S. Petersen, Yu-Ru Su, Sonja I. Berndt, Stephanie A. Bien, Hermann Brenner, Graham Casey, Andrew T. Chan, Jenny Chang-Claude, Jane C. Figueiredo, Steven J. Gallinger, Robert W. Haile, Tabitha A. Harrison, Michael Hoffmeister, Mark A. Jenkins, Amit D. Joshi, Sébastien Küry, Loic Le Marchand, Yi Lin, Noralane M. Lindor, Polly A. Newcomb, John D. Potter, Robert Schoen, Martha L. Slattery, Stephen N. Thibodeau, Emily White, Li Hsu, Ulrike Peters, CCFR, GECCO. Interactions between genetic predictors of gene expression and dietary factors associated with risk of colorectal cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 5268.

Abstract Background Regular use of aspirin and other non-steroidal anti-inflammatory drugs (NSAIDs) is associated with lower risk of colorectal cancer (CRC). Genome-wide interaction analysis (GxE) has identified a few variants that may modify the effects of NSAIDs on CRC risk. However, limited statistical power remains a concern. Restricting analyses by using functional genomic information to aggregate variants into biologically relevant sets can reduce the number of tests, thereby increasing statistical power. We tested the interactions between variant models of gene expression and NSAIDs use on CRC risk. Methods Functional weights of each variant were estimated using PrediXcan based on jointly measured transcriptomes and genomes data of transverse colon tissues from the Genotype-Tissue Expression (GTEx) Project for all genes with sufficient heritability (≥1%). A mixed-effects model was used to assess the GxE effects in a gene among 9,917 incident CRC cases and 10,533 controls from 17 (nested) case-control studies from the Colon Cancer Family Registry (CCFR) and Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). GxE analysis was done by modeling the interaction between predicted gene expression and NSAID use (fixed effects), and residual variant-specific GxE effects that are not accounted for (random effects). Regular user of aspirin and/or non-aspirin NSAIDs was compared to non-regular users. False discovery rate (FDR) was used to account for multiple comparison and FDR≤0.2 was considered genome-wide significant. Secondary analysis was performed on significant genes and duration of use among users. Results Regular use of aspirin was higher among controls (30.3%) compared to cases (24.3%). Among the 4,842 genes tested, SORD significantly modified the effect of regular use of aspirin on CRC risk (p-interaction = 1.45×10-5; FDR=0.07). The SORD gene encodes sorbitol dehydrogenase (SORD), which oxidizes sorbitol to fructose in the polyol pathway. Decreased SORD concentrations were previously observed to dramatically increase in colorectal adenoma cells when compared to normal mucosa cells, suggesting the involvement of dysregulated polyol metabolisms in colorectal tumorigenesis. However, the duration of aspirin use was not statistically significantly associated with SORD gene expression on CRC risk (p=0.149) among 5,560 aspirin users. No significant interactions were observed between genetically determined colon gene expression levels and any NSAID use or non-aspirin NSAID use at FDR<0.2. Conclusions Incorporating functional information, we discovered a novel gene that may interact with aspirin use to confer CRC risk. These findings provide preliminary support for new biological insights that could help understand the chemopreventive mechanisms of aspirin on CRC. We aim to replicate these findings in additional studies. Citation Format: Xiaoliang Wang, Yu-Ru Su, Andrew T. Chan, Stephanie Bien, Sonja I. Bernt, Hermann Brenner, Graham Casey, Jenny Chang-Claude, Steven J. Gallinger, Robert W. Haile, Tabitha A. Harrison, Michael Hoffmeister, Mark A. Jenkins, Amit Joshi, Yi Lin, Noralane M. Lindor, Loic Le Marchand, Hongmei Nan, Polly A. Newcomb, John D. Potter, Martha L. Slattery, Steve N. Thibodeau, Emily White, Li Hsu, Ulrike Peters. Functionally informed genome-wide interaction analysis of nonsteroidal anti-inflammatory drugs on colorectal cancer risk [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 2965.

Abstract Over 50 genetic variants have been associated with colorectal cancer (CRC) risk through genome-wide association studies (GWAS), yet these variants represent only a fraction of the total estimated heritability. CRC is a heterogenous disease with diverse tumor etiology. Assessing genetic risk in molecular subtypes may help to identify novel loci and characterize genetic risk among tumor subtypes. We used microsatellite instability (MSI), an established CRC classifier with etiological and therapeutic relevance, to define CRC subtypes for GWAS analyses. We conducted a case-case analysis to estimate odds ratios (OR) and 95% confidence intervals (CI) for association of genome-wide variants with microsatellite stable (MSS) versus unstable (MSI) carcinomas. We ran an inverse-variance weighted fixed-effects meta-analysis across GWAS in a discovery set of 4,163 population-based CRC cases with harmonized microsatellite instability (MSI) marker and imputed genotype data. For each analysis, we used log-additive logistic regression, adjusting for age, sex, and principal components to account for population substructure. We then followed up with replication of 102 SNPs that reached p-values less than 5x10-6 in 1,698 cases. A total of 845 (20.3%) cancer cases were microsatellite unstable in the discovery population and 174 (10.2%) were unstable in the replication population. No variants reached the genome-wide significance level of 5x10-8 in the discovery set. However, we identified two variants that reached a Bonferroni corrected p-value of 4.0x10-4 in the replication set. This included one variant in MLH1 (Replication: OR=1.74, 95% CI=1.53-1.98, p=1.63x10-5; Discovery+Replication: OR=1.45, 95% CI=1.37-1.54, p=9.76x10-11) and one variant in LOC105377645 (Replication: OR=1.70, 95% CI=1.49-1.94, p=5.13x10-5; Discovery+Replication: OR=1.45, 95% CI=1.37-1.54, p=9.76 x 10-11). The MLH1 gene is a DNA mismatch repair gene implicated in Lynch Syndrome, the hallmark of which is microsatellite instability. This is the first genome-wide scan to identify a common variant in MLH1 that is associated with CRC. This variant (minor allele frequency, MAF = 23% in this all European ancestry population) is located in the 5'-untranslated region of MLH1 and is thought to act as a long-range regulator of DCLK3, a potential tumor driver gene. The second variant, located in LOC105377645 with an MAF of 22%, is in an uncharacterized region of the genome and has not previously been implicated in cancer development. These findings suggest that accounting for molecular heterogeneity is important for discovery and characterization of genetic variants associated with CRC risk. We plan to run polytomous regression analyses, increase our sample size, and further investigate CRC subtypes by CIMP, BRAF mutation, KRAS mutation status. Citation Format: Tabitha A. Harrison, Yiwen Lu, Chenjie Zeng, Flora Qu, Kristin Anderson, Hermann Brenner, Daniel D. Buchanan, Peter T. Campbell, Andrew T. Chan, Jenny Chang-Claude, Graham G. Giles, Bethany Van Guelpen, Michael Hoffmeister, Mark A. Jenkins, Noralane M. Lindor, Roger L. Milne, Polly A. Newcomb, Reiko Nishihara, Michael O. Woods, Shuji Ogino, John D. Potter, Martha L. Slattery, Wei Sun, Stephen N. Thibodeau, Li Hsu, Ulrike Peters. Genome-wide association study by colorectal carcinoma subtype [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 229.

Greater physical activity is associated with a decrease in risk of colorectal cancer for the general population; however, little is known about its relationship with colorectal cancer risk in people with Lynch syndrome, carriers of inherited pathogenic mutations in genes affecting DNA mismatch repair (MMR). We studied a cohort of 2,042 MMR gene mutations carriers (n = 807, diagnosed with colorectal cancer), from the Colon Cancer Family Registry. Self‐reported physical activity in three age‐periods (20–29, 30–49 and ≥50 years) was summarized as average metabolic equivalent of task hours per week (MET‐hr/week) during the age‐period of cancer diagnosis or censoring (near‐term exposure) and across all age‐periods preceding cancer diagnosis or censoring (long‐term exposure). Weighted Cox regression was used to estimate the hazard ratio (HR) and 95% confidence intervals (CI) for the association between physical activity and colorectal cancer risk. Near‐term physical activity was associated with a small reduction in the risk of colorectal cancer (HR ≥35 vs. <3.5 MET‐hr/week, 0.71; 95% CI, 0.53–0.96). The strength and direction of associations were similar for long‐term physical activity, although the associations were not nominally significant. Our results suggest that physical activity is inversely associated with the risk of colorectal cancer for people with Lynch syndrome; however, further confirmation is warranted. The potential modifying effect of physical activity on colorectal cancer risk in people with Lynch syndrome could be useful for risk prediction and support counseling advice for lifestyle modification to reduce cancer risk.

Abstract Background Previous genome-wide association studies (GWAS) have identified 42 loci (P < 5 × 10−8) associated with risk of colorectal cancer (CRC). Expanded consortium efforts facilitating the discovery of additional susceptibility loci may capture unexplained familial risk. Methods We conducted a GWAS in European descent CRC cases and control subjects using a discovery–replication design, followed by examination of novel findings in a multiethnic sample (cumulative n = 163 315). In the discovery stage (36 948 case subjects/30 864 control subjects), we identified genetic variants with a minor allele frequency of 1% or greater associated with risk of CRC using logistic regression followed by a fixed-effects inverse variance weighted meta-analysis. All novel independent variants reaching genome-wide statistical significance (two-sided P < 5 × 10−8) were tested for replication in separate European ancestry samples (12 952 case subjects/48 383 control subjects). Next, we examined the generalizability of discovered variants in East Asians, African Americans, and Hispanics (12 085 case subjects/22 083 control subjects). Finally, we examined the contributions of novel risk variants to familial relative risk and examined the prediction capabilities of a polygenic risk score. All statistical tests were two-sided. Results The discovery GWAS identified 11 variants associated with CRC at P < 5 × 10−8, of which nine (at 4q22.2/5p15.33/5p13.1/6p21.31/6p12.1/10q11.23/12q24.21/16q24.1/20q13.13) independently replicated at a P value of less than .05. Multiethnic follow-up supported the generalizability of discovery findings. These results demonstrated a 14.7% increase in familial relative risk explained by common risk alleles from 10.3% (95% confidence interval [CI] = 7.9% to 13.7%; known variants) to 11.9% (95% CI = 9.2% to 15.5%; known and novel variants). A polygenic risk score identified 4.3% of the population at an odds ratio for developing CRC of at least 2.0. Conclusions This study provides insight into the architecture of common genetic variation contributing to CRC etiology and improves risk prediction for individualized screening.