Modugno, Francesmary

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

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

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

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

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

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

Abstract 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

Abstract Background: The link between modifiable breast cancer risk factors and tumor genomic alterations remains largely unexplored. We evaluated the association of prediagnostic body mass index (BMI), cigarette smoking, and alcohol consumption with somatic copy number variation (SCNV), total somatic mutation burden (TSMB), seven single base substitution (SBS) signatures (SBS1, SBS2, SBS3, SBS5, SBS13, SBS29, and SBS30), and nine driver mutations (CDH1, GATA3, KMT2C, MAP2K4, MAP3K1, NCOR1, PIK3CA, RUNX1, and TP53) in a subset of The Cancer Genome Atlas (TCGA). Methods: Clinical and genomic data were retrieved from the TCGA database. Risk factor information was collected from four TCGA sites (n = 219 women), including BMI (1 year before diagnosis), cigarette smoking (smokers/nonsmokers), and alcohol consumption (current drinkers/nondrinkers). Multivariable regression analyses were conducted in all tumors and stratified according to estrogen receptor (ER) status. Results: Increasing BMI was associated with increasing SCNV in all women (P = 0.039) and among women with ER− tumors (P = 0.031). Smokers had higher SCNV and TSMB versus nonsmokers (P < 0.05 all women). Alcohol drinkers had higher SCNV versus nondrinkers (P < 0.05 all women and among women with ER+ tumors). SBS3 (defective homologous recombination-based repair) was exclusively found in alcohol drinkers with ER− disease. GATA3 mutation was more likely to occur in women with higher BMI. No association was significant after multiple testing correction. Conclusions: This study provides preliminary evidence that BMI, cigarette smoking, and alcohol consumption can influence breast tumor biology, in particular, DNA alterations. Impact: This study demonstrates a link between modifiable breast cancer risk factors and tumor genomic alterations.

As a follow‐up to genome‐wide association analysis of common variants associated with ovarian carcinoma (cancer), our study considers seven well‐known ovarian cancer risk factors and their interactions with 28 genome‐wide significant common genetic variants. The interaction analyses were based on data from 9971 ovarian cancer cases and 15,566 controls from 17 case–control studies. Likelihood ratio and Wald tests for multiplicative interaction and for relative excess risk due to additive interaction were used. The top multiplicative interaction was noted between oral contraceptive pill (OCP) use (ever vs. never) and rs13255292 (p value = 3.48 × 10−4). Among women with the TT genotype for this variant, the odds ratio for OCP use was 0.53 (95% CI = 0.46–0.60) compared to 0.71 (95%CI = 0.66–0.77) for women with the CC genotype. When stratified by duration of OCP use, women with 1–5 years of OCP use exhibited differential protective benefit across genotypes. However, no interaction on either the multiplicative or additive scale was found to be statistically significant after multiple testing correction. The results suggest that OCP use may offer increased benefit for women who are carriers of the T allele in rs13255292. On the other hand, for women carrying the C allele in this variant, longer (5+ years) use of OCP may reduce the impact of carrying the risk allele of this SNP. Replication of this finding is needed. The study presents a comprehensive analytic framework for conducting gene–environment analysis in ovarian cancer.