Cancer stem cells are controlled by developmental networks that are often topologically indistinguishable from normal, healthy stem cells. The question is why cancer stem cells can be both phenotypically distinct and have morphological effects so different from normal stem cells. The difference between cancer stem cells and normal stem cells lies not in differences their network architecture, but rather in the spatial-temporal locality of their activation in the genome and the resulting expression in the body. The metastatic potential cancer stem cells is not based primarily on their network divergence from normal stem cells, but on non-network based genetic changes that enable the evolution of gene-based phenotypic properties of the cell that permit its escape and travel to other parts of the body. Stem cell network theory allows the precise prediction of stem cell behavioral dynamics and a mathematical description of stem cell proliferation for both normal and cancer stem cells. It indicates that the best therapeutic approach is to tackle the highest order stem cells first, otherwise spontaneous remission of so called cured cancers will always be a danger. Stem cell networks poin
Cancer is increasingly perceived as a systems-level, network phenomenon. The major trend of malignant transformation can be described as a two-phase process, where an initial increase of network plasticity is followed by a decrease of plasticity at late stages of tumor development. The fluctuating intensity of stress factors, like hypoxia, inflammation and the either cooperative or hostile interactions of tumor inter-cellular networks, all increase the adaptation potential of cancer cells. This may lead to the bypass of cellular senescence, and to the development of cancer stem cells. We propose that the central tenet of cancer stem cell definition lies exactly in the indefinability of cancer stem cells. Actual properties of cancer stem cells depend on the individual "stress-history" of the given tumor. Cancer stem cells are characterized by an extremely large evolvability (i.e. a capacity to generate heritable phenotypic variation), which corresponds well with the defining hallmarks of cancer stem cells: the possession of the capacity to self-renew and to repeatedly re-build the heterogeneous lineages of cancer cells that comprise a tumor in new environments. Cancer stem cells rep
The homeostasis of epithelial tissue relies on a balance between the self-renewal of stem cell populations, cellular differentiation, and loss. Although this balance needs to be tightly regulated to avoid pathologies, such as tumor growth, the regulatory mechanisms, both cell-intrinsic and collective, which ensure tissue steady-state are still poorly understood. Here, we develop a computational model that incorporates basic assumptions of stem cell renewal into distinct populations and mechanical interactions between cells. We find that the model generates unexpected dynamic features: stem cells repel each other in the bulk tissue and are thus found rather isolated, as in a number of in vivo contexts. By mapping the system onto a gas of passive Brownian particles with effective repulsive interactions, that arise from the generated flows of differentiated cells, we show that we can quantitatively describe such stem cell distribution in tissues. The interaction potential between a pair of stem cells decays exponentially with a characteristic length that spans several cell sizes, corresponding to the volume of cells generated per stem cell division. Our findings may help understanding
In previous theoretical research, we inferred that cancer stem cells (CSCs), the cells that presumably drive tumor growth and resistance to conventional cancer treatments, are not uniformly distributed in the bulk of a tumorsphere. To confirm this theoretical prediction, we cultivated tumorspheres enriched in CSCs, and performed immunofluorecent detection of the stemness marker SOX2 using a confocal microscope. In this article, we present a method developed to process the images that reconstruct the amount and location of the CSCs in the tumorspheres. Its advantage is the use of a statistical criterion to classify the cells in stem and differentiated instead of setting an arbitrary threshold. From the analysis of the results of the methods using graph theory and computational modeling, we concluded that the distribution of Cancer Stem Cells in an experimental tumorsphere is non-homogeneous. This method is independent of the tumorsphere assay being useful for analyzing images in which several different kinds of cells are stained with different markers.
Stem cells are characterized by their ability to self-renew, as well as to differentiate and give rise to new populations of cells. Stem cell divisions are crucial for generative processes that occur during early development, and later in adulthood to support tissue regenerative capabilities. This property of stemness, the ability of self-renewal or tissue-specific differentiation, is also observed in cancer cells facilitating the sustenance of tumor growth, and in bipotent megakaryocytic-erythroid progenitors (MEPs) to produce blood cells. We are interested in modeling the size of the stem cell population required to adequately generate tissues or colonies of cells. We develop a state model that characterizes stem cell divisions and the dynamic changes of the stem cell and differentiated cell populations. In our model, the probabilities of self-renewal and differentiation events that stem cells undergo can vary over time instead of remaining constant throughout the process. We provide an estimation method for the division probabilities and using a simulation study, we show that our method provides good estimates even with a small sample size.
In addition to chemical and mechanical interactions between cells electromagnetic field produced by cells has been considered as another form of signaling for cell-cell communication. The aim of this study is evaluation of electromagnetic effects on viability of Adipose-derived stem cells (ADSCs) without co-culturing. In this study, stem cells were isolated from human adipose tissue enzymatically and proliferated in monolayer culture. Then, 5.(10^4) adipose-derived stem cells were cultured in each well of the test plate. In the first row (4 wells), ADSCs as inducer cells were cultured in DMEM1 with 10 ng/ml Fibroblast growth factor (FGF). In adjacent and the last rows, ADSCs were cultured without FGF (as detector cells). After the three and five days the viability of cells were evaluated. Moreover, ADSCs were cultured in the same conditions but the inducer cells were placed once in the UV-filter tube and once in the quartz tube to see whether there is electromagnetic interaction among cells. Inducer cells caused significant cell proliferation in adjacent row cells (p- value<0.01) in the fifth day. However, using the UV-filter tube and quartz tube both reduced the effect of induc
Human embryonic stem cells (ESC) are undifferentiated and are endowed with the capacities of self renewal and pluripotential differentiation. Adult stem cells renew their own tissue, but whether they can trans-differentiate to other tissues is still controversial. To understand the genetic program that underlies the pluripotency of stem cells, we compared the transcription profile of ESC with that of progenitor/stem cells of human hematopoietic and keratinocytic origins, along with their mature cells to be viewed as snapshots along tissue differentiation. ESC gene profile show higher complexity with significantly more highly expressed genes than adult cells. We hypothesize that ESC use a strategy of expressing genes that represent various differentiation pathways and selection of only a few for continuous expression upon differentiation to a particular target. Such a strategy may be necessary for the pluripotency of ESC. The progenitors of either hematopoietic or keratinocytic cells also follow the same design principle. Using advanced clustering, we show that many of the ESC expressed genes are turned off in the progenitors/stem cells followed by a further downregulation in adult
Cycling tissues such as the intestinal epithelium, germ line, and hair follicles, require a constant flux of differentiated cells. These tissues are maintained by a population of stem cells, which generate differentiated progenies and self-renew. Asymmetric division of each stem cell into one stem cell and one differentiated cell can accomplish both tasks. However, in mammalian cycling tissues, some stem cells divide symmetrically into two differentiated cells and are replaced by a neighbor that divides symmetrically into two stem cells. Besides this heterogeneity in fate (population asymmetry), stem cells also exhibit heterogenous proliferation-rates; in the long run, however, all stem cells proliferate at the same average rate (equipotency). We construct and simulate a mathematical model based on these experimental observations. We show that the complex steady-state dynamics of population-asymmetric stem cells reduces the rate of replicative aging of the tissue --potentially lowering the incidence of somatic mutations and genetics diseases such as cancer. Essentially, slow-dividing stem cells proliferate and purge the population of the fast-dividing --older-- cells which had unde
Embryonic stem cells (ESCs) and induced pluripotent stem cells (iPSCs) derived from somatic cells (SCs) provide promising resources for regenerative medicine and medical research, leading to a daily identification of new cell lines. However, an efficient system to discriminate the cell lines is lacking. Here, we developed a quantitative system to discriminate the three cell types, iPSCs, ESCs and SCs. The system contains DNA-methylation biomarkers and mathematical models, including an artificial neural network and support vector machines. All biomarkers were unbiasedly selected by calculating an eigengene score derived from analysis of genome-wide DNA methylations. With 30 biomarkers, or even with as few as 3 top biomarkers, this system can discriminate SCs from ESCs and iPSCs with almost 100% accuracy, and with approximately 100 biomarkers, the system can distinguish ESCs from iPSCs with an accuracy of 95%. This robust system performs precisely with raw data without normalization as well as with converted data in which the continuous methylation levels are accounted. Strikingly, this system can even accurately predict new samples generated from different microarray platforms and t
Stem cells have the potential to produce lineages of non-stem cell populations (differentiated cells) via a ubiquitous hierarchal division scheme. Differentiation of a stem cell into (partially) differentiated cells can happen either symmetrically or asymmetrically. The selection dynamics of a mutant cancer stem cell should be investigated in the light of a stem cell proliferation hierarchy and presence of a non-stem cell population. By constructing a three-compartment Moran-type model composed of normal stem cells, mutant (cancer) stem cells and differentiated cells, we derive the replicator dynamics of stem cell frequencies where asymmetric differentiation and differentiated cell death rates are included in the model. We determine how these new factors change the conditions for a successful mutant invasion and discuss the variation on the steady state fraction of the population as different model parameters are changed. By including the phenotypic plasticity/dedifferentiation, in which a progenitor/differentiated cell can transform back into a cancer stem cell, we show that the effective fitness of mutant stem cells is not only determined by their proliferation and death rates bu
Rapidly dividing tissues, like intestinal crypts, are frequently chosen to investigate the process of tumor initiation, because of their high rate of mutations. To study the interplay between normal and mutant as well as immortal cells in the human colon or intestinal crypt, we developed a 4-compartmental stochastic model for cell dynamics based on current discoveries. Recent studies of the intestinal crypt have revealed the existence of two stem cell groups. Therefore, our model incorporates two stem cell groups (central stem cells (CeSCs) and border stem cells (BSCs)), plus one compartment for transit amplifying (TA) cells and one compartment of fully differentiated (FD) cells. However, it can be easily modified to have only one stem cell group. We find that the worst-case scenario occurs when CeSCs are mutated, or an immortal cell arises in the TA or FD compartments. The probability that the progeny of a single advantageous CeSC mutant will take over the entire crypt is more than $0.2$, and one immortal cell always causes all FD cells to become immortals.Moreover, when CeSCs are either mutants or wild-type (w.t.) individuals, their progeny will take over the entire crypt in less
In this review, we discuss recent advances on the plasticity of cancer stem cells and highlight their relevance to understand the metastatic process and to guide therapeutic interventions. Recent results suggest that the strict hierarchical structure of cancer cell populations advocated by the cancer stem cell model must be reconsidered since the depletion of cancer stem cells leads the other tumor cells to switch back into the cancer stem cell phenotype. This plasticity has important implications for metastasis since migrating cells do not need to be cancer stem cells in order to seed a metastasis. We also discuss the important role of the immune system and the microenvironment in modulating phenotypic switching and suggest possible avenues to exploit our understanding of this process to develop an effective strategy for precision medicine.
Modern tissue engineering strategies combine living cells and scaffold materials to develop biological substitutes that can restore tissue functions. Both natural and synthetic materials have been fabricated for transplantation of stem cells and their specific differentiation into muscles, bones and cartilages. One of the key objectives for bone regeneration therapy to be successful is to direct stem cells' proliferation and to accelerate their differentiation in a controlled manner through the use of growth factors and osteogenic inducers. Here we show that graphene provides a promising biocompatible scaffold that does not hamper the proliferation of human mesenchymal stem cells (hMSCs) and accelerates their specific differentiation into bone cells. The differentiation rate is comparable to the one achieved with common growth factors, demonstrating graphene's potential for stem cell research.
Cancer stem cells have been shown to be critical to the development of a variety of solid cancers. The precise interplay mechanisms between cancer stem cells and the rest of a tissue are still not elucidated. To shed light on the interactions between stem and non-stem cancer cell populations we develop a two-population mathematical model, which is suitable to describe tumorsphere growth. Both interspecific and intraspecific interactions, mediated by the microenvironment, are included. We show that there is a tipping point, characterized by a transcritical bifurcation, where a purely non-stem cell attractor is replaced by a new attractor that contains both stem and differentiated cancer cells. The model is then applied to describe the outcome of a recent experiment. This description reveals that, while the intraspecific interactions are inhibitory, the interspecific interactions stimulate growth. This can be understood in terms of stem cells needing differentiated cells to reinforce their niches, and phenotypic plasticity favoring the de-differentiation of differentiated cells into cancer stem cells. We posit that this is a consequence of the deregulation of the quorum sensing that
Tumors are defined by their intense proliferation, but sometimes cancer cells turn senescent and stop replicating. In the stochastic cancer model in which all cells are tumorigenic, senescence is seen as the result of random mutations, suggesting that it could represent a barrier to tumor growth. In the hierarchical cancer model a subset of the cells, the cancer stem cells, divide indefinitely while other cells eventually turn senescent. Here we formulate cancer growth in mathematical terms and obtain predictions for the evolution of senescence. We perform experiments in human melanoma cells which are compatible with the hierarchical model and show that senescence is a reversible process controlled by survivin. We conclude that enhancing senescence is unlikely to provide a useful therapeutic strategy to fight cancer, unless the cancer stem cells are specifically targeted
Since the discovery of a cancer initiating side population in solid tumours, studies focussing on the role of so-called cancer stem cells in cancer initiation and progression have abounded. The biological interrogation of these cells has yielded volumes of information about their behaviour, but there has, as of yet, not been many actionable generalised theoretical conclusions. To address this point, we have created a hybrid, discrete/continuous computational cellular automaton model of a generalised stem-cell driven tissue and explored the phenotypic traits inherent in the inciting cell and the resultant tissue growth. We identify the regions in phenotype parameter space where these initiating cells are able to cause a disruption in homeostasis, leading to tissue overgrowth and tumour formation. As our parameters and model are non-specific, they could apply to any tissue cancer stem-cell and do not assume specific genetic mutations. In this way, our model suggests that targeting these phenotypic traits could represent generalizable strategies across cancer types and represents a first attempt to identify the hallmarks of cancer stem cells.
Accumulating evidence suggests that many tumors have a hierarchical organization, with the bulk of the tumor composed of relatively differentiated short-lived progenitor cells that are maintained by a small population of undifferentiated long-lived cancer stem cells. It is unclear, however, whether cancer stem cells originate from normal stem cells or from dedifferentiated progenitor cells. To address this, we mathematically modeled the effect of dedifferentiation on carcinogenesis. We considered a hybrid stochastic-deterministic model of mutation accumulation in both stem cells and progenitors, including dedifferentiation of progenitor cells to a stem cell-like state. We performed exact computer simulations of the emergence of tumor subpopulations with two mutations, and we derived semi-analytical estimates for the waiting time distribution to fixation. Our results suggest that dedifferentiation may play an important role in carcinogenesis, depending on how stem cell homeostasis is maintained. If the stem cell population size is held strictly constant (due to all divisions being asymmetric), we found that dedifferentiation acts like a positive selective force in the stem cell popu
The identification of cancer stem cells in vivo and in vitro relies on specific surface markers that should allow to sort cancer cells in phenotypically distinct subpopulations. Experiments report that sorted cancer cell populations after some time tend to express again all the original markers, leading to the hypothesis of phenotypic switching, according to which cancer cells can transform stochastically into cancer stem cells. Here we explore an alternative explanation based on the hypothesis that markers are not perfect and are thus unable to identify all cancer stem cells. Our analysis is based on a mathematical model for cancer cell proliferation that takes into account phenotypic switching, imperfect markers and error in the sorting process. Our conclusion is that the observation of reversible expression of surface markers after sorting does not provide sufficient evidence in support of phenotypic switching.
MicroRNAs (miRNAs) critically modulate stem cell properties like pluripotency, but the fundamental mechanism remains largely unknown. This study systematically analyzes multiple-omics data and builds a systems physical network including genome-wide interactions between miRNAs and their targets to reveal the systems mechanism of miRNA functions in mouse pluripotent stem cells. Globally, miRNAs directly repress the pluripotent core factors during differentiation state. Surprisingly, during pluripotent state, the top important miRNAs do not directly regulate the pluripotent core factors as thought, but they only directly target the pluripotent signal pathways and directly repress developmental processes. Furthermore, at pluripotent state miRNAs predominately repress DNA methyltransferases, the core enzymes for DNA methylation. The decreasing methylation repressed by miRNAs in turn activates the top miRNAs and pluripotent core factors, creating an active circuit system to modulate pluripotency. MiRNAs vary their functions with different stem cell states. While miRNAs directly repress pluripotent core factors to facilitate the differentiation during cell differentiation, they also help
Statistical and mathematical modeling are crucial to describe, interpret, compare and predict the behavior of complex biological systems including the organization of hematopoietic stem and progenitor cells in the bone marrow environment. The current prominence of high-resolution and live-cell imaging data provides an unprecedented opportunity to study the spatiotemporal dynamics of these cells within their stem cell niche and learn more about aberrant, but also unperturbed, normal hematopoiesis. However, this requires careful quantitative statistical analysis of the spatial and temporal behavior of cells and the interaction with their microenvironment. Moreover, such quantification is a prerequisite for the construction of hypothesis-driven mathematical models that can provide mechanistic explanations by generating spatiotemporal dynamics that can be directly compared to experimental observations. Here, we provide a brief overview of statistical methods in analyzing spatial distribution of cells, cell motility, cell shapes and cellular genealogies. We also describe cell- based modeling formalisms that allow researchers to simulate emergent behavior in a multicellular system based