Phelan, Catherine M.

Abstract Background: Epithelial Ovarian Cancer (EOC) is a lethal gynecologic malignancy. Mitochondria are cellular organelles that regulate energy metabolism and are the primary source of reactive oxygen species (ROS), which damage both the nuclear and mitochondria DNA (mtDNA), contributing to the pathogenesis of EOC. Previously, we identified statistically significant EOC risk associations with 54 single nucleotide polymorphisms (SNPs) from 11 mitochondria-related genes in a study of 1815 EOC cases and 1900 controls. In this replication study we expand our investigation of mitochondrial metabolism genetic variants and associations with EOC risk and overall survival (OS). Methods: In a discovery set of 3761 EOC cases and 2722 controls, we investigated risk associations in 4200 SNPs from 130 mitochondrial metabolism genes. We genotyped 431 of the most significant SNPs from 86 genes in 14,525 EOC cases and 23,447 controls (of European ancestry) from 43 studies in OCAC using an Illumina Infinium iSelect BeadChip. Risk analyses were conducted using unconditional logistic regression under a log-additive model and adjusted for study site and the first five principal components of population substructure. Both invasive cancers combined and the four main histological subtypes were analyzed. Survival analyses were conducted on the same SNPs using Cox hazards models to estimate hazard ratios (HRs) and 95% CIs considering overall survival in 10,084 cases from 28 OCAC studies. We adjusted for multiple tests using False Discovery Rate q<20%. Results: The strongest evidence of an association for risk of invasive cancers combined and the serous subtype was SNP rs3096 in the pyridoxamine 5'-phosphate oxidase (PNPO) gene (OR=1.08, p=2.8x10-6). In clear cell carcinoma, the strongest association was a missense coding SNP rs28379632 (Ile421Met) in the MT-CYB gene (OR=4.0, p=0.003). For endometrioid cancer rs9330440 in the UNQ6494 gene was most significant (OR=0.89, p=0.001) while rs10267797 in the dihydrolipoamide dehydrogenase (DLD) gene was most significant in the mucinous subtype (OR=1.15, p=0.004). The strongest association with overall survival was mitochondrial transcription termination factor 1 (MTERF) rs2540571 (HR=0.92, p= 0.00011). Conclusion: These results suggest that mitochondrial metabolism gene variants appear to be associated with EOC risk and survival. Citation Format: Ganna Chornokur, Hui-Yi Lin, Alvaro Monteiro, Joellen Schildkraut, Ellen L. Goode, Paul DP Pharoah, Simon A. Gayther, Susan Ramus, Thomas A. Sellers, Catherine M. Phelan. Mitochondrial metabolism gene variants are associated with risk and survival in epithelial ovarian cancer (EOC). [abstract]. In: Proceedings of the AACR Special Conference: Metabolism and Cancer; Jun 7-10, 2015; Bellevue, WA. Philadelphia (PA): AACR; Mol Cancer Res 2016;14(1_Suppl):Abstract nr A50.

Abstract Background: Epithelial Ovarian Cancer (EOC) is a lethal gynecologic malignancy and the fifth cause of cancer mortality in women in the US. Normal ovarian physiology is intricately connected to tightly regulated small GTP binding proteins of the Ras superfamily (Ras, Rac, Rho, Rab, Arf, and Ran) which regulate key cellular processes such as signal transduction, cell proliferation, cell motility, and vesicle transport. These proteins function in a highly coordinated manner through signaling cascades and feedback loops within and among the small GTPase subfamilies. We hypothesized that single nucleotide polymorphisms (SNPs) in small GTPase genes are associated with epithelial ovarian cancer (EOC) risk and aberrant expression of these genes correlates with tumor characteristics including overall survival (OS). Methods: In a discovery set of 7931 EOC cases and 9206 controls we investigated 9,356 SNPs from 112 genes, 657 of which showed associations up to the significance level of p<0.05. We genotyped 407 of the most significant SNPs from 112 genes in a combined dataset (discovery and replication) which consisted of 18,736 EOC cases and 23,448 controls (of European ancestry) from 43 studies in the Ovarian Cancer Association Consortium (OCAC), using an Illumina Infinium iSelect BeadChip as part of the Collaborative Oncological Gene-environment Study (COGS). Odds ratios and 95% confidence intervals were calculated using unconditional logistic regression under log-additive models. A False Discovery Rate (FDR) q<0.2 was applied. Pearson's correlation tests were performed on the expression of genes with tumor characteristics and OS in 561 ovarian tumors: serous (n=464); endometrioid (n=52); clear cell carcinoma (n=28) and mucinous (n=17) subtypes. Results: The most significantly associated SNP associations with EOC were: AKAP6 rs1955513, (OR=0.9, p=3.3x10-4) in all invasive; ARGHEF10L rs10788679, (OR=1.05, p=2.6x10-4) in serous; RAB31 rs1166373, (OR=1.25, p=4.0x10-3) in endometrioid; KRAS rs4963872, (OR=1.35, P=4.5x10-4) in mucinous EOC and TNIK rs6780532, (OR=0.88, p=7.7x10-3) in clear cell carcinoma. Increased expression of AKAP6 and TNIK was marginally correlated with advanced EOC stage (p<0.05) while high KRAS expression was correlated with OS (HR=1.18, p=0.038). Conclusions: Genetic variation in the small GTPase genes appears to be associated with ovarian cancer risk and aberrant expression of these genes correlates with tumor characteristics and OS. Citation Format: Hui-Yi Lin, Yin Xiong, Jonathan Tyrer, Douglas C. Marchion, Alvaro NA Monteiro, Andrew Berchuck, Joellen M. Schildkraut, Ellen L. Goode, Susan J. Ramus, Simon A. Gayther, Paul DP Pharoah, Steven A. Narod, Thomas A. Sellers, Catherine M. Phelan. Investigation of small GTPase genes in epithelial ovarian cancer. [abstract]. In: Proceedings of the AACR Special Conference on Advances in Ovarian Cancer Research: Exploiting Vulnerabilities; Oct 17-20, 2015; Orlando, FL. Philadelphia (PA): AACR; Clin Cancer Res 2016;22(2 Suppl):Abstract nr B27.

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: 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 Background Genome-wide association studies (GWAS) previously identified 8 new risk loci for all subtypes of epithelial ovarian carcinoma (EOC) and one specific to invasive serous EOC. One locus, at 2q31, was also associated with invasive mucinous EOC (mEOC) at P=7.3 x 10-7 (Nat Genet 2010:874). Here, we report new risk loci for mEOC, which is thought to develop from a borderline mucinous tumor. These loci may provide clues to the molecular basis of the disease. Methods We used genotype data from 4 GWAS of EOC and from multiple case-control studies in the Ovarian Cancer Association Consortium that had been genotyped using the iCOGS custom SNP panel (210k SNPs) developed for the Collaborative Oncological Gene-environment Study. We imputed genotypes of common variation across the genome using the 1000 Genomes Project reference panel and tested for association using log-additive unconditional regression models. Results Genotypes were available for 1,643 combined invasive and borderline mEOC and for 1,002 invasive-only mEOC of European ancestry. Three loci in new regions, and one locus in a previously published region, with minor allele frequency (MAF) &gt; 1% and P ≤ 10-7 were found. The imputed SNP rs72831838 (imputation r2=0.67, MAF=0.15) is located at 2q13 in an intron of the PAX8 (paired box 8) gene (all mEOC OR=1.35, P=4.2 x 10-8 and invasive-only OR=1.43, P=9.6 x 10-8). PAX8 is a transcription factor and may be a Müllerian lineage marker that distinguishes carcinomas of gynecologic origin from other sites. The imputed variant at chr19:39745695 (imputation r2=0.61, MAF=0.32) is an insertion/deletion polymorphism located at 19q13.2 and falls within the MSRB1P1 (methionine sulfoxide reductase B1 pseudogene 1). This locus was associated with decreased risk (all mEOC OR=0.74, P = 2.1 x 10-9 and invasive-only OR=0.79, P = 8.1 x 10-5) and is enriched with the H3K27Ac histone mark that identifies active enhancers of gene expression. The imputed SNP rs6425050 (imputation r2=0.45, MAF=0.19) is located at 1q25 (all mEOC OR=1.38, P=4.7 x 10-8 and invasive-only OR=1.41, P=1.6x10-6). This SNP is intronic to PLA2G4A (phospholipase A2, group IVA), which encodes an enzyme that hydrolyzes membrane phospholipids for eicosanoid metabolism and signaling. None of these loci were associated with other EOC histological subtypes at P ≤ 10-7. The genotyped SNP, rs711830 (MAF=0.32), is located at 2q31.1 in the 3′ region of the HOXD3 homeobox gene (all mEOC OR=1.30, P=7.5 x 10-12 and invasive-only OR=1.26, P=7.8 x 10-7). This SNP is highly correlated (r2=0.99) with rs2072590, which is located 5.3kb upstream and previously reported by us to be associated with mEOC and invasive serous EOC (Nat Genet 2010:874). The HOXD3 locus is highly conserved and may also be a marker of Müllerian cell type origin. Conclusions The 3 new loci for mEOC require confirmation by genotyping. Fine mapping and functional analysis of all loci is needed in order to determine their likely target genes. Citation Format: Linda E. Kelemen, Jonathan Tyrer, Catherine M. Phelan, Susan J. Ramus, Andrew Berchuck, Simon A. Gayther, Ellen L. Goode, Celeste L. Pearce, Joellen M. Schildkraut, Georiga Chenevix-Trench, Alvaro N. Monteiro, Marc T. Goodman, Thomas A. Sellers, Paul PD Pharoah. GWAS identifies risk variants for mucinous ovarian carcinoma. [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 3283. doi:10.1158/1538-7445.AM2014-3283

Abstract Background:BRCA1 and BRCA2 mutation carriers are at substantially increased risk for developing breast and ovarian cancer. The incomplete penetrance coupled with the variable age at diagnosis in carriers of the same mutation suggests the existence of genetic and nongenetic modifying factors. In this study, we evaluated the putative role of variants in many candidate modifier genes. Methods: Genotyping data from 15,252 BRCA1 and 8,211 BRCA2 mutation carriers, for known variants (n = 3,248) located within or around 445 candidate genes, were available through the iCOGS custom-designed array. Breast and ovarian cancer association analysis was performed within a retrospective cohort approach. Results: The observed P values of association ranged between 0.005 and 1.000. None of the variants was significantly associated with breast or ovarian cancer risk in either BRCA1 or BRCA2 mutation carriers, after multiple testing adjustments. Conclusion: There is little evidence that any of the evaluated candidate variants act as modifiers of breast and/or ovarian cancer risk in BRCA1 or BRCA2 mutation carriers. Impact: Genome-wide association studies have been more successful at identifying genetic modifiers of BRCA1/2 penetrance than candidate gene studies. Cancer Epidemiol Biomarkers Prev; 24(1); 308–16. ©2014 AACR.