Gayther, Simon A.

Abstract BRCA-associated homologous recombination deficiency (HRD) is present in ~ 50% of high-grade serous carcinomas (HGSC) and predicts sensitivity to platinum-based therapy. However, there is little understanding of why some patients with BRCA-deficient tumors experience unexpectedly poor outcomes. We profiled 154 tumors, enriched for patients with BRCA-deficient tumors that experienced short overall survival (≤ 3 years, n = 42), using whole-genome, transcriptome, and methylation analyses. All but one BRCA-deficient tumor exceeded an accepted HRD genomic scarring threshold. However, patients with BRCA1-deficient HGSC with a more elevated HRD score survived significantly longer. Patients with BRCA2-deficient HGSC and loss of NF1 survived twice as long as those without NF1 loss, whereas PIK3CA or RAD21 amplification defined BRCA2-deficient HGSC with exceptionally short survival. BRCA1-deficient tumors in short survivors had evidence of immunosuppressive c-kit signaling and EMT. In a large HGSC cohort (n = 1,389) including 282 individuals with pathogenic germline BRCA variants (gBRCApv), the location of the mutation within functional domains stratified clinical outcomes. Notably, residual disease after primary surgery had limited prognostic effect in gBRCApv-carriers compared to non-carriers. Our findings indicate that tumor HR proficiency in the context of therapy response and survival is not a binary property, and highlight genomic and immune modifiers of outcomes in BRCA-deficient HGSC.

Abstract Background: Ovarian high-grade serous carcinomas (HGSC) comprise four distinct molecular subtypes based on mRNA expression patterns, with differential survival. Understanding risk factor associations is important to elucidate the etiology of HGSC. We investigated associations between different epidemiologic risk factors and HGSC molecular subtypes. Methods: We pooled data from 11 case–control studies with epidemiologic and tumor gene expression data from custom NanoString CodeSets developed through a collaboration within the Ovarian Tumor Tissue Analysis consortium. The PrOTYPE-validated NanoString-based 55-gene classifier was used to assign HGSC gene expression subtypes. We examined associations between epidemiologic factors and HGSC subtypes in 2,070 cases and 16,633 controls using multivariable-adjusted polytomous regression models. Results: Among the 2,070 HGSC cases, 556 (27%) were classified as C1.MES, 340 (16%) as C5.PRO, 538 (26%) as C2.IMM, and 636 (31%) as C4.DIF. The key factors, including oral contraceptive use, parity, breastfeeding, and family history of ovarian cancer, were similarly associated with all subtypes. Heterogeneity was observed for several factors. Former smoking [OR = 1.25; 95% confidence interval (CI) = 1.03, 1.51] and genital powder use (OR = 1.42; 95% CI = 1.08, 1.86) were uniquely associated with C2.IMM. History of endometriosis was associated with C5.PRO (OR = 1.46; 95% CI = 0.98, 2.16) and C4.DIF (OR = 1.27; 95% CI = 0.94, 1.71) only. Family history of breast cancer (OR = 1.44; 95% CI = 1.16, 1.78) and current smoking (OR = 1.40; 95% CI = 1.11, 1.76) were associated with C4.DIF only. Conclusions: This study observed heterogeneous associations of epidemiologic and modifiable factors with HGSC molecular subtypes. Impact: The different patterns of associations may provide key information about the etiology of the four subtypes.

Abstract Purpose: The purpose of this study was to evaluate RB1 expression and survival across ovarian carcinoma histotypes and how co-occurrence of BRCA1 or BRCA2 (BRCA) alterations and RB1 loss influences survival in tubo-ovarian high-grade serous carcinoma (HGSC). Experimental Design: RB1 protein expression was classified by immunohistochemistry in ovarian carcinomas of 7,436 patients from the Ovarian Tumor Tissue Analysis consortium. We examined RB1 expression and germline BRCA status in a subset of 1,134 HGSC, and related genotype to overall survival (OS), tumor-infiltrating CD8+ lymphocytes, and transcriptomic subtypes. Using CRISPR-Cas9, we deleted RB1 in HGSC cells with and without BRCA1 alterations to model co-loss with treatment response. We performed whole-genome and transcriptome data analyses on 126 patients with primary HGSC to characterize tumors with concurrent BRCA deficiency and RB1 loss. Results: RB1 loss was associated with longer OS in HGSC but with poorer prognosis in endometrioid ovarian carcinoma. Patients with HGSC harboring both RB1 loss and pathogenic germline BRCA variants had superior OS compared with patients with either alteration alone, and their median OS was three times longer than those without pathogenic BRCA variants and retained RB1 expression (9.3 vs. 3.1 years). Enhanced sensitivity to cisplatin and paclitaxel was seen in BRCA1-altered cells with RB1 knockout. Combined RB1 loss and BRCA deficiency correlated with transcriptional markers of enhanced IFN response, cell-cycle deregulation, and reduced epithelial–mesenchymal transition. CD8+ lymphocytes were most prevalent in BRCA-deficient HGSC with co-loss of RB1. Conclusions: Co-occurrence of RB1 loss and BRCA deficiency was associated with exceptionally long survival in patients with HGSC, potentially due to better treatment response and immune stimulation.

Abstract Introduction: Identifying dysregulated gene expression in ovarian cancers has been limited by a deficit of available normal tissues. Here, we generated the largest set of high-grade serous ovarian cancer (HGSOC) tumors with normal precursor tissues for transcriptome analyses to identify lncRNA expression patterns, provide insight into their cellular contexts and functions, and assess their relevance to HGSOC pathogenesis and genetic susceptibility. Methods: We performed RNA sequencing on 220 primary HGSOCs and 116 benign epithelia (micro-dissected fallopian tube, ovarian surface, and inclusion cyst epithelia), and combined samples with 428 HGSOCs from TCGA, 60 HGSOCs from a prior study, and 180 bulk ovary tissues from GTEx. Reads were processed with a uniform bioinformatic and quality control pipeline, combined and batch corrected. Normalized TPM expression values were compared using the Wilcoxon rank sum test with >2-fold change and FDR <1% considered statistically significant. We computed weighted gene co-expression networks for each tissue type and tested for enrichment of differentially expressed genes (DEGs) and gene ontologies within co-expression modules. Results: Transcriptomes included 27,700 expressed genes (8,202 lncRNAs) across 706 HGSOCs, 180 bulk ovary, and 88 ovarian epithelia tissues. Most genes were expressed in all tissues, though 5% each showed HGSOC- or normal tissue-specific expression which comprised ≥50% lncRNA. Comparing HGSOCs to normal tissues identified 11,804 DEGs with 4,522 lncRNAs (DElncRNA) of which ~50% were tissue-specific. DEGs included MUC16 and multiple GWAS/TWAS implicated susceptibility genes including RAD51, BRIP1, BNC2, TIPARP-AS1, PRC1, KANSL1, ANKLE1, CHMP4C, ESRP2, and CCNE1. The most highly expressed DElncRNA in HGSOC were upregulated RMRP (P=1.4 × 10−39), SNHG1 (P=3.0 × 10−27), and HAGLR (P=2.0 × 10−24) at the HOXD risk locus. Lower proportions of lncRNA were seen in the HGSOC co-expression modules compared to those in normal tissues and less (3% vs. 15%) HGSOC modules were considered strongly lncRNA regulated (>50% lncRNA in module). Tissue-specific DElncRNA tended toward upregulation in HGSOC compared to ovarian epithelia with enriched modules associated with cell cycle regulation (P=1.5 × 10−38); and toward downregulation in HGSOC compared to bulk ovary with enriched modules associated with chemokine signaling/response (P=1.7 × 10−113). Conclusions: The lower proportion of lncRNA involved in HGSOC co-expression networks indicates a marked loss of regulation by lncRNA. Preliminary interrogation of HGSOC-dysregulated lncRNA expression revealed tissue-specific differences that highlight unique biological pathways in precursor epithelia and the ovarian microenvironment that contribute to HGSOC pathogenesis. Our results provide additional support for previously nominated risk genes. eQTL analysis has recently been completed and integration with GWAS are underway to elucidate lncRNA and lncRNA-regulated susceptibility genes and will be presented at the conference. Citation Format: Brett M. Reid, Ann Chen, Zhihua Chen, Florian A. Karreth, Peter Kanetsky, Jennifer B. Permuth, Ozlen Saglam, Jamie K. Teer, Xiaoqing Yu, Simon Gayther, Ellen Goode, Paul Pharoah, Thomas A. Sellers, Kate Lawrenson, Jonathan Tyrer. Patterns of lncRNA-regulated gene expression in high-grade serous ovarian carcinomas [abstract]. In: Proceedings of the AACR Special Conference on Ovarian Cancer; 2023 Oct 5-7; Boston, Massachusetts. Philadelphia (PA): AACR; Cancer Res 2024;84(5 Suppl_2):Abstract nr A004.

Abstract Purpose: Identifying dysregulated gene expression in ovarian cancers has been limited by a deficit of available normal tissues. Here, we generated the largest set of high-grade serous ovarian cancer (HGSOC) tumors with normal precursor tissues for transcriptome analyses. Methods: We performed RNA sequencing on 220 primary HGSOCs and 116 benign epithelia (micro-dissected fallopian tube, ovarian surface, and inclusion cysts), and combined samples with 428 HGSOCs from TCGA, 60 HGSOCs from a prior study, and 180 bulk ovary tissues from GTEx. Raw reads were processed with a uniform bioinformatic and quality control pipeline; combined data were batch corrected. We tested for differences in median normalized CPM expression values using the Wilcoxon rank sum test with >2-fold change and a false discovery rate <1% considered statistically significant. We also conducted weighted gene co-expression network analysis in each tissue type. The hypergeometric test was used for enrichment of differentially expressed genes (DEGs) and gene ontologies within modules. Results: Transcriptomes comprised 27,700 expressed genes (8,202 lncRNAs) across 706 HGSOCs, 180 bulk ovary, and 88 ovarian epithelia tissues. Most (~90%) genes were expressed in all tissues; 5% each showed HGSOC- or normal tissue-specific expression and ≥50% were lncRNA. Comparing HGSOCs to ovarian epithelia and to bulk ovary identified 11,804 DEGs with 4,522 lncRNAs (DElncRNA) of which ~50% were tissue-specific. DEGs included MUC16 and multiple GWAS/TWAS implicated susceptibility genes including RAD51, BRIP1, BNC2, TIPARP-AS1, PRC1, KANSL1, ANKLE1, CHMP4C, ESRP2, and CCNE1. The most highly expressed DElncRNA in HGSOC were upregulated RMRP (P=1.4x10-39), SNHG1 (P=3.0x10-27), and HAGLR (P=2.0x10-24) at the HOXD risk locus. The highest expressed DElncRNA in precursor tissues was the HGSOC-downregulated MEG3 (P=1.7x10-67). DEGs were enriched in HGSOC co-expression modules associated with immune response, cell motility/localization, cell cycle regulation, angiogenesis, reproductive development, transcriptional regulation, and metabolic processes. Tissue-specific DElncRNA tended toward upregulation compared to ovarian epithelia with enriched modules associated with cell cycle regulation (hub=BUB1B; top DElncRNA TRPM2-AS, P=1.5x10-38); and toward downregulation compared to bulk ovary with enriched modules associated with chemokine signaling/response (hub= GADD45B, top DElncRNA RP11-87P3.1, P=1.7x10-113). Conclusion: HGSOC-dysregulated lncRNA expression revealed tissue-specific differences that highlight unique biological pathways in precursor epithelia and the ovarian microenvironment that contribute to HGSOC pathogenesis. Our results provide additional evidence to support previously nominated risk genes. Integration with eQTL and GWAS are underway to further elucidate novel HGSOC susceptibility genes. Citation Format: Brett M. Reid, Ann Chen, Zhihua Chen, Florian A. Karreth, Peter Kanetsky, Jennifer B. Permuth, Ozlen Saglam, Jamie Teer, Xiaoqing Yu, Simon Gayther, Ellen Goode, Paul Pharoah, Thomas A. Sellers, Kate Lawrenson. Patterns of dysregulated coding and noncoding gene expression in high-grade serous ovarian carcinomas [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 6065.

AbstractOur objective was to test whether p53 expression status is associated with survival for women diagnosed with the most common ovarian carcinoma histotypes (high‐grade serous carcinoma [HGSC], endometrioid carcinoma [EC], and clear cell carcinoma [CCC]) using a large multi‐institutional cohort from the Ovarian Tumor Tissue Analysis (OTTA) consortium. p53 expression was assessed on 6,678 cases represented on tissue microarrays from 25 participating OTTA study sites using a previously validated immunohistochemical (IHC) assay as a surrogate for the presence and functional effect of TP53 mutations. Three abnormal expression patterns (overexpression, complete absence, and cytoplasmic) and the normal (wild type) pattern were recorded. Survival analyses were performed by histotype. The frequency of abnormal p53 expression was 93.4% (4,630/4,957) in HGSC compared to 11.9% (116/973) in EC and 11.5% (86/748) in CCC. In HGSC, there were no differences in overall survival across the abnormal p53 expression patterns. However, in EC and CCC, abnormal p53 expression was associated with an increased risk of death for women diagnosed with EC in multivariate analysis compared to normal p53 as the reference (hazard ratio [HR] = 2.18, 95% confidence interval [CI] 1.36–3.47, p = 0.0011) and with CCC (HR = 1.57, 95% CI 1.11–2.22, p = 0.012). Abnormal p53 was also associated with shorter overall survival in The International Federation of Gynecology and Obstetrics stage I/II EC and CCC. Our study provides further evidence that functional groups of TP53 mutations assessed by abnormal surrogate p53 IHC patterns are not associated with survival in HGSC. In contrast, we validate that abnormal p53 IHC is a strong independent prognostic marker for EC and demonstrate for the first time an independent prognostic association of abnormal p53 IHC with overall survival in patients with CCC.

Abstract Background The role of ovulation in epithelial ovarian cancer (EOC) is supported by the consistent protective effects of parity and oral contraceptive use. Whether these factors protect through anovulation alone remains unclear. We explored the association between lifetime ovulatory years (LOY) and EOC. Methods LOY was calculated using 12 algorithms. Odds ratios (ORs) and 95% confidence intervals (CIs) estimated the association between LOY or LOY components and EOC among 26 204 control participants and 21 267 case patients from 25 studies. To assess whether LOY components act through ovulation suppression alone, we compared beta coefficients obtained from regression models with expected estimates assuming 1 year of ovulation suppression has the same effect regardless of source. Results LOY was associated with increased EOC risk (OR per year increase = 1.014, 95% CI = 1.009 to 1.020 to OR per year increase = 1.044, 95% CI = 1.041 to 1.048). Individual LOY components, except age at menarche, also associated with EOC. The estimated model coefficient for oral contraceptive use and pregnancies were 4.45 times and 12- to 15-fold greater than expected, respectively. LOY was associated with high-grade serous, low-grade serous, endometrioid, and clear cell histotypes (ORs per year increase = 1.054, 1.040, 1.065, and 1.098, respectively) but not mucinous tumors. Estimated coefficients of LOY components were close to expected estimates for high-grade serous but larger than expected for low-grade serous, endometrioid, and clear cell histotypes. Conclusions LOY is positively associated with nonmucinous EOC. Differences between estimated and expected model coefficients for LOY components suggest factors beyond ovulation underlie the associations between LOY components and EOC in general and for non-HGSOC.

AbstractBackgroundCyclin E1 (CCNE1) is a potential predictive marker and therapeutic target in tubo‐ovarian high‐grade serous carcinoma (HGSC). Smaller studies have revealed unfavorable associations for CCNE1 amplification and CCNE1 overexpression with survival, but to date no large‐scale, histotype‐specific validation has been performed. The hypothesis was that high‐level amplification of CCNE1 and CCNE1 overexpression, as well as a combination of the two, are linked to shorter overall survival in HGSC.MethodsWithin the Ovarian Tumor Tissue Analysis consortium, amplification status and protein level in 3029 HGSC cases and mRNA expression in 2419 samples were investigated.ResultsHigh‐level amplification (>8 copies by chromogenic in situ hybridization) was found in 8.6% of HGSC and overexpression (>60% with at least 5% demonstrating strong intensity by immunohistochemistry) was found in 22.4%. CCNE1 high‐level amplification and overexpression both were linked to shorter overall survival in multivariate survival analysis adjusted for age and stage, with hazard stratification by study (hazard ratio [HR], 1.26; 95% CI, 1.08‐1.47, p = .034, and HR, 1.18; 95% CI, 1.05‐1.32, p = .015, respectively). This was also true for cases with combined high‐level amplification/overexpression (HR, 1.26; 95% CI, 1.09‐1.47, p = .033). CCNE1 mRNA expression was not associated with overall survival (HR, 1.00 per 1‐SD increase; 95% CI, 0.94‐1.06; p = .58). CCNE1 high‐level amplification is mutually exclusive with the presence of germline BRCA1/2 pathogenic variants and shows an inverse association to RB1 loss.ConclusionThis study provides large‐scale validation that CCNE1 high‐level amplification is associated with shorter survival, supporting its utility as a prognostic biomarker in HGSC.

AbstractBackgroundKnown risk alleles for epithelial ovarian cancer (EOC) account for approximately 40% of the heritability for EOC. Copy number variants (CNVs) have not been investigated as EOC risk alleles in a large population cohort.MethodsSingle nucleotide polymorphism array data from 13 071 EOC cases and 17 306 controls of White European ancestry were used to identify CNVs associated with EOC risk using a rare admixture maximum likelihood test for gene burden and a by-probe ratio test. We performed enrichment analysis of CNVs at known EOC risk loci and functional biofeatures in ovarian cancer–related cell types.ResultsWe identified statistically significant risk associations with CNVs at known EOC risk genes; BRCA1 (PEOC = 1.60E-21; OREOC = 8.24), RAD51C (Phigh-grade serous ovarian cancer [HGSOC] = 5.5E-4; odds ratio [OR]HGSOC = 5.74 del), and BRCA2 (PHGSOC = 7.0E-4; ORHGSOC = 3.31 deletion). Four suggestive associations (P < .001) were identified for rare CNVs. Risk-associated CNVs were enriched (P < .05) at known EOC risk loci identified by genome-wide association study. Noncoding CNVs were enriched in active promoters and insulators in EOC-related cell types.ConclusionsCNVs in BRCA1 have been previously reported in smaller studies, but their observed frequency in this large population-based cohort, along with the CNVs observed at BRCA2 and RAD51C gene loci in EOC cases, suggests that these CNVs are potentially pathogenic and may contribute to the spectrum of disease-causing mutations in these genes. CNVs are likely to occur in a wider set of susceptibility regions, with potential implications for clinical genetic testing and disease prevention.

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