Goode, Ellen L.

Abstract Purpose: Here we assess whether molecular subtyping identifies biological features of tumors that correlate with survival and surgical outcomes of high-grade serous ovarian cancer (HGSOC). Experimental Design: Consensus clustering of pooled mRNA expression data from over 2,000 HGSOC cases was used to define molecular subtypes of HGSOCs. This de novo classification scheme was then applied to 381 Mayo Clinic HGSOC patients with detailed survival and surgical outcome information. Results: Five molecular subtypes of HGSOC were identified. In the pooled dataset, three subtypes were largely concordant with prior studies describing proliferative, mesenchymal, and immunoreactive tumors (concordance > 70%), and the group of tumors previously described as differentiated type was segregated into two new types, one of which (anti-mesenchymal) had downregulation of genes that were typically upregulated in the mesenchymal subtype. Molecular subtypes were significantly associated with overall survival (P < 0.001) and with rate of optimal surgical debulking (≤1 cm, P = 1.9E−4) in the pooled dataset. Among stage III-C or IV Mayo Clinic patients, molecular subtypes were also significantly associated with overall survival (P = 0.001), as well as rate of complete surgical debulking (no residual disease; 16% in mesenchymal tumors compared with >28% in other subtypes; P = 0.02). Conclusions: HGSOC tumors may be categorized into five molecular subtypes that associate with overall survival and the extent of residual disease following debulking surgery. Because mesenchymal tumors may have features that were associated with less favorable surgical outcome, molecular subtyping may have future utility in guiding neoadjuvant treatment decisions for women with HGSOC. Clin Cancer Res; 23(15); 4077–85. ©2017 AACR.

Abstract Purpose: Recent progress in understanding the molecular biology of epithelial ovarian cancer has not yet translated into individualized treatment for these women or improvements in their disease outcome. Gene expression has been utilized to identify distinct molecular subtypes, but there have been no reports investigating whether or not molecular subtyping is predictive of response to bevacizumab in ovarian cancer. Experimental Design: DASL gene expression arrays were performed on FFPE tissue from patients enrolled on the ICON7 trial. Patients were stratified into four TCGA molecular subtypes. Associations between molecular subtype and the efficacy of randomly assigned therapy with bevacizumab were assessed. Results: Molecular subtypes were assigned as follows: 122 immunoreactive (34%), 96 proliferative (27%), 73 differentiated (20%), and 68 mesenchymal (19%). In univariate analysis patients with tumors of proliferative subtype obtained the greatest benefit from bevacizumab with a median PFS improvement of 10.1 months [HR, 0.55 (95% CI, 0.34–0.90), P = 0.016]. For the mesenchymal subtype, bevacizumab conferred a nonsignificant improvement in PFS of 8.2 months [HR 0.78 (95% CI, 0.44–1.40), P = 0.41]. Bevacizumab conferred modest improvements in PFS for patients with immunoreactive subtype (3.8 months; P = 0.08) or differentiated subtype (3.7 months; P = 0.61). Multivariate analysis demonstrated significant PFS improvement in proliferative subtype patients only [HR, 0.45 (95% CI, 0.27–0.74), P = 0.0015]. Conclusions: Ovarian carcinoma molecular subtypes with the poorest survival (proliferative and mesenchymal) derive a comparably greater benefit from treatment that includes bevacizumab. Validation of our findings in an independent cohort could enable the use of bevacizumab for those patients most likely to benefit, thereby reducing side effects and healthcare cost. Clin Cancer Res; 23(14); 3794–801. ©2017 AACR.

Abstract Background: Epigenetic disturbances are crucial in cancer initiation, potentially with pleiotropic effects, and may be influenced by the genetic background. Methods: In a subsets (ASSET) meta-analytic approach, we investigated associations of genetic variants related to epigenetic mechanisms with risks of breast, lung, colorectal, ovarian and prostate carcinomas using 51,724 cases and 52,001 controls. False discovery rate–corrected P values (q values < 0.05) were considered statistically significant. Results: Among 162,887 imputed or genotyped variants in 555 candidate genes, SNPs in eight genes were associated with risk of more than one cancer type. For example, variants in BABAM1 were confirmed as a susceptibility locus for squamous cell lung, overall breast, estrogen receptor (ER)–negative breast, and overall prostate, and overall serous ovarian cancer; the most significant variant was rs4808076 [OR = 1.14; 95% confidence interval (CI) = 1.10–1.19; q = 6.87 × 10−5]. DPF1 rs12611084 was inversely associated with ER-negative breast, endometrioid ovarian, and overall and aggressive prostate cancer risk (OR = 0.93; 95% CI = 0.91–0.96; q = 0.005). Variants in L3MBTL3 were associated with colorectal, overall breast, ER-negative breast, clear cell ovarian, and overall and aggressive prostate cancer risk (e.g., rs9388766: OR = 1.06; 95% CI = 1.03–1.08; q = 0.02). Variants in TET2 were significantly associated with overall breast, overall prostate, overall ovarian, and endometrioid ovarian cancer risk, with rs62331150 showing bidirectional effects. Analyses of subpathways did not reveal gene subsets that contributed disproportionately to susceptibility. Conclusions: Functional and correlative studies are now needed to elucidate the potential links between germline genotype, epigenetic function, and cancer etiology. Impact: This approach provides novel insight into possible pleiotropic effects of genes involved in epigenetic processes. Cancer Epidemiol Biomarkers Prev; 26(6); 816–25. ©2017 AACR.

Abstract Background: Genome-wide association studies (GWAS) have identified multiple loci associated with epithelial ovarian cancer (EOC) susceptibility, but further progress requires integration of epidemiology and biology to illuminate true risk loci below genome-wide significance levels (P < 5 × 10−8). Most risk SNPs lie within non–protein-encoding regions, and we hypothesize that long noncoding RNA (lncRNA) genes are enriched at EOC risk regions and represent biologically relevant functional targets. Methods: Using imputed GWAS data from about 18,000 invasive EOC cases and 34,000 controls of European ancestry, the GENCODE (v19) lncRNA database was used to annotate SNPs from 13,442 lncRNAs for permutation-based enrichment analysis. Tumor expression quantitative trait locus (eQTL) analysis was performed for sub-genome-wide regions (1 × 10−5 > P > 5 × 10−8) overlapping lncRNAs. Results: Of 5,294 EOC-associated SNPs (P < 1.0 × 10−5), 1,464 (28%) mapped within 53 unique lncRNAs and an additional 3,484 (66%) SNPs were correlated (r2 > 0.2) with SNPs within 115 lncRNAs. EOC-associated SNPs comprised 130 independent regions, of which 72 (55%) overlapped with lncRNAs, representing a significant enrichment (P = 5.0 × 10−4) that was more pronounced among a subset of 5,401 lncRNAs with active epigenetic regulation in normal ovarian tissue. EOC-associated lncRNAs and their putative promoters and transcription factors were enriched for biologically relevant pathways and eQTL analysis identified five novel putative risk regions with allele-specific effects on lncRNA gene expression. Conclusions: lncRNAs are significantly enriched at EOC risk regions, suggesting a mechanistic role for lncRNAs in driving predisposition to EOC. Impact: lncRNAs represent key candidates for integrative epidemiologic and functional studies. Further research on their biologic role in ovarian cancer is indicated. Cancer Epidemiol Biomarkers Prev; 26(1); 116–25. ©2016 AACR.

Abstract Four gene expression subtypes of high-grade serous ovarian cancer (HGSC) have been previously described. In these early studies, a fraction of samples that did not fit well into the four subtype classifications were excluded. Therefore, we sought to systematically determine the concordance of transcriptomic HGSC subtypes across populations without removing any samples. We created a bioinformatics pipeline to independently cluster the five largest mRNA expression datasets using k-means and nonnegative matrix factorization (NMF). We summarized differential expression patterns to compare clusters across studies. While previous studies reported four subtypes, our cross-population comparison does not support four. Because these results contrast with previous reports, we attempted to reproduce analyses performed in those studies. Our results suggest that early results favoring four subtypes may have been driven by the inclusion of serous borderline tumors. In summary, our analysis suggests that either two or three, but not four, gene expression subtypes are most consistent across datasets.

Abstract Background: The precise mechanism by which the immune system is adversely affected in cancer patients remains poorly understood, but the accumulation of immunosuppressive/protumorigenic myeloid-derived suppressor cells (MDSCs) is thought to be a prominent mechanism contributing to immunologic tolerance of malignant cells in epithelial ovarian cancer (EOC). To this end, we hypothesized genetic variation in MDSC pathway genes would be associated with survival after EOC diagnoses. Methods: We measured the hazard of death due to EOC within 10 years of diagnosis, overall and by invasive subtype, attributable to SNPs in 24 genes relevant in the MDSC pathway in 10,751 women diagnosed with invasive EOC. Versatile Gene-based Association Study and the admixture likelihood method were used to test gene and pathway associations with survival. Results: We did not identify individual SNPs that were significantly associated with survival after correction for multiple testing (P < 3.5 × 10−5), nor did we identify significant associations between the MDSC pathway overall, or the 24 individual genes and EOC survival. Conclusions: In this well-powered analysis, we observed no evidence that inherited variations in MDSC-associated SNPs, individual genes, or the collective genetic pathway contributed to EOC survival outcomes. Impact: Common inherited variation in genes relevant to MDSCs was not associated with survival in women diagnosed with invasive EOC. Cancer Epidemiol Biomarkers Prev; 26(3); 420–4. ©2016 AACR.

MicroRNAs (miRNAs) regulate up to one‐third of all protein‐coding genes including genes relevant to cancer. Variants within miRNAs have been reported to be associated with prognosis, survival, response to chemotherapy across cancer types, in vitro parameters of cell growth, and altered risks for development of cancer. Five miRNA variants have been reported to be associated with risk for development of colorectal cancer (CRC). In this study, we evaluated germline genetic variation in 1,123 miRNAs in 899 individuals with CRCs categorized by clinical subtypes and in 204 controls. The role of common miRNA variation in CRC was investigated using single variant and miRNA‐level association tests. Twenty‐nine miRNAs and 30 variants exhibited some marginal association with CRC in at least one subtype of CRC. Previously reported associations were not confirmed (n = 4) or could not be evaluated (n = 1). The variants noted for the CRCs with deficient mismatch repair showed little overlap with the variants noted for CRCs with proficient mismatch repair, consistent with our evolving understanding of the distinct biology underlying these two groups. © 2016 The Authors Genes, Chromosomes & Cancer Published by Wiley Periodicals, Inc.