Citations

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.

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.

Abstract Background: Endometrioid carcinoma (EC) and clear cell carcinoma (CC) histotypes of epithelial ovarian cancer are understudied compared with the more common high-grade serous carcinomas (HGSC). We therefore sought to characterize EC and CC transcriptomes in relation to HGSC. Methods: Following bioinformatics processing and gene abundance normalization, differential expression analysis of RNA sequence data collected on fresh-frozen tumors was completed with nonparametric statistical analysis methods (55 ECs, 19 CCs, 112 HGSCs). Association of gene expression with progression-free survival (PFS) was completed with Cox proportional hazards models. Eight additional multi-histotype expression array datasets (N = 852 patients) were used for replication. Results: In the discovery set, tumors generally clustered together by histotype. Thirty-two protein-coding genes were differentially expressed across histotype (P &lt; 1 × 10−10) and showed similar associations in replication datasets, including MAP2K6, KIAA1324, CDH1, ENTPD5, LAMB1, and DRAM1. Nine genes associated with PFS (P &lt; 0.0001) showed similar associations in replication datasets. In particular, we observed shorter PFS time for CC and EC patients with high gene expression for CCNB2, CORO2A, CSNK1G1, FRMD8, LIN54, LINC00664, PDK1, and PEX6, whereas, the converse was observed for HGSC patients. Conclusions: The results suggest important histotype differences that may aid in the development of treatment options, particularly those for patients with EC or CC. Impact: We present replicated findings on transcriptomic differences and how they relate to clinical outcome for two of the rarer ovarian cancer histotypes of EC and CC, along with comparison with the common histotype of HGSC. Cancer Epidemiol Biomarkers Prev; 27(9); 1101–9. ©2018 AACR.