We have developed a novel protocol for studying heterocellular population dynamics in response to perturbations. This manuscript describes an imaging-based platform that produces quantitative datasets for simultaneous characterization of multiple cellular phenotypes of heterocellular populations in a robust manner.
Cellular processes are complex and result from the interplay between multiple cell types and their environment. Existing cell biology techniques often do not allow for accurate interpretation of this interplay. Using a quantitative imaging-based approach, we present a high-content protocol for characterizing the dynamic phenotypic responses (i.e. morphology changes, proliferation, apoptosis) of heterogeneous cell populations to changes in environmental stimuli. We highlight our ability to distinguish between cell types based upon either fluorescence intensity or inherent morphology features depending on the application. This platform allows for a more comprehensive characterization of subpopulation response to perturbation while utilizing shorter time, smaller amounts of reagents, and lower likelihood of error than traditional cell biology assays. However, in some cases, cell populations may be difficult to identify and quantitate based on complex cellular features and will require additional troubleshooting; we highlight some of these circumstances in the protocol. We demonstrate this application using response to drug in a cancer model; however, it can easily be applied more broadly to other physiological processes. This protocol allows one to identify subpopulations within a co-culture system and characterize the particular response of each to external stimuli.
Cell-based assays have been a workhorse in basic research and drug development settings. However, the limitations of these standard assays have become increasingly apparent with the discordance between in vitro and clinical data and the failure of most drugs to receive FDA approval. Here, we present a novel method for utilizing quantitative imaging to simultaneously analyze heterocellular phenotypes in response to relevant, co-occurring environmental stimuli.
Traditional cell-based assays that are used to measure cell viability include: trypan blue exclusion assays, MTT/MTS, and Annexin V-FITC flow cytometry staining. Trypan blue exclusion assays, while simple and inexpensive, require a large number of cells, are time consuming, and are often influenced by user bias1. MTT and MTS assays indirectly measure cell viability through measurements of mitochondrial metabolic rate. However, the metabolic activity of cells can be affected by different culture conditions (such as media or oxygen concentration), which leads to inaccurate results and prevents standardization across cell types and conditions2,3. Another major disadvantage of these techniques is their inability to distinguish between multiple cell types – most biological systems are heterocellular. While flow cytometry methods have the ability to distinguish between multiple cell populations, cell labels are required, dynamic sampling is challenging, and when using adherent cells, this application becomes time-consuming and error prone.
Other important cellular phenotypes, including morphological changes, occur in response to environmental stimuli but are not captured by traditional cell-based assays. Profiling cell states through morphological characterization and mapping similarities across samples is a powerful, unbiased tool with the ability to provide novel insights into many aspects of basic and translational research, including basic cell biology and drug discovery4. Furthermore, tumor cell morphology has been shown to correlate with tumor subtypes5 and aggressiveness6. Hence, it is of great interest to study these cellular features and how they relate to specific environmental perturbations. Additionally, one can use differences in morphological features to discriminate between subpopulations in co-culture systems. Fluorescently labeling cells has downfalls (i.e. altering inherent cell properties, time-consuming) and therefore additional methods to classify cell types are advantageous.
Microscopy-based imaging is an alternative method for profiling cellular phenotypes in a multiplexed, quantitative, and robust manner. In this manuscript, we apply our quantitative imaging pipeline to highlight the evolutionary dynamics of heterogeneous cell populations within a tumor. We focus on the interaction between Non-Small Cell Lung Cancer (NSCLC) cells and Cancer-Associated Fibroblasts (CAFs), the most prevalent stromal cell type found in tumors. CAFs have been implicated in tumor initiation, progression, and therapeutic response; therefore, performing phenotypic assays on tumor cells in the absence of CAFs can be misleading7,8,9. Specifically, we evaluated effects of CAFs on tumor cells in response to erlotinib, a small molecule targeting the Epidermal Growth Factor Receptor (EGFR) that is often used in clinical treatment of NSCLC. We utilized a high-content screening platform and its accompanying image analysis software for evaluation; however, in an attempt to make this methodology accessible to other researchers we have also developed a comparable downstream protocol using the open-source software: CellProfiler10 and CellProfiler Analyst11. Most image-based high-content screening assays are analyzed with commercialized software specific to a given instrument model. Results are difficult to replicate in other labs with different software because the underlying algorithms are often proprietary. Using this image-based pipeline, cell proliferation, death, and morphology of each subpopulation of a heterocellular culture in response to drug treatment using both fluorescence- and morphology-based classification were measured. The following protocol provides a robust methodology for probing complex cellular processes.
1. Cell Culture
NOTE: Here phenotypic responses of NSCLC cells (H3255) to the EGFR-targeting agent, erlotinib, when co-cultured with lung fibroblasts (CCD-19Lu GFP) were investigated. For the same data set, fluorescence-based and morphology-based classification of the two cell populations (H3255 and CCD-19Lu GFP) to exemplify their concordance was performed. However, only one classification method needs to be used and should be chosen based upon the application.
2. Preparation of Cells
3. Cell Plating
4. Drug Dosing
5. Image Acquisition
NOTE: Images were acquired on days 0, 2 and 3. Depending on cell types and cellular processes being studied, other desired time points may be studied.
6. Image Analysis
NOTE: All image analysis was performed using proprietary software. However, because this is not a publicly available, comparable analyses were also designed on CellProfiler 2.2 and CellProfiler Analyst 2.0, with a brief protocol listed below and detailed protocol provided as supplementary material (test images are already loaded into pipelines to test workflow). The images from this experiment were grouped together on a per-well basis so that each well could be loaded individually. The CellProfiler pipelines below contain a regular-expression that parses metadata information from each image filename and allows the images to be further grouped by channel.
We generated an image set consisting of 25 fields/well, 54 wells/plate (3 cell populations x 6 drug concentrations x 3 replicates), across three plates for a total of 4,050 individual images. The image sets generated over the course of the experiment were analyzed using proprietary software (see table of materials) to extract various quantitative properties of cells (i.e. morphology, fluorescence) which could then be used to classify cell subpopulations. However, because the commercial software used has limited access, comparable downstream pipelines in CellProfiler and CellProfiler Analyst were created.
Heterocellular Classification into Subpopulations
Nuclei were identified and segmented based on the DNA stain (here Hoechst) and cell populations were classified either based upon fluorescence or morphology (Figure 1). For fluorescence-based classification, the fibroblasts (CCD-19Lu) were previously transduced with GFP-lentivirus. The GFP intensity levels were measured for each nucleus, and those that were calculated above the accepted threshold (based on the background signal) were classified as CCD-19Lu while those below were identified as tumor cells (H3255). For morphology-based classification, cells were previously stained with a non-toxic cellular stain (see the table of materials) and this was used to identify and segment the cytoplasm. A machine learning algorithm was trained with ~50-100 cells from each population. Morphological features were identified that were significantly different between the populations, which were then used to design a linear classifier to distinguish between CCD-19Lu and H3255 cells. The fluorescence and morphology classification protocols were 97.4% (n = 1403) concordant at distinguishing between the two cell populations in untreated conditions and 92.5% (n = 916) concordant in drug treated conditions (1 µM erlotinib) (Figure 2).
Phenotypic Analyses of Subpopulations
In addition to discriminating between cell types, we aimed to characterize phenotypic properties of each subpopulation. Multiplexing assays saves time and reagents, adds consistency, and provides additional information regarding the system being studied. There are many potential phenotypic outputs and one should choose them based upon the questions of interest. Here, changes in the cell morphology and viability status in response to erlotinib treatment were investigated. After three days of drug treatment, a decrease in nuclear area and an increase in cellular area of the H3255 cells (Figure 3A) was observed. The mean difference in nuclear area between the "no drug" and "drug" treated populations was found to be statistically significant via a two-sided type-2 (equal variance) t-test (p = 7.92 x 10-16). We hypothesize that this observation is a cellular response to the stress imposed by drug treatment.
It is also of interest to study whether a drug has a cytotoxic (i.e. increase in number of dead cells over time) or cytostatic (i.e. decrease in number of cell births over time) effect on cells, as this has profound clinical impact. For example, a cytostatic drug effect induces growth arrest yet does not eliminate the cells from the tumor, thus there is the potential for cancer cells to reinitiate cell proliferation once the drug is removed. Drug effects can often be context, concentration, and cell type dependent. We previously observed erlotinib eliciting a cytotoxic response in one cell type, while showing a cytostatic response in another13.
Traditional viability assays output relative cell number and therefore, do not discriminate between growth arrests and cell death. Herein, dead cells were identified based upon propidium iodide stain (Figure 3B). Both cytotoxic and cytostatic effects of erlotinib on H3255 cells were observed, with an increase in the number of deaths and a decrease in the number of births following drug treatment (Figure 3C). It is worth noting that the number of dead cells drops following day 1 likely due to cell clearance. CCD-19Lu cells were not affected by the drug. An additional advantage of this platform is the generation of quantitative data. For example, in our co-culture experiment, an initial subpopulation of 1,118 (75.8%) H3255 cells was found to be 2,817 (87.9%) or 396 (57.2%) after three days without or with erlotinib treatment, respectively (Figure 4). Because we can generate actual cell counts instead of relative percentage (as with flow cytometry methods), we conclude that the change in composition during drug treatment is due to a decrease in H3255 cells and not an increase in CCD-19Lu. It is worth nothing that death rates may be underestimated due to cell clearance, which is difficult to assess experimentally and likely differs across cell types.
Figure 1: Overview of the Image Analysis Protocol. Two potential downstream image analysis pipelines to classify heterocellular populations using either morphology-based classification or fluorescence-based classification. Scale bars = 100 µm. Please click here to view a larger version of this figure.
Figure 2: Concordance between Morphology and Fluorescence-based Classification. (A) Concordance plot displaying the overlap of the two classification protocols. The same cells were classified as H3255 using both morphology- and fluorescence-based classification. The two protocols were in agreement, with classification for 97.4% (n = 1403) of untreated cells and 92.5% (n = 916) of cells treated with erlotinib (Note: white area is too small to visualize). (B) 10X images depicting examples of good and poor concordance between fluorescence-based and morphology-based classification. The white arrows point out cells that were inconsistently classified between platforms. Input image: blue – nuclei (Hoechst); green – CCD19Lu (GFP). Classification images: Red – H3255; green – CCD-19Lu. Scale bar = 100 µm. Please click here to view a larger version of this figure.
Figure 3: Multiplexed Phenotypic Measurements from a Single Experimental Setup. (A) Morphological features, such as nuclei and cell area, were calculated on the single cell level in the presence and absence of drug. Note: Cell areas measuring smaller than 100 µm2 were considered debris and excluded from analyses. Box plot depicts median with first and third quartile ranges and 95% confidence interval error bars. (B) H3255 (blue) and CCD-19Lu (green) cells were co-cultured and dead cells were identified based upon the intensity of propidium iodide stain (red) and imaged using a 10X objective. Scale bar = 1 mm (top panel); 100 µm (bottom images). (C) Total number of live and dead cells were calculated over three days with or without drug treatment, with an obvious decrease in number of live cells and increase in dead cells with the addition of erlotinib. Error bars represent standard error of the mean based on of three replicates. Please click here to view a larger version of this figure.
Figure 4: Subpopulation Dynamics over Time. Representative 10X images of wells containing H3255 (blue) and CCD-19Lu (green) on day 0 or day 3 with and without drug. Cells belonging to each subpopulation were counted and proportional pie charts show the actual change in population composition across samples. Scale bars = 1 mm (middle panels, shpwn in ;eftmost panel), 1 mm (top image), 100 µm (bottom images). Please click here to view a larger version of this figure.
The protocol described above improves upon current cell biology assays by providing more comprehensive insights into phenotypic dynamics of multiple cell types in response to environmental perturbations while using reduced reagents and time. A major advantage of this experimental design is the ability to analyze multiple phenotypes with a single setup and generate quantitative data characterizing these phenotypes on a single cell level. One technical advantage to this platform is the ease of initial troubleshooting compared to other assays. Because this method is image-based, one is able to visualize the wells for apparent over/under seeding. It is advisable to have cells in the exponential growth phase for the duration of the experiment and not be limited by nutrient or spatial constraints or confounded by senescence due to scarce seeding. Otherwise, birth and death rates may not be reproducible between experiments. For reference, CellPD is a publicly available program for computation of birth and death rates14. Additionally, one can visualize whether adequate concentrations of dyes were added to each well. Pipetting issues in an individual well could result in missegmentation and skewed data, but can easily be detected with the aforementioned protocol.
Unfortunately, not all cell types may be amenable to this application. It is important to be able to accurately segment the nuclei and cells, therefore analysis of cells that organize in more sphere-like or clumped structures may not be suitable. For some cells, it also may be advantageous to use a cell strainer prior to seeding to ensure initial seeding of single cells. In addition, the linear classifier technique is only applicable to cell types that can be readily distinguished based on morphology features.
For the success of the protocol, it is important to first optimize the imaging conditions, as the validity of downstream analyses is dependent on the quality of the images. While here we performed experiments using a high-content screening platform, image acquisition can also be performed using any fluorescent microscope (although an automated imaging platform is ideal for high-throughput approaches). Test images should be taken prior to each imaging time point to ensure that there are no problems with the microscope or protocol. The signal to noise ratio should be high, especially for the channels that will be used for segmentation (i.e. nuclei, cell stains). Additionally, it is important to image in the optimal plane of focus. If the images are out of focus, segmentation becomes much more difficult and the calculated morphological features will likely be inaccurate. Illumination differences between fields can cause problems with image segmentation as well. Large differences in brightness make the automated selection of threshold values difficult. Additionally, if there are heterogeneous fluorescence intensities between cells, a single threshold value may not sufficiently segment all the cells in an image. In this analysis, these problems were overcome by creating masks around brighter and dimmer cell populations and segmenting each population separately.
While the phenotypes under investigation in this protocol are limited to live, dead, and morphological characterization, they can easily be expanded to investigate other features. For example, functional genetic studies can be added with RNAi, overexpression, or other chemical perturbations.
In this paper, the capabilities of the protocol to measure the response of non-small cell lung cancer cells to erlotinib in the presence and absence of CAFs was demonstrated. However, this is merely one example of the many cell types and microenvironmental parameters that can be tested. We have extended this protocol to be used with other cell types and drug studies, including primary cells isolated from patient tumors13,15.
The authors have nothing to disclose.
This work was funded by National Cancer Institute (NCI) Grants U54CA143798 and U54CA143907 to establish Physical Sciences-Oncology Centers (PS-OCs) at the Dana-Farber Cancer Institute and University of Southern California, respectively. S.M. Mumenthaler received a PS-OC transnetwork award that supported some of this work.
We would like to express our deepest gratitude to our philanthropic supporters, particularly the Stephenson family, Emmet, Toni and Tessa, for their donation of the Operetta HCS platform. We would also like to thank J. Foo for guidance, and the Center for Applied Molecular Medicine team members: D. Agus for clinical guidance and mentorship, K. Patsch for meaningful discussions with experimental design, R. Rawat for aid in image analysis protocols, J. Katz for technical assistance with the Operetta, and P. Macklin and D. Ruderman for helpful discussions and feedback.
RPMI-1640 | Corning | 10-040-CV | cell culture medium |
FBS | Gemini | 100-106 | medium supplement |
Penicillin/streptomycin | Gibco | 15-140-122 | medium supplement |
TC20 | Biorad | 1450102 | cell counter |
TC20 slides with trypan blue | Biorad | 1450003 | cell counter slides |
96-well plates | Corning | 3904 | clear bottom black plates |
Erlotinib | LC Laboratories | E-4007 | |
DMSO | VWR | 317275-100ML | solution to resuspend drug |
Hoechst | Invitrogen | H21491 | nuclear dye |
Propidium Iodide | Invitrogen | P1304MP | dead cell stain |
Cell Tracker Orange CMRA | Life Technologies | C34551 | whole cell stain |
Operetta high content imaging system | Perkin Elmer | ||
CellProfiler | Broad Institue | version 2.2.0 (rev 9969f42) | http://cellprofiler.org/releases/ |
Cell culture incubator | Any cell culture incubator will be suitable – cells were cultured under 37 ºC at 5% CO2. | ||
15 mL Falcon conical tubes | Falcon | 14-959-53A | |
10 cm2 cell culture plates | TPP | 93040 |