Brenner, Hermann

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 Purpose: Expression of the tumor suppressor gene E-cadherin is diminished in lobular breast cancers and has been implicated in epithelial mesenchymal transition. We assessed risk factor associations for breast cancer stratified by low vs. high E-cadherin protein expression in a pooled analysis within the Breast Cancer Association Consortium (BCAC) studies. Methods: E-cadherin tumor tissue staining was performed centrally at the NCI on formalin-fixed paraffin-embedded tissue microarray (TMA) sections representing 6,010 breast cancer patients from 12 US and European BCAC studies. TMAs were digitally scanned and scored using the SlidePath Digital Image Hub (Leica Biosystems, Wetzlar, Germany). For 5,896 cancers with evaluable tumors, E-cadherin was visually scored as estimated percent of positive cells times stain intensity (0, 1+, 2+, 3+) (score range 0-300). E-cadherin low was defined as tumors with a score < 100. Risk factor associations for low vs. high E-cadherin expressing tumors were evaluated by logistic regression, adjusted for age and study site, and stratified by ER status and histologic subtype. To assess the consistency in results by study, meta-analyses were performed using the random-effects model. All statistical tests were two-sided. Results: E-cadherin low cancers comprised 20% of tumors and were associated with lobular histology, well/moderately differentiated cancers, > 2cm in size, and HER2-negative status (χ2, P < 0.003 for all factors). E-cadherin low status was associated with family history of breast cancer (FH) (OR = 1.32, 95% CI = 1.09-1.60, P-het = 0.005) and ever use of menopausal hormones (OR = 1.26, 95% CI = 1.02-1.56, P-het = 0.03). Study specific meta-plots showed consistent effects across studies for menopausal hormone therapy (I2 = 0.0%, p = 0.64); however, E-cadherin loss by FH did show evidence of heterogeneity by study (I2 = 54%, P = 0.03). Differences in E-cadherin expression remained significant for FH, and menopausal hormone use when further adjusted for ER (P < 0.05). We also found relationships with E-cadherin loss to vary by BMI, number of live births, age at first birth, and oral contraceptive use in stratified analysis by ER status and histologic subtype. Conclusion: This large pooled analysis shows that breast cancer risk factor associations may differ by E-cadherin expression independent of ER status, suggesting that it may represent a marker of etiologic heterogeneity. Citation Format: Hisani N. Horne, Hannah Oh, Mark E. Sherman, Maya Palakal, Stephen H. Hewitt, Marjanka Schmidt, Javier Benitez, Roger Milne, Hermann Brenner, Heli Nevanlinna, Arto Mannermaa, Georgia Chenevix-Trench, Fergus Couch, Peter Devilee, Diana Eccles, Maartje Hooning, Anthony J. Swerdlow, Nick Orr, Melissa A. Troester, Renata Cora, Paul D. Pharoah, Montserrat Garcia-Closas, Jonine D. Figueroa. Breast cancer risk factor associations by loss of E-cadherin tumor tissue expression: A pooled analysis of 5,896 cases in 12 studies from the Breast Cancer Association Consortium (BCAC). [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 3451.

Abstract Breast, ovarian, and prostate cancers are hormone-related and may have a shared genetic basis, but this has not been investigated systematically by genome-wide association (GWA) studies. Meta-analyses combining the largest GWA meta-analysis data sets for these cancers totaling 112,349 cases and 116,421 controls of European ancestry, all together and in pairs, identified at P < 10−8 seven new cross-cancer loci: three associated with susceptibility to all three cancers (rs17041869/2q13/BCL2L11; rs7937840/11q12/INCENP; rs1469713/19p13/GATAD2A), two breast and ovarian cancer risk loci (rs200182588/9q31/SMC2; rs8037137/15q26/RCCD1), and two breast and prostate cancer risk loci (rs5013329/1p34/NSUN4; rs9375701/6q23/L3MBTL3). Index variants in five additional regions previously associated with only one cancer also showed clear association with a second cancer type. Cell-type–specific expression quantitative trait locus and enhancer–gene interaction annotations suggested target genes with potential cross-cancer roles at the new loci. Pathway analysis revealed significant enrichment of death receptor signaling genes near loci with P < 10−5 in the three-cancer meta-analysis. Significance: We demonstrate that combining large-scale GWA meta-analysis findings across cancer types can identify completely new risk loci common to breast, ovarian, and prostate cancers. We show that the identification of such cross-cancer risk loci has the potential to shed new light on the shared biology underlying these hormone-related cancers. Cancer Discov; 6(9); 1052–67. ©2016 AACR. This article is highlighted in the In This Issue feature, p. 932

Circulating insulin‐like growth factors (IGFs) and their binding proteins (IGFBPs) are associated with prostate cancer. Using genetic variants as instruments for IGF peptides, we investigated whether these associations are likely to be causal. We identified from the literature 56 single nucleotide polymorphisms (SNPs) in the IGF axis previously associated with biomarker levels (8 from a genome‐wide association study [GWAS] and 48 in reported candidate genes). In ∼700 men without prostate cancer and two replication cohorts (N ∼ 900 and ∼9,000), we examined the properties of these SNPS as instrumental variables (IVs) for IGF‐I, IGF‐II, IGFBP‐2 and IGFBP‐3. Those confirmed as strong IVs were tested for association with prostate cancer risk, low (< 7) vs. high (≥ 7) Gleason grade, localised vs. advanced stage, and mortality, in 22,936 controls and 22,992 cases. IV analysis was used in an attempt to estimate the causal effect of circulating IGF peptides on prostate cancer. Published SNPs in the IGFBP1/IGFBP3 gene region, particularly rs11977526, were strong instruments for IGF‐II and IGFBP‐3, less so for IGF‐I. Rs11977526 was associated with high (vs. low) Gleason grade (OR per IGF‐II/IGFBP‐3 level‐raising allele 1.05; 95% CI: 1.00, 1.10). Using rs11977526 as an IV we estimated the causal effect of a one SD increase in IGF‐II (∼265 ng/mL) on risk of high vs. low grade disease as 1.14 (95% CI: 1.00, 1.31). Because of the potential for pleiotropy of the genetic instruments, these findings can only causally implicate the IGF pathway in general, not any one specific biomarker.

Background BRCA1 interacting protein C-terminal helicase 1 (BRIP1) is one of the Fanconi Anaemia Complementation (FANC) group family of DNA repair proteins. Biallelic mutations in BRIP1 are responsible for FANC group J, and previous studies have also suggested that rare protein truncating variants in BRIP1 are associated with an increased risk of breast cancer. These studies have led to inclusion of BRIP1 on targeted sequencing panels for breast cancer risk prediction. Methods We evaluated a truncating variant, p.Arg798Ter (rs137852986), and 10 missense variants of BRIP1, in 48 144 cases and 43 607 controls of European origin, drawn from 41 studies participating in the Breast Cancer Association Consortium (BCAC). Additionally, we sequenced the coding regions of BRIP1 in 13 213 cases and 5242 controls from the UK, 1313 cases and 1123 controls from three population-based studies as part of the Breast Cancer Family Registry, and 1853 familial cases and 2001 controls from Australia. Results The rare truncating allele of rs137852986 was observed in 23 cases and 18 controls in Europeans in BCAC (OR 1.09, 95% CI 0.58 to 2.03, p=0.79). Truncating variants were found in the sequencing studies in 34 cases (0.21%) and 19 controls (0.23%) (combined OR 0.90, 95% CI 0.48 to 1.70, p=0.75). Conclusions These results suggest that truncating variants in BRIP1, and in particular p.Arg798Ter, are not associated with a substantial increase in breast cancer risk. Such observations have important implications for the reporting of results from breast cancer screening panels.

Abstract Identifying genetic variants with pleiotropic associations can uncover common pathways influencing multiple cancers. We took a two-stage approach to conduct genome-wide association studies for lung, ovary, breast, prostate, and colorectal cancer from the GAME-ON/GECCO Network (61,851 cases, 61,820 controls) to identify pleiotropic loci. Findings were replicated in independent association studies (55,789 cases, 330,490 controls). We identified a novel pleiotropic association at 1q22 involving breast and lung squamous cell carcinoma, with eQTL analysis showing an association with ADAM15/THBS3 gene expression in lung. We also identified a known breast cancer locus CASP8/ALS2CR12 associated with prostate cancer, a known cancer locus at CDKN2B-AS1 with different variants associated with lung adenocarcinoma and prostate cancer, and confirmed the associations of a breast BRCA2 locus with lung and serous ovarian cancer. This is the largest study to date examining pleiotropy across multiple cancer-associated loci, identifying common mechanisms of cancer development and progression. Cancer Res; 76(17); 5103–14. ©2016 AACR.

Previous genome‐wide association studies among women of European ancestry identified two independent breast cancer susceptibility loci represented by single nucleotide polymorphisms (SNPs) rs13281615 and rs11780156 at 8q24. A fine‐mapping study across 2.06 Mb (chr8:127,561,724–129,624,067, hg19) in 55,540 breast cancer cases and 51,168 controls within the Breast Cancer Association Consortium was conducted. Three additional independent association signals in women of European ancestry, represented by rs35961416 (OR = 0.95, 95% CI = 0.93–0.97, conditional p = 5.8 × 10−6), rs7815245 (OR = 0.94, 95% CI = 0.91–0.96, conditional p = 1.1 × 10−6) and rs2033101 (OR = 1.05, 95% CI = 1.02–1.07, conditional p = 1.1 × 10−4) were found. Integrative analysis using functional genomic data from the Roadmap Epigenomics, the Encyclopedia of DNA Elements project, the Cancer Genome Atlas and other public resources implied that SNPs rs7815245 in Signal 3, and rs1121948 in Signal 5 (in linkage disequilibrium with rs11780156, r2 = 0.77), were putatively functional variants for two of the five independent association signals. The results highlighted multiple 8q24 variants associated with breast cancer susceptibility in women of European ancestry.

AbstractGenetic risk scores were used as unconfounded instruments for specific lipid traits (Mendelian randomization) to assess whether circulating lipids causally influence prostate cancer risk. Data from 22,249 prostate cancer cases and 22,133 controls from 22 studies within the international PRACTICAL consortium were analyzed. Allele scores based on single nucleotide polymorphisms (SNPs) previously reported to be uniquely associated with each of low‐density lipoprotein (LDL), high‐density lipoprotein (HDL), and triglyceride (TG) levels, were first validated in an independent dataset, and then entered into logistic regression models to estimate the presence (and direction) of any causal effect of each lipid trait on prostate cancer risk. There was weak evidence for an association between the LDL genetic score and cancer grade: the odds ratio (OR) per genetically instrumented standard deviation (SD) in LDL, comparing high‐ (≥7 Gleason score) versus low‐grade (<7 Gleason score) cancers was 1.50 (95% CI: 0.92, 2.46; P = 0.11). A genetically instrumented SD increase in TGs was weakly associated with stage: the OR for advanced versus localized cancer per unit increase in genetic risk score was 1.68 (95% CI: 0.95, 3.00; P = 0.08). The rs12916‐T variant in 3‐hydroxy‐3‐methylglutaryl‐CoA reductase (HMGCR) was inversely associated with prostate cancer (OR: 0.97; 95% CI: 0.94, 1.00; P = 0.03). In conclusion, circulating lipids, instrumented by our genetic risk scores, did not appear to alter prostate cancer risk. We found weak evidence that higher LDL and TG levels increase aggressive prostate cancer risk, and that a variant in HMGCR (that mimics the LDL lowering effect of statin drugs) reduces risk. However, inferences are limited by sample size and evidence of pleiotropy.

AbstractAutomated methods are needed to facilitate high‐throughput and reproducible scoring of Ki67 and other markers in breast cancer tissue microarrays (TMAs) in large‐scale studies. To address this need, we developed an automated protocol for Ki67 scoring and evaluated its performance in studies from the Breast Cancer Association Consortium. We utilized 166 TMAs containing 16,953 tumour cores representing 9,059 breast cancer cases, from 13 studies, with information on other clinical and pathological characteristics. TMAs were stained for Ki67 using standard immunohistochemical procedures, and scanned and digitized using the Ariol system. An automated algorithm was developed for the scoring of Ki67, and scores were compared to computer assisted visual (CAV) scores in a subset of 15 TMAs in a training set. We also assessed the correlation between automated Ki67 scores and other clinical and pathological characteristics. Overall, we observed good discriminatory accuracy (AUC = 85%) and good agreement (kappa = 0.64) between the automated and CAV scoring methods in the training set. The performance of the automated method varied by TMA (kappa range= 0.37–0.87) and study (kappa range = 0.39–0.69). The automated method performed better in satisfactory cores (kappa = 0.68) than suboptimal (kappa = 0.51) cores (p‐value for comparison = 0.005); and among cores with higher total nuclei counted by the machine (4,000–4,500 cells: kappa = 0.78) than those with lower counts (50–500 cells: kappa = 0.41; p‐value = 0.010). Among the 9,059 cases in this study, the correlations between automated Ki67 and clinical and pathological characteristics were found to be in the expected directions. Our findings indicate that automated scoring of Ki67 can be an efficient method to obtain good quality data across large numbers of TMAs from multicentre studies. However, robust algorithm development and rigorous pre‐ and post‐analytical quality control procedures are necessary in order to ensure satisfactory performance.