García‐Closas, Montserrat

Abstract BackgroundCopy number variants (CNVs) are pervasive in the human genome but potential disease associations with rare CNVs have not been comprehensively assessed in large datasets. We analysed rare CNVs in genes and non-coding regions for 86,788 breast cancer cases and 76,122 controls of European ancestry with genome-wide array data.ResultsGene burden tests detected the strongest association for deletions in BRCA1 (P= 3.7E-18). Nine other genes were associated with a p-value < 0.01 including known susceptibility genes CHEK2 (P= 0.0008), ATM (P= 0.002) and BRCA2 (P= 0.008). Outside the known genes we detected associations with p-values < 0.001 for either overall or subtype-specific breast cancer at nine deletion regions and four duplication regions. Three of the deletion regions were in established common susceptibility loci.ConclusionsThis is the first genome-wide analysis of rare CNVs in a large breast cancer case-control dataset. We detected associations with exonic deletions in established breast cancer susceptibility genes. We also detected suggestive associations with non-coding CNVs in known and novel loci with large effects sizes. Larger sample sizes will be required to reach robust levels of statistical significance.

Abstract Background: Use of menopausal hormone therapy (MHT) is associated with increased risk for breast cancer. However, the relevant mechanisms and its interaction with genetic variants are not fully understood. Methods: We conducted a genome-wide interaction analysis between MHT use and genetic variants for breast cancer risk in 27,585 cases and 34,785 controls from 26 observational studies. All women were post-menopausal and of European ancestry. Multivariable logistic regression models were used to test for multiplicative interactions between genetic variants and current MHT use. We considered interaction p-values<5x10-8 as genome-wide significant, and p-values<1x10-5 as suggestive. Linkage disequilibrium (LD)-based clumping was performed to identify independent candidate variants.Results: None of the 9.7 million genetic variants tested for interactions with MHT use reached genome-wide significance. Only 213 variants, representing 18 independent loci, had p-values<1x105. The strongest evidence was found for rs4674019 (p-value=2.27x10-7), which showed genome-wide significant interaction (p-value=3.8x10-8) with current MHT use when analysis was restricted to population-based studies only. Limiting the analyses to combined estrogen-progesterone MHT use only or to estrogen receptor (ER) positive cases did not identify any genome-wide significant evidence of interactions. Conclusions: In this large genome-wide SNP-MHT interaction study of breast cancer, we found no strong support for common genetic variants modifying the effect of MHT on breast cancer risk. These results suggest that common genetic variation has limited impact on the observed MHT–breast cancer risk association.

Abstract Background Rigorous evaluation of the calibration and discrimination of breast-cancer risk-prediction models in prospective cohorts is critical for applications under clinical guidelines. We comprehensively evaluated an integrated model incorporating classical risk factors and a 313-variant polygenic risk score (PRS) to predict breast-cancer risk. Methods Fifteen prospective cohorts from six countries with 239 340 women (7646 incident breast-cancer cases) of European ancestry aged 19–75 years were included. Calibration of 5-year risk was assessed by comparing expected and observed proportions of cases overall and within risk categories. Risk stratification for women of European ancestry aged 50–70 years in those countries was evaluated by the proportion of women and future cases crossing clinically relevant risk thresholds. Results Among women &lt;50 years old, the median (range) expected-to-observed ratio for the integrated model across 15 cohorts was 0.9 (0.7–1.0) overall and 0.9 (0.7–1.4) at the highest-risk decile; among women ≥50 years old, these were 1.0 (0.7–1.3) and 1.2 (0.7–1.6), respectively. The proportion of women identified above a 3% 5-year risk threshold (used for recommending risk-reducing medications in the USA) ranged from 7.0% in Germany (∼841 000 of 12 million) to 17.7% in the USA (∼5.3 of 30 million). At this threshold, 14.7% of US women were reclassified by adding the PRS to classical risk factors, with identification of 12.2% of additional future cases. Conclusion Integrating a 313-variant PRS with classical risk factors can improve the identification of European-ancestry women at elevated risk who could benefit from targeted risk-reducing strategies under current clinical guidelines.

Abstract Background: It is not known whether modifiable lifestyle factors that predict survival after invasive breast cancer differ by subtype. Methods: We analyzed data for 121,435 women diagnosed with breast cancer from 67 studies in the Breast Cancer Association Consortium with 16,890 deaths (8,554 breast cancer specific) over 10 years. Cox regression was used to estimate associations between risk factors and 10-year all-cause mortality and breast cancer–specific mortality overall, by estrogen receptor (ER) status, and by intrinsic-like subtype. Results: There was no evidence of heterogeneous associations between risk factors and mortality by subtype (Padj &gt; 0.30). The strongest associations were between all-cause mortality and BMI ≥30 versus 18.5–25 kg/m2 [HR (95% confidence interval (CI), 1.19 (1.06–1.34)]; current versus never smoking [1.37 (1.27–1.47)], high versus low physical activity [0.43 (0.21–0.86)], age ≥30 years versus &lt;20 years at first pregnancy [0.79 (0.72–0.86)]; &gt;0–&lt;5 years versus ≥10 years since last full-term birth [1.31 (1.11–1.55)]; ever versus never use of oral contraceptives [0.91 (0.87–0.96)]; ever versus never use of menopausal hormone therapy, including current estrogen–progestin therapy [0.61 (0.54–0.69)]. Similar associations with breast cancer mortality were weaker; for example, 1.11 (1.02–1.21) for current versus never smoking. Conclusions: We confirm associations between modifiable lifestyle factors and 10-year all-cause mortality. There was no strong evidence that associations differed by ER status or intrinsic-like subtype. Impact: Given the large dataset and lack of evidence that associations between modifiable risk factors and 10-year mortality differed by subtype, these associations could be cautiously used in prognostication models to inform patient-centered care.

Abstract Blood lipids have been associated with the development of a range of cancers, including breast, lung and colorectal cancer. For endometrial cancer, observational studies have reported inconsistent associations between blood lipids and cancer risk. To reduce biases from unmeasured confounding, we performed a bidirectional, two‐sample Mendelian randomization analysis to investigate the relationship between levels of three blood lipids (low‐density lipoprotein [LDL] and high‐density lipoprotein [HDL] cholesterol, and triglycerides) and endometrial cancer risk. Genetic variants associated with each of these blood lipid levels ( P < 5 × 10 −8 ) were identified as instrumental variables, and assessed using genome‐wide association study data from the Endometrial Cancer Association Consortium (12 906 cases and 108 979 controls) and the Global Lipids Genetic Consortium (n = 188 578). Mendelian randomization analyses found genetically raised LDL cholesterol levels to be associated with lower risks of endometrial cancer of all histologies combined, and of endometrioid and non‐endometrioid subtypes. Conversely, higher genetically predicted HDL cholesterol levels were associated with increased risk of non‐endometrioid endometrial cancer. After accounting for the potential confounding role of obesity (as measured by genetic variants associated with body mass index), the association between genetically predicted increased LDL cholesterol levels and lower endometrial cancer risk remained significant, especially for non‐endometrioid endometrial cancer. There was no evidence to support a role for triglycerides in endometrial cancer development. Our study supports a role for LDL and HDL cholesterol in the development of non‐endometrioid endometrial cancer. Further studies are required to understand the mechanisms underlying these findings.

Abstract We evaluated the joint associations between a new 313-variant PRS (PRS313) and questionnaire-based breast cancer risk factors for women of European ancestry, using 72 284 cases and 80 354 controls from the Breast Cancer Association Consortium. Interactions were evaluated using standard logistic regression and a newly developed case-only method for breast cancer risk overall and by estrogen receptor status. After accounting for multiple testing, we did not find evidence that per-standard deviation PRS313 odds ratio differed across strata defined by individual risk factors. Goodness-of-fit tests did not reject the assumption of a multiplicative model between PRS313 and each risk factor. Variation in projected absolute lifetime risk of breast cancer associated with classical risk factors was greater for women with higher genetic risk (PRS313 and family history) and, on average, 17.5% higher in the highest vs lowest deciles of genetic risk. These findings have implications for risk prevention for women at increased risk of breast cancer.

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 &lt;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 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 (&lt;3.4), intermediate (3.4-5.4) and high (&gt;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 (&gt;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) &gt;75th vs ≤ 75th percentile = 2.92 (1.77, 4.82); p-value &lt;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 &lt;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 &lt;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 &lt; 0.0001; AIC = 472.3] and visual [4.73 (2.55, 8.76); p-value &lt; 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 Introduction: Genome-wide association studies over the past 10 years have identified over 100 common variants associated with breast cancer risk. The association of these polymorphisms varies according to estrogen receptor status of the tumor, with some associated primarily with ER+ or ER- disease, or both subtypes. However, ER expression in breast cancer cells is continuous. The aim of this study was to establish whether the association between polygenic risk scores and ER expression in breast cancer is continuous in nature. Methods: Data were taken from the international Breast Cancer Association Consortium (BCAC). Germline genotype data for 75 SNPs previously reported to be associated with breast cancer were available for 4,999 breast cancer cases for which ER expression data were also available. Automated ER expression on a continuous scale was measured using two image analysis algorithms: Astrogrid and Ariol. A third form of ER expression score, the manual Allred ordinal, was used. Genotype data were used to calculate polygenic risk scores using the equation: PRS = ∑ βx where β is the log-odds of the odds ratio (OR) of breast cancer associated with that single nucleotide polymorphism (SNP), and x is the number of variant alleles at a SNP locus. Five different PRSs were calculated: i) SNPs associated with both ER+ and ER- disease (PRS-A, n= 74) ii) SNPs associated with ER+ disease at P&lt;0.01 (PRS-B, n= 65) iii) SNPs associated with ER- disease at P&lt;0.01 (PRS-C, n= 43) iv) SNPs associated with ER+ but not ER- disease (PRS-D, n= 30) v) SNPs associated with ER- but not ER+ disease (PRS-E, n = 8) The ER-specific weights for the PRS were taken from previous OR estimates. Linear regression models adjusted for study were used to test for a relationship between the five PRSs and ER expression measured on the Ariol, Astrogrid and Allred scales. Results: There was a strong positive correlation found between PRS-B against Astrogrid (β=0.0176, P=4.06x10-4), a weak negative correlation found with PRS-C (β=-0.012, P=0.019) and a weak positive correlation for PRS-A (β=0.012, P=0.0144). PRS-D and E showed far stronger association with Astrogrid than PRS-B and C, respectively (β=0.0186 and -0.017, P=1.97x10-4 and 8.1x10-4). There was a qualitatively similar association when the PRSs were correlated with Ariol and Allred, apart from the stronger correlations with PRS-B (β=91.2 and 0.0853, P=5.52x10-4 and 1.33x10-3). Conclusion: We have shown a robust association between PRS and ER expression in the directions we predicted. This suggests that the continuous variation in ER expression is biologically plausible and linked to genetic polymorphism. The weakening of association with PRS-A and the stronger associations with PRS-D and E show that having SNPs that are associated with both ER subtypes in the same PRS weakens significance, and suggests that subtype-specific PRSs may be improved on by using SNPs exclusive to the respective subtype. Citation Format: Anamay H. Shetty, Raza Ali, Fiona Blows, Carlos Caldas, Nicholas Inard, Javier Benítez, Jenny Chang-Claude, Georgia Chenevix-Trench, Fergus Couch, Angela Cox, Peter Devilee, Douglas Easton, Montserrat Garcia-Closas, Arto Mannermaa, Heli Nevanlinna, Marjanka Schmidt, Paul Pharoah. Patient level polygenic risk scores and continuous estrogen receptor expression in 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 2211.