Learn more: PMC Disclaimer | PMC Copyright Notice
Leading malignant cells initiate collective epithelial cell invasion in a three-dimensional heterotypic tumor spheroid model
Associated Data
Abstract
Solid tumors consist of genetically and phenotypically diverse subpopulations of cancer cells with unique capacities for growth, differentiation, and invasion. While the molecular and microenvironmental bases for heterogeneity are increasingly appreciated, the outcomes of such intratumor heterogeneity, particularly in the context of tumor invasion and metastasis, remain poorly understood. To study heterotypic cell-cell interactions and elucidate the biological consequences of intratumor heterogeneity, we developed a tissue-engineered multicellular spheroid (MCS) co-culture model that recapitulates the cellular diversity and fully three-dimensional cell-cell and cell-matrix interactions that characterize human carcinomas. We found that “invasion-competent” malignant cells induced the collective invasion of otherwise “invasion-incompetent” epithelial cells, and that these two cell types consistently exhibited distinct leader and follower roles during invasion. Analysis of extracellular matrix microarchitecture revealed that malignant cell invasion was accompanied by extensive extracellular matrix remodeling including matrix alignment and proteolytic track-making. Inhibition of cell contractility- and proteolysis-mediated matrix reorganization prevented leader-follower behavior and malignant cell-induced epithelial cell invasion. These results indicate that heterogeneous subpopulations within a tumor may possess specialized roles during tumor progression and suggest that complex interactions among the various subpopulations of cancer cells within a tumor may regulate critical aspects of tumor biology and affect clinical outcome.
Introduction
The tumor microenvironment contains a complex mixture of normal epithelial cells, highly heterogeneous malignant cells, stromal cells, and extracellular matrix (ECM) whose interactions are critical to cancer progression [1, 2]. Despite advances in in vitro models that have refined our understanding of how the cellular and noncellular components of the tumor microenvironment cooperate to promote or suppress disease, many tissue-engineered tumor models fail to accurately recapitulate the intratumor heterogeneity and three-dimensional tissue architecture that characterize human carcinomas [3–5]. Notably, biophysical and biochemical interactions among dissimilar cell types within the tumor microenvironment enable many key features of developing cancer including abnormal cell growth [6], enhanced angiogenesis [7], and tissue invasion and metastasis [8–10]. However, the importance of interactions among dissimilar cells within the tumor compartment remains to be determined. To better understand how a tumor’s constituents cooperate to regulate critical phases of tumor progression, these processes should be studied using well-controlled in vitro tumor models that permit three-dimensional, in vivo-like interactions among the diverse cells and ECM that comprise the microenvironment.
Since cancer cells continually interact with other cells that are resident within or recruited to the tumor microenvironment, several experimental models have coupled malignant cells with cells such as stromal fibroblasts, endothelial cells, and immune cells to study the underlying mechanisms and clinically-relevant outcomes of heterotypic cell-cell interactions. For example, Gaggioli et al observed that stromal fibroblasts produce cell contractility- and proteolysis-dependent tracks within the ECM that are sufficient to induce squamous cell carcinoma invasion [8]. Additionally, activation of stromal fibroblasts by gastric cancer cells induces stromal expression of pro-angiogenic factors resulting in increased tumor vascularization [7]. Importantly, the cellular diversity of the tumor microenvironment is not limited to these distinct cell types, as cells of epithelial lineage within carcinomas exist along a spectrum ranging from normal epithelial cells to highly variable malignant cells [11]. In addition to the significant heritable diversity that neoplastic cells can exhibit due to the acquisition of genetic and epigenetic alterations during oncogenesis [12], cancer cells exhibit phenotypic plasticity, whereby a cell’s genetic program is integrated with biochemical and biophysical cues from the extracellular microenvironment to regulate cellular phenotype [13–15]. The resulting cellular diversity within the epithelial compartment itself gives rise to substantial morphological, physiological, and behavioral intratumor heterogeneity whose functional consequences are poorly understood.
While it is clear that human carcinomas contain subpopulations of cancer cells with unique capacities for growth, differentiation, tissue invasion, and secondary tumor initiation [11, 12], the significance of interactions among cells in these various subgroups are unknown. Since local tissue invasion is both the first discernible step of metastasis and the basis for histopathological diagnosis of metastatic cancer, the invasive potential of tumor cells is of particular clinical relevance [16, 17]. During invasion, cancer cells physically dissociate from the primary tumor mass, bypass the epithelial basement membrane to move from the epithelial compartment into the stromal compartment, and negotiate the three-dimensional interstitial microenvironment on their way to a secondary site [18]. Migration through the three-dimensional stroma requires specialized cellular mechanisms including proteolysis- and cell contractility-based matrix remodeling to overcome steric hindrance imposed by the extracellular matrix [19–22]. Tissue invasion is dependent upon cells’ invasive fitness as well as cell-cell and cell-matrix interactions [13, 23]. For example, 3D extracellular matrix microarchitecture [24], alignment [25], and mechanics [26] regulate 3D migration efficiency [27]. Biophysical cell-cell and cell-matrix interactions including cell-cell adhesions and proteolytic matrix patterning provide additional control over the nature of invasion [13, 21, 22, 28]. Together, these findings have established important roles for extracellular control of invasion, but there remains a need to investigate cancer invasion in a more physiologically relevant heterogeneous cancer model where only a subset of cancer cells within a heterogeneous tumor may be “invasion-competent”.
Tissue engineering approaches to cancer research have emerged as valuable intermediates between traditional 2D in vitro cell culture techniques, which often fail to accurately represent the microenvironmental complexity of tumors, and in vivo cancer models, whose extensive complexity can limit experimental control and confound findings [5, 29, 30]. Tissue-engineered platforms such as multicellular spheroids (MCS) incorporate cells and ECM in a three-dimensional physiological context, and thus, are able to effectively recapitulate tumor architecture and cancer cell function, which are coupled through regulation of cell-cell and cell-ECM interactions [31–33]. Thus, tumor spheroids have been widely used to investigate tumorigenesis [34], cellular mechanisms of cancer invasion [20, 35], and anticancer drug efficacy [36]. Importantly, MCS can be made to contain multiple cell types in co-culture and can be dynamically and quantitatively analyzed with confocal microscopy, making them an ideal experimental model with which to explore the functional importance of intratumor heterogeneity.
In this study, we investigated cancer invasion in the framework of intratumor heterogeneity using a tissue-engineered co-culture tumor model in which two dissimilar cell types derived from the epithelial compartment were incorporated into heterotypic multicellular spheroids. Co-culture MCS composed of the “invasion-competent” breast adenocarcinoma cell line MDA-MB-231 and the “invasion-incompetent” breast epithelial cell line MCF-10A were embedded within 3D collagen matrices that permitted in vivo-like cell-cell and cell-matrix interactions and, thus, supported physiological in vitro tumor progression. Using this co-culture platform, we found that invasive malignant cells induced and led collective invasion of otherwise non-invasive epithelial cells, and that this leader-follower co-invasive behavior was dependent upon cell contractility- and proteolysis-based ECM remodeling by leading malignant cells.
Materials and Methods
Cell culture and reagents
Malignant MDA-MB-231 breast adenocarcinoma cells (HTB-26; ATCC, Rockville, MD) were maintained in DMEM (Invitrogen, Carlsbad, CA) supplemented with 10% fetal bovine serum (Atlanta Biologicals, Norcross, GA) and 1% penicillin-streptomycin (Invitrogen). MDA-MB-231/GFP cells (AKR-201; Cell Biolabs, San Diego, CA) were maintained in complete MDA-MB-231 media supplemented with 0.1mM MEM Non-Essential Amino Acids (Invitrogen). MCF-10A mammary epithelial cells (CRL-10317; ATCC) and MCF-10CA1a malignant mammary epithelial cells (Barbara Ann Karmanos Cancer Institute, Detroit, MI) were maintained in DMEM:F12 (Invitrogen) supplemented with 5% horse serum (Invitrogen), 0.5 μg ml−1 hydrocortisone (Sigma-Aldrich, St. Louis, MO), 20 ng ml−1 hEGF (Invitrogen), 10 μg ml−1 insulin (Sigma-Aldrich), 100 ng ml−1 cholera toxin (Sigma-Aldrich), and 1% penicillin-streptomycin. PC-3 malignant prostate adenocarcinoma cells (CRL-1435; ATCC) were maintained in Ham’s F-12K Medium (ATCC) supplemented with 10% fetal bovine serum and 1% penicillin-streptomycin. PrEC primary human prostate epithelial cells (CC-2555; Lonza, Walkersville, MD) were maintained in PrEGM prostate epithelial cell growth medium (Lonza) supplemented with SingleQuots (Lonza) according to the manufacturer’s recommended protocol. Mammary spheroid formation media consisted of complete MCF-10A media supplemented with 0.25% methylcellulose (H4110; StemCell Technologies, Vancouver, BC). Mammary spheroid growth media consisted of complete MCF-10A media supplemented with 0.25% methylcellulose and 1% Matrigel (BD; San Jose, CA). Prostate spheroid formation media consisted of complete PrEC media supplemented with 0.25% methylcellulose. Prostate spheroid growth media consisted of complete PrEC media supplemented with 0.25% methylcellulose and 1% Matrigel. All cell culture was maintained at 37°C and 5% CO2.
At least 2 h prior to spheroid generation, epithelial cells (MCF-10A, PrEC) and malignant cells (MDA-MB-231, MCF-10CA1a, PC-3) were labeled with 10 μM CellTracker Orange CMRA (Invitrogen) and 10 μM Cell Tracer Green CMFDA SE (Invitrogen), respectively. GM6001 (EMD Millipore, Billerica, MA) and Y27632 (Sigma-Aldrich) were used at 20 μM and 10 μM, respectively.
Multicellular spheroid generation
To generate multicellular spheroids, fluorescently-labeled cells were trypsinized as for cell passage and resuspended in spheroid formation media at 2.5×104 cells ml−1. To prepare epithelial and malignant monoculture spheroids, 200 μl of the appropriate cell suspension (5×103 cells) was seeded into each well of non-adhesive round-bottom 96-well plates (Corning; Tewksbury, MA). For epithelial/malignant co-culture spheroids, volumetric mixtures of each cell suspension were seeded into each well. Plates were centrifuged at 1000 rpm for 5 min and then incubated at 37°C and 5% CO2 on an orbital shaker at ~60 rpm for 2 h. Spheroid formation media was then replaced with spheroid growth media to promote spheroid compaction.
3D spheroid invasion assay
After 48 h of compaction, spheroids were embedded in 1.5 mg ml−1 and 6.0 mg ml−1 type I collagen gels. Collagen gels were prepared as previously described [24]. Briefly, type I collagen was acid-extracted from rat tail tendons (Pel-Freez Biologicals, Rogers, AR), purified via centrifugation and lyophilization, and reconstituted at 10 mg ml−1 in 0.1% acetic acid. Stock collagen solution was diluted to either 1.5 mg ml−1 or 6.0 mg ml−1 by gently mixing with ice-cold DMEM, and the solution was neutralized to pH 7.0 with 1N NaOH. Spheroids were removed from culture plates and individually embedded within 500 μl collagen gels in glass-bottom 24-well plates (MatTek, Ashland, MA). After 30 minutes of gel polymerization at 37°C, gels were overlaid with 500 μl of complete MCF-10A or PrEC media with or without pharmacological inhibitors.
Microscopy and image analysis
Real-time spheroid compaction and invasion were measured with time-lapse microscopy using a Zeiss AxioObserver Z1 inverted phase contrast microscope equipped with a Hamamatsu ORCA-ER camera. For compaction analysis, spheroids were imaged at 2-h intervals for 24 h in a temperature-, humidity-, and CO2-controlled microscope incubation chamber using AxioVision software (v. 4.8, Carl Zeiss MicroImaging GmbH, Jena, Germany). Spheroid cross-sectional area was measured using ImageJ (v. 1.43u, National Institute of Health, Bethesda, MD) and normalized to initial cross-sectional area for n ≥ 10 spheroids per condition from at least three independent experiments. Simultaneous fluorescence and phase time-lapse microscopy (imaged at 20-min intervals for 24 h) were used to observe invasion of fluorescently-labeled MDA-MB-231 and MCF-10A cells from co-culture mammary spheroids into 1.5 mg ml−1 collagen matrices.Intra-spheroid cell organization and spheroid invasion were imaged using a Zeiss LSM700 confocal microscope operated by ZEN software (v. 2010, Carl Zeiss). Low-power (5× and 10×) fluorescent images of spheroids presented herein are maximum intensity z-projections of confocal image stacks acquired from the bottom surface of the spheroid to the spheroid center. Spheroid invasion was calculated by measuring the projected spheroid area immediately after collagen embedding (A0) and the projected spheroid area following culture within collagen matrix for 48 h (AInv). Invasive Index was defined as (AInv/A0) − 1 for n ≥ 10 spheroids per condition from at least three independent experiments. Maximal invasion distance was determined by measuring the radial distance (in μm) from the spheroid edge to the invasive cells furthest from the spheroid for n ≥ 10 spheroids per condition (four cells measured per spheroid) from at least three independent experiments. To further investigate invasion from co-culture mammary spheroids, the number of epithelial cell-containing invasive strands per spheroid was quantified for n ≥ 10 co-culture spheroids after 48 h of invasion.
For E-cadherin immunofluorescence, mammary spheroids composed of MDA-MB-231/GFP and unlabeled MCF-10A cells were embedded in 1.5 mg ml−1 collagen and either fixed immediately or allowed to invade for 48 h before fixation. Briefly, samples were fixed with 3.7% (v/v) formaldehyde in PBS, permeabilized with 1% (v/v) Triton (JT Baker, Phillipsburg, NJ) in PBS, and blocked with 0.5% (v/v) Tween (JT Baker) with 3% (w/v) bovine serum albumin (Sigma-Aldrich) in PBS. An anti-E-Cadherin primary antibody (sc-7870; Santa Cruz Biotechnology Inc., Santa Cruz, CA) was used at (1:50) in 1% (w/v) BSA in PBS and an Alexa Fluor 568 conjugated secondary antibody (Santa Cruz) was used at (1:100) in 1% (w/v) BSA. Nuclei were stained with 4,6-diamidino-2-phenylindole (DAPI; Sigma-Aldrich).
Cell localization/immunofluorescence and collagen fiber organization before and during spheroid invasion were assessed with high-power (40×) fluorescence and reflectance confocal microscopy, respectively. Confocal reflectance images of collagen fibers are 1-μm thick confocal slices acquired as previously described near the spheroid “equator” [24]. The ImageJ plugin OrientationJ was used as previously described to measure collagen fiber orientation from confocal reflectance images [37]. Grayscale confocal reflectance images of collagen fibers were analyzed and colorized with OrientationJ such that pixel hue corresponded to the angle of local fiber orientation, which could range from −90° to +90° relative to horizontal. To quantify collagen fiber organization, collagen fiber angle (θFiber) relative to the tangent of the original spheroid surface (θSph) was determined for each pixel in the given region of interest using the algorithm |θFiber - θSph|. Using this analysis, the minimum angular difference was 0°, which corresponded to a tangentially-aligned fiber; the maximum angular difference was 90°, which corresponded to a radially-aligned fiber. Results are graphed as frequency distributions of the angular difference, binned in 10-degree increments, between collagen fibers and the spheroid surface for all pixels within the indicated regions of interest.
Statistical analysis
Data was compared by one-way analysis of variance (ANOVA) with post-hoc Tukey HSD tests using JMP Pro software (v. 9.0.2, SAS Institute, Cary, NC). All numeric data is presented as mean ± SE; statistical significance was considered with p < 0.01.
Results
Multicellular spheroid generation and characterization
To experimentally recapitulate intratumor heterogeneity and study heterogeneity-associated phenomena, a multicellular spheroid model was developed. Fluorescently-labeled epithelial cells (labeled red) and malignant cells (labeled green) were seeded in round-bottom 96-well plates (Fig. 1a; 0 h) and allowed to spontaneously coalesce and compact over 24 h to form multicellular spheroids (Fig. 1a; insets). Fluorescence confocal microscopy and transmitted light microscopy were used to probe cellular distribution and spheroid morphology, respectively. MCF-10A epithelial monoculture spheroids were densely-compacted, while MDA-MB-231 malignant monoculture spheroids were more loosely aggregated. Co-culture spheroids (1:1 mixture of MCF-10A cells and MDA-MB-231 cells) exhibited intermediate compaction (Fig. 1a; 24 h). Tracking of spheroid area over 24 h showed that spheroid compaction, and thus, final spheroid size, was dependent on cellular composition (Fig. 1b). Purely epithelial spheroids compacted rapidly and had final cross-sectional areas of 16.9 ± 0.8% of the initial cell pellet. Co-culture spheroids compacted more slowly and had final cross-sectional areas of 19.9 ± 0.5% of the initial cell pellet. Purely malignant spheroids exhibited slower compaction than spheroids containing epithelial cells and had cross-sectional areas of 28.4 ± 2.1 % of the initial cell pellet.
Multicellular spheroid compaction. (a) Representative confocal fluorescence images (5×) show cell compaction during formation of MCF-10A mammary epithelial cell monoculture (red), MDA-MB-231 malignant breast cell monoculture (green), and MCF10A/MDA-MB-231 co-culture spheroids. Phase images illustrate spheroid morphologies after 24 hours of culture. Insets provide cartoon representation of cells (red) compacting within non-adhesive round-bottom wells to form spheroids. (b) Quantification of spheroid cross-sectional area during compaction; data are presented as spheroid cross-sectional area normalized to initial cross-sectional area for n ≥ 11 spheroids from at least three independent experiments. Scale bar = 500 μm
While epithelial and malignant monoculture spheroids contained homogeneous cellular distributions, epithelial/malignant co-culture spheroids contained clusters of each cell type. Confocal microscopy showed that these cell clusters were organized into distinct regions of the spheroid (Fig. 2a). Consistently, a core of malignant cells was encased in a shell of epithelial cells (Fig. 2a; 100–125 μm), which was surrounded by an outer layer of malignant cells (Fig. 2a; all slices). Time-lapse fluorescence microscopy demonstrated that malignant cells peripheral to the epithelial shell were highly motile at the conclusion of compaction (Fig. 2b). E-cadherin was localized to MCF-10A cell membranes (Fig. 2c; white arrowheads), but absent from fluorescently-labeled MDA-MB-231 cells (Fig. 2c; black arrowheads) in co-culture spheroids.
Cell organization within heterotypic co-culture spheroids. (a) Segmented 25-μm thick z-projections (10×) of co-culture tumor spheroid consisting of fluorescently-labeled MCF-10A (red) and MDA-MB-231 (green) cells after 24 h of compaction highlight localization of two cell types through partial thickness of spheroid. (b) Time-lapse imaging at peak spheroid compaction (22–25 h) shows compaction-independent intraspheroid movement of malignant MDA-MB-231 cells. Arrowheads track the movement of a highly-motile cluster of malignant cells that have been excluded to the spheroid periphery. (c) Maximum intensity z-projection of co-culture tumor spheroid consisting of fluorescently-labeled MDA-MB-231 and unlabeled MCF-10A cells. Black arrowheads highlight fluorescently-labeled MDA-MB-231 cells and white arrowheads highlight E-cadherin-positive cell membranes of MCF-10A cells. Scale bars = 50 μm
Three-dimensional multicellular spheroid invasion
To examine invasion from multicellular spheroids in a 3D physiological context, fluorescently-labeled spheroids were fully embedded within collagen gels and imaged immediately and after 48 h. During this culture period, MCF-10A epithelial cells in monoculture spheroids remained non-invasive (Fig. 3a), with these spheroids only undergoing some non-invasive expansive growth, as indicated by a small non-zero Invasive Index of 0.22 ± 0.04 (Fig. 3b) and no detectable Maximal Invasion Distance (Fig. 3c). In contrast, MDA-MB-231 malignant cells at the periphery of both epithelial/malignant co-culture and malignant monoculture spheroids extended exploratory cell protrusions and began to invade the matrix within 3 h after collagen embedding. After 48 h, co-culture and MDA-MB-231 malignant spheroids showed extensive collagen invasion (Fig. 3a), resulting in statistically significant Invasive Indices of 2.3 ± 0.1 and 3.6 ± 0.2, respectively (Fig. 3b). Cells invading from these spheroids exhibited statistically equivalent Maximal Invasion Distances of 276 ± 13 μm and 296 ± 14 μm, respectively (Fig. 3c). MDA-MB-231 malignant cells (labeled green) invaded readily in a disorganized manner from both co-culture and malignant spheroids. Notably, when cultured with malignant cells in co-culture spheroids, normally non-invasive MCF-10A epithelial cells (labeled red) became invasive (Fig. 3a; arrowheads). Examination of cells invading from co-culture spheroids revealed that invading epithelial cells (Fig. 4a; black arrowheads) were consistently organized in invasive strands that were led by one or more leading malignant cells (Fig. 4a; white arrowheads). Another malignant mammary cell type, MCF-10CA1a, was similarly able to induce invasion of MCF-10A epithelial cells (Fig. 4b). Further, PC-3 malignant prostate adenocarcinoma cells induced invasion of otherwise non-invasive (data not shown) PrEC primary prostate epithelial cells (Fig. 4c). In each of these systems, malignant cells (MDA-MB-231, MCF10CA1a, or PC-3) were positioned at the tip of invasive strands (Fig. 4a–c, white arrowheads) and epithelial cells (MCF-10A or PrEC) were trailing behind (Fig. 4a–c, black arrowheads).
Three-dimensional multicellular tumor spheroid invasion. Fluorescently-labeled spheroids were embedded within 1.5 mg ml−1 collagen gels and imaged immediately after embedding and at 48 h. (a) Representative confocal images (5×) of embedded MCF-10A monoculture, MCF-10A/MDA-MB-231 co-culture, and MDA-MB-231 monoculture spheroids at 0 h and 48 h. Arrowheads indicate strands of epithelial MCF-10A cells invading from co-culture spheroids. (b) Quantification of Invasive Index (a.u.) at 48 h; data are presented as spheroid cross-sectional area normalized to initial cross-sectional area at t = 0 h. (c) Quantification of maximal invasion distance (μm) at 48 h; data are presented as the radial distance from the spheroid edge to the invasive cells furthest from the spheroid. Data in (b) and (c) are from n ≥ 10 spheroids from at least three independent experiments. Scale bar = 500 μm
Leader-follower behavior in invasive strands. Confocal images (10×) of fluorescently-labeled cells invading into 1.5 mg ml−1 collagen matrix from co-culture spheroids demonstrate leading malignant cells (green; white arrowheads) and following epithelial cells (red; black arrowheads). Co-culture combinations include: (a) MDA-MB-231 malignant breast adenocarcinoma cells (green) and MCF-10A mammary epithelial cells (red), (b) MCF10CA1a malignant mammary epithelial cells (green) and MCF-10A mammary epithelial cells (red), and (c) PC-3 malignant prostate adenocarcinoma cells (green) and PrEC primary prostate epithelial cells (red). Arrow indicates direction of invasion away from spheroids, which are in the bottom left corner of each image. Representative images from three independent experiments. Scale bar = 100 μm
Extracellular matrix remodeling during spheroid invasion
To more closely examine the differential invasion observed from multicellular spheroids, confocal fluorescence and reflectance microscopy were used to probe cell localization and extracellular matrix organization at the spheroid-matrix interface, respectively. Immediately following collagen embedding, all spheroids were surrounded by uniform collagen matrices containing randomly organized fibers (Fig. 5a; 0 h). After 48 h, collagen fiber reorganization and alignment were observed in all conditions (Fig. 5a; 48 h, white arrowheads) as compared to naïve matrices at 0 h. The non-invasive expansive growth exhibited by MCF-10A epithelial spheroids was accompanied by reorganization of collagen fibers tangential to the spheroid surface such that fibers were wrapped around the spheroid (Fig. 5a; 48 h). Quantitative analysis of fiber alignment using OrientationJ confirmed that the fibers in the indicated region of interest were preferentially aligned parallel to the epithelial spheroid surface, as shown by a small difference between θFiber and θSph (Fig. 5b; MCF-10A). Confocal imaging revealed that cellular invasion from epithelial/malignant co-culture and MDA-MB-231 malignant monoculture spheroids was associated with remodeling of collagen fibers perpendicular to the spheroid surface (Fig. 5a; 48 h). Quantification of ECM alignment in the proximity of invading cells with OrientationJ verified that collagen fibers were oriented orthogonal to the spheroid surface, as indicated by relatively large differences between θFiber and θSph (Fig. 5b; Co-culture and MDA-MB-231). Consistently, leading malignant cells were elongated and radially-oriented. Malignant and epithelial following cells tended to show rounded morphologies, with the latter forming cohesive invasive strands. Notably, invading malignant MDA-MB-231 leader cells left cell-scale microtracks in the matrix behind them, as indicated by regions of high collagen reflectance surrounding regions devoid of extracellular matrix (Fig. 5a; 48 h, black arrowheads in cross-section inset). Time-lapse imaging revealed the time-course and dynamics of invasion from co-culture spheroids and showed that leading malignant MDA-MB-231 cells persisted at the tips of invasive strands (Fig. 6a; white arrowheads). Importantly, co-invasion was maintained even when transient gaps (Fig. 6a; black arrowheads) formed between leading MDA-MB-231 cells and following MCF-10A cells. Together, time-lapse imaging and confocal microscopy (Fig. 6b) showed that trailing epithelial cells exhibited two general induced-invasion strategies: (i) MCF-10A cells fill matrix tunnels (Fig. 6b; black arrowheads) left by invading MDA-MB-231 cells; and (ii) MCF-10A cells invade directly behind and adjacent to leading MDA-MB-231 cells. E-cadherin was consistently detected between following MCF-10A epithelial cells (Fig. 6b,c; white arrowheads), but not in leading MDA-MB-231 cells or between adjacent MDA-MB-231 cells (Fig. 6c,d).
Matrix reorganization associated with malignant cell invasion. (a) Maximum intensity z-projections (40×) of multicellular tumor spheroids composed of epithelial MCF-10A cells (red) and malignant MDA-MB-231 cells (green), embedded within 1.5 mg ml−1 collagen matrix (white fibers imaged by confocal reflectance microscopy) and imaged immediately after embedding (0 h, left panels) and following culture (48 h, middle panels). Regions of interest at 48 h are indicated in middle panels by white boxes and are magnified 3× in right panels. Confocal reflectance images in right panels were colorized with OrientationJ as described in Materials and Methods. White arrowheads highlight areas of collagen fiber alignment. Inset in center panel is cross-sectional view along the invasive strand. Black arrowheads indicate boundaries of matrix microtrack generated by invading malignant leader cells. (b) Quantitative analysis of collagen fiber alignment relative to the spheroid surface. Briefly, local collagen fiber orientations (θFiber) were determined with OrientationJ and compared to the angle given by the original spheroid surface (θSph), resulting in alignment distributions representative of all fibers in the regions of interest. Scale bars = 50 μm
Intercellular interactions during co-invasion. (a) Time-lapse fluorescent and phase microscopy (10×) of a co-culture spheroid invading into 1.5 mg ml−1 collagen matrix. White arrowheads highlight MDA-MB-231 cell (green) leading an invasive strand of MCF-10A cells (red). Black arrowheads identify transient gap between leading MDA-MB-231 cell and following MCF-10A cells. Time shown indicates hours after embedding. (b) Single confocal section (40×) of multicellular spheroid composed of MCF-10A mammary epithelial cells and fluorescently-labeled MDA-MB-231 malignant breast cells (green) after 48 h of invasion into collagen matrix (white fibers). E-cadherin (red) and nuclei (DAPI; blue) were detected by confocal fluorescence microscopy. Inset highlights E-cadherin-positive cell-cell adhesions between MCF-10A cells (white arrowheads). Black arrowheads indicate boundaries of matrix microtracks generated by invading malignant leader cells. Two MCF-10A induced-invasion strategies are demonstrated: (i) MCF-10A cells fill matrix tunnels left by invading MDA-MB-231 cells; (ii) MCF-10A cells remain adjacent to leading MDA-MB-231 cells, but with no detectable E-cadherin-based adhesions between the two cell types. (c, d) Maximum intensity z-projections of invasive cellular strands from different regions of a co-culture spheroid demonstrate (c) an E-cadherin-positive (white arrowheads) cohort of MCF-10A cells following leading MDA-MB-231 cells and (d) lack of E-cadherin-positive cell-cell adhesions between adjacent singly invading MDA-MB-231 cells. Insets highlight E-cadherin distribution. Arrow indicates direction of invasion away from spheroids for panels b-d. Scale bars = 50 μm
Since fiber reorganization and proteolytic track-making were associated with cell invasion and malignant MDA-MB-231 cells consistently led invasive strands of epithelial MCF-10A cells during invasion from co-culture spheroids, we hypothesized that extracellular matrix remodeling by leading malignant cells was required for 3D invasion from multicellular spheroids. To test this hypothesis, co-culture spheroids were embedded within 1.5 mg ml−1 collagen matrices and treated with GM6001 (20 μM) to inhibit MMP activity or Y-27632 (10 μM) to inhibit ROCK-based cell contractility. Confocal reflectance imaging around invading cells at the spheroid-matrix interface showed that collagen fiber alignment and matrix microtracks, which were prominent in control conditions, were absent under MMP and ROCK inhibition (Fig. 7a; black arrowheads). Quantification of ECM alignment indicated that treatment with GM6001 or Y27632 attenuated the fiber alignment observed in control conditions. Briefly, fibers within control ECM tended to be oriented perpendicular to the spheroid surface (large difference between θFiber and θSph) while fibers within GM6001- and Y27632-treated matrices showed reduced preferential orientation (Fig. 7c). MDA-MB-231 invasion past the original spheroid boundary (Fig. 8a; dashed white lines) persisted under both MMP and ROCK inhibition (Fig. 8a). However, the extent of invasion was decreased: MMP and ROCK inhibition significantly reduced Invasive Index from 2.3 ± 0.1 to 1.5 ± 0.2 and 1.4 ± 0.2, respectively (Fig. 8c), and significantly reduced Maximal Invasion Distance from 276 ± 13 μm to 217 ± 19 μm and 146 ± 15 μm, respectively (Fig. 8d). Similar trends of reduced invasion under MMP and ROCK inhibition were observed in MDA-MB-231 monoculture spheroid studies (Supplementary Fig. S1). Notably, co-culture spheroids exhibited 9.8 ± 1.5 MCF-10A+ invasive strands per spheroid (Fig. 8a; white arrowheads), but no such epithelial cell-positive invasive strands were observed under MMP or ROCK inhibition (Fig. 8a,e).
Effects of proteolysis and cell contractility on matrix reorganization during invasion. Confocal reflectance images (40×) of cells invading from co-culture spheroids into (a) 1.5 mg ml−1 or (b) 6.0 mg ml−1 collagen matrices under control conditions or treated with inhibitors of matrix metalloproteinases (GM6001; 20 μM) and ROCK (Y-27632; 10 μM). Images are single slices aligned parallel to the plane of invasion taken at the spheroid-matrix interface. Black arrowheads indicate boundaries of matrix microtracks generated by leading malignant cells. Insets are high-magnification (20μm × 20μm) images of 1.5 and 6.0 mg ml−1 collagen matrices to highlight fiber and pore microarchitecture. Quantification of collagen fiber orientation relative to the spheroid surface at the spheroid-matrix interface within (c) 1.5 mg ml−1 or (d) 6.0 mg ml−1 collagen matrices under control conditions (purple bars), MMP inhibition with GM6001 (orange bars), or ROCK inhibition with Y27632 (blue bars). Scale bar = 50μm
Role of proteolysis and cell contractility in induced collective epithelial cell invasion. Maximum intensity z-projections (10×) of co-culture multicellular spheroids embedded within (a) 1.5 mg ml−1 or (b) 6.0 mg ml−1 collagen gels and treated with GM6001 (20 μM) or Y-27632 (10 μM). Arrowheads indicate epithelial cell-containing invasive strands. Dashed lines indicate original spheroid boundaries. Quantification of (c) Invasive Index, (d) Maximal Invasion Distance, and (e) the incidence of MCF-10A-positive invasive strands under control conditions (black bars), GM6001 treatment (gray bars), and Y27632 treatment (white bars). Dashed gray line in (c) denotes Invasive Index of non-invasive epithelial monoculture spheroids. MCF-10A-positive invasive strands were not detected under GM6001 or Y27632 treatment (n.d.). * indicates statistical significance compared to control 1.5 mg ml−1 condition; ** indicates statistical significance compared to control 6.0 mg ml−1 condition. Scale bar = 50 μm
To investigate the role of matrix microarchitecture in regulating invasion, co-culture spheroids were embedded within high-density 6.0 mg ml −1 collagen matrices. These matrices have a more confined microarchitecture consisting of a tight meshwork of short collagen fibers and small pores (Fig. 7b; inset) as compared to low-density 1.5 mg ml −1 collagen matrices, which contain larger fibers and pores (Fig. 7a; inset). As in low-density gels, matrix remodeling was prominent in control conditions (Fig. 7b; black arrowheads), but was reduced in cultures subjected to MMP or ROCK inhibition, which showed no matrix microtracks and reduced collagen fiber alignment (Fig. 7b,d). Invasion into high-density matrices, as indicated by Invasive Index and Maximal Invasion Distance, was significantly reduced by 55% and 40%, respectively, compared to invasion into low-density matrices under control conditions (Fig. 8c,d); however, matrix density did not change the frequency of MCF-10A+ invasive strands under control conditions (Fig. 8e). Similar to as in low-density matrices, inhibition of MMP or ROCK activity mitigated invasion into high-density matrices (Fig. 8b–d). MMP and ROCK inhibition significantly reduced the Invasive Indices to 0.28 ± 0.04 and 0.35 ± 0.03, respectively, which were close to the Invasive Index exhibited by non-invasive MCF-10A epithelial monoculture spheroids (Fig. 8c; dashed line). Max invasion distance was reduced from 165 ± 10 to 58 ± 7 and 76 ± 6 by MMP and ROCK inhibition, respectively (Fig. 8d). These values are consistent with the observation that malignant cells extended thin protrusions from the spheroid periphery into the matrix, but did not translate completely past the spheroid boundary into the matrix (Fig. 8b; dashed white lines). As in low-density matrix, inhibition of MMP and ROCK activity prevented the development of MCF-10A+ co-invasive strands in 6.0 mg ml−1 matrices (Fig. 8b,e).
Discussion
Solid tumors, which initially arise by clonal expansion of transformed cells, contain diverse subpopulations that result from genetic instability, microenvironmental selective pressures, and the reinforcement of enabling phenotypes [12, 38, 39]. Such heritable diversity, when coupled with the phenotypic plasticity that characterizes cancer cells [15], results in substantial intratumor heterogeneity that can obscure the study, diagnosis, and treatment of cancer [4, 40]. To study the functional consequences of such cellular variety in a controlled in vitro system, we have developed a tissue-engineered model of intratumor heterogeneity using two dissimilar mammary cell types in co-culture multicellular spheroids. Using this platform, we have established a new model for collective cancer cell invasion in which invasive malignant cells can initiate the invasion of otherwise non-invasive cells within the tumor through cell contractility- and proteolysis-dependent matrix remodeling. To our knowledge, this is the first study that links the distinct leader and follower roles inherent in collective cancer cell invasion to differential invasiveness and intratumor heterogeneity. Thus, this study provides new insight into interactions between cell subpopulations within heterogeneous tumors.
Development and characterization of a fundamental in vitro model of intratumor heterogeneity
To generate a well-defined in vitro model of intratumor heterogeneity, co-culture multicellular spheroids were made by mixing equal numbers of fluorescently-labeled cells from two dissimilar mammary cell lines: the malignant breast adenocarcinoma cell line MDA-MB-231 and the breast epithelial cell line MCF-10A. These two cell lines were used throughout this study to represent two discrete subpopulations of cells within a tumor and were chosen because of their distinct invasive/metastatic potentials as well as their considerable historical and scientific significance within the field of breast cancer research. While MCF-10A cells exhibit a more basal phenotype than epithelial cells from normal breast tissue [41], these cells are non-tumorigenic in mice and, because they possess many features of normal mammary epithelial cells, they are extensively used as a model of normal epithelial cells in vitro [42, 43]. Malignant and epithelial monoculture spheroids were used throughout the study to demonstrate these cells’ behavior in monoculture as a reference. To expand upon the co-culture spheroid system, an alternate malignant mammary epithelial cell line, MCF-10CA1a, was used in co-culture with MCF-10A epithelial cells. The MCF-10CA1a cell line was derived from a metastatic subpopulation of Ha-ras-overexpressing MCF-10A cells [44]. Further, PC-3 malignant prostate adenocarcinoma cells and PrEC primary prostate epithelial cells were used to investigate intratumor heterogeneity in other tissue systems.
Using the MDA-MB-231 and MCF-10A cell lines, characterization of MCS compaction showed that the rate and extent of spheroid aggregation was dependent upon cell composition. Epithelial MCF-10A spheroids showed the fastest and most extensive compaction while malignant MDA-MB-231 spheroids showed the slowest and least extensive compaction; as expected, co-culture MCF-10A/MDA-MB-231 spheroids exhibited intermediate compaction. These differences in compaction were likely due to differences in the cell-cell and cell-matrix adhesions that participated in aggregation since epithelial cell aggregation is primarily mediated by E-cadherin and malignant cell aggregation is directed by β1 integrin-ECM interactions [45, 46]. Indeed, E-cadherin was only detected on the cell membranes of epithelial MCF-10A cells, and was absent from malignant MDA-MB-231 cells within co-culture spheroids. Such disparity in intercellular adhesiveness is the basis for the differential adhesion hypothesis and, in part, may underlie the cell sorting that was observed in heterotypic co-culture spheroids [47]. Consistently, intermixed cells segregated to form distinct regions of each cell type: a core of malignant cells was encased within a densely-packed shell of epithelial cells, which was surrounded by an outer rim of highly motile malignant cells. The presence of such mobile malignant cells at the spheroid periphery is supported by previous findings [48] and suggests that these cells were actively remodeling their adhesions to reach the physical configuration that maximized adhesion [47]. Consistent with this hypothesis, peripheral malignant cells eventually settled onto and adhered to the spheroid surface. While the specific mechanisms regulating cellular organization within co-culture spheroids remain to be determined, the time course of spheroid formation suggests that differential adhesion and cell sorting, as opposed to proliferation and apoptosis, are responsible for the rapid establishment of cellular organization [49]. An inclusive characterization of the mechanisms governing, and outcomes of, differential cell adhesion, sorting, and survival within heterotypic spheroids was not the aim of this work, and thus, a more comprehensive and rigorous analysis of cell physiology within co-culture spheroids remains an open challenge. Notably, co-culture spheroids can serve as a physiologically-relevant model system within which to study how microenvironmental factors such as cell-matrix interactions, heterotypic cell-cell interactions, and spatially-defined nutrient and growth signals differentially influence cell behaviors relevant to cancer progression [50].
Physiological 3D invasion from an in vitro model of intratumor heterogeneity
Intratumor heterogeneity is thought to have considerable biological and clinical implications due to vast differences in growth, invasive, and differentiation potential among a tumor’s cells [4, 11]. During tissue invasion, which is one of the earliest rate-limiting steps of metastasis, malignant cells must spatially and temporally coordinate proteolysis with changes in adhesion, contractility, and cell shape to physically disseminate from the primary tumor [16, 18]. Thus, due to the molecular specialization that is required for efficient invasion, we and others have hypothesized that cells within the tumor may possess different capacities for invasion and metastasis based on their phenotypic abilities [17, 51–53]. To test this hypothesis experimentally, malignant and epithelial cell lines were used in co-culture multicellular spheroids to represent cell subpopulations within a tumor having distinct invasive capacities.
When compacted spheroids were embedded into three-dimensional collagen matrices that mimic the in vivo tumor stroma [54], MCF-10A cells in monoculture epithelial spheroids were unable to invade the matrix and instead exhibited some non-invasive expansive growth. While extracellular cues including matrix stiffness [31], cellular density [55], and growth factors [56] can promote an invasive mesenchymal phenotype in MCF-10A cells, cells in monoculture MCF-10A spheroids maintained a non-invasive epithelial phenotype, and thus, these cells represented an “invasion incompetent” subpopulation. Conversely, “invasion-competent” MDA-MB-231 cells in monoculture form malignant spheroids that readily invaded into the surrounding ECM by first extending exploratory cell projections into the matrix and subsequently translating the cell body away from the spheroid. These opposing behaviors, which were expected based on previous work with these cell types in 3D culture [21, 57], validated the physiological relevance of our 3D culture system. When heterotypic co-culture spheroids were embedded into 3D collagen matrix, they initially showed a pattern of invasion similar to malignant monoculture spheroids, with malignant cells promptly invading into the ECM. After 48 h of culture, the extent of invasion as indicated by measurement of the Invasive Index – the normalized area of invasion – was dependent upon the amount of malignant cells in each culture. Assessment of Maximal Invasion Distance revealed that malignant cells in co-culture spheroids are equally as invasive as malignant cells in monoculture malignant spheroids, indicating that epithelial cells did not restrain malignant cells’ invasive capacity in co-culture.
Interestingly, when co-cultured with malignant MDA-MB-231 cells, otherwise non-invasive epithelial MCF-10A cells were triggered to invade into the matrix in cohesive strands following one or more leading malignant cells. This pattern of “co-invasion”, with distinct leader cells and cohesive strands or masses of invading follower cells, is indicative of collective migration, which occurs during development, regeneration, and epithelial cancer invasion [8, 21, 28, 58–60]. An alternate malignant mammary epithelial cell line, MCF-10CA1a, was similarly able to induce and lead invasion of MCF-10A cells and PC-3 malignant prostate cells led invasive strands of normally non-invasive PrEC primary prostate epithelial cells. While there were clear differences in cellular phenotypes and the nature of co-invasion among the cell lines used, the emergence of leader-follower collective invasion and consistent presence of malignant cells at the tips of invasive strands indicate that this phenomenon is not limited to MDA-MB-231/MCF-10A co-cultures. Further, the incidence of leader-follower co-invasion in prostate co-cultures suggests that this behavior extends to other tissue systems and that non-invasive primary cells can be induced to invade by malignant leader cells in vitro. Notably, a defining feature of collective cell migration is the preservation of cell-cell adhesions, and thus, an epithelial phenotype, in migrating cells [13, 61]. Here, we found that epithelial follower cells, but not malignant leader cells, retained such cell-cell cohesion when invading from co-culture spheroids. E-cadherin was maintained at cell-cell junctions between adjacent MCF-10A cells both within the spheroid mass and within invasive strands, indicating retention of an epithelial phenotype in co-invading MCF-10A cells. Thus, using our co-culture system, we demonstrate for the first time that the leader and follower roles characteristic of collective cancer invasion can be assumed by two distinct cellular subpopulations within a tumor.
Collective cell migration through physiological three-dimensional matrices requires that matrices contain adequately large pores through which strands or clusters of cells can pass [13]. However, in the present study as well as in many tissues throughout the body [54], cell diameter can exceed ECM pore diameter, thus providing a barrier to 3D migration that requires cell-mediated extracellular matrix remodeling to permit invasion [62]. The concurrence and interdependence of 3D cell invasion and extracellular matrix remodeling has been extensively documented both in vitro and in vivo [8, 20–22, 25, 63, 64]. Notably, the Weiss group showed that cancer cells use the membrane-anchored proteinase MT1-MMP to cleave collagen and migrate through cross-linked collagen gels in vitro and interstitial ECM in vivo [22]. In defining the spatiotemporal dynamics of pericellular proteolysis during 3D invasion, Wolf et al showed that malignant cells migrating through 3D matrices can proteolytically generate tunnel-like cell-scale microtracks that facilitate further cellular infiltration and proteolysis, eventually leading to the formation of macrotracks and the emergence of collective cell invasion [21]. Using a co-culture tumor model, Gaggioli et al showed that cell-contractility and proteolysis-dependent matrix remodeling by stromal fibroblasts can induce the invasion of squamous cell carcinoma cells and that these cancer cells retain an epithelial phenotype as they collectively advance through matrix tunnels left by fibroblasts [8].
In the current study, confocal reflectance microscopy was used to assess the organization of collagen fibers around spheroids embedded within physiological fibrillar matrices. Collagen fibers around invading malignant and co-culture spheroids were preferentially aligned orthogonal to the spheroid surface while fibers surrounding non-invasive epithelial spheroids were aligned tangential to the spheroid surface. Since naïve collagen matrices and matrix far from spheroids consisted of networks of isotropically oriented fibers, these results indicate that local ECM remodeling at the spheroid-matrix interface occurred in all conditions. The circumferential alignment of fibers around non-invasive MCF-10A epithelial spheroids suggests that the growing spheroid applied outward pressure to compress the matrix [65]. Similar matrix alignment was observed in work by Provenzano et al that identified tangentially organized ECM around non-invasive mammary tumors in vivo [25]. The radial matrix organization associated with invasion in malignant monoculture and epithelial/malignant co-culture spheroids was due to active cell-mediated ECM remodeling by cells at the spheroid-matrix interface. Since cells in 3D physiological matrices exert traction forces on the ECM that can result in long-range physical changes in the surrounding matrix [66] and fibers up to several hundred microns from invading cells showed radial realignment in our model, we concluded that malignant cells caused cell contractility-mediated matrix reorganization in our system. Previous work has shown that such Rho-ROCK-based cell contractility is required for collagen fiber alignment at the tumor-stroma interface [67]. This alignment promotes contact guidance and invasion in vivo [25], and thus, may serve a similar role in our model. Further, invading malignant cells also produced cell-scale microtracks through the matrix as reported elsewhere [21, 22], and in co-culture spheroids, cohorts of epithelial cells infiltrated these microtracks. Such microtracks have been shown to enable non-proteolytic cell infiltration [35, 68]. Notably, E-cadherin was only detected between adjacent following MCF-10A epithelial cells and transient gaps between leading and following cells opened and closed during co-invasion. These results suggest that E-cadherin may be required for cohesive MCF-10A migration but that matrix tunnels, rather than cell-cell adhesions between leading and following cells, enable the induced invasive phenotype. Since we cannot rule out the possibility that malignant and epithelial cells interact using other adhesion molecules [45, 46, 49], future work should further define the mechanisms of leader-follower collective invasion in this system.
The use of cell traction forces and pericellular proteolysis to generate paths of least resistance is a general property of leader cells during three-dimensional collective migration in several systems [21, 60, 69, 70]. To determine if matrix remodeling by leading malignant cells in our co-culture system was required for the development of epithelial cell-positive co-invasive strands, cultures were treated with inhibitors of MMPs (GM6001) and ROCK-mediated cell contractility (Y27632). As expected, these inhibitors prevented the formation of matrix microtracks and quantitatively reduced matrix alignment at the spheroid-matrix interface, validating the efficacy of inhibition. While MMP and ROCK inhibition quantitatively reduced the efficiency of malignant cell invasion in both monoculture MDA-MB-231 and co-culture spheroids, presumably because the repertoire of 3D migration strategies available to cells was limited [71], these treatments completely prevented induced epithelial co-invasion. Together with our matrix remodeling observations, these findings support the hypothesis that matrix remodeling by leading invasive cancer cells was necessary for the induction of collective invasion by otherwise non-invasive cells.
Previously, we showed that 3D migration of MDA-MB-231 malignant cells is sensitive to matrix microarchitecture [24], so we used high-density 6.0 mg ml−1 collagen matrices to investigate the relationship between matrix structure, matrix remodeling, and invasion. Consistent with previous work showing that matrix density negatively correlates with 3D invasion efficiency [21, 24], here, we found that invasion from spheroids was less efficient in high-density matrices. Notably, while MDA-MB-231 cells were able to invade low-density matrix under MMP and ROCK inhibition, malignant cells at the periphery of spheroids embedded in high-density collagen matrix only extended thin protrusions into the surrounding matrix when treated with MMP and ROCK inhibitors, suggesting that high-density collagen matrices, but not low-density matrices, provided a barrier to MMP- and ROCK-independent invasion. The finding that MDA-MB-231 cells within low-density matrices retained their invasive phenotype under MMP and ROCK inhibition is supported by previous work with these cells [67], and suggests that malignant cells were sufficiently deformable to squeeze through the pores in low-density collagen matrices [72]. In support of this hypothesis, we observed that low-density collagen matrices were heterogeneous and contained a wide distribution of pore sizes, including sufficiently large pores that could potentially serve as “pro-invasive” microdomains through which malignant cells could migrate without matrix remodeling. While the absolute requirement for extracellular proteolysis during cell invasion through physiological 3D matrices remains a topic of debate, it is evident that this requirement is dependent upon local physical properties of the matrix (e.g., pore size and the presence of intermolecular crosslinks) that determine the degree to which the matrix serves as a barrier to invasion [20, 22, 73].
ROCK-mediated cell contractility has been shown to serve many functions in three-dimensional cell migration including driving cell deformation for amoeboid movement and the generation of cell tension and strengthening of adhesions for matrix remodeling [67, 74–77]. Interestingly, Sahai and Marshall showed that rounded amoeboid cell motility, but not elongated protrusive movement, has a significant requirement for Rho-ROCK activity [74]. Consistent with this, we found that malignant cells, which adopted an elongated mesenchymal motile phenotype, were still able to invade within permissive low-density matrices under ROCK inhibition. Thus, as previously shown [67], these results suggest that ROCK-mediated contractility was primarily responsible for matrix alignment and remodeling in our system, which acted to enhance, but was not required for, malignant cell invasion. Interestingly, we found that inhibition of force-mediated matrix remodeling prevented the formation of matrix microtracks by malignant cells, and thus, blocked induced MCF-10A invasion in co-cultures. These findings, which are supported by previous work by Gaggioli et al with matrix-remodeling fibroblasts [8], suggest that coordinated proteolytic- and cell force-mediated matrix remodeling were critical for the formation of invasion-inducing microtracks. Indeed, regardless of matrix density, MCF-10A cells were only induced to co-invade when malignant cells were able to both invade into and remodel the surrounding matrix. This co-invasion occurred at the same frequency in low- and high-density matrices, and was similarly associated with extracellular matrix remodeling, suggesting that the ability of leading MDA-MB-231 cells to generate matrix tunnels under control conditions was not sensitive to the matrix densities used in this study.
Our conclusion that proteolytic- and cell force-mediated matrix priming by leading invasive cells promoted the invasion of otherwise non-invasive cells is consistent with previous work in other experimental systems [8, 78]. While the use of global inhibitors of MMPs and ROCK prevents us from definitively ruling out the possibilities that (i) the epithelial follower cells themselves required ROCK and/or MMP activity to invade or that (ii) global MMP and ROCK inhibition indirectly altered off-target aspects of cell physiology [79], our data and others’ suggest that the limiting factor for induced collective epithelial invasion in this model was the efficient generation of matrix tunnels by leading malignant cells. Nonetheless, future work should seek to more selectively down-regulate cell contractility and proteolytic machinery in a cell-type and isoform-specific manner to systematically define the molecular mechanisms required for leader-follower collective invasion and limit off-target effects.
We propose the following model for malignant cell-induced epithelial cell invasion in our novel heterotypic spheroid system. First, Rho-ROCK-based cell contractility enables ECM alignment by malignant cells that extend into the matrix at the spheroid periphery. This matrix priming facilitates efficient 3D migration of malignant cells into the surrounding matrix, where they can couple their migration with pericellular proteolysis to generate cell-scale microtracks through the ECM. Finally, invasion-incompetent epithelial cells infiltrate matrix tunnels as cohesive strands or clusters of cells, resulting in collective invasion. In this paradigm, ROCK and MMP activity are required for distinct but integrated matrix remodeling functions by leading malignant cells that together lead to efficient generation of matrix tunnels and thus, induced epithelial cell invasion. Our results demonstrate a previously unappreciated instance of intratumor cellular cooperation [4] whereby invasion-competent cancer cells could provide ready-made matrix tunnels to be used by other cells that otherwise lack the ability to invade and negotiate the stromal ECM. In support of this concept, Hanahan and Weinberg recently proposed that, since heterogeneous subpopulations of cells within and around the tumor are likely to be significantly involved in all stages of cancer progression, it is unlikely that progression through the invasion-metastasis cascade is a cell-autonomous process [11]. Further, clear genetic differences have been detected among cells from separate but histologically similar regions of single metastatic tumors [12, 40], suggesting that invasive carcinomas can possess tremendous molecular and phenotypic cellular diversity that remains poorly understood. While the cell types used herein were chosen to represent two tumor subpopulations with distinct invasion capabilities, these results may also give an indication of how malignant cells could affect the remaining non-transformed epithelium within the tumor. Considering the critical role of homotypic epithelial cell-cell interactions in maintaining normal epithelial tissue homeostasis [31, 32], it is conceivable that the heterotypic cell-cell interactions characteristic of heterogeneous tumor cell populations could play a role in promoting the dysfunctional cell behavior and tissue structure observed in invasive carcinomas [80]. Ultimately, the relevance of this in vitro model to the metastatic process could be validated in vivo to determine if a subpopulation of invasive cells can similarly enable the local invasion, dissemination, and metastasis of otherwise non-invasive cells. Further characterization of heterotypic cell-cell interactions and their functional outcomes within well-defined, tissue-engineered tumor heterogeneity models should provide additional insight into the biological and clinical consequences of intratumor heterogeneity.
Supplementary Material
10585_2013_9565_MOESM1_ESM
Fig. S1:Effects of MMP and ROCK inhibition on co-culture and malignant spheroid invasion. (a) Maximum intensity z-projections (10×) of malignant MDA-MB-231 multicellular spheroids embedded within 1.5 mg ml−1 collagen gels and treated with inhibitors of matrix metalloproteinases (GM6001; 20 μM) and ROCK (Y-27632; 10 μM) for 48 h. Dashed lines indicate original spheroid boundaries. Quantification of (b) Invasive Index and (c) Maximal Invasion Distance show similar attenuation of invasion by both metrics in co-culture and monoculture malignant spheroids. Scale bar = 50 μm
Acknowledgments
This work was supported in part by the Cornell Center on the Microenvironment & Metastasis through Award Number U54CA143876 from the National Cancer Institute (NCI) to CAR and a National Science Foundation Graduate Research Fellowship to SPC.
Footnotes
Conflict of Interest
The authors declare that they have no conflict of interest.