A Mathematical Framework for Understanding Four-Dimensional Heterogeneous Differentiation of CD4+ T Cells

Multistability and predominant hybrid phenotypes in a four node mutually repressive network of Th1/Th2/Th17/Treg differentiation

Elucidating the emergent dynamics of cellular differentiation networks is crucial to understanding cell-fate decisions. Toggle switch - a network of mutually repressive lineage-specific transcription factors A and B - enables two phenotypes from a common progenitor: (high A, low B) and (low A, high B). However, the dynamics of networks enabling differentiation of more than two phenotypes from a progenitor cell has not been well-studied. Here, we investigate the dynamics of a four-node network A, B, C, and D inhibiting each other, forming a toggle tetrahedron. Our simulations show that this network is multistable and predominantly allows for the co-existence of six hybrid phenotypes where two of the nodes are expressed relatively high as compared to the remaining two, for instance (high A, high B, low C, low D). Finally, we apply our results to understand naïve CD4+ T cell differentiation into Th1, Th2, Th17 and Treg subsets, suggesting Th1/Th2/Th17/Treg decision-making to be a two-step process.

A mathematical model for understanding synergistic regulations and paradoxical feedbacks in the shoot apical meristem

The shoot apical meristem (SAM) is the primary stem cell niche in plant shoots. Stem cells in the SAM are controlled by an intricate regulatory network, including negative feedback between WUSCHEL (WUS) and CLAVATA3 (CLV3). Recently, we identified a group of signals, Epidermal Patterning Factor-Like (EPFL) proteins, that are produced at the peripheral region and are important for SAM homeostasis. Here, we present a mathematical model for the SAM regulatory network. The model revealed that the SAM uses EPFL and signals such as HAIRY MERISTEM from the middle in a synergistic manner to constrain both WUS and CLV3. We found that interconnected negative and positive feedbacks between WUS and CLV3 ensure stable WUS expression in the SAM when facing perturbations, and the positive feedback loop also maintains distinct cell populations containing WUS on and CLV3 on cells in the apical-basal direction. Furthermore, systematic perturbations of the parameters revealed a tradeoff between optimizations of multiple patterning features. Our results provide a holistic view of the regulation of SAM patterning in multiple dimensions. They give insights into how Arabidopsis integrates signals from lateral and apical-basal axes to control the SAM patterning, and they shed light into design principles that may be widely useful for understanding regulatory networks of stem cell niche.

Early Detection of Daylengths with a Feedforward Circuit Coregulated by Circadian and Diurnal Cycles

Light-entrained circadian clocks confer rhythmic dynamics of cellular and molecular activities to animals and plants. These intrinsic clocks allow stable anticipations to light-dark (diel) cycles. Many genes in the model plant Arabidopsis thaliana are regulated by diel cycles via pathways independent of the clock, suggesting that the integration of circadian and light signals is important for the fitness of plants. Previous studies of light-clock signal integrations have focused on moderate phase adjustment of the two signals. However, dynamical features of integrations across a broad range of phases remain elusive. Phosphorylation of ribosomal protein of the small subunit 6 (eS6), a ubiquitous post-translational modification across kingdoms, is influenced by the circadian clock and the light-dark (diel) cycle in an opposite manner. To understand this striking phenomenon and its underlying information processing capabilities, we built a mathematical model for the eS6 phosphorylation (eS6-P) control circuit. We found that the dynamics of eS6-P can be explained by a feedforward circuit with inputs from both circadian and diel cycles. Furthermore, the early day response of this circuit with dual rhythmic inputs is sensitive to the changes in daylength, including both transient and gradual changes observed in realistic light intervals across a year, because of weather and seasons. By analyzing published gene expression data, we found that the dynamics produced by the eS6-P control circuit can be observed in the expression profiles of a large number of genes. Our work provides mechanistic insights into the complex dynamics of a ribosomal protein, and it proposes a previously underappreciated function of the circadian clock, which not only prepares organisms for normal diel cycles but also helps to detect both transient and seasonal changes with a predictive power.

Searching for the Mechanisms of Mammalian Cellular Aging Through Underlying Gene Regulatory Networks

Aging attracts the attention throughout the history of humankind. However, it is still challenging to understand how the internal driving forces, for example, the fundamental building blocks of life, such as genes and proteins, as well as the environments work together to determine longevity in mammals. In this study, we built a gene regulatory network for mammalian cellular aging based on the experimental literature and quantify its underlying driving force for the dynamics as potential and flux landscape. We found three steady-state attractors: a fast-aging state attractor, slow-aging state attractor, and intermediate state attractor. The system can switch from one state attractor to another driven by the intrinsic or external forces through the genetics and the environment. We identified the dominant path from the slow-aging state directly to the fast-aging state. We also identified the dominant path from slow-aging to fast-aging through an intermediate state. We quantified the evolving landscape for revealing the dynamic characteristics of aging through certain regulation changes in time. We also predicted the key genes and regulations for fast-aging and slow-aging through the analysis of the stability for landscape basins. We also found the oscillation dynamics between fast-aging and slow-aging and showed that more energy is required to sustain such oscillations. We found that the flux is the dynamic cause and the entropy production rate the thermodynamic origin of the phase transitions or the bifurcations between the three-state phase and oscillation phase. The landscape quantification provides a global and physical approach to explore the underlying mechanisms of cellular aging in mammals.

Multiscale modeling of vertebrate limb development

We review the current state of mathematical modeling of cartilage pattern formation in vertebrate limbs. We place emphasis on several reaction-diffusion type models that have been proposed in the last few years. These models are grounded in more detailed knowledge of the relevant regulatory processes than previous ones but generally refer to different molecular aspects of these processes. Considering these models in light of comparative phylogenomics permits framing of hypotheses on the evolutionary order of appearance of the respective mechanisms and their roles in the fin-to-limb transition. This article is categorized under: Analytical and Computational Methods > Computational Methods Models of Systems Properties and Processes > Mechanistic Models Developmental Biology > Developmental Processes in Health and Disease Analytical and Computational Methods > Analytical Methods.

Integrative Transcriptomic Analysis Reveals a Multiphasic Epithelial-Mesenchymal Spectrum in Cancer and Non-tumorigenic Cells

Epithelial-mesenchymal transition (EMT), the conversion between rigid epithelial cells and motile mesenchymal cells, is a reversible cellular process involved in tumorigenesis, metastasis, and chemoresistance. Numerous studies have found that several types of tumor cells show a high degree of cell-to-cell heterogeneity in terms of their gene expression signatures and cellular phenotypes related to EMT. Recently, the prevalence and importance of partial or intermediate EMT states have been reported. It is unclear, however, whether there is a general pattern of cancer cell distribution in terms of the overall expression of epithelial-related genes and mesenchymal-related genes, and how this distribution is related to EMT process in normal cells. In this study, we performed integrative transcriptomic analysis that combines cancer cell transcriptomes, time course data of EMT in non-tumorigenic epithelial cells, and epithelial cells with perturbations of key EMT factors. Our statistical analysis shows that cancer cells are widely distributed in the EMT spectrum, and the majority of these cells can be described by an EMT path that connects the epithelial and the mesenchymal states via a hybrid expression region in which both epithelial genes and mesenchymal genes are highly expressed overall. We found that key patterns of this EMT path are observed in EMT progression in non-tumorigenic cells and that transcription factor ZEB1 plays a key role in defining this EMT path via diverse gene regulatory circuits connecting to epithelial genes. We performed Gene Set Variation Analysis to show that the cancer cells at hybrid EMT states also possess hybrid cellular phenotypes with both high migratory and high proliferative potentials. Our results reveal critical patterns of cancer cells in the EMT spectrum and their relationship to the EMT process in normal cells, and provide insights into the mechanistic basis of cancer cell heterogeneity and plasticity.

Combinatorial perturbation analysis reveals divergent regulations of mesenchymal genes during epithelial-to-mesenchymal transition

Epithelial-to-mesenchymal transition (EMT), a fundamental transdifferentiation process in development, produces diverse phenotypes in different physiological or pathological conditions. Many genes involved in EMT have been identified to date, but mechanisms contributing to the phenotypic diversity and those governing the coupling between the dynamics of epithelial (E) genes and that of the mesenchymal (M) genes are unclear. In this study, we employed combinatorial perturbations to mammary epithelial cells to induce a series of EMT phenotypes by manipulating two essential EMT-inducing elements, namely TGF-β and ZEB1. By measuring transcriptional changes in more than 700 E-genes and M-genes, we discovered that the M-genes exhibit a significant diversity in their dependency to these regulatory elements and identified three groups of M-genes that are controlled by different regulatory circuits. Notably, functional differences were detected among the M-gene clusters in motility regulation and in survival of breast cancer patients. We computationally predicted and experimentally confirmed that the reciprocity and reversibility of EMT are jointly regulated by ZEB1. Our integrative analysis reveals the key roles of ZEB1 in coordinating the dynamics of a large number of genes during EMT, and it provides new insights into the mechanisms for the diversity of EMT phenotypes.

An enriched network motif family regulates multistep cell fate transitions with restricted reversibility

Multistep cell fate transitions with stepwise changes of transcriptional profiles are common to many developmental, regenerative and pathological processes. The multiple intermediate cell lineage states can serve as differentiation checkpoints or branching points for channeling cells to more than one lineages. However, mechanisms underlying these transitions remain elusive. Here, we explored gene regulatory circuits that can generate multiple intermediate cellular states with stepwise modulations of transcription factors. With unbiased searching in the network topology space, we found a motif family containing a large set of networks can give rise to four attractors with the stepwise regulations of transcription factors, which limit the reversibility of three consecutive steps of the lineage transition. We found that there is an enrichment of these motifs in a transcriptional network controlling the early T cell development, and a mathematical model based on this network recapitulates multistep transitions in the early T cell lineage commitment. By calculating the energy landscape and minimum action paths for the T cell model, we quantified the stochastic dynamics of the critical factors in response to the differentiation signal with fluctuations. These results are in good agreement with experimental observations and they suggest the stable characteristics of the intermediate states in the T cell differentiation. These dynamical features may help to direct the cells to correct lineages during development. Our findings provide general design principles for multistep cell linage transitions and new insights into the early T cell development. The network motifs containing a large family of topologies can be useful for analyzing diverse biological systems with multistep transitions.

The functions of cells are dynamically controlled in many biological processes including development, regeneration and disease progression. Cell fate transition, or the switch of cellular functions, often involves multiple steps. The intermediate stages of the transition provide the biological systems with the opportunities to regulate the transitions in a precise manner. These transitions are controlled by key regulatory genes of which the expression shows stepwise patterns, but how the interactions of these genes can determine the multistep processes was unclear. Here, we present a comprehensive analysis on the design principles of gene circuits that govern multistep cell fate transition. We found a large network family with common structural features that can generate systems with the ability to control three consecutive steps of the transition. We found that this type of networks is enriched in a gene circuit controlling the development of T lymphocyte, a crucial type of immune cells. We performed mathematical modeling using this gene circuit and we recapitulated the stepwise and irreversible loss of stem cell properties of the developing T lymphocytes. Our findings can be useful to analyze a wide range of gene regulatory networks controlling multistep cell fate transitions.

Modeling the dynamic behavior of biochemical regulatory networks

Strategies for modeling the complex dynamical behavior of gene/protein regulatory networks have evolved over the last 50 years as both the knowledge of these molecular control systems and the power of computing resources have increased. Here, we review a number of common modeling approaches, including Boolean (logical) models, systems of piecewise-linear or fully non-linear ordinary differential equations, and stochastic models (including hybrid deterministic/stochastic approaches). We discuss the pro's and con's of each approach, to help novice modelers choose a modeling strategy suitable to their problem, based on the type and bounty of available experimental information. We illustrate different modeling strategies in terms of some abstract network motifs, and in the specific context of cell cycle regulation.

Exploring intermediate cell states through the lens of single cells

As our catalog of cell states expands, appropriate characterization of these states and the transitions between them is crucial. Here we discuss the roles of intermediate cell states (ICSs) in this growing collection. We begin with definitions and discuss evidence for the existence of ICSs and their relevance in various tissues. We then provide a list of possible functions for ICSs with examples. Finally, we describe means by which ICSs and their functional roles can be identified from single-cell data or predicted from models.

Controlling Stochasticity in Epithelial-Mesenchymal Transition Through Multiple Intermediate Cellular States

Epithelial-mesenchymal transition (EMT) is an instance of cellular plasticity that plays critical roles in development, regeneration and cancer progression. Recent studies indicate that the transition between epithelial and mesenchymal states is a multi-step and reversible process in which several intermediate phenotypes might coexist. These intermediate states correspond to various forms of stem-like cells in the EMT system, but the function of the multi-step transition or the multiple stem cell phenotypes is unclear. Here, we use mathematical models to show that multiple intermediate phenotypes in the EMT system help to attenuate the overall fluctuations of the cell population in terms of phenotypic compositions, thereby stabilizing a heterogeneous cell population in the EMT spectrum. We found that the ability of the system to attenuate noise on the intermediate states depends on the number of intermediate states, indicating the stem-cell population is more stable when it has more sub-states. Our study reveals a novel advantage of multiple intermediate EMT phenotypes in terms of systems design, and it sheds light on the general design principle of heterogeneous stem cell population.

Realization of tristability in a multiplicatively coupled dual-loop genetic network

Multistability is a crucial recurring theme in cell signaling. Multistability is attributed to the presence of positive feedback loops, but the general condition and essential mechanism for realizing multistability remain unclear. Here, we build a generic circuit model comprising two transcription factors and a microRNA, representing a kind of core architecture in gene regulatory networks. The circuit can be decomposed into two positive feedback loops (PFLs) or one PFL and one negative feedback loop (NFL), which are multiplicatively coupled. Bifurcation analyses of the model reveal that the circuit can achieve tristability through four kinds of bifurcation scenarios when parameter values are varied in a wide range. We formulate the general requirement for tristability in terms of logarithmic gain of the circuit. The parameter ranges for tristability and possible transition routes among steady states are determined by the combination of gain features of individual feedback loops. Coupling two PFLs with bistability or one NFL with a bistable PFL is most likely to generate tristability, but the underlying mechanisms are largely different. We also interpret published results and make testable predictions. This work sheds new light on interlinking feedback loops to realize tristability. The proposed theoretical framework can be of wide applicability.