Citations

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.

As a follow‐up to genome‐wide association analysis of common variants associated with ovarian carcinoma (cancer), our study considers seven well‐known ovarian cancer risk factors and their interactions with 28 genome‐wide significant common genetic variants. The interaction analyses were based on data from 9971 ovarian cancer cases and 15,566 controls from 17 case–control studies. Likelihood ratio and Wald tests for multiplicative interaction and for relative excess risk due to additive interaction were used. The top multiplicative interaction was noted between oral contraceptive pill (OCP) use (ever vs. never) and rs13255292 (p value = 3.48 × 10−4). Among women with the TT genotype for this variant, the odds ratio for OCP use was 0.53 (95% CI = 0.46–0.60) compared to 0.71 (95%CI = 0.66–0.77) for women with the CC genotype. When stratified by duration of OCP use, women with 1–5 years of OCP use exhibited differential protective benefit across genotypes. However, no interaction on either the multiplicative or additive scale was found to be statistically significant after multiple testing correction. The results suggest that OCP use may offer increased benefit for women who are carriers of the T allele in rs13255292. On the other hand, for women carrying the C allele in this variant, longer (5+ years) use of OCP may reduce the impact of carrying the risk allele of this SNP. Replication of this finding is needed. The study presents a comprehensive analytic framework for conducting gene–environment analysis in ovarian cancer.

Abstract Genome-wide association studies have identified 40 ovarian cancer risk loci. However, the mechanisms underlying these associations remain elusive. In this study, we conducted a two-pronged approach to identify candidate causal SNPs and assess underlying biological mechanisms at chromosome 9p22.2, the first and most statistically significant associated locus for ovarian cancer susceptibility. Three transcriptional regulatory elements with allele-specific effects and a scaffold/matrix attachment region were characterized and, through physical DNA interactions, BNC2 was established as the most likely target gene. We determined the consensus binding sequence for BNC2 in vitro, verified its enrichment in BNC2 ChIP-seq regions, and validated a set of its downstream target genes. Fine-mapping by dense regional genotyping in over 15,000 ovarian cancer cases and 30,000 controls identified SNPs in the scaffold/matrix attachment region as among the most likely causal variants. This study reveals a comprehensive regulatory landscape at 9p22.2 and proposes a likely mechanism of susceptibility to ovarian cancer. Significance: Mapping the 9p22.2 ovarian cancer risk locus identifies BNC2 as an ovarian cancer risk gene. See related commentary by Choi and Brown, p. 439

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 Background The immunohistochemical (IHC) 4 + clinical treatment score (C) i.e. IHC4+C assay combines IHC assessment of four routinely measured markers in the clinical setting, including estrogen receptor (ER), progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and KI67, with other clinical parameters for aiding decision-making on adjuvant chemotherapy in hormone positive (HR+) breast cancer. It is low cost, widely available and has been shown to provide similar prognostic information to the 21-gene (Oncotype DX) test. Clinical use of IHC4+C is limited by the need for standardized scoring of the IHC markers. Here, we examine the prognostic value of a digital image analysis (DIA) based IHC4+C assay in two study populations and compare this with pathologists' visual assay in a subset of patients. We also evaluated its potential clinical relevance in relation to the Nottingham prognostic index [NPI; i.e. low (<3.4), intermediate (3.4-5.4) and high (>5.4)] and IHC defined breast cancer subtypes (i.e. luminal A-like and B-like). Methods Using DIA data on ER, PR, HER2 and KI67, in addition to clinicopathological characteristics of tumors, we computed DIA based IHC4+C (i.e. DIAIHC4+C) assay for 2,034 women from study populations in Poland (PBCS) and the United Kingdom (SEARCH). A total of 246 deaths occurred over a median follow-up period of 6.7 years. A relevant threshold was determined using quartile analysis and the relationship between dichotomous categories (>75th vs ≤ 75th percentile) of DIAIHC4+C and 10-year breast cancer specific survival (BCSS) was investigated using Kaplan-Meier survival curves and in Cox proportional hazard regression models, overall and stratified by study, NPI and subtype. Results DIAIHC4+C score was associated with worse 10-year BCSS overall [HR (95% CI) >75th vs ≤ 75th percentile = 2.92 (1.77, 4.82); p-value <0.0001] and in both the PBCS [3.48 (1.40, 8.65); p-value = 0.007] and SEARCH [3.87 (2.24, 6.68); p-value <0.0001] study populations. DIAIHC4+C was associated with significantly worse survival outcomes in women with intermediate NPI [1.94 (1.31, 2.88); p-value = 0.001] and luminal A-like subtype of HR+ tumors [5.37 (3.50, 8.26); p-value <0.0001]. DIA and visual IHC4+C results were highly correlated (r = 0.95) and survival associations and model fit parameters were similar for DIA [3.82 (2.08, 7.00); p-value < 0.0001; AIC = 472.3] and visual [4.73 (2.55, 8.76); p-value < 0.0001; AIC = 466.1] IHC4+C scores amongst 444 patients with data from both sources. Conclusion These findings are in support of the prognostic value of DIAIHC4+C assay in HR+ breast cancer and demonstrate its potential clinical relevance as an adjunct to the NPI. Given the advantages of being high-throughput and highly reproducible, DIA based methods could find use in the standardization of the analytical processes for the IHC4+C assay. Citation Format: Mustapha Abubakar, Paul D. Pharoah, Mark E. Sherman, Montserrat Garcia-Closas, Mitch Dowsett. Digital image analysis based IHC4+C assay and prognosis in hormone receptor-positive 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 2971.

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 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.

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.

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 &lt; 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 &lt; 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.