Skip to main content
  • Research article
  • Open access
  • Published:

BHLHE40 confers a pro-survival and pro-metastatic phenotype to breast cancer cells by modulating HBEGF secretion

Abstract

Background

Metastasis is responsible for a significant number of breast cancer-related deaths. Hypoxia, a primary driving force of cancer metastasis, induces the expression of BHLHE40, a transcription regulator. This study aimed to elucidate the function of BHLHE40 in the metastatic process of breast cancer cells.

Methods

To define the role of BHLHE40 in breast cancer, BHLHE40 expression was knocked down by a lentiviral construct expressing a short hairpin RNA against BHLHE40 or knocked out by the CRISPR/Cas9 editing system. Orthotopic xenograft and experimental metastasis (tail vein injection) mouse models were used to analyze the role of BHLHE40 in lung metastasis of breast cancer. Global gene expression analysis and public database mining were performed to identify signaling pathways regulated by BHLHE40 in breast cancer. The action mechanism of BHLHE40 was examined by chromatin immunoprecipitation (ChIP), co-immunoprecipitation (CoIP), exosome analysis, and cell-based assays for metastatic potential.

Results

BHLHE40 knockdown significantly reduced primary tumor growth and lung metastasis in orthotopic xenograft and experimental metastasis models of breast cancer. Gene expression analysis implicated a role of BHLHE40 in transcriptional activation of heparin-binding epidermal growth factor (HBEGF). ChIP and CoIP assays revealed that BHLHE40 induces HBEGF transcription by blocking DNA binding of histone deacetylases (HDAC)1 and HDAC2. Cell-based assays showed that HBEGF is secreted through exosomes and acts to promote cell survival and migration. Public databases provided evidence linking high expression of BHLHE40 and HBEGF to poor prognosis of triple-negative breast cancer.

Conclusion

This study reveals a novel role of BHLHE40 in promoting tumor cell survival and migration by regulating HBEGF secretion.

Background

One in every eight women in the USA will be diagnosed with breast cancer over the course of her lifetime [1]. An estimated 266,120 new cases are expected, and 40,920 women are expected to die from breast cancer in 2018 in the USA [1]. Distant metastasis is the major cause of breast cancer-related deaths. Hypoxia has been recognized as a primary driving force of distant metastasis of breast cancer [2,3,4,5,6,7]. Among hypoxia-responsive genes are both promoting and suppressive factors for malignant progression. It is unclear how the expression and activities of metastasis-promoting factors are preferentially augmented in metastatic tumors. A large body of studies have focused on elucidating the molecular mechanisms by which hypoxia enables cancer cells to survive a variety of stresses imposed by the metastatic process, including nutrient depletion, loss of attachment, and deprivation of growth factors.

Hypoxia-induced exosomic secretion of cytokines and growth factors plays a key role in promoting metastasis through both tumor autonomous and non-autonomous mechanisms [8, 9]. Exosomes are microvesicles (40–130 nm) constitutively released by a variety of cells into the extracellular environment to promote cell-to-cell communication [10]. Tumor cells have been reported to utilize exosomes to transfer nucleotides, lipids, and proteins into surrounding cells or cells in distant metastatic niches [10, 11]. Hypoxia is known to markedly increase the number of secreted exosomes, as well as alter the contents of exosomes [12]. However, our understanding of the regulation of exosome secretion is rudimentary.

Although the cellular response to hypoxia is mainly controlled by two basic helix-loop-helix transcription factors, hypoxia inducible factor (HIF)1A and EPAS1/HIF2A, the outcomes of the hypoxia response are modified by other transcription regulators that are regulated by hypoxia or interact with HIF1A or EPAS1 [13, 14]. This study focuses on the role of a basic helix-loop-helix transcription factor BHLHE40 (also known as DEC1/BHLHB2/SHARP2/STRA13) in metastasis of breast cancer. BHLHE40 expression is directly activated by HIF1A in a variety of tumor cells under hypoxia [15, 16]. High BHLHE40 expression has been linked to activation of a hypoxia-response pathway, elevated metastatic potentials, and poor prognosis of various types of tumors, including hepatocellular carcinoma, pancreatic cancer, and invasive breast cancer [17,18,19]. It was reported that BHLHE40 binds to the E-box elements and regulates the expression of genes associated with circadian rhythm, cell differentiation, cell senescence, lipid metabolism, DNA damage response, and immune response [20,21,22,23]. However, the mechanism of action and downstream targets of BHLHE40 in breast cancer cells is largely unknown. In this study, we provide evidence which suggests that BHLHE40 is a pro-metastasis factor in breast cancer cells which promotes tumor cell survival and migration by modulating exosomic secretion of heparin-binding epidermal growth factor (HBEGF).

Methods

Cell culture

Breast cancer MDA-MB-231 and MCF7 cells were obtained from ATCC (Manassas, VA, USA) and maintained in minimal essential medium (ThermoFisher Scientific, Rockford, IL, USA) supplemented with 10% fetal bovine serum (FBS), 200 U/ml penicillin-streptomycin, and 0.5 μg/ml amphotericin B (Cellgro, Manassa, VA, USA). A lung metastatic derivative of MDA-MB-231 (LM) and a tamoxifen-resistant derivative of MCF7 (TR) were established as described previously [24, 25]. A stable line (BHLHE40-KD) expressing a short hairpin RNA (shRNA) against BHLHE40 (TRCN0000232187, Sigma-Aldrich, St. Louis, MO, USA) was generated by lentiviral transduction and selected in medium containing 2 μg/ml puromycin (Sigma-Aldrich). A colony of BHLHE40 knockout variant (BHLHE40-KO) of MDA-MB-231 was generated by using the CRISPR/Cas9 all-in-one expression system (HCP221270-CG01–1, GeneCopoeia) and selected in medium supplemented with 500 μg/ml gentamicin. Elimination of BHLHE40 expression was examined by immunoblotting and quantitative polymerase chain reaction (qPCR) with a forward primer designed to cover the CRISPR editing site (forward 5’GACGGGGAATAAAGCGGAGC and reverse 5’CCGGTCACGTCTCTTTTTCTC). To knockout HIF1A and EPAS1 by CRISPR/Cas9 editing, gRNA targeting exon 1 of HIF1A or EPAS1 was individually cloned into the pX462-puromycin and pX462-hygromycin vectors (expressing Cas9n, AddGene), respectively. MDA-MB-231 cells were transfected with pX462-puro-HIF1A gRNAs using FuGene HD followed by selection in puromycin. A clonal line that had no HIF1A protein detected by immunoblotting was transfected with pX462-hygromycin-EPAS1 gRNA and selected by hygromycin. Hygromycin-resistant colonies that had no EPAS1 detected by immunoblotting were pooled to generate a HIF1A/EPAS1 double-knockout (HIF-dKO) subline. For all knockdown (KD) or KO sublines, control cells were transfected with corresponding empty vectors (EV) and selected in antibiotics in parallel with cells transfected with shRNAs or gRNAs. Pooled drug-resistant colonies of control cells were used as EV control lines. Cells were exposed to different conditions relevant to solid tumors, including loss of attachment (suspension culture), hypoxia (1% O2), hypoxia in combination with low (1 mg/ml) glucose (1%O2/LG), and hypoxia in combination with glucose depletion (1%O2/GF, a condition that induces rapid apoptosis of tumor cells).

Orthotropic xenograft and experimental lung metastasis models

All in-vivo studies were performed in accordance with the protocols approved by the Institutional Animal Care and Use Committee (IACUC) of the University of Tennessee Health Science Center. NOD.Cg-PrkdcscidIl2rgtm1Wjl/SzJ (NSG) mice were purchased from Jackson Laboratories (Bar Harbor, ME, USA). Orthotropic xenograft and experimental lung metastasis (tail vein injection) models were established using fluorescence-labeled tumor cells, as previously described [24, 26]. Tumor size was monitored and measured weekly using digital calipers. Tumor volume was calculated as: volume = (width2 × length)/2. Lung metastasis was quantified by fluorescent imaging of lungs and qPCR of human Alu DNA repeats (forward primer: 5′: GTCAGGAGATCGAGACCATCCC 3′; reverse primer: 5′: TCCTGCCTCAGCCTCCCAAG 3′). Circulating tumor cells in whole blood collected by cardiac puncture were isolated using the Ficoll-Paque PLUS medium (GE Healthcare Life Sciences, Piscataway, NJ, USA) and counted under fluorescent microscope.

Migration, invasion, and wound healing assays

Transwell membrane inserts with 8-μm pores (BD Biosciences, Bedford, MA, USA), uncoated or coated with Matrigel, were used to determine the migratory and invasive activities of cancer cells, respectively. Cells undergoing migration and invasion were expressed as: percent migration = mean number of cells migrating through the uncoated transwell × 100/mean number of seeded cells; percent invasion = mean number of cells migrating through the Matrigel-coated transwell × 100/mean number of migrating cells through the uncoated pores. Real-time assessment of migratory activity during scratch wound healing was performed using the IncuCyte ZOOM-ImageLock plate system (Essen Bioscience, Michigan, USA). To examine the expression levels of proteins in response to scratch wounds, the EMD Millipore Chemicon Cell Comb Scratch assay kit (Millipore) was used to generate a high-density field of scratches in a confluent cell monolayer to maximize the area of wound edges. To examine the effect of HBEGF on cell migration and invasion, a neutralizing antibody to HBEGF (10 μg/ml; AF-259-NA, R&D systems, Minneapolis, MN, USA) or a HBEGF peptide (20 μg/ml) was added to the medium.

Suspension culture, viable cell counting, and caspase assays

To mimic the loss of attachment, cells were cultured in PolyHEMA (Sigma-Aldrich)-coated plates to prevent adherence. Methylcellulose (1%) was added to the medium to prevent formation of large cell aggregates to accurately measure tumor cell proliferation in the suspension. Viable cells were counted by the Trypan blue exclusion method. Cell apoptosis was determined using the caspase Glo 3/7 assay kit (Promega, Madison, WI, USA) or immunoblotting of cleaved Caspase 9.

Luciferase reporter assay for HIF activity

Luciferase reporter constructs driven by hypoxia responsive elements of LDHA (LDHA-Luc, S721613) and ITGA6 (ITGA6-Luc, S708174) were purchased from SwitchGear Genomics (Carlsbad, CA, USA). Cells (3 × 105) were transfected with 500 ng LDHA-Luc or ITGA6-Luc and 10 ng CMV-β-galactosidase (a control for transfection efficiency), along with 250 ng of an empty vector or a HIF1A-expressing construct (HsCD00444875, DNASU plasmid repository), using Fugene6 (Promega). Forty-eight hours after transfection, the cells were exposed to hypoxia (1% O2) for 6 h. Luciferase and β-galactosidase (β-gal) activities were measured using the LightSwitch Luciferase Assay System and Promega Beta-Glo Assay System, respectively. Relative luciferase activities normalized to β-gal were presented as mean ± SD, n = 6.

Exosome isolation and analysis

To isolate exosomes, 5 × 106 cells were cultured in medium supplemented with exosome-free serum (SystemBio, Palo Alto, CA, USA). Exosomes in the conditioned medium were purified using the ExoQuick-TC solution (SystemBio) and quantified under a fluorescent microscope after being labeled with carboxyfluorescein diacetate succinimidyl ester (CFSE) using the Exo-Glow labeling kit (SystemBio) that is designed to exclude background particles. To analyze the protein content, isolated exosomes were lysed in RIPA buffer (Pierce, ThermoScientific) containing protease inhibitor cocktail (Sigma-Aldrich) and subjected to immunoblot analysis. To examine the effect of purified exosomes on cell migration, exosomes were re-suspended in exosome-free medium and added to cells seeded in transwells.

Gene expression microarray and qPCR analysis

Total RNA from cells exposed to hypoxia (1% O2 for 6 and 48 h) was purified using the RNeasy kit (Qiagen) and submitted to the Molecular Resource Center at the University of Tennessee Health Science Center for labeling and hybridization to the HT-12 expression BeadChips (Illumina, Chicago, IL, USA). Hybridization signals were processed for annotation, background subtraction, quantile normalization, and presence call filtering using the Gene Expression Module of the Genome Studio Software (Illumina). The microarray data can be found in the Gene Expression Omnibus database with accession number GSE107300. Hypoxia-responsive genes were defined as genes whose expression was altered ≥ 1.5-fold (hypoxia versus control) in two independent experiments. BHLHE40 target genes were defined as genes whose expression was altered ≥ 1.5-fold (BHLHE40-KD versus EV) in two independent experiments. To examine the effect of BHLHE40-KD on gene expression in cells exposed to hypoxia (1% O2) in combination with low glucose (1 mM; 1%O2/LG), a condition frequently encountered by cells in solid tumors, total RNA of cells exposed to 1%O2/LG for 4 h from three independent experiments were pooled and analyzed using the GeneChip Human Gene 1.0 ST array (Affymetrix, Santa Clara, CA, USA). The Affymetrix data were extracted, normalized, and summarized with the robust multi-average (RMA) method implemented in the Affymetrix Expression Console. To validate microarray data by qPCR, total RNA from cells in three independent experiments with duplicates was prepared using trizol (Life Technologies, Grand Island, NY, USA). cDNAs were synthesized using iScript cDNA Synthesis kits (Bio-Rad, Hercules, CA, USA) and qPCR was performed on the CFX96TM Real-Time PCR detection system using SYBR green supermix (Bio-Rad). Expression data of mRNA were normalized by the 2–ΔΔCT method to RPL13A and presented as mean ± SD. Primers for qPCR were obtained from the Primerbank [27].

Protein extraction, co-immunoprecipitation, and immunoblotting

For immunoblotting (IB) analysis, whole cell lysates and nuclear proteins were prepared using RIPA buffer supplemented with protease inhibitor cocktails (Sigma-Aldrich) and the NE-PER Nuclear and Cytoplasmic Extraction Reagents (Thermo Scientific), respectively. To detect protein-protein interaction, soluble proteins were extracted using the Pierce IP Lysis Buffer (Thermo Scientific) supplemented with protease inhibitor cocktails and co-immunoprecipitation (CoIP) was performed using the TrueBlot Immunoprecipitation and Western Blot Kit (Rockland Immunochemicals Inc., Limerick, PA, USA). IB signals were developed using the SuperSignal West Dura Extended Duration Substrate and the CL-XPosure Film (Thermo Scientific). Antibodies used in this study were: anti-EGFR-Tyr1110P, anti-EGFR, anti-ERK1/2-Thr202/Tyr204P, anti-ERK1/2, anti-AKT-Tyr416P, anti-AKT, anti-Caspase9, anti-HDAC1, and anti-HDAC2 from Cell Signaling Technologies (Boston, MA, USA), anti-GAPDH from Millipore (Merck, Darmstadt, Germany), anti-TBP from Abcam (Cambridge, MA, USA), anti-CTGF from Abgent (San Diego, CA, USA), anti-HBEGF and anti-CD9 from R&D Biosystems (Minneapolis, MN, USA), anti-BHLHE40, anti-HIF1A, and anti-EPAS1 from Novus Biologicals (Littleton, CO, USA), and anti-ALIX, anti-TSG101, and anti-CD81 from Santa Cruz (Dallas, TX, USA).

Chromatin immunoprecipitation (ChIP)

Protein-DNA crosslink, nuclear fraction extraction, and chromatin fragmentation were performed as described previously [26]. The soluble fraction of sheared chromatin (200–500 bp in length) was pre-cleaned with Magnabind Goat anti-Rabbit IgG (for anti-BHLHE40 and anti-HDAC1) or Magnabind Goat anti-Mouse (for anti-HDAC2) (Life Technologies), followed by immunoprecipitation with control IgG or antibodies against BHLHE40, histone deacetylases (HDAC)1 or HDAC2, and Magnabind beads conjugated with a secondary antibody. DNA in the de-crosslinked immunocomplexes was isolated with the MiniElute PCR purification kit (Qiagen, Germantown, MD, USA). qPCR was performed to detect the presence of the proximal promoter region of HBEGF (−529 to −372 from the transcription start site) using the following primers: forward 5’TGCCTGCAACTTCAACT CTG3′ and reverse 5’CCATCCCTGTCACCCTCTAA3′.

Statistical analysis

Student’s t tests, one-way analysis of variance (ANOVA) with post-hoc Tukey test and correlation significance analyses were performed using the GraphPad Prism 5 software (GraphPad, San Diego, CA, USA); p values < 0.05 were considered statistically significant.

Results

BHLHE40 knockdown leads to decreased primary tumor growth and lung metastases

To define the role of BHLHE40 in breast cancer metastasis, we examined the effect of its knockdown (KD) by a shRNA lentiviral construct on spontaneous lung metastasis of orthotopic xenograft tumors derived from a lung metastasis-enriched subline (LM) of breast cancer MDA-MB-231 cells [28]. The protein levels of BHLHE40 is low in cells under normal growth conditions but is significantly induced by hypoxia (1% O2, 16 h). BHLHE40-shRNA expression effectively reduced both baseline and hypoxia-induced levels of BHLHE40 in LM cells (Fig. 1a). In NSG mice inoculated with 2 × 105 control LM-EV (empty vector) cells in the inguinal mammary gland fat pads, palpable tumors were detected at 2 weeks (Fig. 1b) and lung metastasis became evident at 5 weeks (Fig. 1c) post-inoculation. BHLHE40-KD delayed the onset of primary tumors, which became palpable 3 weeks after inoculation, and reduced the growth rate of primary tumors, coincident with decreased lung metastases (Fig. 1a–c). To further investigate the effect of BHLHE40-KD on lung metastases, primary tumors of EV and BHLHE40-KD cells were surgically removed at 3 and 5 weeks post-inoculation, respectively, when they reached similar size with a diameter of 4–5 mm. Lung metastasis was examined 4 weeks after primary tumor resection (Fig. 1d). BHLHE40-KD substantially reduced lung metastasis in mice with similar primary tumor burdens. Taken together, these results suggest that BHLHE40 plays a role in promoting primary tumor growth and spontaneous distant metastasis of breast cancer cells.

Fig. 1
figure 1

BHLHE40-knockdown (KD) significantly reduced primary tumor size and lung metastatic burden in an orthotopic xenograft model. a BHLHE40-shRNA expression effectively reduced both baseline and hypoxia-induced expression of BHLHE40 protein in the LM cells, as determined by immunoblotting. b Orthotopic xenograft tumors derived from LM-BHLHE40-KD cells exhibited lower growth rate than tumors derived from control LM empty vector (EV) cells. NSG mice were inoculated in the inguinal mammary gland fat pads with 2 × 105 cells. Tumor size was monitored and measured weekly using a digital caliper. Tumor volume was calculated as: volume = (width2 × length)/2. *p <  0.05 (n = 20, KD vs. EV at indicated time points), one-way ANOVA followed by Tukey’s post-hoc tests. c Spontaneous lung metastasis detected by fluorescent imaging of lungs or human ALU repeats qPCR 5 weeks after inoculation of tumor cells in mammary gland fat pads. *p <  0.05 (n = 10, KD vs. EV), Student’s t test. d Lung metastasis in mice after resection of primary tumors. Primary tumors in mammary gland fat pads were resected when they reached a size of 5 × 5 mm and lung metastasis were analyzed 4 weeks post-resection by fluorescent imaging of lungs or human ALU repeats qPCR. *p <  0.05 (n = 10, KD vs. EV), Student’s t test

BHLHE40 knockdown reduces lung colonization of tumor cells inoculated through tail vein

To determine whether BHLHE40 regulates late metastatic events after entry of tumor cells into the blood stream, we examined the effect of BHLHE40-KD on the ability of tumor cells to survive circulation and colonize in the lungs using an experimental metastasis model, in which tumor cells were delivered into the blood stream through tail vain injection to bypass the initial steps of metastasis such as migration and intravasation. LM-EV and LM-BHLHE40-KD cells (5 × 105) were injected into the left lateral tail veins of 5-week-old female NSG mice, and tumor cells in the bloodstream and lung tissues were examined at various times post-injection (Fig. 2). Compared with control LM-EV cells, LM-BHLHE40-KD cells were more rapidly eliminated from the bloodstream (Fig. 2a). LM-EV cells were observed in lung tissues at 72 h and formed large metastatic foci at 4 weeks after tail vein injection (Fig. 2b, c). In contrast, BHLHE40-KD cells were not detected in lung tissue at 72 h and formed less metastatic foci in lungs than EV cells at various time points (Fig. 2b, c). No fluorescent loci of EV or BHLHE40-KD cells were found in other organs (i.e., livers, spleens, and kidneys) within 5 weeks after tail vein inoculation. Together, these results suggest that BHLHE40 is required for tumor cells to survive in the circulation and establish metastatic foci in the lungs.

Fig. 2
figure 2

BHLHE40-knockdown (KD) reduced lung colonization of tumor cells inoculated into the circulation via tail veins. a Less tumor cells were detected in blood of NSG mice inoculated via tail vein injection with LM BHLHE40-KD than mice inoculated with empty vector (EV) control LM EV cells. Tumor cells in whole blood collected by cardiac puncture at indicated times after tail vein injection were isolated using the Ficoll-Paque PLUS medium and counted under fluorescent microscope. *p <  0.05 (n = 3, KD vs. EV), one-way ANOVA followed by Tukey’s post-hoc tests. b Fluorescent imaging of metastatic foci in lungs at different time points post tail-vein injection. c Percentage of areas occupied by metastatic loci in the lungs, as quantified by fluorescent imaging and ImageJ. *p <  0.05 (n = 3, KD vs. EV), one-way ANOVA followed by Tukey’s post-hoc tests

BHLHE40 acts to promote cell migration, invasion and survival

Having established a role for BHLHE40 in distant metastasis of breast cancer cells in vivo, we sought to identify the specific cellular processes that require BHLHE40 activity. Despite the significant effect of BHLHE40-KD on primary tumor growth and lung metastasis of LM cells in vivo, BHLHE40-KD showed no significant effect on proliferation of LM cells under normal two-dimensional growth conditions in vitro. The doubling times, determined by Trypan blue exclusion-based cell counting (daily for 7 days), of LM-EV and LM-BHLHE40-KD cells were 36.37 ± 0.49 h (n = 6) and 38.95 ± 3.61 h (n = 6), respectively. Therefore, we focused on investigating whether BHLHE40 is a downstream effector of HIF1A activation by hypoxia or loss of attachment. Detached breast cells were reported to rely on HIF1A activation to survive under normoxia [29]. In-vitro cell migration and invasion assays showed that BHLHE40-KD reduced the ability of cells to penetrate either uncoated or Matrigel-coated transwells under hypoxia conditions (1% O2; Fig. 3a). Under the nonadherent culture condition for 15 days, in which cells were mixed with growth medium supplemented with 1% methylcellulose to prevent cell aggregation and then seeded in plates coated with polyHEMA to prevent adherence, the number of viable LM BHLHE40-KD cells was significantly lower than LM EV cells (Fig. 3a, lower panel). To examine whether the effects exerted by BHLHE40-KD on LM cells can be extended to the parent MDA-MB-231 cells, we established a BHLHE40 knockout (KO) subline using the CRISPR/Cas9 editing system. BHLHE40 protein depletion in the KO subline under normoxia or hypoxia was confirmed by immunoblotting (Fig. 3b, upper panel). Although residue BHLHE40 protein was detected in the KO subline by immunoblotting, no wild-type mRNA was detected by qPCR with a forward primer designed to cover the CRISPR editing site. BHLHE40-KO resulted in a reduced number of viable cells after a 15-day suspension culture in plates coated with polyHEMA to prevent attachment (Fig. 3b, middle panel), while it showed no significant effect on the number of viable cells after a 5-day adherence culture. In addition, BHLHE40-KO significantly sensitized MDA-MB-231 cells to apoptosis induced by hypoxia in combination with glucose depletion (1%O2/GF, 6 h), as evidenced by the appearance of apoptotic morphology and activation of caspase3/7 (Fig. 3b, lower panel).

Fig. 3
figure 3

BHLHE40 depletion impaired cell migration, invasion, and survival. a BHLHE40-knockdown (KD) by shRNA in LM cells reduced cell migration and invasion, as well as the number of viable cells during suspension culture in comparison to empty vector (EV) control LM EV cells. Migratory and invasive activities were determined using transwells, uncoated or coated with Matrigel, respectively. The results were presented as: percent migration = mean number of cells migrating through the uncoated transwells × 100/mean number of seeded cells; percent invasion = mean number of cells migrating through the Matrigel-coated transwells × 100/mean number of migrating cells through the uncoated transwells. Suspension culture was conducted by seeding cells in medium containing 1% methylcellulose and in dishes coated with PolyHEMA for 15 days. Viable cells were counted under fluorescent microscope. *p <  0.05 (n = 6, KD vs. EV), Student’s t test. b BHLHE40-knockout (KO) by CRISPR/Cas9 editing in MDA-MB-231 cells reduced the number of viable cells after suspension culture and enhanced apoptosis induced by hypoxia combined with glucose depletion (1%O2/GF). Viable cells were determined by Trypan blue exclusion-based cell counting after a 15-day suspension culture. Apoptosis of cells exposed to 1%O2/GF for 6 h was examined by the appearance of apoptotic morphology (as indicated by arrows in the cell images) and Caspase3/7 assays. *p <  0.05 (n = 6, KO vs. EV), one-way ANOVA followed by Tukey’s post-hoc tests. c BHLHE40-KD by shRNA in tamoxifen-resistant subline of MCF7 (TR) reduced the ability of cells to survive 1%O2/GF (6 h) and reduced the number of viable cells after a 15-day suspension culture. Elevated expression of HIF1A and BHLHE40 in TR and fulvestrant-resistant (FR) sublines of MCF7 cells were detected by qPCR. BHLHE40-KD was confirmed by immunoblotting. Apoptosis induced by 1%O2/GF (6 h) was determined by Caspase3/7 assays. The number of viable cells after a 15-day suspension culture was determined by Trypan blue exclusion-based cell counting. **p <  0.05 (n = 6, TR or FR vs. parent MCF7), *p <  0.05 (n = 6, KD vs. EV), one-way ANOVA followed by Tukey’s post-hoc tests

We further examined the function of BHLHE40 in breast tumor cells with elevated baseline activation of HIF1A and BHLHE40 using the tamoxifen-resistant (TR) and fulvestrant-resistant (FR) variants of MCF7 cells [25]. As shown in Fig. 3c, mRNA expression levels of BHLHE40 and HIF1A are significantly elevated in TR and FR cells in comparison with parent MCF7 cells. BHLHE40-KD in TR cells substantially increased apoptosis induced by glucose depletion, under both normoxia and hypoxia conditions, as well as reducing the number of viable cells after a 15-day suspension culture in polyHEMA-coated plates (Fig. 3c). Similarly, BHLHE40-KD reduced the number of viable cells in suspension culture and the ability of FR cells to survive glucose depletion (data not shown). Collectively, these observations provide evidence supporting a role for BHLHE40 in promoting survival and migration.

BHLHE40 is required for transcription activation of a set of cytokines and growth factors

To delineate the molecular pathways regulated by BHLHE40, we performed global gene expression analysis of LM-EV and LM BHLHE40-KD cells exposed to hypoxia (1% O2, 6 h or 48 h). The microarray data can be found in the Gene Expression Omnibus database with accession number GSE107300. Overall, the expression levels of 521 and 646 genes in LM-EV cells were altered (fold-change ≥ 1.5 in two independent experiments) by hypoxia at 6 h and 48 h, respectively. BHLHE40-KD abolished the hypoxia-mediated upregulation of 45 (out of 261, 17.2%) and 98 (out of 361, 27.1%) genes at 6 h and 48 h, respectively. In addition, BHLHE40-KD abolished the hypoxia-mediated downregulation of 30 (out of 260, 10.5%) and 44 (out of 285, 15.4%) genes at 6 h and 48 h, respectively. The hypoxia-induced genes that were affected by BHLHE40-KD were over-represented by genes that encode proteins with cytokine or growth factor activities as defined by Gene Ontology annotation GO:0005125 and GO:0008083 (Fisher’s exact test, p <  0.0001; Fig. 4a). The expression of a subset of these genes was also reduced in BHLHE40-KD cells exposed to 1%O2/LG for 4 h compared with EV cell (Fig. 4b). In contrast, hypoxia-induced expression of a panel of the core hypoxia-responsive genes that are known to be directly targeted by HIF1A [14] was not significantly affected by BHLHE40-KD (Fig. 4c). These observations suggest that BHLHE40-KD preferentially reduced the hypoxia-induced expression of a set of cytokines and growth factors but did not cause a global defect in HIF1A-mediated transcription activation. To confirm this notion, we examined the effect of BHLHE40-KO on HIF-mediated expression of reporter luciferase driven by well-characterized HIF1A-binding sites in the promoter regions of LDHA and ITGA6 [26,27,28]. As shown in Fig. 4d, BHLHE40-KD exhibited no significant effect on hypoxia-induced luciferase activities, in the absence or presence of exogenous HIF1A protein. The effect of BHLHE40-KD on the expression of cytokines and growth factors in cells exposed to 1%O2/LG (4 h) was validated by qPCR (Fig. 4e). Consistent with results in BHLHE40-KD cells cultured in vitro, we detected reduced expression of CTGF and HBEGF, at both the mRNA and protein levels, in the LM-BHLHE40-KD primary tumors in comparison with LM-EV tumors established in mouse mammary gland fat pads (Fig. 4f).

Fig. 4
figure 4

BHLHE40-knockdown reduced hypoxia-induced expression of a panel of cytokines and growth factors. a Heatmaps of cytokines and growth factors whose hypoxia-induced expression (1% O2 at 6 h or 48 h, fold-change ≥ 1.5 in two independent experiments) was diminished by BHLHE40-knockdown (KD) in LM cells. The gene expression levels were determined using the Illumina Human HT-12 expression BeadChips. Normalized (quantile normalization) hybridization signals were log2 transformed and standardized by genes across experiment conditions to generate the heatmap. b Heatmaps of a subset of genes list in a whose expression was affected by BHLHE40-KD in LM cells exposed to hypoxia combined with low (1 mM) glucose (1%O2/LG, 4 h). The gene expression levels were determined using the Affymetrix Human Gene 1.0 ST array. c Heatmaps of hypoxia-induced genes whose expression was not significantly affected by BHLHE40-KD in LM cells as determined by the Illumina Human HT-12 expression BeadChips. d Expression of luciferase reporters driven by hypoxia-responsive elements of ITGA6 or LDHA was not affected by BHLHE40 knockout (KO) by CRISPR/Cas9 editing in MDA-MB-231 cells, in the absence or presence of exogenous HIF1A. Luciferase activities were normalized to co-transfected CMV-β-galactosidase and presented as mean ± SD (n = 6). e Expression of genes in control LM empty vector (EV) and LM BHLHE40-KD cells exposed to 1%O2/LG (4 h). mRNA expression levels were determined by qPCR, normalized to RPL13A, and presented as mean ± SD (n = 6). *p <  0.05 (n = 6, 1%O2/LG vs. untreated control), **p <  0.05 (n = 6, KD vs. EV), one-way ANOVA followed by Tukey’s post-hoc tests. f mRNA and protein expression levels of HBEGF and CTGF in primary xenograft tumors, determined by qPCR and immunoblotting, respectively. *p <  0.05 (n = 6, KD vs. EV), Student’s t test. Representative immunoblotting images of three tumors of KD or EV cells are presented

To determine whether BHLHE40-mediated expression of genes encoding cytokine or growth factors is relevant to clinical samples, we analyzed the mRNA expression data of breast tumors in The Cancer Genome Atlas (TCGA) database [30, 31]. The expression of 71.4% (20 out of 28) of these BHLHE40-dependent genes (as shown in Fig. 4a) was found to be positively correlated with BHLHE40 expression, with statistical significance (p <  0.05) in at least one of the four major subtypes of breast tumors (Table 1). This observation provides supporting evidence of a role for BHLHE40 in the expression of these genes in human breast tumors.

Table 1 Correlated expression of BHLHE40 and its putative targets in breast tumors (The Cancer Genome Atlas)

Since hypoxia-induced cytokines and growth factors are commonly exported to the extracellular space by exosomes [32], we sought to determine whether BHLHE40-KD could affect exosome secretion. As shown in Fig. 5, the number of isolated exosomes was significantly reduced in the conditioned medium from MDA-MB-231 BHLHE40-KO and TR-BHLHE40-KD cells in comparison with the corresponding control cells cultured under both normal conditions or exposed to 1%O2/LG for 6 h. The presence of exosomic markers (i.e., CD9, CD81, ALIX, and TSG101) [33, 34]) in the isolated exosomes was confirmed by immunoblotting (Fig. 5). BHLHE40 depletion reduced the protein levels of HBEGF in the purified exosomes (Fig. 5), reflecting the reduced levels of HBEGF mRNA and HBEGF protein in whole cell extracts of BHLHE40-KD or KO cells (Figs. 4e and 5). These observations suggest that BHLHE40 depletion reduced overall exosome secretion and sorting of HBEGF into exosomes.

Fig. 5
figure 5

Exosomic secretion of HBEGF was reduced by BHLHE40 depletion in MDA-MB231 or TR cells. a BHLHE40-knockout (KO) by CRISPR/Cas9 editing in MDA-MB-231 cells reduced the total number of exosomes and the amount of HBEGF protein in exosomes in comparison with empty vector (EV) control cells. b BHLHE40-knockdown (KD) by shRNA in tamoxifen-resistant (TR) cells reduced the total number of exosomes and the amount of HBEGF protein in exosomes in comparison with TR EV cells. Exosomes were purified from conditioned medium of 5 × 106 cells cultured in medium supplemented with exosome-free serum, either under normal culture condition or exposed to 1% O2/low glucose (LG) for 6 h. Exosomes were purified using the ExoQuick-TC solution and quantified under a fluorescent microscope after being labeled with carboxyfluorescein diacetate succinimidyl ester (CFSE) using the Exo-Glow labeling kit, which is designed to exclude background particles. Exosome number in the bar graph is presented as mean number of exosomes per field ± SD (total of nine fields from three independent experiments were examined). The presence of HBEGF and exosome markers in purified exosomes (3 μg protein/lane) or whole cell extracts (WCE; 30 μg protein/lane) were detected by immunoblotting. **p <  0.05 (n = 9, 1%O2/LG vs. control), *p <  0.05 (n = 9, KO or KD vs. EV), one-way ANOVA followed by Tukey’s post-hoc tests

BHLHE40 activates HBEGF transcription by sequestering HDAC1 and HDAC2 from promoter binding

Among the cytokines and growth factors affected by BHLHE40-KD in LM cells, the expression level of HBEGF mRNA is positively correlated with the expression level of BHLHE40 mRNA in all four major subtypes of breast tumors in the TCGA database (Figs. 4 and 5 and Table 1). HBEGF is a heparin-binding epidermal growth factor (EGF)-like growth factor that promotes cell proliferation and invasion through EGF receptor (EGFR) activation [35]. To examine the molecular mechanism underlying BHLHE40-mediated HBEGF transcription, we performed ChIP analysis. BHLHE40 binding to the proximal promoter region of HBEGF was not affected by 1%O2/LG (data not shown), indicating that HBEGF transcription activation was not caused by increased BHLHE40-DNA binding. However, 1%O2/LG treatment reduced binding of HDAC1 and HDAC2 to the HBEGF promoter (Fig. 6a), which is coincident with increased BHLHE40-HDAC1/2 interaction in the soluble cellular fraction, as detected by reciprocal CoIP followed by IB (Fig. 6b). In cells lacking BHLHE40, HDAC1/HDAC2 remained bound to the promoter region of HBEGF after 1%O2/LG treatment (Fig. 6a). This result suggests that BHLHE40 plays a role in facilitating the dissociation of HDAC1/2 from promoters through protein-protein interaction. To examine whether HDAC1/2-DNA binding plays a key role in suppressing transcription of BHLHE40 target genes, we examined the effect of HDAC inhibitors on the mRNA expression of HBEGF, CTGF, and VEGFC. As shown in Fig. 6c, both HDAC2-specific (BRD6688, 10 μM) and pan-HDAC inhibitor (TSA, 2 μM) increased the expression of BHLHE40 target genes in MDA-MB-231 EV and BHLHE40-KD cells, supporting a role for HDAC1/2 in suppressing transcription of BHLHE40 target genes. Taken together, these observations suggest that sequestering HDAC1/2 from DNA binding contributes to BHLHE40-mediated transcription activation.

Fig. 6
figure 6

BHLHE40 activates gene expression by sequestering histone deacetylase (HDAC)1 and HDAC2 from genome DNA binding in MDA-MB-231 cells exposed to hypoxia and low glucose (1%O2/LG, 4 h). a BHLHE40-knockout (KO) diminished dissociation of HDAC1 and HDAC2 from the promoter region of HBEGF in MDA-MB-231 cells exposed to 1%O2/LG (4 h), as determined by chromatin immunoprecipitation (ChIP) followed by qPCR of the HBEGF promoter region (−529 to −372 from the transcription start site). HBEGF promoter binding activity of HDAC1 or HDAC2 was calculated as: (DNA amount in anti-HADC IP complex – DNA amount in control IgG IP complex)/DNA amount in 1% input. *p <  0.05 (n = 6, 1%O2/LG vs. control), **p <  0.05 (n = 6, KO vs. EV). b 1%O2/LG treatment increased interactions between BHLHE40 and HDAC1/2 in the soluble cellular fraction of MDA-MB-231 empty vector (EV) cells. Protein-protein interaction was detected by reciprocal co-immunoprecipitation (IP)/immunoblotting (IB) analysis. c HDAC inhibition induced expression of BHLHE40 target genes. Cells were exposed to hypoxia (1% O2) or HDAC inhibitors (BRD6688 10 μM or TSA 2 μM) for 24 h. mRNA expression levels were determined by qPCR, normalized to RPL13A, and presented as mean ± SD (n = 6). *p <  0.05 (n = 6, treated vs. untreated control cells), one-way ANOVA followed by Tukey’s post-hoc tests

HBEGF acts to promote cell survival and migration

To examine whether BHLHE40-driven HBEGF expression plays a role in EGFR activation to promote cell survival, we examined the phosphorylation status of EGFR and its downstream targets in MDA-MB-231 and TR sublines exposed to 1%O2/GF for 6 h, a condition known to induce apoptosis as shown in Fig. 3. Compared with cells with intact BHLHE40 activity, MDA-MB-231-BHLHE40-KO and TR-BHLHE40-KD cells expressed lower levels of HBEGF mRNA and protein, which was coincident with reduced levels of phosphorylation of EGFR, AKT, and ERK, and increased caspase 9 cleavage (Fig. 7a, b). Next, we examined whether active HBEGF peptide could rescue MDA-MB-231-BHLHE40-KO cells from apoptosis. As shown in Fig. 7c, the addition of HBEGF peptide into the culture medium of cells exposed to 1%O2/GF significantly reduced activation of caspase 3/7. These observations provide evidence supporting a role of HBEGF in promoting cell survival.

Fig. 7
figure 7

BHLHE40 depletion reduced phosphorylation of epidermal growth factor receptor (EGFR), while it increased Caspase 9 cleavage, in cells exposed to glucose depletion and hypoxia (1%O2/GF). a BHLHE40-knockout (KO) by CRISPR/Cas9 editing in MDA-MB-231 and BHLHE40-knockdown (KD) by shRNA in tamoxifen resistant (TR) cells diminished HBEGF induction by 1%O2/GF (6 h). mRNA expression levels were determined by qPCR, normalized to RPL13A, and presented as mean ± SD (n = 6). *p <  0.05 (n = 6, 1%O2/GF vs. control), **p <  0.05 (n = 6, KO vs. EV), one-way ANOVA followed by Tukey’s post-hoc tests. b BHLHE40 depletion reduced EGFR activation, as indicated by reduced phosphorylation of EGFR and its downstream targets (ERK and AKT), while increasing apoptosis, as indicated by detection of cleaved caspase 9. Data from three independent immunoblotting analyses are presented. c HBEGF peptide (10 μg/ml) reduced apoptosis induced by 1%O2/GF (6 h) in MDA-MB-231 BHLHE40-KO cells. Apoptosis was determined by Caspase 3/7 assays. *p <  0.05 (n = 6, 1%O2/GF vs. control), **p <  0.05 (n = 6, HBEGF vs. untreated with HBEGF), one-way ANOVA followed by Tukey’s post-hoc tests

Monolayer scratch was found to induce the expression of HIF1A, BHLHE40, and HBEGF in MDA-MB-231-EV cells (Fig. 8a), implicating a role for the HIF1A-BHLHE40-HBEGF axis in cell migration during wound healing. Using the IncuCyte ZOOM-ImageLock plate system, we demonstrated that BHLHE40-KO substantially diminished the ability of MDA-MB-231 cells to close the wound gaps, which was restored by the addition of HBEGF peptide (Fig. 8b, c). In contrast, a HBEGF-neutralizing antibody [36] inhibited wound healing of LM-EV cells (Fig. 8b, c). To confirm that exosomic HBEGF plays a key role in promoting cell migration, we examined the migratory activities of MDA-MB-231 BHLHE40-KO cells in the presence of conditioned medium or purified exosomes which were collected from the MDA-MB-231 EV cells at 24 h after extensive wound scratch. The transwell migration assay showed that both conditioned medium and purified exosomes from the wounded EV cells increased the migratory activity of BHLHE40-KO cells (Fig. 8d). Together, these observations suggest that HBEGF act downstream of BHLHE40 to promote cell migration.

Fig. 8
figure 8

The HIF-BHLHE40-HBEGF axis plays a role in promoting cell migration during wound healing. a Monolayer scratch increased protein levels of HIF1A, BHLHE40, and HBEGF in MDA-MB-231 empty vector (EV) cells. Intensive scratch wounds were generated using EMD Millipore’s Cell Comb scratch assay kit and immunoblotting was performed 6 h after cells cultured under normoxia (19% O2) or hypoxia (1% O2). b BHLHE40-knockout (KO) by CRISPR/Cas9 editing reduced the migratory activity of MDA-MB-231 cells, which was restored by the addition of HBEGF peptide into the culture medium. In contrast, a HBEGF-neutralizing antibody reduced the migratory activity of MDA-MB-231-EV cells. Real-time assessment of migratory activity after wound scratch was performed using the IncuCyte ZOOM-ImageLock plate system. *p <  0.05 (n = 6, time points 6–24 h, HBEGF vs. untreated), **p <  0.05 (n = 6, time points 6–24 h, anti-HBEGF vs. untreated), one-way ANOVA followed by Tukey’s post-hoc tests. Representative data from two independent experiments with six replicates are presented. c Images of representative wound fields at 0 and 24 h after wound scratch as described in b. d Conditioned medium or purified exosomes from MDA-MB-231-EV cells (24 h after wound scratch) increased migratory activities of MDA-MB-231 BHLHE40-KO cells, as determined by the transwell migration assays. The migrated cells in six fields were imaged and counted under fluorescent microscopy. The results are presented as: percent migration = mean number of cells migrating through the uncoated transwells × 100/mean number of seeded cells; *p <  0.05 (n = 12, vs. untreated control), one-way ANOVA followed by Tukey’s post-hoc tests

To confirm that BHLHE40 and HBEGF are key downstream effectors of HIFs in promoting cell migration, we examined the effect of BHLHE40 overexpression on a MDA-MB-231 subline (HIF-dKO) in which both HIF isoforms (HIF1A and EPAS1) were knocked out by using the CRISPR/Cas9 editing system. Although the HIF1A mRNA expression level is approximately sixfold higher than EPAS1 mRNA in MDA-MB-231 cells according to reported RNAseq data (GSE73526), compensatory activation of EPAS1 could compromise the effect of HIF1A knockout. Therefore, we used HIF-dKO cells to examine whether BHLHE40 overexpression can rescue molecular and phenotypic changes caused by complete elimination of HIF activities. Gene expression analysis by qPCR showed that HIF-dKO reduced baseline and 1%O2/LG-induced expression of BHLHE40, HBEGF, CTGF, and VEGFC mRNA, which was restored by BHLHE40 overexpression (Fig. 9a). In addition, BHLHE40 overexpression reduced cell-cell contact, as shown by cell imaging, and increased the migratory activity of HIF-dKO cells, as determined by transwell assays (Fig. 9b). Immunoblotting analysis confirmed that BHLHE40 overexpression restored expression levels of HBEGF protein in HIF-dKO cells exposed to 1%O2/LG (Fig. 9c). Together, these observations support the notion that BHLHE40 and HBEGF act as key downstream effectors of HIFs to promote cell migration.

Fig. 9
figure 9

Effect of BHLHE40 overexpression on molecular and phenotypic changes caused by HIF1A/EPAS1 double knockout (HIF-dKO) in MDA-MB-231 cells. a BHLHE40 overexpression restored baseline and hypoxia/low glucose (1%O2/LG (6 h))-induced expression of HBEGF, CTGF, TNFSF12, and VEGFC in HIF-dKO cells. The expression levels of mRNA were determined by qPCR, normalized to RPL13A and presented as mean ± SD (n = 6). *p <  0.05 (n = 6, vs. untreated HIF-dKO), **p <  0.05 (n = 6, vs. HIF-dKO exposed to 1%O2/LG for 6 h), one-way ANOVA followed by Tukey’s post-hoc tests. b BHLHE40 overexpression decreased cell-cell contact (as shown by the cell images) and increased migratory activity of HIF-dKO cells exposed to 1%O2/LG (24 h). Migratory activity was determined by transwell assays and presented as mean percentage of migrating cells ± SD (n = 6). *p <  0.05 (n = 6, HIF-dKO vs. control wild-type cells), **p <  0.05 (n = 6, HIF-dKO/BHLHE40 overexpression vs. HIF-dKO), one-way ANOVA followed by Tukey’s post-hoc tests. c BHLHE40 overexpression restored expression of HBEGF proteins in HIF-dKO cells exposed to 1%O2/GF (6 h). Proteins were detected by immunoblotting. GAPDH and TBP were used as loading control for whole cell extract (WCE) or nuclear extract (NE), respectively

High expression of BHLHE40 and HBEGF is associated with poor prognosis of breast cancer

Having established a role of the BHLHE40-HBEGF axis in enhancing cell survival and migration, we sought to examine the association of BHLHE40 and HBEGF with clinical characteristics of breast tumors using the gene expression data in the Kaplan-Meier plotter database, which contains the Affymetrix microarray expression data of 2178 breast cancer patients [37]. We found that high expression of BHLHE40 or HBEGF is significantly associated with shorter interval of relapse-free survival (RFS) among patients diagnosed with triple-negative breast cancer (TNBC; n = 255) and patients treated with chemotherapy (n = 602) (Fig. 10). However, BHLHE40 and HBEGF are not poor prognostic markers for patients with estrogen receptor-positive tumors or patients treated with endocrine therapy. In addition, we analyzed the association of BHLHE40 and HBEGF with overall survival (OS) of TNBC using the METABRIC dataset in the cBioPortal for Cancer Genomics. High HBEGF expression was found to be associated with a short interval of OS (Fig. 10c). Although TNBC with higher expression of BHLHE40 tends to have a shorter interval of OS, this correlation did not reach statistical significance (Fig. 10c). These findings suggest that activation of the BHLHE40-HBEGF pathway contributes to aggressive behaviors of TNBC and chemoresistance.

Fig. 10
figure 10

High expression of BHLHE40 and HBEGF is associated with poor prognosis of breast cancer. a High expression of BHLHE40 and HBEGF is associated with short interval of relapse-free survival (RFS) of patients diagnosed with triple-negative breast cancer (TNBC; n = 255). The gene expression data and patient information were obtained from the Kaplan-Meier plotter database. b High expression of BHLHE40 and HBEGF is associated with a short interval of RFS among patients treated with chemotherapy (n = 602). The gene expression data and patient information were obtained from the Kaplan-Meier plotter database. c High expression of HBEGF is associated with a short interval of overall survival (OS) of patients diagnosed with TNBC (n = 150). The gene expression data and patient information were obtained from the METABRIC database in the cBioPortal for Cancer Genomics

Discussion

Breast cancer metastasis is the major cause of death in breast cancer patients. Adaptation to hypoxia is a driving force of metastatic progression and drug resistance [3]. Proteins secreted by tumor cells under hypoxia promote metastasis by altering tumor cell behaviors and modifying the tumor microenvironment [2]. Therefore, the regulation of hypoxia-driven protein secretion is currently under intense investigation. In this study, we report a novel role of BHLHE40, a transcription factor directly targeted by HIF1A, in regulating exosomic release of HBEGF. Our results suggest that the HIF-BHLHE40-HBEGF axis constitutes an important signaling mechanism to promote metastasis of breast tumors.

Exosomes are 40- to 100-nm vesicles that originate from the endocytic compartment. Exosomes contain a wide range of proteins, lipids, mRNAs, and microRNAs (miRNAs) that reflect the molecular contents of the parental cells [32]. Compared with normal cells, cancer cells exhibit higher activity of exosome secretion, which is further augmented by stress conditions including TP53 activation, alteration of intracellular calcium levels, senescence, hypoxia, and acidosis [38]. Exosomes released by tumor cells have been reported to contain cytokines and growth factors that promote metastasis and chemoresistance [38,39,40]. However, the precise molecular mechanism governing the release of exosomes remains elusive. This study suggests that BHLHE40 acts as a key downstream effector of HIFs to activate transcription and subsequent exosome secretion of a set of cytokines and growth factors.

BHLHE40 was previously described as a transcriptional repressor that binds to the class B E-box (CACGTG) and recruits HDAC1 and HDAC2 to block transcription [41]. BHLHE40 activation has been linked to cell cycle arrest, senescence, differentiation, and apoptosis [42,43,44]. On the other hand, emerging evidence supports a role for BHLHE40 in transcription activation and promoting cell survival. For instance, BHLHE40 was reported to activate transcription of pro-survival factors in tumor cells, including BIRC5 and DeltaNp63 [45, 46]. In addition, BHLHE40 was reported to activate the transcription of a panel of cytokines required for activation of murine CD4+ T cells [47, 48]. Which factors determine the selectivity of BHLHE40 to suppress or activate transcription remains undefined. The BHLHE40-mediated transcription activation of DeltaNp63 was shown to depend on its direct interaction with HDAC2 [45]. In agreement with this observation, our results suggest that BHLHE40 activates HBEGF transcription by sequestering HDAC1/2 from DNA binding. It remains to be determined whether interfering with HDAC1/2-DNA binding is a general mechanism responsible for BHLHE40-mediated transcription activation.

Elevated EGFR activation is known to promote survival, proliferation, and invasion of tumor cells under hypoxia, and multiple mechanisms have been linked to hypoxia-induced EGFR activation [49]. For example, EPAS1 activation by hypoxia was shown to increase EGFR mRNA translation [50]. Hypoxia-mediated activation of metalloproteases (e.g., ADAM12 and ADAM17) was reported to activate EGFR by increasing ectodomain shedding of HBEGF [51, 52]. Our study provides a novel aspect of EGFR activation through the BHLHE40-HBEGF axis. In addition to autocrine or paracrine effects within tumor cells, exosomic release of HBEGF might exert paracrine effects to remodel tumor stroma or endocrine effects to prime distant metastatic niches [53].

In conclusion, this study provides evidence supporting an essential role of BHLHE40 in exosomic release of HBEGF, a critical pro-survival and pro-metastasis factor. The clinical relevance of our findings is evidenced by the fact that the elevated expression of BHLHE40 and HBEGF in breast tumors is associated with poor prognosis of patients with TNBC and chemoresistance. Therapeutic intervention targeting the BHLHE40-HBEGF axis may represent an effective approach to combat hypoxia-driven drug resistance and metastasis.

Conclusion

Hypoxia-induced activation of BHLHE40 plays a key role in promoting cell survival and metastasis by modulating exosomic secretion of HBEGF.

Abbreviations

BHLHE40:

Basic helix-loop-helix family member E40

ChIP:

Chromatin immunoprecipitation

CoIP:

Co-immunoprecipitation

EGF:

Epidermal growth factor

EGFR:

Epidermal growth factor receptor

EV:

Empty vectors

FR:

Fulvestrant-resistant

GF:

Glucose depletion

HBEGF:

Heparin-binding epidermal growth factor

HDAC:

Histone deacetylases

HIF:

Hypoxia inducible factor

HIF-dKO:

HIF1A and EPAS1 double knockout

IB:

Immunoblotting

KD:

Knockdown

KO:

Knockout

LG:

Low glucose

LM:

Lung metastatic cell line derived from MDA-MB-231

NSG:

NOD.Cg-Prkdcscid Il2rgtm1Wjl/SzJ

OS:

Overall survival

qPCR:

Quantitative polymerase chain reaction

RFS:

Relapse-free survival

shRNA:

Short hairpin RNA

TCGA:

The Cancer Genome Atlas

TNBC:

Triple-negative breast cancer

TR:

Tamoxifen-resistant

References

  1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2018. CA Cancer J Clin. 2018;68(1):7–30.

    Article  Google Scholar 

  2. Gilkes DM, Semenza GL, Wirtz D. Hypoxia and the extracellular matrix: drivers of tumour metastasis. Nat Rev Cancer. 2014;14(6):430–9.

    Article  CAS  Google Scholar 

  3. Rankin EB, Giaccia AJ. Hypoxic control of metastasis. Science. 2016;352(6282):175–80.

    Article  CAS  Google Scholar 

  4. Petrova V, Annicchiarico-Petruzzelli M, Melino G, Amelio I. The hypoxic tumour microenvironment. Oncogenesis. 2018;7(1):10.

    Article  Google Scholar 

  5. Sormendi S, Wielockx B. Hypoxia pathway proteins as central mediators of metabolism in the tumor cells and their microenvironment. Front Immunol. 2018;9:40.

    Article  Google Scholar 

  6. Muz B, de la Puente P, Azab F, Azab AK. The role of hypoxia in cancer progression, angiogenesis, metastasis, and resistance to therapy. Hypoxia. 2015;3:83–92.

    Article  Google Scholar 

  7. Semenza GL. The hypoxic tumor microenvironment: a driving force for breast cancer progression. Biochim Biophys Acta. 2016;1863(3):382–91.

    Article  CAS  Google Scholar 

  8. Weidle HU, Birzele F, Kollmorgen G, RÜGer R. The multiple roles of exosomes in metastasis. Cancer Genomics Proteomics. 2017;14(1):1–16.

    Article  CAS  Google Scholar 

  9. Harris DA, Patel SH, Gucek M, Hendrix A, Westbroek W, Taraska JW. Exosomes released from breast cancer carcinomas stimulate cell movement. PLoS One. 2015;10(3):e0117495.

    Article  Google Scholar 

  10. O'Driscoll L. Expanding on exosomes and ectosomes in cancer. N Engl J Med. 2015;372(24):2359–62.

    Article  CAS  Google Scholar 

  11. Melo SA, Sugimoto H, O'Connell JT, Kato N, Villanueva A, Vidal A, Qiu L, Vitkin E, Perelman LT, Melo CA, et al. Cancer exosomes perform cell-independent microRNA biogenesis and promote tumorigenesis. Cancer Cell. 2014;26(5):707–21.

    Article  CAS  Google Scholar 

  12. King HW, Michael MZ, Gleadle JM. Hypoxic enhancement of exosome release by breast cancer cells. BMC Cancer. 2012;12:421.

    Article  CAS  Google Scholar 

  13. Pawlus MR, Wang L, Hu CJ. STAT3 and HIF1alpha cooperatively activate HIF1 target genes in MDA-MB-231 and RCC4 cells. Oncogene. 2014;33(13):1670–9.

    Article  CAS  Google Scholar 

  14. Villar D, Ortiz-Barahona A, Gomez-Maldonado L, Pescador N, Sanchez-Cabo F, Hackl H, Rodriguez BA, Trajanoski Z, Dopazo A, Huang TH, et al. Cooperativity of stress-responsive transcription factors in core hypoxia-inducible factor binding regions. PLoS One. 2012;7(9):e45708.

    Article  CAS  Google Scholar 

  15. Miyazaki K, Kawamoto T, Tanimoto K, Nishiyama M, Honda H, Kato Y. Identification of functional hypoxia response elements in the promoter region of the DEC1 and DEC2 genes. J Biol Chem. 2002;277(49):47014–21.

    Article  CAS  Google Scholar 

  16. Khurana P, Sugadev R, Jain J, Singh SB. HypoxiaDB: a database of hypoxia-regulated proteins. Database (Oxford). 2013;2013:bat074.

    Article  Google Scholar 

  17. Xiong J, Yang H, Luo W, Shan E, Liu J, Zhang F, Xi T, Yang J. The anti-metastatic effect of 8-MOP on hepatocellular carcinoma is potentiated by the down-regulation of bHLH transcription factor DEC1. Pharmacol Res. 2016;105:121–33.

    Article  CAS  Google Scholar 

  18. Wu Y, Sato F, Yamada T, Bhawal UK, Kawamoto T, Fujimoto K, Noshiro M, Seino H, Morohashi S, Hakamada K, et al. The BHLH transcription factor DEC1 plays an important role in the epithelial-mesenchymal transition of pancreatic cancer. Int J Oncol. 2012;41(4):1337–46.

    Article  CAS  Google Scholar 

  19. Chakrabarti J, Turley H, Campo L, Han C, Harris AL, Gatter KC, Fox SB. The transcription factor DEC1 (stra13, SHARP2) is associated with the hypoxic response and high tumour grade in human breast cancers. Br J Cancer. 2004;91(5):954–8.

    Article  CAS  Google Scholar 

  20. Nakashima A, Kawamoto T, Honda KK, Ueshima T, Noshiro M, Iwata T, Fujimoto K, Kubo H, Honma S, Yorioka N, et al. DEC1 modulates the circadian phase of clock gene expression. Mol Cell Biol. 2008;28(12):4080–92.

    Article  CAS  Google Scholar 

  21. Nishiyama Y, Goda N, Kanai M, Niwa D, Osanai K, Yamamoto Y, Senoo-Matsuda N, Johnson RS, Miura S, Kabe Y, et al. HIF-1alpha induction suppresses excessive lipid accumulation in alcoholic fatty liver in mice. J Hepatol. 2012;56(2):441–7.

    Article  CAS  Google Scholar 

  22. Chung SY, Kao CH, Villarroya F, Chang HY, Chang HC, Hsiao SP, Liou GG, Chen SL. Bhlhe40 represses PGC-1alpha activity on metabolic gene promoters in myogenic cells. Mol Cell Biol. 2015;35(14):2518–29.

    Article  CAS  Google Scholar 

  23. Kanda M, Yamanaka H, Kojo S, Usui Y, Honda H, Sotomaru Y, Harada M, Taniguchi M, Suzuki N, Atsumi T, et al. Transcriptional regulator Bhlhe40 works as a cofactor of T-bet in the regulation of IFN-gamma production in iNKT cells. Proc Natl Acad Sci U S A. 2016;113(24):E3394–402.

    Article  CAS  Google Scholar 

  24. Fan M, Krutilina R, Sun J, Sethuraman A, Yang CH. Wu Z-h, Yue J, Pfeffer LM: comprehensive analysis of microRNA (miRNA) targets in breast cancer cells. J Biol Chem. 2013;288(38):27480–93.

    Article  CAS  Google Scholar 

  25. Fan M, Yan PS, Hartman-Frey C, Chen L, Paik H, Oyer SL, Salisbury JD, Cheng AS, Li L, Abbosh PH, et al. Diverse gene expression and DNA methylation profiles correlate with differential adaptation of breast cancer cells to the antiestrogens tamoxifen and fulvestrant. Cancer Res. 2006;66(24):11954–66.

    Article  CAS  Google Scholar 

  26. Sethuraman A, Brown M, Seagroves TN, Wu ZH, Pfeffer LM, Fan M. SMARCE1 regulates metastatic potential of breast cancer cells through the HIF1A/PTK2 pathway. Breast Cancer Res. 2016;18(1):81.

    Article  Google Scholar 

  27. Patsialou A, Wang Y, Lin J, Whitney K, Goswami S, Kenny PA, Condeelis JS. Selective gene-expression profiling of migratory tumor cells in vivo predicts clinical outcome in breast cancer patients. Breast Cancer Res. 2012;14(5):R139.

    Article  CAS  Google Scholar 

  28. Krutilina R, Sun W, Sethuraman A, Brown M, Seagroves TN, Pfeffer LM, Ignatova T, Fan M. MicroRNA-18a inhibits hypoxia-inducible factor 1alpha activity and lung metastasis in basal breast cancers. Breast Cancer Res. 2014;16(4):R78.

    Article  Google Scholar 

  29. Whelan KA, Schwab LP, Karakashev SV, Franchetti L, Johannes GJ, Seagroves TN, Reginato MJ. The oncogene HER2/neu (ERBB2) requires the hypoxia-inducible factor HIF-1 for mammary tumor growth and anoikis resistance. J Biol Chem. 2013;288(22):15865–77.

    Article  CAS  Google Scholar 

  30. Cancer Genome Atlas N. Comprehensive molecular portraits of human breast tumours. Nature. 2012;490(7418):61–70.

    Article  Google Scholar 

  31. Ciriello G, Gatza ML, Beck AH, Wilkerson MD, Rhie SK, Pastore A, Zhang H, McLellan M, Yau C, Kandoth C, et al. Comprehensive molecular portraits of invasive lobular breast cancer. Cell. 2015;163(2):506–19.

    Article  CAS  Google Scholar 

  32. Keerthikumar S, Chisanga D, Ariyaratne D, Al Saffar H, Anand S, Zhao K, Samuel M, Pathan M, Jois M, Chilamkurti N, et al. ExoCarta: a web-based compendium of exosomal cargo. J Mol Biol. 2016;428(4):688–92.

    Article  CAS  Google Scholar 

  33. Yáñez-Mó M, Siljander PRM, Andreu Z, Zavec AB, Borràs FE, Buzas EI, Buzas K, Casal E, Cappello F, Carvalho J, et al. Biological properties of extracellular vesicles and their physiological functions. J Extracellular Vesicles. 2015;4. https://doi.org/10.3402/jev.v3404.27066.

  34. Yáñez-Mó M, Siljander PR, Andreu Z, Zavec AB, Borràs FE, Buzas EI, Buzas K, Casal E, Cappello F, Carvalho J, et al. Biological properties of extracellular vesicles and their physiological functions. J Extracell Vesicles. 2015;4:27066.

    Article  Google Scholar 

  35. Prenzel N, Zwick E, Daub H, Leserer M, Abraham R, Wallasch C, Ullrich A. EGF receptor transactivation by G-protein-coupled receptors requires metalloproteinase cleavage of proHB-EGF. Nature. 1999;402(6764):884–8.

    Article  CAS  Google Scholar 

  36. Rubin JS, Chan AM, Bottaro DP, Burgess WH, Taylor WG, Cech AC, Hirschfield DW, Wong J, Miki T, Finch PW, et al. A broad-spectrum human lung fibroblast-derived mitogen is a variant of hepatocyte growth factor. Proc Natl Acad Sci U S A. 1991;88(2):415–9.

    Article  CAS  Google Scholar 

  37. Gyorffy B, Lanczky A, Eklund AC, 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 1,809 patients. Breast Cancer Res Treat. 2010;123(3):725–31.

    Article  Google Scholar 

  38. Atretkhany KSN, Drutskaya MS, Nedospasov SA, Grivennikov SI, Kuprash DV. Chemokines, cytokines and exosomes help tumors to shape inflammatory microenvironment. Pharmacol Ther. 2016;168:98–112.

    Article  CAS  Google Scholar 

  39. Kucharzewska P, Christianson HC, Welch JE, Svensson KJ, Fredlund E, Ringner M, Morgelin M, Bourseau-Guilmain E, Bengzon J, Belting M. Exosomes reflect the hypoxic status of glioma cells and mediate hypoxia-dependent activation of vascular cells during tumor development. Proc Natl Acad Sci U S A. 2013;110(18):7312–7.

    Article  CAS  Google Scholar 

  40. Keerthikumar S, Gangoda L, Liem M, Fonseka P, Atukorala I, Ozcitti C, Mechler A, Adda CG, Ang CS, Mathivanan S. Proteogenomic analysis reveals exosomes are more oncogenic than ectosomes. Oncotarget. 2015;6(17):15375–96.

    Article  Google Scholar 

  41. St-Pierre B, Flock G, Zacksenhaus E, Egan SE. Stra13 homodimers repress transcription through class B E-box elements. J Biol Chem. 2002;277(48):46544–51.

    Article  CAS  Google Scholar 

  42. Bhawal UK, Sato F, Arakawa Y, Fujimoto K, Kawamoto T, Tanimoto K, Ito Y, Sasahira T, Sakurai T, Kobayashi M, et al. Basic helix-loop-helix transcription factor DEC1 negatively regulates cyclin D1. J Pathol. 2011;224(3):420–9.

    Article  CAS  Google Scholar 

  43. Seino H, Wu Y, Morohashi S, Kawamoto T, Fujimoto K, Kato Y, Takai Y, Kijima H. Basic helix-loop-helix transcription factor DEC1 regulates the cisplatin-induced apoptotic pathway of human esophageal cancer cells. Biomed Res. 2015;36(2):89–96.

    Article  CAS  Google Scholar 

  44. Bi H, Li S, Qu X, Wang M, Bai X, Xu Z, Ao X, Jia Z, Jiang X, Yang Y, et al. DEC1 regulates breast cancer cell proliferation by stabilizing cyclin E protein and delays the progression of cell cycle S phase. Cell Death Dis. 2015;6:e1891.

    Article  CAS  Google Scholar 

  45. Qian Y, Jung YS, Chen X. DeltaNp63, a target of DEC1 and histone deacetylase 2, modulates the efficacy of histone deacetylase inhibitors in growth suppression and keratinocyte differentiation. J Biol Chem. 2011;286(14):12033–41.

    Article  CAS  Google Scholar 

  46. Li Y, Xie M, Yang J, Yang D, Deng R, Wan Y, Yan B. The expression of antiapoptotic protein survivin is transcriptionally upregulated by DEC1 primarily through multiple sp1 binding sites in the proximal promoter. Oncogene. 2006;25(23):3296–306.

    Article  CAS  Google Scholar 

  47. Martinez-Llordella M, Esensten JH, Bailey-Bucktrout SL, Lipsky RH, Marini A, Chen J, Mughal M, Mattson MP, Taub DD, Bluestone JA. CD28-inducible transcription factor DEC1 is required for efficient autoreactive CD4+ T cell response. J Exp Med. 2013;210(8):1603–19.

    Article  CAS  Google Scholar 

  48. Cowley GS, Weir BA, Vazquez F, Tamayo P, Scott JA, Rusin S, East-Seletsky A, Ali LD, Gerath WF, Pantel SE, et al. Parallel genome-scale loss of function screens in 216 cancer cell lines for the identification of context-specific genetic dependencies. Sci Data. 2014;1:140035.

    Article  CAS  Google Scholar 

  49. Chen Y, Henson ES, Xiao W, Huang D, McMillan-Ward EM, Israels SJ, Gibson SB. Tyrosine kinase receptor EGFR regulates the switch in cancer cells between cell survival and cell death induced by autophagy in hypoxia. Autophagy. 2016;12(6):1029–46.

    Article  CAS  Google Scholar 

  50. Franovic A, Gunaratnam L, Smith K, Robert I, Patten D, Lee S. Translational up-regulation of the EGFR by tumor hypoxia provides a nonmutational explanation for its overexpression in human cancer. Proc Natl Acad Sci U S A. 2007;104(32):13092–7.

    Article  CAS  Google Scholar 

  51. Diaz B, Yuen A, Iizuka S, Higashiyama S, Courtneidge SA. Notch increases the shedding of HB-EGF by ADAM12 to potentiate invadopodia formation in hypoxia. J Cell Biol. 2013;201(2):279–92.

    Article  CAS  Google Scholar 

  52. Wang XJ, Feng CW, Li M. ADAM17 mediates hypoxia-induced drug resistance in hepatocellular carcinoma cells through activation of EGFR/PI3K/Akt pathway. Mol Cell Biochem. 2013;380(1–2):57–66.

    Article  CAS  Google Scholar 

  53. Yotsumoto F, Tokunaga E, Oki E, Maehara Y, Yamada H, Nakajima K, Nam SO, Miyata K, Koyanagi M, Doi K, et al. Molecular hierarchy of heparin-binding EGF-like growth factor-regulated angiogenesis in triple-negative breast cancer. Mol Cancer Res. 2013;11(5):506–17.

    Article  CAS  Google Scholar 

Download references

Funding

This work was supported by a National Institutes of Health grant (CA197206 to MF).

Author information

Authors and Affiliations

Authors

Contributions

AS carried out the majority of the experiments and drafted the manuscript. MB and RK participated in the cell-based assays and manuscript preparation. TNS, Z-HW, and LMP participated in experiment design and manuscript preparation. MF is responsible for the overall experiment design and manuscript preparation. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Meiyun Fan.

Ethics declarations

Ethics approval

All studies with mouse models were performed according to the ethical regulations as prescribed in the Research Compliance guidelines of the University of Tennessee Health Science Center.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sethuraman, A., Brown, M., Krutilina, R. et al. BHLHE40 confers a pro-survival and pro-metastatic phenotype to breast cancer cells by modulating HBEGF secretion. Breast Cancer Res 20, 117 (2018). https://doi.org/10.1186/s13058-018-1046-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13058-018-1046-3

Keywords