Obesity is associated with poor prognosis in triple-negative breast cancer (TNBC). Preclinical models of TNBC were used to test the hypothesis that increased leptin signaling drives obesity-associated TNBC development by promoting cancer stem cell (CSC) enrichment and/or epithelial-to-mesenchymal transition (EMT). MMTV-Wnt-1 transgenic mice, which develop spontaneous basal-like, triple-negative mammary tumors, received either a control diet (10% kcal from fat) or a diet-induced obesity regimen (DIO, 60% kcal from fat) for up to 42 weeks (n = 15/group). Mice were monitored for tumor development and euthanized when tumor diameter reached 1.5 cm. Tumoral gene expression was assessed via RNA sequencing (RNA-seq). DIO mice had greater body weight and percent body fat at termination than controls. DIO mice, versus controls, demonstrated reduced survival, increased systemic metabolic and inflammatory perturbations, upregulated tumoral CSC/EMT gene signature, elevated tumoral aldehyde dehydrogenase activity (a CSC marker), and greater leptin signaling. In cell culture experiments using TNBC cells (murine: E-Wnt and M-Wnt; human: MDA-MB-231), leptin enhanced mammosphere formation, and media supplemented with serum from DIO versus control mice increased cell viability, migration, invasion, and CSC- and EMT-related gene expression, including Foxc2, Twist2, Vim, Akt3, and Sox2. In E-Wnt cells, knockdown of leptin receptor ablated these procancer effects induced by DIO mouse serum. These findings indicate that increased leptin signaling is causally linked to obesity-associated TNBC development by promoting CSC enrichment and EMT.

Implications: Leptin-associated signals impacting CSC and EMT may provide new targets and intervention strategies for decreasing TNBC burden in obese women. Mol Cancer Res; 16(5); 869–79. ©2018 AACR.

This article is featured in Highlights of This Issue, p. 743

The prevalence of obesity, an established risk factor for several breast cancer subtypes, including triple-negative breast cancer (TNBC), is over 40% in U.S. women (1). Obesity is also associated with higher breast cancer–specific and overall mortality in patients with breast cancer, regardless of subtype (2–4). Moreover, obese patients do not respond as well as normal weight patients to several classes of chemotherapy (5–9). Unfortunately, the mechanisms underlying the obesity–TNBC link are incompletely understood, hampering efforts to decrease the burden of breast cancer in obese women.

One possible mediator of obesity-associated TNBC development, progression, and chemotherapy resistance is cancer stem cell (CSC) enrichment. CSCs are an intratumoral subpopulation of cancer cells with stem cell–like properties, including self-renewal and multipotent differentiation, and are posited as primary drivers of tumor initiation and progression (10, 11). CSC enrichment is strongly associated with poor prognosis in patients with breast cancer (10). Obesity in preclinical transplant models of TNBC increases tumor growth, CSC-related gene expression, and markers of epithelial-to-mesenchymal transition (EMT), a key developmental program linked with CSC enrichment and metastasis (11). The procancer effects of obesity in these models are associated with increased circulating leptin levels (12). Leptin, an adipokine produced primarily by white adipose tissue, is positively correlated with body fat levels and regulates food intake, inflammation, cell differentiation, and proliferation. Elevated levels of systemic leptin, intratumoral leptin, and leptin receptor expression confer a worse breast cancer prognosis (13–16).

The MMTV-Wnt-1 transgenic mouse is an attractive, relevant model of basal-like breast cancer, the most common form of TNBC. Although Wnt-1 is not implicated in human breast tumorigenesis, activation of the Wnt/β-catenin pathway is required for the tumorigenic behavior of basal-like breast cancer cells (17), and spontaneous tumors in MMTV-Wnt-1 mice display pathologic and transcriptional characteristics of human basal-like TNBC (11). Using MMTV-Wnt-1 mice and in vitro models of basal-like breast cancer, we tested the hypothesis that increased leptin signaling drives obesity-associated TNBC development by promoting CSC enrichment and EMT.

In vivo MMTV-Wnt-1 transgenic mouse study

Animal studies and procedures were approved and monitored by the University of Texas Institutional Animal Care and Use Committee. Female MMTV-Wnt-1 mice on a C57BL/6 background (n = 30) were obtained from the Hursting laboratory's breeding colony and ovariectomized at 12 weeks of age. Following a 1-week recovery period, the mice were randomized to receive one of two diets, both fed ad libitum: a control regimen or a diet-induced obesity (DIO) regimen (n = 15/diet group) as described previously (18). Mice were followed for up to 42 weeks and euthanized when tumor diameter was >1.5 cm in any direction or at study termination, whichever occurred first. Mice that died or were euthanized for nontumor-related cause prior to study termination were censored from the survival analysis. Mice were palpated weekly for mammary tumors, and tumors were measured (2 dimensions) with calipers twice weekly after detection. Mice were fasted for 4 to 6 hours before blood collection. Tumor (when present), distal mammary fat pad, and blood were collected, processed, and stored as described previously (19). Body fat and lean mass were assessed after euthanization using a Lunar PIXImus Dual Emission X-Ray Absorptiometer (GE Medical Systems). One DIO mouse and 4 control mice died from nontumor-related causes, were censored at the time of death in body weight and survival analyses, and excluded from serum and tissue analyses.

Serum hormone, adipokine, and cytokine analyses

Analyses of serum levels of several obesity-associated hormones, adipokines, and cytokines were performed on 10 randomly selected mice per group. Serum IGF-1 was measured via a Mouse Magnetic Luminex Screening Assay (R&D Systems). Other serum proteins were measured via Bio-Plex Pro Multiplex Immunoassays (Bio-Rad, Inc.). All assays were analyzed on a Bio-Plex MAGPIX Multiplex Reader (Bio-Rad, Inc.).

IHC analyses

Paraffin-embedded tumor tissue from 6 randomly selected mice per group was cut into 4-μm thick sections for IHC analysis. Tumor tissue was stained, processed, and analyzed as described previously (11) with the following primary antibodies: Ki67 (Bethyl Laboratories #IHC-00375), CD31 (BD Pharmingen #550274), E-cadherin (Santa Cruz Biotechnology #sc-7870), and vimentin (Abcam #ab92547).

RNA sequencing and Ingenuity Pathway Analysis

Total RNA was isolated from flash-frozen tumor samples using TRIzol Reagent (Sigma-Aldrich) according to manufacturer's instructions. RNA quality was assessed by Agilent 2100 Bioanalyzer, and randomly selected samples from control (n = 5) and DIO (n = 6) mice were prepared for sequencing using the Illumina TruSeq RNA Library Preparation Kit. RNA libraries were then sequenced on the Illumina HiSeq 2000 instrument (Illumina Inc.). Sequencing reads were quality assessed and trimmed for any remaining sequencing adaptor, and reads smaller than 50 bp were removed. RNA sequencing reads were aligned to Mouse Genome Ensembl GRCm38 using TopHat (version 2.1.0), and the differential levels of transcripts were quantified by Cuffdiff/Cufflink 2.2.1. Normalized fragments per kilobase of exon per million fragments mapped (FPKM) were generated by Cufflinks as representative of gene expression level. With principle component analysis, 3 outlier samples (1 control and 2 DIO) were removed, and 4 samples per group were subjected to further differential expression and pathway enrichment analysis. Genes with an adjusted two-tailed P value of <0.05 and greater than 20.75-fold change (∼1.7-fold) in expression were considered differentially expressed. These genes were then entered into Ingenuity Pathway Analysis (Qiagen) to identify modulated signaling pathways and potential upstream regulators. The RNA sequencing data have been deposited in Sequence Read Archive (accession number SRP129022).

Quantitative RT-PCR analyses

Total RNA was isolated from tumor and distal mammary pad tissues, reverse transcribed, and assayed in triplicate as described previously (19).

Aldehyde dehydrogenase activity assay

A colorimetric aldehyde dehydrogenase (ALDH) Activity Assay Kit (Abcam) was used to assess tumor ALDH enzymatic activity in 6 randomly selected mice per group. Levels of NADH (produced by ALDH oxidation of acetaldehyde) were normalized to protein concentration of the tumor lysate, measured via Bradford protein assay, to determine ALDH activity.

Cell lines

Two mouse mammary tumor cell lines isolated from MMTV-Wnt-1 mice in 2010 by the Hursting laboratory and found to cluster with the basal-like (E-Wnt) and claudin-low (M-Wnt) breast cancer subtypes (11), and human MDA-MB-231 breast cancer cells (ATCC #HTB-26, obtained in 2011), which cluster with the claudin-low subtype, were used in in vitro studies. The cell lines were maintained in RPMI1640 media (GIBCO Life Technologies) supplemented with 10% FBS, 10 mmol/L HEPES buffer, and 2 mmol/L l-glutamine (complete media). All experiments were performed on cells maintained at <30 passages and within 10 passages of cell line removal from liquid nitrogen storage. The cell lines were found to be negative for mycoplasma by the University of North Carolina's Tissue Culture Facility and were authenticated via species identification and karyotyping analysis at the T.C. Shu Molecular Cytogenetics Core at the University of Texas MD Anderson Cancer Center (Houston, TX) in 2016.

An Ob-R shRNA Plasmid (sc-36116-SH, Santa Cruz Biotechnology, Dallas, TX) was transfected into E-Wnt cells using FuGENE 6 to generate Lepr knockdown E-Wnt cell lines. Stably transfected cells were selected based on growth in media with 5 μg/ml puromycin dihydrochloride (Sigma-Aldrich). Two clones with the lowest Lepr expression (EWnt-L1 and EWnt-L2), as determined by quantitative RT-PCR, were maintained in complete media plus 2.5 μg/ml puromycin dihydrochloride. An E-Wnt cell line stably transfected with the scrambled Control shRNA Plasmid-A (sc-108060, Santa Cruz Biotechnology) was also generated (EWnt-S) and maintained via the same methodology. The parental E-Wnt cell line (EWnt-P) used to generate these 3 cell lines was also utilized for comparison in the in vitro experiments.

Mammosphere assays

E-Wnt, M-Wnt, and MDA-MB-231 cells were plated for the mammosphere assays as described previously (12). Serial passage of the mammospheres, a common approach in the CSC field (20), was performed to assess whether leptin treatment further promotes mammary CSC self-renewal in secondary spheres. For propagation 1 (P1), cells were treated with vehicle (sterile H2O) or leptin (100, 200, and 400 ng/mL, final concentrations) for 7 days; then, the number of spheres/well quantified. Previous reports examining the impact of leptin on CSC markers in multiple types of cancer cells have used leptin concentrations in this range (12, 21–23). Propagation 2 (P2) was then initiated by dissociating the mammospheres grown under each treatment condition and replating the cells in new media with the same treatments for an additional 7 days. The number of spheres was quantified at the end of P2, and representative images (×10 magnification) were captured on an EVOS FL Auto 2 Cell Imaging System (Invitrogen). Mammosphere assays of EWnt-P, EWnt-S, EWnt-L1, and EWnt-L2 cells used the same plating methodology and media, but only one propagation period.

Generation and characterization of mouse serum for in vitro studies

Female 8-week-old wild-type C57BL/6 mice (n = 202) were acclimated on the control diet regimen for 1 week, and then randomized to remain on control diet (n = 110) or change to the DIO regimen (n = 92). The mice were maintained on these diets for 14 weeks, and then euthanized. Blood was collected by cardiac puncture and serum isolated. The serum samples were pooled by diet group, and levels of serum hormones, adipokines, and cytokines measured in each pool as described above. The remaining serum was aliquoted and stored at −80°C for use in the in vitro experiments.

In vitro gene expression analyses

To assess basal gene expression in vitro, EWnt-P, EWnt-S, EWnt-L1, and EWnt-L2 cells were maintained in complete media before harvesting. To evaluate DIO-induced gene expression in vitro, E-Wnt, M-Wnt, MDA-MB-231, EWnt-S, and EWnt-L2 cells were serum-starved for 18 hours, and then exposed to media containing 2% control or 2% DIO mouse serum in serum-free media (SFM) for 24 hours. RNA was isolated and reverse transcribed, and quantitative RT-PCR was completed as described above.

In vitro metastatic phenotype assays

For the cell viability assay, all cell lines were seeded at a density of 5 × 103 in 96-well plates. After 24 hours, cells were serum-starved for 18 hours and then continuously exposed to 2% control or 2% DIO mouse serum in SFM for 48 hours. MTT reagent was used to assess cell viability levels as described previously (24), with absorbance read at 570 nm on a Cytation 3 Cell Imaging Multi-Mode Reader (BioTek Instruments Inc.).

To assess cell migration, all cell lines were seeded in 2-Well Culture-Inserts (ibidi) in 24-well plates with 7.5 × 104 cells per insert well. After 24 hours, cells were serum-starved for 18 hours, and then inserts removed and cells continuously exposed to 2% control or 2% DIO mouse serum in SFM for 6 hours. The EVOS FL Auto 2 Cell Imaging System (Invitrogen) was used to capture images (×10 magnification) of the wound diameter at baseline and at 6 hours.

To measure invasive capacity, 1.0 × 104 cells were seeded in SFM in each 24-well Corning Matrigel Invasion Chamber (8.0 μm; Corning), and 2% control or 2% DIO mouse serum in SFM was placed beneath the chambers. After 24 hours, the invading cells were fixed in 100% methanol for 1 minute, and then stained for 30 minutes with 0.5% crystal violet in 50% methanol. The EVOS FL Auto 2 Cell Imaging System (Invitrogen) was used to capture images (×10 magnification) of the stained cells.

Statistical analyses

Animal study data are presented as mean ± SD and in vitro data as mean ± SEM. All in vitro data shown represent the average of at least 3 independent experiments. For all statistical tests, GraphPad Prism software was used (GraphPad Software Inc.), and P < 0.05 was considered significantly different. Differences between animals or cells exposed to two different experimental conditions were analyzed using Student t test. Kaplan–Meier survival curves were plotted, and between-group difference in survival was analyzed using the log-rank (Mantel–Cox) test. Differences between cells exposed to more than two experimental conditions were analyzed using one-way ANOVA, followed by Tukey post hoc test. Data from experiments with more than one independent variable were analyzed using two-way ANOVA, followed by Tukey post hoc test.

DIO regimen promotes excess body weight, excess body fat, and systemic metabolic and inflammatory perturbations

Female ovariectomized MMTV-Wnt-1 transgenic mice fed a DIO versus control diet regimen had significantly greater body weight from week 2 through week 35 on diet (P < 0.05; Supplementary Fig. S1A) and significantly greater body weight and percent body fat at final assessment (week 42 or when tumor diameter was >1.5 cm, whichever occurred first; P < 0.01 for each comparison; Supplementary Fig. S1B and S1C). At final assessment, DIO mice, versus controls, also had significantly elevated serum insulin, leptin, resistin (each P < 0.001), gastric inhibitory peptide (P < 0.01), glucagon-like peptide 1, TNFα, and IL17A (each P < 0.05; Table 1). Serum IGF-1 (P = 0.06), plasminogen activator inhibitor-1 (P = 0.09), and IL6 (P = 0.09) were elevated in DIO mice relative to controls, but these differences did not reach statistical significance.

DIO promotes MMTV-Wnt-1 mammary tumor development

DIO mice had significantly reduced survival compared with controls (P < 0.05; Fig. 1A). One DIO mouse and 4 control mice remained alive and tumor-free at 42 weeks. Mammary tumors from DIO mice, relative to controls, expressed significantly higher levels of the proliferation marker Ki67 (P < 0.01) and endothelial cell marker CD31 (P < 0.05; Fig. 1B). Gross necropsy detected no liver or lung metastases in either diet group.

DIO modulates CSC- and EMT-related tumoral gene and protein expression

Using RNA sequencing and Qiagen's Ingenuity Pathway Analysis (IPA), we characterized differential mammary tumor gene expression in randomly selected tumors from DIO versus control mice. Genes differentially expressed (P < 0.05; Fig. 1C) were evaluated using IPA. This analysis indicated significant upregulation of protumorigenic pathways (associated genes listed parenthetically) characterized as “regulation of EMT” (Akt3, Pik3r1, Twist1, Twist2, Wnt11, Wnt9b; P = 0.003), “human embryonic stem cell pluripotency,” (Akt3, Pik3r1, Wnt11, Wnt9b; P = 0.026), and “Wnt/β-catenin signaling” (Akt3, Pik3r1, Frat1, Wnt11, Wnt9b; P = 0.038), and identified Lep as a highly significant (P < 0.001) upstream master regulator with the highest z-score (z = 2.06).

Using quantitative RT-PCR, we verified diet group–dependent differences in tumoral expression of CSC- and EMT-related genes identified by our RNA sequencing analysis (Fig. 2A) or previous studies (11) in mammary tumor transplant models (Fig. 2B). Tumors from DIO mice versus controls displayed significantly upregulated Aldh1a1 (P < 0.01), Akt3, Pik3r1, Twist1, and Twist2 (each P < 0.05; Fig. 2A) and Pou5f1 (P < 0.01), Foxc2, and Vim (each P < 0.05; Fig. 2B) as well as downregulated Cdh1 (P < 0.05; Fig. 2B). No between-group differences were detected in tumoral expression of Nanog, Notch1, Snai1, Snai2, Tgfb, and Zeb1 (data not shown). Tumor ALDH activity, a measure of CSC enrichment, was greater in DIO mice than controls (P < 0.01; Fig. 2C). In addition, tumors in DIO mice versus controls had significantly elevated protein expression of the mesenchymal marker vimentin (both P < 0.05), but not lower protein expression of the epithelial marker E-cadherin (Fig. 2D), as assessed by IHC.

We examined whether local signaling of leptin (an IPA-identified master regulator) is increased in the tumor microenvironment of DIO mice versus controls. Tumor Lepr and mammary fat pad Lep expression were significantly higher in DIO mice compared with controls (P < 0.05 for both; Fig. 2E). No between-group differences in tumor expression of Lep were found (data not shown).

Leptin signaling promotes CSC phenotype in triple-negative mammary tumor cells in vitro

The in vitro effect of leptin (100, 200, and 400 ng/mL, final concentrations) on mammary CSC enrichment was examined in two murine (E-Wnt and M-Wnt) and one human (MDA-MB-231) triple-negative mammary tumor cell lines by quantifying mammospheres formed during two sequential propagations (P1 and P2). Leptin, compared with vehicle, significantly increased mammosphere formation in E-Wnt cells at 400 ng/mL leptin in P2 (P < 0.05; Fig. 3A), in M-Wnt cells at each leptin concentration in P2 (100 ng/mL, P < 0.05; 200 and 400 ng/mL, P < 0.001; Fig. 3B), and in MDA-MB-231 cells at each leptin concentration in P1 (100 ng/mL, P < 0.01; 200 and 400 ng/mL, P < 0.001) and P2 (100 ng/mL, P < 0.001; 200 and 400 ng/mL, P < 0.0001; Fig. 3C). In each cell line, there were more mammospheres in P2 than P1 for each leptin concentration (E-Wnt: 100 and 200 ng/mL, P < 0.05; 400 ng/mL, P < 0.001; M-Wnt: 100 ng/mL, P < 0.01; 200 and 400 ng/mL, P < 0.05; MDA-MB-231: 400 ng/mL, P < 0.001; 100 and 200 ng/mL, P < 0.0001; Fig. 3A–C).

To establish the causal relationship between obesity-induced leptin signaling and CSC enrichment, we used shRNA to silence Lepr expression in E-Wnt cells. Specifically, we characterized leptin- and CSC/EMT–related gene expression and mammosphere formation capability of two stably transfected Lepr knockdown E-Wnt lines (EWnt-L1 and EWnt-L2) relative to E-Wnt cells stably transfected with a scrambled plasmid (EWnt-S) and parental E-Wnt cells (EWnt-P). Lepr expression in EWnt-L1 and EWnt-L2 was significantly reduced by 58.7% ± 13.1% and 62.3% ± 9.1%, respectively, compared with EWnt-S cells (P < 0.001 for both), whereas there was no significant difference in Lepr expression between EWnt-S and EWnt-P cells (Fig. 3D). In both EWnt-L1 and EWnt-L2 cells compared with EWnt-P and EWnt-S cells, mammosphere formation (each P < 0.001) and expression of the leptin target gene Socs3 (each P < 0.001) were significantly reduced (Fig. 3E and F). In addition, EWnt-L1 and EWnt-L2 cells, relative to EWnt-P and EWnt-S cells, had significantly reduced expression of the CSC/EMT–related genes Foxc2 (P < 0.001 for all) and Twist2 (EWnt-P and EWnt-S vs. EWnt-L1, P < 0.01 and EWnt-L2, P < 0.001), and Vim expression was significantly lower in EWnt-L2 than EWnt-S cells (P < 0.05; Fig. 3G). No significant differences between cell lines were noted in expression of Akt3, Aldh1a1, Cdh1, Cdh2, Pik3r1, Pou5f1, Sox2, Twist1, Wnt9b, and Wnt11 (data not shown).

Obesity-associated circulating factors stimulate metastatic phenotype and CSC/EMT gene signature in triple-negative mammary tumor cells

Because metastasis is closely linked with EMT and CSC enrichment, we examined the impact of obesity-associated circulating factors on in vitro measures of metastatic potential. Pooled serum isolated from female wild-type C57BL/6 mice maintained on the DIO (n = 92) or control (n = 110) diet regimens for 14 weeks was used to model an obese versus nonobese control state in vitro. Body weight in C57BL/6 mice on the DIO diet was significantly greater than controls during weeks 2 to 14 on diet (P < 0.0001 for all weeks; Supplementary Fig. S2). Serum collected from the DIO mice, compared with controls, had significantly elevated levels of several hormones, adipokines, and cytokines, including leptin (P < 0.0001; Supplementary Table S1).

Cell viability was significantly higher in E-Wnt (P < 0.01), M-Wnt, and MDA-MB-231 (both P < 0.05) cells following a 48-hour incubation in media containing serum from DIO versus control mice (Fig. 4A). DIO mouse serum also promoted significantly greater cell migration, relative to control, in E-Wnt, M-Wnt, and MDA-MB-231 cells during a 6-hour incubation (all P < 0.001; Fig. 4B). Furthermore, DIO mouse serum stimulated significantly greater invasive activity in all 3 cell lines, as a higher number of E-Wnt, M-Wnt, and MDA-MB-231 cells invaded over 24 hours through Matrigel membrane toward media containing serum from DIO versus control mice (all P < 0.001, Fig. 4C).

A 24-hour incubation of triple-negative mammary tumor cells in media containing DIO mouse serum, relative to control serum, resulted in significantly increased (or decreased where noted) expression in E-Wnt cells of the CSC- and EMT-related genes Akt3 (P < 0.001), Aldh1a1, Cdh2, Foxc2, Vim (each P < 0.01), Cdh1 (decreased), Pou5f1, Twist1, and Twist2 (each P < 0.05); in M-Wnt cells of Akt3, Twist1, Twist2 (each P < 0.001), Sox2 (P < 0.01), Cdh1 (decreased), Pou5f1, and Vim (each P < 0.05); and in MDA-MB-231 cells of AKT3, TWIST1 (both P < 0.001), TWIST2 (P < 0.01), ALDH1A1, CDH2, FOXC2, POU5F1, SOX2, and VIM (each P < 0.05; Fig. 4D). Aldh1a1 was not detectable in M-Wnt cells. No differences were observed in Wnt9b, Wnt11, or Pik3r1 expression in the 3 cell lines (data not shown).

Leptin signaling contributes to metastatic phenotype and CSC/EMT gene signature in an in vitro model of obesity and TNBC

To examine the causal role of leptin signaling in the metastatic phenotype observed following mammary tumor cell exposure to DIO mouse serum, we assessed cell viability, migration, and invasion capability in Lepr knockdown E-Wnt cell lines following exposure to DIO mouse and control serum. Relative to control, DIO mouse serum promoted significantly greater cell viability in EWnt-S cells (P < 0.05) but not EWnt-L1 or EWnt-L2 cells. Cell viability was significantly lower in EWnt-L2 versus EWnt-S cells exposed to DIO serum (P < 0.05) and did not significantly differ between cell lines exposed to control mouse serum (Fig. 5A). DIO mouse serum, relative to control serum, significantly increased migration of EWnt-P (P < 0.05) and EWnt-S (P < 0.01) cells but not EWnt-L1 or EWnt-L2 cells. Migration was significantly reduced in EWnt-L2 cells, relative to EWnt-P, EWnt-S, and EWnt-L1 cells, when cultured in DIO mouse serum (all P < 0.05), and did not differ between cell lines when cultured in control mouse serum (Fig. 5B). DIO mouse serum, relative to control, significantly increased invasion in Lepr-proficient cells (EWnt-P and EWnt-S, both P < 0.0001) and to a lesser extent Lepr knockdown cells (EWnt-L1 and EWnt-L2, both P < 0.01). When stimulated with DIO mouse serum, the number of invading cells was significantly lower for EWnt-L2 versus EWnt-P (P < 0.0001) and EWnt-S (P < 0.01) cells, whereas EWnt-L1 cells demonstrated intermediate invasion that was significantly lower than the EWnt-P cells (P < 0.05). Invasion did not significantly differ between cell lines when exposed to control mouse serum (Fig. 5C).

Finally, we assessed CSC/EMT gene expression in EWnt-S and EWnt-L2 cell lines when exposed for 24 hours to DIO mouse serum, compared with control serum (Fig. 5D). In EWnt-S cells, exposure to DIO versus control serum significantly modulated Cdh2 (P < 0.001), Cdh1, Sox2, Twist2, Vim (each P < 0.01), Akt3, and Foxc2 (both P < 0.05) expression (as seen previously in parental E-Wnt cells, Fig. 4D). In EWnt-L2 cells, knockdown of Lepr eliminated the difference between the DIO and control mouse serum conditions for Foxc2, Twist2, and Vim expression. These genes also exhibited reduced expression in EWnt-L2 cells compared with EWnt-S cells, particularly in cells exposed to DIO mouse serum (DIO mouse serum-treated EWnt-S vs. EWnt-L2 cells: Foxc2, P < 0.01; Twist2, P < 0.0001; and Vim, P < 0.001). Furthermore, differences between the serum conditions in Akt3 and Sox2 expression were attenuated with Lepr knockdown in the EWnt-L2 cells.

Although obesity is associated with a poor TNBC prognosis (3), mechanisms mediating this link are not fully understood, and the current standard of TNBC care does not differ for obese patients. Our findings suggest a leptin-dependent mechanism involving CSC enrichment and EMT, which could lead to increased tumor development and metastasis. Specifically, we report that (i) obesity decreases tumor-free survival and promotes a CSC/EMT gene signature in the MMTV-Wnt-1 transgenic mouse model of TNBC; (ii) tumors from DIO mice have higher expression of genes that promote stem cell pluripotency and EMT; and (iii) an IPA-predicted upstream regulator of these genes is leptin. Given the significant elevation of systemic leptin in DIO mice, we focused on defining the role of leptin signaling in obesity-induced CSC enrichment and EMT. Utilizing in vitro models of obesity and TNBC, we found that leptin is a key contributor to obesity's effects, specifically promoting a CSC/EMT phenotype and upregulating the expression of multiple CSC/EMT–related genes, including Foxc2, Twist2, Vim, and to a lesser extent Akt3 and Sox2 (Fig. 6). To our knowledge, this is the first study to characterize the impact of DIO on global gene expression in MMTV-Wnt-1 tumors and the first to identify Foxc2, Twist2, and Akt3 as leptin-regulated targets in tumors.

We and others have previously explored leptin's impact on CSC enrichment and/or EMT in breast and other cancers. We established that (i) the growth of orthotopically transplanted MMTV-Wnt-1 tumors is enhanced in genetically obese db/db mice, which have a compensatory elevation in leptin production due to a Lepr gene mutation; and (ii) MMTV-Wnt-1 tumor growth is reduced in the genetically obese ob/ob model, which lacks functional Lep expression (12). These findings suggested that leptin signaling may contribute to obesity's protumor effects on TNBC. In a follow-up study in nonobese C57BL/6 and nude mice, we also showed Lepr silencing in M-Wnt cells inhibits transplanted tumor growth and decreases Nanog expression, whereas LEPR silencing in human MDA-MB-231 cells decreases NANOG, POU5F1, SOX2, and VIM and increases CDH1 expression (25). Chang and colleagues (21) also found that leptin induces an EMT phenotype, mammosphere formation, and LEPR expression via STAT3 and miR-200c in preclinical models of normal breast epithelial cells (MCF12A) and luminal A breast cancer (MCF7). In addition, sorted LEPRHI primary human breast cancer cells have elevated sphere-forming capacity and CSC/EMT markers, including higher NANOG, SOX2, and POU5F1 (21). Others have found similar leptin-mediated increases in CSC/EMT markers in models of both luminal A and TNBC, suggesting that leptin's effects are likely not specific to one breast cancer subtype (26–32). Feldman and colleagues (22) further demonstrated that leptin treatment stimulates liver CSC Sox2 and Pou5f1 (but not Nanog) expression via STAT3, and Lepr expression is induced by SOX2 and OCT4 (the protein encoded by Pou5f1). Leptin treatment also upregulates EMT and/or CSC markers in ovarian (23) and pancreatic (33) cancer cell lines and promotes the self-renewal of normal human breast epithelial stem cells (34). These data suggest that leptin not only promotes CSC enrichment and EMT in multiple types of cancer but also induces CSC Lepr expression to maintain the pluripotent state.

Lepr knockdown in E-Wnt cells indicated that only Foxc2, Twist2, Vim, Akt3, and Sox2 are regulated by leptin signaling in that cell line. The lack of a leptin-mediated effect on Nanog, Pou5f1, and Cdh1 expression contrasted with others' findings (21, 22, 25) and may be related to differences between cell lines used. Zheng and colleagues (25) utilized M-Wnt and MDA-MB-231 cell lines, both of which cluster with claudin-low breast cancer subtype and have a more mesenchymal phenotype, while we knocked down Lepr in the E-Wnt cell line, which clusters genetically with basal-like breast cancer (11). We chose to utilize E-Wnt cells for our mechanistic studies because they are more representative of the majority of the tumor cells that comprise spontaneous basal-like triple-negative MMTV-Wnt-1 tumors (11). In addition, given that they are more epithelial-like and less aggressive than the M-Wnt cell line, we were particularly interested in examining whether leptin and other obesity-associated factors could induce a more mesenchymal phenotype and genotype in these cells.

Our finding that leptin signaling regulates mammary cancer cell Foxc2 and Akt3 expression is novel. In addition, although previous studies have demonstrated that Twist expression is increased with leptin treatment in several breast cancer cell lines (27, 28), we found that leptin regulates the expression of Twist2, but not Twist1. The function and regulation of the third Akt isoform remains poorly characterized in comparison with Akt1/2, but Akt3 signaling is known to stimulate CSC enrichment in breast tumor models (35, 36). Others have demonstrated that Akt3 amplification is more frequent in TNBC versus estrogen receptor–positive tumors and is negatively correlated with recurrence-free survival in TNBC (37, 38). The transcription factors FOXC2 and TWIST2 are both closely linked to stem cell–like properties as well as EMT and metastatic capability in breast and other cancers (39–44). In addition, FOXC2 mediates resistance to multiple forms of chemotherapy, including paclitaxel, epirubicin, 5-fluorouracil, and cisplatin in breast cancer cells, in correlation with an EMT phenotype (41, 45). Furthermore, FOXC2 is enriched in claudin-low breast tumors in comparison with other subtypes, which links this gene with the most aggressive form of breast cancer (32). We previously demonstrated that obesity promotes greater Foxc2 expression in orthotopically transplanted E-Wnt and M-Wnt cells (11), but we are unaware of any studies linking obesity with greater tumoral Twist2 or Akt3 expression. Here, we report that obesity upregulates Foxc2, Twist2, and Akt3 expression in spontaneous MMTV-Wnt-1 tumors in association with greater local and systemic leptin signaling (Fig. 2A and B). Moreover, Lepr knockdown ablates E-Wnt expression of Foxc2, Twist2, and Vim, while partially attenuating Akt3 and Sox2 expression, in response to DIO mouse serum (Fig. 5D). Consequently, Foxc2, Twist2, Vim, Akt3, and Sox2 may be mediators of obesity's pro-CSC/EMT effects that are all induced at least in part by leptin.

This study used global gene expression data from a spontaneous model of TNBC to drive the generation of hypotheses regarding potential mechanisms mediating the reduced survival in obese MMTV-Wnt-1 mice. In addition, the use of mouse sera to model the obese versus nonobese tumor microenvironment in vitro allowed us to examine the combined impact of multiple obesity-associated systemic factors as well as the relative contribution of leptin. However, this in vitro approach to assessing the effects of obesity in mammosphere assays was precluded because serum inhibits sphere formation, so we were limited to treating with leptin alone. Our in vitro model of obesity using DIO mouse serum was also limited by its inability to account for intercellular cross-talk in the tumor microenvironment. Given that TNBC is a highly heterogeneous disease that includes multiple distinct molecular subtypes (46, 47), our use of the basal-like E-Wnt cell line alone for the mechanistic experiments is an additional study limitation. We did not mechanistically link leptin signaling to an EMT phenotype and CSC/EMT–related gene expression in the two claudin-low cell lines used (M-Wnt and MDA-MB-231), hindering our ability to extend our conclusions beyond the basal-like subtype. In addition, we are not able to conclude that the observed effects are limited to TNBC, as we did not assess leptin's effects on models representing other breast cancer subtypes. Finally, although this study provides a broad characterization of leptin's impact on CSC/EMT–related phenotypes and genes in TNBC cells, we did not pursue definition of the exact molecular pathways mediating the observed downstream effects. We can therefore only hypothesize that these effects are likely mediated through one or both of the two major pathways implicated by others in leptin's pro-CSC/EMT effects: JAK2/STAT3 and PI3K/Akt (Fig. 6; refs. 21–23, 25–28, 30).

Given the continuing increase in global obesity rates, it is critical that we improve our understanding of mechanisms mediating the obesity–breast cancer link. Our findings indicate that increased leptin signaling is causally linked to obesity-associated TNBC development by promoting CSC enrichment and EMT. Future studies will further analyze the leptin targets we identified (Foxc2, Twist2, and Akt3) to define their role in obesity-induced CSC enrichment, EMT, and chemotherapy resistance in TNBC. In addition, the contribution of other adipokines and cytokines to this obesity-associated CSC/EMT genotype and phenotype will be examined, as our in vitro analyses indicate that leptin signaling does not fully account for all of the observed effects of DIO on TNBC. Greater understanding of the signals impacting CSC and EMT may provide new targets and intervention strategies for decreasing TNBC burden in obese women.

L.A. deGraffenried has provided expert testimony for Par Pharmaceutical. No potential conflicts of interest were disclosed by the other authors.

Conception and design: L.W. Bowers, E.L. Rossi, L.A. deGraffenried, S.D. Hursting

Development of methodology: L.W. Bowers, S.B. McDonell, S.D. Hursting

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): L.W. Bowers, E.L. Rossi, S.B. McDonell, S.A. Khatib, J.E. Albright, S.D. Hursting

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): L.W. Bowers, S.B. McDonell, S.S. Doerstling, S.A. Khatib, C.G. Lineberger, X. Tang

Writing, review, and/or revision of the manuscript: L.W. Bowers, S.B. McDonell, S.S. Doerstling, S.A. Khatib, X. Tang, L.A. deGraffenried, S.D. Hursting

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): S.B. McDonell, S.S. Doerstling

Study supervision: S.D. Hursting

This study was supported by grants from the Breast Cancer Research Foundation (to S.D. Hursting) and the NIH (R35 CA197627, to S.D. Hursting). L.W. Bowers and E.L. Rossi were supported by fellowships from the NCI (R25 CA057726), and S.B. McDonell was supported by a fellowship from the NIH (T32 LM12420-2).

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Flegal
KM
,
Carroll
MD
,
Kit
BK
,
Ogden
CL
. 
Prevalence of obesity and trends in the distribution of body mass index among US adults, 1999-2010
.
JAMA
2012
;
307
:
491
7
.
2.
Protani
M
,
Coory
M
,
Martin
JH
. 
Effect of obesity on survival of women with breast cancer: systematic review and meta-analysis
.
Breast Cancer Res Treat
2010
;
123
:
627
35
.
3.
Bao
PP
,
Cai
H
,
Peng
P
,
Gu
K
,
Su
Y
,
Shu
XO
, et al
Body mass index and weight change in relation to triple-negative breast cancer survival
.
Cancer Causes Control
2016
;
27
:
229
36
.
4.
Widschwendter
P
,
Friedl
TW
,
Schwentner
L
,
DeGregorio
N
,
Jaeger
B
,
Schramm
A
, et al
The influence of obesity on survival in early, high-risk breast cancer: results from the randomized SUCCESS A trial
.
Breast Cancer Res
2015
;
17
:
129
.
5.
Litton
JK
,
Gonzalez-Angulo
AM
,
Warneke
CL
,
Buzdar
AU
,
Kau
SW
,
Bondy
M
, et al
Relationship between obesity and pathologic response to neoadjuvant chemotherapy among women with operable breast cancer
.
J Clin Oncol
2008
;
26
:
4072
7
.
6.
Chen
S
,
Chen
CM
,
Zhou
Y
,
Zhou
RJ
,
Yu
KD
,
Shao
ZM
. 
Obesity or overweight is associated with worse pathological response to neoadjuvant chemotherapy among Chinese women with breast cancer
.
PLoS One
2012
;
7
:
e41380
.
7.
Iwase
T
,
Sangai
T
,
Nagashima
T
,
Sakakibara
M
,
Sakakibara
J
,
Hayama
S
, et al
Impact of body fat distribution on neoadjuvant chemotherapy outcomes in advanced breast cancer patients
.
Cancer Med
2016
;
5
:
41
8
.
8.
Pajares
B
,
Pollan
M
,
Martin
M
,
Mackey
JR
,
Lluch
A
,
Gavila
J
, et al
Obesity and survival in operable breast cancer patients treated with adjuvant anthracyclines and taxanes according to pathological subtypes: a pooled analysis
.
Breast Cancer Res
2013
;
15
:
R105
.
9.
Del Fabbro
E
,
Parsons
H
,
Warneke
CL
,
Pulivarthi
K
,
Litton
JK
,
Dev
R
, et al
The relationship between body composition and response to neoadjuvant chemotherapy in women with operable breast cancer
.
Oncologist
2012
;
17
:
1240
5
.
10.
Wei
W
,
Lewis
MT
. 
Identifying and targeting tumor-initiating cells in the treatment of breast cancer
.
Endocr Relat Cancer
2015
;
22
:
R135
55
.
11.
Dunlap
SM
,
Chiao
LJ
,
Nogueira
L
,
Usary
J
,
Perou
CM
,
Varticovski
L
, et al
Dietary energy balance modulates epithelial-to-mesenchymal transition and tumor progression in murine claudin-low and basal-like mammary tumor models
.
Cancer Prev Res
2012
;
5
:
930
42
.
12.
Zheng
Q
,
Dunlap
SM
,
Zhu
J
,
Downs-Kelly
E
,
Rich
J
,
Hursting
SD
, et al
Leptin deficiency suppresses MMTV-Wnt-1 mammary tumor growth in obese mice and abrogates tumor initiating cell survival
.
Endocr Relat Cancer
2011
;
18
:
491
503
.
13.
Goodwin
PJ
,
Ennis
M
,
Pritchard
KI
,
Trudeau
ME
,
Koo
J
,
Taylor
SK
, et al
Insulin- and obesity-related variables in early-stage breast cancer: correlations and time course of prognostic associations
.
J Clin Oncol
2012
;
30
:
164
71
.
14.
Ishikawa
M
,
Kitayama
J
,
Nagawa
H
. 
Enhanced expression of leptin and leptin receptor (OB-R) in human breast cancer
.
Clin Cancer Res
2004
;
10
:
4325
31
.
15.
Garofalo
C
,
Koda
M
,
Cascio
S
,
Sulkowska
M
,
Kanczuga-Koda
L
,
Golaszewska
J
, et al
Increased expression of leptin and the leptin receptor as a marker of breast cancer progression: possible role of obesity-related stimuli
.
Clin Cancer Res
2006
;
12
:
1447
53
.
16.
Miyoshi
Y
,
Funahashi
T
,
Tanaka
S
,
Taguchi
T
,
Tamaki
Y
,
Shimomura
I
, et al
High expression of leptin receptor mRNA in breast cancer tissue predicts poor prognosis for patients with high, but not low, serum leptin levels
.
Int J Cancer
2006
;
118
:
1414
9
.
17.
Xu
J
,
Prosperi
JR
,
Choudhury
N
,
Olopade
OI
,
Goss
KH
. 
beta-Catenin is required for the tumorigenic behavior of triple-negative breast cancer cells
.
PLoS One
2015
;
10
:
e0117097
.
18.
O'Flanagan
CH
,
Rossi
EL
,
McDonell
SB
,
Chen
X
,
Tsai
YH
,
Parker
JS
, et al
Metabolic reprogramming underlies metastatic potential in an obesity-responsive murine model of metastatic triple negative breast cancer
.
NPJ Breast Cancer
2017
;
3
:
26
.
19.
Rossi
EL
,
Dunlap
SM
,
Bowers
LW
,
Khatib
SA
,
Doerstling
SS
,
Smith
LA
, et al
Energy balance modulation impacts epigenetic reprogramming, ERalpha and ERbeta expression, and mammary tumor development in MMTV-neu transgenic mice
.
Cancer Res
2017
;
77
:
2500
11
.
20.
Korkaya
H
,
Paulson
A
,
Iovino
F
,
Wicha
MS
. 
HER2 regulates the mammary stem/progenitor cell population driving tumorigenesis and invasion
.
Oncogene
2008
;
27
:
6120
30
.
21.
Chang
CC
,
Wu
MJ
,
Yang
JY
,
Camarillo
IG
,
Chang
CJ
. 
Leptin-STAT3-G9a signaling promotes obesity-mediated breast cancer progression
.
Cancer Res
2015
;
75
:
2375
86
.
22.
Feldman
DE
,
Chen
C
,
Punj
V
,
Tsukamoto
H
,
Machida
K
. 
Pluripotency factor-mediated expression of the leptin receptor (OB-R) links obesity to oncogenesis through tumor-initiating stem cells
.
Proc Natl Acad Sci U S A
2012
;
109
:
829
34
.
23.
Kato
S
,
Abarzua-Catalan
L
,
Trigo
C
,
Delpiano
A
,
Sanhueza
C
,
Garcia
K
, et al
Leptin stimulates migration and invasion and maintains cancer stem-like properties in ovarian cancer cells: an explanation for poor outcomes in obese women
.
Oncotarget
2015
;
6
:
21100
19
.
24.
Bowers
LW
,
Cavazos
DA
,
Maximo
IX
,
Brenner
AJ
,
Hursting
SD
,
deGraffenried
LA
. 
Obesity enhances nongenomic estrogen receptor crosstalk with the PI3K/Akt and MAPK pathways to promote in vitro measures of breast cancer progression
.
Breast Cancer Res
2013
;
15
:
R59
.
25.
Zheng
Q
,
Banaszak
L
,
Fracci
S
,
Basali
D
,
Dunlap
SM
,
Hursting
SD
, et al
Leptin receptor maintains cancer stem-like properties in triple negative breast cancer cells
.
Endocr Relat Cancer
2013
;
20
:
797
808
.
26.
Wang
L
,
Tang
C
,
Cao
H
,
Li
K
,
Pang
X
,
Zhong
L
, et al
Activation of IL-8 via PI3K/Akt-dependent pathway is involved in leptin-mediated epithelial-mesenchymal transition in human breast cancer cells
.
Cancer Biol Ther
2015
;
16
:
1220
30
.
27.
Yan
D
,
Avtanski
D
,
Saxena
NK
,
Sharma
D
. 
Leptin-induced epithelial-mesenchymal transition in breast cancer cells requires beta-catenin activation via Akt/GSK3- and MTA1/Wnt1 protein-dependent pathways
.
J Biol Chem
2012
;
287
:
8598
612
.
28.
Wei
L
,
Li
K
,
Pang
X
,
Guo
B
,
Su
M
,
Huang
Y
, et al
Leptin promotes epithelial-mesenchymal transition of breast cancer via the upregulation of pyruvate kinase M2
.
J Exp Clin Cancer Res
2016
;
35
:
166
.
29.
Thiagarajan
PS
,
Zheng
Q
,
Bhagrath
M
,
Mulkearns-Hubert
EE
,
Myers
MG
,
Lathia
JD
, et al
STAT3 activation by leptin receptor is essential for TNBC stem cell maintenance
.
Endocr Relat Cancer
2017
;
24
:
415
26
.
30.
Mishra
AK
,
Parish
CR
,
Wong
ML
,
Licinio
J
,
Blackburn
AC
. 
Leptin signals via TGFB1 to promote metastatic potential and stemness in breast cancer
.
PLoS One
2017
;
12
:
e0178454
.
31.
Giordano
C
,
Chemi
F
,
Panza
S
,
Barone
I
,
Bonofiglio
D
,
Lanzino
M
, et al
Leptin as a mediator of tumor-stromal interactions promotes breast cancer stem cell activity
.
Oncotarget
2016
;
7
:
1262
75
.
32.
Strong
AL
,
Ohlstein
JF
,
Biagas
BA
,
Rhodes
LV
,
Pei
DT
,
Tucker
HA
, et al
Leptin produced by obese adipose stromal/stem cells enhances proliferation and metastasis of estrogen receptor positive breast cancers
.
Breast Cancer Res
2015
;
17
:
112
.
33.
Harbuzariu
A
,
Rampoldi
A
,
Daley-Brown
DS
,
Candelaria
P
,
Harmon
TL
,
Lipsey
CC
, et al
Leptin-Notch signaling axis is involved in pancreatic cancer progression
.
Oncotarget
2017
;
8
:
7740
52
.
34.
Esper
RM
,
Dame
M
,
McClintock
S
,
Holt
PR
,
Dannenberg
AJ
,
Wicha
MS
, et al
Leptin and adiponectin modulate the self-renewal of normal human breast epithelial stem cells
.
Cancer Prev Res
2015
;
8
:
1174
83
.
35.
Liu
S
,
Patel
SH
,
Ginestier
C
,
Ibarra
I
,
Martin-Trevino
R
,
Bai
S
, et al
MicroRNA93 regulates proliferation and differentiation of normal and malignant breast stem cells
.
PLoS Genet
2012
;
8
:
e1002751
.
36.
Chin
YR
,
Yoshida
T
,
Marusyk
A
,
Beck
AH
,
Polyak
K
,
Toker
A
. 
Targeting Akt3 signaling in triple-negative breast cancer
.
Cancer Res
2014
;
74
:
964
73
.
37.
O'Hurley
G
,
Daly
E
,
O'Grady
A
,
Cummins
R
,
Quinn
C
,
Flanagan
L
, et al
Investigation of molecular alterations of AKT-3 in triple-negative breast cancer
.
Histopathology
2014
;
64
:
660
70
.
38.
Nakatani
K
,
Thompson
DA
,
Barthel
A
,
Sakaue
H
,
Liu
W
,
Weigel
RJ
, et al
Up-regulation of Akt3 in estrogen receptor-deficient breast cancers and androgen-independent prostate cancer lines
.
J Biol Chem
1999
;
274
:
21528
32
.
39.
Mani
SA
,
Yang
J
,
Brooks
M
,
Schwaninger
G
,
Zhou
A
,
Miura
N
, et al
Mesenchyme Forkhead 1 (FOXC2) plays a key role in metastasis and is associated with aggressive basal-like breast cancers
.
Proc Natl Acad Sci U S A
2007
;
104
:
10069
74
.
40.
Paranjape
AN
,
Soundararajan
R
,
Werden
SJ
,
Joseph
R
,
Taube
JH
,
Liu
H
, et al
Inhibition of FOXC2 restores epithelial phenotype and drug sensitivity in prostate cancer cells with stem-cell properties
.
Oncogene
2016
;
35
:
5963
76
.
41.
Hollier
BG
,
Tinnirello
AA
,
Werden
SJ
,
Evans
KW
,
Taube
JH
,
Sarkar
TR
, et al
FOXC2 expression links epithelial-mesenchymal transition and stem cell properties in breast cancer
.
Cancer Res
2013
;
73
:
1981
92
.
42.
Fang
X
,
Cai
Y
,
Liu
J
,
Wang
Z
,
Wu
Q
,
Zhang
Z
, et al
Twist2 contributes to breast cancer progression by promoting an epithelial-mesenchymal transition and cancer stem-like cell self-renewal
.
Oncogene
2011
;
30
:
4707
20
.
43.
Liu
AY
,
Cai
Y
,
Mao
Y
,
Lin
Y
,
Zheng
H
,
Wu
T
Twist2 promotes self-renewal of liver cancer stem-like cells by regulating CD24
.
Carcinogenesis
2014
;
35
:
537
45
.
44.
Wang
T
,
Li
Y
,
Wang
W
,
Tuerhanjiang
A
,
Wu
Z
,
Yang
R
, et al
Twist2, the key Twist isoform related to prognosis, promotes invasion of cervical cancer by inducing epithelial-mesenchymal transition and blocking senescence
.
Hum Pathol
2014
;
45
:
1839
46
.
45.
Cai
J
,
Tian
AX
,
Wang
QS
,
Kong
PZ
,
Du
X
,
Li
XQ
, et al
FOXF2 suppresses the FOXC2-mediated epithelial-mesenchymal transition and multidrug resistance of basal-like breast cancer
.
Cancer Lett
2015
;
367
:
129
37
.
46.
Bianchini
G
,
Balko
JM
,
Mayer
IA
,
Sanders
ME
,
Gianni
L
. 
Triple-negative breast cancer: challenges and opportunities of a heterogeneous disease
.
Nat Rev Clin Oncol
2016
;
13
:
674
90
.
47.
Abramson
VG
,
Mayer
IA
. 
Molecular heterogeneity of triple negative breast cancer
.
Curr Breast Cancer Rep
2014
;
6
:
154
8
.