Abstract
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
Introduction
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.
Materials and Methods
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.
Results
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.
Discussion
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.
Disclosure of Potential Conflicts of Interest
L.A. deGraffenried has provided expert testimony for Par Pharmaceutical. No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
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
Acknowledgments
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.