Spatial pattern dynamics of 3D stem cell loss of pluripotency via rules-based computational modeling

Recognition and reconstruction of cell differentiation patterns with deep learning

Cell lineage decisions occur in three-dimensional spatial patterns that are difficult to identify by eye. There is an ongoing effort to replicate such patterns using mathematical modeling. One approach uses long ranging cell-cell communication to replicate common spatial arrangements like checkerboard and engulfing patterns. In this model, the cell-cell communication has been implemented as a signal that disperses throughout the tissue. On the other hand, machine learning models have been developed for pattern recognition and pattern reconstruction tasks. We combined synthetic data generated by the mathematical model with spatial summary statistics and deep learning algorithms to recognize and reconstruct cell fate patterns in organoids of mouse embryonic stem cells. Application of Moran's index and pair correlation functions for in vitro and synthetic data from the model showed local clustering and radial segregation. To assess the patterns as a whole, a graph neural network was developed and trained on synthetic data from the model. Application to in vitro data predicted a low signal dispersion value. To test this result, we implemented a multilayer perceptron for the prediction of a given cell fate based on the fates of the neighboring cells. The results show a 70% accuracy of cell fate imputation based on the nine nearest neighbors of a cell. Overall, our approach combines deep learning with mathematical modeling to link cell fate patterns with potential underlying mechanisms.

Adjusting the range of cell-cell communication enables fine-tuning of cell fate patterns from checkerboard to engulfing

During development, spatio-temporal patterns ranging from checkerboard to engulfing occur with precise proportions of the respective cell fates. Key developmental regulators are intracellular transcriptional interactions and intercellular signaling. We present an analytically tractable mathematical model based on signaling that reliably generates different cell type patterns with specified proportions. Employing statistical mechanics, We derived a cell fate decision model for two cell types. A detailed steady state analysis on the resulting dynamical system yielded necessary conditions to generate spatially heterogeneous patterns. This allows the cell type proportions to be controlled by a single model parameter. Cell-cell communication is realized by local and global signaling mechanisms. These result in different cell type patterns. A nearest neighbor signal yields checkerboard patterns. Increasing the signal dispersion, cell fate clusters and an engulfing pattern can be generated. Altogether, the presented model allows us to reliably generate heterogeneous cell type patterns of different kinds as well as desired proportions.

Engineering multicellular living systems-a Keystone Symposia report

The ability to engineer complex multicellular systems has enormous potential to inform our understanding of biological processes and disease and alter the drug development process. Engineering living systems to emulate natural processes or to incorporate new functions relies on a detailed understanding of the biochemical, mechanical, and other cues between cells and between cells and their environment that result in the coordinated action of multicellular systems. On April 3-6, 2022, experts in the field met at the Keystone symposium "Engineering Multicellular Living Systems" to discuss recent advances in understanding how cells cooperate within a multicellular system, as well as recent efforts to engineer systems like organ-on-a-chip models, biological robots, and organoids. Given the similarities and common themes, this meeting was held in conjunction with the symposium "Organoids as Tools for Fundamental Discovery and Translation".

Community effects allow bioelectrical reprogramming of cell membrane potentials in multicellular aggregates: Model simulations

Bioelectrical patterns are established by spatiotemporal correlations of cell membrane potentials at the multicellular level, being crucial to development, regeneration, and tumorigenesis. We have conducted multicellular simulations on bioelectrical community effects and intercellular coupling in multicellular aggregates. The simulations aim at establishing under which conditions a local heterogeneity consisting of a small patch of cells can be stabilized against a large aggregate of surrounding identical cells which are in a different bioelectrical state. In this way, instructive bioelectrical information can be persistently encoded in spatiotemporal patterns of separated domains with different cell polarization states. The multicellular community effects obtained are regulated both at the single-cell and intercellular levels, and emerge from a delicate balance between the degrees of intercellular coupling in: (i) the small patch, (ii) the surrounding bulk, and (iii) the interface that separates these two regions. The model is experimentally motivated and consists of two generic voltage-gated ion channels that attempt to establish the depolarized and polarized cell states together with coupling conductances whose individual and intercellular different states permit a dynamic multicellular connectivity. The simulations suggest that community effects may allow the reprogramming of single-cell bioelectrical states, in agreement with recent experimental data. A better understanding of the resulting electrical regionalization can assist the electroceutical correction of abnormally depolarized regions initiated in the bulk of normal tissues as well as suggest new biophysical mechanisms for the establishment of target patterns in multicellular engineering.

Deep neural net tracking of human pluripotent stem cells reveals intrinsic behaviors directing morphogenesis

Lineage tracing is a powerful tool in developmental biology to interrogate the evolution of tissue formation, but the dense, three-dimensional nature of tissue limits the assembly of individual cell trajectories into complete reconstructions of development. Human induced pluripotent stem cells (hiPSCs) can recapitulate aspects of developmental processes, providing an in vitro platform to assess the dynamic collective behaviors directing tissue morphogenesis. Here, we trained an ensemble of neural networks to track individual hiPSCs in time-lapse microscopy, generating longitudinal measures of cell and cellular neighborhood properties on timescales from minutes to days. Our analysis reveals that, while individual cell parameters are not strongly affected by pluripotency maintenance conditions or morphogenic cues, regional changes in cell behavior predict cell fate and colony organization. By generating complete multicellular reconstructions of hiPSC behavior, our tracking pipeline enables fine-grained understanding of morphogenesis by elucidating the role of regional behavior in early tissue formation.

Automated Design of Pluripotent Stem Cell Self-Organization

Human pluripotent stem cells (hPSCs) have the intrinsic ability to self-organize into complex multicellular organoids that recapitulate many aspects of tissue development. However, robustly directing morphogenesis of hPSC-derived organoids requires novel approaches to accurately control self-directed pattern formation. Here, we combined genetic engineering with computational modeling, machine learning, and mathematical pattern optimization to create a data-driven approach to control hPSC self-organization by knock down of genes previously shown to affect stem cell colony organization, CDH1 and ROCK1. Computational replication of the in vitro system in silico using an extended cellular Potts model enabled machine learning-driven optimization of parameters that yielded emergence of desired patterns. Furthermore, in vitro the predicted experimental parameters quantitatively recapitulated the in silico patterns. These results demonstrate that morphogenic dynamics can be accurately predicted through model-driven exploration of hPSC behaviors via machine learning, thereby enabling spatial control of multicellular patterning to engineer human organoids and tissues. A record of this paper's Transparent Peer Review process is included in the Supplemental Information.

Agent-based modeling of morphogenetic systems: Advantages and challenges

The complexity of morphogenesis poses a fundamental challenge to understanding the mechanisms governing the formation of biological patterns and structures. Over the past century, numerous processes have been identified as critically contributing to morphogenetic events, but the interplay between the various components and aspects of pattern formation have been much harder to grasp. The combination of traditional biology with mathematical and computational methods has had a profound effect on our current understanding of morphogenesis and led to significant insights and advancements in the field. In particular, the theoretical concepts of reaction–diffusion systems and positional information, proposed by Alan Turing and Lewis Wolpert, respectively, dramatically influenced our general view of morphogenesis, although typically in isolation from one another. In recent years, agent-based modeling has been emerging as a consolidation and implementation of the two theories within a single framework. Agent-based models (ABMs) are unique in their ability to integrate combinations of heterogeneous processes and investigate their respective dynamics, especially in the context of spatial phenomena. In this review, we highlight the benefits and technical challenges associated with ABMs as tools for examining morphogenetic events. These models display unparalleled flexibility for studying various morphogenetic phenomena at multiple levels and have the important advantage of informing future experimental work, including the targeted engineering of tissues and organs.

Mouse ICM Organoids Reveal Three-Dimensional Cell Fate Clustering

During mammalian preimplantation, cells of the inner cell mass (ICM) adopt either an embryonic or an extraembryonic fate. This process is tightly regulated in space and time and has been studied previously in mouse embryos and embryonic stem cell models. Current research suggests that cell fates are arranged in a salt-and-pepper pattern of random cell positioning or a spatially alternating pattern. However, the details of the three-dimensional patterns of cell fate specification have not been investigated in the embryo nor in in vitro systems. We developed ICM organoids as a, to our knowledge, novel three-dimensional in vitro stem cell system to model mechanisms of fate decisions that occur in the ICM. ICM organoids show similarities to the in vivo system that arise regardless of the differences in geometry and total cell number. Inspecting ICM organoids and mouse embryos, we describe a so far unknown local clustering of cells with identical fates in both systems. These findings are based on the three-dimensional quantitative analysis of spatiotemporal patterns of NANOG and GATA6 expression in combination with computational rule-based modeling. The pattern identified by our analysis is distinct from the current view of a salt-and-pepper pattern. Our investigation of the spatial distributions both in vivo and in vitro dissects the contributions of the different parts of the embryo to cell fate specifications. In perspective, our combination of quantitative in vivo and in vitro analyses can be extended to other mammalian organisms and thus creates a powerful approach to study embryogenesis.

Dynamic intercellular transport modulates the spatial patterning of differentiation during early neural commitment

The initiation of heterogeneity within a population of phenotypically identical progenitors is a critical event for the onset of morphogenesis and differentiation patterning. Gap junction communication within multicellular systems produces complex networks of intercellular connectivity that result in heterogeneous distributions of intracellular signaling molecules. In this study, we investigate emergent systems-level behavior of the intercellular network within embryonic stem cell (ESC) populations and corresponding spatial organization during early neural differentiation. An agent-based model incorporates experimentally-determined parameters to yield complex transport networks for delivery of pro-differentiation cues between neighboring cells, reproducing the morphogenic trajectories during retinoic acid-accelerated mouse ESC differentiation. Furthermore, the model correctly predicts the delayed differentiation and preserved spatial features of the morphogenic trajectory that occurs in response to intercellular perturbation. These findings suggest an integral role of gap junction communication in the temporal coordination of emergent patterning during early differentiation and neural commitment of pluripotent stem cells.

Probing the stochastic property of endoreduplication in cell size determination of Arabidopsis thaliana leaf epidermal tissue

Cell size distribution is highly reproducible, whereas the size of individual cells often varies greatly within a tissue. This is obvious in a population of Arabidopsis thaliana leaf epidermal cells, which ranged from 1,000 to 10,000 μm² in size. Endoreduplication is a specialized cell cycle in which nuclear genome size (ploidy) is doubled in the absence of cell division. Although epidermal cells require endoreduplication to enhance cellular expansion, the issue of whether this mechanism is sufficient for explaining cell size distribution remains unclear due to a lack of quantitative understanding linking the occurrence of endoreduplication with cell size diversity. Here, we addressed this question by quantitatively summarizing ploidy profile and cell size distribution using a simple theoretical framework. We first found that endoreduplication dynamics is a Poisson process through cellular maturation. This finding allowed us to construct a mathematical model to predict the time evolution of a ploidy profile with a single rate constant for endoreduplication occurrence in a given time. We reproduced experimentally measured ploidy profile in both wild-type leaf tissue and endoreduplication-related mutants with this analytical solution, further demonstrating the probabilistic property of endoreduplication. We next extended the mathematical model by incorporating the element that cell size is determined according to ploidy level to examine cell size distribution. This analysis revealed that cell size is exponentially enlarged 1.5 times every endoreduplication round. Because this theoretical simulation successfully recapitulated experimentally observed cell size distributions, we concluded that Poissonian endoreduplication dynamics and exponential size-boosting are the sources of the broad cell size distribution in epidermal tissue. More generally, this study contributes to a quantitative understanding whereby stochastic dynamics generate steady-state biological heterogeneity.

Derivative-Free Optimization of Rate Parameters of Capsid Assembly Models from Bulk in Vitro Data

The assembly of virus capsids proceeds by a complicated cascade of association and dissociation steps, the great majority of which cannot be directly experimentally observed. This has made capsid assembly a rich field for computational models, but there are substantial obstacles to model inference for such systems. Here, we describe progress on fitting kinetic rate constants defining capsid assembly models to experimental data, a difficult data-fitting problem because of the high computational cost of simulating assembly trajectories, the stochastic noise inherent to the models, and the limited and noisy data available for fitting. We evaluate the merits of data-fitting methods based on derivative-free optimization (DFO) relative to gradient-based methods used in prior work. We further explore the advantages of alternative data sources through simulation of a model of time-resolved mass spectrometry data, a technology for monitoring bulk capsid assembly that can be expected to provide much richer data than previously used static light scattering approaches. The results show that advances in both the data and the algorithms can improve model inference. More informative data sources lead to high-quality fits for all methods, but DFO methods show substantial advantages on less informative data sources that better represent current experimental practice.

Mathematical Models of Pluripotent Stem Cells: At the Dawn of Predictive Regenerative Medicine

Regenerative medicine, ranging from stem cell therapy to organ regeneration, is promising to revolutionize treatments of diseases and aging. These approaches require a perfect understanding of cell reprogramming and differentiation. Predictive modeling of cellular systems has the potential to provide insights about the dynamics of cellular processes, and guide their control. Moreover in many cases, it provides alternative to experimental tests, difficult to perform for practical or ethical reasons. The variety and accuracy of biological processes represented in mathematical models grew in-line with the discovery of underlying molecular mechanisms. High-throughput data generation led to the development of models based on data analysis, as an alternative to more established modeling based on prior mechanistic knowledge. In this chapter, we give an overview of existing mathematical models of pluripotency and cell fate, to illustrate the variety of methods and questions. We conclude that current approaches are yet to overcome a number of limitations: Most of the computational models have so far focused solely on understanding the regulation of pluripotency, and the differentiation of selected cell lineages. In addition, models generally interrogate only a few biological processes. However, a better understanding of the reprogramming process leading to ESCs and iPSCs is required to improve stem-cell therapies. One also needs to understand the links between signaling, metabolism, regulation of gene expression, and the epigenetics machinery.

Quantitative multivariate analysis of dynamic multicellular morphogenic trajectories

Interrogating fundamental cell biology principles that govern tissue morphogenesis is critical to better understanding of developmental biology and engineering novel multicellular systems. Recently, functional micro-tissues derived from pluripotent embryonic stem cell (ESC) aggregates have provided novel platforms for experimental investigation; however elucidating the factors directing emergent spatial phenotypic patterns remains a significant challenge. Computational modelling techniques offer a unique complementary approach to probe mechanisms regulating morphogenic processes and provide a wealth of spatio-temporal data, but quantitative analysis of simulations and comparison to experimental data is extremely difficult. Quantitative descriptions of spatial phenomena across multiple systems and scales would enable unprecedented comparisons of computational simulations with experimental systems, thereby leveraging the inherent power of computational methods to interrogate the mechanisms governing emergent properties of multicellular biology. To address these challenges, we developed a portable pattern recognition pipeline consisting of: the conversion of cellular images into networks, extraction of novel features via network analysis, and generation of morphogenic trajectories. This novel methodology enabled the quantitative description of morphogenic pattern trajectories that could be compared across diverse systems: computational modelling of multicellular structures, differentiation of stem cell aggregates, and gastrulation of cichlid fish. Moreover, this method identified novel spatio-temporal features associated with different stages of embryo gastrulation, and elucidated a complex paracrine mechanism capable of explaining spatiotemporal pattern kinetic differences in ESC aggregates of different sizes.

Inverse problem analysis of pluripotent stem cell aggregation dynamics in stirred-suspension cultures

The cultivation of stem cells as aggregates in scalable bioreactor cultures is an appealing modality for the large-scale manufacturing of stem cell products. Aggregation phenomena are central to such bioprocesses affecting the viability, proliferation and differentiation trajectory of stem cells but a quantitative framework is currently lacking. A population balance equation (PBE) model was used to describe the temporal evolution of the embryonic stem cell (ESC) cluster size distribution by considering collision-induced aggregation and cell proliferation in a stirred-suspension vessel. For ESC cultures at different agitation rates, the aggregation kernel representing the aggregation dynamics was successfully recovered as a solution of the inverse problem. The rate of change of the average aggregate size was greater at the intermediate rate tested suggesting a trade-off between increased collisions and agitation-induced shear. Results from forward simulation with obtained aggregation kernels were in agreement with transient aggregate size data from experiments. We conclude that the framework presented here can complement mechanistic studies offering insights into relevant stem cell clustering processes. More importantly from a process development standpoint, this strategy can be employed in the design and control of bioreactors for the generation of stem cell derivatives for drug screening, tissue engineering and regenerative medicine.

Computational modelling of embryonic stem-cell fate control

The maintenance of pluripotency in embryonic stem cells (ESCs), its loss during lineage specification or its re-induction to generate induced pluripotent stem cells are central topics in stem cell biology. To uncover the molecular basis and the design principles of pluripotency control, a multitude of experimental, but also an increasing number of computational, studies have been published. Here, we consider recent reports that apply computational or mathematical modelling approaches to describe the regulatory processes that underlie cell fate decisions in mouse ESCs. We summarise the principles, the strengths and potentials but also the limitations of different computational strategies.

An agent-based model of cancer stem cell initiated avascular tumour growth and metastasis: the effect of seeding frequency and location

It is very important to understand the onset and growth pattern of breast primary tumours as well as their metastatic dissemination. In most cases, it is the metastatic disease that ultimately kills the patient. There is increasing evidence that cancer stem cells are closely linked to the progression of the metastatic tumour. Here, we investigate stem cell seeding to an avascular tumour site using an agent-based stochastic model of breast cancer metastatic seeding. The model includes several important cellular features such as stem cell symmetric and asymmetric division, migration, cellular quiescence, senescence, apoptosis and cell division cycles. It also includes external features such as stem cell seeding frequency and location. Using this model, we find that cell seeding rate and location are important features for tumour growth. We also define conditions in which the tumour growth exhibits decremented and exponential growth patterns. Overall, we find that seeding, senescence and division limit affect not only the number of stem cells, but also their spatial and temporal distribution.

In vitro spatially organizing the differentiation in individual multicellular stem cell aggregates

With significant potential as a robust source to produce specific somatic cells for regenerative medicine, stem cells have attracted increasing attention from both academia and government. In vivo, stem cell differentiation is a process under complicated regulations to precisely build tissue with unique spatial structures. Since multicellular spheroidal aggregates of stem cells, commonly called as embryoid bodies (EBs), are considered to be capable of recapitulating the events in early stage of embryonic development, a variety of methods have been developed to form EBs in vitro for studying differentiation of embryonic stem cells. The regulation of stem cell differentiation is crucial in directing stem cells to build tissue with the correct spatial architecture for specific functions. However, stem cells within the three-dimensional multicellular aggregates undergo differentiation in a less unpredictable and spatially controlled manner in vitro than in vivo. Recently, various microengineering technologies have been developed to manipulate stem cells in vitro in a spatially controlled manner. Herein, we take the spotlight on these technologies and researches that bring us the new potential for manipulation of stem cells for specific purposes.

Stem cells in a three-dimensional scaffold environment

Stem cells have emerged as important players in the generation and maintenance of many tissues. However, the accurate in vitro simulation of the native stem cell niche remains difficult due at least in part to the lack of a comprehensive definition of the critical factors of the stem cell niche based on in vivo models. Three-dimensional (3D) cell culture systems have allowed the development of useful models for investigating stem cell physiology particularly with respect to their ability to sense and generate mechanical force in response to their surrounding environment. We review the use of 3D culture systems for stem cell culture and discuss the relationship between stem cells and 3D growth matrices including the roles of the extracellular matrix, scaffolds, soluble factors, cell-cell interactions and shear stress effects within this environment. We also discuss the potential for novel methods that mimic the native stem cell niche in vitro as well as the current associated challenges.

Single-cell analysis of embryoid body heterogeneity using microfluidic trapping array

The differentiation of pluripotent stem cells as embryoid bodies (EBs) remains a common method for inducing differentiation toward many lineages. However, differentiation via EBs typically yields a significant amount of heterogeneity in the cell population, as most cells differentiate simultaneously toward different lineages, while others remain undifferentiated. Moreover, physical parameters, such as the size of EBs, can modulate the heterogeneity of differentiated phenotypes due to the establishment of nutrient and oxygen gradients. One of the challenges in examining the cellular composition of EBs is the lack of analytical methods that are capable of determining the phenotype of all of the individual cells that comprise a single EB. Therefore, the objective of this work was to examine the ability of a microfluidic cell trapping array to analyze the heterogeneity of cells comprising EBs during the course of early differentiation. The heterogeneity of single cell phenotype on the basis of protein expression of the pluripotent transcription factor OCT-4 was examined for populations of EBs and single EBs of different sizes at distinct stages of differentiation. Results from the cell trap device were compared with flow cytometry and whole mount immunostaining. Additionally, single cells from dissociated pooled EBs or individual EBs were examined separately to discern potential differences in the value or variance of expression between the different methods of analysis. Overall, the analytical method described represents a novel approach for evaluating how heterogeneity is manifested in EB cultures and may be used in the future to assess the kinetics and patterns of differentiation in addition to the loss of pluripotency.

Engineering three-dimensional stem cell morphogenesis for the development of tissue models and scalable regenerative therapeutics

The physiochemical stem cell microenvironment regulates the delicate balance between self-renewal and differentiation. The three-dimensional assembly of stem cells facilitates cellular interactions that promote morphogenesis, analogous to the multicellular, heterotypic tissue organization that accompanies embryogenesis. Therefore, expansion and differentiation of stem cells as multicellular aggregates provides a controlled platform for studying the biological and engineering principles underlying spatiotemporal morphogenesis and tissue patterning. Moreover, three-dimensional stem cell cultures are amenable to translational screening applications and therapies, which underscores the broad utility of scalable suspension cultures across laboratory and clinical scales. In this review, we discuss stem cell morphogenesis in the context of fundamental biophysical principles, including the three-dimensional modulation of adhesions, mechanics, and molecular transport and highlight the opportunities to employ stem cell spheroids for tissue modeling, bioprocessing, and regenerative therapies.

The role of SDF-1α-ECM crosstalk in determining neural stem cell fate

The consequences of central nervous system injury are far-reaching and debilitating and, while an endogenous repair response to neural injury has been observed in recent years, the mechanisms behind this response remain unclear. Neural progenitor/stem cell (NPSC) migration to the site of injury from the neural stem cell niches (e.g. subventricular zone and hippocampus) has been observed to be vasophilic in nature. While the chemotactic stimuli directing NPSC homing to injury is not well established, it is thought to be due in part to an increasing gradient of chemotactic cytokines, such as stromal cell-derived factor 1α (SDF-1α). Based on these recent findings, we hypothesize that critical crosstalk between SDF-1α and the extracellular matrix (ECM) drives injury-induced NPSC behavior. In this study, we investigated the effect of SDF-1α and ECM substrates (Matrigel, laminin, and vitronectin) on the migration, differentiation, and proliferation of NPSCs in vitro using standard assays. The results demonstrated that SDF-1α and laminin-based ECM (Matrigel and laminin) significantly and synergistically enhanced NPSC migration and acute neuronal differentiation. These effects were significantly attenuated with the addition of AMD3100 (an antagonist against the SDF-1α receptor, CXCR4). SDF-1α alone significantly increased NPSC proliferation regardless of ECM substrate, however no synergy was observed between SDF-1α and the ECM. These results serve to elucidate the relationship between adhesive and soluble signaling factors of interest and their effect on NPSC behavior following neural injury. Furthermore, these results better inform the next generation of biomaterials aimed at stimulating endogenous neural regeneration for neural injury and neurodegenerative diseases.

Time scales in epigenetic dynamics and phenotypic heterogeneity of embryonic stem cells

A remarkable feature of the self-renewing population of embryonic stem cells (ESCs) is their phenotypic heterogeneity: Nanog and other marker proteins of ESCs show large cell-to-cell variation in their expression level, which should significantly influence the differentiation process of individual cells. The molecular mechanism and biological implication of this heterogeneity, however, still remain elusive. We address this problem by constructing a model of the core gene-network of mouse ESCs. The model takes account of processes of binding/unbinding of transcription factors, formation/dissolution of transcription apparatus, and modification of histone code at each locus of genes in the network. These processes are hierarchically interrelated to each other forming the dynamical feedback loops. By simulating stochastic dynamics of this model, we show that the phenotypic heterogeneity of ESCs can be explained when the chromatin at the Nanog locus undergoes the large scale reorganization in formation/dissolution of transcription apparatus, which should have the timescale similar to the cell cycle period. With this slow transcriptional switching of Nanog, the simulated ESCs fluctuate among multiple transient states, which can trigger the differentiation into the lineage-specific cell states. From the simulated transitions among cell states, the epigenetic landscape underlying transitions is calculated. The slow Nanog switching gives rise to the wide basin of ESC states in the landscape. The bimodal Nanog distribution arising from the kinetic flow running through this ESC basin prevents transdifferentiation and promotes the definite decision of the cell fate. These results show that the distribution of timescales of the regulatory processes is decisively important to characterize the fluctuation of cells and their differentiation process. The analyses through the epigenetic landscape and the kinetic flow on the landscape should provide a guideline to engineer cell differentiation.

Embryonic stem cells (ESCs) can proliferate indefinitely by keeping pluripotency, i.e., the ability to differentiate into any cell-lineage. ESCs, therefore, have been the focus of intense biological and medical interests. A remarkable feature of ESCs is their phenotypic heterogeneity: ESCs show large cell-to-cell fluctuation in the expression level of Nanog, which is a key factor to maintain pluripotency. Since Nanog regulates many genes in ESCs, this fluctuation should seriously affect individual cells when they start differentiation. In this paper we analyze this phenotypic fluctuation by simulating the stochastic dynamics of gene network in ESCs. The model takes account of the mutually interrelated processes of gene regulation such as binding/unbinding of transcription factors, formation/dissolution of transcription apparatus, and histone-code modification. We show the distribution of timescales of these processes is decisively important to characterize the dynamical behavior of the gene network, and that the slow formation/dissolution of transcription apparatus at the Nanog locus explains the observed large fluctuation of ESCs. The epigenetic landscapes are calculated based on the stochastic simulation, and the role of the phenotypic fluctuation in the differentiation process is analyzed through the landscape picture.

Deconstructing stem cell population heterogeneity: single-cell analysis and modeling approaches

Isogenic stem cell populations display cell-to-cell variations in a multitude of attributes including gene or protein expression, epigenetic state, morphology, proliferation and proclivity for differentiation. The origins of the observed heterogeneity and its roles in the maintenance of pluripotency and the lineage specification of stem cells remain unclear. Addressing pertinent questions will require the employment of single-cell analysis methods as traditional cell biochemical and biomolecular assays yield mostly population-average data. In addition to time-lapse microscopy and flow cytometry, recent advances in single-cell genomic, transcriptomic and proteomic profiling are reviewed. The application of multiple displacement amplification, next generation sequencing, mass cytometry and spectrometry to stem cell systems is expected to provide a wealth of information affording unprecedented levels of multiparametric characterization of cell ensembles under defined conditions promoting pluripotency or commitment. Establishing connections between single-cell analysis information and the observed phenotypes will also require suitable mathematical models. Stem cell self-renewal and differentiation are orchestrated by the coordinated regulation of subcellular, intercellular and niche-wide processes spanning multiple time scales. Here, we discuss different modeling approaches and challenges arising from their application to stem cell populations. Integrating single-cell analysis with computational methods will fill gaps in our knowledge about the functions of heterogeneity in stem cell physiology. This combination will also aid the rational design of efficient differentiation and reprogramming strategies as well as bioprocesses for the production of clinically valuable stem cell derivatives.