Risk Factors for Ovarian Cancer by BRCA Status: A Collaborative Case-Only Analysis

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.

Epidemiologic Factors in Relation to Colorectal Cancer Risk and Survival by Genotoxic Colibactin Mutational Signature

BACKGROUND

The genotoxin colibactin causes a tumor single-base substitution (SBS) mutational signature, SBS88. It is unknown whether epidemiologic factors' association with colorectal cancer risk and survival differs by SBS88.

METHODS

Within the Genetic Epidemiology of Colorectal Cancer Consortium and Colon Cancer Family Registry, we measured SBS88 in 4,308 microsatellite stable/microsatellite instability low tumors. Associations of epidemiologic factors with colorectal cancer risk by SBS88 were assessed using multinomial regression (N = 4,308 cases, 14,192 controls; cohort-only cases N = 1,911), and with colorectal cancer-specific survival using Cox proportional hazards regression (N = 3,465 cases).

RESULTS

392 (9%) tumors were SBS88 positive. Among all cases, the highest quartile of fruit intake was associated with lower risk of SBS88-positive colorectal cancer than SBS88-negative colorectal cancer [odds ratio (OR) = 0.53, 95% confidence interval (CI) 0.37-0.76; OR = 0.75, 95% CI 0.66-0.85, respectively, Pheterogeneity = 0.047]. Among cohort studies, associations of body mass index (BMI), alcohol, and fruit intake with colorectal cancer risk differed by SBS88. BMI ≥30 kg/m2 was associated with worse colorectal cancer-specific survival among those SBS88-positive [hazard ratio (HR) = 3.40, 95% CI 1.47-7.84], but not among those SBS88-negative (HR = 0.97, 95% CI 0.78-1.21, Pheterogeneity = 0.066).

CONCLUSIONS

Most epidemiologic factors did not differ by SBS88 for colorectal cancer risk or survival. Higher BMI may be associated with worse colorectal cancer-specific survival among those SBS88-positive; however, validation is needed in samples with whole-genome or whole-exome sequencing available.

IMPACT

This study highlights the importance of identification of tumor phenotypes related to colorectal cancer and understanding potential heterogeneity for risk and survival.

Genome-wide association analyses of ovarian cancer patients undergoing primary debulking surgery identify candidate genes for residual disease

Survival from ovarian cancer depends on the resection status after primary surgery. We performed genome-wide association analyses for resection status of 7705 ovarian cancer patients, including 4954 with high-grade serous carcinoma (HGSOC), to identify variants associated with residual disease. The most significant association with resection status was observed for rs72845444, upstream of MGMT, in HGSOC (p = 3.9 × 10-8). In gene-based analyses, PPP2R5C was the most strongly associated gene in HGSOC after stage adjustment. In an independent set of 378 ovarian tumours from the AGO-OVAR 11 study, variants near MGMT and PPP2R5C correlated with methylation and transcript levels, and PPP2R5C mRNA levels predicted progression-free survival in patients with residual disease. MGMT encodes a DNA repair enzyme, and PPP2R5C encodes the B56γ subunit of the PP2A tumour suppressor. Our results link heritable variation at these two loci with resection status in HGSOC.

Large-scale genome-wide association study of 398,238 women unveils seven novel loci associated with high-grade serous epithelial ovarian cancer risk

Background

Nineteen genomic regions have been associated with high-grade serous ovarian cancer (HGSOC). We used data from the Ovarian Cancer Association Consortium (OCAC), Consortium of Investigators of Modifiers of BRCA1/BRCA2 (CIMBA), UK Biobank (UKBB), and FinnGen to identify novel HGSOC susceptibility loci and develop polygenic scores (PGS).

Methods

We analyzed >22 million variants for 398,238 women. Associations were assessed separately by consortium and meta-analysed. OCAC and CIMBA data were used to develop PGS which were trained on FinnGen data and validated in UKBB and BioBank Japan.

Results

Eight novel variants were associated with HGSOC risk. An interesting discovery biologically was finding that TP53 3'-UTR SNP rs78378222 was associated with HGSOC (per T allele relative risk (RR)=1.44, 95%CI:1.28-1.62, P=1.76×10-9). The optimal PGS included 64,518 variants and was associated with an odds ratio of 1.46 (95%CI:1.37-1.54) per standard deviation in the UKBB validation (AUROC curve=0.61, 95%CI:0.59-0.62).

Conclusions

This study represents the largest GWAS for HGSOC to date. The results highlight that improvements in imputation reference panels and increased sample sizes can identify HGSOC associated variants that previously went undetected, resulting in improved PGS. The use of updated PGS in cancer risk prediction algorithms will then improve personalized risk prediction for HGSOC.

Genotoxic colibactin mutational signature in colorectal cancer is associated with clinicopathological features, specific genomic alterations and better survival

Background and Aims

The microbiome has long been suspected of a role in colorectal cancer (CRC) tumorigenesis. The mutational signature SBS88 mechanistically links CRC development with the strain of Escherichia coli harboring the pks island that produces the genotoxin colibactin, but the genomic, pathological and survival characteristics associated with SBS88-positive tumors are unknown.

Methods

SBS88-positive CRCs were identified from targeted sequencing data from 5,292 CRCs from 17 studies and tested for their association with clinico-pathological features, oncogenic pathways, genomic characteristics and survival.

Results

In total, 7.5% (398/5,292) of the CRCs were SBS88-positive, of which 98.7% (392/398) were microsatellite stable/microsatellite instability low (MSS/MSI-L), compared with 80% (3916/4894) of SBS88 negative tumors (p=1.5x10-28). Analysis of MSS/MSI-L CRCs demonstrated that SBS88 positive CRCs were associated with the distal colon (OR=1.84, 95% CI=1.40-2.42, p=1x10-5) and rectum (OR=1.90, 95% CI=1.44-2.51, p=6x10-6) tumor sites compared with the proximal colon. The top seven recurrent somatic mutations associated with SBS88-positive CRCs demonstrated mutational contexts associated with colibactin-induced DNA damage, the strongest of which was the APC:c.835-8A>G mutation (OR=65.5, 95%CI=39.0-110.0, p=3x10-80). Large copy number alterations (CNAs) including CNA loss on 14q and gains on 13q, 16q and 20p were significantly enriched in SBS88-positive CRCs. SBS88-positive CRCs were associated with better CRC-specific survival (p=0.007; hazard ratio of 0.69, 95% CI=0.52-0.90) when stratified by age, sex, study, and by stage.

Conclusion

SBS88-positivity, a biomarker of colibactin-induced DNA damage, can identify a novel subtype of CRC characterized by recurrent somatic mutations, copy number alterations and better survival. These findings provide new insights for treatment and prevention strategies for this subtype of CRC.

Concurrent RB1 loss and BRCA-deficiency predicts enhanced immunological response and long-term survival in tubo-ovarian high-grade serous carcinoma

Background

Somatic loss of the tumour suppressor RB1 is a common event in tubo-ovarian high-grade serous carcinoma (HGSC), which frequently co-occurs with alterations in homologous recombination DNA repair genes including BRCA1 and BRCA2 (BRCA). We examined whether tumour expression of RB1 was associated with survival across ovarian cancer histotypes (HGSC, endometrioid (ENOC), clear cell (CCOC), mucinous (MOC), low-grade serous carcinoma (LGSC)), and how co-occurrence of germline BRCA pathogenic variants and RB1 loss influences long-term survival in a large series of HGSC.

Patients and methods

RB1 protein expression patterns were classified by immunohistochemistry in epithelial ovarian carcinomas of 7436 patients from 20 studies participating in the Ovarian Tumor Tissue Analysis consortium and assessed for associations with overall survival (OS), accounting for patient age at diagnosis and FIGO stage. We examined RB1 expression and germline BRCA status in a subset of 1134 HGSC, and related genotype to survival, tumour infiltrating CD8+ lymphocyte counts and transcriptomic subtypes. Using CRISPR-Cas9, we deleted RB1 in HGSC cell lines with and without BRCA1 mutations to model co-loss with treatment response. We also performed genomic analyses on 126 primary HGSC to explore the molecular characteristics of concurrent homologous recombination deficiency and RB1 loss.

Results

RB1 protein loss was most frequent in HGSC (16.4%) and was highly correlated with RB1 mRNA expression. RB1 loss was associated with longer OS in HGSC (hazard ratio [HR] 0.74, 95% confidence interval [CI] 0.66-0.83, P = 6.8 ×10-7), but with poorer prognosis in ENOC (HR 2.17, 95% CI 1.17-4.03, P = 0.0140). Germline BRCA mutations and RB1 loss co-occurred in HGSC (P < 0.0001). Patients with both RB1 loss and germline BRCA mutations had a superior OS (HR 0.38, 95% CI 0.25-0.58, P = 5.2 ×10-6) compared to patients with either alteration alone, and their median OS was three times longer than non-carriers whose tumours retained RB1 expression (9.3 years vs. 3.1 years). Enhanced sensitivity to cisplatin (P < 0.01) and paclitaxel (P < 0.05) was seen in BRCA1 mutated cell lines with RB1 knockout. Among 126 patients with whole-genome and transcriptome sequence data, combined RB1 loss and genomic evidence of homologous recombination deficiency was correlated with transcriptional markers of enhanced interferon response, cell cycle deregulation, and reduced epithelial-mesenchymal transition in primary HGSC. CD8+ lymphocytes were most prevalent in BRCA-deficient HGSC with co-loss of RB1.

Conclusions

Co-occurrence of RB1 loss and BRCA mutation was associated with exceptionally long survival in patients with HGSC, potentially due to better treatment response and immune stimulation.

The Prognostic Effect of Immune Cell Infiltration Depends on Molecular Subtype in Endometrioid Ovarian Carcinomas

PURPOSE

Endometrioid ovarian carcinoma (ENOC) is the second most-common type of ovarian carcinoma, comprising 10%-20% of cases. Recently, the study of ENOC has benefitted from comparisons to endometrial carcinomas including defining ENOC with four prognostic molecular subtypes. Each subtype suggests differential mechanisms of progression, although tumor-initiating events remain elusive. There is evidence that the ovarian microenvironment may be critical to early lesion establishment and progression. However, while immune infiltrates have been well studied in high-grade serous ovarian carcinoma, studies in ENOC are limited.

EXPERIMENTAL DESIGN

We report on 210 ENOC, with clinical follow-up and molecular subtype annotation. Using multiplex IHC and immunofluorescence, we examine the prevalence of T-cell lineage, B-cell lineage, macrophages, and populations with programmed cell death protein 1 or programmed death-ligand 1 across subtypes of ENOC.

RESULTS

Immune cell infiltrates in tumor epithelium and stroma showed higher densities in ENOC subtypes with known high mutation burden (POLEmut and MMRd). While molecular subtypes were prognostically significant, immune infiltrates were not (overall survival P > 0.2). Analysis by molecular subtype revealed that immune cell density was prognostically significant in only the no specific molecular profile (NSMP) subtype, where immune infiltrates lacking B cells (TILB minus) had inferior outcome (disease-specific survival: HR, 4.0; 95% confidence interval, 1.1-14.7; P < 0.05). Similar to endometrial carcinomas, molecular subtype stratification was generally superior to immune response in predicting outcomes.

CONCLUSIONS

Subtype stratification is critical for better understanding of ENOC, in particular the distribution and prognostic significance of immune cell infiltrates. The role of B cells in the immune response within NSMP tumors warrants further study.

Batch-effect correction with sample remeasurement in highly confounded case-control studies

Batch effects are pervasive in biomedical studies. One approach to address the batch effects is repeatedly measuring a subset of samples in each batch. These remeasured samples are used to estimate and correct the batch effects. However, rigorous statistical methods for batch-effect correction with remeasured samples are severely underdeveloped. Here we developed a framework for batch-effect correction using remeasured samples in highly confounded case-control studies. We provided theoretical analyses of the proposed procedure, evaluated its power characteristics and provided a power calculation tool to aid in the study design. We found that the number of samples that need to be remeasured depends strongly on the between-batch correlation. When the correlation is high, remeasuring a small subset of samples is possible to rescue most of the power.

Multi-trait genome-wide association study identifies a novel endometrial cancer risk locus that associates with testosterone levels

To detect novel endometrial cancer risk variants, we leveraged information from endometrial cancer risk factors in a multi-trait GWAS analysis. We first assessed causal relationships between established and suspected endometrial cancer risk factors, and endometrial cancer using Mendelian randomization. Following multivariable analysis, five independent risk factors (waist circumference, testosterone levels, sex hormone binding globulin levels, age at menarche, and age at natural menopause) were included in a multi-trait Bayesian GWAS analysis. We identified three potentially novel loci that associate with endometrial cancer risk, one of which (7q22.1) replicated in an independent endometrial cancer GWAS dataset and was genome-wide significant in a meta-analysis. This locus may affect endometrial cancer risk through altered testosterone levels. Consistent with this, we observed colocalization between the signals for endometrial cancer risk and expression of CYP3A7, a gene involved in testosterone metabolism. Thus, our findings suggest opportunities for hormone therapy to prevent or treat endometrial cancer.

Lifetime ovulatory years and risk of epithelial ovarian cancer: a multinational pooled analysis

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.

Clinical parameters affecting survival outcomes in patients with low-grade serous ovarian carcinoma: an international multicentre analysis

BACKGROUND

Women with low-grade ovarian serous carcinoma (LGSC) benefit from surgical treatment; however, the role of chemotherapy is controversial. We examined an international database through the Ovarian Cancer Association Consortium to identify factors that affect survival in LGSC.

METHODS

We performed a retrospective cohort analysis of patients with LGSC who had had primary surgery and had overall survival data available. We performed univariate and multivariate analyses of progression-free survival and overall survival, and generated Kaplan-Meier survival curves.

RESULTS

Of the 707 patients with LGSC, 680 (96.2%) had available overall survival data. The patients' median age overall was 54 years. Of the 659 patients with International Federation of Obstetrics and Gynecology stage data, 156 (23.7%) had stage I disease, 64 (9.7%) had stage II, 395 (59.9%) had stage III, and 44 (6.7%) had stage IV. Of the 377 patients with surgical data, 200 (53.0%) had no visible residual disease. Of the 361 patients with chemotherapy data, 330 (91.4%) received first-line platinum-based chemotherapy. The median follow-up duration was 5.0 years. The median progression-free survival and overall survival were 43.2 months and 110.4 months, respectively. Multivariate analysis indicated a statistically significant impact of stage and residual disease on progression-free survival and overall survival. Platinum-based chemotherapy was not associated with a survival advantage.

CONCLUSION

This multicentre analysis indicates that complete surgical cytoreduction to no visible residual disease has the most impact on improved survival in LGSC. This finding could immediately inform and change practice.

p53 and ovarian carcinoma survival: an Ovarian Tumor Tissue Analysis consortium study

Our 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.

Methylation Signature Implicated in Immuno-Suppressive Activities in Tubo-Ovarian High-Grade Serous Carcinoma

BACKGROUND

Better understanding of prognostic factors in tubo-ovarian high-grade serous carcinoma (HGSC) is critical, as diagnosis confers an aggressive disease course. Variation in tumor DNA methylation shows promise predicting outcome, yet prior studies were largely platform-specific and unable to evaluate multiple molecular features.

METHODS

We analyzed genome-wide DNA methylation in 1,040 frozen HGSC, including 325 previously reported upon, seeking a multi-platform quantitative methylation signature that we evaluated in relation to clinical features, tumor characteristics, time to recurrence/death, extent of CD8+ tumor-infiltrating lymphocytes (TIL), gene expression molecular subtypes, and gene expression of the ATP-binding cassette transporter TAP1.

RESULTS

Methylation signature was associated with shorter time to recurrence, independent of clinical factors (N = 715 new set, hazard ratio (HR), 1.65; 95% confidence interval (CI), 1.10-2.46; P = 0.015; N = 325 published set HR, 2.87; 95% CI, 2.17-3.81; P = 2.2 × 10-13) and remained prognostic after adjustment for gene expression molecular subtype and TAP1 expression (N = 599; HR, 2.22; 95% CI, 1.66-2.95; P = 4.1 × 10-8). Methylation signature was inversely related to CD8+ TIL levels (P = 2.4 × 10-7) and TAP1 expression (P = 0.0011) and was associated with gene expression molecular subtype (P = 5.9 × 10-4) in covariate-adjusted analysis.

CONCLUSIONS

Multi-center analysis identified a novel quantitative tumor methylation signature of HGSC applicable to numerous commercially available platforms indicative of shorter time to recurrence/death, adjusting for other factors. Along with immune cell composition analysis, these results suggest a role for DNA methylation in the immunosuppressive microenvironment.

IMPACT

This work aids in identification of targetable epigenome processes and stratification of patients for whom tailored treatment may be most beneficial.

CCNE1 and survival of patients with tubo-ovarian high-grade serous carcinoma: An Ovarian Tumor Tissue Analysis consortium study

BACKGROUND

Cyclin 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.

METHODS

Within the Ovarian Tumor Tissue Analysis consortium, amplification status and protein level in 3029 HGSC cases and mRNA expression in 2419 samples were investigated.

RESULTS

High-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.

CONCLUSION

This study provides large-scale validation that CCNE1 high-level amplification is associated with shorter survival, supporting its utility as a prognostic biomarker in HGSC.

Lifestyle and personal factors associated with having macroscopic residual disease after ovarian cancer primary cytoreductive surgery

OBJECTIVE

The presence of macroscopic residual disease after primary cytoreductive surgery (PCS) is an important factor influencing survival for patients with high-grade serous ovarian cancer (HGSC). More research is needed to identify factors associated with having macroscopic residual disease. We analyzed 12 lifestyle and personal exposures known to be related to ovarian cancer risk or inflammation to identify those associated with having residual disease after surgery.

METHODS

This analysis used data on 2054 patients with advanced stage HGSC from the Ovarian Cancer Association Consortium. The exposures were body mass index, breastfeeding, oral contraceptive use, depot-medroxyprogesterone acetate use, endometriosis, first-degree family history of ovarian cancer, incomplete pregnancy, menopausal hormone therapy use, menopausal status, parity, smoking, and tubal ligation. Logistic regression models were fit to assess the association between these exposures and having residual disease following PCS.

RESULTS

Menopausal estrogen-only therapy (ET) use was associated with 33% lower odds of having macroscopic residual disease compared to never use (OR = 0.67, 95%CI 0.46-0.97, p = 0.033). Compared to nulliparous women, parous women who did not breastfeed had 36% lower odds of having residual disease (OR = 0.64, 95%CI 0.43-0.94, p = 0.022), while there was no association among parous women who breastfed (OR = 0.90, 95%CI 0.65-1.25, p = 0.53).

CONCLUSIONS

The association between ET and having no macroscopic residual disease is plausible given a strong underlying biologic hypothesis between this exposure and diagnosis with HGSC. If this or the parity finding is replicated, these factors could be included in risk stratification models to determine whether HGSC patients should receive PCS or neoadjuvant chemotherapy.

Increased FOXJ1 protein expression is associated with improved overall survival in high-grade serous ovarian carcinoma: an Ovarian Tumor Tissue Analysis Consortium Study

BACKGROUND

Recently, we showed a >60% difference in 5-year survival for patients with tubo-ovarian high-grade serous carcinoma (HGSC) when stratified by a 101-gene mRNA expression prognostic signature. Given the varied patient outcomes, this study aimed to translate prognostic mRNA markers into protein expression assays by immunohistochemistry and validate their survival association in HGSC.

METHODS

Two prognostic genes, FOXJ1 and GMNN, were selected based on high-quality antibodies, correlation with protein expression and variation in immunohistochemical scores in a preliminary cohort (n = 134 and n = 80, respectively). Six thousand four hundred and thirty-four (FOXJ1) and 5470 (GMNN) formalin-fixed, paraffin-embedded ovarian neoplasms (4634 and 4185 HGSC, respectively) represented on tissue microarrays from the Ovarian Tumor Tissue Analysis consortium underwent immunohistochemical staining and scoring, then univariate and multivariate survival analysis.

RESULTS

Consistent with mRNA, FOXJ1 protein expression exhibited a linear, increasing association with improved overall survival in HGSC patients. Women with >50% expression had the most favourable outcomes (HR = 0.78, 95% CI 0.67-0.91, p < 0.0001). GMNN protein expression was not significantly associated with overall HSGC patient survival. However, HGSCs with >35% GMNN expression showed a trend for better outcomes, though this was not significant.

CONCLUSION

We provide foundational evidence for the prognostic value of FOXJ1 in HGSC, validating the prior mRNA-based prognostic association by immunohistochemistry.

The genomic and immune landscape of long-term survivors of high-grade serous ovarian cancer

Fewer than half of all patients with advanced-stage high-grade serous ovarian cancers (HGSCs) survive more than five years after diagnosis, but those who have an exceptionally long survival could provide insights into tumor biology and therapeutic approaches. We analyzed 60 patients with advanced-stage HGSC who survived more than 10 years after diagnosis using whole-genome sequencing, transcriptome and methylome profiling of their primary tumor samples, comparing this data to 66 short- or moderate-term survivors. Tumors of long-term survivors were more likely to have multiple alterations in genes associated with DNA repair and more frequent somatic variants resulting in an increased predicted neoantigen load. Patients clustered into survival groups based on genomic and immune cell signatures, including three subsets of patients with BRCA1 alterations with distinctly different outcomes. Specific combinations of germline and somatic gene alterations, tumor cell phenotypes and differential immune responses appear to contribute to long-term survival in HGSC.

Abstract PR008: Evaluation of intra-tumoral pathogenic bacteria pks+ E. coli, enterotoxigenic B. fragilis and fusobacterium nucleatum in 3695 colorectal cancer cases

Background: The microbiome in colorectal cancer (CRC) tumors is now recognized as a potent facilitator of tumorigenesis. Identification of the intra-tumoral microbiome may provide an avenue to enhance prediction of adverse outcomes and ultimately lead to improved survival of CRC, the second leading cause of cancer death. Methods: We performed targeted sequencing of tumor and matched normal DNA samples for 3695 colorectal cancer cases within the Genetics and Colorectal Cancer Consortium (GECCO) and the Colon Cancer Family Registry (CCFR). Tumor DNA was extracted primarily from formalin fixed paraffin embedded (FFPE) or from fresh frozen tissue. FFPE slides were dissected to ensure &gt;70% tumor content. Normal DNA was primarily extracted from blood and otherwise normal surrounding colorectal tissue. We designed probes for bacterial pathogens including Bacteroides fragilis (ETBF), polyketide synthase positive E. coli (pks+ EC), enterotoxigenic and Fusobacterium nucleatum (Fn). Logistic regression was used to assess tumor-associated bacteria and various clinical parameters, tumor molecular characteristics, and tumor stage adjusted for age, sex, and study and a Cox promotional hazard to assess association with CRC-specific survival adjusted for age at diagnosis, sex, mutational burden and microsatellite instability (MSI) and stratified baseline hazards by study. Survival data were available for 2357 patients. Results: The prevalence of patients with bacteria in tumors varied for ETBF (5.2%), pks+ EC (12.9%), and Fn (10.4%). Men had a lower probability of tumors with Fn (odds ratio (OR)=0.8, 95% confidence interval (CI)=0.64-0.99,p=0.04) but higher odds of tumors with pks+ EC (OR=1.38, CI=1.14-1.68,p=0.001). Patients with proximal tumors had higher odds of carrying Fn (OR=1.76, CI=1.38-2.26, p&lt;0.001) as did patients with MSI high tumors (OR 2.79, CI=2.08-3.74) p&lt;0.001). Patients with MSI-high tumors had lower odds of carrying pks+ EC (OR=0.63, CI=0.44-0.88, p=0.008) when compared to patients with microsatellite stable tumors. Among non-hypermutated tumors, we observed that tumors positive for pks+ EC were associated with better CRC-specific survival (OR=0.73, CI=0.57-0.93,p=0.01 whereas in patients with ETBF positive tumors had shorter survival (OR=1.99,CI=1,18-3.36,p=0.009) as did patients with tumors positive for Fn (OR=1.33, CI=1.04-1.72, p=0.03). Conclusions: In this largest investigation of intra-tumoral bacteria in CRC, we demonstrate that inclusion of pathogenic intra-tumoral microbiome may enhance prediction of survival and can identify a subset of patients with CRC that may benefit from targeted therapies that modulate the tumor microbiome.

Citation Format: Meredith A. J. Hullar, Keith R. Curtis, Tabitha A. Harrison, Yi Lin, Robert S. Steinfelder, Sonja I. Berndt, Daniel D. Buchanan, Andrew T. Chan, David Drew, Jane C. Figueiredo, Amy J. French, Ellen L. Goode, Mark A. Jenkins, Yohannes A. Melaku, Victor Moreno, Shuji Ogino, Stephen N. Thibideau, Amanda I. Phipps, Ulrike Peters. Evaluation of intra-tumoral pathogenic bacteria pks+ E. coli, enterotoxigenic B. fragilis and fusobacterium nucleatum in 3695 colorectal cancer cases [abstract]. In: Proceedings of the AACR Special Conference on Colorectal Cancer; 2022 Oct 1-4; Portland, OR. Philadelphia (PA): AACR; Cancer Res 2022;82(23 Suppl_1):Abstract nr PR008.

The proteome of clear cell ovarian carcinoma

Clear cell ovarian carcinoma (CCOC) is the second most common subtype of epithelial ovarian carcinoma. Late-stage CCOC is not responsive to gold-standard chemotherapy and results in suboptimal outcomes for patients. In-depth molecular insight is urgently needed to stratify the disease and drive therapeutic development. We conducted global proteomics for 192 cases of CCOC and compared these with other epithelial ovarian carcinoma subtypes. Our results showed distinct proteomic differences in CCOC compared with other epithelial ovarian cancer subtypes including alterations in lipid and purine metabolism pathways. Furthermore, we report potential clinically significant proteomic subgroups within CCOC, suggesting the biologic plausibility of stratified treatment for this cancer. Taken together, our results provide a comprehensive understanding of the CCOC proteomic landscape to facilitate future understanding and research of this disease. © 2022 The Pathological Society of Great Britain and Ireland.

Molecular Subclasses of Clear Cell Ovarian Carcinoma and Their Impact on Disease Behavior and Outcomes

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.

Copy Number Variants Are Ovarian Cancer Risk Alleles at Known and Novel Risk Loci

BACKGROUND

Known 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.

METHODS

Single 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.

RESULTS

We 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.

CONCLUSIONS

CNVs 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.

Cross-Ancestry Investigation of Venous Thromboembolism Genomic Predictors

BACKGROUND

Venous thromboembolism (VTE) is a life-threatening vascular event with environmental and genetic determinants. Recent VTE genome-wide association studies (GWAS) meta-analyses involved nearly 30 000 VTE cases and identified up to 40 genetic loci associated with VTE risk, including loci not previously suspected to play a role in hemostasis. The aim of our research was to expand discovery of new genetic loci associated with VTE by using cross-ancestry genomic resources.

METHODS

We present new cross-ancestry meta-analyzed GWAS results involving up to 81 669 VTE cases from 30 studies, with replication of novel loci in independent populations and loci characterization through in silico genomic interrogations.

RESULTS

In our genetic discovery effort that included 55 330 participants with VTE (47 822 European, 6320 African, and 1188 Hispanic ancestry), we identified 48 novel associations, of which 34 were replicated after correction for multiple testing. In our combined discovery-replication analysis (81 669 VTE participants) and ancestry-stratified meta-analyses (European, African, and Hispanic), we identified another 44 novel associations, which are new candidate VTE-associated loci requiring replication. In total, across all GWAS meta-analyses, we identified 135 independent genomic loci significantly associated with VTE risk. A genetic risk score of the significantly associated loci in Europeans identified a 6-fold increase in risk for those in the top 1% of scores compared with those with average scores. We also identified 31 novel transcript associations in transcriptome-wide association studies and 8 novel candidate genes with protein quantitative-trait locus Mendelian randomization analyses. In silico interrogations of hemostasis and hematology traits and a large phenome-wide association analysis of the 135 GWAS loci provided insights to biological pathways contributing to VTE, with some loci contributing to VTE through well-characterized coagulation pathways and others providing new data on the role of hematology traits, particularly platelet function. Many of the replicated loci are outside of known or currently hypothesized pathways to thrombosis.

CONCLUSIONS

Our cross-ancestry GWAS meta-analyses identified new loci associated with VTE. These findings highlight new pathways to thrombosis and provide novel molecules that may be useful in the development of improved antithrombosis treatments.

Computational tumor stroma reaction evaluation led to novel prognosis-associated fibrosis and molecular signature discoveries in high-grade serous ovarian carcinoma

Background

As one of the key criteria to differentiate benign vs. malignant tumors in ovarian and other solid cancers, tumor-stroma reaction (TSR) is long observed by pathologists and has been found correlated with patient prognosis. However, paucity of study aims to overcome subjective bias or automate TSR evaluation for enabling association analysis to a large cohort.

Materials and methods

Serving as positive and negative sets of TSR studies, H&E slides of primary tumors of high-grade serous ovarian carcinoma (HGSOC) (n = 291) and serous borderline ovarian tumor (SBOT) (n = 15) were digitally scanned. Three pathologist-defined quantification criteria were used to characterize the extents of TSR. Scores for each criterion were annotated (0/1/2 as none-low/intermediate/high) in the training set consisting of 18,265 H&E patches. Serial of deep learning (DL) models were trained to identify tumor vs. stroma regions and predict TSR scores. After cross-validation and independent validations, the trained models were generalized to the entire HGSOC cohort and correlated with clinical characteristics. In a subset of cases tumor transcriptomes were available, gene- and pathway-level association studies were conducted with TSR scores.

Results

The trained models accurately identified the tumor stroma tissue regions and predicted TSR scores. Within tumor stroma interface region, TSR fibrosis scores were strongly associated with patient prognosis. Cancer signaling aberrations associated 14 KEGG pathways were also found positively correlated with TSR-fibrosis score.

Conclusion

With the aid of DL, TSR evaluation could be generalized to large cohort to enable prognostic association analysis and facilitate discovering novel gene and pathways associated with disease progress.

Into the future: A pilot study combining imaging with molecular profiling to predict resectability in ovarian cancer

OBJECTIVE

We sought to determine the predictive value of combining tumor molecular subtype and computerized tomography (CT) imaging for surgical outcomes after primary cytoreductive surgery in advanced stage high-grade serous ovarian cancer (HGSOC) patients.

METHODS

We identified 129 HGSOC patients who underwent pre-operative CT imaging and post-operative tumor mRNA profiling. A continuous CT-score indicative of overall disease burden was defined based on six imaging measurements of anatomic involvement. Molecular subtypes were derived from mRNA profiling of chemo-naïve tumors and classified as mesenchymal (MES) subtype (36%) or non-MES subtype (64%). Fischer exact tests and multivariate logistic regression examined residual disease and surgical complexity.

RESULTS

Women with higher CT-scores were more likely to have MES subtype tumors (p = 0.014). MES subtypes and a high CT-score were independently predictive of macroscopic disease and high surgical complexity. In multivariate models adjusting for age, stage and American Society of Anesthesiologists (ASA) score, patients with a MES subtype and high CT-score had significantly elevated risk of macroscopic disease (OR = 26.7, 95% CI = [6.42, 187]) and were more likely to undergo high complexity surgery (OR = 9.53, 95% CI = [2.76, 40.6], compared to patients with non-MES tumor and low CT-score.

CONCLUSION

Preoperative CT imaging combined with tumor molecular subtyping can identify a subset of women unlikely to have resectable disease and likely to require high complexity surgery. Along with other clinical factors, these may refine predictive scores for resection and assist treatment planning. Investigating methods for pre-surgical molecular subtyping is an important next step.

Abstract 2248: Environmental factors associated with residual disease after ovarian cancer primary cytoreduction surgery

Background: Ovarian cancer is the deadliest gynecologic cancer. Standard treatments for advanced stage high-grade serous ovarian cancer, the most common type, include (1) primary cytoreductive surgery (PCS) followed by adjuvant chemotherapy or (2) neoadjuvant chemotherapy (NACT) with interval debulking surgery. Patients in whom PCS is unlikely to yield optimal cytoreduction to no visible residual disease (R0) or who have medical contraindications to PCS are recommended to undergo NACT. Previous studies on associations between environmental factors and risk of any macroscopic residual disease have yielded inconsistent results. We aimed to (1) comprehensively examine the associations between demographic, lifestyle and reproductive factors and risk of residual disease after PCS; and (2) develop and internally validate a risk prediction model based on these factors.

Methods: We used pooled data on 3,492 women who had PCS following a diagnosis with advanced stage high-grade serous ovarian, fallopian tube or primary peritoneal cancers from the Ovarian Cancer Association Consortium. Fifteen exposures of interest included age at diagnosis, menopausal status, race/ethnicity, education level, first-degree family history of ovarian cancer, endometriosis, smoking, body mass index, parity, incomplete pregnancy, tubal ligation, use of combined oral contraceptives, depot-medroxyprogesterone acetate, estrogen (ET), and combined estrogen-progestin therapy. We fit logistic regression models to examine each exposure-residual disease association in 80% of the data (n=2,794; random selection of participants). We developed a risk prediction model including the factors with p-values ≤0.2. Area under the receiver operating characteristic curve (AUC) was used to assess the model’s discrimination in the remaining 20% of the data (n=698).

Results: Of the 3,492 participants, 2,003 (57%) had residual disease following PCS. Older age at diagnosis was associated with an increased risk of residual disease (OR=1.07 per five years, 95% CI 1.03-1.11). In contrast, a family history of ovarian cancer (OR=0.63, 95% CI 0.42-0.95) or ET use for 5+ years (OR=0.61, 95% CI 0.39-0.96) were associated with achieving R0. These above factors were included in the risk prediction model as were race/ethnicity, personal history of endometriosis, and smoking, whose p-values ≤0.2. The model showed modest performance in the validation set (AUC=0.64).

Conclusions: Younger age at diagnosis, family history of ovarian cancer, and long-term ET use were associated with achieving R0 following PCS. Future studies incorporating genetic and clinical factors to improve the risk prediction for residual disease following PCS are warranted.

Citation Format: Minh Tung Phung, Andrew Berchuck, Ellen L. Goode, Marc T. Goodman, Gillian E. Hanley, Jean Richardson, Bronwyn Grout, Anne Chase, Cindy McKinnon Deurloo, Beth Y. Karlan, Toon Van Gorp, Keitaro Matsuo, Karen McLean, Malcolm C. Pike, Joellen M. Schildkraut, Kathryn L. Terry, Anna DeFazio, Penelope M. Webb, Paul D. P. Pharoah, Susan J. Ramus, Celeste Leigh Pearce. Environmental factors associated with residual disease after ovarian cancer primary cytoreduction surgery [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 2248.

Germline BRCA variants, lifestyle and ovarian cancer survival

OBJECTIVE

Women with ovarian cancer who have a pathogenic germline variant in BRCA1 or BRCA2 (BRCA) have been shown to have better 5-year survival after diagnosis than women who are BRCA-wildtype (non-carriers). Modifiable lifestyle factors, including smoking, physical activity and body mass index (BMI) have previously been associated with ovarian cancer survival; however, it is unknown whether these associations differ by germline BRCA status.

METHODS

We investigated measures of lifestyle prior to diagnosis in two cohorts of Australian women with invasive epithelial ovarian cancer, using Cox proportional hazards regression to calculate adjusted hazard ratios (HRs) and 95% confidence intervals (CIs).

RESULTS

In the combined studies (n = 1923), there was little association between physical activity, BMI or alcohol intake and survival, and no difference by BRCA status. However, the association between current smoking status before diagnosis and poorer survival was stronger for BRCA variant carriers (HR 1.98; 95% CI 1.20-3.27) than non-carriers (HR 1.18; 95% CI 0.96-1.46; p-interaction 0.02). We saw a similar differential association with smoking when we pooled results from two additional cohorts from the USA and UK (n = 2120). Combining the results from all four studies gave a pooled-HR of 1.94 (95% CI 1.28-2.94) for current smoking among BRCA variant carriers compared to 1.08 (0.90-1.29) for non-carriers.

CONCLUSIONS

Our results suggest that the adverse effect of smoking on survival may be stronger for women with a BRCA variant than those without. Thus, while smoking cessation may improve outcomes for all women with ovarian cancer, it might provide a greater benefit for BRCA variant carriers.

Multi-modal immune profiling of mucinous ovarian carcinoma: analysis from the Ovarian Tumor Tissue Analysis/Multidisciplinary Ovarian Cancer Outcomes Group consortia

Mucinous ovarian carcinoma (MOC) is a rare histotype of ovarian cancer. Women with advanced disease (FIGO stage III/IV) have very poor outcomes, with 5-year survival of &lt;20%. Little is known about the anti-tumor immune response in MOC, limiting potential immunotherapeutic options for patients.

Immune data from 3 platforms were analysed. A NanoString mRNA expression panel with 39 immune-related genes (n=60 MOC, n=4198 other ovarian histotypes), a NanoString mRNA plexset with PD-1, PD-L1 (n=241 MOC, n=115 upper and lower gastrointestinal (GI) cancers), and 2 multicolor immunofluorescence (mcIF) panels (n=121 MOC) for CD68, PD-L1, PD-1, CD8 and FOXP3. Data overlapped all platforms for 22 MOC. Immune marker expression by ovarian cancer histotype was measured. Uni- and multivariable Cox Proportional Hazards assessed overall survival (OS) within MOC. Hierarchical clustering looked for subgroups of patients based on co-expression of immune genes.

Compared to more common high-grade serous ovarian cancer (HGSC) and GI cancers, MOC were less immunologically active. High IGHM expression was associated with poorer OS (HR 1.31 (95% CI 1.02–1.69), p=0.033). High tumor densities of FOXP3 and PD-1 positive cells conferred poorer OS in a model with age, stage and site; however, these and PD-L1 and CD68 were associated with grade. The 39 gene NanoString and mcIF datasets revealed 4 clusters, roughly translating to immune cold, moderate and hot (groups 3,4). mcIF analysis revealed a subset with PD-L1 expression in the absence of T-regulatory cell infiltrates.

This study provides novel insights into the MOC immune landscape. Despite relatively low immune activity overall, immunotherapy could be explored for a subset of patients with MOC.

Supported by NSW Ministry of Health and UNSW Sydney, through the NSW Health PhD Scholarship Program 2017-2022; Translational Cancer Research Network top-up scholarship 2021

Validated biomarker assays confirm that ARID1A loss is confounded with MMR deficiency, CD8+ TIL infiltration, and provides no independent prognostic value in endometriosis-associated ovarian carcinomas

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.

MCM3 is a novel proliferation marker associated with longer survival for patients with tubo-ovarian high-grade serous carcinoma

Tubo-ovarian high-grade serous carcinomas (HGSC) are highly proliferative neoplasms that generally respond well to platinum/taxane chemotherapy. We recently identified minichromosome maintenance complex component 3 (MCM3), which is involved in the initiation of DNA replication and proliferation, as a favorable prognostic marker in HGSC. Our objective was to further validate whether MCM3 mRNA expression and possibly MCM3 protein levels are associated with survival in patients with HGSC. MCM3 mRNA expression was measured using NanoString expression profiling on formalin-fixed and paraffin-embedded tissue (N = 2355 HGSC) and MCM3 protein expression was assessed by immunohistochemistry (N = 522 HGSC) and compared with Ki-67. Kaplan-Meier curves and the Cox proportional hazards model were used to estimate associations with survival. Among chemotherapy-naïve HGSC, higher MCM3 mRNA expression (one standard deviation increase in the score) was associated with longer overall survival (HR = 0.87, 95% CI 0.81-0.92, p < 0.0001, N = 1840) in multivariable analysis. MCM3 mRNA expression was highest in the HGSC C5.PRO molecular subtype, although no interaction was observed between MCM3, survival and molecular subtypes. MCM3 and Ki-67 protein levels were significantly lower after exposure to neoadjuvant chemotherapy compared to chemotherapy-naïve tumors: 37.0% versus 46.4% and 22.9% versus 34.2%, respectively. Among chemotherapy-naïve HGSC, high MCM3 protein levels were also associated with significantly longer disease-specific survival (HR = 0.52, 95% CI 0.36-0.74, p = 0.0003, N = 392) compared to cases with low MCM3 protein levels in multivariable analysis. MCM3 immunohistochemistry is a promising surrogate marker of proliferation in HGSC.

Polygenic risk modeling for prediction of epithelial ovarian cancer risk

Polygenic risk scores (PRS) for epithelial ovarian cancer (EOC) have the potential to improve risk stratification. Joint estimation of Single Nucleotide Polymorphism (SNP) effects in models could improve predictive performance over standard approaches of PRS construction. Here, we implemented computationally efficient, penalized, logistic regression models (lasso, elastic net, stepwise) to individual level genotype data and a Bayesian framework with continuous shrinkage, "select and shrink for summary statistics" (S4), to summary level data for epithelial non-mucinous ovarian cancer risk prediction. We developed the models in a dataset consisting of 23,564 non-mucinous EOC cases and 40,138 controls participating in the Ovarian Cancer Association Consortium (OCAC) and validated the best models in three populations of different ancestries: prospective data from 198,101 women of European ancestries; 7,669 women of East Asian ancestries; 1,072 women of African ancestries, and in 18,915 BRCA1 and 12,337 BRCA2 pathogenic variant carriers of European ancestries. In the external validation data, the model with the strongest association for non-mucinous EOC risk derived from the OCAC model development data was the S4 model (27,240 SNPs) with odds ratios (OR) of 1.38 (95% CI: 1.28-1.48, AUC: 0.588) per unit standard deviation, in women of European ancestries; 1.14 (95% CI: 1.08-1.19, AUC: 0.538) in women of East Asian ancestries; 1.38 (95% CI: 1.21-1.58, AUC: 0.593) in women of African ancestries; hazard ratios of 1.36 (95% CI: 1.29-1.43, AUC: 0.592) in BRCA1 pathogenic variant carriers and 1.49 (95% CI: 1.35-1.64, AUC: 0.624) in BRCA2 pathogenic variant carriers. Incorporation of the S4 PRS in risk prediction models for ovarian cancer may have clinical utility in ovarian cancer prevention programs.

High Prediagnosis Inflammation-Related Risk Score Associated with Decreased Ovarian Cancer Survival

BACKGROUND

There is suggestive evidence that inflammation is related to ovarian cancer survival. However, more research is needed to identify inflammation-related factors that are associated with ovarian cancer survival and to determine their combined effects.

METHODS

This analysis used pooled data on 8,147 women with invasive epithelial ovarian cancer from the Ovarian Cancer Association Consortium. The prediagnosis inflammation-related exposures of interest included alcohol use; aspirin use; other nonsteroidal anti-inflammatory drug use; body mass index; environmental tobacco smoke exposure; history of pelvic inflammatory disease, polycystic ovarian syndrome, and endometriosis; menopausal hormone therapy use; physical inactivity; smoking status; and talc use. Using Cox proportional hazards models, the relationship between each exposure and survival was assessed in 50% of the data. A weighted inflammation-related risk score (IRRS) was developed, and its association with survival was assessed using Cox proportional hazards models in the remaining 50% of the data.

RESULTS

There was a statistically significant trend of increasing risk of death per quartile of the IRRS [HR = 1.09; 95% confidence interval (CI), 1.03-1.14]. Women in the upper quartile of the IRRS had a 31% higher death rate compared with the lowest quartile (95% CI, 1.11-1.54).

CONCLUSIONS

A higher prediagnosis IRRS was associated with an increased mortality risk after an ovarian cancer diagnosis. Further investigation is warranted to evaluate whether postdiagnosis exposures are also associated with survival.

IMPACT

Given that pre- and postdiagnosis exposures are often correlated and many are modifiable, our study results can ultimately motivate the development of behavioral recommendations to enhance survival among patients with ovarian cancer.

Functional analysis of the 1p34.3 risk locus implicates GNL2 in high-grade serous ovarian cancer

The high-grade serous ovarian cancer (HGSOC) risk locus at chromosome 1p34.3 resides within a frequently amplified genomic region signifying the presence of an oncogene. Here, we integrate in silico variant-to-function analysis with functional studies to characterize the oncogenic potential of candidate genes in the 1p34.3 locus. Fine mapping of genome-wide association statistics identified candidate causal SNPs local to H3K27ac-demarcated enhancer regions that exhibit allele-specific binding for CTCF in HGSOC and normal fallopian tube secretory epithelium cells (FTSECs). SNP risk associations colocalized with eQTL for six genes (DNALI1, GNL2, RSPO1, SNIP1, MEAF6, and LINC01137) that are more highly expressed in carriers of the risk allele, and three (DNALI1, GNL2, and RSPO1) were upregulated in HGSOC compared to normal ovarian surface epithelium cells and/or FTSECs. Increased expression of GNL2 and MEAF6 was associated with shorter survival in HGSOC with 1p34.3 amplifications. Despite its activation of β-catenin signaling, RSPO1 overexpression exerted no effects on proliferation or colony formation in our study of ovarian cancer and FTSECs. Instead, GNL2, MEAF6, and SNIP1 silencing impaired in vitro ovarian cancer cell growth. Additionally, GNL2 silencing diminished xenograft tumor formation, whereas overexpression stimulated proliferation and colony formation in FTSECs. GNL2 influences 60S ribosomal subunit maturation and global protein synthesis in ovarian cancer and FTSECs, providing a potential mechanism of how GNL2 upregulation might promote ovarian cancer development and mediate genetic susceptibility of HGSOC.

DNA Methylation Profiles of Ovarian Clear Cell Carcinoma

BACKGROUND

Ovarian clear cell carcinoma (OCCC) is a rare ovarian cancer histotype that tends to be resistant to standard platinum-based chemotherapeutics. We sought to better understand the role of DNA methylation in clinical and biological subclassification of OCCC.

METHODS

We interrogated genome-wide methylation using DNA from fresh frozen tumors from 271 cases, applied nonsmooth nonnegative matrix factorization (nsNMF) clustering, and evaluated clinical associations and biological pathways.

RESULTS

Two approximately equally sized clusters that associated with several clinical features were identified. Compared with Cluster 2 (N = 137), Cluster 1 cases (N = 134) presented at a more advanced stage, were less likely to be of Asian ancestry, and tended to have poorer outcomes including macroscopic residual disease following primary debulking surgery (P < 0.10). Subset analyses of targeted tumor sequencing and IHC data revealed that Cluster 1 tumors showed TP53 mutation and abnormal p53 expression, and Cluster 2 tumors showed aneuploidy and ARID1A/PIK3CA mutation (P < 0.05). Cluster-defining CpGs included 1,388 CpGs residing within 200 bp of the transcription start sites of 977 genes; 38% of these genes (N = 369 genes) were differentially expressed across cluster in transcriptomic subset analysis (P < 10-4). Differentially expressed genes were enriched for six immune-related pathways, including IFNα and IFNγ responses (P < 10-6).

CONCLUSIONS

DNA methylation clusters in OCCC correlate with disease features and gene expression patterns among immune pathways.

IMPACT

This work serves as a foundation for integrative analyses that better understand the complex biology of OCCC in an effort to improve potential for development of targeted therapeutics.

Identification of a Locus Near ULK1 Associated With Progression-Free Survival in Ovarian Cancer

BACKGROUND

Many loci have been found to be associated with risk of epithelial ovarian cancer (EOC). However, although there is considerable variation in progression-free survival (PFS), no loci have been found to be associated with outcome at genome-wide levels of significance.

METHODS

We carried out a genome-wide association study (GWAS) of PFS in 2,352 women with EOC who had undergone cytoreductive surgery and standard carboplatin/paclitaxel chemotherapy.

RESULTS

We found seven SNPs at 12q24.33 associated with PFS (P < 5 × 10-8), the top SNP being rs10794418 (HR = 1.24; 95% CI, 1.15-1.34; P = 1.47 × 10-8). High expression of a nearby gene, ULK1, is associated with shorter PFS in EOC, and with poor prognosis in other cancers. SNP rs10794418 is also associated with expression of ULK1 in ovarian tumors, with the allele associated with shorter PFS being associated with higher expression, and chromatin interactions were detected between the ULK1 promoter and associated SNPs in serous and endometrioid EOC cell lines. ULK1 knockout ovarian cancer cell lines showed significantly increased sensitivity to carboplatin in vitro.

CONCLUSIONS

The locus at 12q24.33 represents one of the first genome-wide significant loci for survival for any cancer. ULK1 is a plausible candidate for the target of this association.

IMPACT

This finding provides insight into genetic markers associated with EOC outcome and potential treatment options.See related commentary by Peres and Monteiro, p. 1604.

Pleiotropy-guided transcriptome imputation from normal and tumor tissues identifies candidate susceptibility genes for breast and ovarian cancer

Familial, sequencing, and genome-wide association studies (GWASs) and genetic correlation analyses have progressively unraveled the shared or pleiotropic germline genetics of breast and ovarian cancer. In this study, we aimed to leverage this shared germline genetics to improve the power of transcriptome-wide association studies (TWASs) to identify candidate breast cancer and ovarian cancer susceptibility genes. We built gene expression prediction models using the PrediXcan method in 681 breast and 295 ovarian tumors from The Cancer Genome Atlas and 211 breast and 99 ovarian normal tissue samples from the Genotype-Tissue Expression project and integrated these with GWAS meta-analysis data from the Breast Cancer Association Consortium (122,977 cases/105,974 controls) and the Ovarian Cancer Association Consortium (22,406 cases/40,941 controls). The integration was achieved through application of a pleiotropy-guided conditional/conjunction false discovery rate (FDR) approach in the setting of a TWASs. This identified 14 candidate breast cancer susceptibility genes spanning 11 genomic regions and 8 candidate ovarian cancer susceptibility genes spanning 5 genomic regions at conjunction FDR < 0.05 that were >1 Mb away from known breast and/or ovarian cancer susceptibility loci. We also identified 38 candidate breast cancer susceptibility genes and 17 candidate ovarian cancer susceptibility genes at conjunction FDR < 0.05 at known breast and/or ovarian susceptibility loci. The 22 genes identified by our cross-cancer analysis represent promising candidates that further elucidate the role of the transcriptome in mediating germline breast and ovarian cancer risk.

Frequent POLE-driven hypermutation in ovarian endometrioid cancer revealed by mutational signatures in RNA sequencing

BACKGROUND

DNA polymerase epsilon (POLE) is encoded by the POLE gene, and POLE-driven tumors are characterized by high mutational rates. POLE-driven tumors are relatively common in endometrial and colorectal cancer, and their presence is increasingly recognized in ovarian cancer (OC) of endometrioid type. POLE-driven cases possess an abundance of TCT > TAT and TCG > TTG somatic mutations characterized by mutational signature 10 from the Catalog of Somatic Mutations in Cancer (COSMIC). By quantifying the contribution of COSMIC mutational signature 10 in RNA sequencing (RNA-seq) we set out to identify POLE-driven tumors in a set of unselected Mayo Clinic OC.

METHODS

Mutational profiles were calculated using expressed single-nucleotide variants (eSNV) in the Mayo Clinic OC tumors (n = 195), The Cancer Genome Atlas (TCGA) OC tumors (n = 419), and the Genotype-Tissue Expression (GTEx) normal ovarian tissues (n = 84). Non-negative Matrix Factorization (NMF) of the mutational profiles inferred the contribution per sample of four distinct mutational signatures, one of which corresponds to COSMIC mutational signature 10.

RESULTS

In the Mayo Clinic OC cohort we identified six tumors with a predicted contribution from COSMIC mutational signature 10 of over five mutations per megabase. These six cases harbored known POLE hotspot mutations (P286R, S297F, V411L, and A456P) and were of endometrioid histotype (P = 5e-04). These six tumors had an early onset (average age of patients at onset, 48.33 years) when compared to non-POLE endometrioid OC cohort (average age at onset, 60.13 years; P = .008). Samples from TCGA and GTEx had a low COSMIC signature 10 contribution (median 0.16 mutations per megabase; maximum 1.78 mutations per megabase) and carried no POLE hotspot mutations.

CONCLUSIONS

From the largest cohort of RNA-seq from endometrioid OC to date (n = 53), we identified six hypermutated samples likely driven by POLE (frequency, 11%). Our result suggests the clinical need to screen for POLE driver mutations in endometrioid OC, which can guide enrollment in immunotherapy clinical trials.

Gene-based mediation analysis in epigenetic studies

Mediation analysis has been a useful tool for investigating the effect of mediators that lie in the path from the independent variable to the outcome. With the increasing dimensionality of mediators such as in (epi)genomics studies, high-dimensional mediation model is needed. In this work, we focus on epigenetic studies with the goal to identify important DNA methylations that act as mediators between an exposure disease outcome. Specifically, we focus on gene-based high-dimensional mediation analysis implemented with kernel principal component analysis to capture potential nonlinear mediation effect. We first review the current high-dimensional mediation models and then propose two gene-based analytical approaches: gene-based high-dimensional mediation analysis based on linearity assumption between mediators and outcome (gHMA-L) and gene-based high-dimensional mediation analysis based on nonlinearity assumption (gHMA-NL). Since the underlying true mediation relationship is unknown in practice, we further propose an omnibus test of gene-based high-dimensional mediation analysis (gHMA-O) by combing gHMA-L and gHMA-NL. Extensive simulation studies show that gHMA-L performs better under the model linear assumption and gHMA-NL does better under the model nonlinear assumption, while gHMA-O is a more powerful and robust method by combining the two. We apply the proposed methods to two datasets to investigate genes whose methylation levels act as important mediators in the relationship: (1) between alcohol consumption and epithelial ovarian cancer risk using data from the Mayo Clinic Ovarian Cancer Case-Control Study and (2) between childhood maltreatment and comorbid post-traumatic stress disorder and depression in adulthood using data from the Gray Trauma Project.

BRCA1 Promoter Methylation and Clinical Outcomes in Ovarian Cancer: An Individual Patient Data Meta-Analysis

BACKGROUND

BRCA1 methylation has been associated with homologous recombination deficiency, a biomarker of platinum sensitivity. Studies evaluating BRCA1-methylated tubal and ovarian cancer (OC) do not consistently support improved survival following platinum chemotherapy. We examine the characteristics of BRCA1-methylated OC in a meta-analysis of individual participant data.

METHODS

Data of 2636 participants across 15 studies were analyzed. BRCA1-methylated tumors were defined according to their original study. Associations between BRCA1 methylation and clinicopathological characteristics were evaluated. The effects of methylation on overall survival (OS) and progression-free survival (PFS) were examined using mixed-effects models. All statistical tests were 2-sided.

RESULTS

430 (16.3%) tumors were BRCA1-methylated. BRCA1 methylation was associated with younger age and advanced-stage, high-grade serous OC. There were no survival differences between BRCA1-methylated and non-BRCA1-methylated OC (median PFS = 20.0 vs 18.5 months, hazard ratio [HR] = 1.01, 95% CI = 0.87 to 1.16; P = .98; median OS = 46.6 vs 48.0 months, HR = 1.02, 95% CI = 0.87 to 1.18; P = .96). Where BRCA1/2 mutations were evaluated (n = 1248), BRCA1 methylation displayed no survival advantage over BRCA1/2-intact (BRCA1/2 wild-type non-BRCA1-methylated) OC. Studies used different methods to define BRCA1 methylation. Where BRCA1 methylation was determined using methylation-specific polymerase chain reaction and gel electrophoresis (n = 834), it was associated with improved survival (PFS: HR = 0.80, 95% CI = 0.66 to 0.97; P = .02; OS: HR = 0.80, 95% CI = 0.63 to 1.00; P = .05) on mixed-effects modeling.

CONCLUSION

BRCA1-methylated OC displays similar clinicopathological features to BRCA1-mutated OC but is not associated with survival. Heterogeneity within BRCA1 methylation assays influences associations. Refining these assays may better identify cases with silenced BRCA1 function and improved patient outcomes.

Expanding Our Understanding of Ovarian Cancer Risk: The Role of Incomplete Pregnancies

BACKGROUND

Parity is associated with decreased risk of invasive ovarian cancer; however, the relationship between incomplete pregnancies and invasive ovarian cancer risk is unclear. This relationship was examined using 15 case-control studies from the Ovarian Cancer Association Consortium (OCAC). Histotype-specific associations, which have not been examined previously with large sample sizes, were also evaluated.

METHODS

A pooled analysis of 10 470 invasive epithelial ovarian cancer cases and 16 942 controls was conducted. Odds ratios (ORs) and 95% confidence intervals (CIs) for the association between incomplete pregnancies and invasive epithelial ovarian cancer were estimated using logistic regression. All models were conditioned on OCAC study, race and ethnicity, age, and education level and adjusted for number of complete pregnancies, oral contraceptive use, and history of breastfeeding. The same approach was used for histotype-specific analyses.

RESULTS

Ever having an incomplete pregnancy was associated with a 16% reduction in ovarian cancer risk (OR = 0.84, 95% CI = 0.79 to 0.89). There was a trend of decreasing risk with increasing number of incomplete pregnancies (2-sided Ptrend < .001). An inverse association was observed for all major histotypes; it was strongest for clear cell ovarian cancer.

CONCLUSIONS

Incomplete pregnancies are associated with a reduced risk of invasive epithelial ovarian cancer. Pregnancy, including incomplete pregnancy, was associated with a greater reduction in risk of clear cell ovarian cancer, but the result was broadly consistent across histotypes. Future work should focus on understanding the mechanisms underlying this reduced risk.

Mendelian randomization analyses suggest a role for cholesterol in the development of endometrial cancer

Blood lipids have been associated with the development of a range of cancers, including breast, lung and colorectal cancer. For endometrial cancer, observational studies have reported inconsistent associations between blood lipids and cancer risk. To reduce biases from unmeasured confounding, we performed a bidirectional, two-sample Mendelian randomization analysis to investigate the relationship between levels of three blood lipids (low-density lipoprotein [LDL] and high-density lipoprotein [HDL] cholesterol, and triglycerides) and endometrial cancer risk. Genetic variants associated with each of these blood lipid levels (P < 5 × 10-8 ) were identified as instrumental variables, and assessed using genome-wide association study data from the Endometrial Cancer Association Consortium (12 906 cases and 108 979 controls) and the Global Lipids Genetic Consortium (n = 188 578). Mendelian randomization analyses found genetically raised LDL cholesterol levels to be associated with lower risks of endometrial cancer of all histologies combined, and of endometrioid and non-endometrioid subtypes. Conversely, higher genetically predicted HDL cholesterol levels were associated with increased risk of non-endometrioid endometrial cancer. After accounting for the potential confounding role of obesity (as measured by genetic variants associated with body mass index), the association between genetically predicted increased LDL cholesterol levels and lower endometrial cancer risk remained significant, especially for non-endometrioid endometrial cancer. There was no evidence to support a role for triglycerides in endometrial cancer development. Our study supports a role for LDL and HDL cholesterol in the development of non-endometrioid endometrial cancer. Further studies are required to understand the mechanisms underlying these findings.

Generating real-world evidence from unstructured clinical notes to examine clinical utility of genetic tests: use case in BRCAness

Background

Next-generation sequencing provides comprehensive information about individuals’ genetic makeup and is commonplace in oncology clinical practice. However, the utility of genetic information in the clinical decision-making process has not been examined extensively from a real-world, data-driven perspective. Through mining real-world data (RWD) from clinical notes, we could extract patients’ genetic information and further associate treatment decisions with genetic information.

Methods

We proposed a real-world evidence (RWE) study framework that incorporates context-based natural language processing (NLP) methods and data quality examination before final association analysis. The framework was demonstrated in a Foundation-tested women cancer cohort (N = 196). Upon retrieval of patients’ genetic information using NLP system, we assessed the completeness of genetic data captured in unstructured clinical notes according to a genetic data-model. We examined the distribution of different topics regarding BRCA1/2 throughout patients’ treatment process, and then analyzed the association between BRCA1/2 mutation status and the discussion/prescription of targeted therapy.

Results

We identified seven topics in the clinical context of genetic mentions including: Information, Evaluation, Insurance, Order, Negative, Positive, and Variants of unknown significance. Our rule-based system achieved a precision of 0.87, recall of 0.93 and F-measure of 0.91. Our machine learning system achieved a precision of 0.901, recall of 0.899 and F-measure of 0.9 for four-topic classification and a precision of 0.833, recall of 0.823 and F-measure of 0.82 for seven-topic classification. We found in result-containing sentences, the capture of BRCA1/2 mutation information was 75%, but detailed variant information (e.g. variant types) is largely missing. Using cleaned RWD, significant associations were found between BRCA1/2 positive mutation and targeted therapies.

Conclusions

In conclusion, we demonstrated a framework to generate RWE using RWD from different clinical sources. Rule-based NLP system achieved the best performance for resolving contextual variability when extracting RWD from unstructured clinical notes. Data quality issues such as incompleteness and discrepancies exist thus manual data cleaning is needed before further analysis can be performed. Finally, we were able to use cleaned RWD to evaluate the real-world utility of genetic information to initiate a prescription of targeted therapy.

Refined cut-off for TP53 immunohistochemistry improves prediction of TP53 mutation status in ovarian mucinous tumors: implications for outcome analyses

TP53 mutations are implicated in the progression of mucinous borderline tumors (MBOT) to mucinous ovarian carcinomas (MOC). Optimized immunohistochemistry (IHC) for TP53 has been established as a proxy for the TP53 mutation status in other ovarian tumor types. We aimed to confirm the ability of TP53 IHC to predict TP53 mutation status in ovarian mucinous tumors and to evaluate the association of TP53 mutation status with survival among patients with MBOT and MOC. Tumor tissue from an initial cohort of 113 women with MBOT/MOC was stained with optimized IHC for TP53 using tissue microarrays (75.2%) or full sections (24.8%) and interpreted using established criteria as normal or abnormal (overexpression, complete absence, or cytoplasmic). Cases were considered concordant if abnormal IHC staining predicted deleterious TP53 mutations. Discordant tissue microarray cases were re-evaluated on full sections and interpretational criteria were refined. The initial cohort was expanded to a total of 165 MBOT and 424 MOC for the examination of the association of survival with TP53 mutation status, assessed either by TP53 IHC and/or sequencing. Initially, 82/113 (72.6%) cases were concordant using the established criteria. Refined criteria for overexpression to account for intratumoral heterogeneity and terminal differentiation improved concordance to 93.8% (106/113). In the expanded cohort, 19.4% (32/165) of MBOT showed evidence for TP53 mutation and this was associated with a higher risk of recurrence, disease-specific death, and all-cause mortality (overall survival: HR = 4.6, 95% CI 1.5-14.3, p = 0.0087). Within MOC, 61.1% (259/424) harbored a TP53 mutation, but this was not associated with survival (overall survival, p = 0.77). TP53 IHC is an accurate proxy for TP53 mutation status with refined interpretation criteria accounting for intratumoral heterogeneity and terminal differentiation in ovarian mucinous tumors. TP53 mutation status is an important biomarker to identify MBOT with a higher risk of mortality.

Development and Validation of the Gene Expression Predictor of High-grade Serous Ovarian Carcinoma Molecular SubTYPE (PrOTYPE)

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 >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.

Circulating CD14+ HLA-DRlo/- monocytic cells as a biomarker for epithelial ovarian cancer progression

PROBLEM

Previous studies identified circulating CD14⁺ HLA-DR^lo/- monocytic cells as an immune suppressive subset in solid malignancies, such as prostate, renal cell carcinoma, and pancreatic cancer. Such monocytic cells have been implicated not only in tumour progression but also as a potential barrier for immunotherapy. This study examined the relationship between the frequency of circulating monocytic cells and epithelial ovarian cancer (EOC) progression pre- and post-frontline chemotherapy, defined by disease stage, which is a leading prognostic factor for this malignancy.

METHOD OF STUDY

Incident cases of 236 women with EOC were recruited and comprehensive flow cytometry was utilized to assess the frequency of peripheral blood CD33⁺ CD11b⁺ HLA-DR^-/low CD14⁺ CD15^- monocytic cells, henceforth termed CD14⁺ HLA-DR^lo/- monocytic cells, prior to and after completion of frontline chemotherapy. Multivariable odds ratios (OR) were used to estimate the association between CD14⁺ HLA-DR^lo/- monocytic cell percentages and disease stage. Wilcoxon signed-rank tests evaluated changes in these monocytic cell levels pre- and post-chemotherapy in a patient subset (n = 70).

RESULTS

Patients with elevated frequencies of circulating CD14⁺ HLA-DR^lo/- monocytic cells at diagnosis were at 3.33-fold greater odds of having advanced stage (III/IV) EOC (CI: 1.04-10.64), with a significant trend in increasing CD14⁺ HLA-DR^lo/- monocytic cell levels (P = .04). There was a 2.02% median decrease of these monocytic cells post-chemotherapy among a subset of patients with advanced stage disease (P < .0001).

CONCLUSION

These findings support the potential clinical relevance of CD14⁺ HLA-DR^lo/- monocytic cells in EOC for prognosis and may indicate a non-invasive biomarker to measure disease progression.

Abstract A19: Aspirin use and ovarian cancer risk by lymphocyte infiltration

Introduction: There is growing evidence to support a modest inverse association between daily aspirin use and ovarian cancer risk. We previously hypothesized that this association may differ by immune cell infiltration. In our prior study evaluating heterogeneity by tumor-associated macrophage (TAM) infiltration, aspirin use was associated with lower risk of ovarian cancer with high protumorigenic M2-type TAM infiltration, but not low M2-type TAM infiltration. The goal of the present study was to evaluate the hypothesis that the association between aspirin use and ovarian cancer risk also differs by infiltration with tumor-infiltrating lymphocytes (TILs).

Methods: We included cases and matched controls from the Nurses’ Health Study (NHS), NHSII, and New England Case Control Study. Study participants self-reported aspirin use including regular use (≥1 x per week), duration of use, and number of tablets. A gynecologic pathologist (JLH), who was blinded to exposure status, reviewed hematoxylin and eosin stained tumor slides to quantify TILs for all cases (n=547). We used polytomous logistic regression, adjusted for ovarian cancer risk factors, to estimate odds ratios (OR) for aspirin use and ovarian cancer risk by TIL infiltration.

Results: Overall, recent aspirin use was associated with lower risk of ovarian cancer compared with no regular aspirin use (OR=0.65, 95%CI=0.47-0.89). Recent aspirin use was associated with a lower risk of ovarian cancer exhibiting TIL infiltration (OR=0.38, 95%CI=0.26-0.56) but not with TIL-negative cancers (OR=1.00, 95%CI=0.71-1.41; p-heterogeneity &lt;0.001). We observed similar results for aspirin duration and tablets. For duration, comparing ≥10 years of aspirin use to &lt;1 year we observed a significantly lower risk of ovarian cancer with TILs (OR=0.47, 95%CI=0.41-0.71; p-trend=0.004), but a nonsignificant reduction for ovarian cancer without TILs (OR=0.86, 95%CI=0.59-1.26; p-trend=0.72; p-heterogeneity=0.03). Comparing ≥6 tablets to &lt;1 tablet per week, there was evidence of a lower risk of ovarian cancer with TILs (OR=0.47, 95%CI=0.30-0.75; p-trend=0.001) but not without TILs (OR=1.02, 95%CI=0.70-1.48; p-trend=0.84; p-heterogeneity=0.003).

Conclusions: Our results suggest that prediagnosis aspirin use may be an important modulator of immune cell infiltration in ovarian tumors, supporting an immunogenic mechanism for aspirin in ovarian carcinogenesis. Further work is under way to evaluate nonaspirin nonsteroidal anti-inflammatory drugs, and to consider which subsets of TAMs and TILs may be driving heterogeneity in the associations between anti-inflammatory drugs and ovarian cancer risk.

Citation Format: Mollie E. Barnard, Britton Trabert, Jonathan L. Hecht, Ellen L. Goode, Naoko Sasamoto, Kathryn L. Terry, Shelley S. Tworoger. Aspirin use and ovarian cancer risk by lymphocyte infiltration [abstract]. In: Proceedings of the AACR Special Conference on Modernizing Population Sciences in the Digital Age; 2019 Feb 19-22; San Diego, CA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2020;29(9 Suppl):Abstract nr A19.

Clinical and pathological associations of PTEN expression in ovarian cancer: a multicentre study from the Ovarian Tumour Tissue Analysis Consortium

BACKGROUND

PTEN loss is a putative driver in histotypes of ovarian cancer (high-grade serous (HGSOC), endometrioid (ENOC), clear cell (CCOC), mucinous (MOC), low-grade serous (LGSOC)). We aimed to characterise PTEN expression as a biomarker in epithelial ovarian cancer in a large population-based study.

METHODS

Tumours from 5400 patients from a multicentre observational, prospective cohort study of the Ovarian Tumour Tissue Analysis Consortium were used to evaluate associations between immunohistochemical PTEN patterns and overall survival time, age, stage, grade, residual tumour, CD8+ tumour-infiltrating lymphocytes (TIL) counts, expression of oestrogen receptor (ER), progesterone receptor (PR) and androgen receptor (AR) by means of Cox proportional hazard models and generalised Cochran-Mantel-Haenszel tests.

RESULTS

Downregulation of cytoplasmic PTEN expression was most frequent in ENOC (most frequently in younger patients; p value = 0.0001) and CCOC and was associated with longer overall survival in HGSOC (hazard ratio: 0.78, 95% CI: 0.65-0.94, p value = 0.022). PTEN expression was associated with ER, PR and AR expression (p values: 0.0008, 0.062 and 0.0002, respectively) in HGSOC and with lower CD8 counts in CCOC (p value < 0.0001). Heterogeneous expression of PTEN was more prevalent in advanced HGSOC (p value = 0.019) and associated with higher CD8 counts (p value = 0.0016).

CONCLUSIONS

PTEN loss is a frequent driver in ovarian carcinoma associating distinctly with expression of hormonal receptors and CD8+ TIL counts in HGSOC and CCOC histotypes.

Menopausal hormone therapy prior to the diagnosis of ovarian cancer is associated with improved survival

PURPOSE

Prior studies of menopausal hormone therapy (MHT) and ovarian cancer survival have been limited by lack of hormone regimen detail and insufficient sample sizes. To address these limitations, a comprehensive analysis of 6419 post-menopausal women with pathologically confirmed ovarian carcinoma was conducted to examine the association between MHT use prior to diagnosis and survival.

METHODS

Data from 15 studies in the Ovarian Cancer Association Consortium were included. MHT use was examined by type (estrogen-only (ET) or estrogen+progestin (EPT)), duration, and recency of use relative to diagnosis. Cox proportional hazards models were used to estimate the association between hormone therapy use and survival. Logistic regression and mediation analysis was used to explore the relationship between MHT use and residual disease following debulking surgery.

RESULTS

Use of ET or EPT for at least five years prior to diagnosis was associated with better ovarian cancer survival (hazard ratio, 0.80; 95% CI, 0.74 to 0.87). Among women with advanced stage, high-grade serous carcinoma, those who used MHT were less likely to have any macroscopic residual disease at the time of primary debulking surgery (p for trend <0.01 for duration of MHT use). Residual disease mediated some (17%) of the relationship between MHT and survival.

CONCLUSIONS

Pre-diagnosis MHT use for 5+ years was a favorable prognostic factor for women with ovarian cancer. This large study is consistent with prior smaller studies, and further work is needed to understand the underlying mechanism.

Prognostic gene expression signature for high-grade serous ovarian cancer

BACKGROUND

Median overall survival (OS) for women with high-grade serous ovarian cancer (HGSOC) is ∼4 years, yet survival varies widely between patients. There are no well-established, gene expression signatures associated with prognosis. The aim of this study was to develop a robust prognostic signature for OS in patients with HGSOC.

PATIENTS AND METHODS

Expression of 513 genes, selected from a meta-analysis of 1455 tumours and other candidates, was measured using NanoString technology from formalin-fixed paraffin-embedded tumour tissue collected from 3769 women with HGSOC from multiple studies. Elastic net regularization for survival analysis was applied to develop a prognostic model for 5-year OS, trained on 2702 tumours from 15 studies and evaluated on an independent set of 1067 tumours from six studies.

RESULTS

Expression levels of 276 genes were associated with OS (false discovery rate < 0.05) in covariate-adjusted single-gene analyses. The top five genes were TAP1, ZFHX4, CXCL9, FBN1 and PTGER3 (P < 0.001). The best performing prognostic signature included 101 genes enriched in pathways with treatment implications. Each gain of one standard deviation in the gene expression score conferred a greater than twofold increase in risk of death [hazard ratio (HR) 2.35, 95% confidence interval (CI) 2.02-2.71; P < 0.001]. Median survival [HR (95% CI)] by gene expression score quintile was 9.5 (8.3 to -), 5.4 (4.6-7.0), 3.8 (3.3-4.6), 3.2 (2.9-3.7) and 2.3 (2.1-2.6) years.

CONCLUSION

The OTTA-SPOT (Ovarian Tumor Tissue Analysis consortium - Stratified Prognosis of Ovarian Tumours) gene expression signature may improve risk stratification in clinical trials by identifying patients who are least likely to achieve 5-year survival. The identified novel genes associated with the outcome may also yield opportunities for the development of targeted therapeutic approaches.

Abstract 2452: Evidence of decreased X chromosome inactivation in primary ovarian tumors

Copy number alterations and loss of heterozygosity (LOH) are common in ovarian tumors, including those on the X chromosome. In females, one of the two copies of the X chromosome is transcriptionally silenced during embryonic development in the epigenetic process of X chromosome inactivation. Specifically, the long non-coding RNA XIST is expressed from the inactive X (Xi) and triggers chromosome-wide DNA methylation. In breast and ovarian cancer cell line studies, Xi deletion with concomitant duplication of the active X (Xa) has been observed. We sought to examine which copy of the X chromosome, Xi or Xa, was gained or lost in a set of primary ovarian tumors, using levels of XIST expression as a marker of X chromosome inactivation. To determine which copy of the X chromosome harbors the alterations, we compared local copy number and LOH in Xq13.2 to the levels of XIST expression in a set of N=134 primary tumors from Mayo Clinic ovarian cancer patients. Associations with XIST expression were examined using ANOVA, with copy number of two and absence of LOH as the reference level. Associations between expression of XIST and expression of other genes were examined with a mixed effects model, to account for the repeated observations across patients, after adjusting for local copy number and LOH. Associations with XIST expression and overall survival were performed with Cox regression. Additional analyses were performed in a subset of 24 patient tumors with X chromosome-wide alterations to simplify interpretation, including a visual examination of chromosome-wide DNA methylation patterns. Analyses were replicated in N=106 ovarian tumors from TCGA with X chromosome-wide alterations. Across N=134 ovarian tumors from Mayo Clinic, abnormal copy number (not equal to two) was associated with reduced XIST expression (p=0.0016). Levels of XIST expression were similar for tumors with one copy and tumors with three or more copies. Similarly, levels of XIST expression were also lower for tumors with LOH (p=5.2E-5). Reduced XIST expression was also associated with poorer survival, after adjusting for age (p=0.02). Results were consistent in the subset of tumors with whole chromosome alterations (pCN=0.0013, pLOH=0.0019). For these tumors, when LOH was present, XIST expression was low and the DNA methylation distribution is bi-modal, suggesting that the inactive X was lost and the active X was gained. Furthermore, after adjusting for copy number and LOH, XIST expression was weakly positively correlated with the median expression of other genes (R2=0.04, p=0.026). In N=106 TCGA tumors with whole X chromosome alterations, XIST expression was also reduced in tumors with abnormal copy number (p=0.006) and with LOH (p=1.5E-10). Taken together, because XIST is expressed from Xi and not Xa, associations of reduced XIST expression with abnormal copy number and presence of LOH suggest that Xi is frequently lost and Xa is gained in ovarian tumors, consistent with cell line studies. These results should be further validated in larger sample sets and examined for association with clinical outcome.

Citation Format: Yilin Song, Bryan M. McCauley, Melissa C. Larson, Sebastian M. Armasu, Kate Lawrenson, Ellen L. Goode, Nicholas B. Larson, Stacey J. Winham. Evidence of decreased X chromosome inactivation in primary ovarian tumors [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 2452.