Cunningham, Julie M.

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 Introduction Chronic lymphocytic leukemia (CLL) is a B-cell malignancy that is known to have a familial component to disease risk. Although 31 loci have been found to be associated with CLL risk, the functional variant(s) driving these associations is mostly unknown. Here we set out to identify rare, highly-penetrant, cosegregating, susceptibility variants within the known GWAS discovered loci using whole exome sequencing (WES) data in CLL families from the Mayo Clinic family study of B-cell malignancies. Methods We performed WES on germ line DNA of 93 CLL families with two or more members with CLL, using Agilent capture kits and Illumina HiSeq2000. Bioinformatics analyses leveraged the following software packages: Novoalign, Picard, The Genome Analysis Toolkit (GATK), and the Biological Reference Repository (bioR). Quality control filters were implemented; subjects with mis-specified relationships were removed, as were variants with &lt;75% call rate, &lt;8X coverage, and those identified as sequencing artifacts. Each GWAS locus was defined by +/- 1Mb of the top GWAS hit within the locus. Linkage disequilibrium (LD) was calculated among the single nucleotide variants (SNVs) located within each locus. Potentially functional SNVs were identified based on: a) uncommon in public databases (&lt; 5%), b) cosegregating in at least two CLL families, c) being highly conserved and in coding regions, and d) functional prediction status of deleterious (SIFT Score), damaging (PolyPhen Score), and a moderate, or high variant impact (SNP Effect). Results In our 93 CLL families, we sequenced 443 individuals: 160 with CLL, 73 with monoclonal B cell lymphocytosis (MBL), and 210 relatives that were not diagnosed with CLL or MBL at the time of sequencing. Median age of CLL diagnosis was 59 years (range 34-87), and 56% were male. Among the MBL individuals and relatives, the median age at recruitment was 55 years (range 18-93), and 40% were male. A total of 317,666 SNVs passed our sequencing quality control filters of which 10,731 were within +/- 1 Mb of known GWAS hits from 31 loci. Of these SNVs, 91% were in coding regions, 18% were reported to have high or moderate impact, 6% were estimated to be damaging and 6% were predicted to be deleterious. From these SNVs, we identified 76 putatively functional variants distributed across 25 GWAS loci that were cosegregating in the individuals with CLL or MBL in multiple CLL pedigrees. These SNVs were all located in coding regions with high or moderate impact and were predicted to be damaging and deleterious. Of these 76 variants, 56 had a frequency of &lt;0.005 in 1000 Genomes’ European population while the remaining 20 had a frequency of 1%. Conclusions Through WES, we identified a number of rare, penetrant and potentially predisposing SNVs located within 25 of the 31 CLL GWAS-discovered loci. These segregating variants provide a list for future validation and functional studies. Citation Format: Sara Beiggi, Daniel R. O'Brien, Sara J. Achenbach, Kari G. Chaffee, Timothy G. Call, Neil E. Kay, Tait D. Shanafelt, Julie Cunningham, James R. Cerhan, Celine M. Vachon, Susan L. Slager. Cosegregating variants in chronic lymphocytic leukemia (CLL) families that are located in loci discovered by genome wide association studies (GWAS). [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 2764. doi:10.1158/1538-7445.AM2015-2764

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