Goode, Ellen L.

AbstractX chromosome inactivation (XCI) is a key epigenetic gene expression regulatory process, which may play a role in women’s cancer. In particular tissues, some genes are known to escape XCI, yet patterns of XCI in ovarian cancer (OC) and their clinical associations are largely unknown. To examine XCI in OC, we integrated germline genotype with tumor copy number, gene expression and DNA methylation information from 99 OC patients. Approximately 10% of genes showed different XCI status (either escaping or being subject to XCI) compared with the studies of other tissues. Many of these genes are known oncogenes or tumor suppressors (e.g. DDX3X, TRAPPC2 and TCEANC). We also observed strong association between cis promoter DNA methylation and allele-specific expression imbalance (P = 2.0 × 10−10). Cluster analyses of the integrated data identified two molecular subgroups of OC patients representing those with regulated (N = 47) and dysregulated (N = 52) XCI. This XCI cluster membership was associated with expression of X inactive specific transcript (P = 0.002), a known driver of XCI, as well as age, grade, stage, tumor histology and extent of rl disease following surgical debulking. Patients with dysregulated XCI (N = 52) had shorter time to recurrence (HR = 2.34, P = 0.001) and overall survival time (HR = 1.87, P = 0.02) than those with regulated XCI, although results were attenuated after covariate adjustment. Similar findings were observed when restricted to high-grade serous tumors. We found evidence of a unique OC XCI profile, suggesting that XCI may play an important role in OC biology. Additional studies to examine somatic changes with paired tumor-normal tissue are needed.

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 Introduction: Mitochondria are intracellular organelles that play an important role in energy metabolism, reactive oxygen species generation, and apoptosis. They possess their own haploid genome distinct from nuclear DNA, and mutations in mitochondrial DNA (mtDNA) have previously been linked to cancer risk, disease progression, and treatment response for multiple tumor types. However, the spectrum of mtDNA variation in epithelial ovarian carcinoma (EOC) and associations with clinical characteristics have not been fully explored. Methods: We extracted mtDNA from fresh frozen EOC tumors of women diagnosed with primary invasive epithelial ovarian, fallopian tube, or peritoneal carcinoma, and performed targeted sequencing using the Illumina MiSeq Reporter mtDNA workflow. Homoplasmic and heteroplasmic SNVs and indels were called using an in-house bioinformatics pipeline. We estimated mtDNA haplogroup using Haplogrep and compared to results from available germline mtDNA SNP genotyping. We derived regional weighted variant burden scores based on mtDNA functional domains, accounting for underlying haplogroup, heteroplasmic fraction, and predicted functional impact. To assess potential differences in variant burden by histotype, we compared burden scores using Kruskal-Wallis testing. Among high-grade serous tubo-ovarian cases, we explored associations with overall survival (OS) using Cox proportional hazards regression. After quality control filtering, 389 mitochondrial genomes were analyzed, including those from 320 tubo-ovarian high-grade serous tumors, 42 endometrioid tumors, and 18 clear cell tumors. Results: Accounting for predicted haplogroup, we identified a median count of two global private heteroplasmies per tumor (range: 0-9). Elevated mutation rates were identified in the regulatory control region as well as among the various mt-tRNAs (>3 per 10kb per sample) relative to the remainder of the mtDNA. The majority of recurrent variants across tumors were also identified in the mtDNA control region. We did not observe an association with overall mtDNA variant burden and patient age at diagnosis. Regional mtDNA burden scores revealed a significant association with histotype (P<0.001) for the oxidative phosphorylation complex IV genes, with higher burden in tubo-ovarian high-grade serous and clear-cell tumors relative to endometrioid tumors. For patients with tubo-ovarian high-grade serous tumors, increased variant burden in MT-CO1 (cytochrome C oxidase I) was associated with improved OS (HR = 0.32; 95% CI = [0.14,0.70], LRT P = 0.001) after adjustment for age, stage, and debulking status. Conclusions: We demonstrate a wide range of tumor mtDNA variant patterns in the largest known collection of EOC patients. This work provides the first extensive evaluation of mtDNA variation in these diseases and highlights the potential prognostic relevance of mitochondrial dysfunction. Citation Format: Nicholas B. Larson, Stacey J. Winham, Chen Wang, Jared M. Evans, Sebastian M. Armasu, Bryan M. McCauley, Sean J. Yoder, Lan Min Zhang, Catherine M. Phelan, Ellen L. Goode. Targeted deep sequencing of mitochondrial DNA in epithelial ovarian carcinomas [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 420.

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.

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.