Beesley, Jonathan

Abstract Genome-wide association studies (GWAS) have identified several susceptibility regions that confer risk of two or more cancers. A region at 19p13 is associated with both breast and ovarian cancer risk in the general population as well as in BRCA1/BRCA2 mutation carriers. Using the iCOGs custom genotyping array, we analyzed 438 SNPs in this region in 7,435 ER-negative breast cancer cases and 42,599 controls from BCAC, and 15,347 ovarian cancer cases and 30,845 controls from OCAC and 15,252 BRCA1 mutation carriers. We identified two significant peaks of association. Peak 1, consisting of 13 highly correlated SNPs (r2>0.95), spans a 19.4kb region harboring the BABAM1, ANKLE1 and ABHD8 genes. Peak 1 SNPs are associated with serous ovarian cancer (top hit rs4808075, P = 9.2×10-20), ER-negative breast cancer (rs61494113, P = 1.1×10-13) and breast cancer in BRCA1 mutation carriers (rs61494113, P = 7.7×10-16). Peak 2 spans over 100kb and contains 7 SNPs associated with breast cancer in BRCA1 mutation carriers only (rs4464206, P = 8.9×10-8). Multiple assays were used to establish the putative functional roles of risk SNPs and candidate susceptibility genes in this region. Expression quantitative trait locus (eQTL) analyses identified significant genotype-gene expression correlations between genotype at the Peak 1 SNP rs480816 and ABHD8 expression in ovarian tumors (P = 3.0×10-5) and normal breast tissues (P = 2.8×10-3). The Peak 2 SNP, rs7246243, was significantly associated with BABAM1 expression (P = 5×10-3) in normal breast tissue only. Overexpression of BABAM1 in normal ovarian and breast cells did not modulate neoplastic phenotypes, whereas overexpression of ABHD8 induced a significant reduction in migration in both cell types in vitro. Differentially expressed genes in ABHD8 models, identified by RNAseq, were associated with cellular movement (P<1.82×10-44). Chromosome conformation capture (3C) identified an interaction between the ABHD8 promoter and a region containing 4 risk SNPs (rs4808075, rs10419397, rs56069439 and rs4808076) in breast and ovarian normal and cancer cells but only rs56069439 coincides with epigenetic enhancer marks in breast and ovarian cells. CRISPR/Cas9 genome editing was used to delete a 57bp region containing rs56069439 and expression of BABAM1, ANKLE1 and ABHD8 were measured. Taken together, these data suggest a common functional mechanism underlying risk of breast and ovarian cancer at this locus, mediated by rs56069439 variants differentially regulating expression of ABHD8, a gene involved in cell motility. Citation Format: Kate Lawrenson, Siddhartha Kar, Karoline Kuckenbaeker, Stacey Edwards, Qiyuan Li, Jonathan Tyrer, Jonathan Beesley, Kyriaki Michailidou, Susan Ramus, Alison Dunning, Fergus Couch, Alvaro Monteiro, Jacques Simard, Matthew Freedman, Douglas Easton, Georgia Chenevix-Trench, Paul Pharoah, Antonis Antoniou, Simon Gayther, OCAC, BCAC and CIMBA. Common functional mechanisms underlying pleiotropy at the 19p13.1 breast and ovarian cancer cusceptibility locus. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 2783. doi:10.1158/1538-7445.AM2015-2783

Abstract Primary and acquired chemotherapy resistance remains a major challenge in the treatment of ovarian cancer. We hypothesize that heritable polymorphisms might be associated with clinical outcome. Through the Ovarian Cancer Association Consortium, we have collected detailed clinical data from a large number of ovarian cancer patients. We have conducted association analyses of progression-free survival (PFS) and overall survival (OS) for 2,901 European women with ovarian carcinoma for ∼8 million imputed SNPs with minor allele frequency ≥ 0.02 and good imputation quality (imputation r2 ≥ 0.3) from the iCOGS experiment (Pharoah et al Nat Genet. 2013). Cox regression models adjusted for known prognostic factors were used to assess the log-additive effects of SNPs on PFS and OS in 2,825 patients who underwent cytoreductive surgery and received any chemotherapy (‘all chemo’) and in the subset of 1,098 patients who underwent cytoreductive surgery followed by ≥ 4 cycles of carboplatin AUC 5 or 6 and paclitaxel 175 or 135 mg/m2 intravenously (‘standard chemo’). We identified five SNPs at four loci that were significantly associated with outcome (p ≤ 10-7) in at least one of four analyses (Table). Analysis of expression data for genes within a 1Mb region flanking each lead SNP indicates significant associations between nearby genes and ovarian cancer survival. Of particular interest is the heat shock protein crystalline αB (CRYAB) at 11q23.1. We have also identified significant associations with KIAA1598 and RAVER2. Further analysis of these loci in a larger genotyped sample is necessary to definitively identify SNP(s) associated with outcome. Our study adds to the growing evidence that genetic variants are associated with ovarian cancer prognosis. Table: SNPs associated with ovarian cancer outcome and gene expression Overall Survival (P)*Progression-Free Survival (P)*Lead SNPChromosome LocusImputation R²Expression AnalysisAll chemoStandard chemoAll chemoStandard chemors1429768671p31.30.742.14E-0512.69E-044.68E-055.11E-073.64E-09rs559375217q32.30.600.198.58E-082.32E-052.80E-031.85E-03rs7885346210q26.110.602.23E-0724.57E-036.57E-081.65E-016.85E-03rs1089121911q23.10.841.53E-0531.51E-059.50E-081.50E-042.32E-04rs3579103211q23.10.811.53E-0531.66E-055.88E-086.46E-051.87E-04*P-values are adjusted for known confounders of ovarian cancer outcomes; FIGO stage, grade, histology, residual disease, and additionally for age at diagnosis for OS.1 RAVER2; 2 KIAA1598; 3 CRYAB Citation Format: Sharon E. Johnatty, Jonathan Tyrer, Jonathan Beesley, Bo Gao, Yi Lu, Stuart MacGregor, Anna deFazio, Paul Pharoah, Ellen Goode, Georgia Chenevix-Trench. Identification of genetic loci associated with ovarian cancer prognosis. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 3286. doi:10.1158/1538-7445.AM2014-3286

Abstract Background: Genome-wide association studies (GWAS) identified variants at 19p13.1 and ZNF365 (10q21.2) as risk factors for breast cancer among BRCA1 and BRCA2 mutation carriers, respectively. We explored associations with ovarian cancer and with breast cancer by tumor histopathology for these variants in mutation carriers from the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). Methods: Genotyping data for 12,599 BRCA1 and 7,132 BRCA2 mutation carriers from 40 studies were combined. Results: We confirmed associations between rs8170 at 19p13.1 and breast cancer risk for BRCA1 mutation carriers [HR, 1.17; 95% confidence interval (CI), 1.07–1.27; P = 7.42 × 10−4] and between rs16917302 at ZNF365 (HR, 0.84; 95% CI, 0.73–0.97; P = 0.017) but not rs311499 at 20q13.3 (HR, 1.11; 95% CI, 0.94–1.31; P = 0.22) and breast cancer risk for BRCA2 mutation carriers. Analyses based on tumor histopathology showed that 19p13 variants were predominantly associated with estrogen receptor (ER)-negative breast cancer for both BRCA1 and BRCA2 mutation carriers, whereas rs16917302 at ZNF365 was mainly associated with ER-positive breast cancer for both BRCA1 and BRCA2 mutation carriers. We also found for the first time that rs67397200 at 19p13.1 was associated with an increased risk of ovarian cancer for BRCA1 (HR, 1.16; 95% CI, 1.05–1.29; P = 3.8 × 10−4) and BRCA2 mutation carriers (HR, 1.30; 95% CI, 1.10–1.52; P = 1.8 × 10−3). Conclusions: 19p13.1 and ZNF365 are susceptibility loci for ovarian cancer and ER subtypes of breast cancer among BRCA1 and BRCA2 mutation carriers. Impact: These findings can lead to an improved understanding of tumor development and may prove useful for breast and ovarian cancer risk prediction for BRCA1 and BRCA2 mutation carriers. Cancer Epidemiol Biomarkers Prev; 21(4); 645–57. ©2012 AACR.

Abstract The 19p13.1 breast cancer susceptibility locus is a modifier of breast cancer risk in BRCA1 mutation carriers and is also associated with the risk of ovarian cancer. Here, we investigated 19p13.1 variation and risk of breast cancer subtypes, defined by estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor-2 (HER2) status, using 48,869 breast cancer cases and 49,787 controls from the Breast Cancer Association Consortium (BCAC). Variants from 19p13.1 were not associated with breast cancer overall or with ER-positive breast cancer but were significantly associated with ER-negative breast cancer risk [rs8170 OR, 1.10; 95% confidence interval (CI), 1.05–1.15; P = 3.49 × 10−5] and triple-negative (ER-, PR-, and HER2-negative) breast cancer (rs8170: OR, 1.22; 95% CI, 1.13–1.31; P = 2.22 × 10−7). However, rs8170 was no longer associated with ER-negative breast cancer risk when triple-negative cases were excluded (OR, 0.98; 95% CI, 0.89–1.07; P = 0.62). In addition, a combined analysis of triple-negative cases from BCAC and the Triple Negative Breast Cancer Consortium (TNBCC; N = 3,566) identified a genome-wide significant association between rs8170 and triple-negative breast cancer risk (OR, 1.25; 95% CI, 1.18–1.33; P = 3.31 × 10−13]. Thus, 19p13.1 is the first triple-negative–specific breast cancer risk locus and the first locus specific to a histologic subtype defined by ER, PR, and HER2 to be identified. These findings provide convincing evidence that genetic susceptibility to breast cancer varies by tumor subtype and that triple-negative tumors and other subtypes likely arise through distinct etiologic pathways. Cancer Res; 72(7); 1795–803. ©2012 AACR.

Abstract Triple-negative breast cancers are an aggressive subtype of breast cancer with poor survival, but there remains little known about the etiologic factors that promote its initiation and development. Commonly inherited breast cancer risk factors identified through genome-wide association studies display heterogeneity of effect among breast cancer subtypes as defined by the status of estrogen and progesterone receptors. In the Triple Negative Breast Cancer Consortium (TNBCC), 22 common breast cancer susceptibility variants were investigated in 2,980 Caucasian women with triple-negative breast cancer and 4,978 healthy controls. We identified six single-nucleotide polymorphisms, including rs2046210 (ESR1), rs12662670 (ESR1), rs3803662 (TOX3), rs999737 (RAD51L1), rs8170 (19p13.1), and rs8100241 (19p13.1), significantly associated with the risk of triple-negative breast cancer. Together, our results provide convincing evidence of genetic susceptibility for triple-negative breast cancer. Cancer Res; 71(19); 6240–9. ©2011 AACR.

Abstract Background: Invasive ovarian cancer is a significant cause of gynecologic cancer mortality. Methods: We examined whether this mortality was associated with inherited variation in approximately 170 candidate genes/regions [993 single-nucleotide polymorphisms (SNPs)] in a multistage analysis based initially on 312 Mayo Clinic cases (172 deaths). Additional analyses used The Cancer Genome Atlas (TCGA; 127 cases, 62 deaths). For the most compelling gene, we immunostained Mayo Clinic tissue microarrays (TMA, 326 cases) and conducted consortium-based SNP replication analysis (2,560 cases, 1,046 deaths). Results: The strongest initial mortality association was in HGF (hepatocyte growth factor) at rs1800793 (HR = 1.7, 95% CI = 1.3–2.2, P = 2.0 × 10−5) and with overall variation in HGF (gene-level test, P = 3.7 × 10−4). Analysis of TCGA data revealed consistent associations [e.g., rs5745709 (r2 = 0.96 with rs1800793): TCGA HR = 2.4, CI = 1.4–4.1, P = 2.2 × 10−3; Mayo Clinic + TCGA HR = 1.6, CI = 1.3–1.9, P = 7.0 × 10−5] and suggested genotype correlation with reduced HGF mRNA levels (P = 0.01). In Mayo Clinic TMAs, protein levels of HGF, its receptor MET (C-MET), and phospho-MET were not associated with genotype and did not serve as an intermediate phenotype; however, phospho-MET was associated with reduced mortality (P = 0.01) likely due to higher expression in early-stage disease. In eight additional ovarian cancer case series, HGF rs5745709 was not associated with mortality (HR = 1.0, CI = 0.9–1.1, P = 0.87). Conclusions: We conclude that although HGF signaling is critical to migration, invasion, and apoptosis, it is unlikely that HGF genetic variation plays a major role in ovarian cancer mortality. Furthermore, any minor role is not related to genetically-determined expression. Impact: Our study shows the utility of multiple data types and multiple data sets in observational studies. Cancer Epidemiol Biomarkers Prev; 20(8); 1638–48. ©2011 AACR.

Abstract As the leading cause of mortality among gynecologic cancers, it is critical to increase understanding of ovarian cancer pathogenesis, which remains largely unknown. Several xenobiotic metabolizing enzymes are responsible for processing and inactivating carcinogens that could initiate ovarian carcinogenesis. Here, we hypothesize that single-nucleotide polymorphisms (SNPs) in genes encoding xenobiotic metabolizing enzymes may be associated with differential risk of ovarian cancer. Cases with epithelial ovarian cancer (N=930) and controls (N=1,037), frequency-matched on age and geographic residence, were enrolled at Mayo Clinic and Duke University, and genotyped at 163 tagSNPs in EPHX1, SULT1E1, ADH5, ADH4, ADH6, ADH1A, ADH1B, ADH1C, NQO2, NAT1, NAT2, GSTP1, CYP19A1, CYP1A1, CYP1A2, and NQO1 using an Illumina GoldenGate™ BeadArray array. Logistic regression was used to estimate odds ratios (ORs) and adjust for population structure and other potential confounders, including oral contraceptive use, parity, and family history. SNP-specific statistical analyses identified several noteworthy results including an association with increased risk of ovarian cancer with a SNP (rs2917666) in NAD(P)H dehydrogenase, quinone 1 (NQO1) (per-allele OR 1.16, 95% confidence interval (CI) 1.01-1.33, p=0.04). The SNP rs1051740 in epoxide hydrolase 1, microsomal (xenobiotic) (EPHX1) was suggestive of increased risk of invasive serous disease (TT v CT OR 0.85, 95% CI 0.66-1.09, TT v CC 1.29, 95% CI 0.89-1.88, p=0.09). In an effort to replicate these observations, a subset of 11 SNPs with a p-value < 0.10 were then genotyped in 727 invasive cancer cases and 1,065 controls from an Australian ovarian cancer study using a Sequenom iPlex assay. Results of combined analyses confirmed an association between EPHX1 rs1051740 and increased invasive serous ovarian cancer risk (TT v CT OR 0.99, 95% CI 0.77-1.29, TT v CC OR 1.61, 95% CI 1.23-2.11, p=0.001; recessive OR 1.61, 95% CI 1.25-2.09, p=0.001). This gene is involved in the activation and detoxification of tobacco-derived carcinogens, and smoking has been associated with the mucinous histologic subtype of ovarian cancer. Analyses by histological subtype and gene-environment interaction analyses considering genotype in conjunction with tobacco use and alcohol intake are underway. Due to the important role of EPHX1 in metabolizing pro-carcinogens and the observed association between rs1051740 and risk of invasive serous disease, an examination of this SNP in a larger ovarian cancer consortium is warranted. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 914.