Clearing the slate: RNA turnover to enable cell state switching?

The distribution of mRNA in tissue is determined by the balance between transcription and decay. Understanding the control of RNA decay during development has been somewhat neglected compared with transcriptional control. Here, we explore the potential for mRNA decay to trigger rapid cell state transitions during development, comparing a bistable switch model of cell state conversion with experimental evidence from different developmental systems. We also consider another potential role for large-scale RNA decay that has emerged from studies of stress-induced cell state transitions, in which removal of mRNA unblocks the translation machinery to prioritise the synthesis of proteins that establish the new cell state.

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.

Zfp281 orchestrates interconversion of pluripotent states by engaging Ehmt1 and Zic2.

Developmental cell fate specification is a unidirectional process that can be reverted in response to injury or experimental reprogramming. Whether differentiation and de-differentiation trajectories intersect mechanistically is unclear. Here, we performed comparative screening in lineage-related mouse naïve embryonic stem cells (ESCs) and primed epiblast stem cells (EpiSCs), and identified the constitutively expressed zinc finger transcription factor (TF) Zfp281 as a bidirectional regulator of cell state interconversion. We showed that subtle chromatin binding changes in differentiated cells translate into activation of the histone H3 lysine 9 (H3K9) methyltransferase Ehmt1 and stabilization of the zinc finger TF Zic2 at enhancers and promoters. Genetic gain-of-function and loss-of-function experiments confirmed a critical role of Ehmt1 and Zic2 downstream of Zfp281 both in driving exit from the ESC state and in restricting reprogramming of EpiSCs. Our study reveals that cell type-invariant chromatin association of Zfp281 provides an interaction platform for remodeling the cis-regulatory network underlying cellular plasticity.

Transition from totipotency to pluripotency in mice: insights into molecular mechanisms.

Totipotency is the ability of a single cell to develop into a full organism and, in mammals, is strictly associated with the early stages of development following fertilization. This unlimited developmental potential becomes quickly restricted as embryonic cells transition into a pluripotent state. The loss of totipotency seems a consequence of the zygotic genome activation (ZGA), a process that determines the switch from maternal to embryonic transcription, which in mice takes place following the first cleavage. ZGA confers to the totipotent cell a transient transcriptional profile characterized by the expression of stage-specific genes and a set of transposable elements that prepares the embryo for subsequent development. The timely silencing of this transcriptional program during the exit from totipotency is required to ensure proper development. Importantly, the molecular mechanisms regulating the transition from totipotency to pluripotency have remained elusive due to the scarcity of embryonic material. However, the development of new in vitro totipotent-like models together with advances in low-input genome-wide technologies, are providing a better mechanistic understanding of how this important transition is achieved. This review summarizes the current knowledge on the molecular determinants that regulate the exit from totipotency.

State transition and intercellular communication of synovial fibroblasts in response to chronic and acute shoulder injuries unveiled by single-cell transcriptomic analyses.

PURPOSE: We aimed to investigate the heterogeneity of synovial fibroblasts and their potential to undergo cell state transitions at the resolution of single cells. MATERIALS AND METHODS: We employed the single-cell RNA sequencing (scRNA-seq) approach to comprehensively map the cellular landscape of the shoulder synovium in individuals with chronic rotator cuff tears (RCTs) and acute proximal humerus fractures (PHFs). Utilizing unbiased clustering analysis, we successfully identified distinct subpopulations of fibroblasts within the synovial environment. We utilized Monocle 3 to delineate the trajectory of synovial fibroblast state transition. And we used CellPhone DB v2.1.0 to predict cell-cell communication patterns within the synovial microenvironment. RESULTS: We identified eight main cell clusters in the shoulder synovium. Unbiased clustering analysis identified four synovial fibroblast subpopulations, with diverse biological functions associated with protein secretion, ECM remodeling, inflammation regulation and cell division. Lining, mesenchymal, pro-inflammatory and proliferative fibroblasts subsets were identified. Combining the results from StemID and characteristic gene features, mesenchymal fibroblasts exhibited characteristics of fibroblast progenitor cells. The trajectory of synovial fibroblast state transition showed a transition from mesenchymal to pro-inflammatory and lining phenotypes. In addition, the cross talk between fibroblast subclusters increased in degenerative shoulder diseases compared to acute trauma. CONCLUSION: We successfully generated the scRNA-seq transcriptomic atlas of the shoulder synovium, which provides a comprehensive understanding of the heterogeneity of synovial fibroblasts, their potential to undergo state transitions, and their intercellular communication in the context of chronic degenerative and acute traumatic shoulder diseases.

spliceJAC: transition genes and state-specific gene regulation from single-cell transcriptome data.

Extracting dynamical information from single-cell transcriptomics is a novel task with the promise to advance our understanding of cell state transition and interactions between genes. Yet, theory-oriented, bottom-up approaches that consider differences among cell states are largely lacking. Here, we present spliceJAC, a method to quantify the multivariate mRNA splicing from single-cell RNA sequencing (scRNA-seq). spliceJAC utilizes the unspliced and spliced mRNA count matrices to constructs cell state-specific gene-gene regulatory interactions and applies stability analysis to predict putative driver genes critical to the transitions between cell states. By applying spliceJAC to biological systems including pancreas endothelium development and epithelial-mesenchymal transition (EMT) in A549 lung cancer cells, we predict genes that serve specific signaling roles in different cell states, recover important differentially expressed genes in agreement with pre-existing analysis, and predict new transition genes that are either exclusive or shared between different cell state transitions.

SOL1 and SOL2 regulate fate transition and cell divisions in the Arabidopsis stomatal lineage.

In the Arabidopsis stomatal lineage, cells transit through several distinct precursor identities, each characterized by unique cell division behaviors. Flexibility in the duration of these precursor phases enables plants to alter leaf size and stomatal density in response to environmental conditions; however, transitions between phases must be complete and unidirectional to produce functional and correctly patterned stomata. Among direct transcriptional targets of the stomatal initiating factor SPEECHLESS, a pair of genes, SOL1 and SOL2, are required for effective transitions in the lineage. We show that these two genes, which are homologs of the LIN54 DNA-binding components of the mammalian DREAM complex, are expressed in a cell cycle-dependent manner and regulate cell fate and division properties in the self-renewing early lineage. In the terminal division of the stomatal lineage, however, these two proteins appear to act in opposition to their closest paralog, TSO1, revealing complexity in the gene family that may enable customization of cell divisions in coordination with development.

Biophysical and biochemical attributes of hybrid epithelial/mesenchymal phenotypes.

The epithelial-mesenchymal transition (EMT) is a biological phenomenon associated with explicit phenotypic and molecular changes in cellular traits. Unlike the earlier-held popular belief of it being a binary process, EMT is now thought of as a landscape including diverse hybrid E/M phenotypes manifested by varying degrees of the transition. These hybrid cells can co-express both epithelial and mesenchymal markers and/or functional traits, and can possess the property of collective cell migration, enhanced tumor-initiating ability, and immune/targeted therapy-evasive features, all of which are often associated with worse patient outcomes. These characteristics of the hybrid E/M cells have led to a surge in studies that map their biophysical and biochemical hallmarks that can be helpful in exploiting their therapeutic vulnerabilities. This review discusses recent advances made in investigating hybrid E/M phenotype(s) from diverse biophysical and biochemical aspects by integrating live cell-imaging, cellular morphology quantification and mathematical modeling, and highlights a set of questions that remain unanswered about the dynamics of hybrid E/M states.

Single-cell analysis of cell fate bifurcation in the chordate Ciona.

BACKGROUND: Inductive signaling interactions between different cell types are a major mechanism for the further diversification of embryonic cell fates. Most blastomeres in the model chordate Ciona robusta become restricted to a single predominant fate between the 64-cell and mid-gastrula stages. The deeply stereotyped and well-characterized Ciona embryonic cell lineages allow the transcriptomic analysis of newly established cell types very early in their divergence from sibling cell states without the pseudotime inference needed in the analysis of less synchronized cell populations. This is the first ascidian study to use droplet scRNAseq with large numbers of analyzed cells as early as the 64-cell stage when major lineages such as primary notochord first become fate restricted. RESULTS AND CONCLUSIONS: We identify 59 distinct cell states, including new subregions of the b-line neural lineage and the early induction of the tail tip epidermis. We find that 34 of these cell states are directly or indirectly dependent on MAPK-mediated signaling critical to early Ciona patterning. Most of the MAPK-dependent bifurcations are canalized with the signal-induced cell fate lost upon MAPK inhibition, but the posterior endoderm is unique in being transformed into a novel state expressing some but not all markers of both endoderm and muscle. Divergent gene expression between newly bifurcated sibling cell types is dominated by upregulation in the induced cell type. The Ets family transcription factor Elk1/3/4 is uniquely upregulated in nearly all the putatively direct inductions. Elk1/3/4 upregulation together with Ets transcription factor binding site enrichment analysis enables inferences about which bifurcations are directly versus indirectly controlled by MAPK signaling. We examine notochord induction in detail and find that the transition between a Zic/Ets-mediated regulatory state and a Brachyury/FoxA-mediated regulatory state is unexpectedly late. This supports a "broad-hourglass" model of cell fate specification in which many early tissue-specific genes are induced in parallel to key tissue-specific transcriptional regulators via the same set of transcriptional inputs.

Drivers of dynamic intratumor heterogeneity and phenotypic plasticity.

Cancer is a clonal disease, i.e., all tumor cells within a malignant lesion trace their lineage back to a precursor somatic cell that acquired oncogenic mutations during development and aging. And yet, those tumor cells tend to have genetic and nongenetic variations among themselves-which is denoted as intratumor heterogeneity. Although some of these variations are inconsequential, others tend to contribute to cell state transition and phenotypic heterogeneity, providing a substrate for somatic evolution. Tumor cell phenotypes can dynamically change under the influence of genetic mutations, epigenetic modifications, and microenvironmental contexts. Although epigenetic and microenvironmental changes are adaptive, genetic mutations are usually considered permanent. Emerging reports suggest that certain classes of genetic alterations show extensive reversibility in tumors in clinically relevant timescales, contributing as major drivers of dynamic intratumor heterogeneity and phenotypic plasticity. Dynamic heterogeneity and phenotypic plasticity can confer resistance to treatment, promote metastasis, and enhance evolvability in cancer. Here, we first highlight recent efforts to characterize intratumor heterogeneity at genetic, epigenetic, and microenvironmental levels. We then discuss phenotypic plasticity and cell state transition by tumor cells, under the influence of genetic and nongenetic determinants and their clinical significance in classification of tumors and therapeutic decision-making.

Cross-species integration of single-cell RNA-seq resolved alveolar-epithelial transitional states in idiopathic pulmonary fibrosis.

Single-cell transcriptomics analyses of the fibrotic lung uncovered two cell states critical to lung injury recovery in the alveolar epithelium-a reparative transitional cell state in the mouse and a disease-specific cell state (KRT5-/KRT17+) in human idiopathic pulmonary fibrosis (IPF). The murine transitional cell state lies between the differentiation from type 2 (AT2) to type 1 pneumocyte (AT1), and the human KRT5-/KRT17+ cell state may arise from the dysregulation of this differentiation process. We review major findings of single-cell transcriptomics analyses of the fibrotic lung and reanalyzed data from seven single-cell RNA sequencing studies of human and murine models of IPF, focusing on the alveolar epithelium. Our comparative and cross-species single-cell transcriptomics analyses allowed us to further delineate the differentiation trajectories from AT2 to AT1 and AT2 to the KRT5-/KRT17+ cell state. We observed AT1 cells in human IPF retain the transcriptional signature of the murine transitional cell state. Using pseudotime analysis, we recapitulated the differentiation trajectories from AT2 to AT1 and from AT2 to KRT5-/KRT17+ cell state in multiple human IPF studies. We further delineated transcriptional programs underlying cell-state transitions and determined the molecular phenotypes at terminal differentiation. We hypothesize that in addition to the reactivation of developmental programs (SOX4, SOX9), senescence (TP63, SOX4) and the Notch pathway (HES1) are predicted to steer intermediate progenitors to the KRT5-/KRT17+ cell state. Our analyses suggest that activation of SMAD3 later in the differentiation process may explain the fibrotic molecular phenotype typical of KRT5-/KRT17+ cells.

Single-cell analysis resolves the cell state transition and signaling dynamics associated with melanoma drug-induced resistance.

Continuous BRAF inhibition of BRAF mutant melanomas triggers a series of cell state changes that lead to therapy resistance and escape from immune control before establishing acquired resistance genetically. We used genome-wide transcriptomics and single-cell phenotyping to explore the response kinetics to BRAF inhibition for a panel of patient-derived BRAFV600 -mutant melanoma cell lines. A subset of plastic cell lines, which followed a trajectory covering multiple known cell state transitions, provided models for more detailed biophysical investigations. Markov modeling revealed that the cell state transitions were reversible and mediated by both Lamarckian induction and nongenetic Darwinian selection of drug-tolerant states. Single-cell functional proteomics revealed activation of certain signaling networks shortly after BRAF inhibition, and before the appearance of drug-resistant phenotypes. Drug targeting those networks, in combination with BRAF inhibition, halted the adaptive transition and led to prolonged growth inhibition in multiple patient-derived cell lines.

High-Throughput Automated Single-Cell Imaging Analysis Reveals Dynamics of Glioblastoma Stem Cell Population During State Transition.

Cancer stem cells (CSCs) are a heterogeneous and dynamic self-renewing population that stands at the top of tumor cellular hierarchy and contribute to tumor recurrence and therapeutic resistance. As methods of CSC isolation and functional interrogation advance, there is a need for a reliable and accessible quantitative approach to assess heterogeneity and state transition dynamics in CSCs. We developed a high-throughput automated single cell imaging analysis (HASCIA) approach for the quantitative assessment of protein expression with single-cell resolution and applied the method to investigate spatiotemporal factors that influence CSC state transition using glioblastoma (GBM) CSCs (GSCs) as a model system. We were able to validate the quantitative nature of this approach through comparison of the protein expression levels determined by HASCIA to those determined by immunoblotting. A virtue of HASCIA was exemplified by detection of a subpopulation of SOX2-low cells, which expanded in fraction size during state transition. HASCIA also revealed that GSCs were committed to loose stem cell state at an earlier time point than the average SOX2 level decreased. Functional assessment of stem cell frequency in combination with the quantification of SOX2 expression by HASCIA defined a stable cutoff of SOX2 expression level for stem cell state. We also developed an approach to assess local cell density and found that denser monolayer areas possess higher average levels of SOX2, higher cell diversity, and a presence of a sub-population of slowly proliferating SOX2-low GSCs. HASCIA is an open source software that facilitates understanding the dynamics of heterogeneous cell population such as that of GSCs and their progeny. It is a powerful and easy-to-use image analysis and statistical analysis tool available at https://hascia.lerner.ccf.org. © 2019 International Society for Advancement of Cytometry.

Pluripotency Deconstructed.

Pluripotency denotes the flexible capacity of single cells to give rise to all somatic lineages and typically also the germline. Mouse ES cells and post-implantation epiblast-derived stem cells (EpiSC) are widely used pluripotent cell culture systems. These two in vitro stem cell types have divergent characteristics. They are considered as representative of distinct developmental stages, distinguished by using the terms "naïve" and "primed". A binary description is an over-simplification, however. Here, we discuss an intermediate stage of pluripotency that we term "formative". Formative pluripotency features a gene regulatory network switch from the naïve state and comprises capacitation of enhancers, signaling pathways and epigenetic machinery in order to install competence for lineage specification.

Stomagenesis versus myogenesis: Parallels in intrinsic and extrinsic regulation of transcription factor mediated specialized cell-type differentiation in plants and animals.

Although the last common unicellular ancestor of plants and animals diverged several billion years ago, and while having developed unique developmental programs that facilitate differentiation and proliferation specific to plant and animal systems, there still exists a high degree of conservation in the logic regulating these developmental processes within these two seemingly diverse kingdoms. Stomatal differentiation in plants involves a series of orchestrated cell division events mediated by a family of closely related bHLH transcription factors (TFs) to create a pair of mature guard cells. These TFs are in turn regulated by a number of upstream signaling components that ultimately function to achieve lineage specific differentiation and organized tissue patterning on the plant epidermis. The logic involved in the specification of the myogenic differentiation program in animals is intriguingly similar to stomatal differentiation in plants: Closely-related myogenic bHLHs, known as MRFs (Myogenic Regulatory Factors) provide lineage specificity essential for cell-fate determination. These MRFs, similar to the bHLHs in plants, are regulated by several upstream signaling cascades that succinctly regulate each differentiation step, leading to the production of mature muscle fibers. This review aims at providing a perspective on the emerging parallels in the logic employed by key bHLH transcription factors and their upstream signaling components that function to precisely regulate key cell-state transition events in the stomatal as well as myogenic cell lineages.

Copper and iron orchestrate cell-state transitions in cancer and immunity.

Whereas genetic mutations can alter cell properties, nongenetic mechanisms can drive rapid and reversible adaptations to changes in their physical environment, a phenomenon termed 'cell-state transition'. Metals, in particular copper and iron, have been shown to be rate-limiting catalysts of cell-state transitions controlling key chemical reactions in mitochondria and the cell nucleus, which govern metabolic and epigenetic changes underlying the acquisition of distinct cell phenotypes. Acquisition of a distinct cell identity, independently of genetic alterations, is an underlying phenomenon of various biological processes, including development, inflammation, erythropoiesis, aging, and cancer. Here, mechanisms that have been uncovered related to the role of these metals in the regulation of cell plasticity are described, illustrating how copper and iron can be exploited for therapeutic intervention.

Molecular cascade reveals sequential milestones underlying hippocampal neural stem cell development into an adult state.

Quiescent adult neural stem cells (NSCs) in the mammalian brain arise from proliferating NSCs during development. Beyond acquisition of quiescence, an adult NSC hallmark, little is known about the process, milestones, and mechanisms underlying the transition of developmental NSCs to an adult NSC state. Here, we performed targeted single-cell RNA-seq analysis to reveal the molecular cascade underlying NSC development in the early postnatal mouse dentate gyrus. We identified two sequential steps, first a transition to quiescence followed by further maturation, each of which involved distinct changes in metabolic gene expression. Direct metabolic analysis uncovered distinct milestones, including an autophagy burst before NSC quiescence acquisition and cellular reactive oxygen species level elevation along NSC maturation. Functionally, autophagy is important for the NSC transition to quiescence during early postnatal development. Together, our study reveals a multi-step process with defined milestones underlying establishment of the adult NSC pool in the mammalian brain.

Sequential drug treatment targeting cell cycle and cell fate regulatory programs blocks non-genetic cancer evolution in acute lymphoblastic leukemia.

BACKGROUND: Targeted therapies exploiting vulnerabilities of cancer cells hold promise for improving patient outcome and reducing side-effects of chemotherapy. However, efficacy of precision therapies is limited in part because of tumor cell heterogeneity. A better mechanistic understanding of how drug effect is linked to cancer cell state diversity is crucial for identifying effective combination therapies that can prevent disease recurrence. RESULTS: Here, we characterize the effect of G2/M checkpoint inhibition in acute lymphoblastic leukemia (ALL) and demonstrate that WEE1 targeted therapy impinges on cell fate decision regulatory circuits. We find the highest inhibition of recovery of proliferation in ALL cells with KMT2A-rearrangements. Single-cell RNA-seq and ATAC-seq of RS4;11 cells harboring KMT2A::AFF1, treated with the WEE1 inhibitor AZD1775, reveal diversification of cell states, with a fraction of cells exhibiting strong activation of p53-driven processes linked to apoptosis and senescence, and disruption of a core KMT2A-RUNX1-MYC regulatory network. In this cell state diversification induced by WEE1 inhibition, a subpopulation transitions to a drug tolerant cell state characterized by activation of transcription factors regulating pre-B cell fate, lipid metabolism, and pre-BCR signaling in a reversible manner. Sequential treatment with BCR-signaling inhibitors dasatinib, ibrutinib, or perturbing metabolism by fatostatin or AZD2014 effectively counteracts drug tolerance by inducing cell death and repressing stemness markers. CONCLUSIONS: Collectively, our findings provide new insights into the tight connectivity of gene regulatory programs associated with cell cycle and cell fate regulation, and a rationale for sequential administration of WEE1 inhibitors with low toxicity inhibitors of pre-BCR signaling or metabolism.

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.

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.