Stratification of Tamoxifen Synergistic Combinations for the Treatment of ER+ Breast Cancer

Abstract

Simple Summary

ER+ breast cancer is the most diagnosed subtype and patient prognosis has improved in recent years thanks largely to the development of endocrine-targeting therapies including tamoxifen. Unfortunately, many tumors will recur as endocrine-therapy-resistant metastases. Therefore, in order to prolong resistance-free survival, combination therapies that mitigate resistance mechanisms should be pursued. In this study, we tested 516 drug combinations with tamoxifen and identified two that were synergistic. These combinations inhibited breast cancer patient-derived xenograft (PDX) growth in animal models better than either drug or tamoxifen alone.

Abstract

Breast cancer alone accounts for the majority of cancer deaths among women, with the most commonly diagnosed subtype being estrogen receptor positive (ER+). Survival has greatly improved for patients with ER+ breast cancer, due in part to the development of antiestrogen compounds, such as tamoxifen. While treatment of the primary disease is often successful, as many as 30% of patients will experience recurrence and metastasis, mainly due to developed endocrine therapy resistance. In this study, we discovered two tamoxifen combination therapies, with simeprevir and VX-680, that reduce the tumor burden in animal models of ER+ breast cancer more than either compound or tamoxifen alone. Additionally, these tamoxifen combinations reduced the expression of HER2, a hallmark of tamoxifen treatment, which can facilitate acquisition of a treatment-resistant phenotype. These combinations could provide clinical benefit by potentiating tamoxifen treatment in ER+ breast cancer.

1. Introduction

1.1. ER+ Breast Cancer and Standard of Care

Breast cancer is the most commonly diagnosed cancer and leading cause of cancer death among women [1]. Breast tumors can be classified by (1) the expression of three receptors, namely, the estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2) and (2) differential transcriptomic signatures such as the PAM50, which, among additional clinical biomarkers, are used to parse breast cancers into one of five intrinsic molecular subtypes (luminal A, luminal B, HER2-enriched, basal-like, and normal-like) [2,3]. ER+ breast cancers are stratified into either the luminal A or B subtype, for which proliferative status and PR expression are differentiating factors. Approximately 70–80% of breast cancers are ER+ at diagnosis and are largely driven by estrogenic activity [4].

In addition to chemotherapy, radiation, and surgical resection, current standard of care regimens for women with ER+ breast cancer focus on the reduction in circulating estrogen or antagonizing ER [5]. These adjuvant endocrine therapy options include selective estrogen-receptor receptor modulators (SERM), the most common of which is tamoxifen, as well as aromatase inhibitors (AI), including letrozole, anastrozole, and exemestane, which block estrogen production [6]. Data from randomized clinical trials suggest five years of adjuvant endocrine therapy reduce recurrence by 50%, and extending treatment to 10 years further improves disease-free survival even beyond the cessation of treatment [7]. AIs are generally prescribed for post-menopausal women, as considerable amounts of estrogen precursors in the ovary limit the therapeutic potential of these compounds in pre-menopausal patients [8]. Tamoxifen is prescribed for pre-menopausal women and any post-menopausal women who may be unable to take AIs due to contraindications or adverse effects [9]. Other treatment options include selective ER degraders (SERD) such as fulvestrant, and oophorectomy. More recently, CDK4/6 inhibitors including palbociclib, ribociclib, and abemaciclib, as well as PI3K targeting compounds, have been introduced to first- and second-line therapy [5,10].

1.2. Etiology of Endocrine Resistance

Though often initially clinically responsive, approximately 30% of patients with ER+ disease will develop long-term endocrine resistance [11]. Considering the loss of ER expression accounts for only 10% of endocrine-resistant breast cancers, understanding the different modalities of the emergence of therapeutic resistance is critical [12]. Much work has been carried out to understand the molecular underpinnings of tamoxifen resistance specifically. Tamoxifen resistance can be the result of direct modulation of ER signaling, the upregulation of receptor tyrosine kinases (HER2, EGFR, FGFR, etc.), aberrations in the PI3K/AKT/mTOR pathway, or activation of NFkB signaling [13]. Targeting these pathways in addition to ER could increase durability, or extend progression-free survival, in patients with ER+ disease, as well as decreasing endocrine therapy resistance.

To overcome resistance, or to potentiate antiestrogen therapies, combination regimens are often prescribed [14]. Many clinical studies evaluating tamoxifen in combination with other compounds and treatments have been performed. For example, trials pairing tamoxifen with EGFR inhibitors (NCT04504331) have been conducted. Additionally, evaluating gonadotropin-releasing hormone agonists (NCT00066703) or the suppression of ovarian function (NCT00066690) with tamoxifen have been performed. As tamoxifen has been used in the treatment of other hormone-sensitive cancers, additional trials have tested the efficacy of tamoxifen with SUBA-Itraconzole, an anti-fungal medication, for the treatment of ovarian cancer (NCT05156892).

1.3. PDX Models/Models of ER+ BC

Patient-derived xenograft (PDX) models are established by directly implanting human tumor cells orthotopically into immune-deficient mice [15]. PDX models are useful in studying translational science, as they faithfully recapitulate facets of human disease such as tumor heterogeneity, therapeutic response, and genomic evolution [16,17]. For these reasons, in vivo and ex vivo PDX models are a useful complement to cell lines [18,19].

1.4. Approach

To overcome drug resistance and increase treatment durability, tamoxifen combination therapies may be pursued. Because novel approaches are necessary, we utilize high-throughput screening (HTS) to discover therapeutic combinations with tamoxifen. To ensure transferability, the combinations that reduced cell viability to the greatest degree were further evaluated, and two were tested in vivo.

We show that tamoxifen and either VX-680, an aurora kinase inhibitor, or simeprevir, a hepatitis C antiviral compound, are synergistic in PDX models of ER+ breast cancer. Additionally, both combinations effectively reduce tumor volume in vivo, more than either compound or tamoxifen alone. We illustrate the capability of these tamoxifen combination therapies to suppress the upregulation of HER2, which can facilitate endocrine resistance. These results highlight the importance of addressing tamoxifen resistance mechanisms in the form of primary treatment to garner a more durable response.

2. Materials and Methods

2.1. Patient-Derived Xenograft (PDX) Models

HCI-011 and HCI-013 were obtained from the Huntsman Cancer Institute, University of Utah; BCM-5097, BCM-15057, and BCM-15034 were obtained from Baylor College of Medicine, Houston, TX, USA. All mouse studies were performed in accordance with the Virginia Commonwealth University (VCU) Institutional Animal Care and Use Committee (IACUC). Xenografts were grown in the fourth mammary fat pad set of 5–7-week-old female non-obese diabetic severe combined immunodeficient gamma (NSG) mice (NOD.Cg-Prkdcscid Il2rgtm1Wjl/SzJ, JAX #005557) [20,21]. Tumor digestion was carried out as previously described. Briefly, tumors were harvested once approximately 300 mm³, finely minced, and digested for one hour at 37 °C with agitation in DMEM/F12 supplemented with 5% FBS, 300 U/mL collagenase (Sigma, St. Louis, MO, USA), and 100 U/mL hyaluronidase (Sigma, St. Louis, MO, USA). Red blood cells were lysed and a single-cell suspension was achieved through trypsin digestion. Cells were resuspended in Hanks’ balanced salt solution (HBSS) with 2% FBS for downstream applications. For in vivo expansion and propagation of PDXs, tumor cells were suspended 1:1 in Cultrex (Bio-Teche, Minneapolis, MN, USA) and injected bilaterally into the fourth set of mammary glands of NSG mice at 500,000 cells per injection. Silastic estrogen pellets (2 mg) were implanted subcutaneously at the time of injection to supplement exogenous estrogen [22]. Tumor growth was monitored longitudinally by weekly caliper measurements.

Development of Endocrine-Resistant PDX Models

For generation of the BCM-15057EI variant, cells were grown as orthotopic tumors in intact female NSG mice with supplemental estrogen via subcutaneous silastic pellet. Once tumors reached approximately 50 mm³, the subcutaneous estrogen pellet was surgically removed. Following tumor regression and subsequent resurgence, tumors were harvested and passaged into ovary-intact mice without exogenous estrogen. Tumors were harvested when they reached 300 mm³ and passaged into ovariectomized mice without supplemental estrogen, wherein tumor growth was observed, and the subline was deemed estrogen independent (EI). To generate a tamoxifen-resistant PDX, BCM-15034 cells were grown orthotopically with supplemental estrogen. Once tumor burden reached 50 mm³, mice were transferred to tamoxifen citrate chow (Envigo, Indianapolis, IN, USA, 600 mg/kg). Tumors were then serially passaged into mice and grown in the presence of tamoxifen, until no difference in growth rate was observed comparing tumors grown under standard breeder chow versus tamoxifen-supplemented chow, wherein this isogenic line was deemed tamoxifen-resistant (BCM-15034TamR).

2.2. Cell Lines

Human ER+ breast cancer cell lines MCF7 and T47D were obtained from the American Type Culture Collection (ATCC, Manassas, VA, USA), and cultured in RPMI-1640 supplemented with GlutaMAX, 10% fetal bovine serum (FBS), and penicillin/streptomycin. UCD12 cells were a generous gift from Carol Sartorius, Ph.D. (University of Colorado Anschutz Medical Campus, Aurora, CO 80045 USA). UCD12 cells were cultured in Dulbecco’s modified minimum essential medium (DMEM)/F12 containing L-glutamine (365  mg/L) buffered with sodium bicarbonate (1200  mg/L) and HEPES (3575  mg/L) with 10% FBS, cholera toxin (100  ng/mL), hydrocortisone (1  μg/mL), insulin (10⁻⁹ M), and penicillin/streptomycin. All cell lines were maintained at 37  °C with 5% CO₂.

2.3. Bulk RNA Sequencing

2.3.1. Bulk RNA Sequencing Analysis

RNA was isolated using Qiagen RNeasy and QIAshredder kits according to the manufacturer’s protocol from flash-frozen tumor fragments. Quality of RNA was assessed using TapeStation (Agilent Technologies, Santa Clara, CA, USA). Library preparation was performed using the NEBNext Ultra II RNA Library Prep Kit for Illumina (New England Biolabs, Ipswich, MA, USA) according to the manufacturer’s instructions. Oligod(T) beads were then used to enrich for mRNA, from which cDNA libraries were generated. After end repair, universal adapters were ligated via PCR with limited cycles. Library validation and quality control were performed by TapeStation and qPCR (KAPA Biosystems, Wilmington, MA, USA).

Illumina HiSeq 4000 was used to perform 2 × 150 bp paired-end sequencing, according to instrument protocols by Azenta Life Sciences. Base calling and image analysis were conducted using HiSeq Control Software (HCS) v3.3.2-v3.4.0. Raw BCL data were converted to FASTQ file format and demultiplexed with Illumina BCL2FASTQ 2.17. One mismatch was allowed for index sequence identification.

2.3.2. Quality Control and Pre-Processing

Preprocessing of bulk RNA-seq data has been described previously [23,24]. Briefly, assessment of sequencing quality was conducted using FastQC v0.11.8 [25]. Adapters and low-quality base pairs were removed using CutAdapt v1.15 [26]. Reads of sufficient quality were aligned using STAR v2.5.2b to a merged human/mouse genome using the command line options: “--outSAMtype BAM Unsorted --outSAMorder Paired --outReadsUnmapped Fastx --quantMode TranscriptomeSAM --outFilterMultimapNmax 1. The Salmon v0.8.2; “quant” algorithm was used to obtain read counts from the aligned BAM files using the “IU” library type [27,28]. For details on merged genome construction, see Alzubi et al. [24]. Log2TPM values were calculated in R and used for gene signature computations and PAM50 subtyping was carried out using genefu v2.11.2 R package [29]. New and previously generated data were included in this study (Supplementary Table S1).

2.3.3. Gene Signatures and Clustering

Previously published gene expression signatures were used to generate signature scores for each PDX by averaging bulk RNA-seq Log2 TPM values from genes included in from each signature [30,31,32,33,34,35,36,37,38,39,40,41,42,43,44]. Unsupervised hierarchical clustering was performed in Morpheus using one minus Pearson’s correlation as the distance metric to cluster both rows and columns of the generated gene signature score matrix.

2.4. Targeted Mutational Profiling

Oncomine Comprehensive Assay v3 was performed as described [45]. Briefly, PDX tumor fragments were flash frozen and embedded in Optimal Cutting Temperature (OCT) compound prior to nucleic acid isolation. Next-generation sequencing (NGS) in accordance with Oncomine Comprehensive Assay v3 (Thermo Fisher, Waltham, MA, USA) was performed by the VCU Pathology Molecular Diagnostics laboratory with clinically validated methodology used for diagnostic tumor profiling.

2.5. Single-Cell RNA Sequencing

Single-Cell RNA Sequencing Quality Control and Preprocessing

Single-cell RNA sequencing (scRNA-seq) was performed on five ER+ breast cancer PDXs (BCM-5097, BCM-15057, BCM-15034, HCI-011, and HCI-013) that were grown orthotopically in NSG mice that were fed normal breeder chow, or in some cases tamoxifen citrate chow. Tumors were dissociated as described previously, and cells were prepared using the Chromium singles cell gene expression kit (10x Genomics, Pleasanton, CA, USA), according to the manufacturer’s instructions. Sequencing was performed at the VCU Genomics core. Reads were aligned to the 10x Genomics human/mouse genome (human being GRCh38, and mouse mm10) in order to remove reads from mouse cells. The remaining human cells were realigned to the human genome, and dead/poor-quality cell filtering was performed using an in-house R v4.1.3 script with the Seurat v.4.3.0 package as described in Boyd et al. [18]. Samples were then normalized using log normalization and merged using Seurat’s merge () function. The merged data were saved as a 10xcounts “.h5” file with all cell annotations saved as .CSV files. The CellRanger “reanalyze” function was then employed to convert the H5 file to a “.cloupe” file so that it could be uploaded Into the 10x Loupe browser for exploration.

2.6. Drug Screening

High-throughput screening was performed on PDX and cell line models of ER+ breast cancer using a 516-drug library at 1 μM concentration, as reported previously [46]. Two biological replicates were averaged, and coefficient of drug interaction (CDI) was calculated for each drug in each model, using the formula CDI = AB/(A × B), where A and B are the viability of tamoxifen and each agent alone, and AB is the viability of the combination [23].

2.6.1. In Vitro

In vitro drug screening of PDX cells has been previously described [23,46,47]. In brief, single-cell suspensions were plated in 96-well plates at a density of 25,000 cells/well in M87 medium and treated for 72 h. For cell line viability assays, 1250 cells per well were plated in 96-well plates in RPMI-1640 supplemented with GlutaMAX, 10% FBS, and penicillin/streptomycin and allowed to adhere for 24 h. Cells were treated for 72 h, and viability was measured relative to vehicle-treated control using the CellTiter-Glo Luminescent Viability Assay (Promega, Madison, WI, USA), according to the manufacturer’s protocol. For 4-hydroxytamoxifen synergism screens, after cells had adhered, media were removed and replaced with RPMI + GlutaMAX with 10% FBS, 1% penicillin/streptomycin, and 4-hydroxytamoxifen at the cell line-specific IC20 dose. Cells were then administered drugs and viability was determined as described above.

2.6.2. In Vivo

Pilot studies were performed to determine tolerable doses of simeprevir (MedChem Express, Junction, NJ, USA, HY-10241) and VX-680 (MedChem Express, HY-10161). For in vivo experiments, tumors were seeded, and, once palpable, mice were randomized into respective treatment groups. VX-680 (70 mg/kg) and simeprevir (40 mg/kg) were suspended in sterile saline and PEG300 at a 50:50 ratio. Doses were administered 5 days per week via intraperitoneal (IP) injection. Tamoxifen treatment was administered via chow (Envigo, 600 mg/kg). Tumors were measured 3 times per week, and tumor area was calculated using the formula A = L/W.

2.7. Western Blot

Tumor fragments were homogenized in a buffer containing 8 M Urea, 1% SDS, 200 mM EPPS pH 8.5 with protease and phosphatase inhibitors. Protein electrophoresis and transfer were performed using BIO-RAD PowerPac Basic and BIO-RAD Tran-Blot Turbo, respectively. Immunoblotting was performed using antibodies against ER (Abcam ab16660, Cambridge, UK), TFF1 (Cell Signaling #15571, Danvers, MA, USA) and GAPDH (Cell Signaling #5174). Detection was achieved using LI-COR Odyssey Fc (987-15227) with Image Studio v5.2.

2.8. Immunohistochemistry (IHC)

PDX tissues were processed and embedded by the Tissue and Data Acquisition and Analysis Core (TDAAC) at VCU. Antigen retrieval was performed using EDTA/TRIS Antigen retrieval buffer in a decloaking chamber (DakoCytomation, Glostrup, Denmark, Pascal). Immunohistochemical staining was performed using antibodies against ER (Abcam ab16660), PR (Cell Signaling #8757), HER2 (Cell Signaling #2242), Ki67 (Cell Signaling #9027), aurora kinase A (Cell Signaling #91590), and aurora kinase B (Cell Signaling #28711). HRP-conjugated anti-rabbit secondary binding was detected using DAB and images were acquired using Zen 2 (blue edition) software.

2.9. Quantification and Statistical Analysis

Densitometry of western blot experiments was quantified using Empiria Studio^® Software v2.2 (LI-COR, Lincoln, NE, USA). Quantification of IHC was conducted using Fiji (ImageJ 2.9.0). All statistical analysis was performed using GraphPad Prism version 9.

3. Results

3.1. Tamoxifen Treatment Response Was Classified in Five ER+ PDX Models In Vivo

Histological analysis of five ER+ PDX models was performed to confirm the expression of hormone receptors, HER2, and the proliferation marker Ki67. Additionally, in vivo dose response curves were collected to establish sensitivity to tamoxifen in ER+ breast cancer PDX models. Analysis of the five ER+ PDX models showed varying expressions of ER and PR. HCI-011 and HCI-013 highly expressed ER and PR, as well as Ki67; however, HCI-013 was extremely responsive to tamoxifen treatment, with the tumor burden being reduced by more than half. Conversely, BCM-5097 was intrinsically resistant to tamoxifen, with the tumor burden being comparable between tamoxifen-treated and untreated groups. HER2 expression was most evident in BCM-15034, which is expected, as this model was generated from a patient with HER2-amplified disease (Figure 1A, Supplementary Figures S1 and S2). Therefore, this model is considered triple positive (ER+ PR+ HER2 amplified) [48]. In all other models interrogated, tumor growth was delayed compared to control; however, tamoxifen treatment did not decrease BCM-15057, BCM-15035 or HCI-011 tumor size at any point (Figure 1B).

3.2. Tumor Heterogeneity Is Conserved in PDX Models

Genomic and transcriptomic analysis was performed to evaluate the intrinsic PAM50 subtypes and mutational landscape of a bank of ER+, HER2-amplified, and TNBC PDXs and cell line models. First, bulk RNA sequencing was performed and integrated with The Cancer Genome Atlas (TCGA) transcriptomic data to demonstrate all 20 PDX models were grouped within the expected subtype and distributed throughout the TCGA data (Figure 2A). The predicted pathogenicity of missense, frameshift, nonsense, and fusion mutations, as well as CNV, splice variants, and the introduction of stop codons in all PDX models and 14 cell lines, was established using the targeted mutational profiling tool Oncomine v3 (Figure 2B) [23]. The most common mutations were found in TP53, with the second most common being PIK3CA, which has been previously reported for many of these models [49]. The ESR1, ARID1A, and FGFR1 mutations were represented within the ER+ PDX samples. Additionally, we discovered that BCM-15057 expressed the ESR1-CCDC170 fusion mutant, which is associated with a poor response to endocrine therapy, and has not previously been reported, to the best of our knowledge [50]. Additionally, BCM-15034 exhibits a copy number variation of 5.54 of ERBB2, the gene that codes for the HER2 protein (Figure 2B, Supplementary Figure S2B). Bulk RNA sequencing data from five ER+ PDX models were used to generate gene signature scores for a selection of gene sets (Figure 2C). The hierarchical clustering of these signatures demonstrated the heterogeneity within the ER+ models. In all, the PDX models in this study were representative of the heterogeneity of human disease with regards to mutational landscape, intrinsic subtyping, and gene expression.

3.3. Genes Associated with Proliferation and ER Signaling Were Profiled at the Single-Cell Level in ER+ PDX Models

Single-cell transcriptomic profiling was performed to establish transcriptome patterns within the ER+ PDX models. After processing, 21,499 cells from the PDX mammary tumors across the five samples were merged, which displayed clusters unique to each model (Figure 3A,B). Each model expressed ESR1, PGR, and ERBB2 to varying degrees, which corresponded to established protein product expression (Figure 1A and Figure 3C). ER-target gene expression also varied among models, with FOS being the most ubiquitously expressed transcript under estrogen-stimulated conditions, and HCI-011 expressing the most robust ER-target gene signature (Figure 3D). Cell-cycle-associated genes RB1 and cyclin-dependent and aurora kinases were most expressed in BCM-15057, HCI-011, and HCI-013 (Figure 3E), corresponding to the tumor growth rate (Figure 1B), as well as to the proportion of cells in the G2M/S phase in each model (Figure 3B). In addition to an amplification of FGFR1 within BCM-15057, which was revealed by Oncomine, transcriptional upregulation of TACC1 was also present (Figure 3F). Interestingly, TACC1 can form a complex with aurora kinase A, which controls translation and cell division in breast cancer [51]. It has been established that HCI-011 tumors harbor a pathogenic mutation to PIK3CA. Here, we demonstrate that HCI-011 displayed a high gene signature score for overexpressed PIK3CA. Clinically, tumors that express pathogenic PIK3CA mutations are often treated with targeted therapy; however, BYL-719 combination therapy exhibited additive to antagonistic effects in HCI-011 in vitro (Figure 3G) [23].

3.4. Tamoxifen Treatment Upregulated Genes Associated with Hormone Therapy Resistance

We next sought to evaluate changes induced by tamoxifen treatment by comparing tamoxifen-treated HCI-013 with untreated cells, as well as BCM-15034 and an isogenic BCM-15034 that has become tamoxifen resistant (BCM-15034TamR). In both cases, tamoxifen reduced the average expression of PR, and increased ERBB2 transcript expression (Figure 4). A comparison of genes that were differentially expressed between the HCI-013 tamoxifen-treated and untreated groups showed increased expression of transcriptional regulators (CITED2, TSC22D1, and PBX1) and decreased expression of genes associated with responses to cAMP and translation (SERBP1, IGFBP5, CKB) (Figure 5A). While the majority of transcripts were expressed at a lower level in BCM-15034TamR when compared to the parental line, only two genes were expressed at a higher level: PIP and TM4SF1 (Figure 5A,B). Interestingly, only three genes were differentially regulated following treatment in both models. TM4SF1 transcripts were increased in both BCM-15034TamR and HCI-013 following treatment, and ATP5ME was expressed at a lower level in both models. However, S100P was more highly expressed in tamoxifen-treated HCI-013 when compared to the untreated PDX and was expressed less in BCM-15034TamR (Figure 5C–E). Regardless of the differences in gene expression, gene signature analysis showed lower early and late estrogen response gene signature scores in tamoxifen-treated cells, as expected (Figure 5F–G). Additionally, gene signatures that show the greatest difference in tamoxifen-treated groups included lower scores for proliferative and pro-survival signatures and higher scores in EGFR- and HER2-associated signatures (Figure 5H). These results validated the previously established upregulation of EGFR and HER2, which occurs as a response to the downregulation of estrogen signaling and can be associated with acquisition of an endocrine-therapy-resistant phenotype [52].

3.5. HTS Identified VX-680 and Simeprevir as Additively Effective Compounds with Tamoxifen

We next sought to discover novel combinatorial agents to be administered with tamoxifen as a more durable form of treatment. First, PDX suspension cultures were screened using 4-OH tamoxifen (4-OH TAM) alone to establish a dose that decreases viability to roughly 90%. While all models remained sensitive to tamoxifen metabolites in vitro, BCM-5097 remained insensitive (Figure 1B and Figure 6A). HTS with a library of 516 single agents alone or in combination with 4-OH TAM was performed on 3 ER+ cell lines and 5 PDX suspension culture models to determine which compounds reduced cell viability below that which could be accounted for by either compound as a single agent (Figure 6B,C). Forty-five combinations were found to be additive in all models based on the coefficient of drug interaction (CDI), with the exception of BCM-5097 which was excluded from analysis due to being insensitive to tamoxifen alone (Figure 6D). All known targets of possible agents were evaluated at the RNA level (Figure 6E). Aurora kinase A and B transcripts were highly expressed in all models, as were ATP1A1 and CKD4. However, some targets of the most synergistic compounds are currently unknown in the context of breast cancer and were therefore excluded from this analysis. After excluding compounds that were contraindicated with tamoxifen, the two most additive combinations, an aurora kinase inhibitor, VX-680, and a Hepatitis C antiviral compound, simeprevir, were selected for further study (Figure 6F) [53,54].

3.6. Tamoxifen Combination Therapy Was Effective in PDX Models of ER+ Breast Cancer

To evaluate tamoxifen and simeprevir or VX-680 combinations in vivo, BCM-15057 and HCI-011 were chosen based on their moderate sensitivity to tamoxifen (Figure 1B). BCM-15034 was excluded from this study as it expressed high levels of HER2 and would be treated with HER2-targeting therapeutics in a clinical setting. In both models, drug combination groups showed a profound reduction in tumor volume when compared to vehicle-treated control mice, as well as a significantly better reduction than the groups treated with either single agent or tamoxifen alone (Figure 7A–D). ER, PR, and TFF1 levels were evaluated as indications that tamoxifen treatment effectively downregulated ER signaling. As expected, ER target genes TFF1 and PR were downregulated, and ER levels were slightly increased (Figure 8A–D, Supplementary Figures S3–S5). Additionally, aurora kinase A and B protein expression was analyzed and reduced expression was seen following treatment with VX-680 (Figure 8C–F). Both combinations were effective in reducing tumor volume beyond that of tamoxifen as a mono-agent or either simeprevir or VX-680 alone, indicating that combination therapy could provide a more durable treatment for patients with ER+ breast cancer.

3.7. Tamoxifen Combination Therapy Abrogated Tamoxifen-Treatment-Induced Upregulation of HER2

Endocrine therapy resistance with respect to tamoxifen has been widely studied. It has been established that estrogen signaling downregulates HER2 expression, which can be reversed by antiestrogens such as tamoxifen [55]. This upregulation of HER2 can drive tamoxifen resistance by facilitating ER-independent proliferative and pro-survival signaling. We established that expression of HER2 is increased in all tamoxifen-sensitive models following tamoxifen treatment (Figure 9A, Supplementary Figure S6). Interestingly, when BCM-15057 and HCI-011 were treated with tamoxifen and simeprevir or VX-680 combination therapy, HER2 was not upregulated to the same degree as in tumors treated with tamoxifen alone (Figure 9B–D, Supplementary Figure S7). This suggests that antiestrogen combination therapies should be pursued that mitigate endocrine therapy resistance mechanisms as a more durable form of treatment.

4. Discussion

Although endocrine therapy has been utilized for the treatment of ER+ breast cancer for decades, as many as 40% of patients will suffer local or distant recurrence, with many cases acquiring endocrine therapy resistance [56]. More durable treatment options may include combination therapy with antiestrogens to prolong remission-free survival and decrease the incidence of endocrine therapy resistance. In this study, we identified two different tamoxifen combination regimens, simeprevir and VX-680, that eliminate the upregulation of HER2, which is often induced by tamoxifen treatment and can facilitate endocrine therapy resistance. Clinically, simeprevir or VX-680 may be added to an existing tamoxifen treatment regimen if HER2-dependent progression occurs. Additionally, this combination may be useful in instances of AI resistance, as evidenced by both HCI-011 and BCM-15057 being sensitive to tamoxifen combination therapy in vivo, even though both models were generated from samples taken from patients who had progressed on AI treatment.

The mechanism by which simeprevir and tamoxifen combination therapy leads to a reduced tumor burden, and reduced HER2 expression, is currently unknown, and may necessitate further study. Additionally, many studies have shown the benefit of aurora kinase inhibitors, as well as other cell cycle inhibitors, in combination with antiestrogens [57,58,59]. Although both tamoxifen and aurora kinase inhibitors are cytostatic, to our knowledge, this is the first study to identify the additional benefit of reduction in HER2 upregulation provided by combination therapy when compared to tamoxifen monotherapy. This may account for the improved anti-tumor activity of this particular combination. Additionally, while this study utilized tamoxifen combination therapies administered concomitantly, it remains unknown if these compounds may provide benefit if prescribed in series, or if pre-treatment with tamoxifen affects the reduction in tumor burden and HER2 expression. These gaps in knowledge may be the subject of further studies. Nevertheless, the present study demonstrates the need for antiestrogen combinations to suppress known mechanisms of endocrine resistance to potentiate their activity.

While tamoxifen has improved patient outcomes by means of overall survival and recurrence-free survival, 30% of patients ultimately progress on tamoxifen therapy or shortly after cessation [60,61,62]. Although this has been established for many years, it is recently that the US FDA has approved the first widely utilized tamoxifen combination therapy [63]. Abemaciclib with tamoxifen was approved for the treatment of ER+ breast cancer with high risk of recurrence in 2021, and as of 2023 has been expanded to include a larger patient population. Therefore, establishing any downstream effects on known mechanisms of endocrine-therapy resistance induced by abemaciclib–tamoxifen combination therapy may also be beneficial.

In this study, we sought to develop novel therapeutic combinations with tamoxifen to enhance anti-tumor efficacy with the hope of developing more durable treatments for ER+ breast cancer. We identified two compounds, simeprevir and VX-680, through HTS that were able to reduce tumor volume in vivo more than either compound or tamoxifen alone. Additionally, tumor cells from the combination groups did not show upregulation of HER2, which is a hallmark of tamoxifen treatment. Our results demonstrate the importance of developing combinations that address the potential mechanisms of endocrine therapy resistance as a more durable form of treatment for ER+ breast cancer patients.

5. Conclusions

The present study identifies two novel tamoxifen combination therapies that increase the efficacy of tamoxifen in vivo. The results reported herein establish the importance of pursuing combinations that reduce the expression or activation of known resistance mechanisms, which can provide more durable treatment in the clinical setting.

Author Contributions

Study design—E.K.Z., J.M.G. and J.C.H.; writing of the manuscript—E.K.Z., D.C.B., A.L.O., J.M.G. and J.C.H.; in vitro studies—E.K.Z. and J.M.G.; in vivo studies—E.K.Z., D.C.B., T.J.L., J.M.G., N.S.H. and C.A.S.; PDX maintenance—E.K.Z., N.S.H. and T.J.L.; Oncomine—D.C.B., A.F.-G., S.A.T.; western blots—E.K.Z. and D.C.B.; sc/bulk RNA prep—J.M.G., M.F.E. and D.C.B.; histology—E.K.Z., D.C.B., A.K.D.; scRNA sequencing, quality control—A.L.O. and M.G.D. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

This study was approved by Virginia Commonwealth University IACUC, protocol number AD10001247 (approved 19 May 2021).

Informed Consent Statement

This study does not have human subjects, and therefore there is no informed consent statement.

Data Availability Statement

The scRNA seq data presented in this manuscript are available online. The data have been deposited on GEO; accession numbers are listed in the Supplementary Materials.

Acknowledgments

The data from Figure 2 were drawn in part from data generated by the TCGA Research Network: https://www.cancer.gov/tcga (accessed on 5 January 2023). Services and products in support of the research project were generated by the Virginia Commonwealth University Tissue and Data Acquisition and Analysis Core Laboratory and the VCU Cancer Mouse Models Core, supported, in part, with funding from the Massey Cancer Center from NIH-NCI Cancer Center Support Grant P30 CA016059. Morpheus (https://software.broadinstitute.org/morpheus (accessed on 3 January 2023)) was used for the RNA sequencing heatmaps and hierarchal clustering. High-performance computing resources provided by the High Performance Research Computing (HPRC) core facility at Virginia Commonwealth University (https://hprc.vcu.edu (accessed on 3 February 2023)) were used for conducting part of the research reported in this work.

Conflicts of Interest

The authors declare no conflict of interest.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers15123179/s1, Figure S1: Representative histograms of signal intensity of IHC images corresponding to Figure 1A; Figure S2: (A) Heatmap representing patient ER, PR and HER2 status from which PDX models were generated, as reported by the Baylor College of Medicine PDX Portal. (B) Copy number of HER2 in PDX model BCM-15034 as established by Oncomine; Figure S3: Complete western blot images of in vivo experiments included in Figure 8; Figure S4: Representative histograms of signal intensity of IHC images corresponding to Figure 8C; Figure S5: Representative histograms of signal intensity of IHC images corresponding to Figure 8D; Figure S6: Representative histograms of signal intensity of IHC images corresponding to Figure 9A; Figure S7: Representative histograms of signal intensity of IHC images corresponding to Figure 9B,C; Table S1: Bulk and scRNA samples used in this study.

References

1. Łukasiewicz, S.; Czeczelewski, M.; Forma, A.; Baj, J.; Sitarz, R.; Stanisławek, A. Breast Cancer-Epidemiology, Risk Factors, Classification, Prognostic Markers, and Current Treatment Strategies—An Updated Review. Cancers; 2021; 13, 4287. [DOI: https://dx.doi.org/10.3390/cancers13174287] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34503097]

2. Perou, C.M.; Sørlie, T.B.; Eisen, M.B.; van de Rijn, M.; Jeffrey, S.S.; Rees, C.A.; Pollack, J.R.; Ross, D.T.; Johnsen, H.; Akslen, L.A. et al. Molecular portraits of human breast tumours. Nature; 2000; 406, pp. 747-752. [DOI: https://dx.doi.org/10.1038/35021093] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10963602]

3. Parker, J.S.; Mullins, M.; Cheang, M.C.; Leung, S.; Voduc, D.; Vickery, T.; Davies, S.; Fauron, C.; He, X.; Hu, Z. et al. Supervised risk predictor of breast cancer based on intrinsic subtypes. J. Clin. Oncol.; 2009; 27, pp. 1160-1167. [DOI: https://dx.doi.org/10.1200/JCO.2008.18.1370] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19204204]

4. DeSantis, C.E.; Ma, J.; Gaudet, M.M.; Newman, L.; Miller, K.; Sauer, A.G.; Jemal, A.; Siegel, R. Breast cancer statistics, 2019. CA Cancer J. Clin.; 2019; 69, pp. 438-451. [DOI: https://dx.doi.org/10.3322/caac.21583] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31577379]

5. Burguin, A.; Diorio, C.; Durocher, F. Breast Cancer Treatments: Updates and New Challenges. J. Pers. Med.; 2021; 11, 808. [DOI: https://dx.doi.org/10.3390/jpm11080808]

6. Haines, C.N.; Wardell, S.E.; McDonnell, D.P. Current and emerging estrogen receptor-targeted therapies for the treatment of breast cancer. Essays Biochem.; 2021; 65, pp. 985-1001. [DOI: https://dx.doi.org/10.1042/EBC20200174]

7. Pistilli, B.; Lohrisch, C.; Sheade, J.; Fleming, G.F. Personalizing Adjuvant Endocrine Therapy for Early-Stage Hormone Receptor–Positive Breast Cancer. Am. Soc. Clin. Oncol. Educ. Book; 2022; 42, pp. 60-72. [DOI: https://dx.doi.org/10.1200/EDBK_350358]

8. Fabian, C.J. The what, why and how of aromatase inhibitors: Hormonal agents for treatment and prevention of breast cancer. Int. J. Clin. Pract.; 2007; 61, pp. 2051-2063. [DOI: https://dx.doi.org/10.1111/j.1742-1241.2007.01587.x]

9. Blondon, M.; Bodmer, A.; Thouvenin, L.; Lecompte, T.; Righini, M.; Fontana, P.; Casini, A. Differential impact of tamoxifen and aromatase inhibitors on thrombin generation: The prospective HEMOBREAST cohort. Blood Adv.; 2022; 6, pp. 2884-2892. [DOI: https://dx.doi.org/10.1182/bloodadvances.2021006623]

10. O’Brien, N.A.; McDermott, M.S.J.; Conklin, D.; Luo, T.; Ayala, R.; Salgar, S.; Chau, K.; DiTomaso, E.; Babbar, N.; Su, F. et al. Targeting activated PI3K/mTOR signaling overcomes acquired resistance to CDK4/6-based therapies in preclinical models of hormone receptor-positive breast cancer. Breast Cancer Res.; 2020; 22, 89. [DOI: https://dx.doi.org/10.1186/s13058-020-01320-8]

11. Achinger-Kawecka, J.; Valdes-Mora, F.; Luu, P.-L.; Giles, K.A.; Caldon, C.E.; Qu, W.; Nair, S.; Soto, S.; Locke, W.J.; Yeo-The, N.S. et al. Epigenetic reprogramming at estrogen-receptor binding sites alters 3D chromatin landscape in endocrine-resistant breast cancer. Nat. Commun.; 2020; 11, 320. [DOI: https://dx.doi.org/10.1038/s41467-019-14098-x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31949157]

12. Shiino, S.; Kinoshita, T.; Yoshida, M.; Jimbo, K.; Asaga, S.; Takayama, S.; Tsuda, H. Prognostic Impact of Discordance in Hormone Receptor Status Between Primary and Recurrent Sites in Patients With Recurrent Breast Cancer. Clin. Breast Cancer; 2016; 16, pp. e133-e140. [DOI: https://dx.doi.org/10.1016/j.clbc.2016.05.014] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27268749]

13. Hosford, S.R.; Miller, T.W. Clinical potential of novel therapeutic targets in breast cancer: CDK4/6, Src, JAK/STAT, PARP, HDAC, and PI3K/AKT/mTOR pathways. Pharm. Pers. Med.; 2014; 7, pp. 203-215.

14. Olson, E. Combination Therapies in Advanced, Hormone Receptor-Positive Breast Cancer. J. Adv. Pract. Oncol.; 2018; 9, pp. 43-54.

15. Murayama, T.; Gotoh, N. Patient-Derived Xenograft Models of Breast Cancer and Their Application. Cells; 2019; 8, 621. [DOI: https://dx.doi.org/10.3390/cells8060621]

16. Woo, X.Y.; Giordano, J.; Srivastava, A.; Zhao, Z.M.; Lloyd, M.W.; de Bruijn, R.; Suh, Y.S.; Patidar, R.; Chen, L.; Scherer, S. et al. Conservation of copy number profiles during engraftment and passaging of patient-derived cancer xenografts. Nat. Genet.; 2021; 53, pp. 86-99. [DOI: https://dx.doi.org/10.1038/s41588-020-00750-6] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33414553]

17. Shi, J.; Li, Y.; Jia, R.; Fan, X. The fidelity of cancer cells in PDX models: Characteristics, mechanism and clinical significance. Int. J. Cancer; 2020; 146, pp. 2078-2088. [DOI: https://dx.doi.org/10.1002/ijc.32662] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31479514]

18. Weeber, F.; Ooft, S.N.; Dijkstra, K.K.; Voest, E.E. Tumor Organoids as a Pre-clinical Cancer Model for Drug Discovery. Cell Chem. Biol.; 2017; 24, pp. 1092-1100. [DOI: https://dx.doi.org/10.1016/j.chembiol.2017.06.012]

19. Williams, J.A. Using PDX for Preclinical Cancer Drug Discovery: The Evolving Field. J. Clin. Med.; 2018; 7, 41. [DOI: https://dx.doi.org/10.3390/jcm7030041]

20. Coughlan, A.M.; Harmon, C.; Whelan, S.; O’Brien, E.C.; O’Reilly, V.P.; Crotty, P.; Kelly, P.; Ryan, M.; Hickey, F.B.; O’Farrelly, C. et al. Myeloid Engraftment in Humanized Mice: Impact of Granulocyte-Colony Stimulating Factor Treatment and Transgenic Mouse Strain. Stem Cells Dev.; 2016; 25, pp. 530-541. [DOI: https://dx.doi.org/10.1089/scd.2015.0289]

21. Shultz, L.D.; Lyons, B.L.; Burzenski, L.M.; Gott, B.; Chen, X.; Chaleff, S.; Kotb, M.; Gillies, S.D.; King, M.; Mangada, J. et al. Human lymphoid and myeloid cell development in NOD/LtSz-scid IL2R gamma null mice engrafted with mobilized human hemopoietic stem cells. J. Immunol.; 2005; 174, pp. 6477-6489. [DOI: https://dx.doi.org/10.4049/jimmunol.174.10.6477] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15879151]

22. Kabos, P.; Finlay-Schultz, J.; Li, C.; Kline, E.; Finlayson, C.; Wisell, J.; Manuel, C.A.; Edgerton, S.M.; Harrell, J.C.; Elias, A. et al. Patient-derived luminal breast cancer xenografts retain hormone receptor heterogeneity and help define unique estrogen-dependent gene signatures. Breast Cancer Res. Treat.; 2012; 135, pp. 415-432. [DOI: https://dx.doi.org/10.1007/s10549-012-2164-8]

23. Boyd, D.C.; Zboril, E.K.; Olex, A.L.; Leftwich, T.J.; Hairr, N.S.; Byers, H.A.; Valentine, A.D.; Altman, J.E.; Alzubi, M.A.; Grible, J.M. et al. Discovering Synergistic Compounds with BYL-719 in PI3K Overactivated Basal-like PDXs. Cancers; 2023; 15, 1582. [DOI: https://dx.doi.org/10.3390/cancers15051582] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36900375]

24. Alzubi, M.A.; Turner, T.H.; Olex, A.L.; Sohal, S.S.; Tobin, N.P.; Recio, S.G.; Bergh, J.; Hatschek, T.; Parker, J.S.; Sartorius, C.A. et al. Separation of breast cancer and organ microenvironment transcriptomes in metastases. Breast Cancer Res.; 2019; 21, 36. [DOI: https://dx.doi.org/10.1186/s13058-019-1123-2]

25. Andrews, S. FastQC: A Quality Control Tool for High throughput Sequence Data. 2010; Volume 5, no. 1 Available online: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 27 January 2023).

26. Martin, M. CUTADAPT removes adapter sequences from high-throughput sequencing reads. EMBnet. J.; 2011; 17, pp. 10-12. [DOI: https://dx.doi.org/10.14806/ej.17.1.200]

27. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics; 2013; 29, pp. 15-21. [DOI: https://dx.doi.org/10.1093/bioinformatics/bts635]

28. Patro, R.; Duggal, G.; Love, M.I.; Irizarry, R.A.; Kingsford, C. Salmon provides fast and bias-aware quantification of transcript expression. Nat. Methods; 2017; 14, pp. 417-419. [DOI: https://dx.doi.org/10.1038/nmeth.4197] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28263959]

29. Genefu|Princess Margaret Bioinformatics and Computational Genomics Laboratory. Available online: https://www.pmgenomics.ca/bhklab/software/genefu (accessed on 27 January 2023).

30. Gatza, M.L.; Lucas, J.E.; Barry, W.T.; Kim, J.W.; Wang, Q.; Crawford, M.D.; Datto, M.B.; Kelley, M.; Mathey-Prevot, B.; Potti, A. et al. A pathway-based classification of human breast cancer. Proc. Natl. Acad. Sci. USA; 2010; 107, pp. 6994-6999. [DOI: https://dx.doi.org/10.1073/pnas.0912708107] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20335537]

31. Hoadley, K.A.; Weigman, V.J.; Fan, C.; Sawyer, L.R.; He, X.; Troester, M.A.; Sartor, C.I.; Rieger-House, T.; Bernard, P.S.; Carey, L.A. et al. EGFR associated expression profiles vary with breast tumor subtype. BMC Genom.; 2007; 8, 258. [DOI: https://dx.doi.org/10.1186/1471-2164-8-258]

32. Sheffer, M.; Lowry, E.; Beelen, N.; Borah, M.; Amara, S.N.; Mader, C.C.; Roth, J.A.; Tsherniak, A.; Freeman, S.S.; Dashevsky, O. et al. Genome-scale screens identify factors regulating tumor cell responses to natural killer cells. Nat. Genet.; 2021; 53, pp. 1196-1206. [DOI: https://dx.doi.org/10.1038/s41588-021-00889-w]

33. Parkes, E.E.; Humphries, M.P.; Gilmore, E.; Sidi, F.A.; Bingham, V.; Phyu, S.M.; Craig, S.; Graham, C.; Miller, J.; Griffin, D. et al. The clinical and molecular significance associated with STING signaling in breast cancer. NPJ Breast Cancer; 2021; 7, 81. [DOI: https://dx.doi.org/10.1038/s41523-021-00283-z] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34172750]

34. Giraddi, R.R.; Chung, C.-Y.; Heinz, R.E.; Balcioglu, O.; Novotny, M.; Trejo, C.L.; Dravis, C.; Hagos, B.M.; Mehrabad, E.M.; Rodewald, L.W. et al. Single-Cell Transcriptomes Distinguish Stem Cell State Changes and Lineage Specification Programs in Early Mammary Gland Development. Cell Rep.; 2018; 24, pp. 1653-1666.e7. [DOI: https://dx.doi.org/10.1016/j.celrep.2018.07.025] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30089273]

35. Witkiewicz, A.K.; Balaji, U.; Knudsen, E.S. Systematically defining single-gene determinants of response to neoadjuvant chemotherapy reveals specific biomarkers. Clin. Cancer Res.; 2014; 20, pp. 4837-4848. [DOI: https://dx.doi.org/10.1158/1078-0432.CCR-14-0885] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25047707]

36. Heng, Y.J.; Lester, S.C.; Tse, G.M.; Factor, R.E.; Allison, K.H.; Collins, L.C.; Chen, Y.Y.; Jensen, K.C.; Johnson, N.B.; Jeong, J.C. et al. The molecular basis of breast cancer pathological phenotypes. J. Pathol.; 2017; 241, pp. 375-391. [DOI: https://dx.doi.org/10.1002/path.4847]

37. Pfefferle, A.D.; Herschkowitz, J.I.; Usary, J.; Harrell, J.C.; Spike, B.T.; Adams, J.R.; Torres-Arzayus, M.I.; Brown, M.; Egan, S.E.; Wahl, G.M. et al. Transcriptomic classification of genetically engineered mouse models of breast cancer identifies human subtype counterparts. Genome Biol.; 2013; 14, R125. [DOI: https://dx.doi.org/10.1186/gb-2013-14-11-r125]

38. Herschkowitz, J.I.; Simin, K.; Weigman, V.J.; Mikaelian, I.; Usary, J.; Hu, Z.; Rasmussen, K.E.; Jones, L.P.; Assefnia, S.; Chandrasekharan, S. et al. Identification of conserved gene expression features between murine mammary carcinoma models and human breast tumors. Genome Biol.; 2007; 8, R76. [DOI: https://dx.doi.org/10.1186/gb-2007-8-5-r76]

39. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S. et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA; 2005; 102, pp. 15545-15550. [DOI: https://dx.doi.org/10.1073/pnas.0506580102]

40. Oh, D.S.; Troester, M.A.; Usary, J.; Hu, Z.; He, X.; Fan, C.; Wu, J.; Carey, L.A.; Perou, C.M. Estrogen-regulated genes predict survival in hormone receptor-positive breast cancers. J. Clin. Oncol.; 2006; 24, pp. 1656-1664. [DOI: https://dx.doi.org/10.1200/JCO.2005.03.2755]

41. Herschkowitz, J.I.; He, X.; Fan, C.; Perou, C.M. The functional loss of the retinoblastoma tumour suppressor is a common event in basal-like and luminal B breast carcinomas. Breast Cancer Res.; 2008; 10, R75. [DOI: https://dx.doi.org/10.1186/bcr2142]

42. Fan, C.; Prat, A.; Parker, J.S.; Liu, Y.; Carey, L.A.; Troester, M.A.; Perou, C.M. Building prognostic models for breast cancer patients using clinical variables and hundreds of gene expression signatures. BMC Med. Genom.; 2011; 4, 3. [DOI: https://dx.doi.org/10.1186/1755-8794-4-3]

43. Lim, E.; Vaillant, F.; Wu, D.; Forrest, N.C.; Pal, B.; Hart, A.H.; Asselin-Labat, M.L.; Gyorki, D.E.; Ward, T.; Partanen, A. et al. Aberrant luminal progenitors as the candidate target population for basal tumor development in BRCA1 mutation carriers. Nat. Med.; 2009; 15, pp. 907-913. [DOI: https://dx.doi.org/10.1038/nm.2000] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19648928]

44. Hennessy, B.T.; Gonzalez-Angulo, A.M.; Stemke-Hale, K.; Gilcrease, M.Z.; Krishnamurthy, S.; Lee, J.S.; Fridlyand, J.; Sahin, A.; Agarwal, R.; Joy, C. et al. Characterization of a naturally occurring breast cancer subset enriched in epithelial-to-mesenchymal transition and stem cell characteristics. Cancer Res.; 2009; 69, pp. 4116-4124. [DOI: https://dx.doi.org/10.1158/0008-5472.CAN-08-3441] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19435916]

45. Finlay-Schultz, J.; Jacobsen, B.M.; Riley, D.; Paul, K.V.; Turner, S.; Ferreira-Gonzalez, A.; Harrell, J.C.; Kabos, P.; Sartorius, C.A. New generation breast cancer cell lines developed from patient-derived xenografts. Breast Cancer Res.; 2020; 22, 68. [DOI: https://dx.doi.org/10.1186/s13058-020-01300-y]

46. Turner, T.H.; Alzubi, M.A.; Harrell, J.C. Identification of synergistic drug combinations using breast cancer patient-derived xenografts. Sci. Rep.; 2020; 10, 1493. [DOI: https://dx.doi.org/10.1038/s41598-020-58438-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32001757]

47. Rashid, N.S.; Hairr, N.S.; Murray, G.; Olex, A.L.; Leftwich, T.J.; Grible, J.M.; Reed, J.; Dozmorov, M.G.; Harrell, J.C. Identification of nuclear export inhibitor-based combination therapies in preclinical models of triple-negative breast cancer. Transl. Oncol.; 2021; 14, 101235. [DOI: https://dx.doi.org/10.1016/j.tranon.2021.101235]

48. Zeng, J.; Edelweiss, M.; Ross, D.S.; Xu, B.; Moo, T.A.; Brogi, E.; D’Alfonso, T.M. Triple-Positive Breast Carcinoma: Histopathologic Features and Response to Neoadjuvant Chemotherapy. Arch. Pathol. Lab Med.; 2021; 145, pp. 728-735. [DOI: https://dx.doi.org/10.5858/arpa.2020-0293-OA] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33112958]

49. Guillen, K.P.; Fujita, M.; Butterfield, A.J.; Scherer, S.D.; Bailey, M.H.; Chu, Z.; DeRose, Y.S.; Zhao, L.; Cortes-Sanchez, E.; Yang, C.H. et al. A human breast cancer-derived xenograft and organoid platform for drug discovery and precision oncology. Nat. Cancer; 2022; 3, pp. 232-250. [DOI: https://dx.doi.org/10.1038/s43018-022-00337-6]

50. Li, L.; Lin, L.; Veeraraghavan, J.; Hu, Y.; Wang, X.; Lee, S.; Tan, Y.; Schiff, R.; Wang, X.-S. Therapeutic role of recurrent ESR1-CCDC170 gene fusions in breast cancer endocrine resistance. Breast Cancer Res.; 2020; 22, 84. [DOI: https://dx.doi.org/10.1186/s13058-020-01325-3]

51. Conte, N.; Delaval, B.; Ginestier, C.; Ferrand, A.; Isnardon, D.; Larroque, C.; Prigent, C.; Séraphin, B.; Jacquemier, J.; Birnbaum, D. TACC1–chTOG–Aurora A protein complex in breast cancer. Oncogene; 2003; 22, pp. 8102-8116. [DOI: https://dx.doi.org/10.1038/sj.onc.1206972]

52. Viedma-Rodríguez, R.; Baiza-Gutman, L.; Salamanca-Gómez, F.; Diaz-Zaragoza, M.; Martínez-Hernández, G.; Esparza-Garrido, R.R.; Velázquez-Flores, M.A.; Arenas-Aranda, D. Mechanisms associated with resistance to tamoxifen in estrogen receptor-positive breast cancer (review). Oncol. Rep.; 2014; 32, pp. 3-15. [DOI: https://dx.doi.org/10.3892/or.2014.3190]

53. Medchemexpress.com. Tozasertib. Available online: https://www.medchemexpress.com/Tozasertib.html (accessed on 10 March 2023).

54. Simeprevir. 2023; Available online: https://www.medchemexpress.com/Simeprevir.html (accessed on 10 March 2023).

55. Mazumder, A.; Shiao, S.; Haricharan, S. HER2 Activation and Endocrine Treatment Resistance in HER2-negative Breast Cancer. Endocrinology; 2021; 162, bqab153. [DOI: https://dx.doi.org/10.1210/endocr/bqab153] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34320193]

56. Pedersen, R.N.; Esen, B.; Mellemkjær, L.; Christiansen, P.; Ejlertsen, B.; Lash, T.L.; Nørgaard, M.; Cronin-Fenton, D. The Incidence of Breast Cancer Recurrence 10–32 Years After Primary Diagnosis. J. Natl. Cancer Inst.; 2022; 114, pp. 391-399. [DOI: https://dx.doi.org/10.1093/jnci/djab202] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34747484]

57. Yde, C.W. Aurora A and Mcl-1: New potential treatment targets in antiestrogen-resistant breast cancer. Mol. Cell. Oncol.; 2021; 8, 998898. [DOI: https://dx.doi.org/10.1080/23723556.2014.998898]

58. Larsen, S.L.; Yde, C.W.; Laenkholm, A.V.; Rasmussen, B.B.; Duun-Henriksen, A.K.; Bak, M.; Lykkesfeldt, A.E.; Kirkegaard, T. Aurora kinase B is important for antiestrogen resistant cell growth and a potential biomarker for tamoxifen resistant breast cancer. BMC Cancer; 2015; 15, 239. [DOI: https://dx.doi.org/10.1186/s12885-015-1210-4] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25885472]

59. Tolaney, S.M.; Beeram, M.; Beck, J.T.; Conlin, A.; Dees, E.C.; Puhalla, S.L.; Rexer, B.N.; Burris, H.A.; Jhaveri, K.; Helsten, T. et al. Abemaciclib in Combination With Endocrine Therapy for Patients With Hormone Receptor-Positive, HER2-Negative Metastatic Breast Cancer: A Phase 1b Study. Front. Oncol.; 2021; 11, 810023. [DOI: https://dx.doi.org/10.3389/fonc.2021.810023]

60. Hwang, K.T.; Kim, E.K.; Jung, S.H.; Lee, E.S.; Kim, S.I.; Lee, S.; Park, H.K.; Kim, J.; Oh, S.; Kim, Y.A. Tamoxifen therapy improves overall survival in luminal A subtype of ductal carcinoma in situ: A study based on nationwide Korean Breast Cancer Registry database. Breast Cancer Res. Treat.; 2018; 169, pp. 311-322. [DOI: https://dx.doi.org/10.1007/s10549-018-4681-6]

61. Yu, N.Y.; Iftimi, A.; Yau, C.; Tobin, N.P.; van ’t Veer, L.; Hoadley, K.A.; Benz, C.C.; Nordenskjöld, B.; Fornander, T.; Stål, O. et al. Assessment of Long-term Distant Recurrence-Free Survival Associated with Tamoxifen Therapy in Postmenopausal Patients With Luminal A or Luminal B Breast Cancer. JAMA Oncol.; 2019; 5, pp. 1304-1309. [DOI: https://dx.doi.org/10.1001/jamaoncol.2019.1856]

62. Pan, H.; Gray, R.; Braybrooke, J.; Davies, C.; Taylor, C.; McGale, P.; Peto, R.; Pritchard, K.I.; Bergh, J.; Dowsett, M. et al. 20-Year Risks of Breast-Cancer Recurrence after Stopping Endocrine Therapy at 5 Years. N. Engl. J. Med.; 2017; 377, pp. 1836-1846. [DOI: https://dx.doi.org/10.1056/NEJMoa1701830]

63. Hamilton, E.; Cortes, J.; Ozyilkan, O.; Chen, S.-C.; Petrakova, K.; Manikhas, A.; Jerusalem, G.; Hegg, R.; Huober, J.; Chapman, S.C. et al. nextMONARCH: Abemaciclib Monotherapy or Combined with Tamoxifen for Metastatic Breast Cancer. Clin. Breast Cancer; 2021; 21, pp. 181-190.e2. [DOI: https://dx.doi.org/10.1016/j.clbc.2020.09.011]