Cytokines-activated nuclear IKKα-FAT10 pathway induces breast cancer tamoxifen-resistance.

Endocrine therapy that blocks estrogen signaling is the most effective treatment for patients with estrogen receptor positive (ER+) breast cancer. However, the efficacy of agents such as tamoxifen (Tam) is often compromised by the development of resistance. Here we report that cytokines-activated nuclear IKKα confers Tam resistance to ER+ breast cancer by inducing the expression of FAT10, and that the expression of FAT10 and nuclear IKKα in primary ER+ human breast cancer was correlated with lymphotoxin ß (LTB) expression and significantly associated with relapse and metastasis in patients treated with adjuvant mono-Tam. IKKα activation or enforced FAT10 expression promotes Tam-resistance while loss of IKKα or FAT10 augments Tam sensitivity. The induction of FAT10 by IKKα is mediated by the transcription factor Pax5, and coordinated via an IKKα-p53-miR-23a circuit in which activation of IKKα attenuates p53-directed repression of FAT10. Thus, our findings establish IKKα-to-FAT10 pathway as a new therapeutic target for the treatment of Tam-resistant ER+ breast cancer.

Nonlinear Mixed-Effects Model of Z-Endoxifen Concentrations in Tamoxifen-Treated Patients from the CEPAM Cohort.

Tamoxifen is widely used in patients with hormone receptor-positive breast cancer. The polymorphic enzyme CYP2D6 is primarily responsible for metabolic activation of tamoxifen, resulting in substantial interindividual variability of plasma concentrations of its most important metabolite, Z-endoxifen. The Z-endoxifen concentration thresholds below which tamoxifen treatment is less efficacious have been proposed but not validated, and prospective trials of individualized tamoxifen treatment to achieve Z-endoxifen concentration thresholds are considered infeasible. Therefore, we aim to validate the association between Z-endoxifen concentration and tamoxifen treatment outcomes, and identify a Z-endoxifen concentration threshold of tamoxifen efficacy, using pharmacometric modeling and simulation. As a first step, the CYP2D6 Endoxifen Percentage Activity Model (CEPAM) cohort was created by pooling data from 28 clinical studies (> 7,000 patients) with measured endoxifen plasma concentrations. After cleaning, data from 6,083 patients were used to develop a nonlinear mixed-effect (NLME) model for tamoxifen and Z-endoxifen pharmacokinetics that includes a conversion factor to allow inclusion of studies that measured total endoxifen but not Z-endoxifen. The final parent-metabolite NLME model confirmed the primary role of CYP2D6, and contributions from body weight, CYP2C9 phenotype, and co-medication with CYP2D6 inhibitors, on Z-endoxifen pharmacokinetics. Future work will use the model to simulate Z-endoxifen concentrations in patients receiving single agent tamoxifen treatment within large prospective clinical trials with long-term survival to identify the Z-endoxifen concentration threshold below which tamoxifen is less efficacious. Identification of this concentration threshold would allow personalized tamoxifen treatment to improve outcomes in patients with hormone receptor-positive breast cancer.

Differences in polygenic score distributions in European ancestry populations: implications for breast cancer risk prediction.

The 313-variant polygenic risk score (PRS313) provides a promising tool for breast cancer risk prediction. However, evaluation of the PRS313 across different European populations which could influence risk estimation has not been performed. Here, we explored the distribution of PRS313 across European populations using genotype data from 94,072 females without breast cancer, of European-ancestry from 21 countries participating in the Breast Cancer Association Consortium (BCAC) and 225,105 female participants from the UK Biobank. The mean PRS313 differed markedly across European countries, being highest in south-eastern Europe and lowest in north-western Europe. Using the overall European PRS313 distribution to categorise individuals leads to overestimation and underestimation of risk in some individuals from south-eastern and north-western countries, respectively. Adjustment for principal components explained most of the observed heterogeneity in mean PRS. Country-specific PRS distributions may be used to calibrate risk categories in individuals from different countries.

Polymorphisms in genes of melatonin biosynthesis and signaling support the light-at-night hypothesis for breast cancer.

Light-at-night triggers the decline of pineal gland melatonin biosynthesis and secretion and is an IARC-classified probable breast-cancer risk factor. We applied a large-scale molecular epidemiology approach to shed light on the putative role of melatonin in breast cancer. We investigated associations between breast-cancer risk and polymorphisms at genes of melatonin biosynthesis/signaling using a study population of 44,405 women from the Breast Cancer Association Consortium (22,992 cases, 21,413 population-based controls). Genotype data of 97 candidate single nucleotide polymorphisms (SNPs) at 18 defined gene regions were investigated for breast-cancer risk effects. We calculated adjusted odds ratios (ORs) and 95% confidence intervals (CI) by logistic regression for the main-effect analysis as well as stratified analyses by estrogen- and progesterone-receptor (ER, PR) status. SNP-SNP interactions were analyzed via a two-step procedure based on logic regression. The Bayesian false-discovery probability (BFDP) was used for all analyses to account for multiple testing. Noteworthy associations (BFDP < 0.8) included 10 linked SNPs in tryptophan hydroxylase 2 (TPH2) (e.g. rs1386492: OR = 1.07, 95% CI 1.02-1.12), and a SNP in the mitogen-activated protein kinase 8 (MAPK8) (rs10857561: OR = 1.11, 95% CI 1.04-1.18). The SNP-SNP interaction analysis revealed noteworthy interaction terms with TPH2- and MAPK-related SNPs (e.g. rs1386483R ∧ rs1473473D ∧ rs3729931D: OR = 1.20, 95% CI 1.09-1.32). In line with the light-at-night hypothesis that links shift work with elevated breast-cancer risks our results point to SNPs in TPH2 and MAPK-genes that may impact the intricate network of circadian regulation.

A genome-wide gene-environment interaction study of breast cancer risk for women of European ancestry.

BACKGROUND: Genome-wide studies of gene-environment interactions (G×E) may identify variants associated with disease risk in conjunction with lifestyle/environmental exposures. We conducted a genome-wide G×E analysis of ~ 7.6 million common variants and seven lifestyle/environmental risk factors for breast cancer risk overall and for estrogen receptor positive (ER +) breast cancer. METHODS: Analyses were conducted using 72,285 breast cancer cases and 80,354 controls of European ancestry from the Breast Cancer Association Consortium. Gene-environment interactions were evaluated using standard unconditional logistic regression models and likelihood ratio tests for breast cancer risk overall and for ER + breast cancer. Bayesian False Discovery Probability was employed to assess the noteworthiness of each SNP-risk factor pairs. RESULTS: Assuming a 1 × 10-5 prior probability of a true association for each SNP-risk factor pairs and a Bayesian False Discovery Probability < 15%, we identified two independent SNP-risk factor pairs: rs80018847(9p13)-LINGO2 and adult height in association with overall breast cancer risk (ORint = 0.94, 95% CI 0.92-0.96), and rs4770552(13q12)-SPATA13 and age at menarche for ER + breast cancer risk (ORint = 0.91, 95% CI 0.88-0.94). CONCLUSIONS: Overall, the contribution of G×E interactions to the heritability of breast cancer is very small. At the population level, multiplicative G×E interactions do not make an important contribution to risk prediction in breast cancer.

Association of the CHEK2 c.1100delC variant, radiotherapy, and systemic treatment with contralateral breast cancer risk and breast cancer-specific survival.

BACKGROUND: Breast cancer (BC) patients with a germline CHEK2 c.1100delC variant have an increased risk of contralateral BC (CBC) and worse BC-specific survival (BCSS) compared to non-carriers. AIM: To assessed the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS. METHODS: Analyses were based on 82,701 women diagnosed with a first primary invasive BC including 963 CHEK2 c.1100delC carriers; median follow-up was 9.1 years. Differential associations with treatment by CHEK2 c.1100delC status were tested by including interaction terms in a multivariable Cox regression model. A multi-state model was used for further insight into the relation between CHEK2 c.1100delC status, treatment, CBC risk and death. RESULTS: There was no evidence for differential associations of therapy with CBC risk by CHEK2 c.1100delC status. The strongest association with reduced CBC risk was observed for the combination of chemotherapy and endocrine therapy [HR (95% CI): 0.66 (0.55-0.78)]. No association was observed with radiotherapy. Results from the multi-state model showed shorter BCSS for CHEK2 c.1100delC carriers versus non-carriers also after accounting for CBC occurrence [HR (95% CI): 1.30 (1.09-1.56)]. CONCLUSION: Systemic therapy was associated with reduced CBC risk irrespective of CHEK2 c.1100delC status. Moreover, CHEK2 c.1100delC carriers had shorter BCSS, which appears not to be fully explained by their CBC risk.

The impact of coding germline variants on contralateral breast cancer risk and survival.

Evidence linking coding germline variants in breast cancer (BC)-susceptibility genes other than BRCA1, BRCA2, and CHEK2 with contralateral breast cancer (CBC) risk and breast cancer-specific survival (BCSS) is scarce. The aim of this study was to assess the association of protein-truncating variants (PTVs) and rare missense variants (MSVs) in nine known (ATM, BARD1, BRCA1, BRCA2, CHEK2, PALB2, RAD51C, RAD51D, and TP53) and 25 suspected BC-susceptibility genes with CBC risk and BCSS. Hazard ratios (HRs) and 95% confidence intervals (CIs) were estimated with Cox regression models. Analyses included 34,401 women of European ancestry diagnosed with BC, including 676 CBCs and 3,449 BC deaths; the median follow-up was 10.9 years. Subtype analyses were based on estrogen receptor (ER) status of the first BC. Combined PTVs and pathogenic/likely pathogenic MSVs in BRCA1, BRCA2, and TP53 and PTVs in CHEK2 and PALB2 were associated with increased CBC risk [HRs (95% CIs): 2.88 (1.70-4.87), 2.31 (1.39-3.85), 8.29 (2.53-27.21), 2.25 (1.55-3.27), and 2.67 (1.33-5.35), respectively]. The strongest evidence of association with BCSS was for PTVs and pathogenic/likely pathogenic MSVs in BRCA2 (ER-positive BC) and TP53 and PTVs in CHEK2 [HRs (95% CIs): 1.53 (1.13-2.07), 2.08 (0.95-4.57), and 1.39 (1.13-1.72), respectively, after adjusting for tumor characteristics and treatment]. HRs were essentially unchanged when censoring for CBC, suggesting that these associations are not completely explained by increased CBC risk, tumor characteristics, or treatment. There was limited evidence of associations of PTVs and/or rare MSVs with CBC risk or BCSS for the 25 suspected BC genes. The CBC findings are relevant to treatment decisions, follow-up, and screening after BC diagnosis.

Association of the CHEK2 c.1100delC variant, radiotherapy, and systemic treatment with contralateral breast cancer risk and breast cancer-specific survival.

Breast cancer (BC) patients with a germline CHEK2 c.1100delC variant have an increased risk of contralateral BC (CBC) and worse BC-specific survival (BCSS) compared to non-carriers. We aimed to assess the associations of CHEK2 c.1100delC, radiotherapy, and systemic treatment with CBC risk and BCSS. Analyses were based on 82,701 women diagnosed with invasive BC including 963 CHEK2 c.1100delC carriers; median follow-up was 9.1 years. Differential associations of treatment by CHEK2 c.1100delC status were tested by including interaction terms in a multivariable Cox regression model. A multi-state model was used for further insight into the relation between CHEK2 c.1100delC status, treatment, CBC risk and death. There was no evidence for differential associations of therapy with CBC risk by CHEK2 c.1100delC status The strongest association with reduced CBC risk was observed for the combination of chemotherapy and endocrine therapy [HR(95%CI): 0.66 (0.55-0.78)]. No association was observed with radiotherapy. Results from the multi-state model showed shorter BCSS for CHEK2 c.1100delC carriers versus non-carriers also after accounting for CBC occurrence [HR(95%CI) :1.30 (1.09-1.56)]. In conclusion, systemic therapy was associated with reduced CBC risk irrespective of CHEK2 c.1100delC status. Moreover, CHEK2 c.1100delC carriers had shorter BCSS, which appears not to be fully explained by their CBC risk. (Main MS: 3201 words).

Cross-Ancestry Genome-Wide Association Study Defines the Extended CYP2D6 Locus as the Principal Genetic Determinant of Endoxifen Plasma Concentrations.

The therapeutic efficacy of tamoxifen is predominantly mediated by its active metabolites 4-hydroxy-tamoxifen and endoxifen, whose formation is catalyzed by the polymorphic cytochrome P450 2D6 (CYP2D6). Yet, known CYP2D6 polymorphisms only partially determine metabolite concentrations in vivo. We performed the first cross-ancestry genome-wide association study with well-characterized patients of European, Middle-Eastern, and Asian descent (n = 497) to identify genetic factors impacting active and parent metabolite formation. Genome-wide significant variants were functionally evaluated in an independent liver cohort (n = 149) and in silico. Metabolite prediction models were validated in two independent European breast cancer cohorts (n = 287, n = 189). Within a single 1-megabase (Mb) region of chromosome 22q13 encompassing the CYP2D6 gene, 589 variants were significantly associated with tamoxifen metabolite concentrations, particularly endoxifen and metabolic ratio (MR) endoxifen/N-desmethyltamoxifen (minimal P = 5.4E-35 and 2.5E-65, respectively). Previously suggested other loci were not confirmed. Functional analyses revealed 66% of associated, mostly intergenic variants to be significantly correlated with hepatic CYP2D6 activity or expression (ρ = 0.35 to -0.52), and six hotspot regions in the extended 22q13 locus impacting gene regulatory function. Machine learning models based on hotspot variants (n = 12) plus CYP2D6 activity score (AS) increased the explained variability (~ 9%) compared with AS alone, explaining up to 49% (median R2 ) and 72% of the variability in endoxifen and MR endoxifen/N-desmethyltamoxifen, respectively. Our findings suggest that the extended CYP2D6 locus at 22q13 is the principal genetic determinant of endoxifen plasma concentration. Long-distance haplotypes connecting CYP2D6 with adjacent regulatory sites and nongenetic factors may account for the unexplained portion of variability.

DNA methylation biomarkers for noninvasive detection of triple-negative breast cancer using liquid biopsy.

Noninvasive detection of aberrant DNA methylation could provide invaluable biomarkers for earlier detection of triple-negative breast cancer (TNBC) which could help clinicians with easier and more efficient treatment options. We evaluated genome-wide DNA methylation data derived from TNBC and normal breast tissues, peripheral blood of TNBC cases and controls and reference samples of sorted blood and mammary cells. Differentially methylated regions (DMRs) between TNBC and normal breast tissues were stringently selected, verified and externally validated. A machine-learning algorithm was applied to select the top DMRs, which then were evaluated on plasma-derived circulating cell-free DNA (cfDNA) samples of TNBC patients and healthy controls. We identified 23 DMRs accounting for the methylation profile of blood cells and reference mammary cells and then selected six top DMRs for cfDNA analysis. We quantified un-/methylated copies of these DMRs by droplet digital PCR analysis in a plasma test set from TNBC patients and healthy controls and confirmed our findings obtained on tissues. Differential cfDNA methylation was confirmed in an independent validation set of plasma samples. A methylation score combining signatures of the top three DMRs overlapping with the SPAG6, LINC10606 and TBCD/ZNF750 genes had the best capability to discriminate TNBC patients from controls (AUC = 0.78 in the test set and AUC = 0.74 in validation set). Our findings demonstrate the usefulness of cfDNA-based methylation signatures as noninvasive liquid biopsy markers for the diagnosis of TNBC.

(Z)-Endoxifen and Early Recurrence of Breast Cancer: An Explorative Analysis in a Prospective Brazilian Study.

Adherence to treatment and use of co-medication, but also molecular factors such as CYP2D6 genotype, affect tamoxifen metabolism, with consequences for early breast cancer prognosis. In a prospective study of 149 tamoxifen-treated early-stage breast cancer patients from Brazil followed up for 5 years, we investigated the association between the active tamoxifen metabolite (Z)-endoxifen at 3 months and event-free survival (EFS) adjusted for clinico-pathological factors. Twenty-five patients (16.8%) had recurred or died at a median follow-up of 52.3 months. When we applied a putative 15 nM threshold used in previous independent studies, (Z)-endoxifen levels below the threshold showed an association with shorter EFS in univariate analysis (p = 0.045) and after adjustment for stage (HR 2.52; 95% CI 1.13-5.65; p = 0.024). However, modeling of plasma concentrations with splines instead of dichotomization did not verify a significant association with EFS (univariate analysis: p = 0.158; adjusted for stage: p = 0.117). Hence, in our small exploratory study, the link between impaired tamoxifen metabolism and early breast cancer recurrence could not be unanimously demonstrated. This inconsistency justifies larger modeling studies backed up by mechanistic pharmacodynamic analyses to shed new light on this suspected association and the stipulation of an appropriate predictive (Z)-endoxifen threshold.

Common variants in breast cancer risk loci predispose to distinct tumor subtypes.

BACKGROUND: Genome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER) status, but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear. METHODS: Among 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes. RESULTS: Eighty-five of 173 variants were associated with at least one tumor feature (false discovery rate < 5%), most commonly ER and grade, followed by PR and HER2. Models for intrinsic-like subtypes found nearly all of these variants (83 of 85) associated at p < 0.05 with risk for at least one luminal-like subtype, and approximately half (41 of 85) of the variants were associated with risk of at least one non-luminal subtype, including 32 variants associated with triple-negative (TN) disease. Ten variants were associated with risk of all subtypes in different magnitude. Five variants were associated with risk of luminal A-like and TN subtypes in opposite directions. CONCLUSION: This report demonstrates a high level of complexity in the etiology heterogeneity of breast cancer susceptibility variants and can inform investigations of subtype-specific risk prediction.

Subunits of BK channels promote breast cancer development and modulate responses to endocrine treatment in preclinical models.

BACKGROUND AND PURPOSE: Pore-forming α subunits of the voltage- and Ca2+ -activated K+ channel with large conductance (BKα) promote malignant phenotypes of breast tumour cells. Auxiliary subunits such as the leucine-rich repeat containing 26 (LRRC26) protein, also termed BKγ1, may be required to permit activation of BK currents at a depolarized resting membrane potential that frequently occur in non-excitable tumour cells. EXPERIMENTAL APPROACH: Anti-tumour effects of BKα loss were investigated in breast tumour-bearing MMTV-PyMT transgenic BKα knockout (KO) mice, primary MMTV-PyMT cell cultures, and in a syngeneic transplantation model of breast cancer derived from these cells. The therapeutic relevance of BK channels in the context of endocrine treatment was assessed in human breast cancer cell lines expressing either low (MCF-7) or high (MDA-MB-453) levels of BKα and BKγ1, as well as in BKα-negative MDA-MB-157. KEY RESULTS: BKα promoted breast cancer onset and overall survival in preclinical models. Conversely, lack of BKα and/or knockdown of BKγ1 attenuated proliferation of murine and human breast cancer cells in vitro. At low concentrations, tamoxifen and its major active metabolites stimulated proliferation of BKα/γ1-positive breast cancer cells, independent of the genomic signalling controlled by the oestrogen receptor. Finally, tamoxifen increased the relative survival time of BKα KO but not of wild-type tumour cell recipient mice. CONCLUSION AND IMPLICATIONS: Breast cancer initiation, progression, and tamoxifen sensitivity depend on functional BK channels thereby providing a rationale for the future exploration of the oncogenic actions of BK channels in clinical outcomes with anti-oestrogen therapy. LINKED ARTICLES: This article is part of a themed issue on New avenues in cancer prevention and treatment (BJP 75th Anniversary). To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v179.12/issuetoc.

Epigenetic quantification of circulating immune cells in peripheral blood of triple-negative breast cancer patients.

BACKGROUND: A shift in the proportions of blood immune cells is a hallmark of cancer development. Here, we investigated whether methylation-derived immune cell type ratios and methylation-derived neutrophil-to-lymphocyte ratios (mdNLRs) are associated with triple-negative breast cancer (TNBC). METHODS: Leukocyte subtype-specific unmethylated/methylated CpG sites were selected, and methylation levels at these sites were used as proxies for immune cell type proportions and mdNLR estimation in 231 TNBC cases and 231 age-matched controls. Data were validated using the Houseman deconvolution method. Additionally, the natural killer (NK) cell ratio was measured in a prospective sample set of 146 TNBC cases and 146 age-matched controls. RESULTS: The mdNLRs were higher in TNBC cases compared with controls and associated with TNBC (odds ratio (OR) range (2.66-4.29), all Padj. < 1e-04). A higher neutrophil ratio and lower ratios of NK cells, CD4 + T cells, CD8 + T cells, monocytes, and B cells were associated with TNBC. The strongest association was observed with decreased NK cell ratio (OR range (1.28-1.42), all Padj. < 1e-04). The NK cell ratio was also significantly lower in pre-diagnostic samples of TNBC cases compared with controls (P = 0.019). CONCLUSION: This immunomethylomic study shows that a shift in the ratios/proportions of leukocyte subtypes is associated with TNBC, with decreased NK cell showing the strongest association. These findings improve our knowledge of the role of the immune system in TNBC and point to the possibility of using NK cell level as a non-invasive molecular marker for TNBC risk assessment, early detection, and prevention.

Mendelian randomisation study of smoking exposure in relation to breast cancer risk.

BACKGROUND: Despite a modest association between tobacco smoking and breast cancer risk reported by recent epidemiological studies, it is still equivocal whether smoking is causally related to breast cancer risk. METHODS: We applied Mendelian randomisation (MR) to evaluate a potential causal effect of cigarette smoking on breast cancer risk. Both individual-level data as well as summary statistics for 164 single-nucleotide polymorphisms (SNPs) reported in genome-wide association studies of lifetime smoking index (LSI) or cigarette per day (CPD) were used to obtain MR effect estimates. Data from 108,420 invasive breast cancer cases and 87,681 controls were used for the LSI analysis and for the CPD analysis conducted among ever-smokers from 26,147 cancer cases and 26,072 controls. Sensitivity analyses were conducted to address pleiotropy. RESULTS: Genetically predicted LSI was associated with increased breast cancer risk (OR 1.18 per SD, 95% CI: 1.07-1.30, P = 0.11 × 10-2), but there was no evidence of association for genetically predicted CPD (OR 1.02, 95% CI: 0.78-1.19, P = 0.85). The sensitivity analyses yielded similar results and showed no strong evidence of pleiotropic effect. CONCLUSION: Our MR study provides supportive evidence for a potential causal association with breast cancer risk for lifetime smoking exposure but not cigarettes per day among smokers.

Author Correction: A case-only study to identify genetic modifiers of breast cancer risk for BRCA1/BRCA2 mutation carriers.

A Correction to this paper has been published: https://doi.org/10.1038/s41467-021-23162-4.

A case-only study to identify genetic modifiers of breast cancer risk for BRCA1/BRCA2 mutation carriers.

Breast cancer (BC) risk for BRCA1 and BRCA2 mutation carriers varies by genetic and familial factors. About 50 common variants have been shown to modify BC risk for mutation carriers. All but three, were identified in general population studies. Other mutation carrier-specific susceptibility variants may exist but studies of mutation carriers have so far been underpowered. We conduct a novel case-only genome-wide association study comparing genotype frequencies between 60,212 general population BC cases and 13,007 cases with BRCA1 or BRCA2 mutations. We identify robust novel associations for 2 variants with BC for BRCA1 and 3 for BRCA2 mutation carriers, P < 10-8, at 5 loci, which are not associated with risk in the general population. They include rs60882887 at 11p11.2 where MADD, SP11 and EIF1, genes previously implicated in BC biology, are predicted as potential targets. These findings will contribute towards customising BC polygenic risk scores for BRCA1 and BRCA2 mutation carriers.

CYP3A7*1C allele: linking premenopausal oestrone and progesterone levels with risk of hormone receptor-positive breast cancers.

BACKGROUND: Epidemiological studies provide strong evidence for a role of endogenous sex hormones in the aetiology of breast cancer. The aim of this analysis was to identify genetic variants that are associated with urinary sex-hormone levels and breast cancer risk. METHODS: We carried out a genome-wide association study of urinary oestrone-3-glucuronide and pregnanediol-3-glucuronide levels in 560 premenopausal women, with additional analysis of progesterone levels in 298 premenopausal women. To test for the association with breast cancer risk, we carried out follow-up genotyping in 90,916 cases and 89,893 controls from the Breast Cancer Association Consortium. All women were of European ancestry. RESULTS: For pregnanediol-3-glucuronide, there were no genome-wide significant associations; for oestrone-3-glucuronide, we identified a single peak mapping to the CYP3A locus, annotated by rs45446698. The minor rs45446698-C allele was associated with lower oestrone-3-glucuronide (-49.2%, 95% CI -56.1% to -41.1%, P = 3.1 × 10-18); in follow-up analyses, rs45446698-C was also associated with lower progesterone (-26.7%, 95% CI -39.4% to -11.6%, P = 0.001) and reduced risk of oestrogen and progesterone receptor-positive breast cancer (OR = 0.86, 95% CI 0.82-0.91, P = 6.9 × 10-8). CONCLUSIONS: The CYP3A7*1C allele is associated with reduced risk of hormone receptor-positive breast cancer possibly mediated via an effect on the metabolism of endogenous sex hormones in premenopausal women.

Breast Cancer Risk Genes - Association Analysis in More than 113,000 Women.

BACKGROUND: Genetic testing for breast cancer susceptibility is widely used, but for many genes, evidence of an association with breast cancer is weak, underlying risk estimates are imprecise, and reliable subtype-specific risk estimates are lacking. METHODS: We used a panel of 34 putative susceptibility genes to perform sequencing on samples from 60,466 women with breast cancer and 53,461 controls. In separate analyses for protein-truncating variants and rare missense variants in these genes, we estimated odds ratios for breast cancer overall and tumor subtypes. We evaluated missense-variant associations according to domain and classification of pathogenicity. RESULTS: Protein-truncating variants in 5 genes (ATM, BRCA1, BRCA2, CHEK2, and PALB2) were associated with a risk of breast cancer overall with a P value of less than 0.0001. Protein-truncating variants in 4 other genes (BARD1, RAD51C, RAD51D, and TP53) were associated with a risk of breast cancer overall with a P value of less than 0.05 and a Bayesian false-discovery probability of less than 0.05. For protein-truncating variants in 19 of the remaining 25 genes, the upper limit of the 95% confidence interval of the odds ratio for breast cancer overall was less than 2.0. For protein-truncating variants in ATM and CHEK2, odds ratios were higher for estrogen receptor (ER)-positive disease than for ER-negative disease; for protein-truncating variants in BARD1, BRCA1, BRCA2, PALB2, RAD51C, and RAD51D, odds ratios were higher for ER-negative disease than for ER-positive disease. Rare missense variants (in aggregate) in ATM, CHEK2, and TP53 were associated with a risk of breast cancer overall with a P value of less than 0.001. For BRCA1, BRCA2, and TP53, missense variants (in aggregate) that would be classified as pathogenic according to standard criteria were associated with a risk of breast cancer overall, with the risk being similar to that of protein-truncating variants. CONCLUSIONS: The results of this study define the genes that are most clinically useful for inclusion on panels for the prediction of breast cancer risk, as well as provide estimates of the risks associated with protein-truncating variants, to guide genetic counseling. (Funded by European Union Horizon 2020 programs and others.).

Combined Associations of a Polygenic Risk Score and Classical Risk Factors With Breast Cancer Risk.

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.