BackgroundCHEK2 c.1100delC was the first moderate-risk breast cancer (BC) susceptibility allele discovered. Despite several genomic, transcriptomic and functional studies, however, it is still... Show moreBackgroundCHEK2 c.1100delC was the first moderate-risk breast cancer (BC) susceptibility allele discovered. Despite several genomic, transcriptomic and functional studies, however, it is still unclear how exactly CHEK2 c.1100delC promotes tumorigenesis. Since the mutational landscape of a tumor reflects the processes that have operated on its development, the aim of this study was to uncover the somatic genomic landscape of CHEK2-associated BC.MethodsWe sequenced primary BC (pBC) and normal genomes of 20 CHEK2 c.1100delC mutation carriers as well as their pBC transcriptomes. Including pre-existing cohorts, we exhaustively compared CHEK2 pBC genomes to those from BRCA1/2 mutation carriers, those that displayed homologous recombination deficiency (HRD) and ER- and ER+ pBCs, totaling to 574 pBC genomes. Findings were validated in 517 metastatic BC genomes subdivided into the same subgroups. Transcriptome data from 168 ER+ pBCs were used to derive a TP53-mutant gene expression signature and perform cluster analysis with CHEK2 BC transcriptomes. Finally, clinical outcome of CHEK2 c.1100delC carriers was compared with BC patients displaying somatic TP53 mutations in two well-described retrospective cohorts totaling to 942 independent pBC cases.ResultsBC genomes from CHEK2 mutation carriers were most similar to ER+ BC genomes and least similar to those of BRCA1/2 mutation carriers in terms of tumor mutational burden as well as mutational signatures. Moreover, CHEK2 BC genomes did not show any evidence of HRD. Somatic TP53 mutation frequency and the size distribution of structural variants (SVs), however, were different compared to ER+ BC. Interestingly, BC genomes with bi-allelic CHEK2 inactivation lacked somatic TP53 mutations and transcriptomic analysis indicated a shared biology with TP53 mutant BC. Moreover, CHEK2 BC genomes had an increased frequency of > 1 Mb deletions, inversions and tandem duplications with peaks at specific sizes. The high chromothripsis frequency among CHEK2 BC genomes appeared, however, not associated with this unique SV size distribution profile.ConclusionsCHEK2 BC genomes are most similar to ER+ BC genomes, but display unique features that may further unravel CHEK2-driven tumorigenesis. Increased insight into this mechanism could explain the shorter survival of CHEK2 mutation carriers that is likely driven by intrinsic tumor aggressiveness rather than endocrine resistance. Show less
Noord, V.E. van der; Stel, W. van der; Louwerens, G.; Verhoeven, D.; Kuiken, H.J.; Lieftink, C.; ... ; Water, B. van de 2023
Triple-negative breast cancer (TNBC) is a subtype of breast cancer with limited treatment options and poor clinical prognosis. Inhibitors of transcriptional CDKs are currently under thorough... Show moreTriple-negative breast cancer (TNBC) is a subtype of breast cancer with limited treatment options and poor clinical prognosis. Inhibitors of transcriptional CDKs are currently under thorough investigation for application in the treatment of multiple cancer types, including breast cancer. These studies have raised interest in combining these inhibitors, including CDK12/13 inhibitor THZ531, with a variety of other anti-cancer agents. However, the full scope of these potential synergistic interactions of transcriptional CDK inhibitors with kinase inhibitors has not been systematically investigated. Moreover, the mechanisms behind these previously described synergistic interactions remain largely elusive.\nKinase inhibitor combination screenings were performed to identify kinase inhibitors that synergize with CDK7 inhibitor THZ1 and CDK12/13 inhibitor THZ531 in TNBC cell lines. CRISPR-Cas9 knockout screening and transcriptomic evaluation of resistant versus sensitive cell lines were performed to identify genes critical for THZ531 resistance. RNA sequencing analysis after treatment with individual and combined synergistic treatments was performed to gain further insights into the mechanism of this synergy. Kinase inhibitor screening in combination with visualization of ABCG2-substrate pheophorbide A was used to identify kinase inhibitors that inhibit ABCG2. Multiple transcriptional CDK inhibitors were evaluated to extend the significance of the found mechanism to other transcriptional CDK inhibitors.\nWe show that a very high number of tyrosine kinase inhibitors synergize with the CDK12/13 inhibitor THZ531. Yet, we identified the multidrug transporter ABCG2 as key determinant of THZ531 resistance in TNBC cells. Mechanistically, we demonstrate that most synergistic kinase inhibitors block ABCG2 function, thereby sensitizing cells to transcriptional CDK inhibitors, including THZ531. Accordingly, these kinase inhibitors potentiate the effects of THZ531, disrupting gene expression and increasing intronic polyadenylation.\nOverall, this study demonstrates the critical role of ABCG2 in limiting the efficacy of transcriptional CDK inhibitors and identifies multiple kinase inhibitors that disrupt ABCG2 transporter function and thereby synergize with these CDK inhibitors. These findings therefore further facilitate the development of new (combination) therapies targeting transcriptional CDKs and highlight the importance of evaluating the role of ABC transporters in synergistic drug-drug interactions in general. Show less
BackgroundHormone receptor (HR)-positive, HER2/neu-negative breast cancers have a sustained risk of recurrence up to 20 years from diagnosis. TEAM (Tamoxifen, Exemestane Adjuvant Multinational) is... Show moreBackgroundHormone receptor (HR)-positive, HER2/neu-negative breast cancers have a sustained risk of recurrence up to 20 years from diagnosis. TEAM (Tamoxifen, Exemestane Adjuvant Multinational) is a large, multi-country, phase III trial that randomized 9776 women for the use of hormonal therapy. Of these 2754 were Dutch patients. The current study aims for the first time to correlate the ten-year clinical outcomes with predictions by CanAssist Breast (CAB)—a prognostic test developed in South East Asia, on a Dutch sub-cohort that participated in the TEAM. The total Dutch TEAM cohort and the current Dutch sub-cohort were almost similar with respect to patient age and tumor anatomical features.MethodsOf the 2754 patients from the Netherlands, which are part of the original TEAM trial, 592 patients’ samples were available with Leiden University Medical Center (LUMC). The risk stratification of CAB was correlated with outcomes of patients using logistic regression approaches entailing Kaplan–Meier survival curves, univariate and multivariate cox-regression hazards model. We used hazard ratios (HRs), the cumulative incidence of distant metastasis/death due to breast cancer (DM), and distant recurrence-free interval (DRFi) for assessment.ResultsOut of 433 patients finally included, the majority, 68.4% had lymph node-positive disease, while only a minority received chemotherapy (20.8%) in addition to endocrine therapy. CAB stratified 67.5% of the total cohort as low-risk [DM = 11.5% (95% CI, 7.6–15.2)] and 32.5% as high-risk [DM = 30.2% (95% CI, 21.9–37.6)] with an HR of 2.90 (95% CI, 1.75–4.80; P < 0.001) at ten years. CAB risk score was an independent prognostic factor in the consideration of clinical parameters in multivariate analysis. At ten years, CAB high-risk had the worst DRFi of 69.8%, CAB low-risk in the exemestane monotherapy arm had the best DRFi of 92.7% [vs CAB high-risk, HR, 0.21 (95% CI, 0.11–0.43), P < 0.001], and CAB low-risk in the sequential arm had a DRFi of 84.2% [vs CAB high-risk, HR, 0.48 (95% CI, 0.28–0.82), P = 0.009].ConclusionsCost-effective CAB is a statistically robust prognostic and predictive tool for ten-year DM for postmenopausal women with HR+/HER2−, early breast cancer. CAB low-risk patients who received exemestane monotherapy had an excellent ten-year DRFi. Show less
Background: Prediction of contralateral breast cancer (CBC) risk is challenging due to moderate performances of the known risk factors. We aimed to improve our previous risk prediction model ... Show moreBackground: Prediction of contralateral breast cancer (CBC) risk is challenging due to moderate performances of the known risk factors. We aimed to improve our previous risk prediction model (PredictCBC) by updated follow-up and including additional risk factors.Methods: We included data from 207,510 invasive breast cancer patients participating in 23 studies. In total, 8225 CBC events occurred over a median follow-up of 10.2 years. In addition to the previously included risk factors, PredictCBC-2.0 included CHEK2 c.1100delC, a 313 variant polygenic risk score (PRS-313), body mass index (BMI), and parity. Fine and Gray regression was used to fit the model. Calibration and a time-dependent area under the curve (AUC) at 5 and 10 years were assessed to determine the performance of the models. Decision curve analysis was performed to evaluate the net benefit of PredictCBC-2.0 and previous PredictCBC models.Results: The discrimination of PredictCBC-2.0 at 10 years was higher than PredictCBC with an AUC of 0.65 (95% prediction intervals (PI) 0.56-0.74) versus 0.63 (95%Pl 0.54-0.71). PredictCBC-2.0 was well calibrated with an observed/ expected ratio at 10 years of 0.92 (95%Pl 0.34-2.54). Decision curve analysis for contralateral preventive mastectomy (CPM) showed the potential clinical utility of PredictCBC-2.0 between thresholds of 4 and 12% 10-year CBC risk for BRCA1/2 mutation carriers and non-carriers.Conclusions: Additional genetic information beyond BRCA1/2 germline mutations improved CBC risk prediction and might help tailor clinical decision-making toward CPM or alternative preventive strategies. Identifying patients who benefit from CPM, especially in the general breast cancer population, remains challenging. Show less
Background: The majority of BRCA1-mutant breast cancers are characterized by a triple-negative phenotype and a basal-like molecular subtype, associated with aggressive clinical behavior. Current... Show moreBackground: The majority of BRCA1-mutant breast cancers are characterized by a triple-negative phenotype and a basal-like molecular subtype, associated with aggressive clinical behavior. Current treatment options are limited, highlighting the need for the development of novel targeted therapies for this tumor subtype. Methods: Our group previously showed that EZH2 is functionally relevant in BRCA1-deficient breast tumors and blocking EZH2 enzymatic activity could be a potent treatment strategy. To validate the role of EZH2 as a therapeutic target and to identify new synergistic drug combinations, we performed a high-throughput drug combination screen in various cell lines derived from BRCA1-deficient and -proficient mouse mammary tumors. Results: We identified the combined inhibition of EZH2 and the proximal DNA damage response kinase ATM as a novel synthetic lethality-based therapy for the treatment of BRCA1-deficient breast tumors. We show that the combined treatment with the EZH2 inhibitor GSK126 and the ATM inhibitor AZD1390 led to reduced colony formation, increased genotoxic stress, and apoptosis-mediated cell death in BRCA1-deficient mammary tumor cells in vitro. These findings were corroborated by in vivo experiments showing that simultaneous inhibition of EZH2 and ATM significantly increased anti-tumor activity in mice bearing BRCA1-deficient mammary tumors. Conclusion: Taken together, we identified a synthetic lethal interaction between EZH2 and ATM and propose this synergistic interaction as a novel molecular combination for the treatment of BRCA1-mutant breast cancer. Show less
Background Genome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor ... Show moreBackground 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. Show less
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... Show moreBackground 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. Show less
Background Breast cancer survivors (BCS) may have increased risk of hypothyroidism, but risk according to treatment modality is unclear. We estimated the incidence of hypothyroidism in women with... Show moreBackground Breast cancer survivors (BCS) may have increased risk of hypothyroidism, but risk according to treatment modality is unclear. We estimated the incidence of hypothyroidism in women with breast cancer, and according to cancer treatment. Methods Using nationwide registries, we identified all Danish women aged >= 35 years diagnosed with non-metastatic breast cancer (1996-2009). We matched up to five cancer-free women (controls) for each BCS. We excluded women with prevalent thyroid disease. Cancer treatment was chemotherapy with or without radiotherapy (RT) targeting the breast/chest wall only, or also the lymph nodes (RTn). We identified hypothyroidism using diagnostic codes, and/or levothyroxine prescriptions. We calculated the cumulative incidence, incidence rates (IR) per 1000 person-years, and used Cox regression to estimate hazard ratios (HR) and associated 95% confidence intervals (CIs) of hypothyroidism, adjusting for comorbidities. Results We included 44,574 BCS and 203,306 matched controls with 2,631,488 person-years of follow-up. BCS had a slightly higher incidence of hypothyroidism than controls [5-year cumulative incidence, 1.8% (95%CI = 1.7-1.9) and 1.6% (95%CI = 1.5-1.6), respectively]. The overall IR was 4.45 (95%CI = 4.25-4.67) and 3.81 (95%CI = 3.73-3.90), corresponding to an adjusted HR = 1.17 (95%CI = 1.11-1.24). BCS who received RTn with chemotherapy (HR = 1.74, 95%CI = 1.50-2.02) or without chemotherapy (HR = 1.31, 95%CI = 1.14-1.51) had an elevated risk of hypothyroidism compared with matched controls and compared with BCS who underwent surgery alone [HR = 1.71, 95%CI = 1.45-2.01 and HR = 1.36, 95%CI = 1.17-1.58, respectively]. Conclusions BCS have an excess risk of hypothyroidism compared with age-matched controls. BCS and those working in cancer survivorship settings ought to be aware that this risk is highest in women treated with radiation therapy to the lymph nodes and chemotherapy. Show less
Background The effect of risk-reducing salpingo-oophorectomy (RRSO) on breast cancer risk for BRCA1 and BRCA2 mutation carriers is uncertain. Retrospective analyses have suggested a protective... Show moreBackground The effect of risk-reducing salpingo-oophorectomy (RRSO) on breast cancer risk for BRCA1 and BRCA2 mutation carriers is uncertain. Retrospective analyses have suggested a protective effect but may be substantially biased. Prospective studies have had limited power, particularly for BRCA2 mutation carriers. Further, previous studies have not considered the effect of RRSO in the context of natural menopause. Methods A multi-centre prospective cohort of 2272 BRCA1 and 1605 BRCA2 mutation carriers was followed for a mean of 5.4 and 4.9 years, respectively; 426 women developed incident breast cancer. RRSO was modelled as a time-dependent covariate in Cox regression, and its effect assessed in premenopausal and postmenopausal women. Results There was no association between RRSO and breast cancer for BRCA1 (HR = 1.23; 95% CI 0.94-1.61) or BRCA2 (HR = 0.88; 95% CI 0.62-1.24) mutation carriers. For BRCA2 mutation carriers, HRs were 0.68 (95% CI 0.40-1.15) and 1.07 (95% CI 0.69-1.64) for RRSO carried out before or after age 45 years, respectively. The HR for BRCA2 mutation carriers decreased with increasing time since RRSO (HR = 0.51; 95% CI 0.26-0.99 for 5 years or longer after RRSO). Estimates for premenopausal women were similar. Conclusion We found no evidence that RRSO reduces breast cancer risk for BRCA1 mutation carriers. A potentially beneficial effect for BRCA2 mutation carriers was observed, particularly after 5 years following RRSO. These results may inform counselling and management of carriers with respect to RRSO. Show less
Background: Breast cancer survivors are at risk for contralateral breast cancer (CBC), with the consequent burden of further treatment and potentially less favorable prognosis. We aimed to develop... Show moreBackground: Breast cancer survivors are at risk for contralateral breast cancer (CBC), with the consequent burden of further treatment and potentially less favorable prognosis. We aimed to develop and validate a CBC risk prediction model and evaluate its applicability for clinical decision-making.Methods: We included data of 132,756 invasive non-metastatic breast cancer patients from 20 studies with 4682 CBC events and a median follow-up of 8.8 years. We developed a multivariable Fine and Gray prediction model (PredictCBC-1A) including patient, primary tumor, and treatment characteristics and BRCA1/2 germline mutation status, accounting for the competing risks of death and distant metastasis. We also developed a model without BRCA1/2 mutation status (PredictCBC-1B) since this information was available for only 6% of patients and is routinely unavailable in the general breast cancer population. Prediction performance was evaluated using calibration and discrimination, calculated by a time-dependent area under the curve (AUC) at 5 and 10 years after diagnosis of primary breast cancer, and an internal-external cross-validation procedure. Decision curve analysis was performed to evaluate the net benefit of the model to quantify clinical utility.Results: In the multivariable model, BRCA1/2 germline mutation status, family history, and systemic adjuvant treatment showed the strongest associations with CBC risk. The AUC of PredictCBC-1A was 0.63 (95% prediction interval (PI) at 5 years, 0.52-0.74; at 10 years, 0.53-0.72). Calibration-in-the-large was -0.13 (95% PI: -1.62-1.37), and the calibration slope was 0.90 (95% PI: 0.73-1.08). The AUC of Predict-1B at 10 years was 0.59 (95% PI: 0.52-0.66); calibration was slightly lower. Decision curve analysis for preventive contralateral mastectomy showed potential clinical utility of PredictCBC-1A between thresholds of 4-10% 10-year CBC risk for BRCA1/2 mutation carriers and non-carriers.Conclusions: We developed a reasonably calibrated model to predict the risk of CBC in women of European-descent; however, prediction accuracy was moderate. Our model shows potential for improved risk counseling, but decision-making regarding contralateral preventive mastectomy, especially in the general breast cancer population where limited information of the mutation status in BRCA1/2 is available, remains challenging. Show less
Ren, J.; Smid, M.; Iaria, J.; Salvatori, D.C.F.; Dam, H. van; Zhu, H.J.; ... ; Dijke, P. ten 2019
Introduction The presence of tumor-infiltrating lymphocytes (TILs) is correlated with good prognosis and outcome after (immuno)therapy in triple-negative and HER2-positive breast cancer. However,... Show moreIntroduction The presence of tumor-infiltrating lymphocytes (TILs) is correlated with good prognosis and outcome after (immuno)therapy in triple-negative and HER2-positive breast cancer. However, the role of TILs in luminal breast cancer is less clear. Emerging evidence has now demonstrated that genetic aberrations in malignant cells influence the immune landscape of tumors. Phosphatidylinositol 3-kinase (PI3K) is the most common altered pathway in ER-positive breast cancer. It is unknown whether changes in the PI3K pathway result in a different composition of the breast tumor microenvironment. Here we present the retrospective analysis of a prospective randomized trial in ER-positive breast cancer on the prognostic and predictive value of specific tumor-associated lymphocytes in the context of PI3K alterations. Methods We included 563 ER-positive tumors from a multicenter trial for stage I to III postmenopausal breast cancer patients, who were randomized to tamoxifen or no adjuvant therapy. The amount of CD8-, CD4-, and FOXP3-positive cells was evaluated by immunohistochemistry and quantified by imaging-analysis software. We analyzed the associations between PIK3CA hotspot mutations, PTEN expression, phosphorylated proteins of the PI3K and MAPK pathway (p-AKT, p-ERK1/2, p-4EBP1, p-p70S6K), and recurrence-free interval after adjuvant tamoxifen or no adjuvant treatment. Results CD8-positive lymphocytes were significantly more abundant in PIK3CA-mutated tumors (OR = 1.65; 95% CI 1.03-2.68). While CD4 and FOXP3 were not significantly associated with prognosis, patients with tumors classified as CD8-high had increased risk of recurrence (HR = 1.98; 95% CI 1.14-3.41; multivariable model including PIK3CA status, treatment arm, and other standard clinicopathological variables). Lymphocytes were more often present in tumors with increased PI3K downstream phosphorylation. This was most pronounced for FOXP3-positive cells. Conclusion These exploratory analyses of a prospective trial in luminal breast cancer suggest high CD8 infiltration is associated with unfavorable outcome and that PI3K pathway alterations might be associated with the composition of the tumor microenvironment. Show less
McLaughlin, R.P.; He, J.; Noord, V.E. van der; Redel, J.; Foekens, J.A.; Martens, J.W.M.; ... ; Water, B. van de 2019
The effective treatment of triple-negative breast cancer (TNBC) remains a profound clinical challenge. Despite frequent epidermal growth factor receptor (EGFR) overexpression and reliance on... Show moreThe effective treatment of triple-negative breast cancer (TNBC) remains a profound clinical challenge. Despite frequent epidermal growth factor receptor (EGFR) overexpression and reliance on downstream signalling pathways in TNBC, resistance to EGFR-tyrosine kinase inhibitors (TKIs) remains endemic. Therefore, the identification of targeted agents, which synergise with current therapeutic options, is paramount.Compound-based, high-throughput, proliferation screening was used to profile the response of TNBC cell lines to EGFR-TKIs, western blotting and siRNA transfection being used to examine the effect of inhibitors on EGFR-mediated signal transduction and cellular dependence on such pathways, respectively. A kinase inhibitor combination screen was used to identify compounds that synergised with EGFR-TKIs in TNBC, utilising sulphorhodamine B (SRB) assay as read-out for proliferation. The impact of drug combinations on cell cycle arrest, apoptosis and signal transduction was assessed using flow cytometry, automated live-cell imaging and western blotting, respectively. RNA sequencing was employed to unravel transcriptomic changes elicited by this synergistic combination and to permit identification of the signalling networks most sensitive to co-inhibition.We demonstrate that a dual cdc7/CDK9 inhibitor, PHA-767491, synergises with multiple EGFR-TKIs (lapatinib, erlotinib and gefitinib) to overcome resistance to EGFR-targeted therapy in various TNBC cell lines. Combined inhibition of EGFR and cdc7/CDK9 resulted in reduced cell proliferation, accompanied by induction of apoptosis, G2-M cell cycle arrest, inhibition of DNA replication and abrogation of CDK9-mediated transcriptional elongation, in contrast to mono-inhibition. Moreover, high expression of cdc7 and RNA polymerase II Subunit A (POLR2A), the direct target of CDK9, is significantly correlated with poor metastasis-free survival in a cohort of breast cancer patients. RNA sequencing revealed marked downregulation of pathways governing proliferation, transcription and cell survival in TNBC cells treated with the combination of an EGFR-TKI and a dual cdc7/CDK9 inhibitor. A number of genes enriched in these downregulated pathways are associated with poor metastasis-free survival in TNBC.Our results highlight that dual inhibition of cdc7 and CDK9 by PHA-767491 is a potential strategy for targeting TNBC resistant to EGFR-TKIs. Show less