Next Article in Journal
Neoantigen Dendritic Cell Vaccination Combined with Anti-CD38 and CpG Elicits Anti-Tumor Immunity against the Immune Checkpoint Therapy-Resistant Murine Lung Cancer Cell Line LLC1
Next Article in Special Issue
Current Understanding of Hypoxia in Glioblastoma Multiforme and Its Response to Immunotherapy
Previous Article in Journal
Targeted Deep Sequencing of Mycosis Fungoides Reveals Intracellular Signaling Pathways Associated with Aggressiveness and Large Cell Transformation
Previous Article in Special Issue
Impact of Hypoxia on Relative Biological Effectiveness and Oxygen Enhancement Ratio for a 62-MeV Therapeutic Proton Beam
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 13
Issue 21
10.3390/cancers13215509
cancers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Post-Hypoxic Cells Promote Metastatic Recurrence after Chemotherapy Treatment in TNBC

by
Inês Godet
1,2,3,
Mahelet Mamo
1,
Andrea Thurnheer
2,
D. Marc Rosen
1 and
Daniele M. Gilkes
1,2,3,4,*
1
Department of Oncology, The Sidney Kimmel Comprehensive Cancer Center, The Johns Hopkins University School of Medicine, Baltimore, MD 21231, USA
2
Department of Chemical and Biomolecular Engineering, The Johns Hopkins University, Baltimore, MD 21218, USA
3
Johns Hopkins Institute for NanoBioTechnology, The Johns Hopkins University, Baltimore, MD 21218, USA
4
Cellular and Molecular Medicine Program, The Johns Hopkins University School of Medicine, Baltimore, MD 21231, USA
*
Author to whom correspondence should be addressed.
Cancers 2021, 13(21), 5509; https://doi.org/10.3390/cancers13215509
Submission received: 22 September 2021 / Revised: 25 October 2021 / Accepted: 29 October 2021 / Published: 2 November 2021
(This article belongs to the Special Issue Hypoxia and Cancer: From Bench to Bedside)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

Intratumoral hypoxia is a negative prognostic factor in breast cancer progression and recurrence. By implementing a hypoxia fate-mapping system, we followed cells that experience intratumoral hypoxia in vivo and determined that these cells have an increased ability to metastasize compared to cells that were never exposed to hypoxia. In this work, we investigate whether cells that experienced intratumoral hypoxia are also resistant to chemotherapy. By utilizing both in vivo and ex vivo models, we conclude that metastatic cells found in the lung and liver, that were exposed to hypoxia in the primary tumor, are less sensitive to doxorubicin and paclitaxel and drive recurrence after treatment. Our studies also suggest that chemoresistance is associated with a cancer stem cell-like phenotype that is maintained in post-hypoxic cells.

Abstract

Hypoxia occurs in 90% of solid tumors and is associated with treatment failure, relapse, and mortality. HIF-1α signaling promotes resistance to chemotherapy in cancer cell lines and murine models via multiple mechanisms including the enrichment of breast cancer stem cells (BCSCs). In this work, we utilize a hypoxia fate-mapping system to determine whether triple-negative breast cancer (TNBC) cells that experience hypoxia in the primary tumor are resistant to chemotherapy at sites of metastasis. Using two orthotopic mouse models of TNBC, we demonstrate that cells that experience intratumoral hypoxia and metastasize to the lung and liver have decreased sensitivity to doxorubicin and paclitaxel but not cisplatin or 5-FU. Resistance to therapy leads to metastatic recurrence caused by post-hypoxic cells. We further determined that the post-hypoxic cells that metastasize are enriched in pathways related to cancer stem cell gene expression. Overall, our results show that even when hypoxic cancer cells are reoxygenated in the bloodstream they retain a hypoxia-induced cancer stem cell-like phenotype that persists and promotes resistance and eventually recurrence.

Graphical Abstract

1. Introduction

One in every eight women will develop breast cancer in their lifetime, and it has been projected that breast cancer will cause over 42,000 deaths this year [1]. The estimated 5-year survival rate for patients with localized breast cancer is 99%; however, this number dramatically drops to 27% when distant metastasis is detected [1]. Breast cancer metastasis occurs when cancer cells spread from the primary tumor through the body, typically to the lung, liver, bone, or brain [2]. It is estimated that 30% of patients with advanced breast cancer will develop distant metastasis [3], which is the main cause of treatment failure and mortality, leading to 90% of cancer-related deaths [4]. Even in patients with tumors that are responsive to radiotherapy or chemotherapy, metastasis is frequently a cause of treatment failure [5]. Metastatic tumors are anatomically less accessible, which makes conventional treatment options such as surgery or radiotherapy extremely challenging [6]. Moreover, metastatic lesions may arise from the clonal evolution of selected aggressive cancer cells that overcome a series of obstacles to thrive at a distant site [7]. In breast cancer, these cells are often pre-exposed to systemic therapy used to treat the primary tumor or to prevent recurrence [8]. Together, these circumstances make metastatic disease nearly incurable [9].
The standard of care for triple-negative metastatic breast cancer (TNMBC) is typically chemotherapy. The main cytotoxic agents used against TNMBC primarily target cell division. For instance, cisplatin binds to the N7 reactive center on purine residues, causing DNA damage in cancer cells, blocking cell division, and resulting in apoptotic cell death [10]. Doxorubicin can act on cancer cells by two main mechanisms: (1) intercalation into DNA and disruption of topoisomerase-II-mediated DNA repair and (2) generation of free radicals that damage cellular membranes, DNA, and proteins [11]. 5-Fluorouracil (5-FU) is converted to fluorodeoxyuridine monophosphate. This forms a stable complex with thymidylate synthase and inhibits deoxythymidine monophosphate (dTMP) production, which is vital for DNA replication and repair, leading to cytotoxicity [12]. Paclitaxel is a microtubule-stabilizing drug that induces mitotic arrest, leading to cell death [13]. The overall contribution of cytotoxic chemotherapy to the 5-year survival rate of adult cancer patients is estimated to be 2.1% in the United States, suggesting that this nontargeted approach has a low success rate [14].
The rapid proliferation of cancer cells along with the lack of efficient vasculature causes regions of hypoxia to develop in 90% of solid tumors [15,16]. Studies on pretreatment oxygenation status of solid tumors revealed that the median partial pressure of oxygen (PO2) in breast tumors is 10 mmHg (1.3% O2) in contrast to 65 mmHg (8% O2) in normal breast tissue [17]. Hypoxic PO2 values (<10 mmHg) have been correlated with increased risk of metastasis and mortality [18,19,20]. The most well-characterized response to oxygen deprivation is via the stabilization of the hypoxia-inducible subunits (HIF-1α and HIF-2α) [21]. Stabilized HIF-1/2α subunits heterodimerize with HIF-1β to form the HIF-1/2 complexes that recognize a 5′-ACGTG-3′ enhancer sequence and transcriptionally regulate more than a thousand genes [22,23]. Hypoxia-induced HIF-1/2 activation has been implicated in multiple mechanisms of chemoresistance, including ABC transporter activity, autophagy, DNA repair activation, enrichment in cancer stem cell-like properties, and repression of apoptosis and senescence [24].
We previously developed a hypoxia fate-mapping system that triggers a permanent fluorescent switch from DsRed to GFP when cells experience intratumoral hypoxia leading to increased HIF-1/2 expression. Our studies using breast cancer murine models revealed that GFP+ cells are 5 times more likely to successfully establish a metastatic lung lesion by adopting a ROS-resistant phenotype [25]. Given that hypoxia has been widely associated with chemoresistance, we sought to investigate whether these potent metastatic cells (hereinafter called post-hypoxic cells) also have decreased sensitivity to standard-of-care chemotherapeutic agents. In this work, we use both ex vivo and in vivo approaches to show that post-hypoxic cells (GFP+) are more resistant to doxorubicin and paclitaxel, but not cisplatin or 5-FU, and are more likely to contribute to recurrence following treatment.

2. Materials and Methods

2.1. Cell Culture

Mycoplasma-free breast cancer cell lines MDA-MB-231 and 4T1 were obtained from the American Type Culture Collection (ATCC) and maintained in DMEM and RPMI (Sigma-Aldrich, St. Louis, MO, USA), respectively, with 10% FBS (Corning, Corning, NY, USA) and 1% penicillin/streptomycin (Invitrogen, Carlsbad, CA, USA).

2.2. Hypoxia Fate-Mapping System

We previously generated and characterized the MDA-MB-231 and 4T1 hypoxia fate-mapping cell lines [25]. The vectors used to generate these cell lines, CMV-loxp-DsRed-loxp-eGFP (#141148) and 4xHRE-MinTK-CRE-ODD (#141147), are available at Addgene.

2.3. Tumor Implantation and Surgical Resection

Female 6- to 8-week-old NOD-SCID Gamma (NSG) mice were utilized according to protocols approved by the Johns Hopkins University Animal Care and Use Committee. Mice were anesthetized with 100 mg/kg ketamine and 16 mg/kg xylazine (Vet One) via intraperitoneal injection. MDA-MB-231 (2 × 106) or 4T1 (500) hypoxia fate-mapping cells were injected into the mammary fat pad under the second left nipple. Tumors were resected 22 to 25 days after tumor implantation. Briefly, mice were anesthetized and transferred to a heating pad. After hair removal, the tumor was carefully detached from adjacent skin, removing adjacent loose tissue and visible lymph nodes, and the wound was closed using 9 mm autoclips (Braintree Scientific, Inc.; Braintree, MA, USA). Betadine was applied to the wound and eye desiccation was prevented by using ophthalmic ointment. Resected tumors were divided in half. Half of the tumor was processed to obtain a cell suspension. Tumor tissue was physically and enzymatically (collagenase 2 mg/mL; MilliporeSigma; Burlington, MA, USA) dissociated, passed through a cell strainer (0.7 μM), and washed. Samples were resuspended in FACS buffer (PBS, 1% BSA, 0.5 mM EDTA, 25 μg/mL DNAse) and analyzed using an SH800 (Sony; San Jose, CA, USA) flow cytometer. The other half of the tumor was formalin-fixed, saturated in 30% sucrose (Sigma-Aldrich; St. Louis, MO, USA) at 4 °C overnight, embedded in OCT media (Fisher Scientific; Hampton, NH, USA), frozen in liquid nitrogen, sectioned via a cryotome CM1100 (Leica; Wetzlar, Germany), and mounted onto Superfrost Plus Microscope Slides (Fisher Scientific). Tumor sections were stained with DAPI (1:1000 for 15 min, RT), mounted with antifade solution, and imaged using an Olympus (UPLFLN 4X) objective on a Cytation 5 instrument (BioTek Instruments; Winooski, VT, USA).

2.4. Generation of Tumor-Derived Cell Lines

Tumor-derived cells were obtained following the previously described tissue dissociation protocol with physical dissociation, followed by enzymatic digestion (collagenase 2 mg/mL (Sigma-Aldrich) and BSA 2 mg/mL (Gemini Bioproducts; West Sacramento, CA, USA) for 1 h at 37 °C and 160 RPM. The resultant cell suspension was washed with PBS, resuspended in sorting buffer (PBS, 1% BSA, 0.5 mM EDTA, 25 mM HEPES pH 8), and sorted into DsRed+/GFP− and DsRed−/GFP+ populations using an SH800 Cell Sorter (Sony Biotechnology). GFP was detected in the FITC channel, and DsRed was detected in the PE channel. Cells were immediately plated in tissue-culture plates with warm DMEM or RPMI media. Cells were mycoplasma-tested and maintained in culture for at most 5–6 passages.

2.5. Ex Vivo Drug Screening

Tumor-derived MDA-MB-231 and 4T1 DsRed+ or GFP+ cells were plated in 96-well plates (5000 cells/well). The following day, fresh media with cisplatin, doxorubicin, 5-FU, or paclitaxel were added. In order to obtain IC50 curves, 7 different drug concentrations and a control condition were evaluated for each chemotherapeutic agent. Following 48 h, the media were aspirated and media with 10% Presto Blue were added. Plates were incubated for 2 h and fluorescence was read using a Cytation5 (BioTek Instruments). The cells were then washed with PBS, fixed with 4% PFA, and stained with DAPI. The entire tissue culture well was imaged in the RFP, GFP, and DAPI channels with an Olympus (UPLFLN 4XPh) phase objective 4× in a montage (4 × 3) (BioTek Instruments). Images were linearly stitched, and the DAPI area was quantified by thresholding as a measure of cellular confluence. IC50 curves based on fluorescence measured by Presto Blue or cell confluence were obtained using the nonlinear fit log(inhibitor) vs. response model available in the GraphPad Prism software. GFP+ tumor-derived cells only retained a resistant phenotype against doxorubicin and paclitaxel for 5–6 passages.

2.6. Ex Vivo Colony Formation Assay

Tumor-derived MDA-MB-231 DsRed+ or GFP+ cells were plated in 96-well plates (5000 cells/well) and treated with the 10 nM of paclitaxel for 48 h, followed by a 30-day recovery period to allow resistant colonies to form. Media were refreshed twice a week, and whole wells were imaged in the RFP and GFP channels with an Olympus (UPLFLN 4XPh) phase objective 4× in a montage (4 × 3) (BioTek Instruments). Images were linearly stitched, and colonies were counted using the images acquired on day 20 of the experiment.

2.7. Chemotherapeutic Treatment of Mice

Chemotherapy treatment was administered to animals via intravenous injection. Briefly, mice were warmed for 5 min with an overhead heat lamp to dilate the veins, while applying friction to the tail. Each animal was restrained, and the dose of chemotherapy (3 mg/kg cisplatin, 2 mg/kg doxorubicin, 20 mg/kg 5-FU, or 16 mg/kg paclitaxel) was injected in appropriate vehicle solution using a 26 G needle. The control group was treated with vehicle solution only.

2.8. Analysis of Metastatic Progression

At the end point of the experiment, animals were sacrificed to assess metastatic burden. Lungs were inflated with an OCT–PBS solution, excised, and immediately formalin-fixed. Livers were excised and fixed. Whole mount fluorescent imaging of lungs and livers was conducted using an Olympus (UPLFLN 1.25XPh) phase objective 1.25× (BioTek Instruments) in the RFP and GFP channels. Images were linearly stitched. The quantification of the metastatic burden was performed using flow cytometry analysis. Briefly, lungs and livers were chopped and digested in 2 mg/mL collagenase (Sigma-Aldrich) and 2 mg/mL BSA (Gemini Bioproducts) for 1 h at 37 °C at 160 RPM. After passing through a 0.70 μm strainer, cells were washed with PBS and resuspended in FACS buffer. GFP and DsRed were detected in the FITC and PE channels, respectively, using an SH800 (Sony Biotechnology). Data were analyzed using the FlowJo V10 software (Tree Star Inc., Ashland, OR, USA).

2.9. Computational Analysis

RNA sequencing was previously performed on DsRed+ and GFP+ cells sorted from tumors and lungs of mice [25]. Gene set enrichment analysis (GSEA) [26] was conducted to determine enrichment in DsRed or GFP cells derived from tumors or lungs against gene sets associated with chemoresistance and resistance to cisplatin, doxorubicin, or 5-FU available in the GSEA database and a paclitaxel-resistance gene signature compiled as described in Table S1. Survival curves were plotted using the Kaplan–Meier Plotter (https://kmplot.com/analysis/; accessed on 26 April 2021) online tool using mRNA gene chip data from breast cancer patients [27]. ROC analysis was conducted using the ROC Plotter (http://www.rocplot.org/; accessed on 9 April 2021) online tool [28].

2.10. Mammosphere Formation Assay

Freshly resected MDA-MB-231 tumors were subjected to physical dissociation followed by enzymatic digestion (collagenase 2 mg/mL (Sigma-Aldrich), BSA 2 mg/mL (Gemini Bioproducts)) for 1 h at 37 °C at 160 RPM. Cells were then passed through a cell strainer (0.7 μM), and 10,000 cells/well were transferred to a 6-well plate previously coated with polyHEMA (12 g/L in 95% EtOH) and air-dried for 48 h with 2 mL of mammosphere formation media per well (MammoCult basal medium (human) with 4 ug/mL heparin, and 0.48 ug/mL hydrocortisone) (STEMCELL Technologies; Vancouver, BC, Canada).

2.11. Statistical Analysis

All data are presented as mean ± SEM, and statistical analysis was performed using GraphPad Prism 9. Differences in metastatic burden were analyzed via ordinary one-way ANOVA with Dunnett (recommended) multiple comparison test or paired one-tailed Student’s t-test when comparing DsRed versus GFP. A comparison of DsRed+ versus GFP+ in matched samples such as the mammosphere assay was performed via paired two-tailed Student’s t-test. Significance levels are reported by displaying exact p-values.

3. Results

3.1. Post-Hypoxic Cells Are Less Sensitive to Doxorubicin and Paclitaxel but Not 5-FU or Cisplatin Treatment Ex Vivo

We recently developed a dual-vector lentiviral system that when expressed in cells causes a permanent fluorescent switch from DsRed to GFP expression under hypoxia. The first vector contains a constitutively active promoter that drives the expression of DsRed. The DsRed gene is flanked by loxp sites and is localized upstream of the GFP gene. The second vector is activated by HIF binding to the hypoxia-responsive elements (HREs) located upstream of a minimal promoter that causes the expression of an altered Cre protein that contains an oxygen-dependent domain (ODD) (Figure 1a). When a cell experiences hypoxia, Cre accumulation drives the cleavage of DsRed and consequently the permanent expression of GFP. By incorporating this system in an MDA-MB-231 orthotopic mouse model of breast cancer metastasis, we previously reported that post-hypoxic cancer cells (GFP+) have a survival advantage in the bloodstream leading to a 5-fold increase in their ability to metastasize to the lung (Figure 1a) [25].
Intratumoral hypoxia has been associated with multiple mechanisms of chemoresistance [24]. We previously performed an RNA sequencing analysis to characterize DsRed+ and GFP+ cancer cells derived from mouse tumor and lung metastasis. Using gene set enrichment analysis (GSEA) [26], we determined that GFP+ cells were enriched in a gene signature associated with chemoresistance [29] (Figure S1a–c). To further explore the impact of hypoxia in resistance, we assessed the effect of four standard-of-care chemotherapies used to treat patients with breast cancer utilizing an ex vivo approach (Figure 1b). DsRed+ and GFP+ MDA-MB-231 cells were sorted from orthotopic tumors harvested 25 days after cancer cell inoculation. After expansion, the IC50 for cisplatin, doxorubicin, 5-fluorouracil, and paclitaxel was determined by measuring fluorescence (RFU) using a Presto Blue proliferation assay, as well as by determining the percentage of cell confluence by analyzing the percentage of DAPI staining. The results showed that cells that experienced intratumoral hypoxia (GFP+) have a 2-fold higher IC50 for doxorubicin and 3-fold higher IC50 for paclitaxel, whereas the IC50 values for cisplatin or 5-FU treatments did not change (Figure 1c–h and Figure S1d–f). GFP+ cells sorted from 4T1 orthotopic tumors also displayed decreased sensitivity compared to DsRed+ cells, particularly to paclitaxel treatment (Figure S2).

3.2. Post-Hypoxic Cells Are Less Sensitive to Doxorubicin and Paclitaxel but Not Cisplatin or 5-FU Treatment In Vivo

To determine whether the sensitivity of GFP+ cells to chemotherapy was also similar in vivo, hypoxia fate-mapping MDA-MB-231 cells were implanted into the mammary fat pad of NSG mice. Tumors were resected 25 days later following a neoadjuvant dose of chemotherapy. We previously demonstrated [25] and confirmed here that at the time of tumor removal about 20% of the cancer cells in the primary tumor had experienced intratumoral hypoxia and were GFP+ (<1% O2) (Figure S3a). Five additional doses of chemotherapy were administered daily after surgery, and the mouse weight was monitored over the course of treatment (Figure 2a and Figure S3b,c). Treatment with cisplatin or 5-FU did not reduce metastatic burden in this model, while doxorubicin and paclitaxel reduced or prevented the outgrowth of lung and liver metastasis by 88–90% and 75%, respectively (Figure 2b,c). However, mice treated with doxorubicin or paclitaxel had almost twice the fraction of GFP+ compared to DsRed+ cells in the lung and 3 to 4 times more GFP+ compared to DsRed+ cells in the liver at the endpoint of the experiment (Figure 2d,e). Taken together the data demonstrated that there are 8 and 14 times more GFP+ than DsRed+ cells present in the metastatic lung and liver lesions, respectively, than there were in the primary tumor at the time of removal.
To ensure that the increase in GFP+ cells was not due to paclitaxel promoting HIF-1α activation [30], we performed a control experiment. We injected hypoxia fate-mapping MDA-MB-231 cells into the tail vein of mice in order to directly seed the mouse lung without forming a primary tumor. The mice were treated with five daily doses of paclitaxel. The colonies that formed in the lung were exclusively DsRed+ (Figure S3d). This confirms that paclitaxel-driven HIF1α expression is not sufficient to promote a DsRed-to-GFP switch in our system. Likewise, the DsRed-to-GFP switch occurs in the primary tumor and is maintained when the cells are reoxygenated in the lung, as we previously reported.
By utilizing GSEA analysis of GFP+ compared to DsRed+ cells, we investigated previously published signatures associated with resistance to cisplatin [31], doxorubicin [31], 5-FU [31], and paclitaxel [32,33,34,35] (Table S1). GFP+ cancer cells sorted from either the primary tumor or lung were significantly enriched in genes associated with resistance to doxorubicin, paclitaxel, and cisplatin, but not 5-FU (Figure 2f and Figure S4). Taken together, these data (1) recapitulate the results obtained using the ex vivo screening assay and (2) demonstrate that cells that experience intratumoral hypoxia (GFP+) and then metastasize have decreased sensitivity to doxorubicin and paclitaxel.

3.3. Post-Hypoxic Cells Are Less Sensitive to Doxorubicin and Paclitaxel

Next, we performed a follow-up experiment with a larger cohort of mice using the same treatment regimen consisting of a neoadjuvant dose of chemotherapy, tumor resection, and four consecutive doses of doxorubicin or paclitaxel (Figure 3a–c and Figure S5a–c). The results replicated our initial findings showing a higher fraction of GFP+ versus DsRed+ cancer cells in the lung and liver of mice treated with chemotherapy versus the vehicle control (Figure 3d–e). To validate our findings in a second model of breast cancer metastasis, we performed the same study using 4T1 hypoxia fate-mapping cells. Tumors were resected at day 22 when they were comprised of about 15% GFP+ cells (Figure S5d). Animals were monitored over the course of the experiment (Figure S5e,f). Paclitaxel treatment reduced metastatic lung burden by 45%. The fraction of GFP+ cells in the lung was 2-fold more than DsRed cells (Figure S5g–i). Together, these data confirm that cells that experience intratumoral hypoxia (GFP+) and metastasize to distant organs are less sensitive to paclitaxel.
In our previous work, we identified a 19-gene signature upregulated in GFP+ cells derived from both tumor and lung metastases that we termed ‘hypoxic memory’ genes [25] (Figure 3f). By utilizing this signature, we compared the distant metastasis-free survival (DMFS) of patients treated with chemotherapy (N = 225) or endocrine therapy (N = 205) [27]. The patients were stratified as either having high or low expression of the 19-gene signature using median expression as the cutoff. Patients treated with chemotherapy that had high expression of the 19-gene signature had a worse outcome (Figure 3g). On the other hand, the expression of this signature was not prognostic for patients treated with endocrine therapy (Figure 3h). These observations suggest that the post-hypoxic phenotype could be associated with reduced sensitivity to chemotherapy.

3.4. Post-Hypoxic Cells Are More Likely to Contribute to Metastatic Recurrence after Treatment

Our initial findings revealed that MDA-MB-231 cells that experience hypoxia in the primary tumor were less sensitive to doxorubicin and paclitaxel both ex vivo and in vivo after metastasizing to the lung and liver. Next, we questioned whether post-hypoxic cells would be more likely to contribute to recurrence after treatment. In a preliminary assessment, we used an ex vivo approach. Tumor-derived DsRed+ or GFP+ cells were treated for 48 h with 10 nM of paclitaxel, followed by a recovery period of 30 days. GFP+ cancer cells formed 5 times more colonies than DsRed+ cells (Figure 4a,b). To test the ability of the cancer cells to recur and cause relapse in vivo, we sacrificed paclitaxel-treated mice on day 42, one week after the last dose of treatment (PAX-II), whereas the vehicle control group was sacrificed on day 35 (Figure 4c). At the time of removal, the tumors contained approximately 20% GFP+ cells (Figure S5j). The mouse body weight was monitored over the course of the experiment (Figure S5k,l). The total metastatic burden in mice treated with paclitaxel was nearly equivalent to that of the control group, suggesting that the population of micrometastatic cells that remain after treatment are sufficient to cause relapse (Figure 4d–f). A higher fraction of GFP+ cancer cells was found at both metastatic sites when compared to the control group (Figure 4g,h).
Taken together, the data demonstrate that even though paclitaxel is effective at drastically reducing metastatic burden in the lung and liver (Figure 2b and Figure 3b), cells that remain at these sites can quickly recover and proliferate once the treatment stops. Importantly, our results demonstrate that even though the primary tumor was composed of only 20% GFP+ cells at the time of tumor resection, GFP+ cells contributed 2–3 times more than DsRed+ cells to metastatic recurrence in the lung and liver.

3.5. Post-Hypoxic Cells Develop a Breast Cancer Stem Cell-like Phenotype

Hypoxia has been implicated in chemoresistance through several HIF-1α dependent mechanisms, including upregulation of ABC transporters, DNA repair mechanisms, autophagy, enrichment in cancer stem cell-like properties, and protection from apoptosis and senescence [24]. Previous studies have demonstrated that HIF-1α is required for CSC enrichment and CSCs are more chemoresistant than cells that do not have a CSC phenotype [30,36,37]. To investigate the cancer stem cell-like potential of cells that experienced intratumoral hypoxia (GFP+), we harvested cancer cells from recently harvested tumors and cultured them in mammosphere formation media following reoxygenation (Figure 5a). After 7 days, we determined that GFP+ cells generated 2 times more mammospheres and were larger in size than DsRed cells (Figure 5b,c).
To further investigate the cancer stem cell-like phenotype in the GFP+ cells harvested from lung metastasis, we used a GSEA analysis. Interestingly, GFP+ cells were significantly enriched for a stemness gene-signature (Zhang et al., 2008 [38]) compared to DsRed+ cells (Figure 5d). The top 10 genes driving this result visibly clustered the DsRed+ versus the GFP+ cancer cells derived from the lung (Figure 5e). By utilizing the top 10 stemness genes, we performed receiver operating characteristic (ROC) analysis to determine if this signature could predict benefit from a specific treatment [28]. Our results showed that patients with breast cancer and high expression of this signature were less likely to respond to taxane therapy (Figure 5f). In contrast, the signature could not predict response to endocrine; anthracycline; or fluorouracil, Adriamycin (doxorubicin), and Cytoxan (FAC) treatment regimens (Figure 5g–i and Figure S6). Together these results suggest that hypoxia enriches the expression of genes related to stemness, which may promote the chemoresistant phenotype of GFP+ cells.

4. Discussion

Metastasis is the cause of 90% of solid cancer-related deaths, and this is in part due to chemoresistance. Furthermore, intratumoral hypoxia has been widely associated with chemoresistance via several mechanisms that depend on HIF-signaling, including HIF-transcriptional activation of ATP-binding cassette (ABC) transporter family genes to promote drug efflux [39,40], hypoxia-dependent autophagy via HIF-1 to promote resistance [41,42], and HIF-1-regulated decreased drug-induced senescence [43]. In our previous work, we designed and implemented a hypoxia fate-mapping system in MDA-MB-231 and 4T1 cells and determined that cells that experience intratumoral hypoxia have a 5-fold increase in their ability to form lung metastasis [25]. In the current study, we ask whether cells that experience hypoxia in the primary tumor and become reoxygenated (termed post-hypoxic cells) as they metastasize to other organs via the bloodstream are resistant to chemotherapy. Gene set enrichment analysis (GSEA) was used as an assessment tool to demonstrate that cells that experienced intratumoral hypoxia showed enrichment of previously published chemoresistance signatures. Moreover, our studies using both ex vivo and in vivo approaches revealed that post-hypoxic MDA-MB-231 cells show decreased sensitivity to paclitaxel and doxorubicin. In line with our observations, in vitro studies have demonstrated that MDA-MB-231 cells exposed to 1% O2 showed increased resistance to paclitaxel [44,45] and doxorubicin [43] via a HIF-1α-dependent mechanism. This suggests that the post-hypoxic resistance observed in GFP+ cells might be linked to HIF-1α signaling, which we have previously demonstrated to be enriched in GFP+ cancer cells isolated from both mouse tumors and lungs [25]. The GFP+ MDA-MB-231 cells did not show resistance to cisplatin or 5-FU. The 4T1 post-hypoxic cells were also less sensitive to paclitaxel in vitro and in vivo. One caveat of our study is that it was conducted using only MDA-MB-231 and 4T1 cell lines. Thus, our results do not rule out that in other cell lines, cancer types, or patients, hypoxia also alters sensitivity to cisplatin or 5-FU. For example, one report demonstrated hypoxia-induced cisplatin resistance in MDA-MB-231 cells in vitro [46], which was not observed in our set-up. This difference in results might suggest that resistance to cisplatin is an effect of active hypoxia as cells in the aforementioned study were exposed and treated for 24 and 48 h, whereas our study is focused on re-oxygenated cells in the lung. In another study, hypoxia-induced resistance to 5-FU was noted in gastric cancer cell lines (AGS and MKN28) [47].
HIF-1α signaling has been implicated in several mechanisms of chemoresistance, particularly enrichment in a cancer stem cell-like phenotype. BCSCs have enhanced resistance to cytotoxic agents, and the percentage of BCSCs following chemotherapy is increased compared with their before treatment status [30]. By utilizing GSEA, our results demonstrate that cells that experience intratumoral hypoxia and form distant metastasis in the lung have enrichment in hallmark hypoxia and cancer stem cell-like signaling pathways. Moreover, our preliminary assessment suggests that tumor-derived GFP+ cells have increased cancer stem cell-like properties and form more mammospheres, a well-accepted phenotype of BCSCs.
In our previous work, we also determined that cells that experience intratumoral hypoxia have a ROS-resistant phenotype that enhances survival in the bloodstream and ultimately drives metastasis. Chemotherapy can promote ROS accumulation, overwhelming redox homeostasis and leading to ROS-induced apoptosis [48]. Both paclitaxel and doxorubicin promote ROS accumulation [49,50]. Thus, it is also plausible that the ROS-resistant phenotype of the GFP+ cells may also contribute to the observed resistance.

5. Conclusions

In this work, we used a hypoxia fate-mapping system to investigate whether intratumoral hypoxia causes chemoresistance in murine models of breast cancer. Overall, our results show that cells that experience hypoxia in the primary tumor are (1) more equipped to metastasize, (2) less sensitive to chemotherapeutics, and (3) more likely to cause recurrence. Biomarkers are urgently needed to identify post-hypoxic cells in primary tumors and at metastatic sites in order to develop a strategy to target them and prevent recurrence.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/cancers13215509/s1, Figure S1: Post-hypoxic tumor cells are less sensitive to chemotherapy ex vivo, Figure S2: 4T1 post-hypoxic tumor cells are resistant to Paclitaxel ex vivo, Figure S3: Ex vivo screening of post-hypoxic tumor cells predicts treatment outcome in vivo, Figure S4: GFP+ cells are enriched in signatures associated with chemoresistance, Figure S5: GFP+ cancer cells are resistant to Doxorubicin and Paclitaxel in vivo, Figure S6: GFP+ cells that metastasize retain breast cancer stem cell phenotype, Table S1: Compiled gene signature of genes associated with resistance to Paclitaxel.

Author Contributions

Conceptualization, I.G. and D.M.G.; methodology, I.G. and D.M.G.; formal analysis, I.G.; investigation, I.G., M.M., A.T. and D.M.R.; data curation, I.G.; writing—original draft preparation, I.G.; writing—review and editing, D.M.G.; supervision, D.M.G.; funding acquisition, D.M.G. All authors have read and agreed to the published version of the manuscript.

Funding

Research in the Gilkes lab is funded by NCI grant number U54-CA210173, The Jayne Koskinas Ted Giovanis Foundation for Health and Policy, Cindy Rosencrans Fund for Metastatic Triple-Negative Breast Cancer, The Emerson Collective, and the NCI/SKCCC Core grant number P50CA006973.

Institutional Review Board Statement

Animal studies were conducted according to protocols approved by the Johns Hopkins University Animal Care and Use Committee (JHMI OLAW Animal Welfare Assurance number: D16-00173).

Informed Consent Statement

Not applicable.

Data Availability Statement

RNA sequencing data are available in the GEO database under the accession codes GSE126609 (tumor) and GSE136372 (lung).

Acknowledgments

Cartoons included in this work’s figures were adapted from Servier Medical Art licensed under a Creative Commons Attribution 3.0 Unported License (smart.servier.com accessed on 26 April 2021).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Siegel, R.L.; Miller, K.D.; Jemal, A. Cancer Statistics, 2020. CA A Cancer J. Clin. 2020, 70, 7–30. [Google Scholar] [CrossRef]
  2. Gilkes, D.M.; Semenza, G.L. Role of hypoxia-inducible factors in breast cancer metastasis. Future Oncol. 2013, 9, 1623–1636. [Google Scholar] [CrossRef] [Green Version]
  3. O’Shaughnessy, J. Extending Survival with Chemotherapy in Metastatic Breast Cancer. Oncologist 2005, 10, 20–29. [Google Scholar] [CrossRef] [Green Version]
  4. Mehlen, P.; Puisieux, A. Metastasis: A question of life or death. Nat. Rev. Cancer 2006, 6, 449–458. [Google Scholar] [CrossRef]
  5. Qiu, W.Z.; Huang, P.Y.; Shi, J.L.; Xia, H.Q.; Zhao, C.; Cao, K.J. Neoadjuvant chemotherapy plus intensity-modulated radiotherapy versus concurrent chemoradiotherapy plus adjuvant chemotherapy for the treatment of locoregionally advanced nasopharyngeal carcinoma: A retrospective controlled study. Chin. J. Cancer 2016, 35, 1–9. [Google Scholar] [CrossRef] [Green Version]
  6. Qian, C.N.; Mei, Y.; Zhang, J. Cancer metastasis: Issues and challenges. Chin. J. Cancer 2017, 36, 1–4. [Google Scholar] [CrossRef]
  7. Turajlic, S.; Swanton, C. Metastasis as an evolutionary process. Science 2016, 352, 169–175. [Google Scholar] [CrossRef]
  8. Ibragimova, M.K.; Tsyganov, M.M.; Litviakov, N.V. Natural and chemotherapy-induced clonal evolution of tumors. Biochem. (Mosc.) 2017, 82, 413–425. [Google Scholar] [CrossRef]
  9. Nicolini, A.; Giardino, R.; Carpi, A.; Ferrari, P.; Anselmi, L.; Colosimo, S.; Conte, M.; Fini, M.; Giavaresi, G.; Berti, P.; et al. Metastatic breast cancer: An updating. Biomed. Pharmacother. 2006, 60, 548–556. [Google Scholar] [CrossRef]
  10. Dasari, S.; Tchounwou, P.B. Cisplatin in cancer therapy: Molecular mechanisms of action. Eur. J. Pharmacol. 2014, 740, 364–378. [Google Scholar] [CrossRef] [Green Version]
  11. Thorn, C.F.; Oshiro, C.; Marsh, S.; Hernandez-Boussard, T.; McLeod, H.; Klein, T.E.; Altman, R.B. Doxorubicin pathways: Pharmacodynamics and adverse effects. Pharmacogenet. Genom. 2011, 21, 440. [Google Scholar] [CrossRef]
  12. Zhang, N.; Yin, Y.; Xu, S.J.; Chen, W.S. 5-Fluorouracil: Mechanisms of resistance and reversal strategies. Molecules 2008, 13, 1551–1569. [Google Scholar] [CrossRef] [Green Version]
  13. Weaver, B.A. How Taxol/paclitaxel kills cancer cells. Mol. Biol. Cell 2014, 25, 2677–2681. [Google Scholar] [CrossRef]
  14. Morgan, G.; Ward, R.; Barton, M. The contribution of cytotoxic chemotherapy to 5-year survival in adult malignancies. Clin. Oncol. 2004, 16, 549–560. [Google Scholar] [CrossRef]
  15. Vaupel, P.; Mayer, A.; Hockel, M. Tumor hypoxia and malignant progression. Methods Enzymol. 2004, 23, 335–354. [Google Scholar]
  16. Cosse, J.-P.; Michiels, C. Tumour hypoxia affects the responsiveness of cancer cells to chemotherapy and promotes cancer progression. Anti-Cancer Agents Med. Chem. 2008, 8, 790–797. [Google Scholar] [CrossRef]
  17. Vaupel, P.; Hockel, M.; Mayer, A. Detection and characterization of tumor hypoxia using pO2 histography. Antioxid. Redox Signal. 2007, 9, 1221–1235. [Google Scholar] [CrossRef]
  18. Vaupel, P.; Mayer, A.; Briest, S.; Höckel, M. Oxygenation Gain Factor: A Novel Parameter Characterizing the Association between Hemoglobin Level and the Oxygenation Status of Breast Cancers. Cancer Res. 2003, 63, 7634–7637. [Google Scholar] [PubMed]
  19. Vaupel, P.; Mayer, A. Hypoxia in cancer: Significance and impact on clinical outcome. Cancer Metastasis Rev. 2007, 26, 225–239. [Google Scholar] [CrossRef]
  20. Walsh, J.C.; Lebedev, A.; Aten, E.; Madsen, K.; Marciano, L.; Kolb, H.C. The Clinical Importance of Assessing Tumor Hypoxia: Relationship of Tumor Hypoxia to Prognosis. Antioxid. Redox Signal. 2014, 21, 1516–1554. [Google Scholar] [CrossRef]
  21. Dengler, V.L.; Galbraith, M.D.; Espinosa, J.M. Transcriptional regulation by hypoxia inducible factors. Crit. Rev. Biochem. Mol. Biol. 2014, 49, 1–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Wang, G.L.; Semenza, G.L. Purification and characterization of hypoxia-inducible factor. J. Biol. Chem. 1995, 270, 1230–1237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Ye, I.C.; Fertig, E.J.; DiGiacomo, J.W.; Considine, M.; Godet, I.; Gilkes, D.M. Molecular portrait of hypoxia in breast cancer: A prognostic signature and novel HIF-regulated genes. Mol. Cancer Res. 2018, 16, 1889–1901. [Google Scholar] [CrossRef] [Green Version]
  24. Minassian, L.M.; Cotechini, T.; Huitema, E.; Graham, C.H. Hypoxia-induced resistance to chemotherapy in cancer. In Hypoxia and Cancer Metastasis; Springer: Cham, Switzerland, 2019; pp. 123–139. [Google Scholar]
  25. Godet, I.; Shin, Y.J.; Ju, J.A.; Ye, I.C.; Wang, G.; Gilkes, D.M. Fate-mapping post-hypoxic tumor cells reveals a ROS-resistant phenotype that promotes metastasis. Nat. Commun. 2019, 10, 1–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. 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. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [Green Version]
  27. Györffy, B.; Lanczky, A.; Eklund, A.C.; Denkert, C.; Budczies, J.; Li, Q.; Szallasi, Z. An online survival analysis tool to rapidly assess the effect of 22,277 genes on breast cancer prognosis using microarray data of 1809 patients. Breast Cancer Res. Treat. 2010, 123, 725–731. [Google Scholar] [CrossRef] [Green Version]
  28. Fekete, T.; Győrffy, B. ROCplot.org: Validating predictive biomarkers of chemotherapy/hormonal therapy/anti-HER2 therapy using transcriptomic data of 3104 breast cancer patients. Int. J. Cancer 2019, 145, 3140–3151. [Google Scholar] [CrossRef]
  29. Yuan, J.; Tan, L.; Yin, Z.; Tao, K.; Wang, G.; Shi, W.; Gao, J. Bioinformatics analysis identifies potential chemoresistance-associated genes across multiple types of cancer. Oncol. Lett. 2019, 18, 2576–2583. [Google Scholar] [CrossRef]
  30. Samanta, D.; Gilkes, M.; Chaturvedia, P.; Xiang, L.; Semenza, G.L. Hypoxia-inducible factors are required for chemotherapy resistance of breast cancer stem cells. Proc. Natl. Acad. Sci. USA 2014, 111, E5429–E5438. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Kang, H.C.; Kim, I.J.; Park, J.H.; Shin, Y.; Ku, J.L.; Jung, M.S.; Yoo, B.C.; Kim, H.K.; Park, J.G. Identification of Genes with Differential Expression in Acquired Drug-Resistant Gastric Cancer Cells Using High-Density Oligonucleotide Microarrays. Clin. Cancer Res. 2004, 10, 272–284. [Google Scholar] [CrossRef] [Green Version]
  32. Murakami, H.; Ito, S.; Tanaka, H.; Kondo, E.; Kodera, Y.; Nakanishi, H. Establishment of new intraperitoneal paclitaxel-resistant gastric cancer cell lines and comprehensive gene expression analysis. Anticancer Res. 2013, 33, 4299–4307. [Google Scholar]
  33. Duan, Z.; Lamendola, D.E.; Duan, Y.; Yusuf, R.Z.; Seiden, M.V. Description of paclitaxel resistance-associated genes in ovarian and breast cancer cell lines. Cancer Chemother. Pharmacol. 2005, 55, 277–285. [Google Scholar] [CrossRef]
  34. Dorman, S.N.; Baranova, K.; Knoll, J.H.; Urquhart, B.L.; Mariani, G.; Carcangiu, M.L.; Rogan, P.K. Genomic signatures for paclitaxel and gemcitabine resistance in breast cancer derived by machine learning. Mol. Oncol. 2016, 10, 85–100. [Google Scholar] [CrossRef]
  35. Jurj, A.; Pop, L.A.; Zanoaga, O.; Ciocan-Cârtiţă, C.A.; Cojocneanu, R.; Moldovan, C.; Raduly, L.; Pop-Bica, C.; Trif, M.; Irimie, A.; et al. New Insights in Gene Expression Alteration as Effect of Paclitaxel Drug Resistance in Triple Negative Breast Cancer Cells. Cell. Physiol. Biochem. 2020, 54, 648–664. [Google Scholar] [PubMed]
  36. Xiong, G.; Stewart, R.L.; Chen, J.; Gao, T.; Scott, T.L.; Samayoa, L.M.; O’Connor, K.; Lane, A.N.; Xu, R. Collagen prolyl 4-hydroxylase 1 is essential for HIF-1α stabilization and TNBC chemoresistance. Nat. Commun. 2018, 9, 1–16. [Google Scholar] [CrossRef] [Green Version]
  37. Cui, C.P.; Wong, C.C.L.; Kai, A.K.L.; Ho, D.W.H.; Lau, E.Y.T.; Tsui, Y.M.; Chan, L.K.; Cheung, T.T.; Chok, K.S.H.; Chan, A.C.; et al. SENP1 promotes hypoxia-induced cancer stemness by HIF-1α deSUMOylation and SENP1/HIF-1α positive feedback loop. Gut 2017, 66, 2149–2159. [Google Scholar] [CrossRef]
  38. Zhang, M.; Behbod, F.; Atkinson, R.L.; Landis, M.D.; Kittrell, F.; Edwards, D.; Medina, D.; Tsimelzon, A.; Hilsenbeck, S.; Green, J.E.; et al. Identification of tumor-initiating cells in a p53-null mouse model of breast cancer. Cancer Res. 2008, 68, 4674–4682. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Krishnamurthy, P.; Ross, D.D.; Nakanishi, T.; Bailey-Dell, K.; Zhou, S.; Mercer, K.E.; Sarkadi, B.; Sorrentino, B.P.; Schuetz, J.D. The stem cell marker Bcrp/ABCG2 enhances hypoxic cell survival through interactions with heme. J. Biol. Chem. 2004, 279, 24218–24225. [Google Scholar] [CrossRef] [Green Version]
  40. Comerford, K.M.; Wallace, T.J.; Karhausen, J.; Louis, N.A.; Montalto, M.C.; Colgan, S.P. Hypoxia-inducible Factor-1-dependent Regulation of the Multidrug Resistance (MDR1) Gene 1. Cancer Res. 2002, 62, 3387–3394. [Google Scholar] [PubMed]
  41. Wu, H.M.; Jiang, Z.F.; Ding, P.S.; Shao, L.J.; Liu, R.Y. Hypoxia-induced autophagy mediates cisplatin resistance in lung cancer cells. Sci. Rep. 2015, 5, 1–5. [Google Scholar] [CrossRef]
  42. Yang, X.; Yin, H.; Zhang, Y.; Li, X.; Tong, H.; Zeng, Y.; Wang, Q.; He, W. Hypoxia-induced autophagy promotes gemcitabine resistance in human bladder cancer cells through hypoxia-inducible factor 1α activation. Int. J. Oncol. 2018, 53, 215–224. [Google Scholar] [CrossRef] [Green Version]
  43. Sullivan, R.; Frederiksen, L.J.; Pare, G.C.; Graham, C.H. Hypoxia-Inducible Factor 1 is required for hypoxia-induced resistance to doxorubicin in breast carcinoma cells. Cancer Res. 2006, 66, 403–404. [Google Scholar]
  44. Flamant, A. Notte, N. Ninane, Raes, M.; Michiels, C. Anti-apoptotic role of HIF-1 and AP-1 in paclitaxel exposed breast cancer cells under hypoxia. Mol. Cancer 2010, 9, 1–15. [Google Scholar] [CrossRef] [Green Version]
  45. Notte, A.; Ninane, N.; Arnould, T.; Michiels, C. Hypoxia counteracts taxol-induced apoptosis in MDA-MB-231 breast cancer cells: Role of autophagy and JNK activation. Cell Death Dis. 2013, 4, e638. [Google Scholar] [CrossRef] [PubMed]
  46. Soleymani Abyaneh, H.; Gupta, N.; Radziwon-Balicka, A.; Jurasz, P.; Seubert, J.; Lai, R.; Lavasanifar, A. STAT3 but not HIF-1α is important in mediating hypoxia-induced chemoresistance in MDA-MB-231, a triple negative breast cancer cell line. Cancers 2017, 9, 137. [Google Scholar] [CrossRef] [Green Version]
  47. Rohwer, N.; Dame, C.; Haugstetter, A.; Wiedenmann, B.; Detjen, K.; Schmitt, C.A.; Cramer, T. Hypoxia-inducible factor 1α determines gastric cancer chemosensitivity via modulation of p53 and Nf-κB. PLoS ONE 2010, 5, e12038. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Perillo, B.; Di Donato, M.; Pezone, A.; Di Zazzo, E.; Giovannelli, P.; Galasso, G.; Castoria, G.; Migliaccio, A. ROS in cancer therapy: The bright side of the moon. Exp. Mol. Med. 2020, 5, 192–203. [Google Scholar] [CrossRef] [PubMed]
  49. Meshkini, A.; Yazdanparast, R. Involvement of oxidative stress in taxol-induced apoptosis in chronic myelogenous leukemia K562 cells. Exp. Toxicol. Pathol. 2012, 64, 357–365. [Google Scholar] [CrossRef] [PubMed]
  50. Asensio-López, M.C.; Soler, F.; Sánchez-Más, J.; Pascual-Figal, D.; Fernández-Belda, F.; Lax, A. Early oxidative damage induced by doxorubicin: Source of production, protection by GKT137831 and effect on Ca2+ transporters in HL-1 cardiomyocytes. Arch. Biochem. Biophys. 2016, 594, 26–36. [Google Scholar] [CrossRef] [PubMed]
Cancers 13 05509 g001 550
Figure 1. Post-hypoxic cells are less sensitive to doxorubicin and paclitaxel but not 5-FU or cisplatin treatment ex vivo. (a) Lentiviral vectors were delivered to MDA-MB-231 cells to generate a hypoxia fate-mapping system (left). A cartoon depicting that cells that experienced intratumoral hypoxia (GFP+) have a higher probability of forming lung metastases (right). (b) Schematic of the experimental set-up for (c–h). Tumors that formed from the orthotopic injection of MDA-MB-231 hypoxia fate-mapping cells in NSG mice were sorted into DsRed+ or GFP+ populations and then cultured in vitro to determine the IC50 of cisplatin (c), doxorubicin (d), 5-fluorouracil (e), or paclitaxel (f) after 48 h of treatment. IC50 curves were generated by assessing cell viability using the area of the cell culture well covered by nuclear DAPI staining (top) or using a Presto Blue assay (bottom) (N = 3, n = 3); GFP versus DsRed. RFU = relative fluorescence units. (g,h) Fluorescent images of tumor-derived DsRed+ and GFP+ sorted cells treated with doxorubicin (g) or paclitaxel (h) using increasing drug concentrations as indicated. Images are quantified in Figure S1f.
Figure 1. Post-hypoxic cells are less sensitive to doxorubicin and paclitaxel but not 5-FU or cisplatin treatment ex vivo. (a) Lentiviral vectors were delivered to MDA-MB-231 cells to generate a hypoxia fate-mapping system (left). A cartoon depicting that cells that experienced intratumoral hypoxia (GFP+) have a higher probability of forming lung metastases (right). (b) Schematic of the experimental set-up for (c–h). Tumors that formed from the orthotopic injection of MDA-MB-231 hypoxia fate-mapping cells in NSG mice were sorted into DsRed+ or GFP+ populations and then cultured in vitro to determine the IC50 of cisplatin (c), doxorubicin (d), 5-fluorouracil (e), or paclitaxel (f) after 48 h of treatment. IC50 curves were generated by assessing cell viability using the area of the cell culture well covered by nuclear DAPI staining (top) or using a Presto Blue assay (bottom) (N = 3, n = 3); GFP versus DsRed. RFU = relative fluorescence units. (g,h) Fluorescent images of tumor-derived DsRed+ and GFP+ sorted cells treated with doxorubicin (g) or paclitaxel (h) using increasing drug concentrations as indicated. Images are quantified in Figure S1f.
Cancers 13 05509 g001
Cancers 13 05509 g002 550
Figure 2. Post-hypoxic cells are less sensitive to doxorubicin and paclitaxel but not cisplatin or 5-FU treatment in vivo. (a) A schematic depicting the experimental set-up. Primary tumors were surgically removed 25 days post implantation following a single neoadjuvant dose of chemotherapy (Ctrl = control; Cis = cisplatin; Dox = doxorubicin; 5FU = 5-fluorouracil; Pax = paclitaxel). Three days after surgery, mice were treated with daily doses of chemotherapy (dose #4 of paclitaxel was suspended). Mice were sacrificed to assess metastatic burden 3 days following the last dose. (b) Whole mount fluorescent image of lung (left) and liver (right) from one representative mouse per treatment group. (c) Percentage of fluorescent (DsRed+ and GFP+) cells measured by flow cytometry analysis (orange arrows indicate decreased metastatic burden compared to the control group). The ratio of mice with detectable metastasis over the total number in the treatment group is displayed for each condition. (d,e) The percentage of DsRed+ or GFP+ cells in the lung (d) and liver (e) as measured by flow cytometry (n = 4–5); p-value is displayed to compare GFP versus DsRed (one-tailed paired t-test). Insets show a shorter y-axis scale for doxorubicin and paclitaxel treatments. (f) Gene set enrichment analysis (GSEA) of gene sets associated with resistance to cisplatin, doxorubicin, 5-FU, and paclitaxel in GFP+ (LG) versus DsRed+ (LR) metastatic cancer cells isolated from the lung. Normalized enrichment score (NES) and p-value (p-val) are displayed.
Figure 2. Post-hypoxic cells are less sensitive to doxorubicin and paclitaxel but not cisplatin or 5-FU treatment in vivo. (a) A schematic depicting the experimental set-up. Primary tumors were surgically removed 25 days post implantation following a single neoadjuvant dose of chemotherapy (Ctrl = control; Cis = cisplatin; Dox = doxorubicin; 5FU = 5-fluorouracil; Pax = paclitaxel). Three days after surgery, mice were treated with daily doses of chemotherapy (dose #4 of paclitaxel was suspended). Mice were sacrificed to assess metastatic burden 3 days following the last dose. (b) Whole mount fluorescent image of lung (left) and liver (right) from one representative mouse per treatment group. (c) Percentage of fluorescent (DsRed+ and GFP+) cells measured by flow cytometry analysis (orange arrows indicate decreased metastatic burden compared to the control group). The ratio of mice with detectable metastasis over the total number in the treatment group is displayed for each condition. (d,e) The percentage of DsRed+ or GFP+ cells in the lung (d) and liver (e) as measured by flow cytometry (n = 4–5); p-value is displayed to compare GFP versus DsRed (one-tailed paired t-test). Insets show a shorter y-axis scale for doxorubicin and paclitaxel treatments. (f) Gene set enrichment analysis (GSEA) of gene sets associated with resistance to cisplatin, doxorubicin, 5-FU, and paclitaxel in GFP+ (LG) versus DsRed+ (LR) metastatic cancer cells isolated from the lung. Normalized enrichment score (NES) and p-value (p-val) are displayed.
Cancers 13 05509 g002
Cancers 13 05509 g003 550
Figure 3. Post-hypoxic cells are less sensitive to doxorubicin and paclitaxel. (a) Schematic depicting the experimental set-up: primary tumors were surgically resected 25 days post implantation following one IV dose of neoadjuvant chemotherapy. Three days after surgery, the mice were treated with four doses of chemotherapy; mice were sacrificed 3 days following treatment. (b) Whole mount fluorescent image of lung (left) and liver (right) from one representative mouse per treatment group. (c) Percentage of fluorescent (DsRed+ and GFP+) cells measured by flow cytometry analysis (orange arrows indicate the decrease in metastatic burden compared to the control group). The ratio of mice with detectable metastasis over the total number in the treatment group is displayed for each condition. (d,e) The percentage of DsRed+ or GFP+ cells in the lung (d) and liver (e) as measured by flow cytometry (n = 8–9); p-value is displayed to compare GFP versus DsRed (one-tailed paired t-test). Insets show a reduced y-axis scale for doxorubicin and paclitaxel treatments. (f) Venn diagram displaying the 19 genes upregulated in both GFP+ (TG) versus DsRed+ (TR) tumor cells and GFP+ (LG) versus DsRed+ (LR) metastatic cells in the lung termed ‘hypoxic memory’ previously reported [25]. (g,h) Microarray expression data were used to perform Kaplan–Meier analysis of distant metastasis-free survival (DMFS) of breast cancer patients stratified by high or low expression of the ’hypoxic memory’ 19-gene signature and treated with either (g) chemotherapy (N = 225) or (h) endocrine therapy (N = 205).
Figure 3. Post-hypoxic cells are less sensitive to doxorubicin and paclitaxel. (a) Schematic depicting the experimental set-up: primary tumors were surgically resected 25 days post implantation following one IV dose of neoadjuvant chemotherapy. Three days after surgery, the mice were treated with four doses of chemotherapy; mice were sacrificed 3 days following treatment. (b) Whole mount fluorescent image of lung (left) and liver (right) from one representative mouse per treatment group. (c) Percentage of fluorescent (DsRed+ and GFP+) cells measured by flow cytometry analysis (orange arrows indicate the decrease in metastatic burden compared to the control group). The ratio of mice with detectable metastasis over the total number in the treatment group is displayed for each condition. (d,e) The percentage of DsRed+ or GFP+ cells in the lung (d) and liver (e) as measured by flow cytometry (n = 8–9); p-value is displayed to compare GFP versus DsRed (one-tailed paired t-test). Insets show a reduced y-axis scale for doxorubicin and paclitaxel treatments. (f) Venn diagram displaying the 19 genes upregulated in both GFP+ (TG) versus DsRed+ (TR) tumor cells and GFP+ (LG) versus DsRed+ (LR) metastatic cells in the lung termed ‘hypoxic memory’ previously reported [25]. (g,h) Microarray expression data were used to perform Kaplan–Meier analysis of distant metastasis-free survival (DMFS) of breast cancer patients stratified by high or low expression of the ’hypoxic memory’ 19-gene signature and treated with either (g) chemotherapy (N = 225) or (h) endocrine therapy (N = 205).
Cancers 13 05509 g003
Cancers 13 05509 g004 550
Figure 4. Post-hypoxic cells are more likely to contribute to metastatic recurrence after treatment. (a) Schematic of the experimental approach. Tumor-derived MDA-MB-231 cells were sorted into DsRed+ or GFP+ populations and treated with 10 nM of paclitaxel for 48 h. Fresh media were added and the cells were cultured for an additional 30 days to determine how many colonies would form. (b) The number of colonies that formed by day 20 of the recovery period was counted (n = 10); p-value is displayed to compare GFP versus DsRed (two-tailed unpaired t-test). (c) Schematic depicting the experimental timeline: primary tumors were surgically resected 25 days post implantation following one IV dose of neoadjuvant chemotherapy. Three days after surgery, the mice were treated with four doses of chemotherapy. Control mice were sacrificed to assess metastatic burden 3 days after treatment, whereas mice treated with paclitaxel were sacrificed 7 days later. (d) Whole mount fluorescent image of lung (top) and liver (bottom) of one representative mouse per treatment group. (e,f) Percentage of fluorescent (DsRed+ and GFP+) cells in (e) lung and (f) liver measured by flow cytometry analysis. The ratio of mice with detectable metastasis over the total number in the treatment group is displayed for each condition. The percentage of metastatic burden in the paclitaxel-treated group (Pax) that was sacrificed 3 days after the last treatment (day 35) presented in Figure 3 is shown as a comparison (N = 2, n = 13–28); treatment versus control (1-way ANOVA with Dunnett (recommended) multiple comparison test). (g,h) The percentage of DsRed+ and GFP+ cells in the lung (g) and liver (h) was quantified by flow cytometry (N = 2, n = 13–28); p-value is displayed for GFP versus DsRed cells (ordinary 2-way ANOVA with Sidak multiple comparison post-test).
Figure 4. Post-hypoxic cells are more likely to contribute to metastatic recurrence after treatment. (a) Schematic of the experimental approach. Tumor-derived MDA-MB-231 cells were sorted into DsRed+ or GFP+ populations and treated with 10 nM of paclitaxel for 48 h. Fresh media were added and the cells were cultured for an additional 30 days to determine how many colonies would form. (b) The number of colonies that formed by day 20 of the recovery period was counted (n = 10); p-value is displayed to compare GFP versus DsRed (two-tailed unpaired t-test). (c) Schematic depicting the experimental timeline: primary tumors were surgically resected 25 days post implantation following one IV dose of neoadjuvant chemotherapy. Three days after surgery, the mice were treated with four doses of chemotherapy. Control mice were sacrificed to assess metastatic burden 3 days after treatment, whereas mice treated with paclitaxel were sacrificed 7 days later. (d) Whole mount fluorescent image of lung (top) and liver (bottom) of one representative mouse per treatment group. (e,f) Percentage of fluorescent (DsRed+ and GFP+) cells in (e) lung and (f) liver measured by flow cytometry analysis. The ratio of mice with detectable metastasis over the total number in the treatment group is displayed for each condition. The percentage of metastatic burden in the paclitaxel-treated group (Pax) that was sacrificed 3 days after the last treatment (day 35) presented in Figure 3 is shown as a comparison (N = 2, n = 13–28); treatment versus control (1-way ANOVA with Dunnett (recommended) multiple comparison test). (g,h) The percentage of DsRed+ and GFP+ cells in the lung (g) and liver (h) was quantified by flow cytometry (N = 2, n = 13–28); p-value is displayed for GFP versus DsRed cells (ordinary 2-way ANOVA with Sidak multiple comparison post-test).
Cancers 13 05509 g004
Cancers 13 05509 g005 550
Figure 5. Post-hypoxic cells develop a breast cancer stem cell-like phenotype. (a) Tumor-derived MDA-MB-231 cells were immediately transferred to mammosphere formation media for a 7-day period and imaged by fluorescent microscopy. (b,c) The number and size of mammospheres that formed were assessed by fluorescent image analysis and normalized by the initial number of DsRed+ and GFP+ cells (N = 3, n = 6); p-value is displayed to compare GFP versus DsRed (two-tailed paired t-test). (d) A GSEA of gene set associated with stemness (Zhang et al., 2008) in DsRed+ (LR) and GFP+ (LG) cancer cells purified from lungs. Normalized enrichment score (NES) and p-value (p-val) are displayed. (e) Heatmap of the expression of the top 10 genes driving the enrichment of stemness in (d) in DsRed+ (LR1 and LR2) and GFP+ (LG1 and LG2) lung-derived cancer cells. (f–i) ROC analysis of the top 10 stemness-associated genes for patients that did or did not have a response to treatment with (f) taxane alone; (g) endocrine therapy; (h) anthracycline alone; or (i) fluorouracil, Adriamycin (doxorubicin), and Cytoxan (FAC) therapies.
Figure 5. Post-hypoxic cells develop a breast cancer stem cell-like phenotype. (a) Tumor-derived MDA-MB-231 cells were immediately transferred to mammosphere formation media for a 7-day period and imaged by fluorescent microscopy. (b,c) The number and size of mammospheres that formed were assessed by fluorescent image analysis and normalized by the initial number of DsRed+ and GFP+ cells (N = 3, n = 6); p-value is displayed to compare GFP versus DsRed (two-tailed paired t-test). (d) A GSEA of gene set associated with stemness (Zhang et al., 2008) in DsRed+ (LR) and GFP+ (LG) cancer cells purified from lungs. Normalized enrichment score (NES) and p-value (p-val) are displayed. (e) Heatmap of the expression of the top 10 genes driving the enrichment of stemness in (d) in DsRed+ (LR1 and LR2) and GFP+ (LG1 and LG2) lung-derived cancer cells. (f–i) ROC analysis of the top 10 stemness-associated genes for patients that did or did not have a response to treatment with (f) taxane alone; (g) endocrine therapy; (h) anthracycline alone; or (i) fluorouracil, Adriamycin (doxorubicin), and Cytoxan (FAC) therapies.
Cancers 13 05509 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Godet, I.; Mamo, M.; Thurnheer, A.; Rosen, D.M.; Gilkes, D.M. Post-Hypoxic Cells Promote Metastatic Recurrence after Chemotherapy Treatment in TNBC. Cancers 2021, 13, 5509. https://doi.org/10.3390/cancers13215509

AMA Style

Godet I, Mamo M, Thurnheer A, Rosen DM, Gilkes DM. Post-Hypoxic Cells Promote Metastatic Recurrence after Chemotherapy Treatment in TNBC. Cancers. 2021; 13(21):5509. https://doi.org/10.3390/cancers13215509

Chicago/Turabian Style

Godet, Inês, Mahelet Mamo, Andrea Thurnheer, D. Marc Rosen, and Daniele M. Gilkes. 2021. "Post-Hypoxic Cells Promote Metastatic Recurrence after Chemotherapy Treatment in TNBC" Cancers 13, no. 21: 5509. https://doi.org/10.3390/cancers13215509

APA Style

Godet, I., Mamo, M., Thurnheer, A., Rosen, D. M., & Gilkes, D. M. (2021). Post-Hypoxic Cells Promote Metastatic Recurrence after Chemotherapy Treatment in TNBC. Cancers, 13(21), 5509. https://doi.org/10.3390/cancers13215509

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop