Abstract
Acquired resistance to HER2-targeted therapies occurs frequently in HER2+ breast tumors and new strategies for overcoming resistance are needed. Here, we report that resistance to trastuzumab is reversible, as resistant cells regained sensitivity to the drug after being cultured in drug-free media. RNA-sequencing analysis showed that cells resistant to trastuzumab or trastuzumab + pertuzumab in combination increased expression of oxidative phosphorylation pathway genes. Despite minimal changes in mitochondrial respiration, these cells exhibited increased expression of ATP synthase genes and selective dependency on ATP synthase function. Resistant cells were sensitive to inhibition of ATP synthase by oligomycin A, and knockdown of ATP5J or ATP5B, components of ATP synthase complex, rendered resistant cells responsive to a low dose of trastuzumab. Furthermore, combining ATP synthase inhibitor oligomycin A with trastuzumab led to regression of trastuzumab-resistant tumors in vivo. In conclusion, we identify a novel vulnerability of cells with acquired resistance to HER2-targeted antibody therapies and reveal a new therapeutic strategy to overcome resistance.
These findings implicate ATP synthase as a novel potential target for tumors resistant to HER2-targeted therapies.
Introduction
About 20%–30% of breast cancers exhibit gene amplification and/or overexpression of HER2 (gene name ERBB2) and are classified as HER2 positive (HER2+; ref. 1). HER2 is a receptor tyrosine kinase that forms dimers with itself or other family members, including EGFR and HER3, to activate downstream signaling. Signaling from HER family kinases controls several important processes that are often deregulated in cancers such as cell growth, proliferation, motility, and survival (reviewed in ref. 2). The standard of care for HER2+ breast cancer includes a HER2-targeted agent in addition to chemotherapy and surgery. The first and most widely used HER2-targeted therapy is a humanized mAb called trastuzumab (trade name Herceptin; Genentech/Roche). Trastuzumab functions by HER2 downregulation (3) and degradation (4), blocking ligand-independent HER2 dimerization and downstream signaling (5–9), and through antibody-dependent cellular cytotoxicity (10, 11). The addition of trastuzumab to a chemotherapy regimen was shown to significantly increase median progression-free survival (PFS) and overall survival (OS) in metastatic breast cancers by 3 and 5 months, respectively (12). Rates of disease-free survival (DFS) and OS were also significantly improved in patients with early-stage breast cancer treated with trastuzumab (13–15). Pertuzumab (trade name Perjeta; Genentech/Roche) is a humanized mAb that can be used in combination with trastuzumab and chemotherapy for the treatment of HER2+ breast cancers at early (16) and late (17) stages. It targets the dimerization domain of HER2 (18) to block ligand-induced dimerization between HER2 and other family members, reducing activation of important downstream signaling pathways such as the PI3K pathway (19). Addition of pertuzumab to a docetaxel and trastuzumab regimen in patients with metastatic breast cancer increased median PFS by 6 months and median OS by 16 months (16). Furthermore, addition of pertuzumab to trastuzumab-based adjuvant chemotherapy increased rates of invasive DFS in patients with operable, early-stage HER2+ breast cancer (17).
Unfortunately, patients with HER2+ breast cancer do not respond uniformly to HER2-targeted antibodies. Fewer than half of patients with HER2+ metastatic breast cancer initially respond to trastuzumab (20). In addition, while trastuzumab significantly extends OS in a fraction of patients with HER2+ metastatic breast cancer, the majority of those patients will experience tumor progression within only one year (12). Many patients with early-stage HER2+ treated with trastuzumab also experience recurrence (13, 15, 21). Several strategies to overcome primary and acquired trastuzumab resistance have been proposed, but at this time, with the exception of pertuzumab, few have translated into clinical improvements (22). Patients with HER2+ breast cancer treated with the combination of chemotherapy, trastuzumab, and pertuzumab may also experience tumor recurrence after adjuvant therapy for early-stage breast cancer and with metastatic breast cancer nearly all inevitably ultimately progress (16, 17). To our knowledge, there are currently no published studies investigating the mechanisms of acquired resistance to this combination therapy. These will become increasingly important, as this regimen is becoming the new standard of care for patients with HER2+ breast cancer diagnosed at any stage.
Several recent studies suggest that drug resistance to anticancer therapies develops initially through a reversible phase. For example, it has been shown in several cancer cell line models that drug-tolerant persister cells can become resensitized to the drug after being cultured without it for some time (23–26). In line with this model of plastic resistance, we show here that trastuzumab resistance was reversed when resistant cells were cultured for several passages in drug-free media. Thus, we set out to uncover novel features of trastuzumab and trastuzumab + pertuzumab resistance by narrowing our focus on changes in gene expression programs. We identified a selective dependency of HER2 antibody–resistant breast cancer cells on mitochondrial oxidative phosphorylation (OXPHOS). This is a major process by which adenosine triphosphate (ATP) is produced in cells. The flow of electrons through the mitochondrial electron transport chain creates an electrochemical proton gradient across the inner mitochondrial membrane. The energy from this gradient is harnessed by the passage of protons through ATP synthase (also called complex V) to phosphorylate adenosine diphosphate (ADP) to produce ATP in the mitochondrial matrix.
Metabolic rewiring to shift reliance toward OXPHOS or glycolysis has been shown to contribute to tumor recurrence in many settings (reviewed in ref. 27). Trastuzumab was shown to hinder glycolysis in breast cancer cells, and the combination of trastuzumab with an inhibitor of glycolysis synergistically decreased tumor cell growth (28). Trastuzumab resistance was also linked to an increase in glycolysis (28, 29). Thus, the glycolytic contribution to trastuzumab resistance has been well characterized in HER2+ breast cancer, but to our knowledge, the contribution from OXPHOS has never been considered. Here, we demonstrate that resistant cells upregulated OXPHOS-related genes and were more sensitive to OXPHOS blockade through inhibition of ATP synthase, even though sensitive and resistant cells respired at a similar rate. Combining trastuzumab with an ATP synthase inhibitor led to the regression of trastuzumab-resistant tumors in vivo. We also show that expression of the ATP synthase gene ATP5B correlated with poor outcome in patients with HER2+ breast cancer. These novel findings may provide a new avenue of treatment for patients with trastuzumab- or trastuzumab + pertuzumab–resistant disease.
Materials and Methods
Chemicals and antibodies
DMSO was purchased from Sigma-Aldrich. Oligomycin A (11342) and rotenone (13995) were purchased from Cayman Chemical. 2-Deoxy-d-glucose (2-DG; 202010) was purchased from Santa Cruz Biotechnology, Inc. Antimycin A (A8674) was purchased from Sigma-Aldrich. Trastuzumab and pertuzumab were provided by or purchased from Genentech/Roche. Anti-phospho-EGFR (D7A5), anti-EGFR (4267), anti-phospho-HER2 (2241), anti-HER2 (4290), anti-phospho-HER3 (2842), anti-HER3 (12708), anti-phospho-AKT S473 (4060), anti-AKT (4691), anti-phospho-ERK (4376), anti-ERK (46955), anti-β-actin (3700), and anti-vinculin (13901) antibodies were purchased from Cell Signaling Technology. Anti-ATP5B (ab14730) and anti-ATP5J (ab224139) antibodies were purchased from Abcam. Anti-cytokeratin (Z0622) and Envision + System HRP-labeled polymer anti-mouse (K4001) antibodies were purchased from Agilent Dako. Goat anti-rabbit Alexa-546 (A11010) and DAPI (D1306) were purchased from Life Technologies. Cyanine-5 tyramide (SAT705A001EA) was purchased from Perkin-Elmer.
Cell culture
BT474 and SKBR3 cells were obtained from the ATCC. BT474 cells were cultured in RPMI 1640. SKBR3 and HEK293T cells were cultured in DMEM. All media were supplemented with 10% FBS, 100 U/mL penicillin, and 100 μg/mL streptomycin. Resistant cell pools were gradually made resistant to trastuzumab or trastuzumab + pertuzumab by being cultured in 10 μg/mL up to 100 μg/mL over the course of approximately three months. Cells were subsequently continuously maintained in 10 μg/mL trastuzumab or 10 μg/mL trastuzumab + 10 μg/mL pertuzumab. Trastuzumab and pertuzumab were used at 10 μg/mL unless explicitly stated otherwise. To allow for direct comparison of parental cells and drug-resistant cells without the influence of drug(s), resistant cells were seeded into media without drug(s) for 2–7 days unless indicated otherwise. MCF10A cells were cultured as described previously (30). Cells were periodically tested for Mycoplasma contamination and authenticated using short tandem repeat profiling.
RNA-sequencing
RNA was isolated from BT474 cells treated for three days with 10 μg/mL trastuzumab, 10 μg/mL trastuzumab + 10 μg/mL pertuzumab, or no treatment. Resistant cells were cultured for 7 days without drugs and then RNA was harvested. RNA was harvested from three replicate 10-cm plates for each cell line using the RNeasy Plus Mini Kit (Qiagen 74136). Samples were evaluated using the Agilent Bioanalyzer before mRNA library preparation. Libraries were prepared according to standard Illumina protocols by the Yale Stem Cell Center Genomics and Bioinformatics Core and run on Illumina HiSeq2000. Bioinformatic analyses were performed as described previously (31). In brief, reads were mapped to the human genome (hg38) with Bowtie2 (32), differential gene expression was determined with DESeq2 (33), and Gene Set Enrichment Analysis (GSEA) software (34) was utilized for pathway analysis. RNA-sequencing (RNA-seq) data were deposited in the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus database under accession number GSE121105.
Oxygen consumption measurements
Agilent Seahorse XF96 plate wells (101085-004) were coated with Cell-Tak (Corning C354240) according to the manufacturer's protocol. BT474 cells ± 10 μg/mL trastuzumab for 3 days, BT-TR cells ± 10 μg/mL trastuzumab for 7 days, and BT-TPR cells ± 10 μg/mL trastuzumab and 10 μg/mL pertuzumab for 7 days were seeded at 25,000 cells/well. Each cell line was seeded into 4–8 wells. The next day, oxygen consumption measurements were obtained on the Agilent Seahorse Bioscience XF96 analyzer with injection ports loaded to attain final concentrations of 1.5 μmol/L oligomycin A, 1 μmol/L carbonyl cyanide-4-(trifluoromethoxy)phenylhydrazone (FCCP), and 1 μmol/L rotenone + 1 μmol/L antimycin A. Afterwards, cells were lysed in high salt lysis buffer and protein content was measured by Bradford assay as the loading control. This assay was performed twice with four or more replicates per condition.
Mitochondrial membrane potential measurements
BT474, BT-TR1, and BT-TR2 cells were seeded into 12-well plates at 250,000 cells/well after being kept out of drug(s) for 7 days. The next day, cells were treated with DMSO or 20 μmol/L FCCP for 10 minutes at 37°C. Subsequently, cells were exposed to 50 nmol/L solution of tetramethylrhodamine ethyl ester (TMRE; Abcam 113852) for 30 minutes at 37°C. Cells were washed with PBS, trypsinized, resuspended, and pushed through a mesh lid to separate cell clumps. Phycoerythrin (PE) signal was measured by flow cytometry. Membrane potential was calculated as the ratio of PE signal of DMSO-treated cells/FCCP-treated cells, as FCCP-treated cells represent staining associated with mitochondrial mass.
Knockdown of ATP5J and ATP5B
pLKO.1 plasmids containing short hairpins against ATP5J (TRCN0000038219-shATP5J-1, TRCN0000038220-shATP5J-2) and ATP5B (TRCN0000300454-shATP5B-1, TRCN0000300455- shATP5B-2) were obtained from Sigma-Aldrich. Scrambled shRNA plasmid was used as control (35). Insert sequences are listed in Supplementary Table S1. HEK-293T cells were transfected with the pLKO.1 plasmids and second-generation psPAX2 and pMD2.G packaging system using Lipofectamine 2000 (Thermo Fisher Scientific 12566014). Virus was sterile filtered and applied to target cells for 16 hours with 8 μg/mL polybrene (hexadimethrine bromide, Sigma H9268). Target cells were selected with 1 μg/mL puromycin and utilized without delay in downstream applications. The knockdown procedure was repeated three times independently to capture differences in cell proliferation and trastuzumab sensitivity before any cellular adaptations.
Animal studies
17β-estradiol pellets (17 mg) were purchased from Huntsman Cancer Institute (Salt Lake City, UT) and implanted subcutaneously into the loose skin at the back of the neck of female NOD scid gamma mice (4–5 weeks old). Five days later, each mouse received two cell injections into the fourth mammary fat pads on opposite flanks. The injections were comprised of 5 × 106 viable cells resuspended in sterile PBS mixed at 1:1 with Matrigel (Corning 356231) in a total volume of 0.1 mL. BT474 cells were injected into the right flank and BT-TR2 cells were injected into the left flank of each animal. After 20 days, mice were randomly assigned to four treatment groups: saline control, trastuzumab [5 mg/kg intraperitoneally (i.p.) every 3 days], oligomycin A [0.5 mg/kg i.p. daily for 3 days followed by 0.25 mg/kg i.p. daily], or the combination of trastuzumab + oligomycin A (each at the same dosage as the monotherapy). Mice in the saline and single treatment arms received control saline injections matching their counterparts receiving drug treatments (1–2 injections per day per mouse). Tumors were measured daily using digital calipers and volume calculated as 0.5 × length × width2. After 13 days of treatment, mice were euthanized and tumors were collected for histologic analyses. This work was performed according to NIH guidelines and approved by the Yale University Institutional Animal Care and Use Committee.
Tissue microarray staining and analysis
The optimal titer of ATP5B antibody was determined by staining the breast cancer standardization tissue microarray (TMA) YTMA263 at five dilutions of primary antibody as described previously (36). ATP5B antibody was used at 0.5 μg/mL for subsequent staining. A tissue microarray was constructed from BT474 cell pellets expressing shControl, shATP5B-1, and shATP5B-2 by the Yale Cancer Center/Pathology Tissue Microarray facility to evaluate ATP5B antibody specificity. Tissue microarrays containing tumor tissue from the Hellenic Oncology Group (HeCOG) 10/05 clinical trial (37) were stained, Automated Quantitative Analysis (AQUA) scores were calculated, and Kaplan–Meier estimates of DFS and OS were calculated as described previously (36). All tissue samples were used in accordance with U.S. Common Rule after approval by the Yale Human Investigation Committee. Informed written consent was obtained from all the patients.
Statistical analysis
Unless indicated otherwise, unpaired, two-tailed Student t test was used to determine significance using GraphPad Prism 7.0 software.
Results
Acquired resistance to trastuzumab is reversible
To study resistance to HER2-targeted antibodies, we used trastuzumab-sensitive BT474 and SKBR3 HER2+ breast cancer cells to generate several derivative pools of BT474 and SKBR3 cells resistant to trastuzumab or trastuzumab + pertuzumab (Fig. 1A; Supplementary Fig. S1A and S1B). We treated parental BT474 and SKBR3 cells with increasing doses of trastuzumab or both antibodies, escalating the dose when the cells appeared to be proliferating normally despite the presence of drug(s) as determined by colony formation assays (Supplementary Fig. S1C). The pools were named with the first two letters of their parental cell line (BT or SK), followed by the designation of trastuzumab-resistant (TR) or trastuzumab + pertuzumab-resistant (TPR) and a number to differentiate between pools generated separately. We first characterized the expression and phosphorylation level of HER family signaling pathway kinases in BT474 and BT-TR2 cells. We found that BT-TR1 and BT-TR2 cells have increased EGFR, HER3, and phosphorylated HER3 compared with parental cells, while they did not show a marked increase in phosphorylation of HER2, EGFR, and ERK (Supplementary Fig. S1D). Short-term trastuzumab treatment decreased phosphorylation of HER3 and ERK while inducing phosphorylation of HER2 and EGFR for both parental and drug-resistant cells (Supplementary Fig. S1D), indicating that these resistant cells respond to trastuzumab similarly to the parental cells for these phosphorylation-mediated signaling events. In contrast to gatekeeper mutations accounting for the majority of acquired resistance cases to targeted small-molecule ATP analogue tyrosine kinase inhibitors, no specific mutation or genetic event has been shown to account for the majority of acquired resistance cases to trastuzumab (38). In line with this, the trastuzumab- and trastuzumab + pertuzumab–resistant pools we generated showed no significant increases in copy number of key receptor tyrosine kinase genes EGFR, ERBB2 (HER2), ERBB3 (HER3), MET, or IGF1R compared with parental cells (Supplementary Fig. S2).
We passaged BT474-derived resistant pools side by side in trastuzumab or drug-free media (referred to as “washout”) and examined their sensitivity to trastuzumab periodically. After 20 doublings (nine passages) in drug-free media, all pools became more sensitive to the drug, and three out of four pools of resistant cells tested regained sensitivity to trastuzumab (Fig. 1B–D; Supplementary Fig. S3A and S3B). To decipher whether the pools regained sensitivity due to clonal selection or flexibility of individual clones, we generated single-cell clones from the resistant pool BT-TR2 and repeated the assay (Fig. 1E). After 23 doublings, two of three resistant clones tested regained sensitivity to trastuzumab (Fig. 1F and G). The third clone demonstrated increased sensitivity after 34 doublings (Supplementary Fig. S3C). Taken together, these results suggested that nongenetic changes may mediate resistance to trastuzumab.
The OXPHOS gene signature is enriched in resistant cells
We hypothesized that alterations in gene expression programs could be the major contributors to resistance. Thus, RNA-seq was performed for sensitive BT474 cells, two pools of BT-TR cells, and two pools of BT-TPR cells cultured in the absence of drug(s) for 7 days to exclude gene expression changes induced by the drug(s) (Supplementary Tables S2–S5). We utilized GSEA to identify differences between resistant pools and BT474 parental cells (Supplementary Tables S6–S13). Several hallmark pathways were positively enriched with nominal P < 0.05 and FDR q-value < 0.1 in each resistant pool compared with BT474 cells. Only one hallmark pathway, protein secretion, was common to both BT-TR pools, but not BT-TPR pools (Fig. 2A). Surprisingly, no pathways were common to both BT-TPR pools without also being enriched in BT-TR pools, highlighting similarities in pools resistant to single and combination therapies. Three GSEA hallmark pathways were positively enriched in all four resistant pools compared with BT474 cells: oxidative phosphorylation, fatty acid metabolism, and MYC targets V1 (Fig. 2A). OXPHOS was the top positively enriched pathway in BT-TR2, BT-TPR1, and BT-TPR2 cells, and third for BT-TR1 (Fig. 2B and C; Supplementary Tables S6–S9).
Among the pathways downregulated in resistant pools compared with BT474 cells with nominal P < 0.05 and FDR q-value < 0.1, four pathways were negatively enriched in all resistant cell pools compared with parental cells: estrogen response early, epithelial-to-mesenchymal transition, estrogen response late, and Hedgehog signaling (Fig. 2D–F). The hallmark gene set for estrogen response early was the top downregulated pathway for both BT-TR pools, while epithelial-to-mesenchymal transition was the top downregulated pathway for both BT-TPR pools (Fig. 2E and F; Supplementary Tables S10–S13). No downregulated gene set was common to BT-TR pools without also being downregulated in a BT-TPR pool (Fig. 2D; Supplementary Tables S10–S13). Interestingly, five hallmark pathways were negatively enriched in both BT-TPR pools, without significant enrichment in BT-TR pools (Fig. 2D). These BT-TPR–specific downregulated pathways include IL6–JAK–STAT3 signaling (Fig. 2G), inflammatory response, IFNγ response, cholesterol homeostasis, and coagulation (Supplementary Tables S12 and S13).
We decided to investigate the role of OXPHOS in resistance to HER2-targeted antibodies as the OXPHOS pathway gene set was highly enriched in all profiled resistant pools compared with parental cells and because OXPHOS is targetable with small-molecule inhibitors, facilitating potential translation to the clinic.
Resistant cells exhibit increased sensitivity to inhibition of OXPHOS, but not increased OXPHOS
To determine whether the increased expression of OXPHOS genes in resistant cells indicated an increased reliance on this pathway, we treated sensitive and resistant cells with oligomycin A, an OXPHOS inhibitor targeting the F0 subunit of ATP synthase. Although BT474 and BT-TR cells were all sensitive to 24-hour exposure to oligomycin A at low nanomolar concentrations, BT-TR1 and BT-TR2 cells were more sensitive than the parental BT474 cells in colony formation and proliferation assays (Fig. 3A and B; Supplementary Fig. S4A). Similarly, BT-TPR1 and BT-TPR2 cells showed more growth inhibition by oligomycin A treatment than BT474 cells in proliferation assays (Fig. 3C and D; Supplementary Fig. S4B). One out of two tested SK-TR cell lines also responded more strongly to oligomycin A than parental SKBR3 cells (Supplementary Fig. S4C). In contrast, MCF10A normal-like mammary epithelial cells did not show any growth inhibition at doses several fold greater than those required to inhibit the breast cancer cell lines (Supplementary Fig. S4D). We then asked whether trastuzumab-resistant cells were also more sensitive to other OXPHOS inhibitors. BT-TR1, BT-TR2, BT-TPR1, and BT-TPR2 cells all showed greater sensitivity to antimycin A, an inhibitor of mitochondrial complex III, and rotenone, an inhibitor of mitochondrial complex I than the parental cells (Supplementary Fig. S4E–S4H).
Given the resistant cells' enrichment in OXPHOS and several other metabolic gene programs compared with BT474 cells (Supplementary Tables S6–S13), it was expected that they were operating at a higher metabolic rate. However, using Seahorse metabolic flux assays, we found that the overall metabolic features of sensitive and resistant cells were very similar, consistent with a previous report (29). While three-day treatment with trastuzumab lowered cellular oxygen consumption rate (OCR) compared with untreated BT474 cells, BT-TR and BT-TPR cells demonstrated basal OCRs similar to or slightly less than parental cells, whether they were kept in drug or cultured in drug-free media for one week beforehand (Fig. 3E–G). Mitochondrial DNA copy number remained very similar between parental BT474 cells and BT-TR cells with or without drug treatments (Fig. 3H). While BT-TPR cells have either decreased or increased mitochondrial DNA copy number compared with parental BT474 cells (Supplementary Fig. S5A), there was no difference in mitochondrial DNA copy number between parental SKBR3 cells and SK-TR cells (Supplementary Fig. S5B). Likewise, ATP production was comparable between BT474 and BT-TR or BT-TPR cells (Fig. 3I; Supplementary Fig. S5C) and between SKBR3 and SK-TR cells (Supplementary Fig. S5D). Furthermore, we did not observe significant difference of mitochondrial membrane potential between parental BT474 and BT-TR cells as determined by TMRE signal normalized to the basal signal in the presence of the uncoupler FCCP (Fig. 3J).
Trastuzumab was reported to shut down glycolysis in HER2+ breast cancer cells (28). Consistent with that, the gene signature for glycolysis is negatively enriched in BT474 cells after 3-day trastuzumab treatment or 3-day trastuzumab + pertuzumab treatment (Supplementary Fig. S6A and S6B). Resistance to trastuzumab has been associated with increased glycolytic capacity (28, 29). In line with this, one of the two BT-TR cell lines profiled with RNA-seq, BT-TR2, showed marginal enrichment (NES = 1.16, nominal P = 0.102, FDR q = 0.237) of the glycolysis gene program (Supplementary Fig. S6C; Supplementary Table S7). This resistant cell line also displayed enhanced ability to increase glycolysis upon inhibition of mitochondrial respiration, as measured by the extracellular acidification rate (ECAR) during the Seahorse metabolic flux assay (Supplementary Fig. S6D and S6E). However, similar to BT-TR1 cells, BT-TR2 cells were less sensitive to 2-DG, a glucose analogue that inhibits glycolysis, compared with parental BT474 cells, indicating that they do not harbor a selective dependence on glycolysis (Supplementary Fig. S6F).
ATP synthase expression is increased in resistant cells, correlates with poor survival of patients with HER2+ breast cancer, and is necessary for maintenance of resistance
To better understand the resistant cells' sensitivity to OXPHOS pathway inhibitors, we looked closer into the genes that comprised the OXPHOS signature and found that several ATP synthase transcripts were increased approximately 1.5 fold in BT-TR and BT-TPR–resistant pools compared with sensitive cells (Fig. 4A). This was validated by RT-qPCR (Fig. 4B). Similar increases were observed in two single-cell clones derived from the resistant pool BT-TR2 cells (Supplementary Fig. S7A). SK-TR2 also demonstrated a moderate increase in ATP5J and ATP5B gene expression compared with parental SKBR3 cells, correlating with its increased sensitivity to oligomycin A (Supplementary Fig. S7B).
To assess the clinical relevance of increased ATP synthase expression in resistant cells, mRNA expression data of tumors collected from patients with HER2+ breast cancer (n = 67) with no prior treatment in The Cancer Genome Atlas dataset were analyzed using the online tool TIMER (39). The majority of these patients were treated with HER2-targeted therapies. We found that patients with high mRNA levels of ATP5B, one of the core components of ATP synthase, had significantly poorer survival compared with the low expression group, despite the small cohort size and limited number of events (Fig. 4C). mRNA levels of ATP5A1, another core component of the complex, followed a similar pattern, although this was not statistically significant by log-rank test due to the limited number of events (Fig. 4C). In contrast, ATP5B and ATPA1 were not significantly associated with survival in basal or luminal subtypes of breast cancer, or in the combined cohorts with all subtypes of breast cancers (Supplementary Fig. S8). The importance of ATP5B protein was explored in patients with HER2+ breast cancer from a clinical trial of adjuvant chemotherapy followed by trastuzumab (HeCOG 10/05; ref. 37). In this unique clinical trial, HER2+ patients were given trastuzumab after (rather than concurrent with) their chemotherapy. Although we observed a trend of association between expression of ATP5B protein with poor overall survival, the associations between expression of ATP5B protein and treatment outcome in this cohort were not statistically significant (Supplementary Fig. S9). However, the analysis was underpowered with only 11 and 7 events for DFS and OS, respectively, for a total of 117 patients included in the analysis (Supplementary Fig. S9).
To examine the relationship between ATP synthase function and trastuzumab sensitivity, we knocked down expression of complex component ATP5J with two independent shRNAs in BT474 and BT-TR cells (Fig. 4D). Notably, knockdown of ATP5J allowed a low dose of trastuzumab to inhibit growth of BT-TR cells, while control BT-TR cells did not respond to trastuzumab at this dose (Fig. 4E). Similarly, ATP5B downregulation with two independent shRNAs also led to reduced growth of BT-TR2 cells by a low dose of trastuzumab (Supplementary Fig. S10). These results confirm a role for ATP synthase in trastuzumab tolerance.
To interrogate these concepts in vivo, we injected BT474 and BT-TR2 cells into the mammary fat pads of female nude mice implanted with pellets to release β-estradiol. After tumor formation, mice were split into four treatment groups: control, trastuzumab, oligomycin A, and trastuzumab + oligomycin A. Despite similar in vitro growth rates, BT474 tumors grew faster than BT-TR2 tumors in vivo with control tumors more than tripling in size for BT474 over approximately 2 weeks, but only doubling for BT-TR2 (Fig. 5A–D). This may be because BT-TR2 tumors have a diminished response to estrogen compared with BT474 tumors, as we noted that the hallmark pathways for estrogen response were significantly downregulated in BT-TR2 cells compared with parental cells (Fig. 2D and E). We analyzed the effect of treatment on each tumor by comparing tumor volume relative to the pretreatment tumor volume. BT474 tumors responded to trastuzumab, but not to oligomycin A. Trastuzumab significantly stagnated tumor growth, while oligomycin A had negligible effects on growth. The combination of the two did not offer further benefit beyond the effects of trastuzumab alone (Fig. 5A). In contrast, BT-TR2 tumors exhibited only a small benefit from trastuzumab or oligomycin A as single agents. The loading dose of oligomycin A appeared to decrease tumor volume for the first days of the experiment, but tumors recovered during the maintenance dose phase. Notably, the combination of the two incited significant tumor regression (Fig. 5B). BT-TR2 tumors treated with trastuzumab + oligomycin A shrank by an average of 60% compared with their starting volumes, demonstrating the potential of this combination for treating trastuzumab-resistant tumors. This combination treatment led to significant apoptosis of tumor cells as measured by cleaved caspase-3 staining (Fig. 5E and F), without affecting overall tumor morphology or tumor cell proliferation as measured by Ki67 staining (Supplementary Fig. S11).
Discussion
In this study, we showed that trastuzumab resistance is reversible, suggesting that nonmutational changes contribute considerably to acquired trastuzumab resistance. Studies in other types of cancer have demonstrated the flexibility of drug tolerance within persister cell populations (23, 25). Here, we demonstrated that fully developed trastuzumab-resistant cell pools and clones are able to revert back to a sensitive state if cultured in drug-free media. This provides rationale for considering a clinical model that includes a drug-free hiatus for patients, followed by retreatment with trastuzumab. It also provides support for identifying new therapeutic targets based on transcriptional reprogramming in resistant cells.
Through transcriptional profiling, we identified several gene expression programs significantly up or downregulated in resistant pools compared with parental cells. Our analysis revealed similarities and differences between the single- or dual therapy–acquired resistance settings. For instance, the protein secretion pathway was exclusively upregulated in trastuzumab-resistant pools, while five pathways were exclusively downregulated in trastuzumab + pertuzumab-resistant pools. Several of the pathways downregulated only in dual-resistant pools were related to inflammation, revealing an interesting finding for future exploration. Four pathways were significantly downregulated by all resistant pools, while three were significantly upregulated. Of note, we found elevated expression of genes involved in OXPHOS, particularly ATP synthase genes, in all four tested trastuzumab- and trastuzumab + pertuzumab–resistant pools. These increases could be due to increased Myc activity (Fig. 2A), which was shown to regulate expression of complex V genes (40), or other posttranscriptional regulatory mechanisms of ATP synthase components that can lead to change of mRNA levels (41). Furthermore, we found that higher expression of ATP synthase component ATP5B significantly correlated with poor overall survival of patients with HER2+ breast cancer. These findings indicate that expression of ATP synthase complex genes, especially ATP5B, may be a good biomarker for predicting response to trastuzumab and clinical outcomes of patients with HER2+ breast cancer.
Treatment with a small-molecule tyrosine kinase inhibitor was previously shown to sensitize leukemia cancer cells to ATP synthase inhibition (42), but it was not known whether antibody-based targeted therapies had a similar consequence or if this occurred in alternate cancer settings. We demonstrated here that long-term HER2-targeted antibody treatment to generate resistant breast cancer cells also created increased sensitivity to ATP synthase inhibition. Resistant cells that demonstrated increased gene expression of ATP synthase demonstrated increased sensitivity to pharmacologic inhibition of ATP synthase compared with parental cells. This was the case for both single- and dual therapy–resistant cells, broadening the potential clinical relevance of these findings to a wider population of patients. Knockdown of ATP5J or ATP5B was sufficient to make trastuzumab-resistant cells respond to a low dose of trastuzumab not effective on the corresponding control cells. Trastuzumab had a cytostatic effect on BT474 tumors in vivo and very little effect on BT-TR2 tumors. However, the combination of oligomycin A and trastuzumab induced regression of trastuzumab-resistant tumors in vivo.
These data suggest that ATP synthase function is required for the maintenance of resistance and propose interference with ATP synthase activity as a method for overcoming resistance. Consistent with our findings, oligomycin A treatment prevented outgrowth of a therapy-resistant melanoma cell subpopulation (43) and prevented tumor recurrence in a mouse model of KRAS-mutant pancreatic cancer (44). However, while the oligomycin A–sensitive cell populations demonstrated increased OCR or ATP production in other models, the resistant cells in this study had similar OCR and ATP production compared with parental cells. These results suggest that ATP synthase in these cells is less efficient in producing ATP and that these cells functionally compensate for this defect by inducing expression of ATP synthase, therefore creating vulnerability to ATP synthase inhibition. The resistant cells showed less sensitivity to 2-DG (Supplementary Fig. S6F) and greater sensitivity to OXPHOS inhibitors than the parental cells (Fig. 3; Supplementary Fig. S4). Interestingly, Myc-positive cells also showed similar trends in sensitivity to these drugs compared with Myc-deficient cells (45). The phenotypic similarity described above is consistent with the idea that increased Myc activity leads to increased expression of OXPHOS genes in drug-resistant cells. In fact, it was reported that breast cancer cells resistant to lapatinib, an inhibitor of HER2 and EGFR, are sensitive to Myc inhibition (46).
Overall, the data presented here suggest that cells resistant to HER2-targeted antibodies have shifted their metabolism in such a way as to make them vulnerable to disruption of ATP synthase function. This approach as a cancer therapy is nuanced because all cells require mitochondrial respiration, but our data suggest that a therapeutic window may exist. The low nanomolar doses of oligomycin A used in vitro selectively slowed the growth of breast cancer cells, but not normal-like breast epithelial cells. Daily oligomycin A injection in vivo did not significantly affect blood counts or markers of renal and hepatic toxicities as assessed by Alvarez-Calderon and colleagues (42) and it was tolerated in experiments performed by us and others (44). Alternative inhibitors of ATP synthase may demonstrate stronger safety profiles and should be explored. For example, FDA-approved drugs such as bedaquiline (47), simvastatin (48), paroxetine (48), and tamoxifen (48) were recently identified as inhibitors of ATP synthase activity and may be better tolerated for this purpose. Taken together, our study revealed that ATP synthase inhibition is a new avenue to defeat resistance to HER2-targeted antibodies in breast cancers.
Disclosure of Potential Conflicts of Interest
M.P. DiGiovanna is a consultant at Merck and has received speakers bureau honoraria from Total Health Information Services LLC. No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
Conception and design: M. Gale, D.F. Stern, Q. Yan
Development of methodology: M. Gale
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): M. Gale, Y. Li, J. Cao, M.A. Holmbeck, M. Zhang, S.M. Lang, M. Do Carmo, S. Gupta, D.L. Rimm
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): M. Gale, Y. Li, J. Cao, Z.Z. Liu, S.M. Lang, S. Gupta, K. Aoshima, M.P. DiGiovanna, G.S. Shadel, Q. Yan
Writing, review, and/or revision of the manuscript: M. Gale, Y. Li, J. Cao, M. Do Carmo, S. Gupta, M.P. DiGiovanna, D.F. Stern, X. Chen, Q. Yan
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): S.M. Lang, L. Wu, D.L. Rimm
Study supervision: D.L. Rimm, X. Chen, Q. Yan
Acknowledgments
We thank the members of the Yan, Shadel, Rimm, Stern, Wajapeyee, and Nguyen laboratories at Yale School of Medicine for their helpful discussions and support. The Illumina sequencing service was conducted by Dr. Mei Zhong at Yale Stem Cell Center Genomics Core facility, which was supported by the Connecticut Regenerative Medicine Research Fund and the Li Ka Shing Foundation. We thank Dr. Wei Wei at Yale School of Public Health for his biostatistical support. This work was supported in part by the NIH (R21 CA187862 to Q. Yan; R01 CA216101 to G.S. Shadel; F32 AG052995 to M.A. Holmbeck; P30 CA016359 to Yale Comprehensive Cancer Center), the Natural Science Foundation of China (Major International Joint Research Program 81620108024 to X. Chen), a National Science Foundation Graduate Research Fellowship (DGE-1122492 to M. Gale), a James Hudson Brown – Alexander Brown Coxe Postdoctoral Fellowship (to M. Zhang), and a Leslie H. Warner Postdoctoral Fellowship (to L. Wu).
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.