Chenevix‐Trench, Georgia

Abstract Introduction: Genome-wide association studies over the past 10 years have identified over 100 common variants associated with breast cancer risk. The association of these polymorphisms varies according to estrogen receptor status of the tumor, with some associated primarily with ER+ or ER- disease, or both subtypes. However, ER expression in breast cancer cells is continuous. The aim of this study was to establish whether the association between polygenic risk scores and ER expression in breast cancer is continuous in nature. Methods: Data were taken from the international Breast Cancer Association Consortium (BCAC). Germline genotype data for 75 SNPs previously reported to be associated with breast cancer were available for 4,999 breast cancer cases for which ER expression data were also available. Automated ER expression on a continuous scale was measured using two image analysis algorithms: Astrogrid and Ariol. A third form of ER expression score, the manual Allred ordinal, was used. Genotype data were used to calculate polygenic risk scores using the equation: PRS = ∑ βx where β is the log-odds of the odds ratio (OR) of breast cancer associated with that single nucleotide polymorphism (SNP), and x is the number of variant alleles at a SNP locus. Five different PRSs were calculated: i) SNPs associated with both ER+ and ER- disease (PRS-A, n= 74) ii) SNPs associated with ER+ disease at P<0.01 (PRS-B, n= 65) iii) SNPs associated with ER- disease at P<0.01 (PRS-C, n= 43) iv) SNPs associated with ER+ but not ER- disease (PRS-D, n= 30) v) SNPs associated with ER- but not ER+ disease (PRS-E, n = 8) The ER-specific weights for the PRS were taken from previous OR estimates. Linear regression models adjusted for study were used to test for a relationship between the five PRSs and ER expression measured on the Ariol, Astrogrid and Allred scales. Results: There was a strong positive correlation found between PRS-B against Astrogrid (β=0.0176, P=4.06x10-4), a weak negative correlation found with PRS-C (β=-0.012, P=0.019) and a weak positive correlation for PRS-A (β=0.012, P=0.0144). PRS-D and E showed far stronger association with Astrogrid than PRS-B and C, respectively (β=0.0186 and -0.017, P=1.97x10-4 and 8.1x10-4). There was a qualitatively similar association when the PRSs were correlated with Ariol and Allred, apart from the stronger correlations with PRS-B (β=91.2 and 0.0853, P=5.52x10-4 and 1.33x10-3). Conclusion: We have shown a robust association between PRS and ER expression in the directions we predicted. This suggests that the continuous variation in ER expression is biologically plausible and linked to genetic polymorphism. The weakening of association with PRS-A and the stronger associations with PRS-D and E show that having SNPs that are associated with both ER subtypes in the same PRS weakens significance, and suggests that subtype-specific PRSs may be improved on by using SNPs exclusive to the respective subtype. Citation Format: Anamay H. Shetty, Raza Ali, Fiona Blows, Carlos Caldas, Nicholas Inard, Javier Benítez, Jenny Chang-Claude, Georgia Chenevix-Trench, Fergus Couch, Angela Cox, Peter Devilee, Douglas Easton, Montserrat Garcia-Closas, Arto Mannermaa, Heli Nevanlinna, Marjanka Schmidt, Paul Pharoah. Patient level polygenic risk scores and continuous estrogen receptor expression in breast 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 2211.

Abstract Large-scale genome-wide association studies (GWAS) have identified approximately 35 loci associated with epithelial ovarian cancer (EOC) risk. The majority of GWAS-identified disease susceptibility variants are located in noncoding regions, and causal genes underlying these associations remain largely unknown. Here, we performed a transcriptome-wide association study to search for novel genetic loci and plausible causal genes at known GWAS loci. We used RNA sequencing data (68 normal ovarian tissue samples from 68 individuals and 6,124 cross-tissue samples from 369 individuals) and high-density genotyping data from European descendants of the Genotype-Tissue Expression (GTEx V6) project to build ovarian and cross-tissue models of genetically regulated expression using elastic net methods. We evaluated 17,121 genes for their cis-predicted gene expression in relation to EOC risk using summary statistics data from GWAS of 97,898 women, including 29,396 EOC cases. With a Bonferroni-corrected significance level of P < 2.2 × 10−6, we identified 35 genes, including FZD4 at 11q14.2 (Z = 5.08, P = 3.83 × 10−7, the cross-tissue model; 1 Mb away from any GWAS-identified EOC risk variant), a potential novel locus for EOC risk. All other 34 significantly associated genes were located within 1 Mb of known GWAS-identified loci, including 23 genes at 6 loci not previously linked to EOC risk. Upon conditioning on nearby known EOC GWAS-identified variants, the associations for 31 genes disappeared and three genes remained (P < 1.47 × 10−3). These data identify one novel locus (FZD4) and 34 genes at 13 known EOC risk loci associated with EOC risk, providing new insights into EOC carcinogenesis. Significance: Transcriptomic analysis of a large cohort confirms earlier GWAS loci and reveals FZD4 as a novel locus associated with EOC risk. Cancer Res; 78(18); 5419–30. ©2018 AACR.

Abstract Introduction:Pathogenic and/or founder germline variants in the ataxia-telangiectasia mutated (ATM) gene confer an increased breast cancer (BC) risk. The protein kinase ATM plays a central role inDNA double-strand break-repair and in the activation of downstream targets such as p53 and BRCA1. We sought to define the repertoire of somatic genetic alterations of BCs from patients with pathogenic germline ATM mutations and whether somatic loss of heterozygosity (LOH) of ATM would be present in these cancers. Methods: 21 BCs from ATM germline mutation carriers were microdissected. Tumor and normal DNA samples were subjected to whole-exome sequencing (WES, n=12) or massively parallel sequencing targeting all coding regions and selected intronic and regulatory regions of 410 key cancer genes (n=9). Somatic mutations, copy number alterations, cancer cell fractions, large-scale state transitions (LSTs) and mutational signatures were defined using state-of-the-art bioinformatics algorithms. ABSOLUTE and FACETS were employed to assess LOH of the wild-type allele of ATM. Results: Of the patients included in this study, 71%, 24% and 5% of cases harbored ATM missense (all but one p.V2424G), frame-shift and nonsense germline mutations, respectively. All tumors were ER-positive and four (19%) were HER2-positive. The median age of the patients was 46 years (32–79 years). Our analyses revealed biallelic inactivation of ATM through LOH of the wild-type allele in 16 of 21 cases (76%), and second somatic ATM mutations were not found. The median number of non-synonymous somatic mutations was 38 (range 15-113) and 2 (range 0-8)in tumors subjected to WES and targeted sequencing, respectively. The repertoire of somatic genetic alterations of ATM-associated BCs was found to be heterogeneous, including clonal PIK3CA mutations (24%), GATA3 mutations (19%), FANCI amplifications (19%) and CCND1 amplifications (14%). Importantly, however, no somatic mutations affecting TP53 were found. Analysis of the WES data revealed that 5 (42%) ATM-associated BCs displayed high LST scores, all of which harbored bi-allelic ATM inactivation. In contrast to BRCA1- and BRCA2-associated BCs, which frequently display the mutational signature 3 associated with defective homologous recombination DNA repair, the ATM-associated BCs studied displayed the ageing mutational signature (i.e. signature 1). Comparison of the mutational profiles of the ATM--associated BCs subjected to WES (n=12) with those of BRCA1- (n=11) and BRCA2-associated (n=10) BCs from The Cancer Genome Atlas revealed that TP53 was more frequently mutated in BCs from BRCA1 germline mutation carriers (0% vs 72%, P<0.001), while no differences with BRCA2-associated BCs were found. Conclusion: ATM-associated BCs frequently display bi-allelic ATM inactivation through LOH of the wild-type allele and a subset of these cases displayed high levels of LSTs. These findings suggest that at least in a subset of ATM-associated BCs, biallelic inactivation of ATM rather than a dominant negative effect of the germline mutation may be the mechanism of inactivation of this tumor suppressor gene. The repertoire of somatic genetic alterations of ATM-associated BCs is heterogeneous, with a noticeable lack of TP53 somatic mutations. Citation Format: Weigelt B, Bi R, Kumar R, James PA, Thorne H, Couch FJ, Eccles DM, Blows F, Geyer FC, Li A, Selenica P, Lim RS, Blecua P, Shen R, Wen H, Robson ME, Reis-Filho JS, Chenevix-Trench G. The landscape of somatic genetic alterations in breast cancers from ATM germline mutation carriers [abstract]. In: Proceedings of the 2017 San Antonio Breast Cancer Symposium; 2017 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2018;78(4 Suppl):Abstract nr PD1-15.