Dissecting heterogeneous cell populations across drug and disease conditions with PopAlign.
Ontology highlight
ABSTRACT:
SUBMITTER: Chen S
PROVIDER: S-EPMC7682438 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
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:The effect of scheduling of targeted therapy combinations on drug resistance is underexplored in triple-negative breast cancer (TNBC). TNBC constitutes heterogeneous cancer cell populations the composition of which can change dynamically during treatment resulting in the selection of resistant clones with a fitness advantage. We evaluated crizotinib (ALK/MET inhibitor) and navitoclax (ABT-263; Bcl-2/Bcl-xL inhibitor) combinations in a large design consisting of 696 two-cycle sequential and concomitant treatment regimens with varying treatment dose, duration, and drug holiday length over a 26-day period in MDA-MB-231 TNBC cells and found that patterns of resistance depend on the schedule and sequence in which the drugs are given. Further, we tracked the clonal dynamics and mechanisms of resistance using DNA-integrated barcodes and single-cell RNA sequencing. Our study suggests that longer formats of treatment schedules in vitro screening assays are required to understand the effects of resistance and guide more realistically in vivo and clinical studies.
| S-EPMC8154902 | 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:BackgroundCell biology research is fundamentally limited by the number of intracellular components, particularly proteins, that can be co-measured in the same cell. Therefore, cell-to-cell heterogeneity in unmeasured proteins can lead to completely different observed relations between the same measured proteins. Attempts to infer such relations in a heterogeneous cell population can yield uninformative average relations if only one underlying biochemical network is assumed. To address this, we developed a method that recursively couples an iterative unmixing process with a Bayesian analysis of each unmixed subpopulation.ResultsOur approach enables to identify the number of distinct cell subpopulations, unmix their corresponding observations and resolve the network structure of each subpopulation. Using simulations of the MAPK pathway upon EGF and NGF stimulations we assess the performance of the method. We demonstrate that the presented method can identify better than clustering approaches the number of subpopulations within a mixture of observations, thus resolving correctly the statistical relations between the proteins.ConclusionsCoupling the unmixing of multiplexed observations with the inference of statistical relations between the measured parameters is essential for the success of both of these processes. Here we present a conceptual and algorithmic solution to achieve such coupling and hence to analyze data obtained from a natural mixture of cell populations. As the technologies and necessity for multiplexed measurements are rising in the systems biology era, this work addresses an important current challenge in the analysis of the derived data.
| S-EPMC4480582 | biostudies-literature
Project description:BackgroundMany biological processes, such as cancer metastasis, organismal development, and acquisition of resistance to cytotoxic therapy, rely on the emergence of rare sub-clones from a larger population. Understanding how the genetic and epigenetic features of diverse clones affect clonal fitness provides insight into molecular mechanisms underlying selective processes. While large-scale barcoding with NGS readout has facilitated cellular fitness assessment at the population level, this approach does not support characterization of clones prior to selection. Single-cell genomics methods provide high biological resolution, but are challenging to scale across large populations to probe rare clones and are destructive, limiting further functional analysis of important clones.ResultsHere, we develop CloneSifter, a methodology for tracking and enriching rare clones throughout their response to selection. CloneSifter utilizes a CRISPR sgRNA-barcode library that facilitates the isolation of viable cells from specific clones within the barcoded population using a sequence-specific retrieval reporter. We demonstrate that CloneSifter can measure clonal fitness of cancer cell models in vitro and retrieve targeted clones at abundance as low as 1 in 1883 in a heterogeneous cell population.ConclusionsCloneSifter provides a means to track and access specific and rare clones of interest across dynamic changes in population structure to comprehensively explore the basis of these changes.
| S-EPMC7687773 | biostudies-literature
Project description:Recent drug discovery and development efforts have created a large arsenal of targeted and chemotherapeutic drugs for precision medicine. However, drug resistance remains a major challenge as minor pre-existing resistant subpopulations are often found to be enriched at relapse. Current drug design has been heavily focused on initial efficacy, and we do not fully understand the effects of drug selective pressure on long-term drug resistance potential. Using a minimal two-population model, taking into account subpopulation proportions and growth/kill rates, we modeled long-term drug treatment and performed parameter sweeps to analyze the effects of each parameter on therapeutic efficacy. We found that drugs with the same overall initial kill may exert differential selective pressures, affecting long-term therapeutic outcome. We validated our conclusions experimentally using a preclinical model of Burkitt's lymphoma. Furthermore, we highlighted an intrinsic tradeoff between drug-imposed overall selective pressure and rate of adaptation. A principled approach in understanding the effects of distinct drug selective pressures on short-term and long-term tumor response enables better design of therapeutics that ultimately minimize relapse.
| S-EPMC5098152 | biostudies-literature
Project description:BackgroundOur group studies the interactions between cells of the brain and the neurotropic parasite Toxoplasma gondii. Using an in vivo system that allows us to permanently mark and identify brain cells injected with Toxoplasma protein, we have identified that Toxoplasma-injected neurons (TINs) are heterogeneously distributed throughout the brain. Unfortunately, standard methods to quantify and map heterogeneous cell populations onto a reference brain atlas are time consuming and prone to user bias.New methodWe developed a novel MATLAB-based semi-automated quantification and mapping program to allow the rapid and consistent mapping of heterogeneously distributed cells on to the Allen Institute Mouse Brain Atlas. The system uses two-threshold background subtraction to identify and quantify cells of interest.ResultsWe demonstrate that we reliably quantify and neuroanatomically localize TINs with low intra- or inter-observer variability. In a follow up experiment, we show that specific regions of the mouse brain are enriched with TINs.Comparison with existing methodsThe procedure we use takes advantage of simple immunohistochemistry labeling techniques, use of a standard microscope with a motorized stage, and low cost computing that can be readily obtained at a research institute. To our knowledge there is no other program that uses such readily available techniques and equipment for mapping heterogeneous populations of cells across the whole mouse brain.ConclusionThe quantification method described here allows reliable visualization, quantification, and mapping of heterogeneous cell populations in immunolabeled sections across whole mouse brains.
| S-EPMC6357967 | biostudies-literature
Project description:Tumor cell migration toward and intravasation into capillaries is an early and key event in cancer metastasis, yet not all cancer cells are imbued with the same capability to do so. This heterogeneity within a tumor is a fundamental property of cancer. Tools to help us understand what molecular characteristics allow a certain subpopulation of cells to spread from the primary tumor are thus critical for overcoming metastasis. Conventional in vitro migration platforms treat populations in aggregate, which leads to a masking of intrinsic differences among cells. Some migration assays reported recently have single-cell resolution, but these platforms do not provide for selective retrieval of the distinct migrating and non-migrating cell populations for further analysis. Thus, to study the intrinsic differences in cells responsible for chemotactic heterogeneity, we developed a single-cell migration platform so that individual cells' migration behavior can be studied and the heterogeneous population sorted based upon chemotactic phenotype. Furthermore, after migration, the highly chemotactic and non-chemotactic cells were retrieved and proved viable for later molecular analysis of their differences. Moreover, we modified the migration channel to resemble lymphatic capillaries to better understand how certain cancer cells are able to move through geometrically confining spaces.
| S-EPMC4435023 | biostudies-literature
Project description:Calcium responses have been observed as spikes of the whole-cell calcium concentration in numerous cell types and are essential for translating extracellular stimuli into cellular responses. While there are several suggestions for how this encoding is achieved, we still lack a comprehensive theory. To achieve this goal it is necessary to reliably predict the temporal evolution of calcium spike sequences for a given stimulus. Here, we propose a modelling framework that allows us to quantitatively describe the timing of calcium spikes. Using a Bayesian approach, we show that Gaussian processes model calcium spike rates with high fidelity and perform better than standard tools such as peri-stimulus time histograms and kernel smoothing. We employ our modelling concept to analyse calcium spike sequences from dynamically-stimulated HEK293T cells. Under these conditions, different cells often experience diverse stimulus time courses, which is a situation likely to occur in vivo. This single cell variability and the concomitant small number of calcium spikes per cell pose a significant modelling challenge, but we demonstrate that Gaussian processes can successfully describe calcium spike rates in these circumstances. Our results therefore pave the way towards a statistical description of heterogeneous calcium oscillations in a dynamic environment.
| S-EPMC5646906 | biostudies-literature
Project description:Models of neocortical networks are increasingly including the diversity of excitatory and inhibitory neuronal classes. Significant variability in cellular properties are also seen within a nominal neuronal class and this heterogeneity can be expected to influence the population response and information processing in networks. Recent studies have examined the population and network effects of variability in a particular neuronal parameter with some plausibly chosen distribution. However, the empirical variability and covariance seen across multiple parameters are rarely included, partly due to the lack of data on parameter correlations in forms convenient for model construction. To addess this we quantify the heterogeneity within and between the neocortical pyramidal-cell classes in layers 2/3, 4, and the slender-tufted and thick-tufted pyramidal cells of layer 5 using a combination of intracellular recordings, single-neuron modelling and statistical analyses. From the response to both square-pulse and naturalistic fluctuating stimuli, we examined the class-dependent variance and covariance of electrophysiological parameters and identify the role of the h current in generating parameter correlations. A byproduct of the dynamic I-V method we employed is the straightforward extraction of reduced neuron models from experiment. Empirically these models took the refractory exponential integrate-and-fire form and provide an accurate fit to the perisomatic voltage responses of the diverse pyramidal-cell populations when the class-dependent statistics of the model parameters were respected. By quantifying the parameter statistics we obtained an algorithm which generates populations of model neurons, for each of the four pyramidal-cell classes, that adhere to experimentally observed marginal distributions and parameter correlations. As well as providing this tool, which we hope will be of use for exploring the effects of heterogeneity in neocortical networks, we also provide the code for the dynamic I-V method and make the full electrophysiological data set available.
| S-EPMC4546387 | biostudies-literature