Citations

AbstractPrevious transcriptome‐wide association studies (TWAS) have identified breast cancer risk genes by integrating data from expression quantitative loci and genome‐wide association studies (GWAS), but analyses of breast cancer subtype‐specific associations have been limited. In this study, we conducted a TWAS using gene expression data from GTEx and summary statistics from the hitherto largest GWAS meta‐analysis conducted for breast cancer overall, and by estrogen receptor subtypes (ER+ and ER−). We further compared associations with ER+ and ER− subtypes, using a case‐only TWAS approach. We also conducted multigene conditional analyses in regions with multiple TWAS associations. Two genes, STXBP4 and HIST2H2BA, were specifically associated with ER+ but not with ER– breast cancer. We further identified 30 TWAS‐significant genes associated with overall breast cancer risk, including four that were not identified in previous studies. Conditional analyses identified single independent breast‐cancer gene in three of six regions harboring multiple TWAS‐significant genes. Our study provides new information on breast cancer genetics and biology, particularly about genomic differences between ER+ and ER− breast cancer.

Abstract Purpose: To determine if the degree of estrogen suppression with aromatase inhibitors (AI: anastrozole, exemestane, letrozole) is associated with efficacy in early-stage breast cancer, and to examine for differences in the mechanism of action between the three AIs. Experimental Design: Matched case–control studies [247 matched sets from MA.27 (anastrozole vs. exemestane) and PreFace (letrozole) trials] were undertaken to assess whether estrone (E1) or estradiol (E2) concentrations after 6 months of adjuvant therapy were associated with risk of an early breast cancer event (EBCE). Preclinical laboratory studies included luciferase activity, cell proliferation, radio-labeled ligand estrogen receptor binding, surface plasmon resonance ligand receptor binding, and nuclear magnetic resonance assays. Results: Women with E1 ≥1.3 pg/mL and E2 ≥0.5 pg/mL after 6 months of AI treatment had a 2.2-fold increase in risk (P = 0.0005) of an EBCE, and in the anastrozole subgroup, the increase in risk of an EBCE was 3.0-fold (P = 0.001). Preclinical laboratory studies examined mechanisms of action in addition to aromatase inhibition and showed that only anastrozole could directly bind to estrogen receptor α (ERα), activate estrogen response element-dependent transcription, and stimulate growth of an aromatase-deficient CYP19A1−/− T47D breast cancer cell line. Conclusions: This matched case–control clinical study revealed that levels of estrone and estradiol above identified thresholds after 6 months of adjuvant anastrozole treatment were associated with increased risk of an EBCE. Preclinical laboratory studies revealed that anastrozole, but not exemestane or letrozole, is a ligand for ERα. These findings represent potential steps towards individualized anastrozole therapy.

1048 Background: Germline BRCA1/2 ( gBRCA) mutations (mt) are some of the few actionable alterations in mBC patients. The PARP inhibitors olaparib/talazoparib are more effective than chemotherapy (ctx) in patients with a gBRCAmt and HER2 negative(-) mBC. In mBC patients the TNT-study suggested a better progression-free survival (PFS) for g BRCA-mt compared to patients with a gBRCA1/2 wildtype (wt) when treated with platinum and not with a taxane. Otherwise little is known about the prognostic effect of g BRCA1/2mt in mBC patients. Methods: PRAEGNANT (NCT02338167) is a prospective mBC registry with a focus on molecular biomarkers. Patients were eligible for this analysis if their mBC was HER2- and treated with ctx for the first time (referred to as first-line ctx). Hormone receptor (HR) positive patients had to have all hormone therapies exhausted. Mutation frequencies and their association with patient and tumor characteristics were analyzed. Multivariable Cox regression models were built with commonly established prognostic factors and g BRCA mutation status as predictors of PFS and overall survival (OS) from first-line ctx. Results: Out of 2932 PRAEGNANT patients, 576 were HER2- and received first-line ctx. Of those 529 patients with g BRCA genotype results and follow up information could be analyzed. 24 patients (4.5%) had a g BRCAmt (11 BRCA1, 13 BRCA2). Mutation rate in HR positive patients was 3.9% (17/432) and 7.2% (7/97) in HR negative patients. Most patients received ctx either as the first treatment in the metastatic setting or after one line of hormone therapy (n=382; 72.2%). Multivariable Cox regression models showed an adjusted hazard ratio for gBRCAmt vs. gBRCAwt patients of 0.70 (95% CI: 0.43-1.15) for PFS and of 0.41 (95% CI: 0.18-0.93) for OS. Most frequent ctx treatments were taxane (52%) or capecitabine based (21%). Additionally, the prevalence of somatic BRCA1/2 mutations in this population will be presented. Conclusions: In this HER2- mBC population under ctx g BRCA mutation rates were within the expected range of about 5%. Within the analyzed population patients with a g BRCA mutation seemed to have a better OS than patients without a mutation. PFS results pointed in the same direction without statistical significance. However, with only 24 mutations replication of these results in additional cohorts is warranted. Clinical trial information: NCT02338167.

Abstract DNA methylation is instrumental for gene regulation. Global changes in the epigenetic landscape have been recognized as a hallmark of cancer. However, the role of DNA methylation in epithelial ovarian cancer (EOC) remains unclear. In this study, high-density genetic and DNA methylation data in white blood cells from the Framingham Heart Study (N = 1,595) were used to build genetic models to predict DNA methylation levels. These prediction models were then applied to the summary statistics of a genome-wide association study (GWAS) of ovarian cancer including 22,406 EOC cases and 40,941 controls to investigate genetically predicted DNA methylation levels in association with EOC risk. Among 62,938 CpG sites investigated, genetically predicted methylation levels at 89 CpG were significantly associated with EOC risk at a Bonferroni-corrected threshold of P < 7.94 × 10−7. Of them, 87 were located at GWAS-identified EOC susceptibility regions and two resided in a genomic region not previously reported to be associated with EOC risk. Integrative analyses of genetic, methylation, and gene expression data identified consistent directions of associations across 12 CpG, five genes, and EOC risk, suggesting that methylation at these 12 CpG may influence EOC risk by regulating expression of these five genes, namely MAPT, HOXB3, ABHD8, ARHGAP27, and SKAP1. We identified novel DNA methylation markers associated with EOC risk and propose that methylation at multiple CpG may affect EOC risk via regulation of gene expression. Significance: Identification of novel DNA methylation markers associated with EOC risk suggests that methylation at multiple CpG may affect EOC risk through regulation of gene expression.

Thymidylate synthase (TYMS) is a crucial enzyme for DNA synthesis. TYMS expression is regulated by its antisense mRNA, ENOSF1. Disrupted regulation may promote uncontrolled DNA synthesis and tumor growth. We sought to replicate our previously reported association between rs495139 in the TYMS-ENOSF1 3′ gene region and increased risk of mucinous ovarian carcinoma (MOC) in an independent sample. Genotypes from 24,351 controls to 15,000 women with invasive OC, including 665 MOC, were available. We estimated per-allele odds ratios (OR) and 95% confidence intervals (CI) using unconditional logistic regression, and meta-analysis when combining these data with our previous report. The association between rs495139 and MOC was not significant in the independent sample (OR = 1.09; 95% CI = 0.97–1.22; p = 0.15; N = 665 cases). Meta-analysis suggested a weak association (OR = 1.13; 95% CI = 1.03–1.24; p = 0.01; N = 1019 cases). No significant association with risk of other OC histologic types was observed (p = 0.05 for tumor heterogeneity). In expression quantitative trait locus (eQTL) analysis, the rs495139 allele was positively associated with ENOSF1 mRNA expression in normal tissues of the gastrointestinal system, particularly esophageal mucosa (r = 0.51, p = 1.7 × 10−28), and nonsignificantly in five MOC tumors. The association results, along with inconclusive tumor eQTL findings, suggest that a true effect of rs495139 might be small.

Abstract Ovarian cancer is the most lethal gynecologic malignancy and the sixth most common cause of death from all cancer in women, with an estimated 140,000 deaths per year worldwide. Endometrial cancer (cancer of the uterine lining) is the most commonly diagnosed gynecologic cancer, ranking sixth in incident cancer in women. Ovarian and endometrial cancer share many epidemiologic, histopathologic, and tumor genetic characteristics. Meta-analyses of genome-wide association study (GWAS) datasets across etiologically related diseases have successfully been used to increase statistical power and identify novel genetic risk regions. We hypothesised that joint meta-analysis of ovarian and endometrial cancer GWAS datasets would identify novel genetic loci for both cancers. Following quality control, summary statistics for >11 million genetic variants were available from the largest genome-wide association studies performed by the Ovarian Cancer Association Consortium (OCAC, Phelan et al., Nature Genetics 2017) and the Endometrial Cancer Association Consortium (ECAC, unpublished). A total of 35,312 cancer cases from all histologic subtypes (22,406 ovarian and 12,906 endometrial cancer cases) and 149,920 controls were included in the analysis. Summary statistics were combined using an inverse-variance, fixed-effects model in METAL and identified ten loci at genome-wide significance (P < 5 x 10-8), of which one locus at 14q23.3 (rs4072776 combined OR 0.94; 95% CI 0.92-0.96; P = 8.0 x 10-9) had not been previously associated with the risk of either cancer. Integration with chromatin conformation data (RNAPolII HiChIP) from the Ishikawa endometrial cancer cell line suggests interaction between the risk signal and the promoter of FUT8. FUT8 encodes a fucosyltransferase and its expression has been implicated in breast, stomach, and melanoma and lung cancer. Interestingly, FUT8 is a downstream target regulated by the loss of PAX2 and mutated p53, which is the earliest known molecular aberration in the progression of the fallopian tissue epithelial to serous ovarian cancer. Further functional evaluation of the 14q23.3 risk region is required to determine the regulatory effect of genetic risk variants on FUT8. Future meta-analyses of datasets, stratifying results by histologic subtypes, will be undertaken. Citation Format: Tracy A. O'Mara, Dylan M. Glubb, Daniel D. Buchanan, Diether Lambrechts, Per Hall, Emma Tham, Jone Trovik, Ellen L. Goode, Peter Fasching, Thilo Dörk, Rodney J. Scott, Paul L. Auer, Roger L. Milne, Graham G. Giles, John Perry, Immaculata de Vivo, Ian Tomlinson, Douglas F. Easton, Deborah J. Thompson, Amanda B. Spurdle, E2C2, BCAC, OCAC, ECAC. Joint genome-wide association study of endometrial cancer and ovarian cancer identifies a novel genetic risk region at 14q23.3 [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 230.

Abstract Tubo-ovarian high-grade serous cancer (HGSC) is the most common histotype of epithelial ovarian cancer (70%-80%) and the leading cause of death from gynecological malignancy. Novel monoclonal antibody and adoptive cell transfer immunotherapies have found success in cancer in the past five years; however success in HGSC treatment has been elusive thus far. We have recently shown that CD8+ T-lymphocytes infiltrating the tumor proper (TILs) associate in a dose response manner with longer patient survival time regardless of age, stage or residual disease. As understanding host factors associated with this clinically relevant tumor phenotype will contribute to the development of targeted therapies, we sought to identify germline genetic factors associated with CD8+ TIL levels in HGSC. We performed a germline genome-wide association study (GWAS) of CD8+ TILs using 1204 HGSC patients of European ancestry from twelve studies with centralized genotyping, immunohistochemistry, and harmonization of clinical data. Germline genotypes were assayed using the Illumina Infinium OncoArray Beadchip and were imputed to over 11.4 million variants based on the 1000 Genomes Project reference panel. CD8 immunohistochemistry was performed on tissue microarrays using the Leica Bond Rx stainer and scored into four levels of CD8+ TILs per x400 magnification of a 0.55-mm diameter field for each tumor hotspot: none (no CD8+ TILs), low (1-2 CD8+ TILs), moderate (3-19 CD8+ TILs), high (20 or more CD8+ TILs). GWAS used multinomial regression modelling of these four CD8+ TIL levels adjusted for study site, age, and European ancestry principal components. On the average, CD8+ TIL levels were 18%, 19%, 44% and 18% for none, low, moderate and high levels, respectively. As expected, CD8+ TIL levels associated with overall survival time (p < 2.6x10-8), and CD8+ TIL levels and age at diagnosis were similarly distributed across study sites. We identified a region on chromosome 22 with multiple, correlated variants (MAF > 20%) associated with CD8+ TIL levels at p-value < 2.0x10-5. The most statistically significant variant was an indel (rs557925408) located in the an intron of the MYO18B gene (Myosin XVIIIB, p-value < 8.7x10-6). Each copy of this variant associated with a reduced extent of CD8+ TIL infiltration (OR low v none = 0.92; OR moderate v none = 0.61). An independent set of 1200 OncoArray genotyped HGSC cases are undergoing CD8+ TIL scoring in December 2017 and will be included in final genome-wide analysis. This will provide improved statistical power and enable consideration of additional genetic models (e.g., linear trend tests). Comprehensive integrative analysis of germline, tumor, and clinical features provides a model for clinical molecular epidemiology studies and will be key to increasing our understanding of this important HGSC immunophenotype. Citation Format: Yanina Natanzon, Martin Köbel, Stacey J. Winham, Sebastian M. Armasu, Bryan M. McCauley, Robert A. Vierkant, Julie M. Cunningham, David Bowtell, Ian G. Campbell, Jenny Chang-Claude, Anna deFazio, Peter A. Fasching, Mark T. Goodman, Beth Y. Karlan, Francesmary Modugno, Kirsten B. Moysich, Roberta B. Ness, Weiva Sieh, Paul D. Pharoah, Susan J. Ramus, Ellen L. Goode. Tumor-infiltrating CD8-positive T-lymphocytes in tubo-ovarian high-grade serous cancer are associated with multiple germline variants in 22q12.1 in a genome-wide association analysis [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 225.

Abstract Background: Radiologically dense area (DA) on a mammogram represents the stromal and epithelial tissue in the breast, while non-dense area (NDA) represents adipose tissue. Quantitative measurements of DA, NDA as well as the proportion of dense area in the breast (percent density, PD) adjusted for age and body mass index (BMI) have all been associated with breast cancer (BC) risk. Previous genome-wide association studies (GWAS) have identified genetic loci that show association with both mammographic density (MD) phenotypes and BC, but little is known about the genome-wide and local genetic correlation between MD phenotypes and BC. Methods: We estimated the genome-wide and local genetic correlation between MD phenotypes adjusted for age and BMI and BC, leveraging GWAS summary statistics for MD phenotypes based on 12,866 women (MODE consortium) and for BC based on 122,977 cases and 105,974 controls (BCAC consortium). We also conducted partitioned heritability analysis to identify specific functional annotations that contribute to MD heritability. Results: We observed a negative genetic correlation between DA and NDA (-0.62, p=0.002). Both DA and PD showed positive genetic correlations with overall BC (DA: rg=0.24, p=0.004; PD: rg=0.27, p=0.0001), estrogen receptor-positive (ER+) BC (DA: rg=0.23, p=0.005; PD: rg=0.23, p=0.001), and ER- BC (DA: rg=0.26, p=0.01; PD: rg=0.32, p=0.0006). In contrast, NDA showed a nonsignificant negative genetic correlation with BC (rg=-0.11, p=0.21). Local genetic correlation analyses between MD phenotypes and BC reflected known findings including positive correlation at 10q21 region between both PD and DA and BC (both p<5x10-7), negative correlation at 22q13.1 between DA and BC (p=4.6x10-6), and positive correlation at 8p11.23 between NDA and BC (p=2.8x10-6). Two known BC regions that have not previously been associated with MD showed suggestive genetic correlations between BC and DA: 1q21.2 (negative genetic correlation, p=1.3x10-4) and 12q21.31 (positive genetic correlation, p=6.7x10-4). In exploratory analyses estimating the genome-wide genetic correlation between MD and hormone-related phenotypes in UK Biobank (n=337,199), we observed nominally significant positive genetic correlations between PD and benign mammary dysplasia (rg=0.45, p=0.01), NDA and endometriosis (rg=0.45, p=0.02), and negative genetic correlation between NDA and menstrual cycle length (rg=-0.30, p=0.05). We identified multiple functional annotations that contribute to MD heritability including histone mark H3K27ac that marks active enhancers (~2,5-fold enrichment in heritability for both DA and NDA) and DNAsel Hypersensitive Sites (DHS, 9.1-fold enrichment in heritability for DA). Conclusion: Our results provide evidence of shared genetic basis between mammographic density and breast cancer. Citation Format: Hongjie Chen, Julie Douglas, William Gordon, Deborah Thompson, Fergus Couch, Peter Fasching, Per Hall, John Hopper, Gertraud Maskarenic, Marina Pollan, Jennifer Stone, Rulla Tamimi, Sara Lindstroem. Evidence for a shared genetic basis between breast cancer and mammographic density phenotypes [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 224.

AbstractWe aimed to validate the prognostic association of p16 expression in ovarian high‐grade serous carcinomas (HGSC) and to explore it in other ovarian carcinoma histotypes. p16 protein expression was assessed by clinical‐grade immunohistochemistry in 6525 ovarian carcinomas including 4334 HGSC using tissue microarrays from 24 studies participating in the Ovarian Tumor Tissue Analysis consortium. p16 expression patterns were interpreted as abnormal (either overexpression referred to as block expression or absence) or normal (heterogeneous). CDKN2A (which encodes p16) mRNA expression was also analyzed in a subset (n = 2280) mostly representing HGSC (n = 2010). Association of p16 expression with overall survival (OS) was determined within histotypes as was CDKN2A expression for HGSC only. p16 block expression was most frequent in HGSC (56%) but neither protein nor mRNA expression was associated with OS. However, relative to heterogeneous expression, block expression was associated with shorter OS in endometriosis‐associated carcinomas, clear cell [hazard ratio (HR): 2.02, 95% confidence (CI) 1.47–2.77, p < 0.001] and endometrioid (HR: 1.88, 95% CI 1.30–2.75, p = 0.004), while absence was associated with shorter OS in low‐grade serous carcinomas (HR: 2.95, 95% CI 1.61–5.38, p = 0.001). Absence was most frequent in mucinous carcinoma (50%), and was not associated with OS in this histotype. The prognostic value of p16 expression is histotype‐specific and pattern dependent. We provide definitive evidence against an association of p16 expression with survival in ovarian HGSC as previously suggested. Block expression of p16 in clear cell and endometrioid carcinoma should be further validated as a prognostic marker, and absence in low‐grade serous carcinoma justifies CDK4 inhibition.