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

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% O₂), hypoxia in combination with low (1 mg/ml) glucose (1%O₂/LG), and hypoxia in combination with glucose depletion (1%O₂/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 = (width² × 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 × 10⁵) 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% O₂) 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 × 10⁶ 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% O₂ 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% O₂) in combination with low glucose (1 mM; 1%O₂/LG), a condition frequently encountered by cells in solid tumors, total RNA of cells exposed to 1%O₂/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.

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% O₂, 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 × 10⁵ 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.

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 × 10⁵) 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.

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% O₂; 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%O₂/GF, 6 h), as evidenced by the appearance of apoptotic morphology and activation of caspase3/7 (Fig. 3b, lower panel).

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% O₂, 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%O₂/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–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%O₂/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).

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%O₂/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.

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%O₂/LG (data not shown), indicating that HBEGF transcription activation was not caused by increased BHLHE40-DNA binding. However, 1%O₂/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%O₂/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.

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%O₂/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%O₂/GF significantly reduced activation of caspase 3/7. These observations provide evidence supporting a role of HBEGF in promoting cell survival.

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.

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%O₂/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%O₂/LG (Fig. 9c). Together, these observations support the notion that BHLHE40 and HBEGF act as key downstream effectors of HIFs to promote cell migration.

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.