Cell Fate Decision as High-Dimensional Critical State Transition

Single-cell analyses reveal increased gene expression variability in human neurodevelopmental conditions

Interindividual variation in phenotypic penetrance and severity is found in many neurodevelopmental conditions, although the underlying mechanisms remain largely unresolved. Within individuals, homogeneous cell types (i.e., genetically identical and in similar environments) can differ in molecule abundance. Here, we investigate the hypothesis that neurodevelopmental conditions can drive increased variability in gene expression, not just differential gene expression. Leveraging independent single-cell and single-nucleus RNA sequencing datasets derived from human brain-relevant cell and tissue types, we identify a significant increase in gene expression variability driven by the autosomal aneuploidy trisomy 21 (T21) as well as autism-associated chromodomain helicase DNA binding protein 8 (CHD8) haploinsufficiency. Our analyses are consistent with a global and, in part, stochastic increase in variability, which is uncoupled from changes in transcript abundance. Highly variable genes tend to be cell-type specific with modest enrichment for repressive H3K27me3, while least variable genes are more likely to be constrained and associated with active histone marks. Our results indicate that human neurodevelopmental conditions can drive increased gene expression variability in brain cell types, with the potential to contribute to diverse phenotypic outcomes. These findings also provide a scaffold for understanding variability in disease, essential for deeper insights into genotype-phenotype relationships.

Exploring cell-to-cell variability and functional insights through differentially variable gene analysis

Single-cell RNA sequencing (scRNA-seq) has revolutionized our understanding of cellular variability by capturing gene expression profiles of individual cells. The importance of cell-to-cell variability in determining and shaping cell function has been widely appreciated. Nevertheless, differential expression (DE) analysis remains a cornerstone method in analytical practice. Current computational analyses overlook the rich information encoded by variability within the single-cell gene expression data by focusing exclusively on mean expression. To offer a deeper understanding of cellular systems, there is a need for approaches to assess data variability rather than just the mean. Here we present spline-DV, a statistical framework for differential variability (DV) analysis using scRNA-seq data. The spline-DV method identifies genes exhibiting significantly increased or decreased expression variability among cells derived from two experimental conditions. Case studies show that DV genes identified using spline-DV are representative and functionally relevant to tested cellular conditions, including obesity, fibrosis, and cancer.

Fokker-Planck diffusion maps of microglial transcriptomes reveal radial differentiation into substates associated with Alzheimer's pathology

The identification of microglia subtypes is important for understanding the role of innate immunity in neurodegenerative diseases. Current methods of unsupervised cell type identification assume a small noise-to-signal ratio of transcriptome measurements to produce well-separated cell clusters. However, identification of subtypes can be obscured by gene expression noise, which diminishes the distances in transcriptome space between distinct cell types, blurs boundaries, and reduces reproducibility. Here we use Fokker-Planck (FP) diffusion maps to model cellular differentiation as a stochastic process whereby cells settle into local minima that correspond to cell subtypes, in a potential landscape constructed from transcriptome data using a nearest neighbor graph approach. By applying critical transition fields, we identify individual cells on the verge of transitioning between subtypes, revealing microglial cells in an inactivated, homeostatic state before radially transitioning into various specialized subtypes. Specifically, we show that cells from Alzheimer's disease patients are enriched in a microglia subtype associated to antigen presentation and T-cell recruitment, and are depleted in an anti-inflammatory subtype.

Computing hematopoiesis plasticity in response to genetic mutations and environmental stimulations

Cell plasticity (CP), describing a dynamic cell state, plays a crucial role in maintaining homeostasis during organ morphogenesis, regeneration, and trauma-to-repair biological process. Single-cell-omics datasets provide an unprecedented resource to empower CP analysis. Hematopoiesis offers fertile opportunities to develop quantitative methods for understanding CP. In this study, we generated high-quality lineage-negative single-cell RNA-sequencing datasets under various conditions and introduced a working pipeline named scPlasticity to interrogate naïve and disturbed plasticity of hematopoietic stem and progenitor cells with mutational or environmental challenges. Using embedding methods UMAP or FA, a continuum of hematopoietic development is visually observed in wild type where the pipeline confirms a low proportion of hybrid cells ( P hc , with bias range: 0.4∼0.6) on a transition trajectory. Upon Tet2 mutation, a driver of leukemia, or treatment of DSS, an inducer of colitis, P hc is increased and plasticity of hematopoietic stem and progenitor cells was enhanced. We prioritized several transcription factors and signaling pathways, which are responsible for P hc alterations. In silico perturbation suggests knocking out EGR regulons or pathways of IL-1R1 and β-adrenoreceptor partially reverses P hc promoted by Tet2 mutation and inflammation.

Attractor Landscape Analysis Reveals a Reversion Switch in the Transition of Colorectal Tumorigenesis

A cell fate change such as tumorigenesis incurs critical transition. It remains a longstanding challenge whether the underlying mechanism can be unraveled and a molecular switch that can reverse such transition is found. Here a systems framework, REVERT, is presented with which can reconstruct the core molecular regulatory network model and a reversion switch based on single-cell transcriptome data over the transition process is identified. The usefulness of REVERT is demonstrated by applying it to single-cell transcriptome of patient-derived matched organoids of colon cancer and normal colon. REVERT is a generic framework that can be applied to investigate various cell fate transition phenomena.

A statistical approach for systematic identification of transition cells from scRNA-seq data

Decoding cellular state transitions is crucial for understanding complex biological processes in development and disease. While recent advancements in single-cell RNA sequencing (scRNA-seq) offer insights into cellular trajectories, existing tools primarily study expressional rather than regulatory state shifts. We present CellTran, a statistical approach utilizing paired-gene expression correlations to detect transition cells from scRNA-seq data without explicitly resolving gene regulatory networks. Applying our approach to various contexts, including tissue regeneration, embryonic development, preinvasive lesions, and humoral responses post-vaccination, reveals transition cells and their distinct gene expression profiles. Our study sheds light on the underlying molecular mechanisms driving cellular state transitions, enhancing our ability to identify therapeutic targets for disease interventions.

Information structure of heterogeneous criticality in a fish school

Integrated information theory (IIT) assesses the degree of consciousness in living organisms from an information-theoretic perspective. This theory can be generalised to other systems, including those exhibiting criticality. In this study, we applied IIT to the collective behaviour of Plecoglossus altivelis and observed that the group integrity (Φ) was maximised at the critical state. Multiple levels of criticality were identified within the group, existing as distinct subgroups. Moreover, these fragmented critical subgroups coexisted alongside the overall criticality of the group. The distribution of high-criticality subgroups was heterogeneous across both time and space. Notably, core fish in the high-criticality subgroups were less affected by internal and external stimuli compared to those in low-criticality subgroups. These findings are consistent with previous interpretations of critical phenomena and offer a new perspective on the dynamics of an empirical critical state.

The network structural entropy for single-cell RNA sequencing data during skin aging

Aging is a complex and heterogeneous biological process at cellular, tissue, and individual levels. Despite extensive effort in scientific research, a comprehensive understanding of aging mechanisms remains lacking. This study analyzed aging-related gene networks, using single-cell RNA sequencing data from >15 000 cells. We constructed a gene correlation network, integrating gene expressions into the weights of network edges, and ranked gene importance using a random walk model to generate a gene importance matrix. This unsupervised method improved the clustering performance of cell types. To further quantify the complexity of gene networks during aging, we introduced network structural entropy. The findings of our study reveal that the overall network structural entropy increases in the aged cells compared to the young cells. However, network entropy changes varied greatly within different cell subtypes. Specifically, the network structural entropy among various cell types may increase, remain unchanged, or decrease. This wide range of changes may be closely related to their individual functions, highlighting the cellular heterogeneity and potential key network reconfigurations. Analyzing gene network entropy provides insights into the molecular mechanisms behind aging. This study offers new scientific evidence and theoretical support for understanding the changes in cell functions during aging.

Multiscale chromatin dynamics and high entropy in plant iPSC ancestors

Plant protoplasts provide starting material for of inducing pluripotent cell masses that are competent for tissue regeneration in vitro, analogous to animal induced pluripotent stem cells (iPSCs). Dedifferentiation is associated with large-scale chromatin reorganisation and massive transcriptome reprogramming, characterised by stochastic gene expression. How this cellular variability reflects on chromatin organisation in individual cells and what factors influence chromatin transitions during culturing are largely unknown. Here, we used high-throughput imaging and a custom supervised image analysis protocol extracting over 100 chromatin features of cultured protoplasts. The analysis revealed rapid, multiscale dynamics of chromatin patterns with a trajectory that strongly depended on nutrient availability. Decreased abundance in H1 (linker histones) is hallmark of chromatin transitions. We measured a high heterogeneity of chromatin patterns indicating intrinsic entropy as a hallmark of the initial cultures. We further measured an entropy decline over time, and an antagonistic influence by external and intrinsic factors, such as phytohormones and epigenetic modifiers, respectively. Collectively, our study benchmarks an approach to understand the variability and evolution of chromatin patterns underlying plant cell reprogramming in vitro.

Inertial effect of cell state velocity on the quiescence-proliferation fate decision

Energy landscapes can provide intuitive depictions of population heterogeneity and dynamics. However, it is unclear whether individual cell behavior, hypothesized to be determined by initial position and noise, is faithfully recapitulated. Using the p21-/Cdk2-dependent quiescence-proliferation decision in breast cancer dormancy as a testbed, we examined single-cell dynamics on the landscape when perturbed by hypoxia, a dormancy-inducing stress. Combining trajectory-based energy landscape generation with single-cell time-lapse microscopy, we found that a combination of initial position and velocity on a p21/Cdk2 landscape, but not position alone, was required to explain the observed cell fate heterogeneity under hypoxia. This is likely due to additional cell state information such as epigenetic features and/or other species encoded in velocity but missing in instantaneous position determined by p21 and Cdk2 levels alone. Here, velocity dependence manifested as inertia: cells with higher cell cycle velocities prior to hypoxia continued progressing along the cell cycle under hypoxia, resisting the change in landscape towards cell cycle exit. Such inertial effects may markedly influence cell fate trajectories in tumors and other dynamically changing microenvironments where cell state transitions are governed by coordination across several biochemical species.

Mitochondrial genetic variations in leukemia: a comprehensive overview

Leukemias are a group of heterogeneous hematological malignancies driven by diverse genetic variations, and the advent of genomic sequencing technologies facilitates the investigation of genetic abnormalities in leukemia. However, these sequencing-based studies mainly focus on nuclear DNAs. Increasing evidence indicates that mitochondrial dysfunction is an important mechanism of leukemia pathogenesis, which is closely related to the mitochondrial genome variations. Here, we provide an overview of current research progress concerning mitochondrial genetic variations in leukemia, encompassing gene mutations and copy number variations. We also summarize currently accessible mitochondrial DNA (mtDNA) sequencing methods. Notably, somatic mtDNA mutations may serve as natural genetic barcodes for lineage tracing and longitudinal assessment of clonal dynamics. Collectively, these findings enhance our understanding of leukemia pathogenesis and foster the identification of novel therapeutic targets and interventions.

Dissecting reversible and irreversible single cell state transitions from gene regulatory networks

Understanding cell state transitions and their governing regulatory mechanisms remains one of the fundamental questions in biology. We develop a computational method, state transition inference using cross-cell correlations (STICCC), for predicting reversible and irreversible cell state transitions at single-cell resolution by using gene expression data and a set of gene regulatory interactions. The method is inspired by the fact that the gene expression time delays between regulators and targets can be exploited to infer past and future gene expression states. From applications to both simulated and experimental single-cell gene expression data, we show that STICCC-inferred vector fields capture basins of attraction and irreversible fluxes. By connecting regulatory information with systems' dynamical behaviors, STICCC reveals how network interactions influence reversible and irreversible state transitions. Compared to existing methods that infer pseudotime and RNA velocity, STICCC provides complementary insights into the gene regulation of cell state transitions.

Form, function, mind: What doesn't compute (and what might)

The applicability of computational and dynamical systems models to organisms is scrutinized, using examples from developmental biology and cognition. Developmental morphogenesis is dependent on the inherent material properties of developing animal (metazoan) tissues, a non-computational modality, but cell differentiation, which utilizes chromatin-based revisable memory banks and program-like function-calling, via the developmental gene co-expression system unique to the metazoans, has a quasi-computational basis. Multi-attractor dynamical models are argued to be misapplied to global properties of development, and it is suggested that along with computationalism, classic forms of dynamicism are similarly unsuitable to accounting for cognitive phenomena. Proposals are made for treating brains and other nervous tissues as novel forms of excitable matter with inherent properties which enable the intensification of cell-based basal cognition capabilities present throughout the tree of life. Finally, some connections are drawn between the viewpoint described here and active inference models of cognition, such as the Free Energy Principle.

Uncovering underlying physical principles and driving forces of cell differentiation and reprogramming from single-cell transcriptomics

Recent advances in single-cell sequencing technology have revolutionized our ability to acquire whole transcriptome data. However, uncovering the underlying transcriptional drivers and nonequilibrium driving forces of cell function directly from these data remains challenging. We address this by learning cell state vector fields from discrete single-cell RNA velocity to quantify the single-cell global nonequilibrium driving forces as landscape and flux. From single-cell data, we quantified the Waddington landscape, showing that optimal paths for differentiation and reprogramming deviate from the naively expected landscape gradient paths and may not pass through landscape saddles at finite fluctuations, challenging conventional transition state estimation of kinetic rate for cell fate decisions due to the presence of the flux. A key insight from our study is that stem/progenitor cells necessitate greater energy dissipation for rapid cell cycles and self-renewal, maintaining pluripotency. We predict optimal developmental pathways and elucidate the nucleation mechanism of cell fate decisions, with transition states as nucleation sites and pioneer genes as nucleation seeds. The concept of loop flux quantifies the contributions of each cycle flux to cell state transitions, facilitating the understanding of cell dynamics and thermodynamic cost, and providing insights into optimizing biological functions. We also infer cell-cell interactions and cell-type-specific gene regulatory networks, encompassing feedback mechanisms and interaction intensities, predicting genetic perturbation effects on cell fate decisions from single-cell omics data. Essentially, our methodology validates the landscape and flux theory, along with its associated quantifications, offering a framework for exploring the physical principles underlying cellular differentiation and reprogramming and broader biological processes through high-throughput single-cell sequencing experiments.

bHLH transcription factors cooperate with chromatin remodelers to regulate cell fate decisions during Arabidopsis stomatal development

The development of multicellular organisms requires coordinated changes in gene expression that are often mediated by the interaction between transcription factors (TFs) and their corresponding cis-regulatory elements (CREs). During development and differentiation, the accessibility of CREs is dynamically modulated by the epigenome. How the epigenome, CREs, and TFs together exert control over cell fate commitment remains to be fully understood. In the Arabidopsis leaf epidermis, meristemoids undergo a series of stereotyped cell divisions, then switch fate to commit to stomatal differentiation. Newly created or reanalyzed scRNA-seq and ChIP-seq data confirm that stomatal development involves distinctive phases of transcriptional regulation and that differentially regulated genes are bound by the stomatal basic helix-loop-helix (bHLH) TFs. Targets of the bHLHs often reside in repressive chromatin before activation. MNase-seq evidence further suggests that the repressive state can be overcome and remodeled upon activation by specific stomatal bHLHs. We propose that chromatin remodeling is mediated through the recruitment of a set of physical interactors that we identified through proximity labeling-the ATPase-dependent chromatin remodeling SWI/SNF complex and the histone acetyltransferase HAC1. The bHLHs and chromatin remodelers localize to overlapping genomic regions in a hierarchical order. Furthermore, plants with stage-specific knockdown of the SWI/SNF components or HAC1 fail to activate specific bHLH targets and display stomatal development defects. Together, these data converge on a model for how stomatal TFs and epigenetic machinery cooperatively regulate transcription and chromatin remodeling during progressive fate specification.

Living in a noisy world-origins of gene expression noise and its impact on cellular decision-making

The expression level of a gene can vary between genetically identical cells under the same environmental condition-a phenomenon referred to as gene expression noise. Several studies have now elucidated a central role of transcription factors in the generation of expression noise. Transcription factors, as the key components of gene regulatory networks, drive many important cellular decisions in response to cellular and environmental signals. Therefore, a very relevant question is how expression noise impacts gene regulation and influences cellular decision-making. In this Review, we summarize the current understanding of the molecular origins of expression noise, highlighting the role of transcription factors in this process, and discuss the ways in which noise can influence cellular decision-making. As advances in single-cell technologies open new avenues for studying expression noise as well as gene regulatory circuits, a better understanding of the influence of noise on cellular decisions will have important implications for many biological processes.

Fokker-Planck diffusion maps of multiple single cell microglial transcriptomes reveals radial differentiation into substates associated with Alzheimer's pathology

The identification of microglia subtypes is important for understanding the role of innate immunity in neurodegenerative diseases. Current methods of unsupervised cell type identification assume a small noise-to-signal ratio of transcriptome measurements that would produce well-separated cell clusters. However, identification of subtypes is obscured by gene expression noise, diminishing the distances in transcriptome space between distinct cell types and blurring boundaries. Here we use Fokker-Planck (FP) diffusion maps to model cellular differentiation as a stochastic process whereby cells settle into local minima, corresponding to cell subtypes, in a potential landscape constructed from transcriptome data using a nearest neighbor graph approach. By applying critical transition fields, we identify individual cells on the verge of transitioning between subtypes, revealing microglial cells in inactivated, homeostatic state before radially transitioning into various specialized subtypes. Specifically, we show that cells from Alzheimer's disease patients are enriched in a microglia subtype associated to antigen presentation and T-cell recruitment.

CPMI: comprehensive neighborhood-based perturbed mutual information for identifying critical states of complex biological processes

BACKGROUND

There exists a critical transition or tipping point during the complex biological process. Such critical transition is usually accompanied by the catastrophic consequences. Therefore, hunting for the tipping point or critical state is of significant importance to prevent or delay the occurrence of catastrophic consequences. However, predicting critical state based on the high-dimensional small sample data is a difficult problem, especially for single-cell expression data.

RESULTS

In this study, we propose the comprehensive neighbourhood-based perturbed mutual information (CPMI) method to detect the critical states of complex biological processes. The CPMI method takes into account the relationship between genes and neighbours, so as to reduce the noise and enhance the robustness. This method is applied to a simulated dataset and six real datasets, including an influenza dataset, two single-cell expression datasets and three bulk datasets. The method can not only successfully detect the tipping points, but also identify their dynamic network biomarkers (DNBs). In addition, the discovery of transcription factors (TFs) which can regulate DNB genes and nondifferential 'dark genes' validates the effectiveness of our method. The numerical simulation verifies that the CPMI method is robust under different noise strengths and is superior to the existing methods on identifying the critical states.

CONCLUSIONS

In conclusion, we propose a robust computational method, i.e., CPMI, which is applicable in both the bulk and single cell datasets. The CPMI method holds great potential in providing the early warning signals for complex biological processes and enabling early disease diagnosis.

Gene regulatory network topology governs resistance and treatment escape in glioma stem-like cells

Poor prognosis and drug resistance in glioblastoma (GBM) can result from cellular heterogeneity and treatment-induced shifts in phenotypic states of tumor cells, including dedifferentiation into glioma stem-like cells (GSCs). This rare tumorigenic cell subpopulation resists temozolomide, undergoes proneural-to-mesenchymal transition (PMT) to evade therapy, and drives recurrence. Through inference of transcriptional regulatory networks (TRNs) of patient-derived GSCs (PD-GSCs) at single-cell resolution, we demonstrate how the topology of transcription factor interaction networks drives distinct trajectories of cell-state transitions in PD-GSCs resistant or susceptible to cytotoxic drug treatment. By experimentally testing predictions based on TRN simulations, we show that drug treatment drives surviving PD-GSCs along a trajectory of intermediate states, exposing vulnerability to potentiated killing by siRNA or a second drug targeting treatment-induced transcriptional programs governing nongenetic cell plasticity. Our findings demonstrate an approach to uncover TRN topology and use it to rationally predict combinatorial treatments that disrupt acquired resistance in GBM.

Combinatorial interpretation of BMP and WNT controls the decision between primitive streak and extraembryonic fates

BMP signaling is essential for mammalian gastrulation, as it initiates a cascade of signals that control self-organized patterning. As development is highly dynamic, it is crucial to understand how time-dependent combinatorial signaling affects cellular differentiation. Here, we show that BMP signaling duration is a crucial control parameter that determines cell fates upon the exit from pluripotency through its interplay with the induced secondary signal WNT. BMP signaling directly converts cells from pluripotent to extraembryonic fates while simultaneously upregulating Wnt signaling, which promotes primitive streak and mesodermal specification. Using live-cell imaging of signaling and cell fate reporters together with a simple mathematical model, we show that this circuit produces a temporal morphogen effect where, once BMP signal duration is above a threshold for differentiation, intermediate and long pulses of BMP signaling produce specification of mesoderm and extraembryonic fates, respectively. Our results provide a systems-level picture of how these signaling pathways control the landscape of early human development.

BASiCS workflow: a step-by-step analysis of expression variability using single cell RNA sequencing data

Cell-to-cell gene expression variability is an inherent feature of complex biological systems, such as immunity and development. Single-cell RNA sequencing is a powerful tool to quantify this heterogeneity, but it is prone to strong technical noise. In this article, we describe a step-by-step computational workflow that uses the BASiCS Bioconductor package to robustly quantify expression variability within and between known groups of cells (such as experimental conditions or cell types). BASiCS uses an integrated framework for data normalisation, technical noise quantification and downstream analyses, propagating statistical uncertainty across these steps. Within a single seemingly homogeneous cell population, BASiCS can identify highly variable genes that exhibit strong heterogeneity as well as lowly variable genes with stable expression. BASiCS also uses a probabilistic decision rule to identify changes in expression variability between cell populations, whilst avoiding confounding effects related to differences in technical noise or in overall abundance. Using a publicly available dataset, we guide users through a complete pipeline that includes preliminary steps for quality control, as well as data exploration using the scater and scran Bioconductor packages. The workflow is accompanied by a Docker image that ensures the reproducibility of our results.

A logic-incorporated gene regulatory network deciphers principles in cell fate decisions

Organisms utilize gene regulatory networks (GRN) to make fate decisions, but the regulatory mechanisms of transcription factors (TF) in GRNs are exceedingly intricate. A longstanding question in this field is how these tangled interactions synergistically contribute to decision-making procedures. To comprehensively understand the role of regulatory logic in cell fate decisions, we constructed a logic-incorporated GRN model and examined its behavior under two distinct driving forces (noise-driven and signal-driven). Under the noise-driven mode, we distilled the relationship among fate bias, regulatory logic, and noise profile. Under the signal-driven mode, we bridged regulatory logic and progression-accuracy trade-off, and uncovered distinctive trajectories of reprogramming influenced by logic motifs. In differentiation, we characterized a special logic-dependent priming stage by the solution landscape. Finally, we applied our findings to decipher three biological instances: hematopoiesis, embryogenesis, and trans-differentiation. Orthogonal to the classical analysis of expression profile, we harnessed noise patterns to construct the GRN corresponding to fate transition. Our work presents a generalizable framework for top-down fate-decision studies and a practical approach to the taxonomy of cell fate decisions.

A senescence restriction point acting on chromatin integrates oncogenic signals

We identify a senescence restriction point (SeRP) as a critical event for cells to commit to senescence. The SeRP integrates the intensity and duration of oncogenic stress, keeps a memory of previous stresses, and combines oncogenic signals acting on different pathways by modulating chromatin accessibility. Chromatin regions opened upon commitment to senescence are enriched in nucleolar-associated domains, which are gene-poor regions enriched in repeated sequences. Once committed to senescence, cells no longer depend on the initial stress signal and exhibit a characteristic transcriptome regulated by a transcription factor network that includes ETV4, RUNX1, OCT1, and MAFB. Consistent with a tumor suppressor role for this network, the levels of ETV4 and RUNX1 are very high in benign lesions of the pancreas but decrease dramatically in pancreatic ductal adenocarcinomas. The discovery of senescence commitment and its chromatin-linked regulation suggests potential strategies for reinstating tumor suppression in human cancers.

Single-cell omics identifies inflammatory signaling as a trans-differentiation trigger in mouse embryos

Trans-differentiation represents a direct lineage conversion; however, insufficient characterization of this process hinders its potential applications. Here, to explore a potential universal principal for trans-differentiation, we performed single-cell transcriptomic analysis of endothelial-to-hematopoietic transition (EHT), endothelial-to-mesenchymal transition, and epithelial-to-mesenchymal transition in mouse embryos. We applied three scoring indexes of entropies, cell-type signature transcription factor expression, and critical transition signals to show common features underpinning the fate plasticity of transition states. Cross-model comparison identified inflammatory-featured transition states and a common trigger role of interleukin-33 in promoting fate conversions. Multimodal profiling (integrative transcriptomic and chromatin accessibility analysis) demonstrated the inflammatory regulation of hematopoietic specification. Furthermore, multimodal omics and fate-mapping analyses showed that endothelium-specific Spi1, as an inflammatory effector, governs appropriate chromatin accessibility and transcriptional programs to safeguard EHT. Overall, our study employs single-cell omics to identify critical transition states/signals and the common trigger role of inflammatory signaling in developmental-stress-induced fate conversions.

Dysfunction of CCT3-associated network signals for the critical state during progression of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the most common malignant tumors and is a serious threat to human health; thus, early diagnosis and adequate treatment are essential. However, there are still great challenges in identifying the tipping point and detecting early warning signals of early HCC. In this study, we aimed to identify the tipping point (critical state) of and key molecules involved in hepatocarcinogenesis based on time series transcriptome expression data of HCC patients. The phase from veHCC (very early HCC) to eHCC (early HCC) was identified as the critical state in HCC progression, with 143 genes identified as key candidate molecules by combining the DDRTree (dimensionality reduction via graph structure learning) and DNB (dynamic network biomarker) methods. Then, we ranked the candidate genes to verify their mRNA levels using the diethylnitrosamine (DEN)-induced HCC mouse model and identified five early warning signals, namely, CCT3, DSTYK, EIF3E, IARS2 and TXNRD1; these signals can be regarded as the potential early warning signals for the critical state of HCC. We identified CCT3 as an independent prognostic factor for HCC, and functions of CCT3 involving in the "MYCtargets_V1" and "E2F-Targets" are closely related to the progression of HCC. The predictive method combining the DDRTree and DNB methods can not only identify the key critical state before cancer but also determine candidate molecules of critical state, thus providing new insight into the early diagnosis and preemptive treatment of HCC.

A neural network-based model framework for cell-fate decisions and development

Gene regulatory networks (GRNs) fulfill the essential function of maintaining the stability of cellular differentiation states by sustaining lineage-specific gene expression, while driving the progression of development. However, accounting for the relative stability of intermediate differentiation stages and their divergent trajectories remains a major challenge for models of developmental biology. Here, we develop an empirical data-based associative GRN model (AGRN) in which regulatory networks store multilineage stage-specific gene expression profiles as associative memory patterns. These networks are capable of responding to multiple instructive signals and, depending on signal timing and identity, can dynamically drive the differentiation of multipotent cells toward different cell state attractors. The AGRN dynamics can thus generate diverse lineage-committed cell populations in a robust yet flexible manner, providing an attractor-based explanation for signal-driven cell fate decisions during differentiation and offering a readily generalizable modelling tool that can be applied to a wide variety of cell specification systems.

Differentiation is accompanied by a progressive loss in transcriptional memory

BACKGROUND

Cell differentiation requires the integration of two opposite processes, a stabilizing cellular memory, especially at the transcriptional scale, and a burst of gene expression variability which follows the differentiation induction. Therefore, the actual capacity of a cell to undergo phenotypic change during a differentiation process relies upon a modification in this balance which favors change-inducing gene expression variability. However, there are no experimental data providing insight on how fast the transcriptomes of identical cells would diverge on the scale of the very first two cell divisions during the differentiation process.

RESULTS

In order to quantitatively address this question, we developed different experimental methods to recover the transcriptomes of related cells, after one and two divisions, while preserving the information about their lineage at the scale of a single cell division. We analyzed the transcriptomes of related cells from two differentiation biological systems (human CD34+ cells and T2EC chicken primary erythrocytic progenitors) using two different single-cell transcriptomics technologies (scRT-qPCR and scRNA-seq).

CONCLUSIONS

We identified that the gene transcription profiles of differentiating sister cells are more similar to each other than to those of non-related cells of the same type, sharing the same environment and undergoing similar biological processes. More importantly, we observed greater discrepancies between differentiating sister cells than between self-renewing sister cells. Furthermore, a progressive increase in this divergence from first generation to second generation was observed when comparing differentiating cousin cells to self renewing cousin cells. Our results are in favor of a gradual erasure of transcriptional memory during the differentiation process.

Retrospective identification of cell-intrinsic factors that mark pluripotency potential in rare somatic cells

Pluripotency can be induced in somatic cells by the expression of OCT4, KLF4, SOX2, and MYC. Usually only a rare subset of cells reprogram, and the molecular characteristics of this subset remain unknown. We apply retrospective clone tracing to identify and characterize the rare human fibroblasts primed for reprogramming. These fibroblasts showed markers of increased cell cycle speed and decreased fibroblast activation. Knockdown of a fibroblast activation factor identified by our analysis increased the reprogramming efficiency. We provide evidence for a unified model in which cells can move into and out of the primed state over time, explaining how reprogramming appears deterministic at short timescales and stochastic at long timescales. Furthermore, inhibiting the activity of LSD1 enlarged the pool of cells that were primed for reprogramming. Thus, even homogeneous cell populations can exhibit heritable molecular variability that can dictate whether individual rare cells will reprogram or not.

Dynamics of cell-type transition mediated by epigenetic modifications

Maintaining tissue homeostasis requires appropriate regulation of stem cell differentiation. The Waddington landscape posits that gene circuits in a cell form a potential landscape of different cell types, wherein cells follow attractors of the probability landscape to develop into distinct cell types. However, how adult stem cells achieve a delicate balance between self-renewal and differentiation remains unclear. We propose that random inheritance of epigenetic states plays a pivotal role in stem cell differentiation and present a hybrid model of stem cell differentiation induced by epigenetic modifications. Our comprehensive model integrates gene regulation networks, epigenetic state inheritance, and cell regeneration, encompassing multi-scale dynamics ranging from transcription regulation to cell population. Through model simulations, we demonstrate that random inheritance of epigenetic states during cell divisions can spontaneously induce cell differentiation, dedifferentiation, and transdifferentiation. Furthermore, we investigate the influences of interfering with epigenetic modifications and introducing additional transcription factors on the probabilities of dedifferentiation and transdifferentiation, revealing the underlying mechanism of cell reprogramming. This in silico model provides valuable insights into the intricate mechanism governing stem cell differentiation and cell reprogramming and offers a promising path to enhance the field of regenerative medicine.

Understanding Cell Differentiation Through Single-Cell Approaches: Conceptual Challenges of the Systemic Approach

The cells of a multicellular organism are derived from a single zygote and genetically almost 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.

From Cell States to Cell Fates: Control of Cell State Transitions

We examine the coordinated behavior of thousands of genes in cell fate transitions through genome expression as an integrated dynamical system using the concepts of self-organized criticality and coherent stochastic behavior. To quantify the effects of the collective behavior of genes, we adopted the flux balance approach and developed it in a new tool termed expression flux analysis (EFA). Here we describe this tool and demonstrate how its application to specific experimental genome-wide expression data provides new insights into the dynamics of the cell-fate transitions. Particularly, we show that in cell fate change, specific stochastic perturbations can spread over the entire system to guide distinct cell fate transitions through switching cyclic flux flow in the genome engine. Utilization of EFA enables us to elucidate a unified genomic mechanism for when and how cell-fate change occurs through critical transitions.

Lhx3/4 initiates a cardiopharyngeal-specific transcriptional program in response to widespread FGF signaling

Individual signaling pathways, such as fibroblast growth factors (FGFs), can regulate a plethora of inductive events. According to current paradigms, signal-dependent transcription factors (TFs), such as FGF/MapK-activated Ets family factors, partner with lineage-determining factors to achieve regulatory specificity. However, many aspects of this model have not been rigorously investigated. One key question relates to whether lineage-determining factors dictate lineage-specific responses to inductive signals or facilitate these responses in collaboration with other inputs. We utilize the chordate model Ciona robusta to investigate mechanisms generating lineage-specific induction. Previous studies in C. robusta have shown that cardiopharyngeal progenitor cells are specified through the combined activity of FGF-activated Ets1/2.b and an inferred ATTA-binding transcriptional cofactor. Here, we show that the homeobox TF Lhx3/4 serves as the lineage-determining TF that dictates cardiopharyngeal-specific transcription in response to pleiotropic FGF signaling. Targeted knockdown of Lhx3/4 leads to loss of cardiopharyngeal gene expression. Strikingly, ectopic expression of Lhx3/4 in a neuroectodermal lineage subject to FGF-dependent specification leads to ectopic cardiopharyngeal gene expression in this lineage. Furthermore, ectopic Lhx3/4 expression disrupts neural plate morphogenesis, generating aberrant cell behaviors associated with execution of incompatible morphogenetic programs. Based on these findings, we propose that combinatorial regulation by signal-dependent and lineage-determinant factors represents a generalizable, previously uncategorized regulatory subcircuit we term "cofactor-dependent induction." Integration of this subcircuit into theoretical models will facilitate accurate predictions regarding the impact of gene regulatory network rewiring on evolutionary diversification and disease ontogeny.

Reversing pathological cell states: the road less travelled can extend the therapeutic horizon

Acquisition of omics data advances at a formidable pace. Yet, our ability to utilize these data to control cell phenotypes and design interventions that reverse pathological states lags behind. Here, we posit that cell states are determined by core networks that control cell-wide networks. To steer cell fate decisions, core networks connecting genotype to phenotype must be reconstructed and understood. A recent method, cell state transition assessment and regulation (cSTAR), applies perturbation biology to quantify causal connections and mechanistically models how core networks influence cell phenotypes. cSTAR models are akin to digital cell twins enabling us to purposefully convert pathological states back to physiologically normal states. While this capability has a range of applications, here we discuss reverting oncogenic transformation.

A two-galectin network establishes mesenchymal condensation phenotype in limb development

Previous work showed that Gal-1A and Gal-8, two proteins belonging to the galactoside-binding galectin family, are the earliest determinants of the patterning of the skeletal elements of embryonic chicken limbs, and further, that their experimentally determined interactions in the embryonic limb bud can be interpreted via a reaction-diffusion-adhesion (2GL: two galectin plus ligands) model. Here, we use an ordinary differential equation-based approach to analyze the intrinsic switching modality of the 2GL network and characterize the network behavior independent of the diffusive and adhesive arms of the patterning mechanism. We identify two states: where the concentrations of both the galectins are respectively, negligible, and very high. This bistable switch-like system arises via a saddle-node bifurcation from a monostable state. For the case of mass-action production terms, we provide an explicit Lyapunov function for the system, which shows that it has no periodic solutions. Our model therefore predicts that the galectin network may exist in low expression and high expression states separated in space or time, without any intermediate states. We test these predictions in experiments performed with high density cultures of chick limb mesenchymal cells and observe that cells inside precartilage protocondensations express Gal-1A at a much higher rate than those outside, for which it was negligible. The Gal-1A and -8-based patterning network is therefore sufficient to partition the mesenchymal cell population into two discrete cell states with different developmental (chondrogenic vs. non-chondrogenic) fates. When incorporated into an adhesion and diffusion-enabled framework this system can generate a spatially patterned limb skeleton.

SGAE: single-cell gene association entropy for revealing critical states of cell transitions during embryonic development

The critical point or pivotal threshold of cell transition occurs in early embryonic development when cell differentiation culminates in its transition to specific cell fates, at which the cell population undergoes an abrupt and qualitative shift. Revealing such critical points of cell transitions can track cellular heterogeneity and shed light on the molecular mechanisms of cell differentiation. However, precise detection of critical state transitions proves challenging when relying on single-cell RNA sequencing data due to their inherent sparsity, noise, and heterogeneity. In this study, diverging from conventional methods like differential gene analysis or static techniques that emphasize classification of cell types, an innovative computational approach, single-cell gene association entropy (SGAE), is designed for the analysis of single-cell RNA-seq data and utilizes gene association information to reveal critical states of cell transitions. More specifically, through the translation of gene expression data into local SGAE scores, the proposed SGAE can serve as an index to quantitatively assess the resilience and critical properties of genetic regulatory networks, consequently detecting the signal of cell transitions. Analyses of five single-cell datasets for embryonic development demonstrate that the SGAE method achieves better performance in facilitating the characterization of a critical phase transition compared with other existing methods. Moreover, the SGAE value can effectively discriminate cellular heterogeneity over time and performs well in the temporal clustering of cells. Besides, biological functional analysis also indicates the effectiveness of the proposed approach.

Differential variability analysis of single-cell gene expression data

The advent of single-cell RNA sequencing (scRNA-seq) technologies has enabled gene expression profiling at the single-cell resolution, thereby enabling the quantification and comparison of transcriptional variability among individual cells. Although alterations in transcriptional variability have been observed in various biological states, statistical methods for quantifying and testing differential variability between groups of cells are still lacking. To identify the best practices in differential variability analysis of single-cell gene expression data, we propose and compare 12 statistical pipelines using different combinations of methods for normalization, feature selection, dimensionality reduction and variability calculation. Using high-quality synthetic scRNA-seq datasets, we benchmarked the proposed pipelines and found that the most powerful and accurate pipeline performs simple library size normalization, retains all genes in analysis and uses denSNE-based distances to cluster medoids as the variability measure. By applying this pipeline to scRNA-seq datasets of COVID-19 and autism patients, we have identified cellular variability changes between patients with different severity status or between patients and healthy controls.

Identifying a developmental transition in honey bees using gene expression data

In many organisms, interactions among genes lead to multiple functional states, and changes to interactions can lead to transitions into new states. These transitions can be related to bifurcations (or critical points) in dynamical systems theory. Characterizing these collective transitions is a major challenge for systems biology. Here, we develop a statistical method for identifying bistability near a continuous transition directly from high-dimensional gene expression data. We apply the method to data from honey bees, where a known developmental transition occurs between bees performing tasks in the nest and leaving the nest to forage. Our method, which makes use of the expected shape of the distribution of gene expression levels near a transition, successfully identifies the emergence of bistability and links it to genes that are known to be involved in the behavioral transition. This proof of concept demonstrates that going beyond correlative analysis to infer the shape of gene expression distributions might be used more generally to identify collective transitions from gene expression data.

Cell Fate Programming by Transcription Factors and Epigenetic Machinery in Stomatal Development

The development of multi-cellular organisms requires coordinated changes in gene expression that are often mediated by the interaction between transcription factors (TFs) and their corresponding cis-regulatory elements (CREs). During development and differentiation, the accessibility of CREs is dynamically modulated by the epigenome. How the epigenome, CREs and TFs together exert control over cell fate commitment remains to be fully understood. In the Arabidopsis leaf epidermis, meristemoids undergo a series of stereotyped cell divisions, then switch fate to commit to stomatal differentiation. Newly created or reanalyzed scRNA-seq and ChIP-seq data confirm that stomatal development involves distinctive phases of transcriptional regulation and that differentially regulated genes are bound by the stomatal basic-helix-loop-helix (bHLH) TFs. Targets of the bHLHs often reside in repressive chromatin before activation. MNase-seq evidence further suggests that the repressive state can be overcome and remodeled upon activation by specific stomatal bHLHs. We propose that chromatin remodeling is mediated through the recruitment of a set of physical interactors that we identified through proximity labeling - the ATPase-dependent chromatin remodeling SWI/SNF complex and the histone acetyltransferase HAC1. The bHLHs and chromatin remodelers localize to overlapping genomic regions in a hierarchical order. Furthermore, plants with stage-specific knock-down of the SWI/SNF components or HAC1 fail to activate specific bHLH targets and display stomatal development defects. Together these data converge on a model for how stomatal TFs and epigenetic machinery cooperatively regulate transcription and chromatin remodeling during progressive fate specification.

An image-guided microfluidic system for single-cell lineage tracking

Cell lineage tracking is a long-standing and unresolved problem in biology. Microfluidic technologies have the potential to address this problem, by virtue of their ability to manipulate and process single-cells in a rapid, controllable and efficient manner. Indeed, when coupled with traditional imaging approaches, microfluidic systems allow the experimentalist to follow single-cell divisions over time. Herein, we present a valve-based microfluidic system able to probe the decision-making processes of single-cells, by tracking their lineage over multiple generations. The system operates by trapping single-cells within growth chambers, allowing the trapped cells to grow and divide, isolating sister cells after a user-defined number of divisions and finally extracting them for downstream transcriptome analysis. The platform incorporates multiple cell manipulation operations, image processing-based automation for cell loading and growth monitoring, reagent addition and device washing. To demonstrate the efficacy of the microfluidic workflow, 6C2 (chicken erythroleukemia) and T2EC (primary chicken erythrocytic progenitors) cells are tracked inside the microfluidic device over two generations, with a cell viability rate in excess of 90%. Sister cells are successfully isolated after division and extracted within a 500 nL volume, which was demonstrated to be compatible with downstream single-cell RNA sequencing analysis.

Systematic analysis and optimization of early warning signals for critical transitions using distribution data

Abrupt shifts between alternative regimes occur in complex systems, from cell regulation to brain functions to ecosystems. Several model-free early warning signals (EWS) have been proposed to detect impending transitions, but failure or poor performance in some systems have called for better investigation of their generic applicability. Notably, there are still ongoing debates whether such signals can be successfully extracted from data in particular from biological experiments. In this work, we systematically investigate properties and performance of dynamical EWS in different deteriorating conditions, and we propose an optimized combination to trigger warnings as early as possible, eventually verified on experimental data from microbiological populations. Our results explain discrepancies observed in the literature between warning signs extracted from simulated models and from real data, provide guidance for EWS selection based on desired systems and suggest an optimized composite indicator to alert for impending critical transitions using distribution data.

Single-cell transcriptional uncertainty landscape of cell differentiation

Background: Single-cell studies have demonstrated the presence of significant cell-to-cell heterogeneity in gene expression. Whether such heterogeneity is only a bystander or has a functional role in the cell differentiation process is still hotly debated. Methods: In this study, we quantified and followed single-cell transcriptional uncertainty - a measure of gene transcriptional stochasticity in single cells - in 10 cell differentiation systems of varying cell lineage progressions, from single to multi-branching trajectories, using the stochastic two-state gene transcription model. Results: By visualizing the transcriptional uncertainty as a landscape over a two-dimensional representation of the single-cell gene expression data, we observed universal features in the cell differentiation trajectories that include: (i) a peak in single-cell uncertainty during transition states, and in systems with bifurcating differentiation trajectories, each branching point represents a state of high transcriptional uncertainty; (ii) a positive correlation of transcriptional uncertainty with transcriptional burst size and frequency; (iii) an increase in RNA velocity preceding the increase in the cell transcriptional uncertainty. Conclusions: Our findings suggest a possible universal mechanism during the cell differentiation process, in which stem cells engage stochastic exploratory dynamics of gene expression at the start of the cell differentiation by increasing gene transcriptional bursts, and disengage such dynamics once cells have decided on a particular terminal cell identity. Notably, the peak of single-cell transcriptional uncertainty signifies the decision-making point in the cell differentiation process.

exFINDER: identify external communication signals using single-cell transcriptomics data

Cells make decisions through their communication with other cells and receiving signals from their environment. Using single-cell transcriptomics, computational tools have been developed to infer cell-cell communication through ligands and receptors. However, the existing methods only deal with signals sent by the measured cells in the data, the received signals from the external system are missing in the inference. Here, we present exFINDER, a method that identifies such external signals received by the cells in the single-cell transcriptomics datasets by utilizing the prior knowledge of signaling pathways. In particular, exFINDER can uncover external signals that activate the given target genes, infer the external signal-target signaling network (exSigNet), and perform quantitative analysis on exSigNets. The applications of exFINDER to scRNA-seq datasets from different species demonstrate the accuracy and robustness of identifying external signals, revealing critical transition-related signaling activities, inferring critical external signals and targets, clustering signal-target paths, and evaluating relevant biological events. Overall, exFINDER can be applied to scRNA-seq data to reveal the external signal-associated activities and maybe novel cells that send such signals.

Macroscopic quorum sensing sustains differentiating embryonic stem cells

Cells can secrete molecules that help each other's replication. In cell cultures, chemical signals might diffuse only within a cell colony or between colonies. A chemical signal's interaction length-how far apart interacting cells are-is often assumed to be some value without rigorous justifications because molecules' invisible paths and complex multicellular geometries pose challenges. Here we present an approach, combining mathematical models and experiments, for determining a chemical signal's interaction length. With murine embryonic stem (ES) cells as a testbed, we found that differentiating ES cells secrete FGF4, among others, to communicate over many millimeters in cell culture dishes and, thereby, form a spatially extended, macroscopic entity that grows only if its centimeter-scale population density is above a threshold value. With this 'macroscopic quorum sensing', an isolated macroscopic, but not isolated microscopic, colony can survive differentiation. Our integrated approach can determine chemical signals' interaction lengths in generic multicellular communities.

Critical transition and reversion of tumorigenesis

Cancer is caused by the accumulation of genetic alterations and therefore has been historically considered to be irreversible. Intriguingly, several studies have reported that cancer cells can be reversed to be normal cells under certain circumstances. Despite these experimental observations, conceptual and theoretical frameworks that explain these phenomena and enable their exploration in a systematic way are lacking. In this review, we provide an overview of cancer reversion studies and describe recent advancements in systems biological approaches based on attractor landscape analysis. We suggest that the critical transition in tumorigenesis is an important clue for achieving cancer reversion. During tumorigenesis, a critical transition may occur at a tipping point, where cells undergo abrupt changes and reach a new equilibrium state that is determined by complex intracellular regulatory events. We introduce a conceptual framework based on attractor landscapes through which we can investigate the critical transition in tumorigenesis and induce its reversion by combining intracellular molecular perturbation and extracellular signaling controls. Finally, we present a cancer reversion therapy approach that may be a paradigm-changing alternative to current cancer cell-killing therapies.

Cell Features Reconstruction from Gene Association Network of Single Cell

Gene expression as an unstable form of cell characterization has been widely used for single-cell analyses. Although there are cell-specific networks (CSN) to explore stable gene associations within a single cell, the amount of information in CSN is huge and there is no method to measure the interaction level between genes. Therefore, this paper presents a two-level approach to reconstructing single-cell features, which transforms the original gene expression feature into the gene ontology feature and gene interaction feature. Specifically, we first squeeze all CSNs into a cell network feature matrix (CNFM) by fusing the global position and neighborhood influence of genes. Next, we propose a computational method of gene gravitation based on CNFM to quantify the extent of gene-gene interaction, and we can construct a gene gravitation network for single cells. Finally, we further design a novel index of gene gravitation entropy to quantitatively evaluate the level of single-cell differentiation. The experiments on eight different scRNA-seq datasets show the effectiveness and broad application prospects of our method.

exFINDER: identify external communication signals using single-cell transcriptomics data

Cells make decisions through their communication with other cells and receiving signals from their environment. Using single-cell transcriptomics, computational tools have been developed to infer cell-cell communication through ligands and receptors. However, the existing methods only deal with signals sent by the measured cells in the data, the received signals from the external system are missing in the inference. Here, we present exFINDER, a method that identifies such external signals received by the cells in the single-cell transcriptomics datasets by utilizing the prior knowledge of signaling pathways. In particular, exFINDER can uncover external signals that activate the given target genes, infer the external signal-target signaling network (exSigNet), and perform quantitative analysis on exSigNets. The applications of exFINDER to scRNA-seq datasets from different species demonstrate the accuracy and robustness of identifying external signals, revealing critical transition-related signaling activities, inferring critical external signals and targets, clustering signal-target paths, and evaluating relevant biological events. Overall, exFINDER can be applied to scRNA-seq data to reveal the external signal-associated activities and maybe novel cells that send such signals.

SPNE: sample-perturbed network entropy for revealing critical states of complex biological systems

Complex biological systems do not always develop smoothly but occasionally undergo a sharp transition; i.e. there exists a critical transition or tipping point at which a drastic qualitative shift occurs. Hunting for such a critical transition is important to prevent or delay the occurrence of catastrophic consequences, such as disease deterioration. However, the identification of the critical state for complex biological systems is still a challenging problem when using high-dimensional small sample data, especially where only a certain sample is available, which often leads to the failure of most traditional statistical approaches. In this study, a novel quantitative method, sample-perturbed network entropy (SPNE), is developed based on the sample-perturbed directed network to reveal the critical state of complex biological systems at the single-sample level. Specifically, the SPNE approach effectively quantifies the perturbation effect caused by a specific sample on the directed network in terms of network entropy and thus captures the criticality of biological systems. This model-free method was applied to both bulk and single-cell expression data. Our approach was validated by successfully detecting the early warning signals of the critical states for six real datasets, including four tumor datasets from The Cancer Genome Atlas (TCGA) and two single-cell datasets of cell differentiation. In addition, the functional analyses of signaling biomarkers demonstrated the effectiveness of the analytical and computational results.

A novel network approach to multiscale biological regulation

Modeling systems at multiple interacting scales is probably the most relevant task for pursuing a physically motivated explanation of biological regulation. In a new study, Smart and Zilman develop a convincing, albeit preliminary, model of the interplay between the cell microscale and the macroscopic tissue organization in biological systems.

Retrospective identification of intrinsic factors that mark pluripotency potential in rare somatic cells

Pluripotency can be induced in somatic cells by the expression of the four "Yamanaka" factors OCT4, KLF4, SOX2, and MYC. However, even in homogeneous conditions, usually only a rare subset of cells admit reprogramming, and the molecular characteristics of this subset remain unknown. Here, we apply retrospective clone tracing to identify and characterize the individual human fibroblast cells that are primed for reprogramming. These fibroblasts showed markers of increased cell cycle speed and decreased fibroblast activation. Knockdown of a fibroblast activation factor identified by our analysis led to increased reprogramming efficiency, identifying it as a barrier to reprogramming. Changing the frequency of reprogramming by inhibiting the activity of LSD1 led to an enlarging of the pool of cells that were primed for reprogramming. Our results show that even homogeneous cell populations can exhibit heritable molecular variability that can dictate whether individual rare cells will reprogram or not.

Distance covariance entropy reveals primed states and bifurcation dynamics in single-cell RNA-Seq data

Cell-fate transitions are fundamental to development and differentiation. Studying them with single-cell omic data is important to advance our understanding of the cell-fate commitment process, yet this remains challenging. Here we present a computational method called DICE, which analyzes the entropy of expression covariation patterns and which is applicable to static and dynamically changing cell populations. Using only single-cell RNA-Seq data, DICE is able to predict multipotent primed states and their regulatory factors, which we subsequently validate with single-cell epigenomic data. DICE reveals that primed states are often defined by epigenetic regulators or pioneer factors alongside lineage-specific transcription factors. In developmental time course single-cell RNA-Seq datasets, DICE can pinpoint the timing of bifurcations more precisely than lineage-trajectory inference algorithms or competing variance-based methods. In summary, by studying the dynamic changes of expression covariation entropy, DICE can help elucidate primed states and bifurcation dynamics without the need for single-cell epigenomic data.