Pharoah, Paul D.P.

PURPOSETo estimate age-specific relative and absolute cancer risks of breast cancer and to estimate risks of ovarian, pancreatic, male breast, prostate, and colorectal cancers associated with germline PALB2 pathogenic variants (PVs) because these risks have not been extensively characterized.METHODSWe analyzed data from 524 families with PALB2 PVs from 21 countries. Complex segregation analysis was used to estimate relative risks (RRs; relative to country-specific population incidences) and absolute risks of cancers. The models allowed for residual familial aggregation of breast and ovarian cancer and were adjusted for the family-specific ascertainment schemes.RESULTSWe found associations between PALB2 PVs and risk of female breast cancer (RR, 7.18; 95% CI, 5.82 to 8.85; P = 6.5 × 10−76), ovarian cancer (RR, 2.91; 95% CI, 1.40 to 6.04; P = 4.1 × 10−3), pancreatic cancer (RR, 2.37; 95% CI, 1.24 to 4.50; P = 8.7 × 10−3), and male breast cancer (RR, 7.34; 95% CI, 1.28 to 42.18; P = 2.6 × 10−2). There was no evidence for increased risks of prostate or colorectal cancer. The breast cancer RRs declined with age ( P for trend = 2.0 × 10−3). After adjusting for family ascertainment, breast cancer risk estimates on the basis of multiple case families were similar to the estimates from families ascertained through population-based studies ( P for difference = .41). On the basis of the combined data, the estimated risks to age 80 years were 53% (95% CI, 44% to 63%) for female breast cancer, 5% (95% CI, 2% to 10%) for ovarian cancer, 2%-3% (95% CI females, 1% to 4%; 95% CI males, 2% to 5%) for pancreatic cancer, and 1% (95% CI, 0.2% to 5%) for male breast cancer.CONCLUSIONThese results confirm PALB2 as a major breast cancer susceptibility gene and establish substantial associations between germline PALB2 PVs and ovarian, pancreatic, and male breast cancers. These findings will facilitate incorporation of PALB2 into risk prediction models and optimize the clinical cancer risk management of PALB2 PV carriers.

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.

AbstractFemale Hodgkin lymphoma (HL) patients treated with chest radiotherapy (RT) have a very high risk of breast cancer. The contribution of genetic factors to this risk is unclear. We therefore examined 211 155 germline single-nucleotide polymorphisms (SNPs) for gene-radiation interaction on breast cancer risk in a case-only analysis including 327 breast cancer patients after chest RT for HL and 4671 first primary breast cancer patients. Nine SNPs showed statistically significant interaction with RT on breast cancer risk (false discovery rate, <20%), of which 1 SNP in the PVT1 oncogene attained the Bonferroni threshold for statistical significance. A polygenic risk score (PRS) composed of these SNPs (RT-interaction-PRS) and a previously published breast cancer PRS (BC-PRS) derived in the general population were evaluated in a case-control analysis comprising the 327 chest-irradiated HL patients with breast cancer and 491 chest-irradiated HL patients without breast cancer. Patients in the highest tertile of the RT-interaction-PRS had a 1.6-fold higher breast cancer risk than those in the lowest tertile. Remarkably, we observed a fourfold increased RT-induced breast cancer risk in the highest compared with the lowest decile of the BC-PRS. On a continuous scale, breast cancer risk increased 1.4-fold per standard deviation of the BC-PRS, similar to the effect size found in the general population. This study demonstrates that genetic factors influence breast cancer risk after chest RT for HL. Given the high absolute breast cancer risk in radiation-exposed women, these results can have important implications for the management of current HL survivors and future patients.

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.

As a follow‐up to genome‐wide association analysis of common variants associated with ovarian carcinoma (cancer), our study considers seven well‐known ovarian cancer risk factors and their interactions with 28 genome‐wide significant common genetic variants. The interaction analyses were based on data from 9971 ovarian cancer cases and 15,566 controls from 17 case–control studies. Likelihood ratio and Wald tests for multiplicative interaction and for relative excess risk due to additive interaction were used. The top multiplicative interaction was noted between oral contraceptive pill (OCP) use (ever vs. never) and rs13255292 (p value = 3.48 × 10−4). Among women with the TT genotype for this variant, the odds ratio for OCP use was 0.53 (95% CI = 0.46–0.60) compared to 0.71 (95%CI = 0.66–0.77) for women with the CC genotype. When stratified by duration of OCP use, women with 1–5 years of OCP use exhibited differential protective benefit across genotypes. However, no interaction on either the multiplicative or additive scale was found to be statistically significant after multiple testing correction. The results suggest that OCP use may offer increased benefit for women who are carriers of the T allele in rs13255292. On the other hand, for women carrying the C allele in this variant, longer (5+ years) use of OCP may reduce the impact of carrying the risk allele of this SNP. Replication of this finding is needed. The study presents a comprehensive analytic framework for conducting gene–environment analysis in ovarian cancer.

Abstract Genome-wide association studies have identified 40 ovarian cancer risk loci. However, the mechanisms underlying these associations remain elusive. In this study, we conducted a two-pronged approach to identify candidate causal SNPs and assess underlying biological mechanisms at chromosome 9p22.2, the first and most statistically significant associated locus for ovarian cancer susceptibility. Three transcriptional regulatory elements with allele-specific effects and a scaffold/matrix attachment region were characterized and, through physical DNA interactions, BNC2 was established as the most likely target gene. We determined the consensus binding sequence for BNC2 in vitro, verified its enrichment in BNC2 ChIP-seq regions, and validated a set of its downstream target genes. Fine-mapping by dense regional genotyping in over 15,000 ovarian cancer cases and 30,000 controls identified SNPs in the scaffold/matrix attachment region as among the most likely causal variants. This study reveals a comprehensive regulatory landscape at 9p22.2 and proposes a likely mechanism of susceptibility to ovarian cancer. Significance: Mapping the 9p22.2 ovarian cancer risk locus identifies BNC2 as an ovarian cancer risk gene. See related commentary by Choi and Brown, p. 439

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 Background The immunohistochemical (IHC) 4 + clinical treatment score (C) i.e. IHC4+C assay combines IHC assessment of four routinely measured markers in the clinical setting, including estrogen receptor (ER), progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and KI67, with other clinical parameters for aiding decision-making on adjuvant chemotherapy in hormone positive (HR+) breast cancer. It is low cost, widely available and has been shown to provide similar prognostic information to the 21-gene (Oncotype DX) test. Clinical use of IHC4+C is limited by the need for standardized scoring of the IHC markers. Here, we examine the prognostic value of a digital image analysis (DIA) based IHC4+C assay in two study populations and compare this with pathologists' visual assay in a subset of patients. We also evaluated its potential clinical relevance in relation to the Nottingham prognostic index [NPI; i.e. low (<3.4), intermediate (3.4-5.4) and high (>5.4)] and IHC defined breast cancer subtypes (i.e. luminal A-like and B-like). Methods Using DIA data on ER, PR, HER2 and KI67, in addition to clinicopathological characteristics of tumors, we computed DIA based IHC4+C (i.e. DIAIHC4+C) assay for 2,034 women from study populations in Poland (PBCS) and the United Kingdom (SEARCH). A total of 246 deaths occurred over a median follow-up period of 6.7 years. A relevant threshold was determined using quartile analysis and the relationship between dichotomous categories (>75th vs ≤ 75th percentile) of DIAIHC4+C and 10-year breast cancer specific survival (BCSS) was investigated using Kaplan-Meier survival curves and in Cox proportional hazard regression models, overall and stratified by study, NPI and subtype. Results DIAIHC4+C score was associated with worse 10-year BCSS overall [HR (95% CI) >75th vs ≤ 75th percentile = 2.92 (1.77, 4.82); p-value <0.0001] and in both the PBCS [3.48 (1.40, 8.65); p-value = 0.007] and SEARCH [3.87 (2.24, 6.68); p-value <0.0001] study populations. DIAIHC4+C was associated with significantly worse survival outcomes in women with intermediate NPI [1.94 (1.31, 2.88); p-value = 0.001] and luminal A-like subtype of HR+ tumors [5.37 (3.50, 8.26); p-value <0.0001]. DIA and visual IHC4+C results were highly correlated (r = 0.95) and survival associations and model fit parameters were similar for DIA [3.82 (2.08, 7.00); p-value < 0.0001; AIC = 472.3] and visual [4.73 (2.55, 8.76); p-value < 0.0001; AIC = 466.1] IHC4+C scores amongst 444 patients with data from both sources. Conclusion These findings are in support of the prognostic value of DIAIHC4+C assay in HR+ breast cancer and demonstrate its potential clinical relevance as an adjunct to the NPI. Given the advantages of being high-throughput and highly reproducible, DIA based methods could find use in the standardization of the analytical processes for the IHC4+C assay. Citation Format: Mustapha Abubakar, Paul D. Pharoah, Mark E. Sherman, Montserrat Garcia-Closas, Mitch Dowsett. Digital image analysis based IHC4+C assay and prognosis in hormone receptor-positive breast cancer [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 2971.

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.