Non-local Parabolic and Hyperbolic Models for Cell Polarisation in Heterogeneous Cancer Cell Populations.
Ontology highlight
ABSTRACT:
SUBMITTER: Bitsouni V
PROVIDER: S-EPMC6153854 | biostudies-literature | 2018 Oct
REPOSITORIES: biostudies-literature
Publications
Non-local Parabolic and Hyperbolic Models for Cell Polarisation in Heterogeneous Cancer Cell Populations.
Bitsouni Vasiliki V Eftimie Raluca R
Bulletin of mathematical biology 20180822 10
Tumours consist of heterogeneous populations of cells. The sub-populations can have different features, including cell motility, proliferation and metastatic potential. The interactions between clonal sub-populations are complex, from stable coexistence to dominant behaviours. The cell-cell interactions, i.e. attraction, repulsion and alignment, processes critical in cancer invasion and metastasis, can be influenced by the mutation of cancer cells. In this study, we develop a mathematical model ...[more]
Similar Datasets
Project description:BackgroundCertain biological processes, such as the development of cancer and immune activation, can be controlled by rare cellular events that are difficult to capture computationally through simulations of individual cells. Information about such rare events can be gleaned from an attractor analysis, for which a variety of methods exist (in particular for Boolean models). However, explicitly simulating a defined mixed population of cells in a way that tracks even the rarest subpopulations remains an open challenge.ResultsHere we show that when cellular states are described using a Boolean network model, one can exactly simulate the dynamics of non-interacting, highly heterogeneous populations directly, without having to model the various subpopulations. This strategy captures even the rarest outcomes of the model with no sampling error. Our method can incorporate heterogeneity in both cell state and, by augmenting the model, the underlying rules of the network as well (e.g., introducing loss-of-function genetic alterations). We demonstrate our method by using it to simulate a heterogeneous population of Boolean networks modeling the T-cell receptor, spanning ∼ 1020 distinct cellular states and mutational profiles.ConclusionsWe have developed a method for using Boolean models to perform a population-level simulation, in which the population consists of non-interacting individuals existing in different states. This approach can be used even when there are far too many distinct subpopulations to model individually.
| S-EPMC5992775 | biostudies-literature
Project description:Computational models are most impactful when they explain and characterize biological phenomena that are non-intuitive, unexpected, or difficult to study experimentally. Countless equation-based models have been built for these purposes, but we have yet to realize the extent to which rules-based models offer an intuitive framework that encourages computational and experimental collaboration. We develop ARCADE, a multi-scale agent-based model to interrogate emergent behavior of heterogeneous cell agents within dynamic microenvironments and demonstrate how complexity of intracellular metabolism and signaling modules impacts emergent dynamics. We perform in silico case studies on context, competition, and heterogeneity to demonstrate the utility of our model for gaining computational and experimental insight. Notably, there exist (i) differences in emergent behavior between colony and tissue contexts, (ii) linear, non-linear, and multimodal consequences of parameter variation on competition in simulated co-cultures, and (iii) variable impact of cell and population heterogeneity on emergent outcomes. Our extensible framework is easily modified to explore numerous biological systems, from tumor microenvironments to microbiomes.
| S-EPMC7301008 | biostudies-literature
Project description:We develop a combined experimental-mathematical framework to investigate heterogeneity in the context of breast cancer treated with doxorubicin. We engineer a cell line to over-express the multi-drug resistance 1 protein (MDR1), an ATP-dependent pump that effluxes intracellular drug. Co-culture experiments mixing the MDR1-overexpressing line with its parental line are evaluated via fluorescence microscopy. To quantify the impact of population heterogeneity on therapy response, these data are analyzed with a coupled pharmacokinetics/pharmacodynamics model. The proliferation and death rates of each line vary with co-culture condition (the relative fraction of each cell line at the time of seeding). For example, the death rate in the parental line under low-dose doxorubicin treatment is increased from 0.64 (± 0.22) × 10-2 to 1.46 (± 0.58) × 10-2 h-1 with increasing fractions of MDR1-overexpressing cells. The growth rate of the MDR1-overexpressing line increases 29% as its relative fraction is decreased. Simulations of the pharmacokinetics/pharmacodynamics model suggest increased efflux from MDR1-overexpressing cells contributes to the increased death rate in the parental cells. Experimentally, the death rate of parental cells is constant across co-culture conditions under co-treatment with an MDR1 inhibitor. These data indicate that intercellular pharmacokinetic variability should be considered in analyzing treatment response in heterogeneous populations.
| S-EPMC5935587 | biostudies-literature
Project description:Certain tumor phenomena, like metabolic heterogeneity and local stable regions of chronic hypoxia, signify a tumor's resistance to therapy. Although recent research has shed light on the intracellular mechanisms of cancer metabolic reprogramming, little is known about how tumors become metabolically heterogeneous or chronically hypoxic, namely the initial conditions and spatiotemporal dynamics that drive these cell population conditions. To study these aspects, we developed a minimal, spatially-resolved simulation framework for modeling tissue-scale mixed populations of cells based on diffusible particles the cells consume and release, the concentrations of which determine their behavior in arbitrarily complex ways, and on stochastic reproduction. We simulate cell populations that self-sort to facilitate metabolic symbiosis, that grow according to tumor-stroma signaling patterns, and that give rise to stable local regions of chronic hypoxia near blood vessels. We raise two novel questions in the context of these results: (1) How will two metabolically symbiotic cell subpopulations self-sort in the presence of glucose, oxygen, and lactate gradients? We observe a robust pattern of alternating striations. (2) What is the proper time scale to observe stable local regions of chronic hypoxia? We observe the stability is a function of the balance of three factors related to O2-diffusion rate, local vessel release rate, and viable and hypoxic tumor cell consumption rate. We anticipate our simulation framework will help researchers design better experiments and generate novel hypotheses to better understand dynamic, emergent whole-tumor behavior.
| S-EPMC5193460 | biostudies-literature
Project description:Estimation of division and death rates of lymphocytes in different conditions is vital for quantitative understanding of the immune system. Deuterium, in the form of deuterated glucose or heavy water, can be used to measure rates of proliferation and death of lymphocytes in vivo. Inferring these rates from labeling and delabeling curves has been subject to considerable debate with different groups suggesting different mathematical models for that purpose. We show that the three most common models, which are based on quite different biological assumptions, actually predict mathematically identical labeling curves with one parameter for the exponential up and down slope, and one parameter defining the maximum labeling level. By extending these previous models, we here propose a novel approach for the analysis of data from deuterium labeling experiments. We construct a model of "kinetic heterogeneity" in which the total cell population consists of many sub-populations with different rates of cell turnover. In this model, for a given distribution of the rates of turnover, the predicted fraction of labeled DNA accumulated and lost can be calculated. Our model reproduces several previously made experimental observations, such as a negative correlation between the length of the labeling period and the rate at which labeled DNA is lost after label cessation. We demonstrate the reliability of the new explicit kinetic heterogeneity model by applying it to artificially generated datasets, and illustrate its usefulness by fitting experimental data. In contrast to previous models, the explicit kinetic heterogeneity model 1) provides a novel way of interpreting labeling data; 2) allows for a non-exponential loss of labeled cells during delabeling, and 3) can be used to describe data with variable labeling length.
| S-EPMC2816685 | biostudies-literature
Project description:Biological systems are endowed with two fundamental but seemingly contradictory properties: robustness, the ability to withstand environmental fluctuations and genetic variability; and evolvability, the ability to acquire selectable and heritable phenotypic changes. Cell populations with heterogeneous genetic makeup, such as those of infectious microbial organisms or cancer, rely on their inherent robustness to maintain viability and fitness, but when encountering environmental insults, such as drug treatment, these populations are also poised for rapid adaptation through evolutionary selection. In this study, we develop a general mathematical model that allows us to explain and quantify this fundamental relationship between robustness and evolvability of heterogeneous cell populations. Our model predicts that robustness is, in fact, essential for evolvability, especially for more adverse environments, a trend we observe in aneuploid budding yeast and breast cancer cells. Robustness also compensates for the negative impact of the systems' complexity on their evolvability. Our model also provides a mathematical means to estimate the number of independent processes underlying a system's performance and identify the most generally adapted subpopulation, which may resemble the multi-drug-resistant "persister" cells observed in cancer.
| S-EPMC5994894 | biostudies-other
Project description:Recent studies suggest trends toward more mastectomies for primary breast cancer treatment. We assessed survival after mastectomy and breast-conserving surgery (BCS) with radiation for early-stage breast cancer among non-selected populations of women and among women similar to those in clinical trials. Using population-based data from Surveillance Epidemiology, and End Results cancer registries linked with Medicare administrative data from 1992 to 2005, we conducted propensity score analysis of survival following primary therapy for early-stage breast cancer, including BCS with radiation, BCS without radiation, mastectomy with radiation, and mastectomy without radiation. Adjusted survival was greatest among women who had BCS with radiation (median survival = 10.98 years). Compared with this group, mortality was higher among women who had mastectomy without radiation (median survival 10.04 years, adjusted hazard ratio (HR) = 1.19, 95% confidence interval (CI) = 1.14-1.23), mastectomy with radiation (median survival 10.02 years, HR = 1.20, 95% CI = 1.14-1.27), and BCS without radiation (median survival 7.63 years, HR = 1.81, 95% CI = 1.70-1.92). Among women representative of those eligible for clinical trials (age ≤70 years, Charlson comorbidity score = 0/1, and stage 1 tumors), there were no differences in survival for women who underwent BCS with radiation or mastectomy. In conclusion, after careful adjustment for differences in patient, physician, and hospital characteristics, we found better survival for BCS with radiation versus mastectomy among older early-stage breast cancer patients, with no difference in survival for BCS with radiation versus mastectomy among women representative of those in clinical trials. These findings are reassuring in light of recent trends towards more aggressive primary breast cancer therapy.
| S-EPMC2924956 | biostudies-literature
Project description:Circulating tumor cells (CTCs) are commonly isolated from the blood by targeting the epithelial cell adhesion molecule (EpCAM) through positive selection. However, EpCAM can be downregulated during metastatic progression, or it can be initially not present. We designed the present prospective trial to characterize CTCs as well as other circulating cell populations in blood samples from women with metastatic breast cancer without EpCAM-dependent enrichment and/or isolation technology.A total of 32 patients with metastatic breast cancer were enrolled, and blood samples were processed using a previously described negative depletion immunomagnetic methodology. Samples from healthy volunteers were run as controls (n?=?5). Multistep sequential labeling was performed to label and fix cell-surface markers followed by permeabilization for cytokeratins (CK) 8, 18 and 19. Multiparametric flow cytometry (FCM) analysis was conducted using a BD LSR II flow cytometer or a BD FACSAria II or FACSAria III cell sorter. Immunocytochemical staining on postenrichment specimens for DAPI, EpCAM, CD45, CK, epidermal growth factor receptor and vimentin was performed. Expression of these markers was visualized using confocal microscopy (CM).CD45-negative/CK-positive (CD45- CK+) populations with EpCAM?+?and EpCAM?-?expression were identified with both FCM and CM from the negatively enriched patient samples. In addition, EpCAM?+?and EpCAM?-?populations that were CK?+?and coexpressing the pan-hematopoietic marker CD45 were also noted. There were more CK?+?EpCAM?-?events/ml than CK?+?EpCAM?+?events/ml in both the CD45- and CD45+ fractions (both statistically significant at P???0.0005). The number of CK?+?CD45- and CK?+?CD45+ events per milliliter in blood samples (regardless of EpCAM status) was higher in patient samples than in normal control samples (P???0.0005 and P???0.026, respectively). Further, a significant fraction of the CK?+?CD45+ events also expressed CD68, a marker associated with tumor-associated macrophages. Higher levels of CD45-CK?+?EpCAM?-?were associated with worse overall survival (P?=?0.0292).Metastatic breast cancer patients have atypical cells that are CK?+?EpCAM?-?circulating in their blood. Because a substantial number of these patients do not have EpCAM?+?CTCs, additional studies are needed to evaluate the role of EpCAM?-?circulating cells as a prognostic and predictive marker.
| S-EPMC4053256 | biostudies-literature
Project description:Intratumoral heterogeneity challenges existing paradigms for anti-cancer therapy. We have previously demonstrated that the human embryonic stem cells (hESC)-derived cellular microenvironment in immunocompromised mice, enables functional distinction of heterogeneous tumor cells, including cells which do not grow into a tumor in a conventional direct tumor xenograft platform. We have identified and characterized six cancer cell subpopulations each clonally expanded from a single cell, derived from human ovarian clear cell carcinoma of a single tumor, to demonstrate striking intratumoral phenotypic heterogeneity that is dynamically dependent on the tumor growth microenvironment. These cancer cell subpopulations, characterized as cancer stem cell subpopulations, faithfully recapitulate the full spectrum of histological phenotypic heterogeneity known for human ovarian clear cell carcinoma. Each of the six subpopulations displays a different level of morphologic and tumorigenic differentiation wherein growth in the hESC-derived microenvironment favors growth of CD44+/aldehyde dehydrogenase positive pockets of self-renewing cells that sustain tumor growth through a process of tumorigenic differentiation into CD44-/aldehyde dehydrogenase negative derivatives. Strikingly, these derivative cells display microenvironment-dependent plasticity with the capacity to restore self-renewal markers and CD44 expression. In the current study, we delineate the distinct gene expression and epigenetic profiles of two such subpopulations, representing extremes of phenotypic heterogeneity in terms of niche-dependent self-renewal and tumorigenic differentiation. By combining Gene Set Enrichment, Gene Ontology and Pathway-focused array analyses with methylation status, we propose a suite of robust differences in tumor self-renewal and differentiation pathways that underlie the striking intratumoral phenotypic heterogeneity which characterize this and other solid tumor malignancies.
| S-EPMC3866276 | biostudies-literature
Project description:Local adaptation is determined by the strength of selection and the level of gene flow within heterogeneous landscapes. The presence of benign habitat can act as an evolutionary stepping stone for local adaptation to challenging environments by providing the necessary genetic variation. At the same time, migration load from benign habitats will hinder adaptation. In a community context, interspecific competition is expected to select against maladapted migrants, hence reducing migration load and facilitating adaptation. As the interplay between competition and spatial heterogeneity on the joint ecological and evolutionary dynamics of populations is poorly understood, we performed an evolutionary experiment using the herbivore spider mite Tetranychus urticae as a model. We studied the species's demography and local adaptation in a challenging environment that consisted of an initial sink (pepper plants) and/or a more benign environment (cucumber plants). Half of the experimental populations were exposed to a competitor, the congeneric T. ludeni. We show that while spider mites only adapted to the challenging pepper environment when it was spatially interspersed with benign cucumber habitat, this adaptation was only temporary and disappeared when the populations in the benign cucumber environment were expanding and spilling-over to the challenging pepper environment. Although the focal species outcompeted the competitor after about two months, a negative effect of competition on the focal species's performance persisted in the benign environment. Adaptation to challenging habitat in heterogeneous landscapes thus highly depends on demography and source-sink dynamics, but also on competitive interactions with other species, even if they are only present for a short time span.
| S-EPMC6599998 | biostudies-literature