Brenner, Hermann

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.

AbstractMotivationThe use of single nucleotide polymorphism (SNP) interactions to predict complex diseases is getting more attention during the past decade, but related statistical methods are still immature. We previously proposed the SNP Interaction Pattern Identifier (SIPI) approach to evaluate 45 SNP interaction patterns/patterns. SIPI is statistically powerful but suffers from a large computation burden. For large-scale studies, it is necessary to use a powerful and computation-efficient method. The objective of this study is to develop an evidence-based mini-version of SIPI as the screening tool or solitary use and to evaluate the impact of inheritance mode and model structure on detecting SNP–SNP interactions.ResultsWe tested two candidate approaches: the ‘Five-Full’ and ‘AA9int’ method. The Five-Full approach is composed of the five full interaction models considering three inheritance modes (additive, dominant and recessive). The AA9int approach is composed of nine interaction models by considering non-hierarchical model structure and the additive mode. Our simulation results show that AA9int has similar statistical power compared to SIPI and is superior to the Five-Full approach, and the impact of the non-hierarchical model structure is greater than that of the inheritance mode in detecting SNP–SNP interactions. In summary, it is recommended that AA9int is a powerful tool to be used either alone or as the screening stage of a two-stage approach (AA9int+SIPI) for detecting SNP–SNP interactions in large-scale studies.Availability and implementationThe ‘AA9int’ and ‘parAA9int’ functions (standard and parallel computing version) are added in the SIPI R package, which is freely available at https://linhuiyi.github.io/LinHY_Software/.Supplementary informationSupplementary data are available at Bioinformatics online.

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.

AbstractEpidemiological, biological, and molecular data suggest links between endometriosis and endometrial cancer, with recent epidemiological studies providing evidence for an association between a previous diagnosis of endometriosis and risk of endometrial cancer. We used genetic data as an alternative approach to investigate shared biological etiology of these two diseases. Genetic correlation analysis of summary level statistics from genomewide association studies (GWAS) using LD Score regression revealed moderate but significant genetic correlation (rg = 0.23, P = 9.3 × 10−3), and SNP effect concordance analysis provided evidence for significant SNP pleiotropy (P = 6.0 × 10−3) and concordance in effect direction (P = 2.0 × 10−3) between the two diseases. Cross‐disease GWAS meta‐analysis highlighted 13 distinct loci associated at P ≤ 10−5 with both endometriosis and endometrial cancer, with one locus (SNP rs2475335) located within PTPRD associated at a genomewide significant level (P = 4.9 × 10−8, OR = 1.11, 95% CI = 1.07–1.15). PTPRD acts in the STAT3 pathway, which has been implicated in both endometriosis and endometrial cancer. This study demonstrates the value of cross‐disease genetic analysis to support epidemiological observations and to identify biological pathways of relevance to multiple diseases.

Limited epidemiological evidence suggests that the etiology of hormone receptor positive (HR+) breast cancer may differ by levels of histologic grade and proliferation. We pooled risk factor and pathology data on 5,905 HR+ breast cancer cases and 26,281 controls from 11 epidemiological studies. Proliferation was determined by centralized automated measures of KI67 in tissue microarrays. Odds ratios (OR), 95% confidence intervals (CI) and p‐values for case–case and case–control comparisons for risk factors in relation to levels of grade and quartiles (Q1–Q4) of KI67 were estimated using polytomous logistic regression models. Case–case comparisons showed associations between nulliparity and high KI67 [OR (95% CI) for Q4 vs. Q1 = 1.54 (1.22, 1.95)]; obesity and high grade [grade 3 vs. 1 = 1.68 (1.31, 2.16)] and current use of combined hormone therapy (HT) and low grade [grade 3 vs. 1 = 0.27 (0.16, 0.44)] tumors. In case–control comparisons, nulliparity was associated with elevated risk of tumors with high but not low levels of proliferation [1.43 (1.14, 1.81) for KI67 Q4 vs. 0.83 (0.60, 1.14) for KI67 Q1]; obesity among women ≥50 years with high but not low grade tumors [1.55 (1.17, 2.06) for grade 3 vs. 0.88 (0.66, 1.16) for grade 1] and HT with low but not high grade tumors [3.07 (2.22, 4.23) for grade 1 vs. 0.85 (0.55, 1.30) for grade 3]. Menarcheal age and family history were similarly associated with HR+ tumors of different grade or KI67 levels. These findings provide insights into the etiologic heterogeneity of HR+ tumors.