Validation of the BOADICEA model and a 313-variant polygenic risk score for breast cancer risk prediction in a Dutch prospective cohort

Purpose: We evaluated the performance of the recently extended Breast and Ovarian Analysis of Disease Incidence and Carrier Estimation Algorithm (BOADICEA version 5) in a Dutch prospective cohort, using a polygenic risk score (PRS) based on 313 breast cancer (BC)–associated variants (PRS313) and other, nongenetic risk factors. Methods: Since 1989, 6522 women without BC aged 45 or older of European descent have been included in the Rotterdam Study. The PRS313 was calculated per 1 SD in controls from the Breast Cancer Association Consortium (BCAC). Cox regression analysis was performed to estimate the association between the PRS313 and incident BC risk. Cumulative 10-year risks were calculated with BOADICEA including different sets of variables (age, risk factors and PRS313). C-statistics were used to evaluate discriminative ability. Results: In total, 320 women developed BC. The PRS313 was significantly associated with BC (hazard ratio [HR] per SD of 1.56, 95% confidence interval [CI] [1.40–1.73]). Using 10-year risk estimates including age and the PRS313, other risk factors improved the discriminatory ability of the BOADICEA model marginally, from a C-statistic of 0.636 to 0.653. Conclusions: The effect size of the PRS313 is highly reproducible in the Dutch population. Our results validate the BOADICEA v5 model for BC risk assessment in the Dutch general population

Association of germline variation with the survival of women with BRCA1/2 pathogenic variants and breast cancer

Germline genetic variation has been suggested to influence the survival of breast cancer patients independently of tumor pathology. We have studied survival associations of genetic variants in two etiologically unique groups of breast cancer patients, the carriers of germline pathogenic variants in BRCA1 or BRCA2 genes. We found that rs57025206 was significantly associated with the overall survival, predicting higher mortality of BRCA1 carrier patients with estrogen receptor-negative breast cancer, with a hazard ratio 4.37 (95% confidence interval 3.03-6.30, P = 3.1 × 10-9). Multivariable analysis adjusted for tumor characteristics suggested that rs57025206 was an independent survival marker. In addition, our exploratory analyses suggest that the associations between genetic variants and breast cancer patient survival may depend on tumor biological subgroup and clinical patient characteristics

Tamoxifen and risk of contralateral breast cancer for BRCA1 and BRCA2 mutation carriers

Purpose To determine whether adjuvant tamoxifen treatment for breast cancer (BC) is associated with reduced contralateral breast cancer (CBC) risk for BRCA1 and/or BRCA2 mutation carriers. Methods Analysis of pooled observational cohort data, self-reported at enrollment and at follow-up from the International BRCA1, and BRCA2 Carrier Cohort Study, Kathleen Cuningham Foundation Consortium for Research into Familial Breast Cancer, and Breast Cancer Family Registry. Eligible women were BRCA1 and BRCA2 mutation carriers diagnosed with unilateral BC since 1970 and no other invasive cancer or tamoxifen use before first BC. Hazard ratios (HRs) for CBC associated with tamoxifen use were estimated using Cox regression, adjusting for year and age of diagnosis, country, and bilateral oophorectomy and censoring at contralateral mastectomy, death, or loss to follow-up. Results Of 1,583 BRCA1 and 881 BRCA2 mutation carriers, 383 (24%) and 454 (52%), respectively, took tamoxifen after first BC d

Inheritance of deleterious mutations at both BRCA1 and BRCA2 in an international sample of 32,295 women

Background: Most BRCA1 or BRCA2 mutation carriers have inherited a single (heterozygous) mutation. Transheterozygotes (TH) who have inherited deleterious mutations in both BRCA1 and BRCA2 are rare, and the consequences of transheterozygosity are poorly understood. Methods: From 32,295 female BRCA1/2 mutation carriers, we identified 93 TH (0.3 %). "Cases" were defined as TH, and "controls" were single mutations at BRCA1 (SH1) or BRCA2 (SH2). Matched SH1 "controls" carried a BRCA1 mutation found in the TH "case". Matched SH2 "controls" carried a BRCA2 mutation found in the TH "case". After matching the TH carriers with SH1 or SH2, 91 TH were matched to 9316 SH1, and 89 TH were matched to 3370 SH2. Results: The majority of TH (45.2 %) involved the three common Jewish mutations. TH were more likely than SH1 and SH2 women to have been ever diagnosed with breast cancer (BC; p = 0.002). TH were more likely to be diagnosed with ovarian cancer (OC) than SH2 (p = 0.017), but not SH1. Age at BC diagnosis was the same in TH vs. SH1 (p = 0.231), but was on average 4.5 years younger in TH than in SH2 (p < 0.001). BC in TH was more likely to be estrogen receptor (ER) positive (p = 0.010) or progesterone receptor (PR) positive (p = 0.013) than in SH1, but less likely to be ER positive (p < 0.001) or PR positive (p = 0.012) than SH2. Among 15 tumors from TH patients, there was no clear pattern of loss of heterozygosity (LOH) for BRCA1 or BRCA2 in either BC or OC. Conclusions: Our observations suggest that clinical TH phenotypes resemble SH1. However, TH breast tumor marker characteristics are phenotypically intermediate to SH1 and SH2

Associations of common breast cancer susceptibility alleles with risk of breast cancer subtypes in BRCA1 and BRCA2 mutation carriers

Introduction: More than 70 common alleles are known to be involved in breast cancer (BC) susceptibility, and several exhibit significant heterogeneity in their associations with different BC subtypes. Although there are differences in the association patterns between BRCA1 and BRCA2 mutation carriers and the general population for several loci, no study has comprehensively evaluated the associations of all known BC susceptibility alleles with risk of BC subtypes in BRCA1 and BRCA2 carriers. Methods: We used data from 15,252 BRCA1 and 8,211 BRCA2 carriers to analyze the associations between approximately 200,000 genetic variants on the iCOGS array and risk of BC subtypes defined by estrogen receptor (ER), progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and triple-negative- (TN) status; morphologic subtypes; histological grade; and nodal involvement. Results: The estimated BC hazard ratios (HRs) for the 74 known BC alleles in BRCA1 carriers exhibited moderate correlations with the corresponding odds ratios from the general population. However, their associations with ER-positive BC in BRCA1 carriers were more consistent with the ER-positive as

Assessing associations between the AURKAHMMR-TPX2-TUBG1 functional module and breast cancer risk in BRCA1/2 mutation carriers

While interplay between BRCA1 and AURKA-RHAMM-TPX2-TUBG1 regulates mammary epithelial polarization, common genetic variation in HMMR (gene product RHAMM) may be associated with risk of breast cancer in BRCA1 mutation carriers. Following on these observations, we further assessed the link between the AURKA-HMMR-TPX2-TUBG1 functional module and risk of breast cancer in BRCA1 or BRCA2 mutation carriers. Forty-one single nucleotide polymorphisms (SNPs) were genotyped in 15,252 BRCA1 and 8,211 BRCA2 mutation carriers and subsequently analyzed using a retrospective likelihood appr

Fine-Scale Mapping of the 4q24 Locus Identifies Two Independent Loci Associated with Breast Cancer Risk

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

Functional mechanisms underlying pleiotropic risk alleles at the 19p13.1 breast-ovarian cancer susceptibility locus

A locus at 19p13 is associated with breast cancer (BC) and ovarian cancer (OC) risk. Here we analyse 438 SNPs in this region in 46,451 BC and 15,438 OC cases, 15,252 BRCA1 mutation carriers and 73,444 controls and identify 13 candidate causal SNPs associated with serous OC (P=9.2 × 10-20), ER-negative BC (P=1.1 × 10-13), BRCA1-associated BC (P=7.7 × 10-16) and triple negative BC (P-diff=2 × 10-5). Genotype-gene expression associations are identified for candidate target genes ANKLE1 (P=2 × 10-3) and ABHD8 (P<2 × 10-3). Chromosome conformation capture identifies interactions between four candidate SNPs and ABHD8, and luciferase assays indicate six risk alleles increased transactivation of the ADHD8 promoter. Targeted deletion of a region containing risk SNP rs56069439 in a putative enhancer induces ANKLE1 downregulation; and mRNA stability assays indicate functional effects for an ANKLE1 3′-UTR SNP. Altogether, these data suggest that multiple SNPs at 19p13 regulate ABHD8 and perhaps ANKLE1 expression, and indicate common mechanisms underlying breast and ovarian cancer risk

Refined histopathological predictors of BRCA1 and BRCA2 mutation status: A large-scale analysis of breast cancer characteristics from the BCAC, CIMBA, and ENIGMA consortia

Introduction: The distribution of histopathological features of invasive breast tumors in BRCA1 or BRCA2 germline mutation carriers differs from that of individuals with no known mutation. Histopathological features thus have utility for mutation prediction, including statistical modeling to assess pathogenicity of BRCA1 or BRCA2 variants of uncertain clinical significance. We analyzed large pathology datasets accrued by the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) and the Breast Cancer Association Consortium (BCAC) to reassess histopathological predictors of BRCA1 and BRCA2 mutation status, and provide robust likelihood ratio (LR) estimates for statistical modeling. Methods: Selection criteria for study/center inclusion were estrogen receptor (ER) status or grade data available for invasive breast cancer diagnosed younger than 70 years. The dataset included 4,477 BRCA1 mutation carriers, 2,565 BRCA2 mutation carriers, and 47,565 BCAC breast cancer cases. Country-stratified estimates of the

Genome-Wide Association Study in BRCA1 Mutation Carriers Identifies Novel Loci Associated with Breast and Ovarian Cancer Risk

BRCA1-associated breast and ovarian cancer risks can be modified by common genetic variants. To identify further cancer risk-modifying loci, we performed a multi-stage GWAS of 11,705 BRCA1 carriers (of whom 5,920 were diagnosed with breast and 1,839 were diagnosed with ovarian cancer), with a further replication in an additional sample of 2,646 BRCA1 carriers. We identified a novel breast cancer risk modifier locus at 1q32 for BRCA1 carriers (rs2290854, P = 2.7×10-8, HR = 1.14, 95% CI: 1.09-1.20). In addition, we identified two novel ovarian cancer risk modifier loci: 17q21.31 (rs17631303, P = 1.4×10-8, HR = 1.27, 95% CI: 1.17-1.38) and 4q32.3 (rs4691139, P = 3.4×10-8, HR = 1.20, 95% CI: 1.17-1.38). The 4q32.3 locus was not associated with ovarian cancer risk in the general population or BRCA2 carriers, suggesting a BRCA1-specific associat