Bien, Stephanie A.

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.

AbstractBackgroundBody mass index (BMI) and diabetes are established risk factors for colorectal cancer (CRC), likely through perturbations in metabolic traits (e.g. insulin resistance and glucose homeostasis). Identification of interactions between variation in genes and these metabolic risk factors may identify novel biologic insights into CRC etiology.MethodsTo improve statistical power and interpretation for gene‐environment interaction (G × E) testing, we tested genetic variants that regulate expression of a gene together for interaction with BMI (kg/m2) and diabetes on CRC risk among 26 017 cases and 20 692 controls. Each variant was weighted based on PrediXcan analysis of gene expression data from colon tissue generated in the Genotype‐Tissue Expression Project for all genes with heritability ≥1%. We used a mixed‐effects model to jointly measure the G × E interaction in a gene by partitioning the interactions into the predicted gene expression levels (fixed effects), and residual G × E effects (random effects). G × BMI analyses were stratified by sex as BMI‐CRC associations differ by sex. We used false discovery rates to account for multiple comparisons and reported all results with FDR <0.2.ResultsAmong 4839 genes tested, genetically predicted expressions of FOXA1 (P = 3.15 × 10−5), PSMC5 (P = 4.51 × 10−4) and CD33 (P = 2.71 × 10−4) modified the association of BMI on CRC risk for men; KIAA0753 (P = 2.29 × 10−5) and SCN1B (P = 2.76 × 10−4) modified the association of BMI on CRC risk for women; and PTPN2 modified the association between diabetes and CRC risk in both sexes (P = 2.31 × 10−5).ConclusionsAggregating G × E interactions and incorporating functional information, we discovered novel genes that may interact with BMI and diabetes on CRC risk.

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 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) &lt; 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 &lt; 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&lt;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 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.