skip to main content


Title: A mathematical model for phenotypic heterogeneity in breast cancer with implications for therapeutic strategies
Inevitably, almost all cancer patients develop resistance to targeted therapy. Intratumour heterogeneity is a major cause of drug resistance. Mathematical models that explain experiments quantitatively are useful in understanding the origin of intratumour heterogeneity, which then could be used to explore scenarios for efficacious therapy. Here, we develop a mathematical model to investigate intratumour heterogeneity in breast cancer by exploiting the observation that HER2+ and HER2− cells could divide symmetrically or asymmetrically. Our predictions for the evolution of cell fractions are in quantitative agreement with single-cell experiments. Remarkably, the colony size of HER2+ cells emerging from a single HER2− cell (or vice versa), which occurs in about four cell doublings, also agrees with experimental results, without tweaking any parameter in the model. The theory explains experimental data on the responses of breast tumours under different treatment protocols. We then used the model to predict that, not only the order of two drugs, but also the treatment period for each drug and the tumour cell plasticity could be manipulated to improve the treatment efficacy. Mathematical models, when integrated with data on patients, make possible exploration of a broad range of parameters readily, which might provide insights in devising effective therapies.  more » « less
Award ID(s):
1708128
NSF-PAR ID:
10382586
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Journal of The Royal Society Interface
Volume:
19
Issue:
186
ISSN:
1742-5662
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The CDK4/6 inhibitor palbociclib blocks cell cycle progression in Estrogen receptor–positive, human epidermal growth factor 2 receptor–negative (ER+/HER2−) breast tumor cells. Despite the drug’s success in improving patient outcomes, a small percentage of tumor cells continues to divide in the presence of palbociclib—a phenomenon we refer to as fractional resistance. It is critical to understand the cellular mechanisms underlying fractional resistance because the precise percentage of resistant cells in patient tissue is a strong predictor of clinical outcomes. Here, we hypothesize that fractional resistance arises from cell-to-cell differences in core cell cycle regulators that allow a subset of cells to escape CDK4/6 inhibitor therapy. We used multiplex, single-cell imaging to identify fractionally resistant cells in both cultured and primary breast tumor samples resected from patients. Resistant cells showed premature accumulation of multiple G1 regulators including E2F1, retinoblastoma protein, and CDK2, as well as enhanced sensitivity to pharmacological inhibition of CDK2 activity. Using trajectory inference approaches, we show how plasticity among cell cycle regulators gives rise to alternate cell cycle “paths” that allow individual tumor cells to escape palbociclib treatment. Understanding drivers of cell cycle plasticity, and how to eliminate resistant cell cycle paths, could lead to improved cancer therapies targeting fractionally resistant cells to improve patient outcomes.

     
    more » « less
  2. null (Ed.)
    Heterogeneity is a hallmark of all cancers. Tumor heterogeneity is found at different levels — interpatient, intrapatient, and intratumor heterogeneity. All of them pose challenges for clinical treatments. The latter two scenarios can also increase the risk of developing drug resistance. Although the existence of tumor heterogeneity has been known for two centuries, a clear understanding of its origin is still elusive, especially at the level of intratumor heterogeneity (ITH). The coexistence of different subpopulations within a single tumor has been shown to play crucial roles during all stages of carcinogenesis. Here, using concepts from evolutionary game theory and public goods game, often invoked in the context of the tragedy of commons, we explore how the interactions among subclone populations influence the establishment of ITH. By using an evolutionary model, which unifies several experimental results in distinct cancer types, we develop quantitative theoretical models for explaining data from in vitro experiments involving pancreatic cancer as well as in vivo data in glioblastoma multiforme. Such physical and mathematical models complement experimental studies, and could optimistically provide new ideas for the design of efficacious therapies for cancer patients. 
    more » « less
  3. Abstract

    Intercellular mechanisms by which the stromal microenvironment contributes to solid tumor progression and targeted therapy resistance remain poorly understood, presenting significant clinical hurdles. PEAK1 (Pseudopodium-Enriched Atypical Kinase One) is an actin cytoskeleton- and focal adhesion-associated pseudokinase that promotes cell state plasticity and cancer metastasis by mediating growth factor-integrin signaling crosstalk. Here, we determined that stromal PEAK1 expression predicts poor outcomes in HER2-positive breast cancers high in SNAI2 expression and enriched for MSC content. Specifically, we identified that the fibroblastic stroma in HER2-positive breast cancer patient tissue stains positive for both nuclear SNAI2 and cytoplasmic PEAK1. Furthermore, mesenchymal stem cells (MSCs) and cancer-associated fibroblasts (CAFs) express high PEAK1 protein levels and potentiate tumorigenesis, lapatinib resistance and metastasis of HER2-positive breast cancer cells in a PEAK1-dependent manner. Analysis of PEAK1-dependent secreted factors from MSCs revealed INHBA/activin-A as a necessary factor in the conditioned media of PEAK1-expressing MSCs that promotes lapatinib resistance. Single-cell CycIF analysis of MSC-breast cancer cell co-cultures identified enrichment of p-Akthigh/p-gH2AXlow, MCL1high/p-gH2AXlowand GRP78high/VIMhighbreast cancer cell subpopulations by the presence of PEAK1-expressing MSCs and lapatinib treatment. Bioinformatic analyses on a PEAK1-centric stroma-tumor cell gene set and follow-up immunostaining of co-cultures predict targeting antiapoptotic and stress pathways as a means to improve targeted therapy responses and patient outcomes in HER2-positive breast cancer and other stroma-rich malignancies. These data provide the first evidence that PEAK1 promotes tumorigenic phenotypes through a previously unrecognized SNAI2-PEAK1-INHBA stromal cell axis.

     
    more » « less
  4. Abstract

    Drug resistance poses a significant challenge in cancer treatment. Despite the initial effectiveness of therapies such as chemotherapy, targeted therapy and immunotherapy, many patients eventually develop resistance. To gain deep insights into the underlying mechanisms, single-cell profiling has been performed to interrogate drug resistance at cell level. Herein, we have built the DRMref database (https://ccsm.uth.edu/DRMref/) to provide comprehensive characterization of drug resistance using single-cell data from drug treatment settings. The current version of DRMref includes 42 single-cell datasets from 30 studies, covering 382 samples, 13 major cancer types, 26 cancer subtypes, 35 treatment regimens and 42 drugs. All datasets in DRMref are browsable and searchable, with detailed annotations provided. Meanwhile, DRMref includes analyses of cellular composition, intratumoral heterogeneity, epithelial–mesenchymal transition, cell–cell interaction and differentially expressed genes in resistant cells. Notably, DRMref investigates the drug resistance mechanisms (e.g. Aberration of Drug’s Therapeutic Target, Drug Inactivation by Structure Modification, etc.) in resistant cells. Additional enrichment analysis of hallmark/KEGG (Kyoto Encyclopedia of Genes and Genomes)/GO (Gene Ontology) pathways, as well as the identification of microRNA, motif and transcription factors involved in resistant cells, is provided in DRMref for user’s exploration. Overall, DRMref serves as a unique single-cell-based resource for studying drug resistance, drug combination therapy and discovering novel drug targets.

     
    more » « less
  5. Tumor microenvironment is a complex niche consisting of cancer cells and stromal cells in a network of extracellular matrix proteins and various soluble factors. Dynamic interactions among cellular and non-cellular components of the tumor microenvironment regulate tumor initiation and progression. Fibroblasts are the most abundant stromal cell type and dynamically interact with cancer cells both in primary tumors and in metastases. Cancer cells activate resident fibroblasts to produce and secrete soluble signaling molecules that support proliferation, migration, matrix invasion, and drug resistance of cancer cell and tumor angiogenesis. In recent years, various forms of three-dimensional tumor models have been developed to study tumor–stromal interactions and to identify anti-cancer drugs that block these interactions. There is currently a technological gap in development of tumor models that are physiologically relevant, scalable, and allow convenient, on-demand addition of desired components of the tumor microenvironment. In this review, we discuss three studies from our group that focus on developing bioengineered models to study tumor-stromal signaling. We will present these studies chronologically and based on their increasing complexity. We will discuss the validation of the models using a CXCL12-CXCR4 chemokine-receptor signaling present among activated fibroblasts and breast cancer cells in solid tumors, highlight the advantages and shortcomings of the models, and conclude with our perspectives on their applications. Impact statement Tumor stroma plays an important role in progression of cancers to a fatal metastatic disease. Modern treatment strategies are considering targeting tumor stroma to improve outcomes for cancer patients. A current challenge to develop stroma-targeting therapeutics is the lack of preclinical physiologic tumor models. Animal models widely used in cancer research lack human stroma and are not amenable to screening of chemical compounds for cancer drug discovery. In this review, we outline in vitro three-dimensional tumor models that we have developed to study the interactions among cancer cells and stromal cells. We describe development of the tumor models in a modular fashion, from a spheroid model to a sophisticated organotypic model, and discuss the importance of using correct physiologic models to recapitulate tumor-stromal signaling. These biomimetic tumor models will facilitate understanding of tumor-stromal signaling biology and provide a scalable approach for testing and discovery of cancer drugs. 
    more » « less