Permuth‐Wey, Jennifer

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 Background: Recent bioinformatic analyses and high-throughput sequencing efforts have revealed that RNA editing is a widespread post-transcriptional modification that affects non-coding repetitive elements throughout the genome. This modification, which involves the conversion of adenosine into inosine, is mediated by the action of the adenosine deaminases acting on RNA (ADAR) family of enzymes. Based on reports of dysfunctional ADAR expression in epithelial ovarian cancers (EOC), we hypothesized that single nucleotide polymorphisms (SNPs) in five ADAR genes (ADAD1, ADAR1, ADAR2, ADAR3, SND1) may contribute to EOC susceptibility, potentially by altering the expression of corresponding genes in ovarian tissues. Methods: Using directly genotyped and imputed data from 10,891 invasive EOC cases and 21,693 healthy controls of European ancestry represented in the international Collaborative Oncological Gene-environment Study, we evaluated the frequency of 7,133 SNPs in 5 RNA editing genes. Unconditional logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI) between genotypes and case status under a log-additive genetic model, with adjustment for the first five principal components representing European ancestry. Subgroup analysis was conducted for the four main histologic sub-types of EOC. We also aggregated SNP-level data and conducted gene-level analyses using the Admixture Maximum Likelihood (AML) test and the Sequence-Kernel Association test for the combined effect of common and rare variants (SKAT-CR). Results: Association analysis revealed a top risk-associated SNP (OR (95% CI) = 1.39 (1.17-1.64)) in ADAR3 with P = 1.3 x10−4, 17 additional SNPs in ADAR1, ADAR3, or SND1 with P = 10−4, and 51 SNPs with P = 10−3; this observed number of statistically-significant SNPs was greater than expected. When restricting to the 6,500 invasive serous adenocarcinomas, the magnitude of association strengthened for numerous SNPs. Exploratory analysis for the less common histologic subtypes (endometrioid (n = 1,439), mucinous (n = 696), and clear cell (n = 660)) revealed SNP-level associations (P&lt;10−3) unique to each sub-type. Gene-level analyses based on AML and SKAT-CR revealed that ADAR1 was globally associated (P&lt;0.10) with susceptibility to invasive EOC, with P = 0.024 and P = 0.023, respectively. Expression quantitative trait locus analysis is underway to evaluate associations between genotype and gene expression in EOC and precursor tissues using data from The Cancer Genome Atlas Project and other sources. Conclusions: Preliminary data from this large-scale collaboration suggest that germline variation involving RNA editing genes may influence susceptibility to EOC. This novel finding warrants further investigation, as inherited and acquired alterations affecting RNA editing may serve as important biological mechanisms that promote the development of EOC and possibly other malignancies. Citation Format: Jennifer Permuth-Wey, Brett M. Reid, Y. Ann Chen, Hui-Yi Lin, Alvaro Monteiro, Zhihua Chen, Andrew Berchuck, Georgia Chenevix-Trench, Jennifer Doherty, Simon Gayther, Ellen Goode, Edwin Iversen, Leigh Pearce, Paul D. P. Pharoah, Catherine Phelan, Susan Ramus, Mary Anne Rossing, Joellen Schildkraut, Thomas Sellers, on behalf of the Ovarian Cancer Association Consortium. Inherited variants affecting RNA editing may contribute to ovarian cancer susceptibility. [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 4635. doi:10.1158/1538-7445.AM2015-4635

Abstract Background: Cell-type-specific super enhancers, or large clusters of mediator-rich enhancers, are highly transcribed regions of the genome with key functions in the maintenance of cell identity and control. These regulatory elements are easily perturbed and enriched for disease-associated sequence variation. Identification of disease-associated super-enhancers and investigation of their sequence variants may help identify important, biologically relevant associations previously overlooked by genome-wide studies. Thus, we sought to identify variants within ovarian tissue specific super-enhancers and investigate their association with epithelial ovarian cancer (EOC) susceptibility. Methods: We utilized genotype data from a European population of ∼11,000 invasive EOC cases and 22,000 healthy controls from the Ovarian Cancer Association Consortium. A public catalog of 478 super-enhancers, identified by peak histone H3K27ac ChIP-seq signal in ovarian tissue, was used to identify 72,116 variants located within 344 super-enhancer regions. Unconditional logistic regression (log-additive model) was used to assess individual SNP-EOC risk associations. Additionally, we reviewed associations within 500 kb of super-enhancer regions to identify proximal associated variants. When LD structure indicated, we employed conditional analysis to determine independence of the super-enhancer signal. Gene-EOC risk associations were assessed with the Admixture Maximum Likelihood (AML) test and Sequence-Kernel Association Test (SKAT) for the combined effect of common and rare (MAF&lt;0.01) variants. Models were adjusted for population stratification using the first five principal components. To adjust for multiple comparisons, a false discovery rate of 10% was used. Results: The most significant SNP observed was rs147487657 (p = 7.7×10−6) on locus 2q37.3 in the HDAC4 gene and there were an additional 12 SNPs with p&lt;1.0×10−4, although not significant at 10% FDR. Gene-level analyses identified two super-enhancers significantly associated with EOC risk. An enhancer at the previously reported 2q31 locus, containing HOXD8 and HOXD-AS2, associated (SKAT p = 4.3×10−5; AML p&lt;0.001) with the top SNP located within the 5′UTR of HOXD8 transcription factor (rs847163 (T&gt;C); p&lt;2.08×10−5). Conditional analysis did not identify an independent association within the super-enhancer region. The second EOC-associated super-enhancer identified was at locus 1q32 (SKAT p = 0.00052, AML p = 0.001) containing the CSRP1 gene. There were 6 SNPs with p&lt;0.001 and 16 SNPs with p&lt;0.01, all were inversely associated with risk and consistent with the top hit at rs927901 T&gt;C (OR = 0.94; 95% CI: 0.90-0.97; p = 2.1×10−4). Associations varied by histologic subtype, possibly due to variance in disease relevant super-enhancer cell-type.Conclusion: These findings suggest germline variation within super-enhancer regions may contribute to EOC susceptibility Citation Format: Brett M. Reid, Jennifer Permuth-Wey, Y. Ann Chen, Hui-Yi Lin, Alvaro Monteiro, Zhihua Chen, Andrew Berchuck, Georgia Chenevix-Trench, Jennifer Doherty, Simon Gayther, Ellen L. Goode, Edwin Iversen, Leigh Pearce, Paul Pharoah, Catherine Phelan, Susan Ramus, Mary Anne Rossing, Joellen Schildkraut, Thomas Sellers, on behalf of the Ovarian Cancer Association Consortium. Variants within super-enhancer regulatory elements associate with epithelial ovarian cancer risk. [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 4634. doi:10.1158/1538-7445.AM2015-4634

Abstract Background: Genome wide association studies (GWAS) have identified 20 loci associated with epithelial ovarian cancer (EOC) risk. Additional EOC risk loci await identification, and biological approaches may represent a useful strategy for overcoming sample size limitations. The ENCODE project recently concluded that only ∼1.2% of the genome encodes proteins, but at least 20% exhibits biological function and over 80% exhibits biochemical indices of function. Post-GWAS follow-up studies have firmly established that some associations are due to inter-individual differences in non-coding RNAs (ncRNAs). To estimate the magnitude of impact of variants in long non-coding RNAs (lncRNAs), we assessed enrichment of EOC-associated SNPs in lncRNA regions by comparing the density of EOC-associated loci between lncRNA regions, non-lncRNA regions, and the whole genome. Methods: The study population included ∼18,000 invasive EOC cases and 34,000 healthy controls of European ancestry from the Ovarian Cancer Association Consortium with GWAS data imputed to 1000 Genomes Project (1KGP) density. Density was compared with two measures: average chromosome length required for one EOC-associated SNP (kb/locus) and average number of tested SNPs containing one EOC-associated SNP (SNPs/locus). Coordinates and annotation for 13,870 lncRNA genes were downloaded from the publicly available GENCODE (v19) database and were used to annotate SNPs from the 1KGP imputed GWAS data, yielding 13,192 lncRNA genes with biallelic variants imputed at rsquare ≥0.25. A SNP was considered EOC-associated if it reached a significance level of P&lt;10−5. Results: Genome-wide, 11.6% of the 15,123,646 tested SNPs fell within lncRNA regions. Of 5,294 EOC-associated SNPs that were identified, 27.7% were located within lncRNAs, a higher frequency than expected (&lt;v:shape style = “WIDTH: 16.5pt; HEIGHT: 18pt” id = _x0000_i1026 equationxml = ‘14pχ2′ type = “#_x0000_t75”&gt; = 9.35×10−22). Additionally, EOC-associated SNPs were located on average every 237.6 kb (or 1204 tested SNPs) in lncRNA regions versus every 657.7 kb (or 3489 SNPS) in non-lncRNA regions and every 541.5 kb (or 2857 SNPs) in the whole genome, clearly demonstrating an enrichment of EOC-associated loci in lncRNA regions. A total of 1,464 lncRNA SNPs from 53 unique lncRNAs associated with EOC risk (P&lt;10−5); 70 of the identified lncRNA SNPs had P-values between 10−20 and 10−32, and 221 had P-values between 10−10 and 10−15, representing 21 unique regions &gt; 500kb apart. Conclusions: These results justify efforts to integrate functional and population level genetic data on lncRNA regions to explore a new paradigm that addresses the problem of the missing heritability for cancer where restrictions in sample size and statistical power limit the ability to discover additional risk loci. Citation Format: Thomas A. Sellers, Brett M. Reid, Y. Ann Chen, Hui-Yi Lin, Edward Richards, Jamie Teer, Alvaro Monteiro, Zhihua Chen, Andrew Berchuck, Georgia Chenevix-Trench, Jennifer Doherty, Ellen Goode, Edwin Iverson, Leigh Pearce, Paul Pharoah, Catherine Phelan, Susan Ramus, Mary Anne Rossing, Joellen Schildkraut, Jin Cheng, Simon Gayther, Jennifer Permuth-Wey, on behalf of the Ovarian Cancer Association Consortium. Evidence that long non-coding RNA variants associate with epithelial ovarian cancer risk. [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 4633. doi:10.1158/1538-7445.AM2015-4633

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 Background: Common genetic variants identified by genome-wide association studies and follow-up genotyping initiatives only explain a small proportion of the heritability of epithelial ovarian cancer (EOC). We hypothesized that some of the missing heritability may be explained by rare (minor allele frequency (MAF) &lt;0.5%) and low-frequency (MAF 0.5-5%) protein-coding variants that are not well captured by genotyping arrays that focus on common variants. Methods: To examine the association between uncommon variants and EOC susceptibility, we genotyped 7,060 EOC cases and 6,712 controls from the Ovarian Cancer Association Consortium using the Affymetrix Axiom Exome Genotyping Array. This platform contains 300,785 putative functional coding variants selected from exome sequences of 12,000+ individuals representing multiple ethnicities and complex traits. After performing quality control (QC), we used unconditional logistic regression treating the number of minor alleles carried as an ordinal variable (log-additive model) to evaluate SNP-EOC risk associations. Genotype cluster plots were visually inspected for top-ranking variants to determine assay success, and allele frequencies of samples from the 1000 Genomes Project and other cohorts were evaluated for QC purposes. Results: A total of 13,318 samples (97%) passed QC, and 12,779 eligible samples (6,293 invasive serous and endometrioid cases and 6,486 controls) were included in the analyses. 280,737 (93%) of the 300,785 attempted markers were successfully genotyped and had call rates &gt;90%. After removing 31,436 monomorphic markers, the 249,301 remaining markers were included in association analyses. We confirmed variants at previously reported EOC susceptibility loci (9p22, 3q25, 17cen-q21.3, and 8q24), with p-values between 10-5 and 10-11. Several correlated low-frequency non-synonymous variants (MAF=5%) residing in a novel EOC susceptibility gene at 3q25.31 were risk-associated. The most significant variant at 3q25.31 had an OR=1.36 and p = 5.58x10-8. Also noteworthy is a rare variant at 15q15.1-q21.1 (MAF=0.5%) that introduces a stop codon in a gene belonging to the tyrosine kinase family (OR=0.39, CI=0.26-0.58, p = 4.52x10-6). We also identified a more common synonymous coding variant in a phosphatase at 6q21.3 that was associated with increased risk (OR=1.22, CI=1.13-1.33, p = 1.84x10-6, MAF=9.6%). Population stratification and gene-based analyses are underway to further examine the impact of the newly-identified variants on EOC risk. A meta-analysis to include data from an additional 2104 EOC cases and 2516 controls genotyped with Illumina's HumanExome Beadchip rare variant array is also underway. Functional analysis of the most promising variants will follow. Conclusions: Preliminary data from this large-scale study reveals novel rare and low-frequency variants that may influence susceptibility to epithelial ovarian cancer. Citation Format: Jennifer Permuth-Wey, Y. Ann Chen, Zhihua Chen, Andrew Berchuck, Georgia Chenevix-Trench, Jennifer Doherty, Simon Gayther, Ellen L. Goode, Edwin Iversen, Alvaro N.A. Monteiro, Leigh Pearce, Paul D.P. Pharoah, Catherine M. Phelan, Ailith Pirie, Susan Ramus, Mary Ann Rossing, Joellen M. Schildkraut, Thomas A. Sellers, on behalf of the Ovarian Cancer Association Consortium. Exome genotyping array identifies rare and low-frequency variants that may be associated with ovarian cancer risk. [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 946. doi:10.1158/1538-7445.AM2014-946

Abstract Prostate cancer is the second most common malignancy among men worldwide. Genome-wide association studies have identified 100 risk variants for prostate cancer, which can explain approximately 33% of the familial risk of the disease. We hypothesized that a comprehensive analysis of genetic variations found within the 3′ untranslated region of genes predicted to affect miRNA binding (miRSNP) can identify additional prostate cancer risk variants. We investigated the association between 2,169 miRSNPs and prostate cancer risk in a large-scale analysis of 22,301 cases and 22,320 controls of European ancestry from 23 participating studies. Twenty-two miRSNPs were associated (P &lt; 2.3 × 10−5) with risk of prostate cancer, 10 of which were within 7 genes previously not mapped by GWAS studies. Further, using miRNA mimics and reporter gene assays, we showed that miR-3162-5p has specific affinity for the KLK3 rs1058205 miRSNP T-allele, whereas miR-370 has greater affinity for the VAMP8 rs1010 miRSNP A-allele, validating their functional role. Significance: Findings from this large association study suggest that a focus on miRSNPs, including functional evaluation, can identify candidate risk loci below currently accepted statistical levels of genome-wide significance. Studies of miRNAs and their interactions with SNPs could provide further insights into the mechanisms of prostate cancer risk. Cancer Discov; 5(4); 368–79. ©2015 AACR. See related commentary by Yousef, p. 351 This article is highlighted in the In This Issue feature, p. 333