Transition states and cell fate decisions in epigenetic landscapes

scMaui: a widely applicable deep learning framework for single-cell multiomics integration in the presence of batch effects and missing data

The recent advances in high-throughput single-cell sequencing have created an urgent demand for computational models which can address the high complexity of single-cell multiomics data. Meticulous single-cell multiomics integration models are required to avoid biases towards a specific modality and overcome sparsity. Batch effects obfuscating biological signals must also be taken into account. Here, we introduce a new single-cell multiomics integration model, Single-cell Multiomics Autoencoder Integration (scMaui) based on variational product-of-experts autoencoders and adversarial learning. scMaui calculates a joint representation of multiple marginal distributions based on a product-of-experts approach which is especially effective for missing values in the modalities. Furthermore, it overcomes limitations seen in previous VAE-based integration methods with regard to batch effect correction and restricted applicable assays. It handles multiple batch effects independently accepting both discrete and continuous values, as well as provides varied reconstruction loss functions to cover all possible assays and preprocessing pipelines. We demonstrate that scMaui achieves superior performance in many tasks compared to other methods. Further downstream analyses also demonstrate its potential in identifying relations between assays and discovering hidden subpopulations.

Cell cycle-driven transcriptome maturation confers multilineage competence to cardiopharyngeal progenitors

During development, stem and progenitor cells divide and transition through germ layer- and lineage-specific multipotent states to generate the diverse cell types that compose an animal. Defined changes in biomolecular composition underlie the progressive loss of potency and acquisition of lineage-specific characteristics. For example, multipotent cardiopharyngeal progenitors display multilineage transcriptional priming, whereby both the cardiac and pharyngeal muscle programs are partially active and coexist in the same progenitor cells, while their daughter cells engage in a cardiac or pharyngeal muscle differentiation path only after cell division. Here, using the tunicate Ciona, we studied the acquisition of multilineage competence and the coupling between fate decisions and cell cycle progression. We showed that multipotent cardiopharyngeal progenitors acquire the competence to produce distinct Tbx1/10 (+) and (-) daughter cells shortly before mitosis, which is necessary for Tbx1/10 activation. By combining transgene-based sample barcoding with single cell RNA-seq (scRNA-seq), we uncovered transcriptome-wide dynamics in migrating cardiopharyngeal progenitors as cells progress through G1, S and G2 phases. We termed this process "transcriptome maturation", and identified candidate "mature genes", including the Rho GAP-coding gene Depdc1 , which peak in late G2. Functional assays indicated that transcriptome maturation fosters cardiopharyngeal competence, in part through multilineage priming and proper oriented and asymmetric division that influences subsequent fate decisions, illustrating the concept of "behavioral competence". Both classic feedforward circuits and coupling with cell cycle progression drive transcriptome maturation, uncovering distinct levels of coupling between cell cycle progression and fateful molecular transitions. We propose that coupling competence and fate decision with the G2 and G1 phases, respectively, ensures the timely deployment of lineage-specific programs.

The Molecular Basis of Differentiation Wave Activity in Embryogenesis

As development varies greatly across the tree of life, it may seem difficult to suggest a model that proposes a single mechanism for understanding collective cell behaviors and the coordination of tissue formation. Here we propose a mechanism called differentiation waves, which unify many disparate results involving developmental systems from across the tree of life. We demonstrate how a relatively simple model of differentiation proceeds not from function-related molecular mechanisms, but from so-called differentiation waves. A phenotypic model of differentiation waves is introduced, and its relation to molecular mechanisms is proposed. These waves contribute to a differentiation tree, which is an alternate way of viewing cell lineage and local action of the molecular factors. We construct a model of differentiation wave-related molecular mechanisms (genome, epigenome, and proteome) based on bioinformatic data from the nematode Caenorhabditis elegans. To validate this approach across different modes of development, we evaluate protein expression across different types of development by comparing Caenorhabditis elegans with several model organisms: fruit flies (Drosophila melanogaster), yeast (Saccharomyces cerevisiae), and mouse (Mus musculus). Inspired by gene regulatory networks, two Models of Interactive Contributions (fully-connected MICs and ordered MICs) are used to suggest potential genomic contributions to differentiation wave-related proteins. This, in turn, provides a framework for understanding differentiation and development.

Approximate Bayesian computation for inferring Waddington landscapes from single-cell data

Single-cell technologies allow us to gain insights into cellular processes at unprecedented resolution. In stem cell and developmental biology snapshot data allow us to characterize how the transcriptional states of cells change between successive cell types. Here, we show how approximate Bayesian computation (ABC) can be employed to calibrate mathematical models against single-cell data. In our simulation study, we demonstrate the pivotal role of the adequate choice of distance measures appropriate for single-cell data. We show that for good distance measures, notably optimal transport with the Sinkhorn divergence, we can infer parameters for mathematical models from simulated single-cell data. We show that the ABC posteriors can be used (i) to characterize parameter sensitivity and identify dependencies between different parameters and (ii) to construct representations of the Waddington or epigenetic landscape, which forms a popular and interpretable representation of the developmental dynamics. In summary, these results pave the way for fitting mechanistic models of stem cell differentiation to single-cell data.

Cell State Transition Models Stratify Breast Cancer Cell Phenotypes and Reveal New Therapeutic Targets

Understanding signaling patterns of transformation and controlling cell phenotypes is a challenge of current biology. Here we applied a cell State Transition Assessment and Regulation (cSTAR) approach to a perturbation dataset of single cell phosphoproteomic patterns of multiple breast cancer (BC) and normal breast tissue-derived cell lines. Following a separation of luminal, basal, and normal cell states, we identified signaling nodes within core control networks, delineated causal connections, and determined the primary drivers underlying oncogenic transformation and transitions across distinct BC subtypes. Whereas cell lines within the same BC subtype have different mutational and expression profiles, the architecture of the core network was similar for all luminal BC cells, and mTOR was a main oncogenic driver. In contrast, core networks of basal BC were heterogeneous and segregated into roughly four major subclasses with distinct oncogenic and BC subtype drivers. Likewise, normal breast tissue cells were separated into two different subclasses. Based on the data and quantified network topologies, we derived mechanistic cSTAR models that serve as digital cell twins and allow the deliberate control of cell movements within a Waddington landscape across different cell states. These cSTAR models suggested strategies of normalizing phosphorylation networks of BC cell lines using small molecule inhibitors.

Alteration in DNA methylation patterns: Epigenetic signatures in gastrointestinal cancers

Abnormalities in epigenetic modifications can cause malignant transformations in cells, leading to cancers of the gastrointestinal (GI) tract, which accounts for 20% of all cancers worldwide. Among the epigenetic alterations, DNA hypomethylation is associated with genomic instability. In addition, CpG methylation and promoter hypermethylation have been recognized as biomarkers for different malignancies. In GI cancers, epigenetic alterations affect genes responsible for cell cycle control, DNA repair, apoptosis, and tumorigenic-specific signaling pathways. Understanding the pattern of alterations in DNA methylation in GI cancers could help scientists discover new molecular-based pharmaceutical treatments. This study highlights alterations in DNA methylation in GI cancers. Understanding epigenetic differences among GI cancers may improve targeted therapies and lead to the discovery of new diagnostic biomarkers.

Establishing a conceptual framework for holistic cell states and state transitions

Cell states were traditionally defined by how they looked, where they were located, and what functions they performed. In this post-genomic era, the field is largely focused on a molecular view of cell state. Moving forward, we anticipate that the observables used to define cell states will evolve again as single-cell imaging and analytics are advancing at a breakneck pace via the collection of large-scale, systematic cell image datasets and the application of quantitative image-based data science methods. This is, therefore, a key moment in the arc of cell biological research to develop approaches that integrate the spatiotemporal observables of the physical structure and organization of the cell with molecular observables toward the concept of a holistic cell state. In this perspective, we propose a conceptual framework for holistic cell states and state transitions that is data-driven, practical, and useful to enable integrative analyses and modeling across many data types.

Decoding the principle of cell-fate determination for its reverse control

Understanding and manipulating cell fate determination is pivotal in biology. Cell fate is determined by intricate and nonlinear interactions among molecules, making mathematical model-based quantitative analysis indispensable for its elucidation. Nevertheless, obtaining the essential dynamic experimental data for model development has been a significant obstacle. However, recent advancements in large-scale omics data technology are providing the necessary foundation for developing such models. Based on accumulated experimental evidence, we can postulate that cell fate is governed by a limited number of core regulatory circuits. Following this concept, we present a conceptual control framework that leverages single-cell RNA-seq data for dynamic molecular regulatory network modeling, aiming to identify and manipulate core regulatory circuits and their master regulators to drive desired cellular state transitions. We illustrate the proposed framework by applying it to the reversion of lung cancer cell states, although it is more broadly applicable to understanding and controlling a wide range of cell-fate determination processes.

What can we learn when fitting a simple telegraph model to a complex gene expression model?

In experiments, the distributions of mRNA or protein numbers in single cells are often fitted to the random telegraph model which includes synthesis and decay of mRNA or protein, and switching of the gene between active and inactive states. While commonly used, this model does not describe how fluctuations are influenced by crucial biological mechanisms such as feedback regulation, non-exponential gene inactivation durations, and multiple gene activation pathways. Here we investigate the dynamical properties of four relatively complex gene expression models by fitting their steady-state mRNA or protein number distributions to the simple telegraph model. We show that despite the underlying complex biological mechanisms, the telegraph model with three effective parameters can accurately capture the steady-state gene product distributions, as well as the conditional distributions in the active gene state, of the complex models. Some effective parameters are reliable and can reflect realistic dynamic behaviors of the complex models, while others may deviate significantly from their real values in the complex models. The effective parameters can also be applied to characterize the capability for a complex model to exhibit multimodality. Using additional information such as single-cell data at multiple time points, we provide an effective method of distinguishing the complex models from the telegraph model. Furthermore, using measurements under varying experimental conditions, we show that fitting the mRNA or protein number distributions to the telegraph model may even reveal the underlying gene regulation mechanisms of the complex models. The effectiveness of these methods is confirmed by analysis of single-cell data for E. coli and mammalian cells. All these results are robust with respect to cooperative transcriptional regulation and extrinsic noise. In particular, we find that faster relaxation speed to the steady state results in more precise parameter inference under large extrinsic noise.

Cellular remodeling and JAK inhibition promote zygotic gene expression in the Ciona germline

Transcription control is a major determinant of cell fate decisions in somatic tissues. By contrast, early germline fate specification in numerous vertebrate and invertebrate species relies extensively on RNA-level regulation, exerted on asymmetrically inherited maternal supplies, with little-to-no zygotic transcription. However delayed, a maternal-to-zygotic transition is nevertheless poised to complete the deployment of pre-gametic programs in the germline. Here, we focus on early germline specification in the tunicate Ciona to study zygotic genome activation. We first demonstrate that a peculiar cellular remodeling event excludes localized postplasmic Pem-1 mRNA, which encodes the general inhibitor of transcription. Subsequently, zygotic transcription begins in Pem-1-negative primordial germ cells (PGCs), as revealed by histochemical detection of elongating RNA Polymerase II, and nascent Mef2 transcripts. In addition, we uncover a provisional antagonism between JAK and MEK/BMPRI/GSK3 signaling, which controls the onset of zygotic gene expression, following cellular remodeling of PGCs. We propose a 2-step model for the onset of zygotic transcription in the Ciona germline and discuss the significance of germ plasm dislocation and remodeling in the context of developmental fate specification.

PhyloVelo enhances transcriptomic velocity field mapping using monotonically expressed genes

Single-cell RNA sequencing (scRNA-seq) is a powerful approach for studying cellular differentiation, but accurately tracking cell fate transitions can be challenging, especially in disease conditions. Here we introduce PhyloVelo, a computational framework that estimates the velocity of transcriptomic dynamics by using monotonically expressed genes (MEGs) or genes with expression patterns that either increase or decrease, but do not cycle, through phylogenetic time. Through integration of scRNA-seq data with lineage information, PhyloVelo identifies MEGs and reconstructs a transcriptomic velocity field. We validate PhyloVelo using simulated data and Caenorhabditis elegans ground truth data, successfully recovering linear, bifurcated and convergent differentiations. Applying PhyloVelo to seven lineage-traced scRNA-seq datasets, generated using CRISPR-Cas9 editing, lentiviral barcoding or immune repertoire profiling, demonstrates its high accuracy and robustness in inferring complex lineage trajectories while outperforming RNA velocity. Additionally, we discovered that MEGs across tissues and organisms share similar functions in translation and ribosome biogenesis.

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.

Contributions of transcriptional noise to leukaemia evolution: KAT2A as a case-study

Transcriptional noise is proposed to participate in cell fate changes, but contributions to mammalian cell differentiation systems, including cancer, remain associative. Cancer evolution is driven by genetic variability, with modulatory or contributory participation of epigenetic variants. Accumulation of epigenetic variants enhances transcriptional noise, which can facilitate cancer cell fate transitions. Acute myeloid leukaemia (AML) is an aggressive cancer with strong epigenetic dependencies, characterized by blocked differentiation. It constitutes an attractive model to probe links between transcriptional noise and malignant cell fate regulation. Gcn5/KAT2A is a classical epigenetic transcriptional noise regulator. Its loss increases transcriptional noise and modifies cell fates in stem and AML cells. By reviewing the analysis of KAT2A-depleted pre-leukaemia and leukaemia models, I discuss that the net result of transcriptional noise is diversification of cell fates secondary to alternative transcriptional programmes. Cellular diversification can enable or hinder AML progression, respectively, by differentiation of cell types responsive to mutations, or by maladaptation of leukaemia stem cells. KAT2A-dependent noise-responsive genes participate in ribosome biogenesis and KAT2A loss destabilizes translational activity. I discuss putative contributions of perturbed translation to AML biology, and propose KAT2A loss as a model for mechanistic integration of transcriptional and translational control of noise and fate decisions. This article is part of a discussion meeting issue 'Causes and consequences of stochastic processes in development and disease'.

Transcriptome free energy can serve as a dynamic patient-specific biomarker in acute myeloid leukemia

Acute myeloid leukemia (AML) is prevalent in both adult and pediatric patients. Despite advances in patient categorization, the heterogeneity of AML remains a challenge. Recent studies have explored the use of gene expression data to enhance AML diagnosis and prognosis, however, alternative approaches rooted in physics and chemistry may provide another level of insight into AML transformation. Utilizing publicly available databases, we analyze 884 human and mouse blood and bone marrow samples. We employ a personalized medicine strategy, combining state-transition theory and surprisal analysis, to assess the RNA transcriptome of individual patients. The transcriptome is transformed into physical parameters that represent each sample's steady state and the free energy change (FEC) from that steady state, which is the state with the lowest free energy.We found the transcriptome steady state was invariant across normal and AML samples. FEC, representing active molecular processes, varied significantly between samples and was used to create patient-specific barcodes to characterize the biology of the disease. We discovered that AML samples that were in a transition state had the highest FEC. This disease state may be characterized as the most unstable and hence the most therapeutically targetable since a change in free energy is a thermodynamic requirement for disease progression. We also found that distinct sets of ongoing processes may be at the root of otherwise similar clinical phenotypes, implying that our integrated analysis of transcriptome profiles may facilitate a personalized medicine approach to cure AML and restore a steady state in each patient.

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.

Nano-CUT&Tag for multimodal chromatin profiling at single-cell resolution

The ability to comprehensively analyze the chromatin state with single-cell resolution is crucial for understanding gene regulatory principles in heterogenous tissues or during development. Recently, we developed a nanobody-based single-cell CUT&Tag (nano-CT) protocol to simultaneously profile three epigenetic modalities-two histone marks and open chromatin state-from the same single cell. Nano-CT implements a new set of secondary nanobody-Tn5 fusion proteins to direct barcoded tagmentation by Tn5 transposase to genomic targets labeled by primary antibodies raised in different species. Such nanobody-Tn5 fusion proteins are currently not commercially available, and their in-house production and purification can be completed in 3-4 d by following our detailed protocol. The single-cell indexing in nano-CT is performed on a commercially available platform, making it widely accessible to the community. In comparison to other multimodal methods, nano-CT stands out in data complexity, low sample requirements and the flexibility to choose two of the three modalities. In addition, nano-CT works efficiently with fresh brain samples, generating multimodal epigenomic profiles for thousands of brain cells at single-cell resolution. The nano-CT protocol can be completed in just 3 d by users with basic skills in standard molecular biology and bioinformatics, although previous experience with single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) is beneficial for more in-depth data analysis. As a multimodal assay, nano-CT holds immense potential to reveal interactions of various chromatin modalities, to explore epigenetic heterogeneity and to increase our understanding of the role and interplay that chromatin dynamics has in cellular development.

Epigenetics in T-cell driven inflammation and cancer

For decades, scientists have been investigating how processes such as gene expression, stem cell plasticity, and cell differentiation can be modulated. The discovery of epigenetics helped unravel these processes and enabled the identification of major underlying mechanisms that, for example, are central for T cell maturation. T cells go through various stages in their development evolving from progenitor cells into double positive CD4/CD8 T cells that finally leave the thymus as naïve T cells. One major mechanism driving T cell maturation is the modulation of gene activity by temporally sequenced transcription of spatially exposed gene loci. DNA methylation, demethylation, and acetylation are key processes that enable a sequenced gene expression required for T cell differentiation. In vivo, differentiated T cells are subjected to enormous pressures originating from the microenvironment. Signals from this environment, particularly from an inflammatory or a tumor microenvironment, can push T cells to differentiate into specific effector and memory T cells, and even prompt T cells to adopt a state of dysfunctional exhaustion, en route of an epigenetically controlled mechanism. Fundamentals of these processes will be discussed in this review highlighting potential therapeutic interventions, in particular those beneficial to revive exhausted T cells.

Using steady-state formula to estimate time-dependent parameters of stochastic gene transcription models

When studying stochastic gene transcription, it is important to understand how system parameters are temporally modulated in response to varying environments. Experimentally, the dynamic distribution data of RNA copy numbers measured at multiple time points are often fitted to stochastic transcription models to estimate time-dependent parameters. However, current methods require determining which parameters are time-dependent, as well as their analytical formulas, before the optimal fit. In this study, we developed a method to estimate time-dependent parameters in a classical two-state model without prior assumptions regarding the system parameters. At each measured time point, the method fitted the dynamic distribution data using a steady-state distribution formula, in which the estimated constant parameters were approximated as time-dependent parameter values at the measured time point. The accuracy of this method can be guaranteed for RNA molecules with relatively high degradation rates and genes with relatively slow responses to induction. We quantify the accuracy of the method and implemented this method on two sets of dynamic distribution data from prokaryotic and eukaryotic cells, and revealed the temporal modulation of transcription burst size in response to environmental changes.

Studying temporal dynamics of single cells: expression, lineage and regulatory networks

Learning how multicellular organs are developed from single cells to different cell types is a fundamental problem in biology. With the high-throughput scRNA-seq technology, computational methods have been developed to reveal the temporal dynamics of single cells from transcriptomic data, from phenomena on cell trajectories to the underlying mechanism that formed the trajectory. There are several distinct families of computational methods including Trajectory Inference (TI), Lineage Tracing (LT), and Gene Regulatory Network (GRN) Inference which are involved in such studies. This review summarizes these computational approaches which use scRNA-seq data to study cell differentiation and cell fate specification as well as the advantages and limitations of different methods. We further discuss how GRNs can potentially affect cell fate decisions and trajectory structures.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-023-01090-5.

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.

Repeated Decision Stumping Distils Simple Rules from Single-Cell Data

Single-cell data afford unprecedented insights into molecular processes. But the complexity and size of these data sets have proved challenging and given rise to a large armory of statistical and machine learning approaches. The majority of approaches focuses on either describing features of these data, or making predictions and classifying unlabeled samples. In this study, we introduce repeated decision stumping (ReDX) as a method to distill simple models from single-cell data. We develop decision trees of depth one-hence "stumps"-to identify in an inductive manner, gene products involved in driving cell fate transitions, and in applications to published data we are able to discover the key players involved in these processes in an unbiased manner without prior knowledge. Our algorithm is deliberately targeting the simplest possible candidate hypotheses that can be extracted from complex high-dimensional data. There are three reasons for this: (1) the predictions become straightforwardly testable hypotheses; (2) the identified candidates form the basis for further mechanistic model development, for example, for engineering and synthetic biology interventions; and (3) this approach complements existing descriptive modeling approaches and frameworks. The approach is computationally efficient, has remarkable predictive power, including in simulation studies where the ground truth is known, and yields robust and statistically stable predictors; the same set of candidates is generated by applying the algorithm to different subsamples of experimental data.

Rethinking cancer initiation: The role of large-scale mutational events

Cancer initiation is revisited in light of recent discoveries in cancer pathogenesis. Of note is the detection of mutated cancer genes in benign conditions. More significantly, somatic clones, which harbor mutations in cancer genes, arise in normal tissues from early development through adulthood, but seldom do they transform into cancer. Further, clustered mutational events-kataegis, chromothripsis and chromoplexy-are widespread in cancer, generating point mutations and chromosomal rearrangements in a single cellular catastrophe. These observations are contrary to the prevailing somatic mutation theory, which states that a cancer is caused by the gradual accumulation of mutations over time. A different perspective is proposed within the framework of Waddington's epigenetic landscape wherein tumorigenesis is viewed primarily as a disruption of cell development. Cell types are defined by their specific gene-expression profiles, determined by the gene regulatory network, and can be regarded as attractor states of the network dynamics: they represent specific, self-stabilizing patterns of gene activities across the genome. However, large-scale mutational events reshape the landscape topology, creating abnormal "unphysiological" attractors. This is the crux of the process of initiation. Initiation primes the cell for conversion into a tumor phenotype by oncogenes and tumor suppressor genes, which drive cell proliferation and clonal diversification. This view of tumorigenesis calls for a different approach to therapy.

Deciphering driver regulators of cell fate decisions from single-cell transcriptomics data with CEFCON

Single-cell technologies enable the dynamic analyses of cell fate mapping. However, capturing the gene regulatory relationships and identifying the driver factors that control cell fate decisions are still challenging. We present CEFCON, a network-based framework that first uses a graph neural network with attention mechanism to infer a cell-lineage-specific gene regulatory network (GRN) from single-cell RNA-sequencing data, and then models cell fate dynamics through network control theory to identify driver regulators and the associated gene modules, revealing their critical biological processes related to cell states. Extensive benchmarking tests consistently demonstrated the superiority of CEFCON in GRN construction, driver regulator identification, and gene module identification over baseline methods. When applied to the mouse hematopoietic stem cell differentiation data, CEFCON successfully identified driver regulators for three developmental lineages, which offered useful insights into their differentiation from a network control perspective. Overall, CEFCON provides a valuable tool for studying the underlying mechanisms of cell fate decisions from single-cell RNA-seq data.

Haematopoietic development and HSC formation in vitro: promise and limitations of gastruloid models

Haematopoietic stem cells (HSCs) are the most extensively studied adult stem cells. Yet, six decades after their first description, reproducible and translatable generation of HSC in vitro remains an unmet challenge. HSC production in vitro is confounded by the multi-stage nature of blood production during development. Specification of HSC is a late event in embryonic blood production and depends on physical and chemical cues which remain incompletely characterised. The precise molecular composition of the HSC themselves is incompletely understood, limiting approaches to track their origin in situ in the appropriate cellular, chemical and mechanical context. Embryonic material at the point of HSC emergence is limiting, highlighting the need for an in vitro model of embryonic haematopoietic development in which current knowledge gaps can be addressed and exploited to enable HSC production. Gastruloids are pluripotent stem cell-derived 3-dimensional (3D) cellular aggregates which recapitulate developmental events in gastrulation and early organogenesis with spatial and temporal precision. Gastruloids self-organise multi-tissue structures upon minimal and controlled external cues, and are amenable to live imaging, screening, scaling and physicochemical manipulation to understand and translate tissue formation. In this review, we consider the haematopoietic potential of gastruloids and review early strategies to enhance blood progenitor and HSC production. We highlight possible strategies to achieve HSC production from gastruloids, and discuss the potential of gastruloid systems in illuminating current knowledge gaps in HSC specification.

ENTRAIN: integrating trajectory inference and gene regulatory networks with spatial data to co-localize the receptor-ligand interactions that specify cell fate

MOTIVATION

Cell fate is commonly studied by profiling the gene expression of single cells to infer developmental trajectories based on expression similarity, RNA velocity, or statistical mechanical properties. However, current approaches do not recover microenvironmental signals from the cellular niche that drive a differentiation trajectory.

RESULTS

We resolve this with environment-aware trajectory inference (ENTRAIN), a computational method that integrates trajectory inference methods with ligand-receptor pair gene regulatory networks to identify extracellular signals and evaluate their relative contribution towards a differentiation trajectory. The output from ENTRAIN can be superimposed on spatial data to co-localize cells and molecules in space and time to map cell fate potentials to cell-cell interactions. We validate and benchmark our approach on single-cell bone marrow and spatially resolved embryonic neurogenesis datasets to identify known and novel environmental drivers of cellular differentiation.

AVAILABILITY AND IMPLEMENTATION

ENTRAIN is available as a public package at https://github.com/theimagelab/entrain and can be used on both single-cell and spatially resolved datasets.

Collective signalling drives rapid jumping between cell states

Development can proceed in 'fits and starts', with rapid transitions between cell states involving concerted transcriptome-wide changes in gene expression. However, it is not clear how these transitions are regulated in complex cell populations, in which cells receive multiple inputs. We address this issue using Dictyostelium cells undergoing development in their physiological niche. A continuous single cell transcriptomics time series identifies a sharp 'jump' in global gene expression marking functionally different cell states. By simultaneously imaging the physiological dynamics of transcription and signalling, we show the jump coincides with the onset of collective oscillations of cAMP. Optogenetic control of cAMP pulses shows that different jump genes respond to distinct dynamic features of signalling. Late jump gene expression changes are almost completely dependent on cAMP, whereas transcript changes at the onset of the jump require additional input. The coupling of collective signalling with gene expression is a potentially powerful strategy to drive robust cell state transitions in heterogeneous signalling environments. Based on the context of the jump, we also conclude that sharp gene expression transitions may not be sufficient for commitment.

Single-Cell Measurements and Modeling and Computation of Decision-Making Errors in a Molecular Signaling System with Two Output Molecules

A cell constantly receives signals and takes different fates accordingly. Given the uncertainty rendered by signal transduction noise, a cell may incorrectly perceive these signals. It may mistakenly behave as if there is a signal, although there is none, or may miss the presence of a signal that actually exists. In this paper, we consider a signaling system with two outputs, and introduce and develop methods to model and compute key cell decision-making parameters based on the two outputs and in response to the input signal. In the considered system, the tumor necrosis factor (TNF) regulates the two transcription factors, the nuclear factor κB (NFκB) and the activating transcription factor-2 (ATF-2). These two system outputs are involved in important physiological functions such as cell death and survival, viral replication, and pathological conditions, such as autoimmune diseases and different types of cancer. Using the introduced methods, we compute and show what the decision thresholds are, based on the single-cell measured concentration levels of NFκB and ATF-2. We also define and compute the decision error probabilities, i.e., false alarm and miss probabilities, based on the concentration levels of the two outputs. By considering the joint response of the two outputs of the signaling system, one can learn more about complex cellular decision-making processes, the corresponding decision error rates, and their possible involvement in the development of some pathological conditions.

Molecular Time Capsules Enable Transcriptomic Recording in Living Cells

Live-cell transcriptomic recording can help reveal hidden cellular states that precede phenotypic transformation. Here we demonstrate the use of protein-based encapsulation for preserving samples of cytoplasmic RNAs inside living cells. These molecular time capsules (MTCs) can be induced to create time-stamped transcriptome snapshots, preserve RNAs after cellular transitions, and enable retrospective investigations of gene expression programs that drive distinct developmental trajectories. MTCs also open the possibility to uncover transcriptomes in difficult-to-reach conditions.

Single-cell sequencing technology applied to epigenetics for the study of tumor heterogeneity

BACKGROUND

Previous studies have traditionally attributed the initiation of cancer cells to genetic mutations, considering them as the fundamental drivers of carcinogenesis. However, recent research has shed light on the crucial role of epigenomic alterations in various cell types present within the tumor microenvironment, suggesting their potential contribution to tumor formation and progression. Despite these significant findings, the progress in understanding the epigenetic mechanisms regulating tumor heterogeneity has been impeded over the past few years due to the lack of appropriate technical tools and methodologies.

RESULTS

The emergence of single-cell sequencing has enhanced our understanding of the epigenetic mechanisms governing tumor heterogeneity by revealing the distinct epigenetic layers of individual cells (chromatin accessibility, DNA/RNA methylation, histone modifications, nucleosome localization) and the diverse omics (transcriptomics, genomics, multi-omics) at the single-cell level. These technologies provide us with new insights into the molecular basis of intratumoral heterogeneity and help uncover key molecular events and driving mechanisms in tumor development.

CONCLUSION

This paper provides a comprehensive review of the emerging analytical and experimental approaches of single-cell sequencing in various omics, focusing specifically on epigenomics. These approaches have the potential to capture and integrate multiple dimensions of individual cancer cells, thereby revealing tumor heterogeneity and epigenetic features. Additionally, this paper outlines the future trends of these technologies and their current technical limitations.

PI3K signaling through a biochemical systems lens

Following 3 decades of extensive research into PI3K signaling, it is now evidently clear that the underlying network does not equate to a simple ON/OFF switch. This is best illustrated by the multifaceted nature of the many diseases associated with aberrant PI3K signaling, including common cancers, metabolic disease, and rare developmental disorders. However, we are still far from a complete understanding of the fundamental control principles that govern the numerous phenotypic outputs that are elicited by activation of this well-characterized biochemical signaling network, downstream of an equally diverse set of extrinsic inputs. At its core, this is a question on the role of PI3K signaling in cellular information processing and decision making. Here, we review the determinants of accurate encoding and decoding of growth factor signals and discuss outstanding questions in the PI3K signal relay network. We emphasize the importance of quantitative biochemistry, in close integration with advances in single-cell time-resolved signaling measurements and mathematical modeling.

The complex biology of aryl hydrocarbon receptor activation in cancer and beyond

The aryl hydrocarbon receptor (AHR) signaling pathway is a complex regulatory network that plays a critical role in various biological processes, including cellular metabolism, development, and immune responses. The complexity of AHR signaling arises from multiple factors, including the diverse ligands that activate the receptor, the expression level of AHR itself, and its interaction with the AHR nuclear translocator (ARNT). Additionally, the AHR crosstalks with the AHR repressor (AHRR) or other transcription factors and signaling pathways and it can also mediate non-genomic effects. Finally, posttranslational modifications of the AHR and its interaction partners, epigenetic regulation of AHR and its target genes, as well as AHR-mediated induction of enzymes that degrade AHR-activating ligands may contribute to the context-specificity of AHR activation. Understanding the complexity of AHR signaling is crucial for deciphering its physiological and pathological roles and developing therapeutic strategies targeting this pathway. Ongoing research continues to unravel the intricacies of AHR signaling, shedding light on the regulatory mechanisms controlling its diverse functions.

Unbalanced response to growth variations reshapes the cell fate decision landscape

The global regulation of cell growth rate on gene expression perturbs the performance of gene networks, which would impose complex variations on the cell-fate decision landscape. Here we use a simple synthetic circuit of mutual repression that allows a bistable landscape to examine how such global regulation would affect the stability of phenotypic landscape and the accompanying dynamics of cell-fate determination. We show that the landscape experiences a growth-rate-induced bifurcation between monostability and bistability. Theoretical and experimental analyses reveal that this bifurcating deformation of landscape arises from the unbalanced response of gene expression to growth variations. The path of growth transition across the bifurcation would reshape cell-fate decisions. These results demonstrate the importance of growth regulation on cell-fate determination processes, regardless of specific molecular signaling or regulation.

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.

Collective behavior and self-organization in neural rosette morphogenesis

Neural rosettes develop from the self-organization of differentiating human pluripotent stem cells. This process mimics the emergence of the embryonic central nervous system primordium, i.e., the neural tube, whose formation is under close investigation as errors during such process result in severe diseases like spina bifida and anencephaly. While neural tube formation is recognized as an example of self-organization, we still do not understand the fundamental mechanisms guiding the process. Here, we discuss the different theoretical frameworks that have been proposed to explain self-organization in morphogenesis. We show that an explanation based exclusively on stem cell differentiation cannot describe the emergence of spatial organization, and an explanation based on patterning models cannot explain how different groups of cells can collectively migrate and produce the mechanical transformations required to generate the neural tube. We conclude that neural rosette development is a relevant experimental 2D in-vitro model of morphogenesis because it is a multi-scale self-organization process that involves both cell differentiation and tissue development. Ultimately, to understand rosette formation, we first need to fully understand the complex interplay between growth, migration, cytoarchitecture organization, and cell type evolution.

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.

Dynamic network-guided CRISPRi screen identifies CTCF-loop-constrained nonlinear enhancer gene regulatory activity during cell state transitions

Comprehensive enhancer discovery is challenging because most enhancers, especially those contributing to complex diseases, have weak effects on gene expression. Our gene regulatory network modeling identified that nonlinear enhancer gene regulation during cell state transitions can be leveraged to improve the sensitivity of enhancer discovery. Using human embryonic stem cell definitive endoderm differentiation as a dynamic transition system, we conducted a mid-transition CRISPRi-based enhancer screen. We discovered a comprehensive set of enhancers for each of the core endoderm-specifying transcription factors. Many enhancers had strong effects mid-transition but weak effects post-transition, consistent with the nonlinear temporal responses to enhancer perturbation predicted by the modeling. Integrating three-dimensional genomic information, we were able to develop a CTCF-loop-constrained Interaction Activity model that can better predict functional enhancers compared to models that rely on Hi-C-based enhancer-promoter contact frequency. Our study provides generalizable strategies for sensitive and systematic enhancer discovery in both normal and pathological cell state transitions.

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.

ETV2 primes hematoendothelial gene enhancers prior to hematoendothelial fate commitment

Mechanisms underlying distinct specification, commitment, and differentiation phases of cell fate determination remain undefined due to difficulties capturing these processes. Here, we interrogate the activity of ETV2, a transcription factor necessary and sufficient for hematoendothelial differentiation, within isolated fate intermediates. We observe transcriptional upregulation of Etv2 and opening of ETV2-binding sites, indicating new ETV2 binding, in a common cardiac-hematoendothelial progenitor population. Accessible ETV2-binding sites are active at the Etv2 locus but not at other hematoendothelial regulator genes. Hematoendothelial commitment coincides with the activation of a small repertoire of previously accessible ETV2-binding sites at hematoendothelial regulators. Hematoendothelial differentiation accompanies activation of a large repertoire of new ETV2-binding sites and upregulation of hematopoietic and endothelial gene regulatory networks. This work distinguishes specification, commitment, and sublineage differentiation phases of ETV2-dependent transcription and suggests that the shift from ETV2 binding to ETV2-bound enhancer activation, not ETV2 binding to target enhancers, drives hematoendothelial fate commitment.

Time-aligned hourglass gastrulation models in rabbit and mouse

The hourglass model describes the convergence of species within the same phylum to a similar body plan during development; however, the molecular mechanisms underlying this phenomenon in mammals remain poorly described. Here, we compare rabbit and mouse time-resolved differentiation trajectories to revisit this model at single-cell resolution. We modeled gastrulation dynamics using hundreds of embryos sampled between gestation days 6.0 and 8.5 and compared the species using a framework for time-resolved single-cell differentiation-flows analysis. We find convergence toward similar cell-state compositions at E7.5, supported by the quantitatively conserved expression of 76 transcription factors, despite divergence in surrounding trophoblast and hypoblast signaling. However, we observed noticeable changes in specification timing of some lineages and divergence of primordial germ cell programs, which in the rabbit do not activate mesoderm genes. Comparative analysis of temporal differentiation models provides a basis for studying the evolution of gastrulation dynamics across mammals.

Regulation of chromatin organization during animal regeneration

Activation of regeneration upon tissue damages requires the activation of many developmental genes responsible for cell proliferation, migration, differentiation, and tissue patterning. Ample evidence revealed that the regulation of chromatin organization functions as a crucial mechanism for establishing and maintaining cellular identity through precise control of gene transcription. The alteration of chromatin organization can lead to changes in chromatin accessibility and/or enhancer-promoter interactions. Like embryogenesis, each stage of tissue regeneration is accompanied by dynamic changes of chromatin organization in regeneration-responsive cells. In the past decade, many studies have been conducted to investigate the contribution of chromatin organization during regeneration in various tissues, organs, and organisms. A collection of chromatin regulators were demonstrated to play critical roles in regeneration. In this review, we will summarize the progress in the understanding of chromatin organization during regeneration in different research organisms and discuss potential common mechanisms responsible for the activation of regeneration response program.

Self-Regulated Symmetry Breaking Model for Stem Cell Differentiation

In conventional disorder-order phase transitions, a system shifts from a highly symmetric state, where all states are equally accessible (disorder) to a less symmetric state with a limited number of available states (order). This transition may occur by varying a control parameter that represents the intrinsic noise of the system. It has been suggested that stem cell differentiation can be considered as a sequence of such symmetry-breaking events. Pluripotent stem cells, with their capacity to develop into any specialized cell type, are considered highly symmetric systems. In contrast, differentiated cells have lower symmetry, as they can only carry out a limited number of functions. For this hypothesis to be valid, differentiation should emerge collectively in stem cell populations. Additionally, such populations must have the ability to self-regulate intrinsic noise and navigate through a critical point where spontaneous symmetry breaking (differentiation) occurs. This study presents a mean-field model for stem cell populations that considers the interplay of cell-cell cooperativity, cell-to-cell variability, and finite-size effects. By introducing a feedback mechanism to control intrinsic noise, the model can self-tune through different bifurcation points, facilitating spontaneous symmetry breaking. Standard stability analysis showed that the system can potentially differentiate into several cell types mathematically expressed as stable nodes and limit cycles. The existence of a Hopf bifurcation in our model is discussed in light of stem cell differentiation.

Structure and function in artificial, zebrafish and human neural networks

Network science provides a set of tools for the characterization of the structure and functional behavior of complex systems. Yet a major problem is to quantify how the structural domain is related to the dynamical one. In other words, how the diversity of dynamical states of a system can be predicted from the static network structure? Or the reverse problem: starting from a set of signals derived from experimental recordings, how can one discover the network connections or the causal relations behind the observed dynamics? Despite the advances achieved over the last two decades, many challenges remain concerning the study of the structure-dynamics interplay of complex systems. In neuroscience, progress is typically constrained by the low spatio-temporal resolution of experiments and by the lack of a universal inferring framework for empirical systems. To address these issues, applications of network science and artificial intelligence to neural data have been rapidly growing. In this article, we review important recent applications of methods from those fields to the study of the interplay between structure and functional dynamics of human and zebrafish brain. We cover the selection of topological features for the characterization of brain networks, inference of functional connections, dynamical modeling, and close with applications to both the human and zebrafish brain. This review is intended to neuroscientists who want to become acquainted with techniques from network science, as well as to researchers from the latter field who are interested in exploring novel application scenarios in neuroscience.

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.

The diverging epigenomic landscapes of honeybee queens and workers revealed by multiomic sequencing

The role of the epigenome in phenotypic plasticity is unclear presently. Here we used a multiomics approach to explore the nature of the epigenome in developing honey bee (Apis mellifera) workers and queens. Our data clearly showed distinct queen and worker epigenomic landscapes during the developmental process. Differences in gene expression between workers and queens become more extensive and more layered during the process of development. Genes known to be important for caste differentiation were more likely to be regulated by multiple epigenomic systems than other differentially expressed genes. We confirmed the importance of two candidate genes for caste differentiation by using RNAi to manipulate the expression of two genes that differed in expression between workers and queens were regulated by multiple epigenomic systems. For both genes the RNAi manipulation resulted in a decrease in weight and fewer ovarioles of newly emerged queens compared to controls. Our data show that the distinct epigenomic landscapes of worker and queen bees differentiate during the course of larval development.

Dynamic network-guided CRISPRi screen reveals CTCF loop-constrained nonlinear enhancer-gene regulatory activity in cell state transitions

Comprehensive enhancer discovery is challenging because most enhancers, especially those affected in complex diseases, have weak effects on gene expression. Our network modeling revealed that nonlinear enhancer-gene regulation during cell state transitions can be leveraged to improve the sensitivity of enhancer discovery. Utilizing hESC definitive endoderm differentiation as a dynamic transition system, we conducted a mid-transition CRISPRi-based enhancer screen. The screen discovered a comprehensive set of enhancers (4 to 9 per locus) for each of the core endoderm lineage-specifying transcription factors, and many enhancers had strong effects mid-transition but weak effects post-transition. Through integrating enhancer activity measurements and three-dimensional enhancer-promoter interaction information, we were able to develop a CTCF loop-constrained Interaction Activity (CIA) model that can better predict functional enhancers compared to models that rely on Hi-C-based enhancer-promoter contact frequency. Our study provides generalizable strategies for sensitive and more comprehensive enhancer discovery in both normal and pathological cell state transitions.

Well-TEMP-seq as a microwell-based strategy for massively parallel profiling of single-cell temporal RNA dynamics

Single-cell RNA sequencing (scRNA-seq) reveals the transcriptional heterogeneity of cells, but the static snapshots fail to reveal the time-resolved dynamics of transcription. Herein, we develop Well-TEMP-seq, a high-throughput, cost-effective, accurate, and efficient method for massively parallel profiling the temporal dynamics of single-cell gene expression. Well-TEMP-seq combines metabolic RNA labeling with scRNA-seq method Well-paired-seq to distinguish newly transcribed RNAs marked by T-to-C substitutions from pre-existing RNAs in each of thousands of single cells. The Well-paired-seq chip ensures a high single cell/barcoded bead pairing rate (~80%) and the improved alkylation chemistry on beads greatly alleviates chemical conversion-induced cell loss (~67.5% recovery). We further apply Well-TEMP-seq to profile the transcriptional dynamics of colorectal cancer cells exposed to 5-AZA-CdR, a DNA-demethylating drug. Well-TEMP-seq unbiasedly captures the RNA dynamics and outperforms the splicing-based RNA velocity method. We anticipate that Well-TEMP-seq will be broadly applicable to unveil the dynamics of single-cell gene expression in diverse biological processes.

Cell state-directed therapy - epigenetic modulation of gene transcription demonstrated with a quantitative systems pharmacology model of temozolomide

Cancer therapy continues to be plagued by modest therapeutic advances. This is particularly evident in glioblastoma multiforme (GBM) wherein treatment failures are attributed to intratumoral heterogeneity (ITH), a dynamic process of cell state transitions or plasticity. To address ITH, we introduce the concept of cell state-directed (CSD) therapy through a quantitative systems pharmacology model of temozolomide (TMZ), a cornerstone of GBM drug therapy. The model consisting of multiple modules incorporated an epigenetic-based gene transcription-translation module that enabled CSD therapy. Numerous model simulations were conducted to demonstrate the potential impact of CSD therapy on TMZ activity. The simulations included those based on global sensitivity analyses to identify fragile nodes - MDM2 and XIAP - in the network, and also how an epigenetic modifier (birabresib) could overcome a mechanism of TMZ resistance. The positive results of CSD therapy on TMZ activity supports continued efforts to develop CSD therapy as a new anticancer approach.

LSM1-mediated Major Satellite RNA decay is required for nonequilibrium histone H3.3 incorporation into parental pronuclei

Epigenetic reprogramming of the parental genome is essential for zygotic genome activation and subsequent embryo development in mammals. Asymmetric incorporation of histone H3 variants into the parental genome has been observed previously, but the underlying mechanism remains elusive. In this study, we discover that RNA-binding protein LSM1-mediated major satellite RNA decay plays a central role in the preferential incorporation of histone variant H3.3 into the male pronucleus. Knockdown of Lsm1 disrupts nonequilibrium pronucleus histone incorporation and asymmetric H3K9me3 modification. Subsequently, we find that LSM1 mainly targets major satellite repeat RNA (MajSat RNA) for decay and that accumulated MajSat RNA in Lsm1-depleted oocytes leads to abnormal incorporation of H3.1 into the male pronucleus. Knockdown of MajSat RNA reverses the anomalous histone incorporation and modifications in Lsm1-knockdown zygotes. Our study therefore reveals that accurate histone variant incorporation and incidental modifications in parental pronuclei are specified by LSM1-dependent pericentromeric RNA decay.