High-resolution deconstruction of evolution induced by chemotherapy treatments in breast cancer xenografts

The Power and Promise of Patient-Derived Xenografts of Human Breast Cancer

In 2016, a group of researchers engaged in the development of patient-derived xenografts (PDXs) of human breast cancer provided a comprehensive review of the state of the field. In that review, they summarized the clinical problem that PDXs might address, the technical approaches to their generation (including a discussion of host animals and transplant conditions tested), and presented transplantation success (take) rates across groups and across transplantation conditions. At the time, there were just over 500 unique PDX models created by these investigators representing all three clinically defined subtypes (ER+, HER2+, and TNBC). Today, many of these PDX resources have at least doubled in size, and several more PDX development groups now exist, such that there may be well upward of 1000 PDX models of human breast cancer in existence worldwide. They also presented a series of open questions for the field. Many of these questions have been addressed. However, several remain open, or only partially addressed. Herein, we revisit these questions, and recount the progress that has been made in a number of areas with respect to generation, characterization, and use of PDXs in translational research, and re-present questions that remain open. These open questions, and others, are now being addressed not only by individual investigators, but also large, well-funded consortia including the PDXNet program of the National Cancer Institute in the United States, and the EuroPDX Consortium, an organization of PDX developers across Europe. Finally, we discuss the new opportunities in PDX-based research.

High-throughput deconvolution of 3D organoid dynamics at cellular resolution for cancer pharmacology with Cellos

Three-dimensional (3D) organoid cultures are flexible systems to interrogate cellular growth, morphology, multicellular spatial architecture, and cellular interactions in response to treatment. However, computational methods for analysis of 3D organoids with sufficiently high-throughput and cellular resolution are needed. Here we report Cellos, an accurate, high-throughput pipeline for 3D organoid segmentation using classical algorithms and nuclear segmentation using a trained Stardist-3D convolutional neural network. To evaluate Cellos, we analyze ~100,000 organoids with ~2.35 million cells from multiple treatment experiments. Cellos segments dye-stained or fluorescently-labeled nuclei and accurately distinguishes distinct labeled cell populations within organoids. Cellos can recapitulate traditional luminescence-based drug response of cells with complex drug sensitivities, while also quantifying changes in organoid and nuclear morphologies caused by treatment as well as cell-cell spatial relationships that reflect ecological affinity. Cellos provides powerful tools to perform high-throughput analysis for pharmacological testing and biological investigation of organoids based on 3D imaging.

Cellos: High-throughput deconvolution of 3D organoid dynamics at cellular resolution for cancer pharmacology

Three-dimensional (3D) culture models, such as organoids, are flexible systems to interrogate cellular growth and morphology, multicellular spatial architecture, and cell interactions in response to drug treatment. However, new computational methods to segment and analyze 3D models at cellular resolution with sufficiently high throughput are needed to realize these possibilities. Here we report Cellos (Cell and Organoid Segmentation), an accurate, high throughput image analysis pipeline for 3D organoid and nuclear segmentation analysis. Cellos segments organoids in 3D using classical algorithms and segments nuclei using a Stardist-3D convolutional neural network which we trained on a manually annotated dataset of 3,862 cells from 36 organoids confocally imaged at 5 μm z-resolution. To evaluate the capabilities of Cellos we then analyzed 74,450 organoids with 1.65 million cells, from multiple experiments on triple negative breast cancer organoids containing clonal mixtures with complex cisplatin sensitivities. Cellos was able to accurately distinguish ratios of distinct fluorescently labelled cell populations in organoids, with ≤3% deviation from the seeding ratios in each well and was effective for both fluorescently labelled nuclei and independent DAPI stained datasets. Cellos was able to recapitulate traditional luminescence-based drug response quantifications by analyzing 3D images, including parallel analysis of multiple cancer clones in the same well. Moreover, Cellos was able to identify organoid and nuclear morphology feature changes associated with treatment. Finally, Cellos enables 3D analysis of cell spatial relationships, which we used to detect ecological affinity between cancer cells beyond what arises from local cell division or organoid composition. Cellos provides powerful tools to perform high throughput analysis for pharmacological testing and biological investigation of organoids based on 3D imaging.

Genomic and epigenomic BRCA alterations predict adaptive resistance and response to platinum-based therapy in patients with triple-negative breast and ovarian carcinomas

Triple-negative breast cancer (TNBC) and ovarian carcinomas (OvCas) with BRCA1 promoter methylation ( BRCA1 meth) respond more poorly to alkylating agents compared to those bearing mutations in BRCA1 and BRCA2 ( BRCA mut). This is a conundrum given the biologically equivalent homologous recombination deficiency (HRD) induced by these genetic and epigenetic BRCA perturbations. We dissected this problem through detailed genomic analyses of TNBC and OvCa cohorts and experimentation with patient-derived xenografts and genetically engineered cell lines. We found that despite identical downstream genomic mutational signatures associated with BRCA1 meth and BRCA mut states, BRCA1 meth uniformly associates with poor outcomes. Exposure of BRCA1 meth TNBCs to platinum chemotherapy, either as clinical treatment of a patient or as experimental in vivo exposure of preclinical patient derived xenografts, resulted in allelic loss of BRCA1 methylation and increased BRCA1 expression and platinum resistance. These data suggest that, unlike BRCA mut cancers, where BRCA loss is a genetically “fixed” deficiency state, BRCA1 meth cancers are highly adaptive to genotoxin exposure and, through reversal of promoter methylation, recover BRCA1 expression and become resistant to therapy. We further found a specific augmented immune transcriptional signal associated with enhanced response to platinum chemotherapy but only in patients with BRCA-proficient cancers. We showed how integrating both this cancer immune signature and the presence of BRCA mutations results in more accurate predictions of patient response when compared to either HRD status or BRCA status alone. This underscores the importance of defining BRCA heterogeneity in optimizing the predictive precision of assigning response probabilities in TNBC and OvCa.

The impact of the spatial heterogeneity of resistant cells and fibroblasts on treatment response

A long-standing practice in the treatment of cancer is that of hitting hard with the maximum tolerated dose to eradicate tumors. This continuous therapy, however, selects for resistant cells, leading to the failure of the treatment. A different type of treatment strategy, adaptive therapy, has recently been shown to have a degree of success in both preclinical xenograft experiments and clinical trials. Adaptive therapy is used to maintain a tumor's volume by exploiting the competition between drug-sensitive and drug-resistant cells with minimum effective drug doses or timed drug holidays. To further understand the role of competition in the outcomes of adaptive therapy, we developed a 2D on-lattice agent-based model. Our simulations show that the superiority of the adaptive strategy over continuous therapy depends on the local competition shaped by the spatial distribution of resistant cells. Intratumor competition can also be affected by fibroblasts, which produce microenvironmental factors that promote cancer cell growth. To this end, we simulated the impact of different fibroblast distributions on treatment outcomes. As a proof of principle, we focused on five types of distribution of fibroblasts characterized by different locations, shapes, and orientations of the fibroblast region with respect to the resistant cells. Our simulation shows that the spatial architecture of fibroblasts modulates tumor progression in both continuous and adaptive therapy. Finally, as a proof of concept, we simulated the outcomes of adaptive therapy of a virtual patient with four metastatic sites composed of different spatial distributions of fibroblasts and drug-resistant cell populations. Our simulation highlights the importance of undetected metastatic lesions on adaptive therapy outcomes.

Adaptive Therapy for Metastatic Melanoma: Predictions from Patient Calibrated Mathematical Models

Adaptive therapy is an evolution-based treatment approach that aims to maintain tumor volume by employing minimum effective drug doses or timed drug holidays. For successful adaptive therapy outcomes, it is critical to find the optimal timing of treatment switch points in a patient-specific manner. Here we develop a combination of mathematical models that examine interactions between drug-sensitive and resistant cells to facilitate melanoma adaptive therapy dosing and switch time points. The first model assumes genetically fixed drug-sensitive and -resistant popul tions that compete for limited resources. The second model considers phenotypic switching between drug-sensitive and -resistant cells. We calibrated each model to fit melanoma patient biomarker changes over time and predicted patient-specific adaptive therapy schedules. Overall, the models predict that adaptive therapy would have delayed time to progression by 6-25 months compared to continuous therapy with dose rates of 6-74% relative to continuous therapy. We identified predictive factors driving the clinical time gained by adaptive therapy, such as the number of initial sensitive cells, competitive effect, switching rate from resistant to sensitive cells, and sensitive cell growth rate. This study highlights that there is a range of potential patient-specific benefits of adaptive therapy and identifies parameters that modulate this benefit.

Conservation of copy number profiles during engraftment and passaging of patient-derived cancer xenografts

Patient-derived xenografts (PDXs) are resected human tumors engrafted into mice for preclinical studies and therapeutic testing. It has been proposed that the mouse host affects tumor evolution during PDX engraftment and propagation, affecting the accuracy of PDX modeling of human cancer. Here, we exhaustively analyze copy number alterations (CNAs) in 1,451 PDX and matched patient tumor (PT) samples from 509 PDX models. CNA inferences based on DNA sequencing and microarray data displayed substantially higher resolution and dynamic range than gene expression-based inferences, and they also showed strong CNA conservation from PTs through late-passage PDXs. CNA recurrence analysis of 130 colorectal and breast PT/PDX-early/PDX-late trios confirmed high-resolution CNA retention. We observed no significant enrichment of cancer-related genes in PDX-specific CNAs across models. Moreover, CNA differences between patient and PDX tumors were comparable to variations in multiregion samples within patients. Our study demonstrates the lack of systematic copy number evolution driven by the PDX mouse host.

High accuracy gene expression profiling of sorted cell subpopulations from breast cancer PDX model tissue

Intratumor Heterogeneity (ITH) is a functionally important property of tumor tissue and may be involved in drug resistance mechanisms. Although descriptions of ITH can be traced back to very early reports about cancer tissue, mechanistic investigations are still limited by the precision of analysis methods and access to relevant tissue sources. PDX models have provided a reproducible source of tissue with at least a partial representation of naturally occurring ITH. We investigated the properties of phenotypically distinct cell populations by Fluorescence activated cell sorting (FACS) tissue derived cells from multiple tumors from a triple negative breast cancer patient derived xenograft (PDX) model. We subsequently subjected each population to in depth gene expression analysis. Our findings suggest that process related gene expression changes (caused by tissue dissociation and FACS sorting) are restricted to Immediate Early Genes (IEGs). This allowed us to discover highly reproducible gene expression profiles of distinct cellular compartments identifiable by cell surface markers in this particular tumor model. Within the context of data from a previously published model our work suggests that gene expression profiles associated with hypoxia, stemness and drug resistance may reside in tumor subpopulations predictably growing in PDX models. This approach provides a novel opportunity for prospective mechanistic studies of ITH.

Unveiling functional heterogeneity in breast cancer multicellular tumor spheroids through single-cell RNA-seq

Heterogeneity is an intrinsic characteristic of cancer. Even in isogenic tumors, cell populations exhibit differential cellular programs that overall supply malignancy and decrease treatment efficiency. In this study, we investigated the functional relationship among cell subtypes and how this interdependency can promote tumor development in a cancer cell line. To do so, we performed single-cell RNA-seq of MCF7 Multicellular Tumor Spheroids as a tumor model. Analysis of single-cell transcriptomes at two-time points of the spheroid growth, allowed us to dissect their functional relationship. As a result, three major robust cellular clusters, with a non-redundant complementary composition, were found. Meanwhile, one cluster promotes proliferation, others mainly activate mechanisms to invade other tissues and serve as a reservoir population conserved over time. Our results provide evidence to see cancer as a systemic unit that has cell populations with task stratification with the ultimate goal of preserving the hallmarks in tumors.

Targeting MYCN-expressing triple-negative breast cancer with BET and MEK inhibitors

Triple-negative breast cancer (TNBC) is an aggressive form of breast cancer that does not respond to endocrine therapy or human epidermal growth factor receptor 2 (HER2)-targeted therapies. Individuals with TNBC experience higher rates of relapse and shorter overall survival compared to patients with receptor-positive breast cancer subtypes. Preclinical discoveries are needed to identify, develop, and advance new drug targets to improve outcomes for patients with TNBC. Here, we report that MYCN, an oncogene typically overexpressed in tumors of the nervous system or with neuroendocrine features, is heterogeneously expressed within a substantial fraction of primary and recurrent TNBC and is expressed in an even higher fraction of TNBCs that do not display a pathological complete response after neoadjuvant chemotherapy. We performed high-throughput chemical screens on TNBC cell lines with varying amounts of MYCN expression and determined that cells with higher expression of MYCN were more sensitive to bromodomain and extraterminal motif (BET) inhibitors. Combined BET and MEK inhibition resulted in a synergistic decrease in viability, both in vitro and in vivo, using cell lines and patient-derived xenograft (PDX) models. Our preclinical data provide a rationale to advance a combination of BET and MEK inhibitors to clinical investigation for patients with advanced MYCN-expressing TNBC.

Molecular Biology and Evolution of Cancer: From Discovery to Action

Cancer progression is an evolutionary process. During this process, evolving cancer cell populations encounter restrictive ecological niches within the body, such as the primary tumor, circulatory system, and diverse metastatic sites. Efforts to prevent or delay cancer evolution-and progression-require a deep understanding of the underlying molecular evolutionary processes. Herein we discuss a suite of concepts and tools from evolutionary and ecological theory that can inform cancer biology in new and meaningful ways. We also highlight current challenges to applying these concepts, and propose ways in which incorporating these concepts could identify new therapeutic modes and vulnerabilities in cancer.

Treating Cancer as an Invasive Species

To cure a patient's cancer is to eradicate invasive cells from the ecosystem of the body. However, the ecologic complexity of this challenge is not well understood. Here we show how results from eradications of invasive mammalian species from islands-one of the few contexts in which invasive species have been regularly cleared-inform new research directions for treating cancer. We first summarize the epidemiologic characteristics of island invader eradications and cancer treatments by analyzing recent datasets from the Database of Invasive Island Species Eradications and The Cancer Genome Atlas, detailing the superior successes of island eradication projects. Next, we compare how genetic and environmental factors impact success in each system. These comparisons illuminate a number of promising cancer research and treatment directions, such as heterogeneity engineering as motivated by gene drives and adaptive therapy; multiscale analyses of how population heterogeneity potentiates treatment resistance; and application of ecological data mining techniques to high-throughput cancer data. We anticipate that interdisciplinary comparisons between tumor progression and invasive species would inspire development of novel paradigms to cure cancer.