Citations

AbstractPrevious transcriptome‐wide association studies (TWAS) have identified breast cancer risk genes by integrating data from expression quantitative loci and genome‐wide association studies (GWAS), but analyses of breast cancer subtype‐specific associations have been limited. In this study, we conducted a TWAS using gene expression data from GTEx and summary statistics from the hitherto largest GWAS meta‐analysis conducted for breast cancer overall, and by estrogen receptor subtypes (ER+ and ER−). We further compared associations with ER+ and ER− subtypes, using a case‐only TWAS approach. We also conducted multigene conditional analyses in regions with multiple TWAS associations. Two genes, STXBP4 and HIST2H2BA, were specifically associated with ER+ but not with ER– breast cancer. We further identified 30 TWAS‐significant genes associated with overall breast cancer risk, including four that were not identified in previous studies. Conditional analyses identified single independent breast‐cancer gene in three of six regions harboring multiple TWAS‐significant genes. Our study provides new information on breast cancer genetics and biology, particularly about genomic differences between ER+ and ER− breast cancer.

Abstract Pathogenic sequence variants (PSV) in BRCA1 or BRCA2 (BRCA1/2) are associated with increased risk and severity of prostate cancer. We evaluated whether PSVs in BRCA1/2 were associated with risk of overall prostate cancer or high grade (Gleason 8+) prostate cancer using an international sample of 65 BRCA1 and 171 BRCA2 male PSV carriers with prostate cancer, and 3,388 BRCA1 and 2,880 BRCA2 male PSV carriers without prostate cancer. PSVs in the 3′ region of BRCA2 (c.7914+) were significantly associated with elevated risk of prostate cancer compared with reference bin c.1001-c.7913 [HR = 1.78; 95% confidence interval (CI), 1.25–2.52; P = 0.001], as well as elevated risk of Gleason 8+ prostate cancer (HR = 3.11; 95% CI, 1.63–5.95; P = 0.001). c.756-c.1000 was also associated with elevated prostate cancer risk (HR = 2.83; 95% CI, 1.71–4.68; P = 0.00004) and elevated risk of Gleason 8+ prostate cancer (HR = 4.95; 95% CI, 2.12–11.54; P = 0.0002). No genotype–phenotype associations were detected for PSVs in BRCA1. These results demonstrate that specific BRCA2 PSVs may be associated with elevated risk of developing aggressive prostate cancer. Significance: Aggressive prostate cancer risk in BRCA2 mutation carriers may vary according to the specific BRCA2 mutation inherited by the at-risk individual.

AbstractAn association between genetic variants in the vitamin D receptor (VDR) gene and epithelial ovarian cancer (EOC) was previously reported in women of African ancestry (AA). We sought to examine associations between genetic variants in VDR and additional genes from vitamin D biosynthesis and pathway targets (EGFR, UGT1A, UGT2A1/2, UGT2B, CYP3A4/5, CYP2R1, CYP27B1, CYP24A1, CYP11A1, and GC). Genotyping was performed using the custom‐designed 533,631 SNP Illumina OncoArray with imputation to the 1,000 Genomes Phase 3 v5 reference set in 755 EOC cases, including 537 high‐grade serous (HGSOC), and 1,235 controls. All subjects are of African ancestry (AA). Logistic regression was performed to estimate odds ratios (OR) and 95% confidence intervals (CI). We further evaluated statistical significance of selected SNPs using the Bayesian False Discovery Probability (BFDP). A significant association with EOC was identified in the UGT2A1/2 region for the SNP rs10017134 (per allele OR = 1.4, 95% CI = 1.2‐1.7, P = 1.2 × 10−6, BFDP = 0.02); and an association with HGSOC was identified in the EGFR region for the SNP rs114972508 (per allele OR = 2.3, 95% CI = 1.6‐3.4, P = 1.6 × 10−5, BFDP = 0.29) and in the UGT2A1/2 region again for rs1017134 (per allele OR = 1.4, 95% CI = 1.2‐1.7, P = 2.3 × 10−5, BFDP = 0.23). Genetic variants in the EGFR and UGT2A1/2 may increase susceptibility of EOC in AA women. Future studies to validate these findings are warranted. Alterations in EGFR and UGT2A1/2 could perturb enzyme efficacy, proliferation in ovaries, impact and mark susceptibility to EOC.

Abstract DNA methylation is instrumental for gene regulation. Global changes in the epigenetic landscape have been recognized as a hallmark of cancer. However, the role of DNA methylation in epithelial ovarian cancer (EOC) remains unclear. In this study, high-density genetic and DNA methylation data in white blood cells from the Framingham Heart Study (N = 1,595) were used to build genetic models to predict DNA methylation levels. These prediction models were then applied to the summary statistics of a genome-wide association study (GWAS) of ovarian cancer including 22,406 EOC cases and 40,941 controls to investigate genetically predicted DNA methylation levels in association with EOC risk. Among 62,938 CpG sites investigated, genetically predicted methylation levels at 89 CpG were significantly associated with EOC risk at a Bonferroni-corrected threshold of P < 7.94 × 10−7. Of them, 87 were located at GWAS-identified EOC susceptibility regions and two resided in a genomic region not previously reported to be associated with EOC risk. Integrative analyses of genetic, methylation, and gene expression data identified consistent directions of associations across 12 CpG, five genes, and EOC risk, suggesting that methylation at these 12 CpG may influence EOC risk by regulating expression of these five genes, namely MAPT, HOXB3, ABHD8, ARHGAP27, and SKAP1. We identified novel DNA methylation markers associated with EOC risk and propose that methylation at multiple CpG may affect EOC risk via regulation of gene expression. Significance: Identification of novel DNA methylation markers associated with EOC risk suggests that methylation at multiple CpG may affect EOC risk through regulation of gene expression.

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.

Abstract Tubo-ovarian high-grade serous cancer (HGSC) is the most common histotype of epithelial ovarian cancer (70%-80%) and the leading cause of death from gynecological malignancy. Novel monoclonal antibody and adoptive cell transfer immunotherapies have found success in cancer in the past five years; however success in HGSC treatment has been elusive thus far. We have recently shown that CD8+ T-lymphocytes infiltrating the tumor proper (TILs) associate in a dose response manner with longer patient survival time regardless of age, stage or residual disease. As understanding host factors associated with this clinically relevant tumor phenotype will contribute to the development of targeted therapies, we sought to identify germline genetic factors associated with CD8+ TIL levels in HGSC. We performed a germline genome-wide association study (GWAS) of CD8+ TILs using 1204 HGSC patients of European ancestry from twelve studies with centralized genotyping, immunohistochemistry, and harmonization of clinical data. Germline genotypes were assayed using the Illumina Infinium OncoArray Beadchip and were imputed to over 11.4 million variants based on the 1000 Genomes Project reference panel. CD8 immunohistochemistry was performed on tissue microarrays using the Leica Bond Rx stainer and scored into four levels of CD8+ TILs per x400 magnification of a 0.55-mm diameter field for each tumor hotspot: none (no CD8+ TILs), low (1-2 CD8+ TILs), moderate (3-19 CD8+ TILs), high (20 or more CD8+ TILs). GWAS used multinomial regression modelling of these four CD8+ TIL levels adjusted for study site, age, and European ancestry principal components. On the average, CD8+ TIL levels were 18%, 19%, 44% and 18% for none, low, moderate and high levels, respectively. As expected, CD8+ TIL levels associated with overall survival time (p < 2.6x10-8), and CD8+ TIL levels and age at diagnosis were similarly distributed across study sites. We identified a region on chromosome 22 with multiple, correlated variants (MAF > 20%) associated with CD8+ TIL levels at p-value < 2.0x10-5. The most statistically significant variant was an indel (rs557925408) located in the an intron of the MYO18B gene (Myosin XVIIIB, p-value < 8.7x10-6). Each copy of this variant associated with a reduced extent of CD8+ TIL infiltration (OR low v none = 0.92; OR moderate v none = 0.61). An independent set of 1200 OncoArray genotyped HGSC cases are undergoing CD8+ TIL scoring in December 2017 and will be included in final genome-wide analysis. This will provide improved statistical power and enable consideration of additional genetic models (e.g., linear trend tests). Comprehensive integrative analysis of germline, tumor, and clinical features provides a model for clinical molecular epidemiology studies and will be key to increasing our understanding of this important HGSC immunophenotype. Citation Format: Yanina Natanzon, Martin Köbel, Stacey J. Winham, Sebastian M. Armasu, Bryan M. McCauley, Robert A. Vierkant, Julie M. Cunningham, David Bowtell, Ian G. Campbell, Jenny Chang-Claude, Anna deFazio, Peter A. Fasching, Mark T. Goodman, Beth Y. Karlan, Francesmary Modugno, Kirsten B. Moysich, Roberta B. Ness, Weiva Sieh, Paul D. Pharoah, Susan J. Ramus, Ellen L. Goode. Tumor-infiltrating CD8-positive T-lymphocytes in tubo-ovarian high-grade serous cancer are associated with multiple germline variants in 22q12.1 in a genome-wide association analysis [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 225.

Abstract Large-scale genome-wide association studies (GWAS) have identified approximately 35 loci associated with epithelial ovarian cancer (EOC) risk. The majority of GWAS-identified disease susceptibility variants are located in noncoding regions, and causal genes underlying these associations remain largely unknown. Here, we performed a transcriptome-wide association study to search for novel genetic loci and plausible causal genes at known GWAS loci. We used RNA sequencing data (68 normal ovarian tissue samples from 68 individuals and 6,124 cross-tissue samples from 369 individuals) and high-density genotyping data from European descendants of the Genotype-Tissue Expression (GTEx V6) project to build ovarian and cross-tissue models of genetically regulated expression using elastic net methods. We evaluated 17,121 genes for their cis-predicted gene expression in relation to EOC risk using summary statistics data from GWAS of 97,898 women, including 29,396 EOC cases. With a Bonferroni-corrected significance level of P < 2.2 × 10−6, we identified 35 genes, including FZD4 at 11q14.2 (Z = 5.08, P = 3.83 × 10−7, the cross-tissue model; 1 Mb away from any GWAS-identified EOC risk variant), a potential novel locus for EOC risk. All other 34 significantly associated genes were located within 1 Mb of known GWAS-identified loci, including 23 genes at 6 loci not previously linked to EOC risk. Upon conditioning on nearby known EOC GWAS-identified variants, the associations for 31 genes disappeared and three genes remained (P < 1.47 × 10−3). These data identify one novel locus (FZD4) and 34 genes at 13 known EOC risk loci associated with EOC risk, providing new insights into EOC carcinogenesis. Significance: Transcriptomic analysis of a large cohort confirms earlier GWAS loci and reveals FZD4 as a novel locus associated with EOC risk. Cancer Res; 78(18); 5419–30. ©2018 AACR.