Modeling signaling-dependent pluripotency with Boolean logic to predict cell fate transitions

Tracking the gene expression programs and clonal relationships that underlie mast, myeloid, and T lineage specification from stem cells

T cells develop from hematopoietic progenitors in the thymus and protect against pathogens and cancer. However, the emergence of human T cell-competent blood progenitors and their subsequent specification to the T lineage have been challenging to capture in real time. Here, we leveraged a pluripotent stem cell differentiation system to understand the transcriptional dynamics and cell fate restriction events that underlie this critical developmental process. Time-resolved single-cell RNA sequencing revealed that downregulation of the multipotent hematopoietic program, upregulation of >90 lineage-associated transcription factors, and cell-cycle exit all occur within a highly coordinated developmental window. Gene-regulatory network inference uncovered a role for YBX1 in T lineage specification. We mapped the differentiation cell fate hierarchy using transcribed lineage barcoding and discovered that mast and myeloid potential bifurcate from each other early in hematopoiesis, upstream of T lineage restriction. Our systems-level analyses provide a quantitative, time-resolved model of human T cell fate specification. A record of this paper's transparent peer review process is included in the supplemental information.

OneSC: a computational platform for recapitulating cell state transitions

MOTIVATION

Computational modeling of cell state transitions has been a great interest of many in the field of developmental biology, cancer biology, and cell fate engineering because it enables performing perturbation experiments in silico more rapidly and cheaply than could be achieved in a lab. Recent advancements in single-cell RNA-sequencing (scRNA-seq) allow the capture of high-resolution snapshots of cell states as they transition along temporal trajectories. Using these high-throughput datasets, we can train computational models to generate in silico "synthetic" cells that faithfully mimic the temporal trajectories.

RESULTS

Here we present OneSC, a platform that can simulate cell state transitions using systems of stochastic differential equations govern by a regulatory network of core transcription factors (TFs). Different from many current network inference methods, OneSC prioritizes on generating Boolean network that produces faithful cell state transitions and terminal cell states that mimic real biological systems. Applying OneSC to real data, we inferred a core TF network using a mouse myeloid progenitor scRNA-seq dataset and showed that the dynamical simulations of that network generate synthetic single-cell expression profiles that faithfully recapitulate the four myeloid differentiation trajectories going into differentiated cell states (erythrocytes, megakaryocytes, granulocytes, and monocytes). Finally, through the in silico perturbations of the mouse myeloid progenitor core network, we showed that OneSC can accurately predict cell fate decision biases of TF perturbations that closely match with previous experimental observations.

AVAILABILITY AND IMPLEMENTATION

OneSC is implemented as a Python package on GitHub (https://github.com/CahanLab/oneSC) and on Zenodo (https://zenodo.org/records/14052421).

Canalizing kernel for cell fate determination

The tendency for cell fate to be robust to most perturbations, yet sensitive to certain perturbations raises intriguing questions about the existence of a key path within the underlying molecular network that critically determines distinct cell fates. Reprogramming and trans-differentiation clearly show examples of cell fate change by regulating only a few or even a single molecular switch. However, it is still unknown how to identify such a switch, called a master regulator, and how cell fate is determined by its regulation. Here, we present CAESAR, a computational framework that can systematically identify master regulators and unravel the resulting canalizing kernel, a key substructure of interconnected feedbacks that is critical for cell fate determination. We demonstrate that CAESAR can successfully predict reprogramming factors for de-differentiation into mouse embryonic stem cells and trans-differentiation of hematopoietic stem cells, while unveiling the underlying essential mechanism through the canalizing kernel. CAESAR provides a system-level understanding of how complex molecular networks determine cell fates.

OneSC: A computational platform for recapitulating cell state transitions

Computational modelling of cell state transitions has been a great interest of many in the field of developmental biology, cancer biology and cell fate engineering because it enables performing perturbation experiments in silico more rapidly and cheaply than could be achieved in a wet lab. Recent advancements in single-cell RNA sequencing (scRNA-seq) allow the capture of high-resolution snapshots of cell states as they transition along temporal trajectories. Using these high-throughput datasets, we can train computational models to generate in silico 'synthetic' cells that faithfully mimic the temporal trajectories. Here we present OneSC, a platform that can simulate synthetic cells across developmental trajectories using systems of stochastic differential equations govern by a core transcription factors (TFs) regulatory network. Different from the current network inference methods, OneSC prioritizes on generating Boolean network that produces faithful cell state transitions and steady cell states that mimic real biological systems. Applying OneSC to real data, we inferred a core TF network using a mouse myeloid progenitor scRNA-seq dataset and showed that the dynamical simulations of that network generate synthetic single-cell expression profiles that faithfully recapitulate the four myeloid differentiation trajectories going into differentiated cell states (erythrocytes, megakaryocytes, granulocytes and monocytes). Finally, through the in-silico perturbations of the mouse myeloid progenitor core network, we showed that OneSC can accurately predict cell fate decision biases of TF perturbations that closely match with previous experimental observations.

A gray box framework that optimizes a white box logical model using a black box optimizer for simulating cellular responses to perturbations

Predicting cellular responses to perturbations requires interpretable insights into molecular regulatory dynamics to perform reliable cell fate control, despite the confounding non-linearity of the underlying interactions. There is a growing interest in developing machine learning-based perturbation response prediction models to handle the non-linearity of perturbation data, but their interpretation in terms of molecular regulatory dynamics remains a challenge. Alternatively, for meaningful biological interpretation, logical network models such as Boolean networks are widely used in systems biology to represent intracellular molecular regulation. However, determining the appropriate regulatory logic of large-scale networks remains an obstacle due to the high-dimensional and discontinuous search space. To tackle these challenges, we present a scalable derivative-free optimizer trained by meta-reinforcement learning for Boolean network models. The logical network model optimized by the trained optimizer successfully predicts anti-cancer drug responses of cancer cell lines, while simultaneously providing insight into their underlying molecular regulatory mechanisms.

Human Stem Cells for Ophthalmology: Recent Advances in Diagnostic Image Analysis and Computational Modelling

Purpose of Review

To explore the advances and future research directions in image analysis and computational modelling of human stem cells (hSCs) for ophthalmological applications.

Recent Findings

hSCs hold great potential in ocular regenerative medicine due to their application in cell-based therapies and in disease modelling and drug discovery using state-of-the-art 2D and 3D organoid models. However, a deeper characterisation of their complex, multi-scale properties is required to optimise their translation to clinical practice. Image analysis combined with computational modelling is a powerful tool to explore mechanisms of hSC behaviour and aid clinical diagnosis and therapy.

Summary

Many computational models draw on a variety of techniques, often blending continuum and discrete approaches, and have been used to describe cell differentiation and self-organisation. Machine learning tools are having a significant impact in model development and improving image classification processes for clinical diagnosis and treatment and will be the focus of much future research.

Modulating the dynamics of NFκB and PI3K enhances the ensemble-level TNFR1 signaling mediated apoptotic response

Cell-to-cell variability during TNFα stimulated Tumor Necrosis Factor Receptor 1 (TNFR1) signaling can lead to single-cell level pro-survival and apoptotic responses. This variability stems from the heterogeneity in signal flow through intracellular signaling entities that regulate the balance between these two phenotypes. Using systematic Boolean dynamic modeling of a TNFR1 signaling network, we demonstrate that the signal flow path variability can be modulated to enable cells favour apoptosis. We developed a computationally efficient approach "Boolean Modeling based Prediction of Steady-state probability of Phenotype Reachability (BM-ProSPR)" to accurately predict the network's ability to settle into different phenotypes. Model analysis juxtaposed with the experimental observations revealed that NFκB and PI3K transient responses guide the XIAP behaviour to coordinate the crucial dynamic cross-talk between the pro-survival and apoptotic arms at the single-cell level. Model predicted the experimental observations that ~31% apoptosis increase can be achieved by arresting Comp1 - IKK* activity which regulates the NFκB and PI3K dynamics. Arresting Comp1 - IKK* activity causes signal flow path re-wiring towards apoptosis without significantly compromising NFκB levels, which govern adequate cell survival. Priming an ensemble of cancerous cells with inhibitors targeting the specific interaction involving Comp1 and IKK* prior to TNFα exposure could enable driving them towards apoptosis.

Modeling unveils sex differences of signaling networks in mouse embryonic stem cells

For a short period during early development of mammalian embryos, both X chromosomes in females are active, before dosage compensation is ensured through X-chromosome inactivation. In female mouse embryonic stem cells (mESCs), which carry two active X chromosomes, increased X-dosage affects cell signaling and impairs differentiation. The underlying mechanisms, however, remain poorly understood. To dissect X-dosage effects on the signaling network in mESCs, we combine systematic perturbation experiments with mathematical modeling. We quantify the response to a variety of inhibitors and growth factors for cells with one (XO) or two X chromosomes (XX). We then build models of the signaling networks in XX and XO cells through a semi-quantitative modeling approach based on modular response analysis. We identify a novel negative feedback in the PI3K/AKT pathway through GSK3. Moreover, the presence of a single active X makes mESCs more sensitive to the differentiation-promoting Activin A signal and leads to a stronger RAF1-mediated negative feedback in the FGF-triggered MAPK pathway. The differential response to these differentiation-promoting pathways can explain the impaired differentiation propensity of female mESCs.

On the History of Ecosystem Dynamical Modeling: The Rise and Promises of Qualitative Models

Ecosystem modeling is a complex and multidisciplinary modeling problem which emerged in the 1950s. It takes advantage of the computational turn in sciences to better understand anthropogenic impacts and improve ecosystem management. For that purpose, ecosystem simulation models based on difference or differential equations were built. These models were relevant for studying dynamical phenomena and still are. However, they face important limitations in data-poor situations. As a response, several formal and non-formal qualitative dynamical modeling approaches were independently developed to overcome some limitations of the existing methods. Qualitative approaches allow studying qualitative dynamics as relevant abstractions of those provided by quantitative models (e.g., response to press perturbations). Each modeling framework can be viewed as a different assemblage of properties (e.g., determinism, stochasticity or synchronous update of variable values) designed to satisfy some scientific objectives. Based on four stated objectives commonly found in complex environmental sciences ((1) grasping qualitative dynamics, (2) making as few assumptions as possible about parameter values, (3) being explanatory and (4) being predictive), our objectives were guided by the wish to model complex and multidisciplinary issues commonly found in ecosystem modeling. We then discussed the relevance of existing modeling approaches and proposed the ecological discrete-event networks (EDEN) modeling framework for this purpose. The EDEN models propose a qualitative, discrete-event, partially synchronous and possibilistic view of ecosystem dynamics. We discussed each of these properties through ecological examples and existing analysis techniques for such models and showed how relevant they are for environmental science studies.

Quantifying cancer cell plasticity with gene regulatory networks and single-cell dynamics

Phenotypic plasticity of cancer cells can lead to complex cell state dynamics during tumor progression and acquired resistance. Highly plastic stem-like states may be inherently drug-resistant. Moreover, cell state dynamics in response to therapy allow a tumor to evade treatment. In both scenarios, quantifying plasticity is essential for identifying high-plasticity states or elucidating transition paths between states. Currently, methods to quantify plasticity tend to focus on 1) quantification of quasi-potential based on the underlying gene regulatory network dynamics of the system; or 2) inference of cell potency based on trajectory inference or lineage tracing in single-cell dynamics. Here, we explore both of these approaches and associated computational tools. We then discuss implications of each approach to plasticity metrics, and relevance to cancer treatment strategies.

Leveraging developmental landscapes for model selection in Boolean gene regulatory networks

Boolean models are a well-established framework to model developmental gene regulatory networks (DGRNs) for acquisition of cellular identities. During the reconstruction of Boolean DGRNs, even if the network structure is given, there is generally a large number of combinations of Boolean functions that will reproduce the different cell fates (biological attractors). Here we leverage the developmental landscape to enable model selection on such ensembles using the relative stability of the attractors. First we show that previously proposed measures of relative stability are strongly correlated and we stress the usefulness of the one that captures best the cell state transitions via the mean first passage time (MFPT) as it also allows the construction of a cellular lineage tree. A property of great computational importance is the insensitivity of the different stability measures to changes in noise intensities. That allows us to use stochastic approaches to estimate the MFPT and thereby scale up the computations to large networks. Given this methodology, we revisit different Boolean models of Arabidopsis thaliana root development, showing that a most recent one does not respect the biologically expected hierarchy of cell states based on relative stabilities. We therefore developed an iterative greedy algorithm that searches for models which satisfy the expected hierarchy of cell states and found that its application to the root development model yields many models that meet this expectation. Our methodology thus provides new tools that can enable reconstruction of more realistic and accurate Boolean models of DGRNs.

Target Control of Asynchronous Boolean Networks

We study the target control of asynchronous Boolean networks, to identify interventions that can drive the dynamics of a given Boolean network from any initial state to the desired target attractor. Based on the application time, the control can be realised with three types of perturbations, including instantaneous, temporary and permanent perturbations. We develop efficient methods to compute the target control for a given target attractor with these three types of perturbations. We compare our methods with the stable motif-based control method on a variety of real-life biological networks to evaluate their performance. We show that our methods scale well for large Boolean networks and they are able to identify a rich set of solutions with a small number of perturbations.

Virtual cells in a virtual microenvironment recapitulate early development-like patterns in human pluripotent stem cell colonies

The mechanism by which morphogenetic signals engage the regulatory networks responsible for early embryonic tissue patterning is incompletely understood. Here, we developed a minimal gene regulatory network (GRN) model of human pluripotent stem cell (hPSC) lineage commitment and embedded it into "cellular" agents that respond to a dynamic morphogenetic signaling microenvironment. Simulations demonstrated that GRN wiring had significant non-intuitive effects on tissue pattern order, composition, and dynamics. Experimental perturbation of GRN connectivities supported model predictions and demonstrated the role of OCT4 as a master regulator of peri-gastrulation fates. Our so-called GARMEN strategy provides a multiscale computational platform to understand how single-cell-based regulatory interactions scale to tissue domains. This foundation provides new opportunities to simulate the impact of network motifs on normal and aberrant tissue development.

NETISCE: a network-based tool for cell fate reprogramming

The search for effective therapeutic targets in fields like regenerative medicine and cancer research has generated interest in cell fate reprogramming. This cellular reprogramming paradigm can drive cells to a desired target state from any initial state. However, methods for identifying reprogramming targets remain limited for biological systems that lack large sets of experimental data or a dynamical characterization. We present NETISCE, a novel computational tool for identifying cell fate reprogramming targets in static networks. In combination with machine learning algorithms, NETISCE estimates the attractor landscape and predicts reprogramming targets using signal flow analysis and feedback vertex set control, respectively. Through validations in studies of cell fate reprogramming from developmental, stem cell, and cancer biology, we show that NETISCE can predict previously identified cell fate reprogramming targets and identify potentially novel combinations of targets. NETISCE extends cell fate reprogramming studies to larger-scale biological networks without the need for full model parameterization and can be implemented by experimental and computational biologists to identify parts of a biological system relevant to the desired reprogramming task.

Langley M, Fisher CL, Aguilar-Hidalgo D, Shukla S, Yachie-Kinoshita A, Hughes M, M. McNagny K, Zandstra PW. IQCELL: A platform for predicting the effect of gene perturbations on developmental trajectories using single-cell RNA-seq data

The increasing availability of single-cell RNA-sequencing (scRNA-seq) data from various developmental systems provides the opportunity to infer gene regulatory networks (GRNs) directly from data. Herein we describe IQCELL, a platform to infer, simulate, and study executable logical GRNs directly from scRNA-seq data. Such executable GRNs allow simulation of fundamental hypotheses governing developmental programs and help accelerate the design of strategies to control stem cell fate. We first describe the architecture of IQCELL. Next, we apply IQCELL to scRNA-seq datasets from early mouse T-cell and red blood cell development, and show that the platform can infer overall over 74% of causal gene interactions previously reported from decades of research. We will also show that dynamic simulations of the generated GRN qualitatively recapitulate the effects of known gene perturbations. Finally, we implement an IQCELL gene selection pipeline that allows us to identify candidate genes, without prior knowledge. We demonstrate that GRN simulations based on the inferred set yield results similar to the original curated lists. In summary, the IQCELL platform offers a versatile tool to infer, simulate, and study executable GRNs in dynamic biological systems.

The Clinical Trials of Mesenchymal Stromal Cells Therapy

Mesenchymal stromal cells (MSCs) are a heterogeneous population of adult stem cells, which are multipotent and possess the ability to differentiate/transdifferentiate into mesodermal and nonmesodermal cell lineages. MSCs display broad immunomodulatory properties since they are capable of secreting growth factors and chemotactic cytokines. Safety, accessibility, and isolation from patients without ethical concern make MSCs valuable sources for cell therapy approaches in autoimmune, inflammatory, and degenerative diseases. Many studies have been conducted on the application of MSCs as a new therapy, but it seems that a low percentage of them is related to clinical trials, especially completed clinical trials. Considering the importance of clinical trials to develop this type of therapy as a new treatment, the current paper is aimed at describing characteristics of MSCs and reviewing relevant clinical studies registered on the NIH database during 2016-2020 to discuss recent advances on MSC-based therapeutic approaches being used in different diseases.

Context-aware synthetic biology by controller design: Engineering the mammalian cell

The rise of systems biology has ushered a new paradigm: the view of the cell as a system that processes environmental inputs to drive phenotypic outputs. Synthetic biology provides a complementary approach, allowing us to program cell behavior through the addition of synthetic genetic devices into the cellular processor. These devices, and the complex genetic circuits they compose, are engineered using a design-prototype-test cycle, allowing for predictable device performance to be achieved in a context-dependent manner. Within mammalian cells, context effects impact device performance at multiple scales, including the genetic, cellular, and extracellular levels. In order for synthetic genetic devices to achieve predictable behaviors, approaches to overcome context dependence are necessary. Here, we describe control systems approaches for achieving context-aware devices that are robust to context effects. We then consider cell fate programing as a case study to explore the potential impact of context-aware devices for regenerative medicine applications.

Unsupervised logic-based mechanism inference for network-driven biological processes

Modern analytical techniques enable researchers to collect data about cellular states, before and after perturbations. These states can be characterized using analytical techniques, but the inference of regulatory interactions that explain and predict changes in these states remains a challenge. Here we present a generalizable, unsupervised approach to generate parameter-free, logic-based models of cellular processes, described by multiple discrete states. Our algorithm employs a Hamming-distance based approach to formulate, test, and identify optimized logic rules that link two states. Our approach comprises two steps. First, a model with no prior knowledge except for the mapping between initial and attractor states is built. We then employ biological constraints to improve model fidelity. Our algorithm automatically recovers the relevant dynamics for the explored models and recapitulates key aspects of the biochemical species concentration dynamics in the original model. We present the advantages and limitations of our work and discuss how our approach could be used to infer logic-based mechanisms of signaling, gene-regulatory, or other input-output processes describable by the Boolean formalism.

Computational Stem Cell Biology: Open Questions and Guiding Principles

Computational biology is enabling an explosive growth in our understanding of stem cells and our ability to use them for disease modeling, regenerative medicine, and drug discovery. We discuss four topics that exemplify applications of computation to stem cell biology: cell typing, lineage tracing, trajectory inference, and regulatory networks. We use these examples to articulate principles that have guided computational biology broadly and call for renewed attention to these principles as computation becomes increasingly important in stem cell biology. We also discuss important challenges for this field with the hope that it will inspire more to join this exciting area.

Evaluating totipotency using criteria of increasing stringency

Totipotency is the ability of a single cell to give rise to all of the differentiated cell types that build the conceptus, yet how to capture this property in vitro remains incompletely understood. Defining totipotency relies on a variety of assays of variable stringency. Here, we describe criteria to define totipotency. We explain how distinct criteria of increasing stringency can be used to judge totipotency by evaluating candidate totipotent cell types in mice, including early blastomeres and expanded or extended pluripotent stem cells. Our data challenge the notion that expanded or extended pluripotent states harbour increased totipotent potential relative to conventional embryonic stem cells under in vitro and in vivo conditions.

Evaluating totipotency using criteria of increasing stringency

Totipotency is the ability of a single cell to give rise to all of the differentiated cell types that build the conceptus, yet how to capture this property in vitro remains incompletely understood. Defining totipotency relies on a variety of assays of variable stringency. Here, we describe criteria to define totipotency. We explain how distinct criteria of increasing stringency can be used to judge totipotency by evaluating candidate totipotent cell types in mice, including early blastomeres and expanded or extended pluripotent stem cells. Our data challenge the notion that expanded or extended pluripotent states harbour increased totipotent potential relative to conventional embryonic stem cells under in vitro and in vivo conditions.

Identification of potential transcription factors that enhance human iPSC generation

Although many factors have been identified and used to enhance the iPSC reprogramming process, its efficiency remains quite low. In addition, reprogramming efficacy has been evidenced to be affected by disease mutations that are present in patient samples. In this study, using RNA-seq platform we have identified and validated the differential gene expression of five transcription factors (TFs) (GBX2, NANOGP8, SP8, PEG3, and ZIC1) that were associated with a remarkable increase in the number of iPSC colonies generated from a patient with Parkinson's disease. We have applied different bioinformatics tools (Gene ontology, protein-protein interaction, and signaling pathways analyses) to investigate the possible roles of these TFs in pluripotency and developmental process. Interestingly, GBX2, NANOGP8, SP8, PEG3, and ZIC1 were found to play a role in maintaining pluripotency, regulating self-renewal stages, and interacting with other factors that are involved in pluripotency regulation including OCT4, SOX2, NANOG, and KLF4. Therefore, the TFs identified in this study could be used as additional transcription factors that enhance reprogramming efficiency to boost iPSC generation technology.

Development of a kinetic model expressing anomalous phenomena in human induced pluripotent stem cell culture

During culture with feeder cells, deviation from the undifferentiated state of human induced pluripotent stem cells (hiPSCs) occurs at a very low frequency. Anomalous cell migration in central and peripheral regions of hiPSC colonies has been suggested to be the trigger for this phenomenon. To confirm this hypothesis, sequential cell migration prior to deviation must be demonstrated. This has been difficult using in vitro methods. We therefore developed a kinetic model with a proposed definition of anomalous cell migration as continuous relatively fast or slow cell migration. The developed model was validated via in silico reproduction of deviation phenomenon observed in vitro, such as the positions of deviated cells in a colony and the frequency of deviation in culture. This model suggests that anomalous cell migration-driven hiPSC deviation can be explained by two factors: a mechanical stimulus, represented by cell migration, and duration of the mechanical stimulus. The factor "duration of mechanical stimulus" sets our model apart from others, and helps to realize the ultra-rare trigger (approximately 10^-5) of deviation from the undifferentiated state in hiPSC culture.

The Importance of Computational Modeling in Stem Cell Research

The generation of large amounts of omics data is increasingly enabling not only the processing and analysis of large data sets but also the development of computational models in the field of stem cell research. Although computational models have been proposed in recent decades, we believe that the stem cell community is not fully aware of the potentiality of computational modeling in guiding their experimental research. In this regard, we discuss how single-cell technologies provide the right framework for computational modeling at different scales of biological organization in order to address challenges in the stem cell field and to guide experimentalists in the design of new strategies for stem cell therapies and treatment of congenital disorders.

An integrative method to predict signalling perturbations for cellular transitions

Induction of specific cellular transitions is of clinical importance, as it allows to revert disease cellular phenotype, or induce cellular reprogramming and differentiation for regenerative medicine. Signalling is a convenient way to accomplish such transitions without transfer of genetic material. Here we present the first general computational method that systematically predicts signalling molecules, whose perturbations induce desired cellular transitions. This probabilistic method integrates gene regulatory networks (GRNs) with manually-curated signalling pathways obtained from MetaCore from Clarivate Analytics, to model how signalling cues are received and processed in the GRN. The method was applied to 219 cellular transition examples, including cell type transitions, and overall correctly predicted experimentally validated signalling molecules, consistently outperforming other well-established approaches, such as differential gene expression and pathway enrichment analyses. Further, we validated our method predictions in the case of rat cirrhotic liver, and identified the activation of angiopoietins receptor Tie2 as a potential target for reverting the disease phenotype. Experimental results indicated that this perturbation induced desired changes in the gene expression of key TFs involved in fibrosis and angiogenesis. Importantly, this method only requires gene expression data of the initial and desired cell states, and therefore is suited for the discovery of signalling interventions for disease treatments and cellular therapies.

Inferring the transcriptional regulatory mechanism of signal-dependent gene expression via an integrative computational approach

Eukaryotic transcription factors (TFs) coordinate different upstream signals to regulate the expression of their target genes. To unveil this regulatory network in B-cell receptor signaling, we developed a computational pipeline to systematically analyze the extracellular signal-regulated kinase (ERK)- and IκB kinase (IKK)-dependent transcriptome responses. We combined a bilinear regression method and kinetic modeling to identify the signal-to-TF and TF-to-gene dynamics, respectively. We input a set of time-course experimental data for B cells and concentrated on transcriptional activators. The results show that the combination of TFs differentially controlled by ERK and IKK could contribute divergent expression dynamics in orchestrating the B-cell response. Our findings provide insights into the regulatory mechanisms underlying signal-dependent gene expression in eukaryotic cells.

Decoding pluripotency: Genetic screens to interrogate the acquisition, maintenance, and exit of pluripotency

Pluripotent stem cells have the ability to unlimitedly self-renew and differentiate to any somatic cell lineage. A number of systems biology approaches have been used to define this pluripotent state. Complementary to systems level characterization, genetic screens offer a unique avenue to functionally interrogate the pluripotent state and identify the key players in pluripotency acquisition and maintenance, exit of pluripotency, and lineage differentiation. Here we review how genetic screens have helped us decode pluripotency regulation. We will summarize results from RNA interference (RNAi) based screens, discuss recent advances in CRISPR/Cas-based genetic perturbation methods, and how these advances have made it possible to more comprehensively interrogate pluripotency and differentiation through genetic screens. Such investigations will not only provide a better understanding of this unique developmental state, but may enhance our ability to use pluripotent stem cells as an experimental model to study human development and disease progression. Functional interrogation of pluripotency also provides a valuable roadmap for utilizing genetic perturbation to gain systems level understanding of additional cellular states, from later stages of development to pathological disease states. This article is categorized under: Developmental Biology > Stem Cell Biology and Regeneration Developmental Biology > Developmental Processes in Health and Disease Biological Mechanisms > Cell Fates.

Automated Segmentation of Fluorescence Microscopy Images for 3D Cell Detection in human-derived Cardiospheres

The ‘cardiosphere’ is a 3D cluster of cardiac progenitor cells recapitulating a stem cell niche-like microenvironment with a potential for disease and regeneration modelling of the failing human myocardium. In this multicellular 3D context, it is extremely important to decrypt the spatial distribution of cell markers for dissecting the evolution of cellular phenotypes by direct quantification of fluorescent signals in confocal microscopy. In this study, we present a fully automated method, named CARE (‘CARdiosphere Evaluation’), for the segmentation of membranes and cell nuclei in human-derived cardiospheres. The proposed method is tested on twenty 3D-stacks of cardiospheres, for a total of 1160 images. Automatic results are compared with manual annotations and two open-source software designed for fluorescence microscopy. CARE performance was excellent in cardiospheres membrane segmentation and, in cell nuclei detection, the algorithm achieved the same performance as two expert operators. To the best of our knowledge, CARE is the first fully automated algorithm for segmentation inside in vitro 3D cell spheroids, including cardiospheres. The proposed approach will provide, in the future, automated quantitative analysis of markers distribution within the cardiac niche-like environment, enabling predictive associations between cell mechanical stresses and dynamic phenotypic changes.

Boolean network modeling of β-cell apoptosis and insulin resistance in type 2 diabetes mellitus

BACKGROUND

Major alteration in lifestyle of human population has promoted Type 2 diabetes mellitus (T2DM) to the level of an epidemic. This metabolic disorder is characterized by insulin resistance and pancreatic β-cell dysfunction and apoptosis, triggered by endoplasmic reticulum (ER) stress, oxidative stress and cytokines. Computational modeling is necessary to consolidate information from various sources in order to obtain a comprehensive understanding of the pathogenesis of T2DM and to investigate possible interventions by performing in silico simulations.

RESULTS

In this paper, we propose a Boolean network model integrating the insulin resistance pathway with pancreatic β-cell apoptosis pathway which are responsible for T2DM. The model has five input signals, i.e. ER stress, oxidative stress, tumor necrosis factor α (TNF α), Fas ligand (FasL), and interleukin-6 (IL-6). We performed dynamical simulations using random order asynchronous update and with different combinations of the input signals. From the results, we observed that the proposed model made predictions that closely resemble the expression levels of genes in T2DM as reported in the literature.

CONCLUSION

The proposed model can make predictions about expression levels of genes in T2DM that are in concordance with literature. Although experimental validation of the model is beyond the scope of this study, the model can be useful for understanding the aetiology of T2DM and discovery of therapeutic intervention for this prevalent complex disease. The files of our model and results are available at https://github.com/JieZheng-ShanghaiTech/boolean-t2dm .

Ensembles, dynamics, and cell types: Revisiting the statistical mechanics perspective on cellular regulation

Genetic regulatory networks control ontogeny. For fifty years Boolean networks have served as models of such systems, ranging from ensembles of random Boolean networks as models for generic properties of gene regulation to working dynamical models of a growing number of sub-networks of real cells. At the same time, their statistical mechanics has been thoroughly studied. Here we recapitulate their original motivation in the context of current theoretical and empirical research. We discuss ensembles of random Boolean networks whose dynamical attractors model cell types. A sub-ensemble is the critical ensemble. There is now strong evidence that genetic regulatory networks are dynamically critical, and that evolution is exploring the critical sub-ensemble. The generic properties of this sub-ensemble predict essential features of cell differentiation. In particular, the number of attractors in such networks scales as the DNA content raised to the 0.63 power. Data on the number of cell types as a function of the DNA content per cell shows a scaling relationship of 0.88. Thus, the theory correctly predicts a power law relationship between the number of cell types and the DNA contents per cell, and a comparable slope. We discuss these new scaling values and show prospects for new research lines for Boolean networks as a base model for systems biology.

A common molecular logic determines embryonic stem cell self-renewal and reprogramming

During differentiation and reprogramming, new cell identities are generated by reconfiguration of gene regulatory networks. Here, we combined automated formal reasoning with experimentation to expose the logic of network activation during induction of naïve pluripotency. We find that a Boolean network architecture defined for maintenance of naïve state embryonic stem cells (ESC) also explains transcription factor behaviour and potency during resetting from primed pluripotency. Computationally identified gene activation trajectories were experimentally substantiated at single‐cell resolution by RT–qPCR. Contingency of factor availability explains the counterintuitive observation that Klf2, which is dispensable for ESC maintenance, is required during resetting. We tested 124 predictions formulated by the dynamic network, finding a predictive accuracy of 77.4%. Finally, we show that this network explains and predicts experimental observations of somatic cell reprogramming. We conclude that a common deterministic program of gene regulation is sufficient to govern maintenance and induction of naïve pluripotency. The tools exemplified here could be broadly applied to delineate dynamic networks underlying cell fate transitions.

Regulatory network characterization in development: challenges and opportunities

Embryonic development and stem cell differentiation, during which coordinated cell fate specification takes place in a spatial and temporal context, serve as a paradigm for studying the orderly assembly of gene regulatory networks (GRNs) and the fundamental mechanism of GRNs in driving lineage determination. However, knowledge of reliable GRN annotation for dynamic development regulation, particularly for unveiling the complex temporal and spatial architecture of tissue stem cells, remains inadequate. With the advent of single-cell RNA sequencing technology, elucidating GRNs in development and stem cell processes poses both new challenges and unprecedented opportunities. This review takes a snapshot of some of this work and its implication in the regulative nature of early mammalian development and specification of the distinct cell types during embryogenesis.

Inference of Genome-Scale Gene Regulatory Networks: Are There Differences in Biological and Clinical Validations?

Causal networks, e.g., gene regulatory networks (GRNs) inferred from gene expression data, contain a wealth of information but are defying simple, straightforward and low-budget experimental validations. In this paper, we elaborate on this problem and discuss distinctions between biological and clinical validations. As a result, validation differences for GRNs reflect known differences between basic biological and clinical research questions making the validations context specific. Hence, the meaning of biologically and clinically meaningful GRNs can be very different. For a concerted approach to a problem of this size, we suggest the establishment of the HUMAN GENE REGULATORY NETWORK PROJECT which provides the information required for biological and clinical validations alike.

A stochastic and dynamical view of pluripotency in mouse embryonic stem cells

Pluripotent embryonic stem cells are of paramount importance for biomedical sciences because of their innate ability for self-renewal and differentiation into all major cell lines. The fateful decision to exit or remain in the pluripotent state is regulated by complex genetic regulatory networks. The rapid growth of single-cell sequencing data has greatly stimulated applications of statistical and machine learning methods for inferring topologies of pluripotency regulating genetic networks. The inferred network topologies, however, often only encode Boolean information while remaining silent about the roles of dynamics and molecular stochasticity inherent in gene expression. Herein we develop a framework for systematically extending Boolean-level network topologies into higher resolution models of networks which explicitly account for the promoter architectures and gene state switching dynamics. We show the framework to be useful for disentangling the various contributions that gene switching, external signaling, and network topology make to the global heterogeneity and dynamics of transcription factor populations. We find the pluripotent state of the network to be a steady state which is robust to global variations of gene switching rates which we argue are a good proxy for epigenetic states of individual promoters. The temporal dynamics of exiting the pluripotent state, on the other hand, is significantly influenced by the rates of genetic switching which makes cells more responsive to changes in extracellular signals.

Modeling signaling-dependent pluripotency with Boolean logic to predict cell fate transitions

Pluripotent stem cells (PSCs) exist in multiple stable states, each with specific cellular properties and molecular signatures. The mechanisms that maintain pluripotency, or that cause its destabilization to initiate development, are complex and incompletely understood. We have developed a model to predict stabilized PSC gene regulatory network (GRN) states in response to input signals. Our strategy used random asynchronous Boolean simulations (R-ABS) to simulate single-cell fate transitions and strongly connected components (SCCs) strategy to represent population heterogeneity. This framework was applied to a reverse-engineered and curated core GRN for mouse embryonic stem cells (mESCs) and used to simulate cellular responses to combinations of five signaling pathways. Our simulations predicted experimentally verified cell population compositions and input signal combinations controlling specific cell fate transitions. Extending the model to PSC differentiation, we predicted a combination of signaling activators and inhibitors that efficiently and robustly generated a Cdx2⁺Oct4^- cells from naïve mESCs. Overall, this platform provides new strategies to simulate cell fate transitions and the heterogeneity that typically occurs during development and differentiation.