Fridley, Brooke L.

Abstract An ovarian cancer Genome Wide Association Study (GWAS) identified 9p22.2 as a novel susceptibility locus with the most statistically significant single nucleotide polymorphisms (SNPs) located in an intergenic region near the Basonuclin 2 (BNC2) gene, which codes for a putative transcription regulator containing three pairs of zinc finger (ZF) domains. The minor alleles are protective in terms of ovarian cancer susceptibility. However, the molecular mechanisms by which these SNPs modify susceptibility remain largely unknown. The significant SNPs in the 9p22.2 locus lie in non-coding regions and therefore are hypothesized to affect the activity of regulatory elements and modify the expression of a target gene(s). In order to test this hypothesis we conducted fine mapping for the locus which delimited a ∼74kb region containing SNPs with a p-value less than 10-8. FAIRE-Seq and ChIP-Seq experiments for histone markers conducted in immortalized ovarian surface epithelial cells (IOSE) and fallopian tube epithelial cells (IFTE) were used to prioritize functional SNP candidates which overlap with regulatory elements. Luciferase assays tested the ability of these regulatory elements to activate transcription. Chromosome conformation capture (3C) experiments demonstrate a physical interaction between the BNC2 promoter and candidate regulatory elements containing risk-associated SNPs. Methylation Quantification at Trait Loci (mQTL) revealed an association between decreased methylation at the BNC2 promoter and the protective minor alleles. Expression analysis shows decreased expression of BNC2 in cancer versus normal tissue implicating tumor suppressor function of BNC2. Therefore the protective minor allele likely increases expression of BNC2 which in turn contains tumor suppressing properties that decrease ovarian cancer risk. BNC2 protein is in complex with transcriptional regulatory proteins, in particular the NURD complex components. Moreover, BNC2 acts as a transcriptional repressor in in vitro transfection assays indicating that it functions in transcriptional repression. We then used protein binding microarrays and CHIP-Seq experiments to identify its putative DNA binding motifs and its downstream target genes. Analysis of this data suggests that BNC2 functions in a regulatory transcription network that impacts on genes implicated in ovarian cancer with enrichment for genes in the TGF-beta response pathway. Citation Format: Melissa A. Buckley, Howard C. Shen, Gustavo A. Mendoza-Fandino, Nicholas T. Woods, Anxhela Gjyshi, Juliet French, Kate Lawrenson, Honglin Song, Jonathan Tyrer, Renato S. Carvalho, Alexandra Valle, Ann Chen, Sean Yoder, Gregory Bloom, Ya-Yu Tsai, Ally Yang, Timothy R. Hughes, Xiaotao Qu, Mine Cicek, Melissa Larson, Ellen Goode, Brooke Fridley, Susan Ramus, Georgia Chenevix-Trench, Paul Pharoah, Thomas A. Sellers, Simon Gayther, Alvaro N.A. Monteiro, Ovarian Cancer Association Consortium. Functional analysis of the 9p22 locus implicates the transcriptional regulation of BNC2 as a mechanism in ovarian cancer predisposition. [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 3285. doi:10.1158/1538-7445.AM2014-3285

Abstract Purpose: To analyze the effect of germline mutations in BRCA1 and BRCA2 on mortality in patients with ovarian cancer up to 10 years after diagnosis. Experimental Design: We used unpublished survival time data for 2,242 patients from two case–control studies and extended survival time data for 4,314 patients from previously reported studies. All participants had been screened for deleterious germline mutations in BRCA1 and BRCA2. Survival time was analyzed for the combined data using Cox proportional hazard models with BRCA1 and BRCA2 as time-varying covariates. Competing risks were analyzed using Fine and Gray model. Results: The combined 10-year overall survival rate was 30% [95% confidence interval (CI), 28%–31%] for non-carriers, 25% (95% CI, 22%–28%) for BRCA1 carriers, and 35% (95% CI, 30%–41%) for BRCA2 carriers. The HR for BRCA1 was 0.53 at time zero and increased over time becoming greater than one at 4.8 years. For BRCA2, the HR was 0.42 at time zero and increased over time (predicted to become greater than 1 at 10.5 years). The results were similar when restricted to 3,202 patients with high-grade serous tumors and to ovarian cancer–specific mortality. Conclusions: BRCA1/2 mutations are associated with better short-term survival, but this advantage decreases over time and in BRCA1 carriers is eventually reversed. This may have important implications for therapy of both primary and relapsed disease and for analysis of long-term survival in clinical trials of new agents, particularly those that are effective in BRCA1/2 mutation carriers. Clin Cancer Res; 21(3); 652–7. ©2014 AACR.

Abstract Tumor infiltrating regulatory T cells (Tregs) appear to play a role in survival in a number of solid tumors, including ovarian cancer. We assessed genetic variation via tag SNPs (N=749) in 25 Treg associated genes (CD28, CTLA4, FOXP3, IDO1, IL10, IL10RA, IL15, 1L17RA, IL23A, IL23R, IL2RA, IL6, IL6R, IL8, LGALS1, LGALS9, MAP3K8, STAT5A, STAT5B, TGFB1, TGFB2, TGFB3, TGFBR1, TGRBR2, and TGFBR3) in relation to ovarian cancer survival. SNPs were selected based on r2 ≥ 0.8, MAF ≥ 0.05, 5 kb upstream and downstream from the gene. We had SNP and outcome data for ovarian cancer from 10,084 combined ovarian cancer cases, including 5,248 high grade serous, 1,452 endometrioid, 661 mucinous, and 795 clear cell cases of European descent across multiple studies from the Ovarian Cancer Association Consortium (OCAC). Cox proportional hazards regression modeling was used to estimate per-allele hazard ratios (HRs) and 95% confidence intervals (CIs) for associations with overall survival, adjusting for age at diagnosis, population substructure PCs, study site, histology (for the all case analysis), tumor stage, tumor grade and oral contraceptive use. The strongest associations were found for endometrioid histologic subtype and IL2RA SNPs (maximum r2=0.67) rs11256497 [HR=1.42, 95% CI: 1.22-1.64; p=5.7 x 10−6], rs791587 [HR=1.36, 95% CI: 1.17-1.57; p=6.2 x 10−5], rs2476491 [HR=1.40, 95% CI: 1.19-1.64; p=5.6 x 10−5], and rs10795763 [HR=1.35, 95% CI: 1.17-1.57; p=7.9 x 10−5], as well as clear cell ovarian cancer and CTLA4 SNP rs231775 [HR=0.67, 95% CI: 0.54-0.82; p=9.3 x 10−5]. No other associations were observed at p<1 x 10−4. Overall, this study provided evidence that variation in genes related to Tregs, particularly IL2RA and CTLA4, influence ovarian cancer survival. Citation Format: Bridget Charbonneau, Keith Knutson, Robert A. Vierkant, Julie M. Cunningham, Zachary C. Fogarty, Kimberly R. Kalli, Matthew J. Maurer, Ann L. Oberg, Brooke L. Fridley, Claudia C. Preston, Paul D.P. Pharoah, Susan J. Ramus, Catherine M. Phelan, Kunle Odunsi, Kirsten Moysich, Ellen L. Goode. Polymorphisms in regulatory T cell related genes and ovarian cancer survival. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 4844. doi:10.1158/1538-7445.AM2013-4844

Abstract BACKGROUND Metformin, a widely used anti-diabetic drug, is being considered as a highly promising agent for treatment and prevention of many types of cancer, including breast cancer. Although the exact mechanism of action of metformin is unknown, it is thought that metformin is an AMPK activator. To better individualize metformin therapy and explore additional underlying molecular mechanism associated with metformin action, we conducted a metformin pharmacogenomic study using 266 genomic data-enriched lymphoblastoid cell lines (LCLs). METHODS Genome-wide microarray data were generated for all LCLs, including 1.3 million SNPs, 485K DNA methylation probes and 54K mRNA expression probe sets. Metformin cytotoxicity assay was performed using MTS assay. Genome-wide association (GWA) analyses and integrated analyses using all the genomic data and metformin cytotoxicity (IC50 value) were performed to identify top candidate genes for functional validation in breast cancer cell lines, MDA-MB 231 and Hs578T, using siRNA screening approach. Further mechanistic study was performed to investigate the role of candidates in regulation of AMPK activity. RESULTS GWA analyses in LCLs identified 198 mRNA expression probe sets and 210 DNA methylation probes associated with metformin IC50 with p-value <10-4, respectively. Analysis of 1.3 million SNPs found 12 loci (a region containing at least 1 SNP of p-value <10-5 and 1 SNP of p-value <10-3 within 50kb) that were associated with metformin IC50. Integrated SNP loci-mRNA expression-IC50 analysis indicated that the SNP loci on chromosome 16 was associated with metformin IC50 through a trans-regulation of expression of 7 genes with p-value <10-4. Integrated methylation-mRNA expression-IC50 analysis showed that 15 DNA methylation probes in 6 genes were associated with metformin IC50 through both cis- and trans-regulation of expression of 48 genes with p-value <10-4, including 4 methylation probes in 3 of those 6 genes that were cis-correlated with its own gene expression. Ingenuity pathway analysis of top 62 candidate genes identified four major pathways and 55 genes within those pathways were selected for functional validation in breast cancer cell lines. The knockdown experiment showed that down regulation of 25 genes significantly altered metformin sensitivity in MDA-MB231, and 14 of them also showed the same effect in Hs578T. Preliminary mechanistic experiment indicated that knockdown of 7 genes significantly reduced AMPK activity in MDA-MB231 and 3 of them had same effect in Hs578T. CONCLUSIONS GWAS using a genomic data-enriched LCL model system, together with functional validation in breast cancer cell lines, could help us to identify novel genetic and epigenetic biomarkers involved in metformin response and help us to better understand the mechanisms of metformin in cancer treatment. Citation Format: Nifang Niu, Xianglin Tan, Brooke L. Fridley, Daniel J. Schaid, Ryan P. Abo, Anthony Batzler, Erin E. Carlson, Gregory Jenkins, Sebastian Moran, Holger Andrea Heyn, Manel Esteller Badosa, Liewei Wang. Metformin pharmacogenomics: A genome-wide associate study to identify genetic and epigenetic biomarkers involved in metformin response. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 2271. doi:10.1158/1538-7445.AM2013-2271

Abstract The FKBP5 gene encodes the FKBP51 immunophilin, a negative regulator of AKT. Variability in level of FKBP5 expression is a major factor contributing to variation in response to several classes of chemotherapeutic agents including gemcitabine, first line treatment for pancreatic cancer. Genetic variation of FKBP5 could influence its expression, function, and, ultimately, treatment response of pancreatic cancer. We set out to comprehensively study the role of genetic variation in FKBP5 identified by Next Generation DNA resequencing in response to gemcitabine treatment of pancreatic cancer by utilizing both tumor and germline DNA samples from 43 pancreatic cancer patients, including 19 paired tumor-normal samples. Genotype-phenotype association studies with overall survival (in tumor only) as well as with FKBP5 gene expression were performed using the same samples in which Next Generation DNA resequencing had been performed, followed by functional genomic studies of top selected single nucleotide polymorphisms (SNPs). In-depth resequencing identified 404 FKBP5 SNPs in both normal and tumor DNA patient samples. SNPs with the smallest p-values for association with FKBP5 survival or expression were subjected to functional genomic studies. Electromobility shift assay (EMSA) showed that the rs73748206 “T” allele altered DNA-protein binding patterns, consistent with a significantly increased reporter gene activity. The effect of rs73748206 was also confirmed by mRNA expression studies using lymphoblastoid cell lines (LCL) derived from the same patients from whom DNA used for resequencing had been obtained. This comprehensive FKBP5 resequencing study provides insight into the role of inheritance in variation in response to gemcitabine therapy of pancreatic cancer. Citation Format: Katarzyna A. Ellsworth, Bruce W. Eckloff, Liang Li, Irene Moon, Brooke L. Fridley, Greg D. Jenkins, Erin Carlson, Abra Brisbin, Ryan Abo, William Bamlet, Gloria Petersen, Eric D. Wieben, Liewei Wang. Contribution of FKBP5 genetic variation to gemcitabine treatment and survival in pancreatic cancer. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 2209. doi:10.1158/1538-7445.AM2013-2209

ScopeWe reevaluated previously reported associations between variants in pathways of one‐carbon (1‐C) (folate) transfer genes and ovarian carcinoma (OC) risk, and in related pathways of purine and pyrimidine metabolism, and assessed interactions with folate intake.Methods and resultsOdds ratios (OR) for 446 genetic variants were estimated among 13 410 OC cases and 22 635 controls, and among 2281 cases and 3444 controls with folate information. Following multiple testing correction, the most significant main effect associations were for dihydropyrimidine dehydrogenase (DPYD) variants rs11587873 (OR = 0.92; p = 6 × 10−5) and rs828054 (OR = 1.06; p = 1 × 10−4). Thirteen variants in the pyrimidine metabolism genes, DPYD, DPYS, PPAT, and TYMS, also interacted significantly with folate in a multivariant analysis (corrected p = 9.9 × 10−6) but collectively explained only 0.2% of OC risk. Although no other associations were significant after multiple testing correction, variants in SHMT1 in 1‐C transfer, previously reported with OC, suggested lower risk at higher folate (pinteraction = 0.03–0.006).ConclusionVariation in pyrimidine metabolism genes, particularly DPYD, which was previously reported to be associated with OC, may influence risk; however, stratification by folate intake is unlikely to modify disease risk appreciably in these women. SHMT1 SNP‐by‐folate interactions are plausible but require further validation. Polymorphisms in selected genes in purine metabolism were not associated with OC.

Abstract Survival in epithelial ovarian cancer (EOC) is influenced by the host immune response, yet the key genetic determinants of inflammation and immunity that affect prognosis are not known. The nuclear factor-κB (NF-κB) transcription factor family plays an important role in many immune and inflammatory responses, including the response to cancer. We studied common inherited variation in 210 genes in the NF-κB family in 10,084 patients with invasive EOC (5,248 high-grade serous, 1,452 endometrioid, 795 clear cell, and 661 mucinous) from the Ovarian Cancer Association Consortium. Associations between genotype and overall survival were assessed using Cox regression for all patients and by major histology, adjusting for known prognostic factors and correcting for multiple testing (threshold for statistical significance, P < 2.5 × 10−5). Results were statistically significant when assessed for patients of a single histology. Key associations were with caspase recruitment domain family, member 11 (CARD11) rs41324349 in patients with mucinous EOC [HR, 1.82; 95% confidence interval (CI), 1.41–2.35; P = 4.13 × 10−6] and tumor necrosis factor receptor superfamily, member 13B (TNFRSF13B) rs7501462 in patients with endometrioid EOC (HR, 0.68; 95% CI, 0.56–0.82; P = 2.33 × 10−5). Other associations of note included TNF receptor–associated factor 2 (TRAF2) rs17250239 in patients with high-grade serous EOC (HR, 0.84; 95% CI, 0.77–0.92; P = 6.49 × 10−5) and phospholipase C, gamma 1 (PLCG1) rs11696662 in patients with clear cell EOC (HR, 0.43; 95% CI, 0.26–0.73; P = 4.56 × 10−4). These associations highlight the potential importance of genes associated with host inflammation and immunity in modulating clinical outcomes in distinct EOC histologies. Cancer Epidemiol Biomarkers Prev; 23(7); 1421–7. ©2014 AACR.

Abstract To reveal biologic mechanisms underlying clinical outcome of high-grade serous (HGS) epithelial ovarian carcinomas (EOC), we evaluated the association between tumor epigenetic changes and time to recurrence (TTR). We assessed methylation at approximately 450,000 genome-wide CpGs in tumors of 337 Mayo Clinic (Rochester, MN) patients. Semi-supervised clustering of discovery (n = 168) and validation (n = 169) sets was used to determine clinically relevant methylation classes. Clustering identified two methylation classes based on 60 informative CpGs, which differed in TTR in the validation set [R vs. L class, P = 2.9 × 10−3, HR = 0.52; 95% confidence interval (CI), 0.34–0.80]. Follow-up analyses considered genome-wide tumor mRNA expression (n = 104) and CD8 T-cell infiltration (n = 89) in patient subsets. Hypomethylation of CpGs located in 6p21.3 in the R class associated with cis upregulation of genes enriched in immune response processes (TAP1, PSMB8, PSMB9, HLA-DQB1, HLA-DQB2, HLA-DMA, and HLA-DOA), increased CD8 T-cell tumor infiltration (P = 7.6 × 10−5), and trans-regulation of genes in immune-related pathways (P = 1.6 × 10−32). This is the most comprehensive assessment of clinical outcomes with regard to epithelial ovarian carcinoma tumor methylation to date. Collectively, these results suggest that an epigenetically mediated immune response is a predictor of recurrence and, possibly, treatment response for HGS EOC. Cancer Res; 74(11); 3084–91. ©2014 AACR.