French, Amy J.

Abstract Motivation Interpreting genetic variation in noncoding regions of the genome is an important challenge for personal genome analysis. One mechanism by which noncoding single nucleotide variants (SNVs) influence downstream phenotypes is through the regulation of gene expression. Methods to predict whether or not individual SNVs are likely to regulate gene expression would aid interpretation of variants of unknown significance identified in whole-genome sequencing studies. Results We developed FIRE (Functional Inference of Regulators of Expression), a tool to score both noncoding and coding SNVs based on their potential to regulate the expression levels of nearby genes. FIRE consists of 23 random forests trained to recognize SNVs in cis-expression quantitative trait loci (cis-eQTLs) using a set of 92 genomic annotations as predictive features. FIRE scores discriminate cis-eQTL SNVs from non-eQTL SNVs in the training set with a cross-validated area under the receiver operating characteristic curve (AUC) of 0.807, and discriminate cis-eQTL SNVs shared across six populations of different ancestry from non-eQTL SNVs with an AUC of 0.939. FIRE scores are also predictive of cis-eQTL SNVs across a variety of tissue types. Availability and implementation FIRE scores for genome-wide SNVs in hg19/GRCh37 are available for download at https://sites.google.com/site/fireregulatoryvariation/. Supplementary information Supplementary data are available at Bioinformatics online.

Abstract The Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO) in collaboration with the Colorectal Cancer Family Registry (CCFR) aims to identify genetic variants and environmental risk factors that impact colorectal cancer (CRC). Over 30 studies from North America, Europe, and Australia participate in GECCO. These studies have collected clinical, epidemiological, and survival data, as well as blood and tumor biospecimens, for over 80,000 CRC cases and controls. The current study aims to conduct targeted deep sequencing of tumors and matching normal DNA to identify recurrent and novel somatic and germline variants in 4,200 CRC cases. To achieve this goal, an AmpliSeq targeted sequencing panel of 1.12 Mbp was constructed to sequence the coding regions of 190 significantly mutated genes identified from whole exome sequencing datasets generated by the Nurses’ Health Study and Health Professional’s Follow-up Study, and The Cancer Genome Atlas. The panel also covers coding regions of 15 genes with germline high penetrance mutations in CRC, 54 regions to detect CRC-related copy number alterations (CNAs), and microsatellite and homopolymer repeat regions to identify defective DNA mismatch repair. Primers were also included to detect Fusobacterium nucleatum DNA in tumor biopsies, as F. nucleatum is thought to promote CRC carcinogenesis. Sequencing of the DNA libraries on Illumina HiSeq 2500 produced a mean coverage of greater than 500X for tumor DNA and 100X for normal DNA, with >85% of the bases covered at the target at 50x. So far, targeted sequencing of >1,500 DNA samples from CRC tumors and normal tissues has identified recurrent and novel somatic mutations, germline genetic variants, and hypermutation status of the tumors due to defective DNA mismatch repair or pathogenic mutations in the POLE gene. Targeted sequencing has also allowed quantification of the F. nucleatum DNA in tumor biopsies; the results were validated by a multiplex QPCR assay. At the AACR annual meeting, we will present targeted sequencing results generated from the first two GECCO-participating studies (n=1,300 cases). These data will be valuable for future association testing of somatic mutations, CNAs, hypermutation status, and F. nucleatum with germline genetic variants, lifestyle, and environmental risk factors and survival. This large study will allow development of better strategies for diagnosis, treatment, and prevention of CRC. Citation Format: Syed H. Zaidi, Wei Sun, Jeroen Huyghe, Catherine S. Grasso, Quang Trinh, Charles Connolly, Amy French, Jasmine Mu, Marios Giannakis, Eve Shinbrot, Ivan Borozan, Michael J. Quist, Hermann Brenner, Daniel Buchanan, Peter Campbell, Andrew Chan, Jenny Chang-Claude, Vincent Ferretti, Charles Fuchs, Andrea Gsur, Marc Gunter, Tabitha Harrison, Michael Hoffmeister, Wen-Yi Huang, Paul Krzyzanowski, Stephen Lee, Mathieu Lemire, Jessica Miller, Danielle Pasternack, Cherie Teney, Elaine Mardis, Polly Newcomb, Lincoln Stein, Lee Timms, David Wheeler, Christina Yung, Niha Zubair, Levi Garraway, Shuji Ogino, Li Hsu, Steven Gallinger, Stephen Thibodeau, Thomas Hudson, Ulrike Peters. Targeted deep sequencing of colorectal tumor tissues to study associations of tumor subtypes with germline genetic, lifestyle, and environmental risk factors [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1286. doi:10.1158/1538-7445.AM2017-1286

Abstract The completion of The Cancer Genome Atlas (TCGA) project for colorectal cancer (CRC) has enabled a new, focused phase of sequencing tumor samples for which genome-wide genetic, epidemiological, clinical and lifestyle data have been collected. By combining somatic mutational profiles with these aforementioned data, we may identify and report effective prevention and treatment approaches for a broader population of individuals. The advent of targeted deep sequencing approaches using DNA obtained from formalin fixed paraffin embedded tissues has made possible the genetic characterization of the large numbers of patients needed to make such an effort relevant. As a first step, we describe a custom gene panel designed from large-scale studies for targeted deep sequencing, and its application to over 4,200 CR tumors, collected by the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). This study is designed to identify recurrent somatic mutations and copy number alterations (CNAs) for association testing with germline genetic and lifestyle and environmental risk factors for CRC, and thereby identify relevant approaches to impact cancer prevention. The targeted CRC panel includes 205 genes. Genes were primarily selected as significantly mutated genes identified from the Nurses’ Health Study and Health Professional's Follow-up Study (n = 700), and TCGA (n = 525). We also included 15 genes associated with high penetrance germline mutations and augmented the list to include genes in pathways of somatically altered genes, identified by literature review, from public databases and known to be associated with loss of heterozygosity. For these 205 genes, amplification primers were designed to include all coding regions of transcripts that are listed in the UCSC Genome Browser. For regions with CNAs, the TCGA dataset was analyzed to include regions with greater than or equal to 2 copy focal gains or losses that were found in more than 4 or 3 tumors, respectively. Candidate genes in regions with significant CNAs from the TCGA CRC publication (Nature 2012) were also included. For CNAs, 6 to 12 amplicons were designed for each of the 56 selected regions (32 gains and 24 losses). Additional target regions include: 1) 25 microsatellite loci previously used to identify defective DNA mismatch repair and 212 homoploymer repeats. These include microsatellite loci recommended by the NCI Consensus Panel for identifying MSI; 2) amelogenin (for gender); and 3) Fusobacterium to detect a putative CRC-associated pathogen in tumor biopsies. At the AACR annual meeting, we expect to present results for deep sequencing (∼1,000X) of the first 1,000 CR tumors, including any preliminarily identified pathways and subtypes that may provide the basis for association testing with germline genetic and lifestyle and environmental risk factors needed for inferring better approaches to prevention and treatment of CRC. Citation Format: Syed H. Zaidi, Catie Grasso, Jasmine Mu, Eve Shinbrot, Marios Giannakis, Charles Connolly, Ivan Borozan, Hermann Brenner, Peter Campbell, Andrew Chan, Jenny Chang-Claude, Mengmeng Du, Vincent Ferretti, Amy French, Charles Fuchs, Steven Gallinger, Levi Garraway, Andrea Gsur, Marc Gunter, Tabitha Harrison, Michael Hoffmeister, Li Hsu, Wen-Yi Huang, Jeroen Huyghe, Mathieu Lemire, Elaine Mardis, John McPherson, Polly Newcomb, Lincoln Stein, Wei Sun, Lee Timms, Quang Trinh, David Wheeler, Christina Yung, Niha Zubair, Shuji Ogino, Stephen Thibodeau, Ulrike Peters, Thomas Hudson. Linking the molecular profile of colorectal tumors to germline genetic and environmental risk factors. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 5221.

Abstract A genome-wide association study (GWAS) of prostate cancer in Kaiser Permanente health plan members (7,783 cases, 38,595 controls; 80.3% non-Hispanic white, 4.9% African-American, 7.0% East Asian, and 7.8% Latino) revealed a new independent risk indel rs4646284 at the previously identified locus 6q25.3 that replicated in PEGASUS (N = 7,539) and the Multiethnic Cohort (N = 4,679) with an overall P = 1.0 × 10−19 (OR, 1.18). Across the 6q25.3 locus, rs4646284 exhibited the strongest association with expression of SLC22A1 (P = 1.3 × 10−23) and SLC22A3 (P = 3.2 × 10−52). At the known 19q13.33 locus, rs2659124 (P = 1.3 × 10−13; OR, 1.18) nominally replicated in PEGASUS. A risk score of 105 known risk SNPs was strongly associated with prostate cancer (P < 1.0 × 10−8). Comparing the highest to lowest risk score deciles, the OR was 6.22 for non-Hispanic whites, 5.82 for Latinos, 3.77 for African-Americans, and 3.38 for East Asians. In non-Hispanic whites, the 105 risk SNPs explained approximately 7.6% of disease heritability. The entire GWAS array explained approximately 33.4% of heritability, with a 4.3-fold enrichment within DNaseI hypersensitivity sites (P = 0.004). Significance: Taken together, our findings of independent risk variants, ethnic variation in existing SNP replication, and remaining unexplained heritability have important implications for further clarifying the genetic risk of prostate cancer. Our findings also suggest that there may be much promise in evaluating understudied variation, such as indels and ethnically diverse populations. Cancer Discov; 5(8); 878–91. ©2015 AACR. This article is highlighted in the In This Issue feature, p. 783

Abstract Background: Cigarette smoking (smoking), hormone therapy (MHT), and folate intake (folate) are each thought to influence colorectal cancer risk, but the underlying molecular mechanisms remain incompletely defined. Expression of estrogen receptor β (ESR2) has been associated with colorectal cancer stage and survival. Methods: In this prospective cohort study, we examined smoking, MHT, and folate-associated colorectal cancer risks by ESR2 protein expression level among participants in the Iowa Women's Health Study (IWHS). Self-reported exposure variables were assessed at baseline. Archived, paraffin-embedded colorectal cancer tissue specimens were collected and evaluated for ESR2 protein expression by IHC. Multivariate Cox regression models were fit to estimate relative risks (RR) and 95% confidence intervals (CI) for associations between smoking, MHT, or folate and ESR2-defined colorectal cancer subtypes. Results: Informative environmental exposure and protein expression data were available for 491 incident colorectal cancer cases. Positive associations between ESR2-low and -high tumors and several smoking-related variables were noted, most prominently with average number of cigarettes per day (RR, 4.24; 95% CI, 1.81–9.91 for ESR2-low and RR, 2.15; 95% CI, 1.05–4.41 for ESR2-high for ≥40 cigarettes compared with nonsmokers). For MHT, a statistically significant association with ESR2-low tumors was observed with longer duration of exposure (RR, 0.54; 95% CI, 0.26–1.13 for >5 years compared with never use). No associations were found for folate. Conclusions: In this study, smoking and MHT were associated with ESR2 expression patterns. Impact: These data support possible heterogeneous effects from smoking and MHT on ERβ-related pathways of colorectal carcinogenesis in older women. Cancer Epidemiol Biomarkers Prev; 24(4); 713–9. ©2015 AACR.

AbstractGenetic changes in colon cancer are known to parallel the tissue abnormalities associated with the disease, namely adenoma and adenocarcinoma. The role of microRNA dysregulation in dysplastic progression, however, is not well understood. Here, we show that miR‐182 and miR‐503 undergo sequential up‐regulation and drive the progression of colon adenoma to adenocarcinoma by cooperatively down‐regulating the tumour suppressor FBXW7. We identified that increased expression of miR‐182 is a feature of adenomas. A subsequent increase in miR‐503 expression works cooperatively with miR‐182 to induce transformation of an adenoma to adenocarcinoma. We show that introducing miR‐503 into AAC1 cells, which are derived from a benign adenoma, confers tumourigenic potential. We also demonstrated that blocking both miR‐182 and miR‐503 in HCT116 colon cancer cells resulted in increased FBXW7 expression and significantly reduced tumour size in xenograft models. We confirmed relevance of these results in patients by examining the expression levels of miR‐182 and miR‐503 in over 200 colon cancer patients with 12 year survival outcome data. Decreased patient survival was correlated with elevated expression of both miRNAs, suggesting that elevated levels of both miR‐182 and miR‐503 define a novel prognostic biomarker for colon cancer patients. In conclusion, we show that a sequential expression of miR‐182 and miR‐503 in benign adenoma cooperatively regulates the tumour suppressor FBXW7, contributing to the malignant transformation of colon adenoma to adenocarcinoma and miR‐182 and miR‐503 may prove to be novel therapeutic targets. Array data are available at: http://www.oncomir.umn.edu/ Copyright © 2014 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd

Abstract Background: Hormone therapy (HT) may prevent colorectal cancer (CRC), but the molecular mechanisms have not been fully defined. Expression of estrogen receptor beta (ERβ) has been associated with CRC stage and survival and may influence HT-related prevention pathways. In this prospective cohort study, we examined HT-associated CRC risk by ERβ protein expression levels in the population-based Iowa Women's Health Study (IWHS). Methods: The IWHS recruited 41,836 randomly selected Iowa women, ages 55-69 years, with a valid driver's license at study entry (1986). Exposure variables were assessed at baseline, by self-report. Incident CRC cases were ascertained by annual linkage with the Iowa Cancer Registry. Archived, paraffin-embedded tissue specimens were collected from 732 CRC cases, diagnosed through 12/31/2002, for molecular epidemiology studies. HT was categorized as never, former, and current and with durations of <=5 years or >5 years. ERβ protein expression was determined by immunohistochemistry on tissue microarrays. For each tumor core, both staining intensity (0-3) and proportion of cells stained (0-5) were scored. These two scores were added to determine negative (0), low (1-5) or high (6-8) ERβ protein expression levels. Multivariate Cox regression models were fit to estimate relative risks (RRs) and 95% confidence intervals (CIs). Results: Complete HT and ERβ data were available from 484 CRC cases. Associations between HT and CRC subtypes defined by ERβ expression level are shown in the Table. Conclusions: Data from the IWHS cohort suggest that current and longer-duration HT exposure favorably influence the risk for ERβ-neg or ERβ-low CRC cases to a greater degree than ERβ-high CRC cases. Further investigation is required to confirm these observations, and to clarify the underlying mechanisms. ERβ Neg or Low ERβ High Hormone Therapy Person years Events, N RR* (95% CI) Events, N RR* (95% CI) Never 341377 128 1.00 (Ref.) 196 1.00 (Ref.) Ever 212696 62 0.74 (0.54,1.03) 98 0.82 (0.63,1.07) Former 151535 48 0.79 (0.56,1.12) 67 0.73 (0.54,0.99) Current 61161 14 0.60 (0.33,1.09) 31 1.11 (0.74,1.64) P-trend 0.05 0.52 Duration < = 5yrs 148704 47 0.84 (0.59,1.19) 68 0.80 (0.60,1.08) >5yrs 60064 14 0.5 (0.26,0.93) 28 0.89 (0.59,1.34) P-trend 0.03 0.25 Citation Format: Lori S. Tillmans, Robert A. Vierkant, Alice H. Wang, Charles F. Lynch, Kristin E. Anderson, Amy J. French, Robert W. Haile, John D. Potter, Susan L. Slager, Thomas C. Smyrk, Stephen N. Thibodeau, James R. Cerhan, Paul J. Limburg. Hormone therapy and incident colorectal cancer risk by ER beta protein expression level in a population-based cohort of older women. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 3652. doi:10.1158/1538-7445.AM2013-3652

Abstract Cigarette smoking (CS), hormone therapy (HT), and folate intake (FI) are each thought to influence colorectal cancer risk, but the underlying molecular mechanisms remain incompletely defined. The TP53 (p53) protein, encoded by the TP53 tumor-suppressor gene that is commonly mutated in colorectal cancer, can be readily assessed to differentiate biologically distinct colorectal cancer subtypes. In this prospective cohort study, we examined CS-, HT-, and FI-associated colorectal cancer risks by TP53 protein expression level among Iowa Women's Health Study (IWHS) participants. The IWHS recruited 41,836 randomly selected Iowa women, ages 55 to 69 years, with a valid driver's license at study entry in 1986. Self-reported exposure variables were assessed at baseline. Incident colorectal cancer cases were ascertained by annual linkage with the Iowa Cancer Registry. Archived, paraffin-embedded tissue specimens were collected and evaluated for TP53 protein expression by immunohistochemistry. Multivariate Cox regression models were fit to estimate relative risks (RR) and 95% confidence intervals (CI) for associations between CS, HT, or FI and TP53-defined colorectal cancer subtypes. Informative environmental exposure and protein expression data were available for 492 incident colorectal cancer cases: 222 (45.1%) TP53 negative, 72 (14.6%) TP53 low, and 198 (40.2%) TP53 high. Longer duration (>5 years) of HT was inversely associated with TP53 high colorectal cancers (RR, 0.50; 95% CI, 0.27–0.94). No other statistically significant associations were observed. These data support possible heterogeneous effects from HT on TP53-related pathways of colorectal carcinogenesis in older women. Cancer Epidemiol Biomarkers Prev; 23(2); 350–5. ©2013 AACR.

Abstract Background: Cigarette smoking (CS) is an established colorectal cancer (CRC) risk factor. The p53 protein, encoded by the TP53 tumor suppressor gene that is commonly mutated in CRC, is important in cell cycle arrest needed to repair DNA damage or induce apoptosis to prevent tumor propagation. In this prospective cohort study, we examined CS-associated CRC risks by p53 protein expression level in the population-based Iowa Women's Health Study (IWHS). Methods: The IWHS recruited 41,836 randomly selected Iowa women, ages 55-69 years, with a valid driver's license at study entry (1986). CS and other exposure variables were assessed at baseline, by self-report. Incident CRC cases were ascertained by annual linkage with the Iowa Cancer Registry. Archived, paraffin-embedded tissue specimens were collected from 732 CRC cases diagnosed through 12/31/2002. CS was categorized by smoking status (never, current, former), average number of cigarettes per day (1-19, 20-39, > 40) and cumulative pack-years (1-19, 20-39, > 40). P53 protein expression was determined by immunohistochemistry on tissue microarrays. Tumor cores were scored by combining both staining intensity and percent of cells stained to determine negative, low or high p53 protein expression levels. Multivariate Cox regression models were fit to estimate relative risks (RRs) and 95% confidence intervals (CIs). Results: Of the 732 available CRC cases, 486 (66%) had complete CS and p53 data. Associations between CS and CRC subtypes defined by p53 expression level are shown in the Table. Conclusions: These data demonstrate no clear heterogeneity of CS with CRC subtypes defined by p53 expression levels among older women. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 1677. doi:1538-7445.AM2012-1677