Huntsman, David G.

Abstract Purpose: To identify molecular subclasses of clear cell ovarian carcinoma (CCOC) and assess their impact on clinical presentation and outcomes. Experimental Design: We profiled 421 primary CCOCs that passed quality control using a targeted deep sequencing panel of 163 putative CCOC driver genes and whole transcriptome sequencing of 211 of these tumors. Molecularly defined subgroups were identified and tested for association with clinical characteristics and overall survival. Results: We detected a putative somatic driver mutation in at least one candidate gene in 95% (401/421) of CCOC tumors including ARID1A (in 49% of tumors), PIK3CA (49%), TERT (20%), and TP53 (16%). Clustering of cancer driver mutations and RNA expression converged upon two distinct subclasses of CCOC. The first was dominated by ARID1A-mutated tumors with enriched expression of canonical CCOC genes and markers of platinum resistance; the second was largely comprised of tumors with TP53 mutations and enriched for the expression of genes involved in extracellular matrix organization and mesenchymal differentiation. Compared with the ARID1A-mutated group, women with TP53-mutated tumors were more likely to have advanced-stage disease, no antecedent history of endometriosis, and poorer survival, driven by their advanced stage at presentation. In women with ARID1A-mutated tumors, there was a trend toward a lower rate of response to first-line platinum-based therapy. Conclusions: Our study suggests that CCOC consists of two distinct molecular subclasses with distinct clinical presentation and outcomes, with potential relevance to both traditional and experimental therapy responsiveness. See related commentary by Lheureux, p. 4838

Abstract ARID1A (BAF250a) is a component of the SWI/SNF chromatin modifying complex, plays an important tumour suppressor role, and is considered prognostic in several malignancies. However, in ovarian carcinomas there are contradictory reports on its relationship to outcome, immune response, and correlation with clinicopathological features. We assembled a series of 1623 endometriosis‐associated ovarian carcinomas, including 1078 endometrioid (ENOC) and 545 clear cell (CCOC) ovarian carcinomas, through combining resources of the Ovarian Tumor Tissue Analysis (OTTA) Consortium, the Canadian Ovarian Unified Experimental Resource (COEUR), local, and collaborative networks. Validated immunohistochemical surrogate assays for ARID1A mutations were applied to all samples. We investigated associations between ARID1A loss/mutation, clinical features, outcome, CD8 + tumour‐infiltrating lymphocytes (CD8 + TILs), and DNA mismatch repair deficiency (MMRd). ARID1A loss was observed in 42% of CCOCs and 25% of ENOCs. We found no associations between ARID1A loss and outcomes, stage, age, or CD8 + TIL status in CCOC. Similarly, we found no association with outcome or stage in endometrioid cases. In ENOC, ARID1A loss was more prevalent in younger patients ( p = 0.012) and was associated with MMRd ( p < 0.001) and the presence of CD8 + TILs ( p = 0.008). Consistent with MMRd being causative of ARID1A mutations, in a subset of ENOCs we also observed an association with ARID1A loss‐of‐function mutation as a result of small indels ( p = 0.035, versus single nucleotide variants). In ENOC, the association with ARID1A loss, CD8 + TILs, and age appears confounded by MMRd status. Although this observation does not explicitly rule out a role for ARID1A influence on CD8 + TIL infiltration in ENOC, given current knowledge regarding MMRd, it seems more likely that effects are dominated by the hypermutation phenotype. This large dataset with consistently applied biomarker assessment now provides a benchmark for the prevalence of ARID1A loss‐of‐function mutations in endometriosis‐associated ovarian cancers and brings clarity to the prognostic significance. © 2021 The Pathological Society of Great Britain and Ireland.

Abstract Background: Accumulating evidence suggests a relationship between endometrial cancer and ovarian cancer. Independent genome-wide association studies (GWAS) for endometrial cancer and ovarian cancer have identified 16 and 27 risk regions, respectively, four of which overlap between the two cancers. We aimed to identify joint endometrial and ovarian cancer risk loci by performing a meta-analysis of GWAS summary statistics from these two cancers. Methods: Using LDScore regression, we explored the genetic correlation between endometrial cancer and ovarian cancer. To identify loci associated with the risk of both cancers, we implemented a pipeline of statistical genetic analyses (i.e., inverse-variance meta-analysis, colocalization, and M-values) and performed analyses stratified by subtype. Candidate target genes were then prioritized using functional genomic data. Results: Genetic correlation analysis revealed significant genetic correlation between the two cancers (rG = 0.43, P = 2.66 × 10−5). We found seven loci associated with risk for both cancers (PBonferroni &lt; 2.4 × 10−9). In addition, four novel subgenome-wide regions at 7p22.2, 7q22.1, 9p12, and 11q13.3 were identified (P &lt; 5 × 10−7). Promoter-associated HiChIP chromatin loops from immortalized endometrium and ovarian cell lines and expression quantitative trait loci data highlighted candidate target genes for further investigation. Conclusions: Using cross-cancer GWAS meta-analysis, we have identified several joint endometrial and ovarian cancer risk loci and candidate target genes for future functional analysis. Impact: Our research highlights the shared genetic relationship between endometrial cancer and ovarian cancer. Further studies in larger sample sets are required to confirm our findings.

Abstract Background: High-grade serous ovarian cancer (HGSOC) is the most common and aggressive epithelial ovarian cancer. HGSOC has poor survival, with five-year survival of less than 40%, as it is often diagnosed at a late stage. The median overall survival (OS) is approximately four years. Prognostic gene expression signatures for HGSOC have been inconsistent and have not been translated into clinical practice. The aim of this study was to develop a robust prognostic signature for HGSOC, in order to identify high-risk patients requiring alternative treatments. Methods: A meta-analysis of 1,455 tumors was performed and 200 prognostic genes were selected for expression profiling. An additional 315 candidate genes were selected from the literature, drug-targetable pathways, and other hypotheses. NanoString expression analysis of the 513 genes was performed on formalin-fixed, paraffin-embedded (FFPE) tumor tissue from 3,769 women with HGSOC from the Ovarian Tumor Tissue Analysis (OTTA) Consortium. A prognostic model for 5-year survival was developed using regression-based and machine learning methods. The model was trained on 2,702 tumors and evaluated on an independent set of 1,067 tumors. A NanoString expression analysis was performed on 2,692 additional tumors, for a subset of 218 of the genes, including the prognostic signature. Results: Expression of 276 genes was associated with OS (false discovery rate (FDR) &lt; 0.05) in single-gene analyses adjusted for age, stage, and race. The top five genes were TAP1, ZFHX4, CXCL9, FBN1, and PTGER3 (p &lt;&lt; 0.001). A 101-gene prognostic signature was established, where a gain of one standard deviation in the gene expression score conferred a greater than two-fold increase in risk of death (HR (hazard ratio) = 2.35 [2.02, 2.71]; p &lt;&lt; 0.001). Median survival by quintile group was 9.5, 5.4, 3.8, 3.2, and 2.3 years. Sensitivity analyses, applying the signature to specific patient groups, showed that the prognostic power of the signature is not explained by residual disease, treatment, BRCA status, or CD8 score. The prognostic signature was enriched in pathways with treatment implications, such as PI3K, GPCR, and immune. This signature is suitable for use in FFPE tumor material and future cases can be assigned an HR from the gene expression score. The signature was developed in cases that did not receive neoadjuvant treatment (NACT); however, application to NACT cases was assessed in the second NanoString expression analysis. Conclusion: The OTTA-SPOT (Ovarian Tumor Tissue Analysis consortium – Stratified Prognosis of Ovarian Tumors) gene expression signature may improve risk stratification in clinical trials by identifying patients who are least likely to achieve 5-year survival. Citation Format: Joshua Millstein, Timothy Budden, Ellen L. Goode, Michael S. Anglesio, Aline Talhouk, OTTA consortium, David G. Huntsman, David D. Bowtell, James D. Brenton, Jennifer A. Doherty, Paul P.D. Pharoah, Susan J. Ramus. The Ovarian Tumour Tissue Analysis Consortium: Stratified Prognosis of Ovarian Tumors (OTTA-SPOT) signature for high-grade serous ovarian cancer [abstract]. In: Proceedings of the AACR Special Conference on Advances in Ovarian Cancer Research; 2019 Sep 13-16, 2019; Atlanta, GA. Philadelphia (PA): AACR; Clin Cancer Res 2020;26(13_Suppl):Abstract nr IA12.

Abstract Purpose: Gene expression–based molecular subtypes of high-grade serous tubo-ovarian cancer (HGSOC), demonstrated across multiple studies, may provide improved stratification for molecularly targeted trials. However, evaluation of clinical utility has been hindered by nonstandardized methods, which are not applicable in a clinical setting. We sought to generate a clinical grade minimal gene set assay for classification of individual tumor specimens into HGSOC subtypes and confirm previously reported subtype-associated features. Experimental Design: Adopting two independent approaches, we derived and internally validated algorithms for subtype prediction using published gene expression data from 1,650 tumors. We applied resulting models to NanoString data on 3,829 HGSOCs from the Ovarian Tumor Tissue Analysis consortium. We further developed, confirmed, and validated a reduced, minimal gene set predictor, with methods suitable for a single-patient setting. Results: Gene expression data were used to derive the predictor of high-grade serous ovarian carcinoma molecular subtype (PrOTYPE) assay. We established a de facto standard as a consensus of two parallel approaches. PrOTYPE subtypes are significantly associated with age, stage, residual disease, tumor-infiltrating lymphocytes, and outcome. The locked-down clinical grade PrOTYPE test includes a model with 55 genes that predicted gene expression subtype with &gt;95% accuracy that was maintained in all analytic and biological validations. Conclusions: We validated the PrOTYPE assay following the Institute of Medicine guidelines for the development of omics-based tests. This fully defined and locked-down clinical grade assay will enable trial design with molecular subtype stratification and allow for objective assessment of the predictive value of HGSOC molecular subtypes in precision medicine applications. See related commentary by McMullen et al., p. 5271

3639 Background: The first phase of the multi-center prospectively specified retrospective study Validating Indicators To Associate Recurrence (VITAR), assessing the relationship between GCC gene expression in formalin fixed (FFPE) LNs and time to recurrence (TTR) in stage II CC pts not treated with adjuvant chemotherapy (Sargent, Annals Surg Onc 2011), showed promising initial results. Here we report a validation set of 463 new stage II CC pts. Methods: GCC mRNA was quantified by RT-qPCR using FFPE LNs from untreated T3N0 CC pts diagnosed from 1999-2008 with at least 12 LNs examined , blinded to clinical outcomes. Patients were classified by GCC LN ratio (LNR) (high risk: LNR > 0.1; low risk: LNR ≤ 0.1), with LNR defined as ratio of GCC positive to GCC informative LNs. Cox regression models tested the relationship between GCC and the primary endpoint of TTR, adjusted for age, tumor grade, number of LN examined pathologically, and lymphovascular invasion. Mismatch repair (MMR) status was also assessed. All primary analyses and cut-points were pre-specified. Results: 46pts (10%) recurred (rec), median follow-up was 65 months, median LNs examined was 20, and 42% (195/463) were classified high risk. Overall, TTR was not significantly associated with binary GCC LNR risk class (HR=1.47, p=.208) or DFS (HR= 1.39, p=.097). One site’s (n=97) tissue grossing method precluded appropriate LN assessment with existing GCC qualification methods. Excluding this site resulted in a TTR HR=1.91, p=0.051 (multivariate). In a post-hocanalysis excluding this site and using a 3-level GCC risk group of high (LNR > 0.20), intermediate (0.10 < LNR < 0.20) and low (LNR < 0.10), high risk group pts had a 5-yr rec risk of 22% versus 8% in low risk (HR 2.72, p=0.006). MMR status was not significantly associated with TTR (multivariate p=0.30). Conclusions: GCC status is a promising prognostic factor in appropriately staged stage II CC pts not treated with adjuvant therapy independent of traditional histopathology risk factors, but GCC determination must be performed with methodology adapted to the tissue procurement and fixation technique. Outcome associations were strengthened when considering a 3-level GCC categorization.

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

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.