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.

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.

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).

Abstract P1-07-08: Endoxifen predicts early recurrence in breast cancer: A Brazilian prospective study

Tamoxifen is widely used in adjuvant therapy for breast cancer patients, but up to one-third of patients recur within 15 years after surgery even with its regular use. Low plasma levels of (Z)-endoxifen, the most active tamoxifen metabolite, have been studied as a possible predictor of treatment failure. Here, we evaluated plasma (Z)-endoxifen as a predictor of response to adjuvant treatment in breast cancer patients. Methods: One hundred forty nine patients were selected from a prospective cohort of consecutive newly diagnosed breast cancer patients recruited at a national referral center for treatment of cancer in Curitiba, Southern Brazil. Inclusion criteria were breast cancer women on stage I - III, positive hormone receptor tumors, operable, undergoing or not (neo)adjuvant chemotherapy, and signed for adjuvant tamoxifen for at least five years. At diagnosis, peripheral blood samples were collected for CYP2D6 genotyping for *1 *2, *3, *4, *5, *6, *7, *8, *9, *10, *17, *29 and *41 and subsequent metabolic activity scores according to phenotypes assigned for each allele (poor metabolizer PM; intermediate, IM; extensive, EM and ultrarapid metabolizer, UM). After the tamoxifen initiation, patients were evaluated at months 3, 6, and 12 for adherence and drug interactions, and plasma (Z)-endoxifen quantification. They were followed from April 2014 to April 2019. The study's primary objective was to evaluate whether (Z)-endoxifen plasma levels predict early recurrence in patients on adjuvant tamoxifen. The secondary aim was to evaluate an association between (Z)-endoxifen levels and clinical and phenotypic CYP2D6 metabolism variables. A Cox model estimated hazard ratios (HR) comparing recurrence free survival (RFS) according to (Z)-endoxifen plasma levels. The survival was compared using the Kaplan-Meier method and the log-rank test. Results: The median follow-up period was 52.3 months and at that point, there were 20 recurrences of 149 patients. At 3 months there were no cases of low adherence to tamoxifen. The (Z)-endoxifen levels were established by quartiles, with the first lower quartile being used as a cutoff (14.72nM). High (Z)-endoxifen levels at 3 months (&gt;14.72nM) were detected in 75.8% (n = 113) of cases and low levels (≤14.72nM) in 24.2% (n = 36). There was no statistically significant difference between the groups with low or high (Z)-endoxifen levels regarding clinical and pathological variables, but there was a higher CYP2D6 PM/IM frequency in the low (Z)-endoxifen levels group (p&lt;0.001). High (Z)-endoxifen levels were associated with lower RFS compared to low (Z)-endoxifen levels (HR 2.5, 95% CI 1.00-6.22, p=0.041). Low (Z)-endoxifen levels increased the recurrence risk even when adjusted for staging and molecular subtypes (HR 2.55, 95% CI 1.00-6.48, p&lt;0.001) or staging and CYP2D6 phenotype (UM/EM vs IM/PM; HR 4.38, 95% CI 1.66-11.59, p=0.003), or prior chemotherapy exposure (HR 2.65, 95% CI 1.06-6.61, p=0.036). Conclusion: This is the first Brazilian prospective trial evaluating the association of plasma (Z)-endoxifen levels, with the clinical outcome in breast cancer women treated with adjuvant tamoxifen, controlled by predictable confounding factors of plasma (Z)-endoxifen levels, such as tamoxifen adherence, and drug interaction with CYP2D6 inhibitors. In our study, a low plasma level of (Z)-endoxifen (≤ 14.72nM) after 3 months of tamoxifen treatment was a predictor of early recurrence in patients under adjuvant tamoxifen (p=0.041). A longer follow-up and a larger cohort may be needed to confirm the (Z)-endoxifen cutoff value that is a predictor of long-term recurrence.

Citation Format: Thais Almeida, Werner Schroth, Jeanine Marie Nardin, Thomas Mürdter, Solane Picolotto, Reiner Hoppe, Jenifer Primon Kogin, Elisa Daniele Gaio, Angela Dasenbrock, Raquel Cristina Skrsypcsak, Lucia Noronha, Hiltrud Brauch, José Claudio Casali-da-Rocha. Endoxifen predicts early recurrence in breast cancer: A Brazilian prospective study [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P1-07-08.

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.

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 P_adj. < 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 P_adj. < 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.

Association of germline genetic variants with breast cancer-specific survival in patient subgroups defined by clinic-pathological variables related to tumor biology and type of systemic treatment

BACKGROUND

Given the high heterogeneity among breast tumors, associations between common germline genetic variants and survival that may exist within specific subgroups could go undetected in an unstratified set of breast cancer patients.

METHODS

We performed genome-wide association analyses within 15 subgroups of breast cancer patients based on prognostic factors, including hormone receptors, tumor grade, age, and type of systemic treatment. Analyses were based on 91,686 female patients of European ancestry from the Breast Cancer Association Consortium, including 7531 breast cancer-specific deaths over a median follow-up of 8.1 years. Cox regression was used to assess associations of common germline variants with 15-year and 5-year breast cancer-specific survival. We assessed the probability of these associations being true positives via the Bayesian false discovery probability (BFDP < 0.15).

RESULTS

Evidence of associations with breast cancer-specific survival was observed in three patient subgroups, with variant rs5934618 in patients with grade 3 tumors (15-year-hazard ratio (HR) [95% confidence interval (CI)] 1.32 [1.20, 1.45], P = 1.4E-08, BFDP = 0.01, per G allele); variant rs4679741 in patients with ER-positive tumors treated with endocrine therapy (15-year-HR [95% CI] 1.18 [1.11, 1.26], P = 1.6E-07, BFDP = 0.09, per G allele); variants rs1106333 (15-year-HR [95% CI] 1.68 [1.39,2.03], P = 5.6E-08, BFDP = 0.12, per A allele) and rs78754389 (5-year-HR [95% CI] 1.79 [1.46,2.20], P = 1.7E-08, BFDP = 0.07, per A allele), in patients with ER-negative tumors treated with chemotherapy.

CONCLUSIONS

We found evidence of four loci associated with breast cancer-specific survival within three patient subgroups. There was limited evidence for the existence of associations in other patient subgroups. However, the power for many subgroups is limited due to the low number of events. Even so, our results suggest that the impact of common germline genetic variants on breast cancer-specific survival might be limited.

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

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.

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.

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.).

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 Polygenic Risk Score and Contralateral Breast Cancer Risk

Previous research has shown that polygenic risk scores (PRSs) can be used to stratify women according to their risk of developing primary invasive breast cancer. This study aimed to evaluate the association between a recently validated PRS of 313 germline variants (PRS313) and contralateral breast cancer (CBC) risk. We included 56,068 women of European ancestry diagnosed with first invasive breast cancer from 1990 onward with follow-up from the Breast Cancer Association Consortium. Metachronous CBC risk (N = 1,027) according to the distribution of PRS313 was quantified using Cox regression analyses. We assessed PRS313 interaction with age at first diagnosis, family history, morphology, ER status, PR status, and HER2 status, and (neo)adjuvant therapy. In studies of Asian women, with limited follow-up, CBC risk associated with PRS313 was assessed using logistic regression for 340 women with CBC compared with 12,133 women with unilateral breast cancer. Higher PRS313 was associated with increased CBC risk: hazard ratio per standard deviation (SD) = 1.25 (95%CI = 1.18-1.33) for Europeans, and an OR per SD = 1.15 (95%CI = 1.02-1.29) for Asians. The absolute lifetime risks of CBC, accounting for death as competing risk, were 12.4% for European women at the 10th percentile and 20.5% at the 90th percentile of PRS313. We found no evidence of confounding by or interaction with individual characteristics, characteristics of the primary tumor, or treatment. The C-index for the PRS313 alone was 0.563 (95%CI = 0.547-0.586). In conclusion, PRS313 is an independent factor associated with CBC risk and can be incorporated into CBC risk prediction models to help improve stratification and optimize surveillance and treatment strategies.

Downregulation of 15-hydroxyprostaglandin dehydrogenase during acquired tamoxifen resistance and association with poor prognosis in ERα-positive breast cancer

Aim:

Tamoxifen (TAM) resistance remains a clinical issue in breast cancer. The authors previously reported that 15-hydroxyprostaglandin dehydrogenase (HPGD) was significantly downregulated in tamoxifen-resistant (TAMr) breast cancer cell lines. Here, the authors investigated the relationship between HPGD expression, TAM resistance and prediction of outcome in breast cancer.

Methods:

HPGD overexpression and silencing studies were performed in isogenic TAMr and parental human breast cancer cell lines to establish the impact of HPGD expression on TAM resistance. HPGD expression and clinical outcome relationships were explored using immunohistochemistry and in silico analysis.

Results:

Restoration of HPGD expression and activity sensitised TAMr MCF-7 cells to TAM and 17β-oestradiol, whilst HPGD silencing in parental MCF-7 cells reduced TAM sensitivity. TAMr cells released more prostaglandin E₂ (PGE₂) than controls, which was reduced in TAMr cells stably transfected with HPGD. Exogenous PGE₂ signalled through the EP4 receptor to reduce breast cancer cell sensitivity to TAM. Decreased HPGD expression was associated with decreased overall survival in ERα-positive breast cancer patients.

Conclusions:

HPGD downregulation in breast cancer is associated with reduced response to TAM therapy via PGE₂-EP4 signalling and decreases patient survival. The data offer a potential target to develop combination therapies that may overcome acquired tamoxifen resistance.

Gene Expression Signatures of BRCAness and Tumor Inflammation Define Subgroups of Early-Stage Hormone Receptor-Positive Breast Cancer Patients

PURPOSE

Patients with estrogen receptor- and/or progesterone receptor-positive, early breast cancer benefit from hormonal treatment, yet high global death burdens due to high prevalence and long-term recurrence risk call for biomarkers to guide additional treatment approaches.

EXPERIMENTAL DESIGN

From a prospective, observational study of postmenopausal early breast cancer patients treated with tamoxifen or aromatase inhibitors, gene expression analyses of 612 tumors was performed using the NanoString Breast Cancer 360 panel to interrogate 23 breast cancer pathways. Candidate signatures associated with disease subtype and event-free survival (EFS) were obtained by cluster analysis, Cox modeling, and conditional inference trees, and were independently tested in 613 patients from BreastMark. Tumor-infiltrating lymphocytes (TIL) were assessed on tissue sections, and mutational burden was assessed in 36 tumors by whole-exome sequencing.

RESULTS

PAM50-derived classification distinguished lower-risk (Luminal A) from higher-risk subtypes (Luminal B, P = 0.04; HER2, P = 0.006; Basal, P = 0.008). In higher-risk patients, shorter EFS was associated with low androgen receptor [HR = 3.61; 95% confidence interval (CI), 1.72-7.56; P = 0.001] or high BRCAness signature expression (HR = 3.58; 95% CI, 1.19-10.7; P = 0.023). BRCAness was independently confirmed as a predictor of shorter EFS (HR = 2.64; 95% CI, 1.31-5.34; P = 0.007). About 13%-15% of patients, enriched for high-grade, higher-risk subtypes (P ≤ 0.0001), had strong expression of the Tumor Inflammation Signature (TIS) suggestive of an inhibited antitumor immune response. TIS scores were strongly associated with TIL numbers (P < 1e-30) but not with tumor mutation status.

CONCLUSIONS

BRCA-related DNA repair deficiency and suppressed tumor immune responses may be clinically relevant predictors of endocrine therapy complementing treatment options in subgroups of hormone-sensitive early breast cancer.

Obesity Alters Endoxifen Plasma Levels in Young Breast Cancer Patients: A Pharmacometric Simulation Approach

Endoxifen is one of the most important metabolites of the prodrug tamoxifen. High interindividual variability in endoxifen steady-state concentrations (CSS,min ENDX ) is observed under tamoxifen standard dosing and patients with breast cancer who do not reach endoxifen concentrations above a proposed therapeutic threshold of 5.97 ng/mL may be at a 26% higher recurrence risk compared with patients with endoxifen concentrations exceeding this value. In this investigation, 10 clinical tamoxifen studies were pooled (1,388 patients) to investigate influential factors on CSS,min ENDX using nonlinear mixed-effects modeling. Age and body weight were found to significantly impact CSS,min ENDX in addition to CYP2D6 phenotype. Compared with postmenopausal patients, premenopausal patients had a 30% higher risk for subtarget CSS,min ENDX at tamoxifen 20 mg per day. In treatment simulations for distinct patient subpopulations, young overweight patients had a 3.1-13.8-fold higher risk for subtarget CSS,min ENDX compared with elderly low-weight patients. Considering ever-rising obesity rates and the clinical importance of tamoxifen for premenopausal patients, this subpopulation may benefit most from individualized tamoxifen dosing.

Transcriptome-wide association study of breast cancer risk by estrogen-receptor status

Previous 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.

Germline HOXB13 mutations p.G84E and p.R217C do not confer an increased breast cancer risk

In breast cancer, high levels of homeobox protein Hox-B13 (HOXB13) have been associated with disease progression of ER-positive breast cancer patients and resistance to tamoxifen treatment. Since HOXB13 p.G84E is a prostate cancer risk allele, we evaluated the association between HOXB13 germline mutations and breast cancer risk in a previous study consisting of 3,270 familial non-BRCA1/2 breast cancer cases and 2,327 controls from the Netherlands. Although both recurrent HOXB13 mutations p.G84E and p.R217C were not associated with breast cancer risk, the risk estimation for p.R217C was not very precise. To provide more conclusive evidence regarding the role of HOXB13 in breast cancer susceptibility, we here evaluated the association between HOXB13 mutations and increased breast cancer risk within 81 studies of the international Breast Cancer Association Consortium containing 68,521 invasive breast cancer patients and 54,865 controls. Both HOXB13 p.G84E and p.R217C did not associate with the development of breast cancer in European women, neither in the overall analysis (OR = 1.035, 95% CI = 0.859-1.246, P = 0.718 and OR = 0.798, 95% CI = 0.482-1.322, P = 0.381 respectively), nor in specific high-risk subgroups or breast cancer subtypes. Thus, although involved in breast cancer progression, HOXB13 is not a material breast cancer susceptibility gene.

Genome-wide association study identifies 32 novel breast cancer susceptibility loci from overall and subtype-specific analyses

Breast cancer susceptibility variants frequently show heterogeneity in associations by tumor subtype1-3. To identify novel loci, we performed a genome-wide association study including 133,384 breast cancer cases and 113,789 controls, plus 18,908 BRCA1 mutation carriers (9,414 with breast cancer) of European ancestry, using both standard and novel methodologies that account for underlying tumor heterogeneity by estrogen receptor, progesterone receptor and human epidermal growth factor receptor 2 status and tumor grade. We identified 32 novel susceptibility loci (P < 5.0 × 10-8), 15 of which showed evidence for associations with at least one tumor feature (false discovery rate < 0.05). Five loci showed associations (P < 0.05) in opposite directions between luminal and non-luminal subtypes. In silico analyses showed that these five loci contained cell-specific enhancers that differed between normal luminal and basal mammary cells. The genetic correlations between five intrinsic-like subtypes ranged from 0.35 to 0.80. The proportion of genome-wide chip heritability explained by all known susceptibility loci was 54.2% for luminal A-like disease and 37.6% for triple-negative disease. The odds ratios of polygenic risk scores, which included 330 variants, for the highest 1% of quantiles compared with middle quantiles were 5.63 and 3.02 for luminal A-like and triple-negative disease, respectively. These findings provide an improved understanding of genetic predisposition to breast cancer subtypes and will inform the development of subtype-specific polygenic risk scores.

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 Ca²⁺ -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.

Abstract P2-11-12: The influences of adherence to tamoxifen and CYP2D6 pharmacogenetics on plasma concentrations of the active metabolite (Z)-endoxifen in breast cancer

One-third of early breast cancer patients with estrogen receptor (ER) positive tumors treated with the selective estrogen receptor modulator tamoxifen either relapse or die from the disease in the following decade. Improvement of tamoxifen efficacy requires a better knowledge of factors determining outcome of which the medication-taking behavior, i.e. adherence, is a strong suspect. Discontinuation and non-adherence to tamoxifen are frequent and result in increased mortality. Commonly used methods to assess drug adherence such as prescription-refill patterns infer that the medication is taken exactly as prescribed yet, partial adherence or actual use of the medication cannot be controlled. An objective surrogate of tamoxifen adherence is the monitoring of drug concentrations in the patients’ plasma including the parent drug and its active metabolite (Z)-endoxifen (endoxifen), however individual concentrations depend on pharmacogene polymorphisms particularly those of the liver enzyme cytochrome P450 2D6 (CYP2D6). Patients without CYP2D6 variants have normal enzyme function (extensive metabolizer, EM) and achieve high endoxifen concentrations compared to patients with variants conferring absent activity (poor metabolizer, PM) or reduced activity (intermediate metabolizer, IM) who have significantly lower plasma concentrations of the active metabolite. A putative threshold of 5.9 ng/mL endoxifen has been proposed suggesting that patients with endoxifen concentration above or below threshold have lower or higher risk for recurrence, assuming that low endoxifen exposure in the tumor likely results in incomplete inhibition of ER-dependent growth signaling. During tamoxifen treatment it is therefore important that patients achieve clinically relevant endoxifen levels while at the same time avoid confounders of which adherence is a prime candidate. To this end, we investigated patients’ ability to achieve relevant plasma concentrations during tamoxifen treatment in relation to their CYP2D6 metabolizer status and their tamoxifen adherence behavior. Our prospective study included 192 Brazilian breast cancer patients treated with 20 mg tamoxifen daily. Two-thirds of the patients were premenopausal. Plasma levels of tamoxifen and its metabolites were measured by LC-MS/MS at 3, 6 and 12 months of treatment. Adherence to tamoxifen was evaluated by the questionnaire-based Morisky medication adherence scale. CYP2D6 polymorphisms were determined by MALDI-TOF mass spectrometry and real-time PCR. Adherence to tamoxifen decreased from 74% at month 3 to 64% at month 12, and explained 47% of the variability of tamoxifen plasma concentrations (P&lt;.001). Accordingly, plasma monitoring of tamoxifen is useful to inform on tamoxifen adherence. Endoxifen concentrations were associated with CYP2D6 status and adherence. While CYP2D6 explained 26.4% of endoxifen variability at 12 months, the combination with adherence explained 40% (P&lt;.001). The influence of low adherence to not achieving therapeutic endoxifen levels was highest in patients with non-compromised CYP2D6 function (EM phenotype, RR 3.65; 95% CI: 1.48-8.99). we conclude that for patients with intact CYP2D6 function adherence to tamoxifen is critical to achieve therapeutic endoxifen plasma levels. CYP2D6 phenotyping and endoxifen plasma monitoring may help to guide EM (and IM) patients to adhere to tamoxifen which is particularly relevant for premenopausal women who are at higher risk to not adhere to tamoxifen.

Citation Format: Jose Claudio Casali da Rocha, Jeanine M Nardin, Werner Schroth, Thais A Almeida, Thomas Mürdter, Solange Picolotto, Evelyn CL Vendramini, Reiner Hoppe, Jenifer P Kogin, Diandra Miqueleto, Silvia DR Moraes, Matthias Schwab, Roberto F Pecoits-Filho, Hiltrud Brauch. The influences of adherence to tamoxifen and CYP2D6 pharmacogenetics on plasma concentrations of the active metabolite (Z)-endoxifen in breast cancer [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr P2-11-12.

A network analysis to identify mediators of germline-driven differences in breast cancer prognosis

Identifying the underlying genetic drivers of the heritability of breast cancer prognosis remains elusive. We adapt a network-based approach to handle underpowered complex datasets to provide new insights into the potential function of germline variants in breast cancer prognosis. This network-based analysis studies ~7.3 million variants in 84,457 breast cancer patients in relation to breast cancer survival and confirms the results on 12,381 independent patients. Aggregating the prognostic effects of genetic variants across multiple genes, we identify four gene modules associated with survival in estrogen receptor (ER)-negative and one in ER-positive disease. The modules show biological enrichment for cancer-related processes such as G-alpha signaling, circadian clock, angiogenesis, and Rho-GTPases in apoptosis.

Fine-mapping of 150 breast cancer risk regions identifies 191 likely target genes

Genome-wide association studies have identified breast cancer risk variants in over 150 genomic regions, but the mechanisms underlying risk remain largely unknown. These regions were explored by combining association analysis with in silico genomic feature annotations. We defined 205 independent risk-associated signals with the set of credible causal variants in each one. In parallel, we used a Bayesian approach (PAINTOR) that combines genetic association, linkage disequilibrium and enriched genomic features to determine variants with high posterior probabilities of being causal. Potentially causal variants were significantly over-represented in active gene regulatory regions and transcription factor binding sites. We applied our INQUSIT pipeline for prioritizing genes as targets of those potentially causal variants, using gene expression (expression quantitative trait loci), chromatin interaction and functional annotations. Known cancer drivers, transcription factors and genes in the developmental, apoptosis, immune system and DNA integrity checkpoint gene ontology pathways were over-represented among the highest-confidence target genes.

The Influences of Adherence to Tamoxifen and CYP2D6 Pharmacogenetics on Plasma Concentrations of the Active Metabolite (Z)-Endoxifen in Breast Cancer

Tamoxifen efficacy in breast cancer is suspected to depend on adherence and intact drug metabolism. We evaluated the role of adherence behavior and pharmacogenetics on the formation rate of (Z)-endoxifen. In 192 Brazilian patients, we assessed plasma levels of tamoxifen and its metabolites at 3, 6, and 12 months of treatment (liquid-chromatography tandem mass spectrometry), adherence behavior (Morisky, Green, and Levine medication adherence scale), and cytochrome P450 2D6 (CYP2D6) and other pharmacogene polymorphisms (matrix-assisted laser-desorption-ionization time of flight) mass spectrometry, real-time polymerase chain reaction). Adherence explained 47% of the variability of tamoxifen plasma concentrations (P < 0.001). Although CYP2D6 alone explained 26.4%, the combination with adherence explained 40% of (Z)-endoxifen variability at 12 months (P < 0.001). The influence of low adherence to not achieving relevant (Z)-endoxifen levels was highest in patients with noncompromised CYP2D6 function (relative risk 3.65; 95% confidence interval 1.48-8.99). As a proof-of-concept, we demonstrated that (Z)-endoxifen levels are influenced both by patient adherence to tamoxifen and CYP2D6, which is particularly relevant for patients with full CYP2D6 function.

The FANCM:p.Arg658* truncating variant is associated with risk of triple-negative breast cancer

Breast cancer is a common disease partially caused by genetic risk factors. Germline pathogenic variants in DNA repair genes BRCA1, BRCA2, PALB2, ATM, and CHEK2 are associated with breast cancer risk. FANCM, which encodes for a DNA translocase, has been proposed as a breast cancer predisposition gene, with greater effects for the ER-negative and triple-negative breast cancer (TNBC) subtypes. We tested the three recurrent protein-truncating variants FANCM:p.Arg658*, p.Gln1701*, and p.Arg1931* for association with breast cancer risk in 67,112 cases, 53,766 controls, and 26,662 carriers of pathogenic variants of BRCA1 or BRCA2. These three variants were also studied functionally by measuring survival and chromosome fragility in FANCM -/- patient-derived immortalized fibroblasts treated with diepoxybutane or olaparib. We observed that FANCM:p.Arg658* was associated with increased risk of ER-negative disease and TNBC (OR = 2.44, P = 0.034 and OR = 3.79; P = 0.009, respectively). In a country-restricted analysis, we confirmed the associations detected for FANCM:p.Arg658* and found that also FANCM:p.Arg1931* was associated with ER-negative breast cancer risk (OR = 1.96; P = 0.006). The functional results indicated that all three variants were deleterious affecting cell survival and chromosome stability with FANCM:p.Arg658* causing more severe phenotypes. In conclusion, we confirmed that the two rare FANCM deleterious variants p.Arg658* and p.Arg1931* are risk factors for ER-negative and TNBC subtypes. Overall our data suggest that the effect of truncating variants on breast cancer risk may depend on their position in the gene. Cell sensitivity to olaparib exposure, identifies a possible therapeutic option to treat FANCM-associated tumors.

Publisher Correction: Shared heritability and functional enrichment across six solid cancers

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Two truncating variants in FANCC and breast cancer risk

Fanconi anemia (FA) is a genetically heterogeneous disorder with 22 disease-causing genes reported to date. In some FA genes, monoallelic mutations have been found to be associated with breast cancer risk, while the risk associations of others remain unknown. The gene for FA type C, FANCC, has been proposed as a breast cancer susceptibility gene based on epidemiological and sequencing studies. We used the Oncoarray project to genotype two truncating FANCC variants (p.R185X and p.R548X) in 64,760 breast cancer cases and 49,793 controls of European descent. FANCC mutations were observed in 25 cases (14 with p.R185X, 11 with p.R548X) and 26 controls (18 with p.R185X, 8 with p.R548X). There was no evidence of an association with the risk of breast cancer, neither overall (odds ratio 0.77, 95%CI 0.44-1.33, p = 0.4) nor by histology, hormone receptor status, age or family history. We conclude that the breast cancer risk association of these two FANCC variants, if any, is much smaller than for BRCA1, BRCA2 or PALB2 mutations. If this applies to all truncating variants in FANCC it would suggest there are differences between FA genes in their roles on breast cancer risk and demonstrates the merit of large consortia for clarifying risk associations of rare variants.

Genome-wide association and transcriptome studies identify target genes and risk loci for breast cancer

Genome-wide association studies (GWAS) have identified more than 170 breast cancer susceptibility loci. Here we hypothesize that some risk-associated variants might act in non-breast tissues, specifically adipose tissue and immune cells from blood and spleen. Using expression quantitative trait loci (eQTL) reported in these tissues, we identify 26 previously unreported, likely target genes of overall breast cancer risk variants, and 17 for estrogen receptor (ER)-negative breast cancer, several with a known immune function. We determine the directional effect of gene expression on disease risk measured based on single and multiple eQTL. In addition, using a gene-based test of association that considers eQTL from multiple tissues, we identify seven (and four) regions with variants associated with overall (and ER-negative) breast cancer risk, which were not reported in previous GWAS. Further investigation of the function of the implicated genes in breast and immune cells may provide insights into the etiology of breast cancer.

Genome-wide association study of germline variants and breast cancer-specific mortality

BACKGROUND

We examined the associations between germline variants and breast cancer mortality using a large meta-analysis of women of European ancestry.

METHODS

Meta-analyses included summary estimates based on Cox models of twelve datasets using ~10.4 million variants for 96,661 women with breast cancer and 7697 events (breast cancer-specific deaths). Oestrogen receptor (ER)-specific analyses were based on 64,171 ER-positive (4116) and 16,172 ER-negative (2125) patients. We evaluated the probability of a signal to be a true positive using the Bayesian false discovery probability (BFDP).

RESULTS

We did not find any variant associated with breast cancer-specific mortality at P < 5 × 10-8. For ER-positive disease, the most significantly associated variant was chr7:rs4717568 (BFDP = 7%, P = 1.28 × 10-7, hazard ratio [HR] = 0.88, 95% confidence interval [CI] = 0.84-0.92); the closest gene is AUTS2. For ER-negative disease, the most significant variant was chr7:rs67918676 (BFDP = 11%, P = 1.38 × 10-7, HR = 1.27, 95% CI = 1.16-1.39); located within a long intergenic non-coding RNA gene (AC004009.3), close to the HOXA gene cluster.

CONCLUSIONS

We uncovered germline variants on chromosome 7 at BFDP < 15% close to genes for which there is biological evidence related to breast cancer outcome. However, the paucity of variants associated with mortality at genome-wide significance underpins the challenge in providing genetic-based individualised prognostic information for breast cancer patients.

Shared heritability and functional enrichment across six solid cancers

Quantifying the genetic correlation between cancers can provide important insights into the mechanisms driving cancer etiology. Using genome-wide association study summary statistics across six cancer types based on a total of 296,215 cases and 301,319 controls of European ancestry, here we estimate the pair-wise genetic correlations between breast, colorectal, head/neck, lung, ovary and prostate cancer, and between cancers and 38 other diseases. We observed statistically significant genetic correlations between lung and head/neck cancer (rg = 0.57, p = 4.6 × 10-8), breast and ovarian cancer (rg = 0.24, p = 7 × 10-5), breast and lung cancer (rg = 0.18, p =1.5 × 10-6) and breast and colorectal cancer (rg = 0.15, p = 1.1 × 10-4). We also found that multiple cancers are genetically correlated with non-cancer traits including smoking, psychiatric diseases and metabolic characteristics. Functional enrichment analysis revealed a significant excess contribution of conserved and regulatory regions to cancer heritability. Our comprehensive analysis of cross-cancer heritability suggests that solid tumors arising across tissues share in part a common germline genetic basis.

Polygenic Risk Scores for Prediction of Breast Cancer and Breast Cancer Subtypes

Stratification of women according to their risk of breast cancer based on polygenic risk scores (PRSs) could improve screening and prevention strategies. Our aim was to develop PRSs, optimized for prediction of estrogen receptor (ER)-specific disease, from the largest available genome-wide association dataset and to empirically validate the PRSs in prospective studies. The development dataset comprised 94,075 case subjects and 75,017 control subjects of European ancestry from 69 studies, divided into training and validation sets. Samples were genotyped using genome-wide arrays, and single-nucleotide polymorphisms (SNPs) were selected by stepwise regression or lasso penalized regression. The best performing PRSs were validated in an independent test set comprising 11,428 case subjects and 18,323 control subjects from 10 prospective studies and 190,040 women from UK Biobank (3,215 incident breast cancers). For the best PRSs (313 SNPs), the odds ratio for overall disease per 1 standard deviation in ten prospective studies was 1.61 (95%CI: 1.57-1.65) with area under receiver-operator curve (AUC) = 0.630 (95%CI: 0.628-0.651). The lifetime risk of overall breast cancer in the top centile of the PRSs was 32.6%. Compared with women in the middle quintile, those in the highest 1% of risk had 4.37- and 2.78-fold risks, and those in the lowest 1% of risk had 0.16- and 0.27-fold risks, of developing ER-positive and ER-negative disease, respectively. Goodness-of-fit tests indicated that this PRS was well calibrated and predicts disease risk accurately in the tails of the distribution. This PRS is a powerful and reliable predictor of breast cancer risk that may improve breast cancer prevention programs.

Long intergenic noncoding RNA 299 methylation in peripheral blood is a biomarker for triple-negative breast cancer

Aim: To identify DNA methylation biomarkers in peripheral blood samples from triple-negative breast cancer (TNBC) patients. Materials & methods: We conducted an epigenome-wide association study (EWAS): the most promising markers were identified in 233 TNBC case–control pairs (discovery set) and subsequently validated in an independent validation set (57 TNBC patients and 124 controls). Results: cg06588802 (LINC00299/ID2) showed a higher methylation in TNBC patients compared with controls (discovery set: 3% increase, p-value = 0.0009; validation set: 2% increase, p-value = 0.01). Consistent results at four neighboring methylation probes and the strong negative correlation (rho = -0.93) with LINC00299 expression add plausibility to this result. Conclusion: Hypermethylation of LINC00299 in peripheral blood may constitute a useful circulating biomarker for TNBC.

Identification of nine new susceptibility loci for endometrial cancer

Endometrial cancer is the most commonly diagnosed cancer of the female reproductive tract in developed countries. Through genome-wide association studies (GWAS), we have previously identified eight risk loci for endometrial cancer. Here, we present an expanded meta-analysis of 12,906 endometrial cancer cases and 108,979 controls (including new genotype data for 5624 cases) and identify nine novel genome-wide significant loci, including a locus on 12q24.12 previously identified by meta-GWAS of endometrial and colorectal cancer. At five loci, expression quantitative trait locus (eQTL) analyses identify candidate causal genes; risk alleles at two of these loci associate with decreased expression of genes, which encode negative regulators of oncogenic signal transduction proteins (SH2B3 (12q24.12) and NF1 (17q11.2)). In summary, this study has doubled the number of known endometrial cancer risk loci and revealed candidate causal genes for future study.

Triple-Negative Breast Cancer Risk Genes Identified by Multigene Hereditary Cancer Panel Testing

BACKGROUND

Germline genetic testing with hereditary cancer gene panels can identify women at increased risk of breast cancer. However, those at increased risk of triple-negative (estrogen receptor-negative, progesterone receptor-negative, human epidermal growth factor receptor-negative) breast cancer (TNBC) cannot be identified because predisposition genes for TNBC, other than BRCA1, have not been established. The aim of this study was to define the cancer panel genes associated with increased risk of TNBC.

METHODS

Multigene panel testing for 21 genes in 8753 TNBC patients was performed by a clinical testing laboratory, and testing for 17 genes in 2148 patients was conducted by a Triple Negative Breast Cancer Consortium (TNBCC) of research studies. Associations between deleterious mutations in cancer predisposition genes and TNBC were evaluated using results from TNBC patients and reference controls.

RESULTS

Germline pathogenic variants in BARD1, BRCA1, BRCA2, PALB2, and RAD51D were associated with high risk (odds ratio > 5.0) of TNBC and greater than 20% lifetime risk for overall breast cancer among Caucasians. Pathogenic variants in BRIP1, RAD51C, and TP53 were associated with moderate risk (odds ratio > 2) of TNBC. Similar trends were observed for the African American population. Pathogenic variants in these TNBC genes were detected in 12.0% (3.7% non-BRCA1/2) of all participants.

CONCLUSIONS

Multigene hereditary cancer panel testing can identify women with elevated risk of TNBC due to mutations in BARD1, BRCA1, BRCA2, PALB2, and RAD51D. These women can potentially benefit from improved screening, risk management, and cancer prevention strategies. Patients with mutations may also benefit from specific targeted therapeutic strategies.

Abstract 3425: Exome-sequencing derived mutations of endocrine treated ER-positive early breast cancer

Anti-hormonal treatment of ER-positive early breast cancer with tamoxifen and/or aromatase inhibitors (AI) significantly improves patient outcome, however there is a lack of markers that predict the one third of treated patients who relapse or die. Although ER positive (luminal) tumors are known to be genetically highly diverse, the relevance of somatic mutations with regards to endocrine treatment outcome and course of the patients' disease is not fully understood. We conducted a prospective observational clinical trial for the investigation of treatment outcome determinants including 1219 tamoxifen and/or AI-treated postmenopausal ER-positive early breast cancer patients (DRKS00000605, recruitment between 2005-2010, median follow up 5.1 years). Of these, 40 tumors divided into 20 patients each with and without disease recurrence matched by prognostic variables (age, tumor size, grade, nodal status, treatment) were analyzed for their mutation pattern. We performed whole-exome sequencing of tumor and normal DNA (Agilent SureSelect, Illumina) for the identification of single nucleotide variants and small indels. Somatic mutations were identified following subtraction of germline single nucleotide polymorphisms (SNP), filtering based on raw-read quality parameters, in silico predicted biological function, and disease relatedness using the Ingenuity variant analysis software package (Qiagen). Clinical outcome associations were tested either on gene or on variant level. Somatic mutations with predicted deleterious effects were identified in 18 out of 40 known mutation-driver genes including PIK3CA, GATA3, MAP3K1, and TP53 (45%). Thirteen patients carried a mutation in the most frequently affected PIK3CA (32.5%). Of note, three out of four patients with a E542K mutation belonged to the breast cancer recurrence group, while 5 of 6 tumors with the H1047R mutation originated from patients without recurrence suggesting differential effects of these PIK3CA hotspot mutations. Common germline SNPs at IL17RB, ANK2, POM121, PGR, PTPRB, RAD51D, RFPLA4, and SUMO3 provided first hints for a possible association with relapse status (p&lt;0.05). We will report on the extended genotype-outcome-correlations to further evaluate the relevance of somatic mutations for endocrine treatment outcome of early ER-positive breast cancer.

Citation Format: Werner Schroth, Siarhei Kandabarau, Peter Fritz, Thomas Mürdter, Liza Bacchus, German Tamoxifen Study Group, Matthias Schwab, Hiltrud Brauch. Exome-sequencing derived mutations of endocrine treated ER-positive early 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 3425.

A transcriptome-wide association study of 229,000 women identifies new candidate susceptibility genes for breast cancer

The breast cancer risk variants identified in genome-wide association studies explain only a small fraction of the familial relative risk, and the genes responsible for these associations remain largely unknown. To identify novel risk loci and likely causal genes, we performed a transcriptome-wide association study evaluating associations of genetically predicted gene expression with breast cancer risk in 122,977 cases and 105,974 controls of European ancestry. We used data from the Genotype-Tissue Expression Project to establish genetic models to predict gene expression in breast tissue and evaluated model performance using data from The Cancer Genome Atlas. Of the 8,597 genes evaluated, significant associations were identified for 48 at a Bonferroni-corrected threshold of P < 5.82 × 10-6, including 14 genes at loci not yet reported for breast cancer. We silenced 13 genes and showed an effect for 11 on cell proliferation and/or colony-forming efficiency. Our study provides new insights into breast cancer genetics and biology.

Genetic overlap between endometriosis and endometrial cancer: evidence from cross-disease genetic correlation and GWAS meta-analyses

Epidemiological, biological, and molecular data suggest links between endometriosis and endometrial cancer, with recent epidemiological studies providing evidence for an association between a previous diagnosis of endometriosis and risk of endometrial cancer. We used genetic data as an alternative approach to investigate shared biological etiology of these two diseases. Genetic correlation analysis of summary level statistics from genomewide association studies (GWAS) using LD Score regression revealed moderate but significant genetic correlation (rg  = 0.23, P = 9.3 × 10-3 ), and SNP effect concordance analysis provided evidence for significant SNP pleiotropy (P = 6.0 × 10-3 ) and concordance in effect direction (P = 2.0 × 10-3 ) between the two diseases. Cross-disease GWAS meta-analysis highlighted 13 distinct loci associated at P ≤ 10-5 with both endometriosis and endometrial cancer, with one locus (SNP rs2475335) located within PTPRD associated at a genomewide significant level (P = 4.9 × 10-8 , OR = 1.11, 95% CI = 1.07-1.15). PTPRD acts in the STAT3 pathway, which has been implicated in both endometriosis and endometrial cancer. This study demonstrates the value of cross-disease genetic analysis to support epidemiological observations and to identify biological pathways of relevance to multiple diseases.

The BRCA2 c.68-7T > A variant is not pathogenic: A model for clinical calibration of spliceogenicity

Although the spliceogenic nature of the BRCA2 c.68-7T > A variant has been demonstrated, its association with cancer risk remains controversial. In this study, we accurately quantified by real-time PCR and digital PCR (dPCR), the BRCA2 isoforms retaining or missing exon 3. In addition, the combined odds ratio for causality of the variant was estimated using genetic and clinical data, and its associated cancer risk was estimated by case-control analysis in 83,636 individuals. Co-occurrence in trans with pathogenic BRCA2 variants was assessed in 5,382 families. Exon 3 exclusion rate was 4.5-fold higher in variant carriers (13%) than controls (3%), indicating an exclusion rate for the c.68-7T > A allele of approximately 20%. The posterior probability of pathogenicity was 7.44 × 10-115 . There was neither evidence for increased risk of breast cancer (OR 1.03; 95% CI 0.86-1.24) nor for a deleterious effect of the variant when co-occurring with pathogenic variants. Our data provide for the first time robust evidence of the nonpathogenicity of the BRCA2 c.68-7T > A. Genetic and quantitative transcript analyses together inform the threshold for the ratio between functional and altered BRCA2 isoforms compatible with normal cell function. These findings might be exploited to assess the relevance for cancer risk of other BRCA2 spliceogenic variants.

Publisher Correction: Evidence that breast cancer risk at the 2q35 locus is mediated through IGFBP5 regulation

This corrects the article DOI: 10.1038/ncomms5999.

Joint associations of a polygenic risk score and environmental risk factors for breast cancer in the Breast Cancer Association Consortium

Background

Polygenic risk scores (PRS) for breast cancer can be used to stratify the population into groups at substantially different levels of risk. Combining PRS and environmental risk factors will improve risk prediction; however, integrating PRS into risk prediction models requires evaluation of their joint association with known environmental risk factors.

Methods

Analyses were based on data from 20 studies; datasets analysed ranged from 3453 to 23 104 invasive breast cancer cases and similar numbers of controls, depending on the analysed environmental risk factor. We evaluated joint associations of a 77-single nucleotide polymorphism (SNP) PRS with reproductive history, alcohol consumption, menopausal hormone therapy (MHT), height and body mass index (BMI). We tested the null hypothesis of multiplicative joint associations for PRS and each of the environmental factors, and performed global and tail-based goodness-of-fit tests in logistic regression models. The outcomes were breast cancer overall and by estrogen receptor (ER) status.

Results

The strongest evidence for a non-multiplicative joint associations with the 77-SNP PRS was for alcohol consumption (P-interaction = 0.009), adult height (P-interaction = 0.025) and current use of combined MHT (P-interaction = 0.038) in ER-positive disease. Risk associations for these factors by percentiles of PRS did not follow a clear dose-response. In addition, global and tail-based goodness of fit tests showed little evidence for departures from a multiplicative risk model, with alcohol consumption showing the strongest evidence for ER-positive disease (P = 0.013 for global and 0.18 for tail-based tests).

Conclusions

The combined effects of the 77-SNP PRS and environmental risk factors for breast cancer are generally well described by a multiplicative model. Larger studies are required to confirm possible departures from the multiplicative model for individual risk factors, and assess models specific for ER-negative disease.

Is Mistletoe Treatment Beneficial in Invasive Breast Cancer? A New Approach to an Unresolved Problem

BACKGROUND/AIM

In this retrospective study, we compared breast cancer patients treated with and without mistletoe lectin I (ML-I) in addition to standard breast cancer treatment in order to determine a possible effect of this complementary treatment.

PATIENTS AND METHODS

This study included 18,528 patients with invasive breast cancer. Data on additional ML-I treatments were reported for 164 patients. We developed a "similar case" method with a distance measure retrieved from the beta variable in Cox regression to compare these patients, after stage adjustment, with their non-ML-1 treated counterparts in order to answer three hypotheses concerning overall survival, recurrence free survival and life quality.

RESULTS

Raw data analysis of an additional ML-I treatment yielded a worse outcome (p=0.02) for patients with ML treatment, possibly due to a bias inherent in the ML-I-treated patients. Using the "similar case" method (a case-based reasoning approach) we could not confirm this harm for patients using ML-I. Analysis of life quality data did not demonstrate reliable differences between patients treated with ML-I treatment and those without proven ML-I treatment.

CONCLUSION

Based on a "similar case" model we did not observe any differences in the overall survival (OS), recurrence-free survival (RFS), and quality of life data between breast cancer patients with standard treatment and those who in addition to standard treatment received ML-I treatment.

Profiles of miRNAs matched to biology in aromatase inhibitor resistant breast cancer

Aromatase inhibitor (AI) resistance during breast cancer treatment is mimicked by MCF-7:5C (5C) and MCF-7:2A (2A) cell lines that grow spontaneously. Survival signaling is reconfigured but cells are vulnerable to estradiol (E₂)-inducible apoptosis. These model systems have alterations of stress related pathways including the accumulation of endoplasmic reticulum, oxidative, and inflammatory stress that occur prior to E₂-induced apoptosis. We investigated miRNA expression profiles of 5C and 2A to characterize their AI resistance phenotypes. Affymetrix GeneChip miRNA2.0 arrays identified 184 miRNAs differentially expressed in 2A and 5C compared to E₂-free wild-type MCF-7:WS8. In 5C, 34 miRNAs of the DLK1-DIO3 locus and miR-31 were overexpressed, whereas miR-222 was low. TCGA data revealed poor and favorable overall survival for low miR-31 and miR-222 levels, respectively (HR=3.0, 95% CI:1.9-4.8; HR=0.3, 95% CI:0.1-0.6). Targets of deregulated miRNAs were identified using CLIP-confirmed TargetScan predictions. KEGG enrichment analyses for 5C- and 2A-specific target gene sets revealed pathways associated with cell proliferation including insulin, mTOR, and ErbB signaling as well as immune response and metabolism. Key genes overrepresented in 5C- and 2A-specific pathway interaction networks including EGFR, IGF1R and PIK3R1 had lower protein levels in 5C compared to 2A and were found to be differentially modulated by respective miRNA sets. Distinct up-regulated miRNAs from the DLK1-DIO3 locus may cause these attenuative effects as they are predicted to interact with corresponding 3′ untranslated regions. These new miRNA profiles become an important regulatory database to explore E₂-induced apoptotic mechanisms of clinical relevance for the treatment of resistant breast cancer.