Ramus, Susan J.

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: Many epithelial ovarian cancer (EOC) risk factors relate to hormone exposure and elevated estrogen levels are associated with obesity in postmenopausal women. Therefore, we hypothesized that gene–environment interactions related to hormone-related risk factors could differ between obese and non-obese women. Methods: We considered interactions between 11,441 SNPs within 80 candidate genes related to hormone biosynthesis and metabolism and insulin-like growth factors with six hormone-related factors (oral contraceptive use, parity, endometriosis, tubal ligation, hormone replacement therapy, and estrogen use) and assessed whether these interactions differed between obese and non-obese women. Interactions were assessed using logistic regression models and data from 14 case–control studies (6,247 cases; 10,379 controls). Histotype-specific analyses were also completed. Results: SNPs in the following candidate genes showed notable interaction: IGF1R (rs41497346, estrogen plus progesterone hormone therapy, histology = all, P = 4.9 × 10−6) and ESR1 (rs12661437, endometriosis, histology = all, P = 1.5 × 10−5). The most notable obesity–gene–hormone risk factor interaction was within INSR (rs113759408, parity, histology = endometrioid, P = 8.8 × 10−6). Conclusions: We have demonstrated the feasibility of assessing multifactor interactions in large genetic epidemiology studies. Follow-up studies are necessary to assess the robustness of our findings for ESR1, CYP11A1, IGF1R, CYP11B1, INSR, and IGFBP2. Future work is needed to develop powerful statistical methods able to detect these complex interactions. Impact: Assessment of multifactor interaction is feasible, and, here, suggests that the relationship between genetic variants within candidate genes and hormone-related risk factors may vary EOC susceptibility. Cancer Epidemiol Biomarkers Prev; 25(5); 780–90. ©2016 AACR.

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 The known genetic variance for epithelial ovarian cancer (EOC) accounts for less than half the familial risks of the disease. BRCA1 and BRCA2 are the major high-penetrance EOC susceptibility genes. More recently genes that interact with BRCA1 and BRCA2 in the same biological pathways, but conferring more moderate penetrance (e.g. RAD51C, RAD51D) have also been identified. Mutations in these genes result in defects in homologous recombination and double strand break repair (HR/DSB). We hypothesized that other genes operating in these and other DNA repair pathways may also confer EOC susceptibility. Whole exome sequencing data, available from more than 400 high-grade serous ovarian cancers, were used to identify 24 candidate genes harboring germline truncating mutations in at least 2 cases, followed by replication analysis using targeted resequencing in ~3,000 cases and 3,000 controls. Subjects were also sequenced for the known susceptibility genes listed above. The data confirmed the previously reported contributions of BRCA1, BRCA2, RAD51C and RAD51D to EOC. In addition, we identified BRIP1, and FANCM as candidate EOC susceptibility genes. Surprisingly, for several genes we found truncating mutations were equally prevalent in cases as they were in controls. These mutations may still represent rare low-penetrance susceptibility alleles for EOC, which we are too limited in sample size and power to confirm; or more likely, it indicates that cells can tolerate heterozygous loss-of-function mutations in apparently critical DNA repair genes without any impact on neoplastic development. These findings may have implications for risk prediction and prevention strategies for EOC and possibly for the identification of novel therapeutic targets. Citation Format: Susan J. Ramus, Honglin Song, Ed Dicks, Maria Intermaggio, Jane Hayward, Jonathan P. Tyrer, Jacek Gronwald, Susanne Krüger Kjær, Alice Wittemore, Beth Karlan, Kunle Odunsi, Usha Menon, David D. Bowtell, Ellen L. Goode, Paul D.P. Pharoah, Simon A. Gayther. Whole exome and targeted resequencing, of population based ovarian cancer cases and controls, identifies susceptibility genes for ovarian cancer [abstract]. In: Proceedings of the 10th Biennial Ovarian Cancer Research Symposium; Sep 8-9, 2014; Seattle, WA. Philadelphia (PA): AACR; Clin Cancer Res 2015;21(16 Suppl):Abstract nr POSTER-CTRL-1213.

Abstract Background: Epithelial ovarian cancer (EOC) is composed of five major histologic types: 1) high-grade serous carcinoma (HGSC), accounting for most cases (∼70%); and the rarer 2) clear cell, 3) endometrioid, 4) mucinous, and 5) low- grade serous carcinoma (LGSC). As EOC risk factors are histology specific, as are the site of precursor lesions, accurate subtyping is critical to understanding EOC prognostic factors and etiology. Our aim was to compare two currently employed histology assignment strategies using gene expression, DNA methylation, and clinical outcome data Methods: Histology assignment based on: a) pathologist, and b) an integrated pathologist and IHC prediction algorithm (pathIHC) (using ARID1A, CDKN2A, DKK1, HNF1B, MDM2, PGR, TP53, TFF3, VIM, and WT1 staining patterns), was compared for tumors of all histologies at the Mayo Clinic, Rochester. The 500 most variable probes from Illumina Methylation450 BeadChips (N = 259) and Agilent 4×44K expression arrays (N = 245) were used to perform unsupervised hierarchical clustering. Fisher's Exact test was used to test the association of clusters with histology. Cox proportional hazards regression analysis was used to test the association of clusters and histology with time to progression (TTP) and time to death (TTD). Results: Eighty-two percent of tumors were concordant between histology assignment strategies. Clustering based on methylation data produced two distinct clusters. PathIHC produced more homogenous clusters (p-value = 2.8×10-16) than pathology alone (p = 1.3×10-9). Cluster M1, characterized by high levels of methylation across almost all probes, was enriched for the rarer EOC subtypes; cluster M2, characterized by moderate levels of methylation across 50%-75% of probes, was predominantly HGSC. These patterns were observed irrespective of histology assignment strategy; however, when using histology by pathologist, cluster M2 had more endometrioid and LGSC tumors interspersed with HGSC tumors. Clustering based on expression data also produced two clusters. PathIHC produced slightly more homogenous clusters (p = 9.0×10-11, pathology alone, p = 2.3×10-10). Cluster E1 was enriched for the rarer EOC subtypes; cluster E2 was primarily serous lineage tumours (HGSC and LGSC), particularly when using histology by pathIHC. Considering only clinical outcomes, histologies by pathology were more significantly different in terms of TTP (p = 1.1×10-6) and TTD (p = 1.7×10-3) than histologies by pathIHC (TTP, p = 5.4×10-6; TTD, p = 5.7×10-3). Conclusions: Histology by pathIHC produced more homogeneous clusters in both the methylation and expression data; thus, molecular data supports this strategy. Differences in clinical outcomes (TTP and TTD) between histology groups were more pronounced when using assignment by pathologist; thus, clinical data supports this strategy. Citation Format: M A. Earp, S J. Winham, S M. Armasu, B L. Fridley, M C. Larson, Z C. Fogarty, K R. Kalli, C Wang, G L. Keeney, J M. Cunningham, S Ramus, M Kobel, E L. Goode. Comparison of pathology versus IHC-based ovarian carcinoma histology assignment using gene expression, DNA methylation, and clinical outcome data. [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 4246. doi:10.1158/1538-7445.AM2015-4246

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