Cell Fate Decision as High-Dimensional Critical State Transition

Dual dynamics of mitochondrial permeability transition pore opening

Mitochondria play an essential role in bioenergetics and cellular Ca\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}^{2+}$$\end{document}2+ handling. The mitochondrial permeability transition pore (mPTP) is a non-specific channel located in the inner mitochondrial membrane. Long-lasting openings of the pore allow the rapid passage of ions and large molecules, which can result in cell death. The mPTP also exhibits transient, low conductance openings that contribute to Ca\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}^{2+}$$\end{document}2+ homeostasis. Although many regulators of the pore have been identified, none of them uniquely governs the passage between the two operating modes, which thus probably relies on a still unidentified network of interactions. By developing a core computational model for mPTP opening under the control of mitochondrial voltage and Ca\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}^{2+}$$\end{document}2+, we uncovered the existence of a positive feedback loop leading to bistability. The characteristics of the two stable steady-states correspond to those of the two opening states. When inserted in a full model of Ca\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}^{2+}$$\end{document}2+ handling by mitochondria, our description of the pore reproduces observations in mitochondrial suspensions. Moreover, the model predicted the occurrence of hysteresis in the switching between the two modes, upon addition and removal of free Ca\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}^{2+}$$\end{document}2+ in the extra-mitochondrial medium. Stochastic simulations then confirmed that the pore can undergo transient openings resembling those observed in intact cells.

Analysis of Single-Cell Gene Pair Coexpression Landscapes by Stochastic Kinetic Modeling Reveals Gene-Pair Interactions in Development

Single-cell transcriptomics is advancing discovery of the molecular determinants of cell identity, while spurring development of novel data analysis methods. Stochastic mathematical models of gene regulatory networks help unravel the dynamic, molecular mechanisms underlying cell-to-cell heterogeneity, and can thus aid interpretation of heterogeneous cell-states revealed by single-cell measurements. However, integrating stochastic gene network models with single cell data is challenging. Here, we present a method for analyzing single-cell gene-pair coexpression patterns, based on biophysical models of stochastic gene expression and interaction dynamics. We first developed a high-computational-throughput approach to stochastic modeling of gene-pair coexpression landscapes, based on numerical solution of gene network Master Equations. We then comprehensively catalogued coexpression patterns arising from tens of thousands of gene-gene interaction models with different biochemical kinetic parameters and regulatory interactions. From the computed landscapes, we obtain a low-dimensional "shape-space" describing distinct types of coexpression patterns. We applied the theoretical results to analysis of published single cell RNA sequencing data and uncovered complex dynamics of coexpression among gene pairs during embryonic development. Our approach provides a generalizable framework for inferring evolution of gene-gene interactions during critical cell-state transitions.

Early-warning signals using dynamical network markers selected by covariance

It is an important issue, particularly in the context of sustainable society, to predict critical transitions across which a system state abruptly shifts toward a contrasting state. In this study, we propose an indicator of critical transitions in multivariate dynamical systems, based on the concept of the dynamical network marker (DNM). The DNM is originally defined based on the eigendecomposition of the Jacobian matrix of a nonlinear system and corresponds to large-magnitude components of the dominant eigenvector, which contributes primarily to transitions. Our DNM-based indicator is derived from the sample covariance matrix of state variables in a target system. Simulation results to predict transitions in complex network systems consisting of a harvesting model consistently show the superiority of our indicator as a precursor of transitions regardless of network structure characteristics, as compared to a conventional indicator.

Resolution of Complex Issues in Genome Regulation and Cancer Requires Non-Linear and Network-Based Thermodynamics

The apparent lack of success in curing cancer that was evidenced in the last four decades of molecular medicine indicates the need for a global re-thinking both its nature and the biological approaches that we are taking in its solution. The reductionist, one gene/one protein method that has served us well until now, and that still dominates in biomedicine, requires complementation with a more systemic/holistic approach, to address the huge problem of cross-talk between more than 20,000 protein-coding genes, about 100,000 protein types, and the multiple layers of biological organization. In this perspective, the relationship between the chromatin network organization and gene expression regulation plays a fundamental role. The elucidation of such a relationship requires a non-linear thermodynamics approach to these biological systems. This change of perspective is a necessary step for developing successful 'tumour-reversion' therapeutic strategies.

Revealing Dynamic Mechanisms of Cell Fate Decisions From Single-Cell Transcriptomic Data

Cell fate decisions play a pivotal role in development, but technologies for dissecting them are limited. We developed a multifunction new method, Topographer, to construct a "quantitative" Waddington's landscape of single-cell transcriptomic data. This method is able to identify complex cell-state transition trajectories and to estimate complex cell-type dynamics characterized by fate and transition probabilities. It also infers both marker gene networks and their dynamic changes as well as dynamic characteristics of transcriptional bursting along the cell-state transition trajectories. Applying this method to single-cell RNA-seq data on the differentiation of primary human myoblasts, we not only identified three known cell types, but also estimated both their fate probabilities and transition probabilities among them. We found that the percent of genes expressed in a bursty manner is significantly higher at (or near) the branch point (~97%) than before or after branch (below 80%), and that both gene-gene and cell-cell correlation degrees are apparently lower near the branch point than away from the branching. Topographer allows revealing of cell fate mechanisms in a coherent way at three scales: cell lineage (macroscopic), gene network (mesoscopic), and gene expression (microscopic).

Phenotypic transitions enacted by simulated microgravity do not alter coherence in gene transcription profile

Cells in simulated microgravity undergo a reversible morphology switch, causing the appearance of two distinct phenotypes. Despite the dramatic splitting into an adherent-fusiform and a floating-spherical population, when looking at the gene-expression phase space, cell transition ends up in a largely invariant gene transcription profile characterized by only mild modifications in the respective Pearson's correlation coefficients. Functional changes among the different phenotypes emerging in simulated microgravity using random positioning machine are adaptive modifications-as cells promptly recover their native phenotype when placed again into normal gravity-and do not alter the internal gene coherence. However, biophysical constraints are required to drive phenotypic commitment in an appropriate way, compatible with physiological requirements, given that absence of gravity foster cells to oscillate between different attractor states, thus preventing them to acquire a exclusive phenotype. This is a proof-of-concept of the adaptive properties of gene-expression networks supporting very different phenotypes by coordinated 'profile preserving' modifications.

Drugs modulating stochastic gene expression affect the erythroid differentiation process

To better understand the mechanisms behind cells decision-making to differentiate, we assessed the influence of stochastic gene expression (SGE) modulation on the erythroid differentiation process. It has been suggested that stochastic gene expression has a role in cell fate decision-making which is revealed by single-cell analyses but studies dedicated to demonstrate the consistency of this link are still lacking. Recent observations showed that SGE significantly increased during differentiation and a few showed that an increase of the level of SGE is accompanied by an increase in the differentiation process. However, a consistent relation in both increasing and decreasing directions has never been shown in the same cellular system. Such demonstration would require to be able to experimentally manipulate simultaneously the level of SGE and cell differentiation in order to observe if cell behavior matches with the current theory. We identified three drugs that modulate SGE in primary erythroid progenitor cells. Both Artemisinin and Indomethacin decreased SGE and reduced the amount of differentiated cells. On the contrary, a third component called MB-3 simultaneously increased the level of SGE and the amount of differentiated cells. We then used a dynamical modelling approach which confirmed that differentiation rates were indeed affected by the drug treatment. Using single-cell analysis and modeling tools, we provide experimental evidence that, in a physiologically relevant cellular system, SGE is linked to differentiation.

Enhanced Waddington landscape model with cell-cell communication can explain molecular mechanisms of self-organization

MOTIVATION

The molecular mechanisms of self-organization that orchestrate embryonic cells to create astonishing patterns have been among major questions of developmental biology. It is recently shown that embryonic stem cells (ESCs), when cultured in particular micropatterns, can self-organize and mimic the early steps of pre-implantation embryogenesis. A systems-biology model to address this observation from a dynamical systems perspective is essential and can enhance understanding of the phenomenon.

RESULTS

Here, we propose a multicellular mathematical model for pattern formation during in vitro gastrulation of human ESCs. This model enhances the basic principles of Waddington epigenetic landscape with cell-cell communication, in order to enable pattern and tissue formation. We have shown the sufficiency of a simple mechanism by using a minimal number of parameters in the model, in order to address a variety of experimental observations such as the formation of three germ layers and trophectoderm, responses to altered culture conditions and micropattern diameters and unexpected spotted forms of the germ layers under certain conditions. Moreover, we have tested different boundary conditions as well as various shapes, observing that the pattern is initiated from the boundary and gradually spreads towards the center. This model provides a basis for in-silico modeling of self-organization.

AVAILABILITY AND IMPLEMENTATION

https://github.com/HFooladi/Self_Organization.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Single-sample landscape entropy reveals the imminent phase transition during disease progression

Motivation

The time evolution or dynamic change of many biological systems during disease progression is not always smooth but occasionally abrupt, that is, there is a tipping point during such a process at which the system state shifts from the normal state to a disease state. It is challenging to predict such disease state with the measured omics data, in particular when only a single sample is available.

Results

In this study, we developed a novel approach, i.e. single-sample landscape entropy (SLE) method, to identify the tipping point during disease progression with only one sample data. Specifically, by evaluating the disorder of a network projected from a single-sample data, SLE effectively characterizes the criticality of this single sample network in terms of network entropy, thereby capturing not only the signals of the impending transition but also its leading network, i.e. dynamic network biomarkers. Using this method, we can characterize sample-specific state during disease progression and thus achieve the disease prediction of each individual by only one sample. Our method was validated by successfully identifying the tipping points just before the serious disease symptoms from four real datasets of individuals or subjects, including influenza virus infection, lung cancer metastasis, prostate cancer and acute lung injury.

Availability and implementation

https://github.com/rabbitpei/SLE.

Supplementary information

Supplementary data are available at Bioinformatics online.

Taking Systems Medicine to Heart

Systems medicine is a holistic approach to deciphering the complexity of human physiology in health and disease. In essence, a living body is constituted of networks of dynamically interacting units (molecules, cells, organs, etc) that underlie its collective functions. Declining resilience because of aging and other chronic environmental exposures drives the system to transition from a health state to a disease state; these transitions, triggered by acute perturbations or chronic disturbance, manifest as qualitative shifts in the interactions and dynamics of the disease-perturbed networks. Understanding health-to-disease transitions poses a high-dimensional nonlinear reconstruction problem that requires deep understanding of biology and innovation in study design, technology, and data analysis. With a focus on the principles of systems medicine, this Review discusses approaches for deciphering this biological complexity from a novel perspective, namely, understanding how disease-perturbed networks function; their study provides insights into fundamental disease mechanisms. The immediate goals for systems medicine are to identify early transitions to cardiovascular (and other chronic) diseases and to accelerate the translation of new preventive, diagnostic, or therapeutic targets into clinical practice, a critical step in the development of personalized, predictive, preventive, and participatory (P4) medicine.

Critical Transitions in Intensive Care Units: A Sepsis Case Study

The progression of complex human diseases is associated with critical transitions across dynamical regimes. These transitions often spawn early-warning signals and provide insights into the underlying disease-driving mechanisms. In this paper, we propose a computational method based on surprise loss (SL) to discover data-driven indicators of such transitions in a multivariate time series dataset of septic shock and non-sepsis patient cohorts (MIMIC-III database). The core idea of SL is to train a mathematical model on time series in an unsupervised fashion and to quantify the deterioration of the model’s forecast (out-of-sample) performance relative to its past (in-sample) performance. Considering the highest value of the moving average of SL as a critical transition, our retrospective analysis revealed that critical transitions occurred at a median of over 35 hours before the onset of septic shock, which suggests the applicability of our method as an early-warning indicator. Furthermore, we show that clinical variables at critical-transition regions are significantly different between septic shock and non-sepsis cohorts. Therefore, our paper contributes a critical-transition-based data-sampling strategy that can be utilized for further analysis, such as patient classification. Moreover, our method outperformed other indicators of critical transition in complex systems, such as temporal autocorrelation and variance.

Lateral Roots: Random Diversity in Adversity

Lateral roots are essential for soil foraging and uptake of minerals and water. They feature a large morphological diversity that results from divergent primordia or root growth and development patterns. Besides a structured diversity, resulting from the hierarchical and developmental organization of root systems, there exists a random diversity, occurring between roots of similar age, of the same hierarchical order, and exposed to uniform conditions. The physiological bases and functional consequences of this random diversity are largely ignored. Here we review the evidence for such random diversity throughout the plant kingdom, present innovative approaches based on statistical modeling to account for such diversity, and set the list of its potential benefits in front of a variable and unpredictable soil environment.

Challenges in measuring and understanding biological noise

Biochemical reactions are intrinsically stochastic, leading to variation in the production of mRNAs and proteins within cells. In the scientific literature, this source of variation is typically referred to as 'noise'. The observed variability in molecular phenotypes arises from a combination of processes that amplify and attenuate noise. Our ability to quantify cell-to-cell variability in numerous biological contexts has been revolutionized by recent advances in single-cell technology, from imaging approaches through to 'omics' strategies. However, defining, accurately measuring and disentangling the stochastic and deterministic components of cell-to-cell variability is challenging. In this Review, we discuss the sources, impact and function of molecular phenotypic variability and highlight future directions to understand its role.

Embracing Systems Toxicology at Single-Cell Resolution

As systems biology expands its multi-omic spectrum to increasing resolutions, distinguishing cells based on single-cell profiles becomes feasible. Unlike traditional bulk assays that average cellular responses and blur the distinct identities of responsive cells, single-cell technologies enable sensitive detection of small cellular changes and precise identification of those cells perturbed by toxicants. Among the suite of omic technologies that continue to expand and become affordable, single-cell RNA sequencing (scRNA-seq) is at the cutting edge and leading the way to transform systems toxicology. Single-cell systems toxicology can provide a wealth of information to elucidate cell-specific alterations and response trajectories, detect points-of-departure, map and develop dynamical models of toxicity pathways.

Stem Cell Differentiation as a Non-Markov Stochastic Process

Pluripotent stem cells can self-renew in culture and differentiate along all somatic lineages in vivo. While much is known about the molecular basis of pluripotency, the mechanisms of differentiation remain unclear. Here, we profile individual mouse embryonic stem cells as they progress along the neuronal lineage. We observe that cells pass from the pluripotent state to the neuronal state via an intermediate epiblast-like state. However, analysis of the rate at which cells enter and exit these observed cell states using a hidden Markov model indicates the presence of a chain of unobserved molecular states that each cell transits through stochastically in sequence. This chain of hidden states allows individual cells to record their position on the differentiation trajectory, thereby encoding a simple form of cellular memory. We suggest a statistical mechanics interpretation of these results that distinguishes between functionally distinct cellular “macrostates” and functionally similar molecular “microstates” and propose a model of stem cell differentiation as a non-Markov stochastic process.

•

We profile individual stem cells as they differentiate along the neural lineage

•

Regulatory network changes and increased cell variability accompany differentiation

•

Analysis of dynamics with a hidden Markov model reveals unobserved molecular states

•

We propose a model of stem cell differentiation as a non-Markov stochastic process

Morphological State Transition Dynamics in EGF-Induced Epithelial to Mesenchymal Transition

Epithelial to Mesenchymal Transition (EMT) is a multi-state process. Here, we investigated phenotypic state transition dynamics of Epidermal Growth Factor (EGF)-induced EMT in a breast cancer cell line MDA-MB-468. We have defined phenotypic states of these cells in terms of their morphologies and have shown that these cells have three distinct morphological states-cobble, spindle, and circular. The spindle and circular states are the migratory phenotypes. Using quantitative image analysis and mathematical modeling, we have deciphered state transition trajectories in different experimental conditions. This analysis shows that the phenotypic state transition during EGF-induced EMT in these cells is reversible, and depends upon the dose of EGF and level of phosphorylation of the EGF receptor (EGFR). The dominant reversible state transition trajectory in this system was cobble to circular to spindle to cobble. We have observed that there exists an ultrasensitive on/off switch involving phospho-EGFR that decides the transition of cells in and out of the circular state. In general, our observations can be explained by the conventional quasi-potential landscape model for phenotypic state transition. As an alternative to this model, we have proposed a simpler discretized energy-level model to explain the observed state transition dynamics.

Inherency of Form and Function in Animal Development and Evolution

I discuss recent work on the origins of morphology and cell-type diversification in Metazoa – collectively the animals – and propose a scenario for how these two properties became integrated, with the help of a third set of processes, cellular pattern formation, into the developmental programs seen in present-day metazoans. Inherent propensities to generate familiar forms and cell types, in essence a parts kit for the animals, are exhibited by present-day organisms and were likely more prominent in primitive ones. The structural motifs of animal bodies and organs, e.g., multilayered, hollow, elongated and segmented tissues, internal and external appendages, branched tubes, and modular endoskeletons, can be accounted for by the properties of mesoscale masses of metazoan cells. These material properties, in turn, resulted from the recruitment of “generic” physical forces and mechanisms – adhesion, contraction, polarity, chemical oscillation, diffusion – by toolkit molecules that were partly conserved from unicellular holozoan antecedents and partly novel, distributed in the different metazoan phyla in a fashion correlated with morphological complexity. The specialized functions of the terminally differentiated cell types in animals, e.g., contraction, excitability, barrier function, detoxification, excretion, were already present in ancestral unicellular organisms. These functions were implemented in metazoan differentiation in some cases using the same transcription factors as in single-celled ancestors, although controlled by regulatory mechanisms that were hybrids between earlier-evolved processes and regulatory innovations, such as enhancers. Cellular pattern formation, mediated by released morphogens interacting with biochemically responsive and excitable tissues, drew on inherent self-organizing processes in proto-metazoans to transform clusters of holozoan cells into animal embryos and organs.

The Cognitive Lens: a primer on conceptual tools for analysing information processing in developmental and regenerative morphogenesis

Brains exhibit plasticity, multi-scale integration of information, computation and memory, having evolved by specialization of non-neural cells that already possessed many of the same molecular components and functions. The emerging field of basal cognition provides many examples of decision-making throughout a wide range of non-neural systems. How can biological information processing across scales of size and complexity be quantitatively characterized and exploited in biomedical settings? We use pattern regulation as a context in which to introduce the Cognitive Lens-a strategy using well-established concepts from cognitive and computer science to complement mechanistic investigation in biology. To facilitate the assimilation and application of these approaches across biology, we review tools from various quantitative disciplines, including dynamical systems, information theory and least-action principles. We propose that these tools can be extended beyond neural settings to predict and control systems-level outcomes, and to understand biological patterning as a form of primitive cognition. We hypothesize that a cognitive-level information-processing view of the functions of living systems can complement reductive perspectives, improving efficient top-down control of organism-level outcomes. Exploration of the deep parallels across diverse quantitative paradigms will drive integrative advances in evolutionary biology, regenerative medicine, synthetic bioengineering, cognitive neuroscience and artificial intelligence. This article is part of the theme issue 'Liquid brains, solid brains: How distributed cognitive architectures process information'.

Phenotypic heterogeneity and evolution of melanoma cells associated with targeted therapy resistance

Phenotypic plasticity is associated with non-genetic drug tolerance in several cancers. Such plasticity can arise from chromatin remodeling, transcriptomic reprogramming, and/or protein signaling rewiring, and is characterized as a cell state transition in response to molecular or physical perturbations. This, in turn, can confound interpretations of drug responses and resistance development. Using BRAF-mutant melanoma cell lines as the prototype, we report on a joint theoretical and experimental investigation of the cell-state transition dynamics associated with BRAF inhibitor drug tolerance. Thermodynamically motivated surprisal analysis of transcriptome data was used to treat the cell population as an entropy maximizing system under the influence of time-dependent constraints. This permits the extraction of an epigenetic potential landscape for drug-induced phenotypic evolution. Single-cell flow cytometry data of the same system were modeled with a modified Fokker-Planck-type kinetic model. The two approaches yield a consistent picture that accounts for the phenotypic heterogeneity observed over the course of drug tolerance development. The results reveal that, in certain plastic cancers, the population heterogeneity and evolution of cell phenotypes may be understood by accounting for the competing interactions of the epigenetic potential landscape and state-dependent cell proliferation. Accounting for such competition permits accurate, experimentally verifiable predictions that can potentially guide the design of effective treatment strategies.

Cancer cells exhibit varied degrees of phenotypic heterogeneity. These phenotypes, each of them with unique molecular and functional profiles, display dynamic interconversion in response to drug perturbations, and can evolve to form new drug-tolerant phenotypes. Such phenotypic plasticity, in turn, renders tumor cells extremely difficult to treat. To get a quantitative biophysical understanding of the origins of the phenotypic equilibrium and evolution associated with drug tolerance development in highly plastic patient-derived melanoma cells, we employed joint experimental and computational approaches, using either bulk or single cell measurements as input, to interrogate the epigenetic landscape of the phenotypic evolution. We found that the observed phenotypic equilibria were established via competition between state-dependent net proliferation rates and landscape potential. The results reveal how the tumor cells maintain a phenotypic heterogeneity that facilitates appropriate responses to external cues. They implicate that, in certain phenotypically plastic tumor cells, drug targeting the driver oncogenes may not have sustained efficacy unless the phenotypic plasticity of the tumor is co-targeted.

The Hidden Memory of Differentiating Cells

An integrated approach unifies experimental observations and mathematical modeling to represent differentiation dynamics as discrete transition events underpinned by stochastic transitions between hidden states.

Identification of an early cell fate regulator by detecting dynamics in transcriptional heterogeneity and co-regulation during astrocyte differentiation

There are an increasing number of reports that characterize the temporal behavior of gene expression at the single-cell level during cell differentiation. Despite accumulation of data describing the heterogeneity of biological responses, the dynamics of gene expression heterogeneity and its regulation during the differentiation process have not been studied systematically. To understand transcriptional heterogeneity during astrocyte differentiation, we analyzed single-cell transcriptional data from cells representing the different stages of astrocyte differentiation. When we compared the transcriptional variability of co-expressed genes between the undifferentiated and differentiated states, we found that there was significant increase in transcriptional variability in the undifferentiated state. The genes showing large changes in both "variability" and "correlation" between neural stem cells (NSCs) and astrocytes were found to be functionally involved in astrocyte differentiation. We determined that these genes are potentially regulated by Ascl1, a previously known oscillatory gene in NSCs. Pharmacological blockade of Ntsr2, which is transcriptionally co-regulated with Ascl1, showed that Ntsr2 may play an important role in the differentiation from NSCs to astrocytes. This study shows the importance of characterizing transcriptional heterogeneity and rearrangement of the co-regulation network between different cell states. It also highlights the potential for identifying novel regulators of cell differentiation that will further increase our understanding of the molecular mechanisms underlying the differentiation process.

What Differentiates Poor and Good Outcome Psychotherapy? A Statistical-Mechanics-Inspired Approach to Psychotherapy Research

Statistical mechanics investigates how emergent properties of macroscopic systems (such as temperature and pressure) relate to microscopic state fluctuations. The underlying idea is that global statistical descriptors of order and variability can monitor the relevant dynamics of the whole system at hand. Here we test the possibility of extending such an approach to psychotherapy research investigating the possibility of predicting the outcome of psychotherapy on the sole basis of coarse-grained empirical macro-parameters. Four good-outcome and four poor-outcome brief psychotherapies were recorded, and their transcripts coded in terms of standard psychological categories (abstract, positive emotional and negative emotional language pertaining to patient and therapist). Each patient-therapist interaction is considered as a discrete multivariate time series made of subsequent word-blocks of 150-word length, defined in terms of the above categories. “Static analyses” (Principal Component Analysis) highlighted a substantial difference between good-outcome and poor-outcome cases in terms of mutual correlations among those descriptors. In the former, the patient’s use of abstract language correlated with therapist’s emotional negative language, while in the latter it co-varied with therapist’s emotional positive language, thus showing the different judgment of the therapists regarding the same variable (abstract language) in poor and good outcome cases. On the other hand, the “dynamic analyses”, based on five coarse-grained descriptors related to variability, the degree of order and complexity of the series, demonstrated a relevant case-specific effect, pointing to the possibility of deriving a consistent picture of any single psychotherapeutic process. Overall, the results showed that the systemic approach to psychotherapy (an old tenet of psychology) is mature enough to shift from a metaphorical to a fully quantitative status.

Transition state dynamics during a stochastic fate choice

The generation of multiple fates from a uniform cell population via self-organisation is a recurring feature in development and regeneration. However, for most self-organising systems, we have little understanding of the processes that allow cells to become different. One of the clearest examples of developmental self-organisation is shown by Dictyostelium, with cells segregating into two major fates, stalk and spore, within multicellular aggregates. To characterise the gene expression decisions that underlie this cell fate bifurcation, we carried out single cell transcriptomics on Dictyostelium aggregates. Our data show the transition of progenitors into prespore and prestalk cells occurs via distinct developmental intermediates. Few cells were captured switching between states, with minimal overlap in fate marker expression between cell types, suggesting states are discrete and transitions rapid. Surprisingly, fate-specific transcript dynamics were a small proportion of overall gene expression changes, with transcript divergence coinciding precisely with large-scale remodelling of the transcriptome shared by prestalk and prespore cells. These observations suggest the stepwise separation of cell identity is temporally coupled to global expression transitions common to both fates.

Identifying Critical State of Complex Diseases by Single-Sample-Based Hidden Markov Model

The progression of complex diseases is generally divided as a normal state, a pre-disease state or tipping point, and a disease state. Developing individual-specific method that can identify the pre-disease state just before a catastrophic deterioration, is critical for patients with complex diseases. However, with only a case sample, it is challenging to detect a pre-disease state which has little significant differences comparing with a normal state in terms of phenotypes and gene expressions. In this study, by regarding the tipping point as the end point of a stationary Markov process, we proposed a single-sample-based hidden Markov model (HMM) approach to explore the dynamical differences between a normal and a pre-disease states, and thus can signal the upcoming critical transition immediately after a pre-disease state. Using this method, we identified the pre-disease state or tipping point in a numerical simulation and two real datasets including stomach adenocarcinoma and influenza infection, which demonstrate the effectiveness of the method.

Noise-induced tipping under periodic forcing: Preferred tipping phase in a non-adiabatic forcing regime

We consider a periodically forced 1D Langevin equation that possesses two stable periodic solutions in the absence of noise. We ask the question: is there a most likely noise-induced transition path between these periodic solutions that allows us to identify a preferred phase of the forcing when tipping occurs? The quasistatic regime, where the forcing period is long compared to the adiabatic relaxation time, has been well studied; our work instead explores the case when these time scales are comparable. We compute optimal paths using the path integral method incorporating the Onsager-Machlup functional and validate results with Monte Carlo simulations. Results for the preferred tipping phase are compared with the deterministic aspects of the problem. We identify parameter regimes where nullclines, associated with the deterministic problem in a 2D extended phase space, form passageways through which the optimal paths transit. As the nullclines are independent of the relaxation time and the noise strength, this leads to a robust deterministic predictor of the preferred tipping phase in a regime where forcing is neither too fast nor too slow.

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.

The Impact of Heterogeneity on Single-Cell Sequencing

The importance of diversity and cellular specialization is clear for many reasons, from population-level diversification, to improved resiliency to unforeseen stresses, to unique functions within metazoan organisms during development and differentiation. However, the level of cellular heterogeneity is just now becoming clear through the integration of genome-wide analyses and more cost effective Next Generation Sequencing (NGS). With easy access to single-cell NGS (scNGS), new opportunities exist to examine different levels of gene expression and somatic mutational heterogeneity, but these assays can generate yottabyte scale data. Here, we model the importance of heterogeneity for large-scale analysis of scNGS data, with a focus on the utilization in oncology and other diseases, providing a guide to aid in sample size and experimental design.

Towards a Quantitative Understanding of Cell Identity

Cells have traditionally been characterized using expression levels of a few proteins that are thought to specify phenotype. This requires a priori selection of proteins, which can introduce descriptor bias, and neglects the wealth of additional molecular information nested within each cell in a population, which often makes these sparse descriptors qualitative. Recently, more unbiased and quantitative cell characterization has been made possible by new high-throughput, information-dense experimental approaches and data-driven computational methods. This review discusses such quantitative descriptors in the context of three central concepts of cell identity: definition, creation, and stability. Collectively, these concepts are essential for constructing quantitative phenotypic landscapes, which will enhance our understanding of cell biology and facilitate cell engineering for research and clinical applications.

Single cell RNA-seq and ATAC-seq analysis of cardiac progenitor cell transition states and lineage settlement

Formation and segregation of cell lineages forming the heart have been studied extensively but the underlying gene regulatory networks and epigenetic changes driving cell fate transitions during early cardiogenesis are still only partially understood. Here, we comprehensively characterize mouse cardiac progenitor cells (CPCs) marked by Nkx2-5 and Isl1 expression from E7.5 to E9.5 using single-cell RNA sequencing and transposase-accessible chromatin profiling (ATAC-seq). By leveraging on cell-to-cell transcriptome and chromatin accessibility heterogeneity, we identify different previously unknown cardiac subpopulations. Reconstruction of developmental trajectories reveal that multipotent Isl1+ CPC pass through an attractor state before separating into different developmental branches, whereas extended expression of Nkx2-5 commits CPC to an unidirectional cardiomyocyte fate. Furthermore, we show that CPC fate transitions are associated with distinct open chromatin states critically depending on Isl1 and Nkx2-5. Our data provide a model of transcriptional and epigenetic regulations during cardiac progenitor cell fate decisions at single-cell resolution.

Computational Modeling and Reverse Engineering to Reveal Dominant Regulatory Interactions Controlling Osteochondral Differentiation: Potential for Regenerative Medicine

The specialization of cartilage cells, or chondrogenic differentiation, is an intricate and meticulously regulated process that plays a vital role in both bone formation and cartilage regeneration. Understanding the molecular regulation of this process might help to identify key regulatory factors that can serve as potential therapeutic targets, or that might improve the development of qualitative and robust skeletal tissue engineering approaches. However, each gene involved in this process is influenced by a myriad of feedback mechanisms that keep its expression in a desirable range, making the prediction of what will happen if one of these genes defaults or is targeted with drugs, challenging. Computer modeling provides a tool to simulate this intricate interplay from a network perspective. This paper aims to give an overview of the current methodologies employed to analyze cell differentiation in the context of skeletal tissue engineering in general and osteochondral differentiation in particular. In network modeling, a network can either be derived from mechanisms and pathways that have been reported in the literature (knowledge-based approach) or it can be inferred directly from the data (data-driven approach). Combinatory approaches allow further optimization of the network. Once a network is established, several modeling technologies are available to interpret dynamically the relationships that have been put forward in the network graph (implication of the activation or inhibition of certain pathways on the evolution of the system over time) and to simulate the possible outcomes of the established network such as a given cell state. This review provides for each of the aforementioned steps (building, optimizing, and modeling the network) a brief theoretical perspective, followed by a concise overview of published works, focusing solely on applications related to cell fate decisions, cartilage differentiation and growth plate biology. Particular attention is paid to an in-house developed example of gene regulatory network modeling of growth plate chondrocyte differentiation as all the aforementioned steps can be illustrated. In summary, this paper discusses and explores a series of tools that form a first step toward a rigorous and systems-level modeling of osteochondral differentiation in the context of regenerative medicine.

Detecting early-warning signals of influenza outbreak based on dynamic network marker

The seasonal outbreaks of influenza infection cause globally respiratory illness, or even death in all age groups. Given early‐warning signals preceding the influenza outbreak, timely intervention such as vaccination and isolation management effectively decrease the morbidity. However, it is usually a difficult task to achieve the real‐time prediction of influenza outbreak due to its complexity intertwining both biological systems and social systems. By exploring rich dynamical and high‐dimensional information, our dynamic network marker/biomarker (DNM/DNB) method opens a new way to identify the tipping point prior to the catastrophic transition into an influenza pandemics. In order to detect the early‐warning signals before the influenza outbreak by applying DNM method, the historical information of clinic hospitalization caused by influenza infection between years 2009 and 2016 were extracted and assembled from public records of Tokyo and Hokkaido, Japan. The early‐warning signal, with an average of 4‐week window lead prior to each seasonal outbreak of influenza, was provided by DNM‐based on the hospitalization records, providing an opportunity to apply proactive strategies to prevent or delay the onset of influenza outbreak. Moreover, the study on the dynamical changes of hospitalization in local district networks unveils the influenza transmission dynamics or landscape in network level.

Histone Acetyltransferase KAT2A Stabilizes Pluripotency with Control of Transcriptional Heterogeneity

Cell fate transitions in mammalian stem cell systems have often been associated with transcriptional heterogeneity; however, existing data have failed to establish a functional or mechanistic link between the two phenomena. Experiments in unicellular organisms support the notion that transcriptional heterogeneity can be used to facilitate adaptability to environmental changes and have identified conserved chromatin‐associated factors that modulate levels of transcriptional noise. Herein, we show destabilization of pluripotency‐associated gene regulatory networks through increased transcriptional heterogeneity of mouse embryonic stem cells in which paradigmatic histone acetyl‐transferase, and candidate noise modulator, Kat2a (yeast orthologue Gcn5), have been inhibited. Functionally, network destabilization associates with reduced pluripotency and accelerated mesendodermal differentiation, with increased probability of transitions into lineage commitment. Thus, we show evidence of a relationship between transcriptional heterogeneity and cell fate transitions through manipulation of the histone acetylation landscape of mouse embryonic stem cells, suggesting a general principle that could be exploited in other normal and malignant stem cell fate transitions. stem cells2018;36:1828–11

Thermodynamic Analysis of Time Evolving Networks

The problem of how to represent networks, and from this representation, derive succinct characterizations of network structure and in particular how this structure evolves with time, is of central importance in complex network analysis. This paper tackles the problem by proposing a thermodynamic framework to represent the structure of time-varying complex networks. More importantly, such a framework provides a powerful tool for better understanding the network time evolution. Specifically, the method uses a recently-developed approximation of the network von Neumann entropy and interprets it as the thermodynamic entropy for networks. With an appropriately-defined internal energy in hand, the temperature between networks at consecutive time points can be readily derived, which is computed as the ratio of change of entropy and change in energy. It is critical to emphasize that one of the main advantages of the proposed method is that all these thermodynamic variables can be computed in terms of simple network statistics, such as network size and degree statistics. To demonstrate the usefulness of the thermodynamic framework, the paper uses real-world network data, which are extracted from time-evolving complex systems in the financial and biological domains. The experimental results successfully illustrate that critical events, including abrupt changes and distinct periods in the evolution of complex networks, can be effectively characterized.

Correcting the Mean-Variance Dependency for Differential Variability Testing Using Single-Cell RNA Sequencing Data

Cell-to-cell transcriptional variability in otherwise homogeneous cell populations plays an important role in tissue function and development. Single-cell RNA sequencing can characterize this variability in a transcriptome-wide manner. However, technical variation and the confounding between variability and mean expression estimates hinder meaningful comparison of expression variability between cell populations. To address this problem, we introduce an analysis approach that extends the BASiCS statistical framework to derive a residual measure of variability that is not confounded by mean expression. This includes a robust procedure for quantifying technical noise in experiments where technical spike-in molecules are not available. We illustrate how our method provides biological insight into the dynamics of cell-to-cell expression variability, highlighting a synchronization of biosynthetic machinery components in immune cells upon activation. In contrast to the uniform up-regulation of the biosynthetic machinery, CD4+ T cells show heterogeneous up-regulation of immune-related and lineage-defining genes during activation and differentiation.

Bridging the Data Gap From in vitro Toxicity Testing to Chemical Safety Assessment Through Computational Modeling

Chemical toxicity testing is moving steadily toward a human cell and organoid-based in vitro approach for reasons including scientific relevancy, efficiency, cost, and ethical rightfulness. Inferring human health risk from chemical exposure based on in vitro testing data is a challenging task, facing various data gaps along the way. This review identifies these gaps and makes a case for the in silico approach of computational dose-response and extrapolation modeling to address many of the challenges. Mathematical models that can mechanistically describe chemical toxicokinetics (TK) and toxicodynamics (TD), for both in vitro and in vivo conditions, are the founding pieces in this regard. Identifying toxicity pathways and in vitro point of departure (PoD) associated with adverse health outcomes requires an understanding of the molecular key events in the interacting transcriptome, proteome, and metabolome. Such an understanding will in turn help determine the sets of sensitive biomarkers to be measured in vitro and the scope of toxicity pathways to be modeled in silico. In vitro data reporting both pathway perturbation and chemical biokinetics in the culture medium serve to calibrate the toxicity pathway and virtual tissue models, which can then help predict PoDs in response to chemical dosimetry experienced by cells in vivo. Two types of in vitro to in vivo extrapolation (IVIVE) are needed. (1) For toxic effects involving systemic regulations, such as endocrine disruption, organism-level adverse outcome pathway (AOP) models are needed to extrapolate in vitro toxicity pathway perturbation to in vivo PoD. (2) Physiologically-based toxicokinetic (PBTK) modeling is needed to extrapolate in vitro PoD dose metrics into external doses for expected exposure scenarios. Linked PBTK and TD models can explore the parameter space to recapitulate human population variability in response to chemical insults. While challenges remain for applying these modeling tools to support in vitro toxicity testing, they open the door toward population-stratified and personalized risk assessment.

Transition state characteristics during cell differentiation

Models describing the process of stem-cell differentiation are plentiful, and may offer insights into the underlying mechanisms and experimentally observed behaviour. Waddington's epigenetic landscape has been providing a conceptual framework for differentiation processes since its inception. It also allows, however, for detailed mathematical and quantitative analyses, as the landscape can, at least in principle, be related to mathematical models of dynamical systems. Here we focus on a set of dynamical systems features that are intimately linked to cell differentiation, by considering exemplar dynamical models that capture important aspects of stem cell differentiation dynamics. These models allow us to map the paths that cells take through gene expression space as they move from one fate to another, e.g. from a stem-cell to a more specialized cell type. Our analysis highlights the role of the transition state (TS) that separates distinct cell fates, and how the nature of the TS changes as the underlying landscape changes-change that can be induced by e.g. cellular signaling. We demonstrate that models for stem cell differentiation may be interpreted in terms of either a static or transitory landscape. For the static case the TS represents a particular transcriptional profile that all cells approach during differentiation. Alternatively, the TS may refer to the commonly observed period of heterogeneity as cells undergo stochastic transitions.

Identification of Biologically Essential Nodes via Determinative Power in Logical Models of Cellular Processes

A variety of biological networks can be modeled as logical or Boolean networks. However, a simplification of the reality to binary states of the nodes does not ease the difficulty of analyzing the dynamics of large, complex networks, such as signal transduction networks, due to the exponential dependence of the state space on the number of nodes. This paper considers a recently introduced method for finding a fairly small subnetwork, representing a collection of nodes that determine the states of most other nodes with a reasonable level of entropy. The subnetwork contains the most determinative nodes that yield the highest information gain. One of the goals of this paper is to propose an algorithm for finding a suitable subnetwork size. The information gain is quantified by the so-called determinative power of the nodes, which is obtained via the mutual information, a concept originating in information theory. We find the most determinative nodes for 36 network models available in the online database Cell Collective (http://cellcollective.org). We provide statistical information that indicates a weak correlation between the subnetwork size and other variables, such as network size, or maximum and average determinative power of nodes. We observe that the proportion represented by the subnetwork in comparison to the whole network shows a weak tendency to decrease for larger networks. The determinative power of nodes is weakly correlated to the number of outputs of a node, and it appears to be independent of other topological measures such as closeness or betweenness centrality. Once the subnetwork of the most determinative nodes is identified, we generate a biological function analysis of its nodes for some of the 36 networks. The analysis shows that a large fraction of the most determinative nodes are essential and involved in crucial biological functions. The biological pathway analysis of the most determinative nodes shows that they are involved in important disease pathways.

Determining Relative Dynamic Stability of Cell States Using Boolean Network Model

Cell state transition is at the core of biological processes in metazoan, which includes cell differentiation, epithelial-to-mesenchymal transition (EMT) and cell reprogramming. In these cases, it is important to understand the molecular mechanism of cellular stability and how the transitions happen between different cell states, which is controlled by a gene regulatory network (GRN) hard-wired in the genome. Here we use Boolean modeling of GRN to study the cell state transition of EMT and systematically compare four available methods to calculate the cellular stability of three cell states in EMT in both normal and genetically mutated cases. The results produced from four methods generally agree but do not totally agree with each other. We show that distribution of one-degree neighborhood of cell states, which are the nearest states by Hamming distance, causes the difference among the methods. From that, we propose a new method based on one-degree neighborhood, which is the simplest one and agrees with other methods to estimate the cellular stability in all scenarios of our EMT model. This new method will help the researchers in the field of cell differentiation and cell reprogramming to calculate cellular stability using Boolean model, and then rationally design their experimental protocols to manipulate the cell state transition.

Conceptual Challenges of the Systemic Approach in Understanding Cell Differentiation

The cells of a multicellular organism are derived from a single zygote and genetically identical. Yet, they are phenotypically very different. This difference is the result of a process commonly called cell differentiation. How the phenotypic diversity emerges during ontogenesis or regeneration is a central and intensely studied but still unresolved issue in biology. Cell biology is facing conceptual challenges that are frequently confused with methodological difficulties. How to define a cell type? What stability or change means in the context of cell differentiation and how to deal with the ubiquitous molecular variations seen in the living cells? What are the driving forces of the change? We propose to reframe the problem of cell differentiation in a systemic way by incorporating different theoretical approaches. The new conceptual framework is able to capture the insights made at different levels of cellular organization and considered previously as contradictory. It also provides a formal strategy for further experimental studies.

scEpath: energy landscape-based inference of transition probabilities and cellular trajectories from single-cell transcriptomic data

Motivation

Single-cell RNA-sequencing (scRNA-seq) offers unprecedented resolution for studying cellular decision-making processes. Robust inference of cell state transition paths and probabilities is an important yet challenging step in the analysis of these data.

Results

Here we present scEpath, an algorithm that calculates energy landscapes and probabilistic directed graphs in order to reconstruct developmental trajectories. We quantify the energy landscape using 'single-cell energy' and distance-based measures, and find that the combination of these enables robust inference of the transition probabilities and lineage relationships between cell states. We also identify marker genes and gene expression patterns associated with cell state transitions. Our approach produces pseudotemporal orderings that are-in combination-more robust and accurate than current methods, and offers higher resolution dynamics of the cell state transitions, leading to new insight into key transition events during differentiation and development. Moreover, scEpath is robust to variation in the size of the input gene set, and is broadly unsupervised, requiring few parameters to be set by the user. Applications of scEpath led to the identification of a cell-cell communication network implicated in early human embryo development, and novel transcription factors important for myoblast differentiation. scEpath allows us to identify common and specific temporal dynamics and transcriptional factor programs along branched lineages, as well as the transition probabilities that control cell fates.

Availability and implementation

A MATLAB package of scEpath is available at https://github.com/sqjin/scEpath.

Supplementary information

Supplementary data are available at Bioinformatics online.

Aging in a Relativistic Biological Space-Time

Here we present a theoretical and mathematical perspective on the process of aging. We extend the concepts of physical space and time to an abstract, mathematically-defined space, which we associate with a concept of “biological space-time” in which biological dynamics may be represented. We hypothesize that biological dynamics, represented as trajectories in biological space-time, may be used to model and study different rates of biological aging. As a consequence of this hypothesis, we show how dilation or contraction of time analogous to relativistic corrections of physical time resulting from accelerated or decelerated biological dynamics may be used to study precipitous or protracted aging. We show specific examples of how these principles may be used to model different rates of aging, with an emphasis on cancer in aging. We discuss how this theory may be tested or falsified, as well as novel concepts and implications of this theory that may improve our interpretation of biological aging.

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.

Network-Based Predictors of Progression in Head and Neck Squamous Cell Carcinoma

The heterogeneity in head and neck squamous cell carcinoma (HNSCC) has made reliable stratification extremely challenging. Behavioral risk factors such as smoking and alcohol consumption contribute to this heterogeneity. To help elucidate potential mechanisms of progression in HNSCC, we focused on elucidating patterns of gene interactions associated with tumor progression. We performed de-novo gene co-expression network inference utilizing 229 patient samples from The Cancer Genome Atlas (TCGA) previously annotated by Bornstein et al. (2016). Differential network analysis allowed us to contrast progressor and non-progressor cohorts. Beyond standard differential expression (DE) analysis, this approach evaluates changes in gene expression variance (differential variability DV) and changes in covariance, which we denote as differential wiring (DW). The set of affected genes was overlaid onto the co-expression network, identifying 12 modules significantly enriched in DE, DV, and/or DW genes. Additionally, we identified modules correlated with behavioral measures such as alcohol consumption and smoking status. In the module enriched for differentially wired genes, we identified network hubs including IL10RA, DOK2, APBB1IP, UBASH3A, SASH3, CELF2, TRAF3IP3, GIMAP6, MYO1F, NCKAP1L, WAS, FERMT3, SLA, SELPLG, TNFRSF1B, WIPF1, AMICA1, PTPN22; the network centrality and progression specificity of these genes suggest a potential role in tumor evolution mechanisms related to inflammation and microenvironment. The identification of this network-based gene signature could be further developed to guide progression stratification, highlighting how network approaches may help improve clinical research end points and ultimately aid in clinical utility.

Modeling of cytometry data in logarithmic space: When is a bimodal distribution not bimodal?

Recent efforts in systems immunology lead researchers to build quantitative models of cell activation and differentiation. One goal is to account for the distributions of proteins from single-cell measurements by flow cytometry or mass cytometry as readout of biological regulation. In that context, large cell-to-cell variability is often observed in biological quantities. We show here that these readouts, viewed in logarithmic scale may result in two easily-distinguishable modes, while the underlying distribution (in linear scale) is unimodal. We introduce a simple mathematical test to highlight this mismatch. We then dissect the flow of influence of cell-to-cell variability proposing a graphical model which motivates higher-dimensional analysis of the data. Finally we show how acquiring additional biological information can be used to reduce uncertainty introduced by cell-to-cell variability, helping to clarify whether the data is uni- or bimodal. This communication has cautionary implications for manual and automatic gating strategies, as well as clustering and modeling of single-cell measurements. © 2018 International Society for Advancement of Cytometry.

Dynamic network biomarker indicates pulmonary metastasis at the tipping point of hepatocellular carcinoma

Developing predictive biomarkers that can detect the tipping point before metastasis of hepatocellular carcinoma (HCC), is critical to prevent further irreversible deterioration. To discover such early-warning signals or biomarkers of pulmonary metastasis in HCC, we analyse time-series gene expression data in spontaneous pulmonary metastasis mice HCCLM3-RFP model with our dynamic network biomarker (DNB) method, and identify CALML3 as a core DNB member. All experimental results of gain-of-function and loss-of-function studies show that CALML3 could indicate metastasis initiation and act as a suppressor of metastasis. We also reveal the biological role of CALML3 in metastasis initiation at a network level, including proximal regulation and cascading influences in dysfunctional pathways. Our further experiments and clinical samples show that DNB with CALML3 reduced pulmonary metastasis in liver cancer. Actually, loss of CALML3 predicts shorter overall and relapse-free survival in postoperative HCC patients, thus providing a prognostic biomarker and therapy target in HCC.

Reprogramming cell fate with artificial transcription factors

Transcription factors (TFs) reprogram cell states by exerting control over gene regulatory networks and the epigenetic landscape of a cell. Artificial transcription factors (ATFs) are designer regulatory proteins comprised of modular units that can be customized to overcome challenges faced by natural TFs in establishing and maintaining desired cell states. Decades of research on DNA-binding proteins and synthetic molecules has provided a molecular toolkit for ATF design and the construction of genome-scale libraries of ATFs capable of phenotypic manipulation and reprogramming of cell states. Here, we compare the unique strengths and limitations of different ATF platforms, highlight the advantages of cooperative assembly, and present the potential of ATF libraries in revealing gene regulatory networks that govern cell fate choices.

Automated cell cycle and cell size measurements for single-cell gene expression studies

OBJECTIVES

Recent rise of single-cell studies revealed the importance of understanding the role of cell-to-cell variability, especially at the transcriptomic level. One of the numerous sources of cell-to-cell variation in gene expression is the heterogeneity in cell proliferation state. In order to identify how cell cycle and cell size influences gene expression variability at the single-cell level, we provide an universal and automatic toxic-free label method, compatible with single-cell high-throughput RT-qPCR. The method consists of isolating cells after a double-stained, analyzing their morphological parameters and performing a transcriptomic analysis on the same identified cells.

RESULTS

This led to an unbiased gene expression analysis and could be also used for improving single-cell tracking and imaging when combined with cell isolation. As an application for this technique, we showed that cell-to-cell variability in chicken erythroid progenitors was negligibly influenced by cell size nor cell cycle.