Anton‐Culver, Hoda

Abstract Background: A recent association study identified a common variant (rs9790517) at 4q24 to be associated with breast cancer risk. Independent association signals and potential functional variants in this locus have not been explored. Methods: We conducted a fine-mapping analysis in 55,540 breast cancer cases and 51,168 controls from the Breast Cancer Association Consortium. Results: Conditional analyses identified two independent association signals among women of European ancestry, represented by rs9790517 [conditional P = 2.51 × 10−4; OR, 1.04; 95% confidence interval (CI), 1.02–1.07] and rs77928427 (P = 1.86 × 10−4; OR, 1.04; 95% CI, 1.02–1.07). Functional annotation using data from the Encyclopedia of DNA Elements (ENCODE) project revealed two putative functional variants, rs62331150 and rs73838678 in linkage disequilibrium (LD) with rs9790517 (r2 ≥ 0.90) residing in the active promoter or enhancer, respectively, of the nearest gene, TET2. Both variants are located in DNase I hypersensitivity and transcription factor–binding sites. Using data from both The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC), we showed that rs62331150 was associated with level of expression of TET2 in breast normal and tumor tissue. Conclusion: Our study identified two independent association signals at 4q24 in relation to breast cancer risk and suggested that observed association in this locus may be mediated through the regulation of TET2. Impact: Fine-mapping study with large sample size warranted for identification of independent loci for breast cancer risk. Cancer Epidemiol Biomarkers Prev; 24(11); 1680–91. ©2015 AACR.

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 < 0.05 and FDR < 0.05). These results were replicated (P < 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 Recent studies have shown that associations between common genetic variation and gene expression in breast tumors provide insight into the functional mechanisms underlying breast cancer risk loci. However, the relationship between genetic variation, gene expression, and risk of triple negative (TN) breast cancer remains largely unexplored. We performed a genome-wide expression quantitative trait locus (eQTL) analysis of TN breast cancer using 668 formalin-fixed paraffin embedded TN tumors from the Triple Negative Breast Cancer Consortium (TNBCC) to identify novel genes relevant to TN breast cancer risk and to further explore the biology underlying known TN risk loci. Cis-eQTLs were defined as correlations between single nucleotide polymorphisms (SNP) and expression of genes within 1Mb, and trans-eQTLs were defined as the remainder of all SNP-gene correlations. In a genome-wide analysis, we identified 68,012 cis-eQTLs, representing correlations between 42,193 SNPs and 2,092 genes that were significant at a 10% false discovery rate (FDR). We also identified 70 trans-eQTLs across two unique genes at the same FDR threshold. Five cis-eQTLs involved SNPs previously associated with TN breast cancer risk (p<5.0x10-5) in 3,700 TN cases and 4,700 controls from the TNBCC. Three of these, located on chromosomes 5p (cis-eQTL p=4.8x10-34; risk p=1.1x10-7), 7p (cis-eQTL p=6.1x10-10; risk p=1.3x10-7), and 17q (cis-eQTL p=1.8x10-15; risk p=2.9x10-7), have previously been reported in public databases of normal tissue. The remaining two TN breast cancer cis-eQTLs, located on chromosomes 8q (cis-eQTL p=6.3x10-6, risk p=8.3x10-6) and 14q (cis-eQTL p=3.1x10-6, risk p=1.2x10-5), have not been reported as eQTLs in other tissues. In addition, 12 of 25 known TN breast cancer risk loci (PEX14, MDM4, 2q31.1, ESR1, 11q13.1, 11q24.3, 12p13.1, PTHLH, NTN4, 12q24, 19p13.1, MLK1), contained at least one significant cis-eQTL. In summary, we have identified five cis-eQTLs that may identify novel risk loci for TN breast cancer. Citation Format: Kristen S. Purrington, Drakoulis Yannoukakos, Jane Carpenter, Heli Nevanlinna, Angela Cox, Gianluca Severi, Christine Ambrosone, Amanda Ewart Toland, Andrew K. Godwin, Hiltrud Brauch, Peter A. Fasching, Penelope Miron, Jenny Chang-Claude, Nicholas G. Martin, Grant W. Montgomery, Vessela Kristensen, Hoda Anton-Culver, Paul Goodfellow, Janet E. Olson, Hugues Sicotte, Naresh Prodduturi, Daniel W. Visscher, Jeanette E. Eckel-Passow, S. Keith Anderson, Seth Slettedahl, Curtis Olswold, Xianshu Wang, V. Shane Pankratz, Susan Slager, Wei Zheng, Arto Mannermaa, Ute Hamann, Diana M. Eccles, Celine M. Vachon, Fergus J. Couch. Expression quantitative trait locus analysis of triple negative breast 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 3266. doi:10.1158/1538-7445.AM2014-3266

ScopeWe reevaluated previously reported associations between variants in pathways of one‐carbon (1‐C) (folate) transfer genes and ovarian carcinoma (OC) risk, and in related pathways of purine and pyrimidine metabolism, and assessed interactions with folate intake.Methods and resultsOdds ratios (OR) for 446 genetic variants were estimated among 13 410 OC cases and 22 635 controls, and among 2281 cases and 3444 controls with folate information. Following multiple testing correction, the most significant main effect associations were for dihydropyrimidine dehydrogenase (DPYD) variants rs11587873 (OR = 0.92; p = 6 × 10−5) and rs828054 (OR = 1.06; p = 1 × 10−4). Thirteen variants in the pyrimidine metabolism genes, DPYD, DPYS, PPAT, and TYMS, also interacted significantly with folate in a multivariant analysis (corrected p = 9.9 × 10−6) but collectively explained only 0.2% of OC risk. Although no other associations were significant after multiple testing correction, variants in SHMT1 in 1‐C transfer, previously reported with OC, suggested lower risk at higher folate (pinteraction = 0.03–0.006).ConclusionVariation in pyrimidine metabolism genes, particularly DPYD, which was previously reported to be associated with OC, may influence risk; however, stratification by folate intake is unlikely to modify disease risk appreciably in these women. SHMT1 SNP‐by‐folate interactions are plausible but require further validation. Polymorphisms in selected genes in purine metabolism were not associated with OC.

Abstract Part of the substantial unexplained familial aggregation of breast cancer may be due to interactions between common variants, but few studies have had adequate statistical power to detect interactions of realistic magnitude. We aimed to assess all two-way interactions in breast cancer susceptibility between 70 917 single nucleotide polymorphisms (SNPs) selected primarily based on prior evidence of a marginal effect. Thirty-eight international studies contributed data for 46 450 breast cancer cases and 42 461 controls of European origin as part of a multi-consortium project (COGS). First, SNPs were preselected based on evidence (P < 0.01) of a per-allele main effect, and all two-way combinations of those were evaluated by a per-allele (1 d.f.) test for interaction using logistic regression. Second, all 2.5 billion possible two-SNP combinations were evaluated using Boolean operation-based screening and testing, and SNP pairs with the strongest evidence of interaction (P < 10−4) were selected for more careful assessment by logistic regression. Under the first approach, 3277 SNPs were preselected, but an evaluation of all possible two-SNP combinations (1 d.f.) identified no interactions at P < 10−8. Results from the second analytic approach were consistent with those from the first (P > 10−10). In summary, we observed little evidence of two-way SNP interactions in breast cancer susceptibility, despite the large number of SNPs with potential marginal effects considered and the very large sample size. This finding may have important implications for risk prediction, simplifying the modelling required. Further comprehensive, large-scale genome-wide interaction studies may identify novel interacting loci if the inherent logistic and computational challenges can be overcome.

Abstract Background: There are several well-established environmental risk factors for ovarian cancer, and recent genome-wide association studies have also identified six variants that influence disease risk. However, the interplay between such risk factors and susceptibility loci has not been studied. Methods: Data from 14 ovarian cancer case–control studies were pooled, and stratified analyses by each environmental risk factor with tests for heterogeneity were conducted to determine the presence of interactions for all histologic subtypes. A genetic “risk score” was created to consider the effects of all six variants simultaneously. A multivariate model was fit to examine the association between all environmental risk factors and genetic risk score on ovarian cancer risk. Results: Among 7,374 controls and 5,566 cases, there was no statistical evidence of interaction between the six SNPs or genetic risk score and the environmental risk factors on ovarian cancer risk. In a main effects model, women in the highest genetic risk score quartile had a 65% increased risk of ovarian cancer compared with women in the lowest [95% confidence interval (CI), 1.48–1.84]. Analyses by histologic subtype yielded risk differences across subtype for endometriosis (Phet < 0.001), parity (Phet < 0.01), and tubal ligation (Phet = 0.041). Conclusions: The lack of interactions suggests that a multiplicative model is the best fit for these data. Under such a model, we provide a robust estimate of the effect of each risk factor that sets the stage for absolute risk prediction modeling that considers both environmental and genetic risk factors. Further research into the observed differences in risk across histologic subtype is warranted. Cancer Epidemiol Biomarkers Prev; 22(5); 880–90. ©2013 AACR.