Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPhenotypic heterogeneity in modeling cancer evolution
63
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The unwelcome evolution of malignancy during cancer progression emerges through a selection process in a complex heterogeneous population structure. In the present work, we investigate evolutionary dynamics in a phenotypically heterogeneous population of stem cells (SCs) and their associated progenitors. The fate of a malignant mutation is determined not only by overall stem cell and differentiated cell growth rates but also differentiation and dedifferentiation rates. We investigate the effect of such a complex population structure on the evolution of malignant mutations. We derive exact analytic results for the fixation probability of a mutant arising in each of the subpopulations. The analytic results are in almost perfect agreement with the numerical simulations. Moreover, a condition for evolutionary advantage of a mutant cell versus the wild type population is given in the present study. We also show that microenvironment-induced plasticity in invading mutants leads to more aggressive mutants with higher fixation probability. Our model predicts that decreasing polarity between stem and differentiated cells turnover would raise the survivability of non-plastic mutants; while it would suppress the development of malignancy for plastic mutants. We discuss our model in the context of colorectal/intestinal cancer (at the epithelium). This novel mathematical framework can be applied more generally to a variety of problems concerning selection in heterogeneous populations, in other contexts such as population genetics, and ecology.
rate research
Read More
A tumor can be thought of as an ecosystem, which critically means that we cannot just consider it as a collection of mutated cells but more as a complex system of many interacting cellular and microenvironmental elements. At its simplest, a growing tumor with increased proliferation capacity must compete for space as a limited resource. Hypercellularity leads to a contact-inhibited core with a competitive proliferating rim. Evolution and selection occurs, and an individual cells capacity to survive and propagate is determined by its combination of traits and interaction with the environment. With heterogeneity in phenotypes, the clone that will dominate is not always obvious as there are both local interactions and global pressures. Several combinations of phenotypes can coexist, changing the fitness of the whole. To understand some aspects of heterogeneity in a growing tumor we build an off-lattice agent based model consisting of individual cells with assigned trait values for proliferation and migration rates. We represent heterogeneity in these traits with frequency distributions and combinations of traits with density maps. How the distributions change over time is dependent on how traits are passed on to progeny cells, which is our main inquiry. We bypass the translation of genetics to behavior by focussing on the functional end result of inheritance of the phenotype combined with the environmental influence of limited space.
The paradigm of phenotypic plasticity indicates reversible relations of different cancer cell phenotypes, which extends the cellular hierarchy proposed by the classical cancer stem cell (CSC) theory. Since it is still question able if the phenotypic plasticity is a crucial improvement to the hierarchical model or just a minor extension to it, it is worthwhile to explore the dynamic behavior characterizing the reversible phenotypic plasticity. In this study we compare the hierarchical model and the reversible model in predicting the cell-state dynamics observed in biological experiments. Our results show that the hierarchical model shows significant disadvantages over the reversible model in describing both long-term stability (phenotypic equilibrium) and short-term transient dynamics (overshoot) of cancer cells. In a very specific case in which the total growth of population due to each cell type is identical, the hierarchical model predicts neither phenotypic equilibrium nor overshoot, whereas thereversible model succeeds in predicting both of them. Even though the performance of the hierarchical model can be improved by relaxing the specific assumption, its prediction to the phenotypic equilibrium strongly depends on a precondition that may be unrealistic in biological experiments, and it also fails to capture the overshoot of CSCs. By comparison, it is more likely for the reversible model to correctly describe the stability of the phenotypic mixture and various types of overshoot behavior.
Background: Tumours are diverse ecosystems with persistent heterogeneity in various cancer hallmarks like self-sufficiency of growth factor production for angiogenesis and reprogramming of energy-metabolism for aerobic glycolysis. This heterogeneity has consequences for diagnosis, treatment, and disease progression. Methods: We introduce the double goods game to study the dynamics of these traits using evolutionary game theory. We model glycolytic acid production as a public good for all tumour cells and oxygen from vascularization via VEGF production as a club good benefiting non-glycolytic tumour cells. This results in three viable phenotypic strategies: glycolytic, angiogenic, and aerobic non-angiogenic. Results: We classify the dynamics into three qualitatively distinct regimes: (1) fully glycolytic, (2) fully angiogenic, or (3) polyclonal in all three cell types. The third regime allows for dynamic heterogeneity even with linear goods, something that was not possible in prior public good models that considered glycolysis or growth-factor production in isolation. Conclusion: The cyclic dynamics of the polyclonal regime stress the importance of timing for anti-glycolysis treatments like lonidamine. The existence of qualitatively different dynamic regimes highlights the order effects of treatments. In particular, we consider the potential of vascular renormalization as a neoadjuvant therapy before follow up with interventions like buffer therapy.
We present a mathematical study of the emergence of phenotypic heterogeneity in vascularised tumours. Our study is based on formal asymptotic analysis and numerical simulations of a system of non-local parabolic equations that describes the phenotypic evolution of tumour cells and their nonlinear dynamic interactions with the oxygen, which is released from the intratumoural vascular network. Numerical simulations are carried out both in the case of arbitrary distributions of intratumour blood vessels and in the case where the intratumoural vascular network is reconstructed from clinical images obtained using dynamic optical coherence tomography. The results obtained support a more in-depth theoretical understanding of the eco-evolutionary process which underpins the emergence of phenotypic heterogeneity in vascularised tumours. In particular, our results offer a theoretical basis for empirical evidence indicating that the phenotypic properties of cancer cells in vascularised tumours vary with the distance from the blood vessels, and establish a relation between the degree of tumour tissue vascularisation and the level of intratumour phenotypic heterogeneity.
Inevitably, almost all cancer patients develop resistance to targeted therapy. Intratumor heterogeneity (ITH) is a major cause of drug resistance. Mathematical models that explain experiments quantitatively is useful in understanding the origin of ITH, which then could be used to explore scenarios for efficacious therapy. Here, we develop a mathematical model to investigate ITH 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, agrees perfectly with experimental results, without tweaking any parameter in the model. The theory quantitatively explains experimental data on the responses of breast cancer tumor 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 tumor 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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
