Pharoah, Paul D.P.

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 Despite homologous recombination (HR) deficiency being an important predictor of survival in high-grade serous ovarian carcinoma (HGSC), its status is rarely considered in large-scale studies of prognostic signatures and biomarkers. We sought to develop a gene expression-based classifier of HR deficient (HRD) HGSCs, suitable for archival FFPE tumor samples, to determine whether established prognostic markers were dependent on HR status. Using 361 HGSCs from The Cancer Genome Atlas (TCGA) where HR status was defined by the genomic gold standard, a HRD score, we trained and tested a ridge penalized logistic regression model. This 218-gene mRNA expression signature was then applied in 6002 HGSCs from the Ovarian Tumor Tissue Analysis (OTTA) consortium to predict the probability of a sample being HRD. In this large cohort of FFPE samples, we sought to understand the relationship between HR and overall survival. Applying Kaplan-Meier analyses and multivariate Cox proportional hazard regression modelling, adjusting for age and stratified by stage and study, we examined the relationship between HR status and the PrOTYPE subtype, SPOT prognostic score and other known prognostic markers in HGSC. Our machine learning model was 84% accurate in a held-out subset of the TCGA (n=73). Of the 6002 OTTA HGSCs, 40.4% were predicted to be HRD. These cases were significantly younger, were enriched for early stage (FIGO I/II) and had improved outcomes when compared to HR proficient (HRP) cases (hazard ratio [HZR]=0.71, 95% confidence interval [CI]=0.66-0.76, p<2x10-16). A linear increase in survival time was observed as the probability of a sample being HRD increased. HGSCs with >75% likelihood of being HRD had a median 5-year survival of 49.3% (95%CI=46.2-52.6%) which was a 17.7% increase from those with <25% likelihood. Both PrOTYPE and SPOT were significantly associated with survival, independent of HR status (p<2x10-16 and p<2x10-16 respectively). Notably the probability of being HRD was negatively correlated with an unfavorable SPOT score (R=-0.32, p<2x10-16), suggesting that SPOT partially captures the survival advantage of HR deficiency. Several individual biomarkers, assessed in 1248 to 2858 cases, showed HR dependent survival associations. HRD HGSCs with high GMNN expression had significantly improved, long-term outcomes (4-10years, HZR=0.39, 95%CI=0.20-0.78, p=0.01), and may be a marker of exceptional survival, while AR negative HRD cases had poorer outcomes, earlier in disease progression (0-4years, HZR=1.43, 95%CI=1.03-1.98, p=0.03). Significantly poorer prognoses were observed in HRP cases expressing PTEN (HZR=1.39, 95%CI=1.04-1.87, p=0.03), and in HRD cases with high expression of MyD88 (HZR=1.21, 95%CI=1.01-1.45, p=0.04). RB1 and FOXJ1 protein expression, and CD8+ T cell abundance were prognostic irrespective of HR status. Our study contributes to understanding the association between HR and HGSC survival in a statistically powered cohort and demonstrates heterogeneity in outcomes through clinically relevant prognostic signatures and biomarkers. Citation Format: Ashley L. Weir, Derek Chiu, Aline Talhouk, Mike S. Anglesio, Dale W. Garsed, Anna DeFazio, Holly H. Harris, Joellen M. Schildkraut, Andrew Berchuck, Ellen L. Goode, Paul D.P. Pharoah, David Huntsman, David Bowtell, Jen A. Doherty, Nadia M. Davidson, Martin Köbel, Susan J. Ramus. Homologous recombination DNA repair dependent survival associations in high-grade serous ovarian carcinoma: an Ovarian Tumor Tissue Analysis consortium study [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Advances in Ovarian Cancer Research; 2025 Sep 19-21; Denver, CO. Philadelphia (PA): AACR; Cancer Res 2025;85(18_Suppl):Abstract nr B046.

BACKGROUND Despite an overall poor prognosis, about 15% of patients with advanced-stage tubo-ovarian high-grade serous carcinoma (HGSC) survive 10 or more years after standard treatment. METHODS We evaluated the tumor microenvironment of this exceptional, understudied group using a large international cohort enriched for long-term survivors (LTS; 10+ years; n = 374) compared with mid-term (MTS; 5–7.99 years; n = 433) and short-term survivors (STS; 2–4.99 years; n = 416). Primary tumor samples were immunostained and scored for intraepithelial and intrastromal densities of 10 immune-cell subsets (including T cells, B cells, plasma cells, myeloid cells, PD-1 + cells, and PD-L1 + cells) and epithelial content. RESULTS Positive associations with LTS compared with STS were seen for 9 of 10 immune-cell subsets. In particular, the combination of intraepithelial CD8 + T cells and intrastromal B cells showed near 5-fold increased odds of LTS compared with STS. All of these associations were stronger in tumors with high epithelial content and/or the C4/Differentiated molecular subtype, despite immune-cell densities generally being higher in tumors with low epithelial content and/or the C2/Immunoreactive molecular subtype. CONCLUSION The tumor microenvironment of HGSC LTS is distinguished by the intersection of T and B cell coinfiltration, high epithelial content, and C4/differentiated molecular subtype, features which may inspire new approaches to immunotherapy. FUNDING Ovarian Cancer Research Program (OCRP) of the Congressionally Directed Medical Research Program (CDMRP), U.S. Department of Defense (DOD); American Cancer Society; BC Cancer Foundation; Canada’s Networks of Centres of Excellence; Canadian Cancer Society; Canadian Institutes of Health Research; Cancer Councils of New South Wales, Victoria, Queensland, South Australia, and Tasmania, Cancer Foundation of Western Australia; Cancer Institute NSW; Cancer Research UK; Deutsche Forschungsgesellschaft; ELAN Funds of the University of Erlangen-Nuremberg; Fred C. and Katherine B. Andersen Foundation; Genome BC; German Cancer Research Center; German Federal Ministry of Education and Research, Programme of Clinical Biomedical Research; Instituto de Salud Carlos III; Mayo Foundation; Minnesota Ovarian Cancer Alliance; Ministerio de Economía y Competitividad; Medical Research Council (MRC); National Center for Advancing Translational Sciences; National Health and Medical Research Council of Australia (NHMRC); Ovarian Cancer Australia; Peter MacCallum Foundation; Sydney West Translational Cancer Research Centre; Terry Fox Research Institute; The Eve Appeal (The Oak Foundation); UK National Institute for Health Research Biomedical Research Centres at the University of Cambridge; University of Pittsburgh School of Medicine; U.S. National Cancer Institute of the National Institutes of Health; VGH & UBC Hospital Foundation; Victorian Cancer Agency.

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 Background: Women with an inherited pathogenic variant in BRCA1 or BRCA2 have a greatly increased risk of developing ovarian cancer, but the importance of behavioral factors is less clear. We used a case-only design to compare the magnitude of associations with established reproductive, hormonal, and lifestyle risk factors between BRCA mutation carriers and noncarriers. Methods: We pooled data from five studies from the Ovarian Cancer Association Consortium including 637 BRCA carriers and 4,289 noncarriers. Covariate-adjusted generalized linear mixed models were used to estimate interaction risk ratios (IRR) and 95% confidence intervals (CI), with BRCA (carrier vs. noncarrier) as the response variable. Results: IRRs were above 1.0 for known protective factors including ever being pregnant (IRR = 1.29, 95% CI; 1.00–1.67) and ever using the oral contraceptive pill (1.30, 95% CI; 1.07–1.60), suggesting the protective effects of these factors may be reduced in carriers compared with noncarriers. Conversely, the IRRs for risk factors including endometriosis and menopausal hormone therapy were below 1.0, suggesting weaker positive associations among BRCA carriers. In contrast, associations with lifestyle factors including smoking, physical inactivity, body mass index, and aspirin use did not appear to differ by BRCA status. Conclusions: Our results suggest that associations with hormonal and reproductive factors are generally weaker for those with a pathogenic BRCA variant than those without, while associations with modifiable lifestyle factors are similar for carriers and noncarriers. Impact: Advice to maintain a healthy weight, be physically active, and refrain from smoking will therefore benefit BRCA carriers as well as noncarriers.

Background Polygenic risk score (PRS), calculated based on genome-wide association studies (GWASs), can improve breast cancer (BC) risk assessment. To date, most BC GWASs have been performed in individuals of European (EUR) ancestry, and the generalisation of EUR-based PRS to other populations is a major challenge. In this study, we examined the performance of EUR-based BC PRS models in Ashkenazi Jewish (AJ) women. Methods We generated PRSs based on data on EUR women from the Breast Cancer Association Consortium (BCAC). We tested the performance of the PRSs in a cohort of 2161 AJ women from Israel (1437 cases and 724 controls) from BCAC (BCAC cohort from Israel (BCAC-IL)). In addition, we tested the performance of these EUR-based BC PRSs, as well as the established 313-SNP EUR BC PRS, in an independent cohort of 181 AJ women from Hadassah Medical Center (HMC) in Israel. Results In the BCAC-IL cohort, the highest OR per 1 SD was 1.56 (±0.09). The OR for AJ women at the top 10% of the PRS distribution compared with the middle quintile was 2.10 (±0.24). In the HMC cohort, the OR per 1 SD of the EUR-based PRS that performed best in the BCAC-IL cohort was 1.58±0.27. The OR per 1 SD of the commonly used 313-SNP BC PRS was 1.64 (±0.28). Conclusions Extant EUR GWAS data can be used for generating PRSs that identify AJ women with markedly elevated risk of BC and therefore hold promise for improving BC risk assessment in AJ women.

Abstract Background Breast cancer (BC) patients with a germline CHEK2 c.1100delC variant have an increased risk of contralateral BC (CBC) and worse BC‐specific survival (BCSS) compared to non‐carriers. Aim To assessed the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS. Methods Analyses were based on 82,701 women diagnosed with a first primary invasive BC including 963 CHEK2 c.1100delC carriers; median follow‐up was 9.1 years. Differential associations with treatment by CHEK2 c.1100delC status were tested by including interaction terms in a multivariable Cox regression model. A multi‐state model was used for further insight into the relation between CHEK2 c.1100delC status, treatment, CBC risk and death. Results There was no evidence for differential associations of therapy with CBC risk by CHEK2 c.1100delC status. The strongest association with reduced CBC risk was observed for the combination of chemotherapy and endocrine therapy [HR (95% CI): 0.66 (0.55–0.78)]. No association was observed with radiotherapy. Results from the multi‐state model showed shorter BCSS for CHEK2 c.1100delC carriers versus non‐carriers also after accounting for CBC occurrence [HR (95% CI): 1.30 (1.09–1.56)]. Conclusion Systemic therapy was associated with reduced CBC risk irrespective of CHEK2 c.1100delC status. Moreover, CHEK2 c.1100delC carriers had shorter BCSS, which appears not to be fully explained by their CBC risk.

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.