Phelan, Catherine

Thymidylate synthase (TYMS) is a crucial enzyme for DNA synthesis. TYMS expression is regulated by its antisense mRNA, ENOSF1. Disrupted regulation may promote uncontrolled DNA synthesis and tumor growth. We sought to replicate our previously reported association between rs495139 in the TYMS-ENOSF1 3′ gene region and increased risk of mucinous ovarian carcinoma (MOC) in an independent sample. Genotypes from 24,351 controls to 15,000 women with invasive OC, including 665 MOC, were available. We estimated per-allele odds ratios (OR) and 95% confidence intervals (CI) using unconditional logistic regression, and meta-analysis when combining these data with our previous report. The association between rs495139 and MOC was not significant in the independent sample (OR = 1.09; 95% CI = 0.97–1.22; p = 0.15; N = 665 cases). Meta-analysis suggested a weak association (OR = 1.13; 95% CI = 1.03–1.24; p = 0.01; N = 1019 cases). No significant association with risk of other OC histologic types was observed (p = 0.05 for tumor heterogeneity). In expression quantitative trait locus (eQTL) analysis, the rs495139 allele was positively associated with ENOSF1 mRNA expression in normal tissues of the gastrointestinal system, particularly esophageal mucosa (r = 0.51, p = 1.7 × 10−28), and nonsignificantly in five MOC tumors. The association results, along with inconclusive tumor eQTL findings, suggest that a true effect of rs495139 might be small.

Telomeres cap chromosome ends, protecting them from degradation, double‐strand breaks, and end‐to‐end fusions. Telomeres are maintained by telomerase, a reverse transcriptase encoded by TERT, and an RNA template encoded by TERC. Loci in the TERT and adjoining CLPTM1L region are associated with risk of multiple cancers. We therefore investigated associations between variants in 22 telomere structure and maintenance gene regions and colorectal, breast, prostate, ovarian, and lung cancer risk. We performed subset‐based meta‐analyses of 204,993 directly‐measured and imputed SNPs among 61,851 cancer cases and 74,457 controls of European descent. Independent associations for SNP minor alleles were identified using sequential conditional analysis (with gene‐level p value cutoffs ≤3.08 × 10−5). Of the thirteen independent SNPs observed to be associated with cancer risk, novel findings were observed for seven loci. Across the DCLRE1B region, rs974494 and rs12144215 were inversely associated with prostate and lung cancers, and colorectal, breast, and prostate cancers, respectively. Across the TERC region, rs75316749 was positively associated with colorectal, breast, ovarian, and lung cancers. Across the DCLRE1B region, rs974404 and rs12144215 were inversely associated with prostate and lung cancers, and colorectal, breast, and prostate cancers, respectively. Near POT1, rs116895242 was inversely associated with colorectal, ovarian, and lung cancers, and RTEL1 rs34978822 was inversely associated with prostate and lung cancers. The complex association patterns in telomere‐related genes across cancer types may provide insight into mechanisms through which telomere dysfunction in different tissues influences cancer risk.

Abstract Breast, ovarian, and prostate cancers are hormone-related and may have a shared genetic basis, but this has not been investigated systematically by genome-wide association (GWA) studies. Meta-analyses combining the largest GWA meta-analysis data sets for these cancers totaling 112,349 cases and 116,421 controls of European ancestry, all together and in pairs, identified at P < 10−8 seven new cross-cancer loci: three associated with susceptibility to all three cancers (rs17041869/2q13/BCL2L11; rs7937840/11q12/INCENP; rs1469713/19p13/GATAD2A), two breast and ovarian cancer risk loci (rs200182588/9q31/SMC2; rs8037137/15q26/RCCD1), and two breast and prostate cancer risk loci (rs5013329/1p34/NSUN4; rs9375701/6q23/L3MBTL3). Index variants in five additional regions previously associated with only one cancer also showed clear association with a second cancer type. Cell-type–specific expression quantitative trait locus and enhancer–gene interaction annotations suggested target genes with potential cross-cancer roles at the new loci. Pathway analysis revealed significant enrichment of death receptor signaling genes near loci with P < 10−5 in the three-cancer meta-analysis. Significance: We demonstrate that combining large-scale GWA meta-analysis findings across cancer types can identify completely new risk loci common to breast, ovarian, and prostate cancers. We show that the identification of such cross-cancer risk loci has the potential to shed new light on the shared biology underlying these hormone-related cancers. Cancer Discov; 6(9); 1052–67. ©2016 AACR. This article is highlighted in the In This Issue feature, p. 932

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<10−3) unique to each sub-type. Gene-level analyses based on AML and SKAT-CR revealed that ADAR1 was globally associated (P<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<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<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<0.001) with the top SNP located within the 5′UTR of HOXD8 transcription factor (rs847163 (T>C); p<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<0.001 and 16 SNPs with p<0.01, all were inversely associated with risk and consistent with the top hit at rs927901 T>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<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 (<v:shape style = “WIDTH: 16.5pt; HEIGHT: 18pt” id = _x0000_i1026 equationxml = ‘14pχ2′ type = “#_x0000_t75”> = 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<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 > 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 Introduction The epithelial-mesenchymal transition (EMT) pathway contributes to epithelial ovarian cancer (EOC) progression. Previously we identified associations of single nucleotide polymorphisms (SNPs) in EMT-related genes and EOC risk (data not shown). To further investigate the role of EMT-related gene variants in EOC and to derive more reliable risk estimates, we evaluated associations for 795 SNPs from 278 EMT-related genes in a replication study of 43 studies within the Ovarian Cancer Association Consortium (OCAC). Methods The study population included 14,736 cases and 23,448 controls of European ancestry, 89 cases and 200 controls of African ancestry, and 249 cases and 1574 controls of Asian ancestry. The 795 SNPs were genotyped using an Illumina Infinium iSelect BeadChip as part of the Collaborative Oncological Gene-environment Study (COGS). For women of European ancestry, both invasive cancers combined and the four main histological subtypes (serous [n=8,372], endometroid [n=2,068], clear cell [n=1,025] and mucinous [n=943] were analyzed, while for women of African and Asian ancestries only the serous subtype was analyzed. SNP analyses were conducted using unconditional logistic regression under a log-additive model separately for women of each ancestry. All analyses were adjusted for study site and population substructure within each ancestry. Results For women of European ancestry, the strongest evidence of an association for invasive cancers combined was observed for IGF1R rs10794486 (OR=1.05, 95%CI=1.02-1.09, P=0.002), which was also associated with serous (P=0.019), but not the other histological subtypes. In the histological subtype-specific analysis, the most significant association for serous was observed at NRP2 rs3771044 (OR=1.09, 95%CI=1.03-1.15, P=0.00014), for endometroid at rs1770247 (13q31.3, OR=1.16, 95%CI=1.07-1.25, P=0.00028), for clear cell at rs4848300 (2q13, OR=0.84, 95%CI=0.76-0.93, P=0.00099) and for mucinous at SEMA4B rs8030039 (OR=1.09, 95%CI=1.01-1.17, P=0.00019). For women of African ancestry, the strongest association for serous histological subtype was observed at F9 rs6048 (OR=1.71, 95%CI=1.17-2.50, P=0.006). This SNP was not significant in women of European or Asian ancestry. The strongest association for women of Asian ancestry was IGF1R rs10794486 (OR=1.68, 95%CI=1.16-2.45, P=0.007), which was also significant in women of European, but not African, ancestry. Conclusion Findings from this large study provide additional evidence that variants in EMT-related genes may be associated with ovarian cancer risk and that EMT susceptibility loci may differ by ovarian cancer histological subtypes. Furthermore, the findings suggest that susceptibility locus for the serous subtype may be different for women of African ancestry. Future studies are warranted to confirm these findings. Citation Format: Ernest K. Amankwah, Jonathan Tyrer, Hui-Yi Lin, Ya-Yu Tsai, Zhihua Chen, Gang Han, Xiaotao Qu, Ellen Goode, Julie Cunninghan, Edward Iverson, Susan Ramus, Andrew Berchuck, Joellen Schildkraut, Alvaro Monteiro, Simon Gayther, Steven Narod, Paul Pharoah, Thomas A. Sellers, Catherine Phelan. Epithelial-mesenchymal transition (EMT) gene variants influence epithelial ovarian cancer risk in women of European, African and Asian ancestry. [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 1352. doi:10.1158/1538-7445.AM2013-1352

Abstract As the leading cause of mortality among gynecologic cancers, it is critical to increase understanding of ovarian cancer pathogenesis, which remains largely unknown. Several xenobiotic metabolizing enzymes are responsible for processing and inactivating carcinogens that could initiate ovarian carcinogenesis. Here, we hypothesize that single-nucleotide polymorphisms (SNPs) in genes encoding xenobiotic metabolizing enzymes may be associated with differential risk of ovarian cancer. Cases with epithelial ovarian cancer (N=930) and controls (N=1,037), frequency-matched on age and geographic residence, were enrolled at Mayo Clinic and Duke University, and genotyped at 163 tagSNPs in EPHX1, SULT1E1, ADH5, ADH4, ADH6, ADH1A, ADH1B, ADH1C, NQO2, NAT1, NAT2, GSTP1, CYP19A1, CYP1A1, CYP1A2, and NQO1 using an Illumina GoldenGate™ BeadArray array. Logistic regression was used to estimate odds ratios (ORs) and adjust for population structure and other potential confounders, including oral contraceptive use, parity, and family history. SNP-specific statistical analyses identified several noteworthy results including an association with increased risk of ovarian cancer with a SNP (rs2917666) in NAD(P)H dehydrogenase, quinone 1 (NQO1) (per-allele OR 1.16, 95% confidence interval (CI) 1.01-1.33, p=0.04). The SNP rs1051740 in epoxide hydrolase 1, microsomal (xenobiotic) (EPHX1) was suggestive of increased risk of invasive serous disease (TT v CT OR 0.85, 95% CI 0.66-1.09, TT v CC 1.29, 95% CI 0.89-1.88, p=0.09). In an effort to replicate these observations, a subset of 11 SNPs with a p-value < 0.10 were then genotyped in 727 invasive cancer cases and 1,065 controls from an Australian ovarian cancer study using a Sequenom iPlex assay. Results of combined analyses confirmed an association between EPHX1 rs1051740 and increased invasive serous ovarian cancer risk (TT v CT OR 0.99, 95% CI 0.77-1.29, TT v CC OR 1.61, 95% CI 1.23-2.11, p=0.001; recessive OR 1.61, 95% CI 1.25-2.09, p=0.001). This gene is involved in the activation and detoxification of tobacco-derived carcinogens, and smoking has been associated with the mucinous histologic subtype of ovarian cancer. Analyses by histological subtype and gene-environment interaction analyses considering genotype in conjunction with tobacco use and alcohol intake are underway. Due to the important role of EPHX1 in metabolizing pro-carcinogens and the observed association between rs1051740 and risk of invasive serous disease, an examination of this SNP in a larger ovarian cancer consortium is warranted. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 914.

Abstract The p53 protein is critical for multiple cellular functions including cell growth and DNA repair. We assessed whether polymorphisms in the region encoding TP53 were associated with risk of invasive ovarian cancer. The study population includes a total of 5,206 invasive ovarian cancer cases (2,829 of which were serous) and 8,790 controls from 13 case-control or nested case-control studies participating in the Ovarian Cancer Association Consortium (OCAC). Three of the studies performed independent discovery investigations involving genotyping of up to 23 single nucleotide polymorphisms (SNP) in the TP53 region. Significant findings from this discovery phase were followed up for replication in the other OCAC studies. Mixed effects logistic regression was used to generate posterior median per allele odds ratios (OR), 95% probability intervals (PI), and Bayes factors (BF) for genotype associations. Five SNPs showed significant associations with risk in one or more of the discovery investigations and were followed up by OCAC. Mixed effects analysis confirmed associations with serous invasive cancers for two correlated (r2 = 0.62) SNPs: rs2287498 (median per allele OR, 1.30; 95% PI, 1.07–1.57) and rs12951053 (median per allele OR, 1.19; 95% PI, 1.01–1.38). Analyses of other histologic subtypes suggested similar associations with endometrioid but not with mucinous or clear cell cancers. This large study provides statistical evidence for a small increase in risk of ovarian cancer associated with common variants in the TP53 region. [Cancer Res 2009;69(6):2349–57]

Abstract Studies suggest that underglycosylation of the cell membrane mucin MUC1 may be associated with epithelial ovarian cancer. We identified 26 genes involved in glycosylation and examined 93 single nucleotide polymorphisms (SNP) with a minor allele frequency of ≥0.05 in relation to incident ovarian cancer. Cases were ascertained at the Mayo Clinic, Rochester, MN (n = 396) or a 48-county region in North Carolina (Duke University; n = 534). Ovarian cancer–free controls (n = 1,037) were frequency matched to the cases on age, race, and residence. Subjects were interviewed to obtain data on risk factors and a sample of blood for DNA and genotyped using the Illumina GoldenGate assay. We excluded subjects and individual SNPs with genotype call rates of <90%. Data were analyzed using logistic regression, with adjustment for age and residence. We fitted dominant, log additive, and recessive genetic models. Among Caucasians, nine SNPs in eight genes were associated with risk at P < 0.05 under at least one genetic model before adjusting for multiple testing. A SNP in GALNT1 (rs17647532) was the only one that remained statistically significant after Bonferroni adjustment for multiple testing but was not statistically significant in Hardy-Weinberg equilibrium among controls. Haplotype analyses revealed a global association of GALNT1 with risk (P = 0.038, under a recessive genetic model), which largely reflected a decreased risk of one haplotype (0.10 frequency; odds ratio, 0.07; P = 0.01) compared with the most common haplotype (0.39 frequency). These results suggest that genetic polymorphisms in the glycoslyation process may be novel risk factors for ovarian cancer. (Cancer Epidemiol Biomarkers Prev 2008;17(2):397–404)