Fridley, Brooke L.

Abstract Background: While numerous susceptibility loci for epithelial ovarian cancer (EOC) have been identified, few associations have been reported with overall survival. In the absence of common prognostic genetic markers, we hypothesize that rare coding variants may be associated with overall EOC survival and assessed their contribution in two exome-based genotyping projects of the Ovarian Cancer Association Consortium (OCAC). Methods: The primary patient set (Set 1) included 14 independent EOC studies (4,293 patients) and 227,892 variants, and a secondary patient set (Set 2) included six additional EOC studies (1,744 patients) and 114,620 variants. Because power to detect rare variants individually is reduced, gene-level tests were conducted. Sets were analyzed separately at individual variants and by gene, and then combined with meta-analyses (73,203 variants and 13,163 genes overlapped). Results: No individual variant reached genome-wide statistical significance. A SNP previously implicated to be associated with EOC risk and, to a lesser extent, survival, rs8170, showed the strongest evidence of association with survival and similar effect size estimates across sets (Pmeta = 1.1E−6, HRSet1 = 1.17, HRSet2 = 1.14). Rare variants in ATG2B, an autophagy gene important for apoptosis, were significantly associated with survival after multiple testing correction (Pmeta = 1.1E−6; Pcorrected = 0.01). Conclusions: Common variant rs8170 and rare variants in ATG2B may be associated with EOC overall survival, although further study is needed. Impact: This study represents the first exome-wide association study of EOC survival to include rare variant analyses, and suggests that complementary single variant and gene-level analyses in large studies are needed to identify rare variants that warrant follow-up study. Cancer Epidemiol Biomarkers Prev; 25(3); 446–54. ©2016 AACR.

ABSTRACTEpithelial‐mesenchymal transition (EMT) is a process whereby epithelial cells assume mesenchymal characteristics to facilitate cancer metastasis. However, EMT also contributes to the initiation and development of primary tumors. Prior studies that explored the hypothesis that EMT gene variants contribute to epithelial ovarian carcinoma (EOC) risk have been based on small sample sizes and none have sought replication in an independent population. We screened 15,816 single‐nucleotide polymorphisms (SNPs) in 296 genes in a discovery phase using data from a genome‐wide association study of EOC among women of European ancestry (1,947 cases and 2,009 controls) and identified 793 variants in 278 EMT‐related genes that were nominally (P < 0.05) associated with invasive EOC. These SNPs were then genotyped in a larger study of 14,525 invasive‐cancer patients and 23,447 controls. A P‐value <0.05 and a false discovery rate (FDR) <0.2 were considered statistically significant. In the larger dataset, GPC6/GPC5 rs17702471 was associated with the endometrioid subtype among Caucasians (odds ratio (OR) = 1.16, 95% CI = 1.07–1.25, P = 0.0003, FDR = 0.19), whereas F8 rs7053448 (OR = 1.69, 95% CI = 1.27–2.24, P = 0.0003, FDR = 0.12), F8 rs7058826 (OR = 1.69, 95% CI = 1.27–2.24, P = 0.0003, FDR = 0.12), and CAPN13 rs1983383 (OR = 0.79, 95% CI = 0.69–0.90, P = 0.0005, FDR = 0.12) were associated with combined invasive EOC among Asians. In silico functional analyses revealed that GPC6/GPC5 rs17702471 coincided with DNA regulatory elements. These results suggest that EMT gene variants do not appear to play a significant role in the susceptibility to EOC.

Abstract Epithelial ovarian cancer (EOC) is the fifth leading cause of cancer death among women in the United States (5% of cancer deaths). The standard treatment for patients with advanced EOC is initial debulking surgery followed by carboplatin-paclitaxel combination chemotherapy. Unfortunately, even with modern chemotherapy, most patients relapse and die with the five-year overall survival around 45%. In addition, those patients that initially respond to taxane-platinum therapy eventually develop platinum-resistant tumors and relapse. Thus, finding novel therapeutics for treating EOC is essential. One approach that has been used widely in cancer drug discovery is Connectivity Mapping (CMAP) using gene expression data and drug phenotype data assessed on cancer cell lines (CCL) (i.e. “sensitive” or “resistant” based on dose-response drug data). However, recent publications have highlighted some issues with the use of these standard CCL sets that are often used in these drug screening studies. Therefore, instead of CMAP analysis based on gene expression signatures developed on drug response data collected on a small set of EOC cell lines, we took a novel CMAP approach based on information collected on EOC patients and clinical endpoints. That is, we determined the gene expression signature (e.g., set of genes) that were associated with time to recurrence in data collected on EOC patients within TCGA (n = 518) and a Mayo Clinic study (n = 474). The set of genes associated with time to recurrence (with and without adjustment for additional clinical covariates) was completed separately for TCGA and Mayo Clinic studies, restricting all analyses to EOC tumor of high-grade serous histology. Each of these 4 sets of genes, where genes with hazard ratios (HR) &gt; 1 coded as “positively” associated and genes with HR &lt; 1 coded as “negatively” associated with the phenotype, were inputted into CMAP software (Broad Institute) to determine a set of drugs for which our signature “matches” the “reference” signature. The drugs that overlapped between the CMAP analyses on the two studies were carried forward for validation studies involving drug screens on a set of 10 EOC cell lines. Of the 11 drugs carried forward, we a priori hypothesized that 5 of these drugs would kill EOC cells (mitoxantrone, podophyllotoxin, wortmannin, doxorubicin, and 17-AAG), for which all 5 showed an ability to kill EOC cells in vitro (cell lines were treated using serial dilutions of the drugs (20-0 μM) for 72 h followed by measuring cell viability using the CellTiter-Blue assay). Further research is on-going to understand the functional mechanism by which these candidate drugs might treat EOC and to apply this bioinformatics approach to other cancer drug screening studies. Note: This abstract was not presented at the meeting. Citation Format: Brooke L. Fridley, Rama Raghavan, Gottfried E. Konecny, Chen Wang, Ellen L. Goode, Harsh B. Pathak, Stephen Hyter. A bioinformatics approach to drug discovery: patient based connectivity mapping. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4853. doi:10.1158/1538-7445.AM2015-4853

Abstract Background: We investigated the role of multifactor interactions and epithelial ovarian cancer (EOC) risk using data collected within the Collaborative Oncological Gene-environment Study and the Ovarian Cancer Association Consortium containing 18,000 EOC cases and 26,000 controls. To date, researchers have identified 18 EOC susceptibility loci. However, it has been estimated that many more common variant loci exist; we hypothesized that some of the unexplained variation is due to gene-environment interactions. Similar to breast and endometrial cancers, many EOC risk factors relate to estrogen exposure, and increased levels of estrogen has been associated with obesity in post-menopausal women. Therefore, we hypothesized that gene-environment interactions dealing with hormone-related risk factors could differ between obese and non-obese women. Methods: We considered multifactor interactions between single nucleotide polymorphisms (SNPs) and six hormone-related factors: oral contraceptive use; parity; endometriosis; tubal ligation; hormone replacement therapy (HRT); and estrogen use; and assessed whether these GE interactions differed between obese and non-obese women. The SNPs included in this analysis included the top associated markers within the 18 confirmed risk loci, as well as, 36,811 SNPs that lie within 84 candidate gene regions related to hormone biosynthesis & metabolism and insulin-like growth factors. For the candidate SNP analyses, we assessed multifactor interactions across environmental factors using logistic regression models including age, study site, and population substructure. For the candidate gene analyses, a two-stage screening and testing procedure was utilized. All analyses were completed including: (1) all histologies; and (2) only serous histology. Results: For the 18 known EOC risk loci, no hormone-related risk factors showed statistically significant interaction at a threshold of 10-4. In contrast, the two-stage analysis of the candidate genes found several significant gene-environment-environment effects, after adjusting for multiple testing. For analysis of all cases, the SNPs in the following candidate genes showed a multifactor interaction with a hormone-related risk factor and obesity: IGFBP5 (HRT, rs729597, p = 1.3×10-13), CGA (HRT, rs58124219, p = 1.0×10-5), and SSTR1 (parity, rs77266093, p = 1.1×10-5). Similarly, the serous histology analyses detected a multifactor interaction with HRT, obesity and IGFBP5 (rs729597, p = 8.8×10-5). Additionally, serous only analyses also detected an interaction involving ESR1 (tubal ligation, Chr6:152246911(indel), p = 8.0×10-6). Conclusions: In a very large case-control collection, assessment of multifactor gene-environment interaction is feasible, and, here, suggest that the relationship between SNPs in candidate genes and hormone-related risk factors in EOC susceptibility may vary. Citation Format: Joseph Usset, Brooke Fridley, Ellen Goode, Joellen Schildkraut, Paul Pharoah. Assessment of multifactor gene-environment interactions and ovarian cancer risk: SNPs, obesity, and hormone-related risk factors. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4684. doi:10.1158/1538-7445.AM2015-4684

Abstract Many germline genetic variants have been found to be associated with susceptibility to epithelial ovarian cancer (EOC), but currently few associations have been identified for EOC outcomes such as overall survival. In the absence of evidence for common germline markers, it is possible that rare variants not captured in GWAS arrays may be associated with overall survival. We hypothesize that rare variants, either individually or combined across a gene, may be associated with overall survival in EOC. We utilized data from two genotyping projects of the Ovarian Cancer Association Consortium (OCAC) that used commercial exome-based arrays. The primary sample, genotyped on the Affymetrix Axiom Exome Array at 227,892 standard and custom variants, consisted of 14 independent EOC studies (4293 cases; 2257 deaths). The secondary sample, genotyped on the Illumina Infinium HumanExome BeadChip at 114,620 variants, consisted of six additional EOC studies (1744 cases; 1027 deaths). Both sets were analyzed using Cox proportional hazards regression to determine the association of each variant with overall survival time. Meta-analysis was conducted across samples at 73,203 overlapping variants. Gene-level tests of over 18,000 genes were also conducted using burden tests and Sequencing Kernel Association Tests (SKAT) adapted for Cox proportional hazards models. Gene-level tests were conducted separately for the primary and secondary cases and then were meta-analyzed. Models were adjusted for age, site, and three principle components. Six individual variants had meta-analysis p-values &lt; 1.0E-4, five of which were rare (MAF &lt; 0.003). Common variant rs8170 (MAF = 0.2) in BABAM1 encoding BRISC and BRCA1 A complex member 1 - previously reported to be associated with EOC risk_showed evidence of association with overall survival (Pmeta = 1.5E-5) and consistent estimates of effect size across both samples (HR1 = 1.18, HR2 = 1.19). In imputed genotype data from over 18,000 subjects from OCAC genotyped on the COGS array, this association reaches genome-wide significance (HR = 1.15, P = 4.7E-9). Five genes had gene-level meta-analysis p-values &lt; 1.0E-4, including top ranked ATG2B (p = 1.1E-6), an autophagy gene important for apoptosis; zinc finger protein PEG3 (p = 1.8E-5), playing a role in autophagy and down-regulated in EOC; helicase-like transcription factor HLTF (p = 1.9E-5), involved in post-replication DNA repair; mitochondrial fission regulator 1-like gene FAM54B (p = 3.1E-5); and T-cell surface glycoprotein CD1E (p = 8.8E-5). With over 6000 EOC cases in the primary and secondary samples, this is the largest genome-wide assessment of rare variant genetic associations with EOC overall survival. Variant rs8170 in BABAM1, a necessary component of the BRCA1 complex, is associated with overall survival. Single variant and gene-level analyses provide complementary approaches, suggesting additional candidates that warrant follow-up study. Citation Format: Stacey J. Winham, Brooke L. Fridley, Melissa C. Larson, Zachary Fogarty, Andrew Berchuck, Yian Ann Chen, Hui-Yi Lin, Georgia Chenevix-Trench, Jenny Permuth-Wey, Thomas A. Sellers, Ailith Pirie, Ellen L. Goode, Ovarian Cancer Association Consortium. Investigation of exome variants associated with overall survival in ovarian cancer. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4636. doi:10.1158/1538-7445.AM2015-4636

Abstract Background: Genome-wide association studies (GWAS) have so far reported 12 loci associated with serous epithelial ovarian cancer (EOC) risk. We hypothesized that some of these loci function through nearby transcription factor (TF) genes and that putative target genes of these TFs as identified by coexpression may also be enriched for additional EOC risk associations. Methods: We selected TF genes within 1 Mb of the top signal at the 12 genome-wide significant risk loci. Mutual information, a form of correlation, was used to build networks of genes strongly coexpressed with each selected TF gene in the unified microarray dataset of 489 serous EOC tumors from The Cancer Genome Atlas. Genes represented in this dataset were subsequently ranked using a gene-level test based on results for germline SNPs from a serous EOC GWAS meta-analysis (2,196 cases/4,396 controls). Results: Gene set enrichment analysis identified six networks centered on TF genes (HOXB2, HOXB5, HOXB6, HOXB7 at 17q21.32 and HOXD1, HOXD3 at 2q31) that were significantly enriched for genes from the risk-associated end of the ranked list (P &lt; 0.05 and FDR &lt; 0.05). These results were replicated (P &lt; 0.05) using an independent association study (7,035 cases/21,693 controls). Genes underlying enrichment in the six networks were pooled into a combined network. Conclusion: We identified a HOX-centric network associated with serous EOC risk containing several genes with known or emerging roles in serous EOC development. Impact: Network analysis integrating large, context-specific datasets has the potential to offer mechanistic insights into cancer susceptibility and prioritize genes for experimental characterization. Cancer Epidemiol Biomarkers Prev; 24(10); 1574–84. ©2015 AACR.

Abstract Purpose: Chemotherapy resistance remains a major challenge in the treatment of ovarian cancer. We hypothesize that germline polymorphisms might be associated with clinical outcome. Experimental Design: We analyzed approximately 2.8 million genotyped and imputed SNPs from the iCOGS experiment for progression-free survival (PFS) and overall survival (OS) in 2,901 European epithelial ovarian cancer (EOC) patients who underwent first-line treatment of cytoreductive surgery and chemotherapy regardless of regimen, and in a subset of 1,098 patients treated with ≥4 cycles of paclitaxel and carboplatin at standard doses. We evaluated the top SNPs in 4,434 EOC patients, including patients from The Cancer Genome Atlas. In addition, we conducted pathway analysis of all intragenic SNPs and tested their association with PFS and OS using gene set enrichment analysis. Results: Five SNPs were significantly associated (P ≤ 1.0 × 10−5) with poorer outcomes in at least one of the four analyses, three of which, rs4910232 (11p15.3), rs2549714 (16q23), and rs6674079 (1q22), were located in long noncoding RNAs (lncRNAs) RP11-179A10.1, RP11-314O13.1, and RP11-284F21.8, respectively (P ≤ 7.1 × 10−6). ENCODE ChIP-seq data at 1q22 for normal ovary show evidence of histone modification around RP11-284F21.8, and rs6674079 is perfectly correlated with another SNP within the super-enhancer MEF2D, expression levels of which were reportedly associated with prognosis in another solid tumor. YAP1- and WWTR1 (TAZ)-stimulated gene expression and high-density lipoprotein (HDL)-mediated lipid transport pathways were associated with PFS and OS, respectively, in the cohort who had standard chemotherapy (pGSEA ≤6 × 10−3). Conclusions: We have identified SNPs in three lncRNAs that might be important targets for novel EOC therapies. Clin Cancer Res; 21(23); 5264–76. ©2015 AACR.

Abstract Background: The standard treatment for epithelial ovarian cancer (EOC) patients with advanced disease is initial debulking surgery followed by carboplatin-paclitaxel combination therapy. Unfortunately, even with modern chemotherapy, most patients with advanced disease relapse and die of EOC. One of the biggest challenges in the treatment of EOC is the large inter-patient variation in clinical response and toxicity observed for taxane and platinum therapies. Purpose: To determine pharmacogenomic markers related to carboplatin-paclitaxel therapy with 33 Epstein-Barr virus (EBV) transformed lymphoblastoid cell lines (LCLs) derived from EOC patients for which genome-wide genetic information exists. Methods: The 33 LCLs were treated with increasing concentrations of carboplatin and/or paclitaxel, followed by MTT cytotoxicity and Caspase3/7 assays. Using the measurements of MTT absorbance or Caspase3/7 fluorescence, dose response curves were fit to estimate the in vitro drug response (MTT or caspase3/7) for each LCL and treatment (paclitaxel, carboplatin, and combination) using the inflection point of the fitted curve (i.e.., the MTT IC50 or effective dose that kills 50% of the cells). The six in vitro phenotypes where then associated with genome-wide SNP data, imputed up to the markers in the 1000 Genomes Project. Results: Genetic association analysis of carboplatin, paclitaxel and combination therapy cell viability, as measured by MTT, found one SNP, rs202167095 within HS3ST4, in common with all three phenotypes (p &lt; 0.0001). HS3ST4 is regulated by TRF2 and activates natural killer cells. Furthermore, TRF2 maintains telomere integrity and is involved in activating DNA damage response (DDR) pathway.. Similarly, three SNPs were found to be associated with the three Caspase3/7 activities: SNPs on chromosome 3 (rs13098237 and rs13098239), and a SNP on chromosome 10 (rs138004421 within LINC00858). Lastly, SNPs in the following genes were associated with carboplatin IC50 with p &lt; 1x10-6 (i.e., the top associated SNPs with any drug response phenotype): JAKMIP3 (p=2.3x10-7), TSHZ3 (p = 4.9x10-7) and TGFA (p=3.6x10-7). JAKMIP3 has been implicated in some cancers and interacts with other janus kinases and JAK-STAT pathways, while TSHZ3 inhibits the expression of CASP4, thus limiting cell apoptosis. Lastly, TGFA codes for a growth factor that is a ligand for epidermal growth factor receptor, which activates a signaling pathway for cell proliferation, differentiation and development. Conclusions: Identifying genetic loci associated with response to platinum-taxane therapies will ultimately improve clinicians' ability to tailor therapy decisions for EOC patients. This research discovered novel candidate genes for which further research is needed to understand the mechanism by which they impact clinical response. Citation Format: Brooke L. Fridley, Taraswi M. Ghosh, Alice Wang, Ellen L. Goode, Jatinder Lamba. In vitro pharmacogenomic study of ovarian cancer finds novel genetic variants associated with response to platinum and taxane therapies. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 5565. doi:10.1158/1538-7445.AM2014-5565

Abstract Both genetic and epigenetic factors influence the development and progression of epithelial ovarian cancer (EOC). However, there is an incomplete understanding of the interrelationship between these factors and whether they interact to impact disease risk. Given that the development of ovarian cancer is complex, we aimed to identify DNA methylation marks that are candidates for mediating ovarian cancer genetic risk. We used 214 cases and 214 age-matched controls from the Mayo Clinic Ovarian Cancer Study. Pretreatment, blood-derived DNA was profiled for genome-wide methylation (Illumina Infinium HumanMethylation27 BeadArray) and single nucleotide polymorphisms (SNPs, Illumina Infinium HD Human610-Quad BeadArray). We employed a three-step filtering procedure, followed by the Causal Inference Test (CIT), to distinguish CpG sites that mediate genetic risk, from those that are consequential or independently acted on by genotype. Controlling for the estimated distribution of immune cells and other key covariates, we identified 1,993 out of 25,926 (7.7%) CpGs that were significantly differentially methylated between cases and controls (FDR, q &lt; 0.05). The relationship between methylation and case-control status among these 1,993 CpGs was found to be highly consistent with the results of an independent study. Implementation of the CIT test revealed 13 CpG/SNP pairs, comprising 13 unique CpGs and 17 unique SNPs, which represent potential methylation-mediated relationships between genotype and EOC risk. These findings provide additional insight into EOC etiology and may serve as novel biomarkers for EOC susceptibility. Citation Format: Devin C. Koestler, Prabhakar Chalise, Mine S. Cicek, Julie M. Cunningham, Sebastian Armasu, Melissa C. Larson, Jeremy Chien, Matthew Block, Kimberly R. Kalli, Thomas A. Sellers, Ellen L. Goode, Brooke L. Fridley. Integrative genomic analysis identifies epigenetic marks that mediate genetic risk for epithelial ovarian cancer. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 260. doi:10.1158/1538-7445.AM2014-260