Tumor progression results from a complex interplay between cellular heterogeneity, treatment response, microenvironment and heterocellular interactions. Existing approaches to characterize this interplay suffer from an inability to distinguish between multiple cell types, often lack environmental context, and are unable to perform multiplex phenotypic profiling of cell populations. Here we present a high-throughput platform for characterizing, with single-cell resolution, the dynamic phenotypic responses (i.e. morphology changes, proliferation, apoptosis) of heterogeneous cell populations both during standard growth and in response to multiple, co-occurring selective pressures. The speed of this platform enables a thorough investigation of the impacts of diverse selective pressures including genetic alterations, therapeutic interventions, heterocellular components and microenvironmental factors. The platform has been applied to both 2D and 3D culture systems and readily distinguishes between (1) cytotoxic versus cytostatic cellular responses; and (2) changes in morphological features over time and in response to perturbation. These important features can directly influence tumor evolution and clinical outcome. Our image-based approach provides a deeper insight into the cellular dynamics and heterogeneity of tumors (or other complex systems), with reduced reagents and time, offering advantages over traditional biological assays.

Targeted agents have improved the efficacy of chemotherapy for cancer patients, however, there remains a lack of understanding of how these therapies affect the unsuspecting bystanders of the stromal microenvironment. Cetuximab, a monoclonal antibody therapy targeting the epidermal growth factor receptor (EGFR), is given in combination with chemotherapy as the standard of care for a subset of metastatic colorectal cancer patients. The overall response to this treatment is underwhelming and, while genetic mutations that confer resistance have been identified, it is still not known why this drug is ineffective for some patients. We discovered that cancer-associated fibroblasts (CAFs), a major cellular subset of the tumor stroma, can provide a source of cancer cell resistance. Specifically, we observed that upon treatment with cetuximab, CAFs increased their secretion of EGF, which was sufficient to render neighboring cancer cells resistant to cetuximab treatment through sustained mitogen-activated protein kinases (MAPK) signaling. Furthermore, we show the cetuximab-induced EGF secretion to be specific to CAFs and not to cancer cells or normal fibroblasts. Altogether, this work emphasizes the importance of the tumor microenvironment and considering the potential unintended consequences of therapeutically targeting cancer-driving proteins on non-tumorigenic cell types.

Experimental studies have shown that one key factor in driving the emergence of drug resistance in solid tumors is tumor hypoxia, which leads to the formation of localized environmental niches where drug-resistant cell populations can evolve and survive. Hypoxia-activated prodrugs (HAPs) are compounds designed to penetrate to hypoxic regions of a tumor and release cytotoxic or cytostatic agents; several of these HAPs are currently in clinical trial. However, preliminary results have not shown a survival benefit in several of these trials. We hypothesize that the efficacy of treatments involving these prodrugs depends heavily on identifying the correct treatment schedule, and that mathematical modeling can be used to help design potential therapeutic strategies combining HAPs with standard therapies to achieve long-term tumor control or eradication. We develop this framework in the specific context of EGFR-driven non-small cell lung cancer, which is commonly treated with the tyrosine kinase inhibitor erlotinib. We develop a stochastic mathematical model, parametrized using clinical and experimental data, to explore a spectrum of treatment regimens combining a HAP, evofosfamide, with erlotinib. We design combination toxicity constraint models and optimize treatment strategies over the space of tolerated schedules to identify specific combination schedules that lead to optimal tumor control. We find that (i) combining these therapies delays resistance longer than any monotherapy schedule with either evofosfamide or erlotinib alone, (ii) sequentially alternating single doses of each drug leads to minimal tumor burden and maximal reduction in probability of developing resistance, and (iii) strategies minimizing the length of time after an evofosfamide dose and before erlotinib confer further benefits in reduction of tumor burden. These results provide insights into how hypoxia-activated prodrugs may be used to enhance therapeutic effectiveness in the clinic.

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.