Schoemaker, Minouk J.

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 Background: Prospective validation of breast cancer risk models integrating classical risk factors and genetic variants is required for risk-stratified prevention and screening strategies. The objective of this study was to validate a breast cancer risk model integrating classical risk factors and a 313-variant polygenic risk score (PRS) in multiple prospective cohort studies, and to project five-year risk of breast cancer in six different countries. Methods: The study population included 7,529 cases and 230,103 controls from 14 prospective cohort studies in Australia, Germany, the Netherlands, Sweden, UK, and USA. We used the Individualized Coherent Absolute Risk Estimator (iCARE) tool for risk model building, validation, and risk projection. Expected five-year risk of invasive or in situ breast cancer was compared to observed risk, overall and within deciles of expected risk using goodness of fit statistics. We evaluated calibration of the relative risk through meta-analysis across cohorts, and of the absolute risk within each cohort. Model discrimination was evaluated using the area under the curve (AUC), and percentages of women crossing risk thresholds. Projections of five-year risk distributions were estimated for women of European ancestry aged 50-70 years in the general populations of these six countries. Results: Analysis showed overall good calibration of the integrated iCARE-based model relative risk for both women younger than 50 years (χ2=14.9, P=0.09) and aged 50 years or older (χ2=14.6, P=0.10), with a small overestimation of risk for women in the highest decile of expected risk (RR = 3.5 expected vs 2.3 (95% CI 1.6 to 3.2) observed for women <50 years; and 2.8 expected vs 2.3 (95% CI 2.0 to 2.7) observed for women 50+ years). The age-adjusted AUCs for the integrated model were 63.1 (95% CI 60.9 to 65.3) and 62.9 (95% CI 61.8 to 64.0), for the two age groups respectively. The calibration of absolute risk showed substantial variation across cohorts, particularly for the older group, but had no systematic bias. Model based projections in the general populations showed that compared to the population average, women in the 1st and 99th percentiles of the integrated risk score had relative risks 0.19 and 3.56 respectively. The proportion of women of European ancestry aged 50-70 years with a five-year risk greater than 3% (threshold for consideration of risk-lowering drugs by U.S. Preventive Services Task Force) ranged from 7.1% in Germany to 18.2% in the US, which corresponds to ~5.5 million women in the US. Conclusions: Five-year risk predictions from a model with classical risk factors and PRS are well calibrated and provide substantial risk stratification across multiple cohorts in six different countries. Further studies are needed to evaluate the clinical utility of the validated model for risk stratified screening and prevention of breast cancer. Citation Format: Parichoy Pal Choudhury, Amber Wilcox, Chi Gao, Brian Carter, Anika Husing, Mark Brook, Mikael Eriksson, Kara Martin, Chris Scott, Min Shi, Thomas Ahearn, Michael Jones, Nick Orr, Minouk Schoemaker, Kamila Czene, Jenny Chang-Claude, Jacques Simard, Doug Easton, Marjanka K. Schmidt, Dale Sandler, Clarice R. Weinberg, Celine Vachon, Roger Milne, Per Hall, Anthony Swerdlow, Rudolph Kaaks, Myrto Barrdahl, Mia Gaudet, Antonis Antoniou, Peter Kraft, Montserrat Garcia-Closas, Nilanjan Chatterjee. Validation of breast cancer risk model incorporating classical risk factors and polygenic risk scores in 14 prospective cohort studies in 6 countries [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 962.

Abstract Genome-wide association studies (GWAS) have so far identified 25 loci associated with glioma risk, with most showing specificity for either glioblastoma (GBM) or non-GBM tumors. The majority of these GWAS susceptibility variants reside in noncoding regions and the causal genes underlying the associations are largely unknown. Here we performed a transcriptome-wide association study to search for novel risk loci and candidate causal genes at known GWAS loci using Genotype-Tissue Expression Project (GTEx) data to predict cis-predicted gene expression in relation to GBM and non-GBM risk in conjunction with GWAS summary statistics on 12,488 glioma cases (6,183 GBM and 5,820 non-GBM) and 18,169 controls. Imposing a Bonferroni-corrected significance level of P < 5.69 × 10−6, we identified 31 genes, including GALNT6 at 12q13.33, as a candidate novel risk locus for GBM (mean Z = 4.43; P = 5.68 × 10−6). GALNT6 resides at least 55 Mb away from any previously identified glioma risk variant, while all other 30 significantly associated genes were located within 1 Mb of known GWAS-identified loci and were not significant after conditioning on the known GWAS-identified variants. These data identify a novel locus (GALNT6 at 12q13.33) and 30 genes at 12 known glioma risk loci associated with glioma risk, providing further insights into glioma tumorigenesis. Significance: This study identifies new genes associated with glioma risk, increasing understanding of how these tumors develop.

Limited epidemiological evidence suggests that the etiology of hormone receptor positive (HR+) breast cancer may differ by levels of histologic grade and proliferation. We pooled risk factor and pathology data on 5,905 HR+ breast cancer cases and 26,281 controls from 11 epidemiological studies. Proliferation was determined by centralized automated measures of KI67 in tissue microarrays. Odds ratios (OR), 95% confidence intervals (CI) and p‐values for case–case and case–control comparisons for risk factors in relation to levels of grade and quartiles (Q1–Q4) of KI67 were estimated using polytomous logistic regression models. Case–case comparisons showed associations between nulliparity and high KI67 [OR (95% CI) for Q4 vs. Q1 = 1.54 (1.22, 1.95)]; obesity and high grade [grade 3 vs. 1 = 1.68 (1.31, 2.16)] and current use of combined hormone therapy (HT) and low grade [grade 3 vs. 1 = 0.27 (0.16, 0.44)] tumors. In case–control comparisons, nulliparity was associated with elevated risk of tumors with high but not low levels of proliferation [1.43 (1.14, 1.81) for KI67 Q4 vs. 0.83 (0.60, 1.14) for KI67 Q1]; obesity among women ≥50 years with high but not low grade tumors [1.55 (1.17, 2.06) for grade 3 vs. 0.88 (0.66, 1.16) for grade 1] and HT with low but not high grade tumors [3.07 (2.22, 4.23) for grade 1 vs. 0.85 (0.55, 1.30) for grade 3]. Menarcheal age and family history were similarly associated with HR+ tumors of different grade or KI67 levels. These findings provide insights into the etiologic heterogeneity of HR+ tumors.

Abstract BACKGROUND: Glioma accounts for ~27% of all primary brain tumors and is responsible for ~13,000 cancer-related deaths in the US each year. Glioma tumors can be broadly classified into glioblastoma (GBM) and lower-grade non-GBM. Typically gliomas have a poor prognosis irrespective of medical care, with the most common form, GBM, having a five-year survival rate of only 5%. While genome-wide association studies (GWAS) have transformed our understanding of glioma susceptibility, individual studies have had limited power to identify risk loci. METHODS: We performed the largest glioma GWAS to date, comprising a meta-analysis of six existing GWAS (6,405 cases, 14,100 controls) as well as new GWAS from the Glioma International Case Control Consortium (GICC; 4,572 cases and 3,286 controls) and University of California, San Francisco (UCSF)-Mayo (1,519 cases, 804 controls), totaling 12,496 cases (6,191 classified as GBM, 5,819 as non-GBM) and 18,190 controls. RESULTS: We identified five new risk loci for GBM at 1p31.3 (rs12752552; near JAK1, P=2.04×10-9, odds ratio (OR)=1.22), 11q14.1 (rs11233250; P=9.95×10-10, OR=1.24), 16p13.3 (rs2562152; near MPG, P=1.93x10-8, OR=1.21), 16q12.1 (rs10852606; HEATR3, P=1.29×10-11, OR=1.18), 22q13.1 (rs2235573; P=1.76×10-10, OR=1.15) and eight for non-GBM at 1q32.1 (rs4252707; MDM4, P=3.34×10-9, OR=1.19), 1q44 (rs12076373; AKT3, P=2.63×10-10, OR=1.23), 2q33.3 (rs7572263; near IDH1, P=2.18×10-10, OR=1.20), 3p14.1 (rs11706832; LRIG1, P=7.66×10-9, OR=1.15), 10q24.33 (rs11598018; OBFC1, P=3.39×10-8, OR=1.14), 11q21 (rs7107785; P=3.87×10-10, OR=1.16), 14q12 (rs10131032; P=5.07x10-11, OR=1.33) and 16p13.3 (rs3751667; P=2.61×10-9, OR=1.18). Case-only analyses confirmed the specificity of 11q14.1, 16p13.3 and 22q13.1 associations for GBM and 1q44, 2q33.3, 3p14.1, 11q21 and 14q12 for non-GBM tumors. In the combined meta-analysis, among previously published glioma risk SNPs, those for all glioma at 17p13.1 (TP53), GBM at 5p15.33 (TERT), 7p11.2 (EGFR), 9p21.3 (CDKN2B-AS1) and 20q13.33 (RTEL1) and for non-GBM at 8q24.21 (CCDC26), 11q23.2, 11q23.3 (PHLDB1) and 15q24.2 (ETFA) showed even greater evidence for association. SNPs at 10q25.2 and 12q12.1 for non-GBM tumors retained genome-wide significance (i.e. P<5.0x10-8). Associations at the previously reported loci for GBM at 3q26.2 (near TERC) and 12q23.33 (POLR3B) did not retain statistical significance. CONCLUSIONS: Our findings substantiate genetic susceptibility to GBM and non-GBM glioma being highly distinct, consistent with their distinctive molecular profiles presumably resulting from different etiological pathways. Functional analyses should lead to further insights into the biological basis of the different glioma histologies. Such information can inform gene discovery initiatives and therefore have a measurable impact on the successful development of new therapeutic agents. Citation Format: Ben Kinnersley, Beatrice S. Melin, Jill S. Barnholtz-Sloan, Margaret R. Wrensch, Christoffer Johansen, Dora Il’yasova, Quinn Ostrom, and members of GICC, Karim Labreche, Jeanette E. Eckel-Passow, Paul A. Decker, Marianne Labussière, Ahmed Idbaih, Khe Hoang-Xuan, Anna-Luisa Di Stefano, Karima Mokhtari, Jean-Yves Delattre, Pilar Galan, Konstantinos Gousias, Johannes Schramm, Minouk J. Schoemaker, Sarah J. Fleming, Stefan Herms, Stefanie Heilmann, Marcus M. Nöthen, Heinz-Erich Wichmann, Stefan Schreiber, Anthony Swerdlow, Mark Lathrop, Matthias Simon, Marc Sanson, Preetha Rajaraman, Stephen Chanock, Martha Linet, Zhaoming Wang, Meredith Yeager, Rose K. Lai, Elizabeth B. Claus, Sara H. Olson, Robert B. Jenkins, Richard S. Houlston, Melissa L. Bondy. Genome-wide association study of glioma reveals specific differences in genetic susceptibility to glioblastoma and non-glioblastoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1302. doi:10.1158/1538-7445.AM2017-1302

Investigating the most likely causal variants identified by fine‐mapping analyses may improve the power to detect gene–environment interactions. We assessed the interplay between 70 single nucleotide polymorphisms identified by genetic fine‐scale mapping of susceptibility loci and 11 epidemiological breast cancer risk factors in relation to breast cancer. Analyses were conducted on up to 58,573 subjects (26,968 cases and 31,605 controls) from the Breast Cancer Association Consortium, in one of the largest studies of its kind. Analyses were carried out separately for estrogen receptor (ER) positive (ER+) and ER negative (ER–) disease. The Bayesian False Discovery Probability (BFDP) was computed to assess the noteworthiness of the results. Four potential gene–environment interactions were identified as noteworthy (BFDP < 0.80) when assuming a true prior interaction probability of 0.01. The strongest interaction result in relation to overall breast cancer risk was found between CFLAR‐rs7558475 and current smoking (ORint = 0.77, 95% CI: 0.67–0.88, pint = 1.8 × 10−4). The interaction with the strongest statistical evidence was found between 5q14‐rs7707921 and alcohol consumption (ORint =1.36, 95% CI: 1.16–1.59, pint = 1.9 × 10−5) in relation to ER– disease risk. The remaining two gene–environment interactions were also identified in relation to ER– breast cancer risk and were found between 3p21‐rs6796502 and age at menarche (ORint = 1.26, 95% CI: 1.12–1.43, pint =1.8 × 10−4) and between 8q23‐rs13267382 and age at first full‐term pregnancy (ORint = 0.89, 95% CI: 0.83–0.95, pint = 5.2 × 10−4). While these results do not suggest any strong gene–environment interactions, our results may still be useful to inform experimental studies. These may in turn, shed light on the potential interactions observed.

Abstract Breast cancer risks conferred by many germline missense variants in the BRCA1 and BRCA2 genes, often referred to as variants of uncertain significance (VUS), have not been established. In this study, associations between 19 BRCA1 and 33 BRCA2 missense substitution variants and breast cancer risk were investigated through a breast cancer case–control study using genotyping data from 38 studies of predominantly European ancestry (41,890 cases and 41,607 controls) and nine studies of Asian ancestry (6,269 cases and 6,624 controls). The BRCA2 c.9104A&gt;C, p.Tyr3035Ser (OR = 2.52; P = 0.04), and BRCA1 c.5096G&gt;A, p.Arg1699Gln (OR = 4.29; P = 0.009) variant were associated with moderately increased risks of breast cancer among Europeans, whereas BRCA2 c.7522G&gt;A, p.Gly2508Ser (OR = 2.68; P = 0.004), and c.8187G&gt;T, p.Lys2729Asn (OR = 1.4; P = 0.004) were associated with moderate and low risks of breast cancer among Asians. Functional characterization of the BRCA2 variants using four quantitative assays showed reduced BRCA2 activity for p.Tyr3035Ser compared with wild-type. Overall, our results show how BRCA2 missense variants that influence protein function can confer clinically relevant, moderately increased risks of breast cancer, with potential implications for risk management guidelines in women with these specific variants. Cancer Res; 77(11); 2789–99. ©2017 AACR.

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 &lt; 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 &lt; 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

AbstractAutomated methods are needed to facilitate high‐throughput and reproducible scoring of Ki67 and other markers in breast cancer tissue microarrays (TMAs) in large‐scale studies. To address this need, we developed an automated protocol for Ki67 scoring and evaluated its performance in studies from the Breast Cancer Association Consortium. We utilized 166 TMAs containing 16,953 tumour cores representing 9,059 breast cancer cases, from 13 studies, with information on other clinical and pathological characteristics. TMAs were stained for Ki67 using standard immunohistochemical procedures, and scanned and digitized using the Ariol system. An automated algorithm was developed for the scoring of Ki67, and scores were compared to computer assisted visual (CAV) scores in a subset of 15 TMAs in a training set. We also assessed the correlation between automated Ki67 scores and other clinical and pathological characteristics. Overall, we observed good discriminatory accuracy (AUC = 85%) and good agreement (kappa = 0.64) between the automated and CAV scoring methods in the training set. The performance of the automated method varied by TMA (kappa range= 0.37–0.87) and study (kappa range = 0.39–0.69). The automated method performed better in satisfactory cores (kappa = 0.68) than suboptimal (kappa = 0.51) cores (p‐value for comparison = 0.005); and among cores with higher total nuclei counted by the machine (4,000–4,500 cells: kappa = 0.78) than those with lower counts (50–500 cells: kappa = 0.41; p‐value = 0.010). Among the 9,059 cases in this study, the correlations between automated Ki67 and clinical and pathological characteristics were found to be in the expected directions. Our findings indicate that automated scoring of Ki67 can be an efficient method to obtain good quality data across large numbers of TMAs from multicentre studies. However, robust algorithm development and rigorous pre‐ and post‐analytical quality control procedures are necessary in order to ensure satisfactory performance.