Multilineage gene expression precedes commitment in the hemopoietic system

Metabolic adaptation pilots the differentiation of human hematopoietic cells

A continuous supply of energy is an essential prerequisite for survival and represents the highest priority for the cell. We hypothesize that cell differentiation is a process of optimization of energy flow in a changing environment through phenotypic adaptation. The mechanistic basis of this hypothesis is provided by the established link between core energy metabolism and epigenetic covalent modifications of chromatin. This theory predicts that early metabolic perturbations impact subsequent differentiation. To test this, we induced transient metabolic perturbations in undifferentiated human hematopoietic cells using pharmacological inhibitors targeting key metabolic reactions. We recorded changes in chromatin structure and gene expression, as well as phenotypic alterations by single-cell ATAC and RNA sequencing, time-lapse microscopy, and flow cytometry. Our observations suggest that these metabolic perturbations are shortly followed by alterations in chromatin structure, leading to changes in gene expression. We also show that these transient fluctuations alter the differentiation potential of the cells.

Early life imprinting of intestinal immune tolerance and tissue homeostasis

Besides its canonical role in protecting the host from pathogens, the immune system plays an arguably equally important role in maintaining tissue homeostasis. Within barrier tissues that interface with the external microenvironment, induction of immune tolerance to innocuous antigens, such as commensal, dietary, and environmental antigens, is key to establishing immune homeostasis. The early postnatal period represents a critical window of opportunity in which parallel development of the tissue, immune cells, and microbiota allows for reciprocal regulation that shapes the long-term immunological tone of the tissue and subsequent risk of immune-mediated diseases. During early infancy, the immune system appears to sacrifice pro-inflammatory functions, prioritizing the establishment of tissue tolerance. In this review, we discuss mechanisms underlying early life windows for intestinal tolerance with a focus on newly identified RORγt+ antigen-presenting cells-Thetis cells-and highlight the role of the intestinal microenvironment in shaping intestinal immune system development and tolerance.

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.

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'.

Mapping functional to morphological variation reveals the basis of regional extracellular matrix subversion and nerve invasion in pancreatic cancer

Intratumor morphological heterogeneity of pancreatic ductal adenocarcinoma (PDAC) predicts clinical outcomes but is only partially understood at the molecular level. To elucidate the gene expression programs underpinning intratumor morphological variation in PDAC, we investigated and deconvoluted at single cell level the molecular profiles of histologically distinct clusters of PDAC cells. We identified three major morphological and functional variants that co-exist in varying proportions in all PDACs, display limited genetic diversity, and are associated with a distinct organization of the extracellular matrix: a glandular variant with classical ductal features; a transitional variant displaying abortive ductal structures and mixed endodermal and myofibroblast-like gene expression; and a poorly differentiated variant lacking ductal features and basement membrane, and showing neuronal lineage priming. Ex vivo and in vitro evidence supports the occurrence of dynamic transitions among these variants in part influenced by extracellular matrix composition and stiffness and associated with local, specifically neural, invasion.

Differentiation is accompanied by a progressive loss in transcriptional memory

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Neutrophil Extracellular DNA Traps in Response to Infection or Inflammation, and the Roles of Platelet Interactions

Neutrophils present the host's first line of defense against bacterial infections. These immune effector cells are mobilized rapidly to destroy invading pathogens by (a) reactive oxygen species (ROS)-mediated oxidative bursts and (b) via phagocytosis. In addition, their antimicrobial service is capped via a distinct cell death mechanism, by the release of their own decondensed nuclear DNA, supplemented with a variety of embedded proteins and enzymes. The extracellular DNA meshwork ensnares the pathogenic bacteria and neutralizes them. Such neutrophil extracellular DNA traps (NETs) have the potential to trigger a hemostatic response to pathogenic infections. The web-like chromatin serves as a prothrombotic scaffold for platelet adhesion and activation. What is less obvious is that platelets can also be involved during the initial release of NETs, forming heterotypic interactions with neutrophils and facilitating their responses to pathogens. Together, the platelet and neutrophil responses can effectively localize an infection until it is cleared. However, not all microbial infections are easily cleared. Certain pathogenic organisms may trigger dysregulated platelet-neutrophil interactions, with a potential to subsequently propagate thromboinflammatory processes. These may also include the release of some NETs. Therefore, in order to make rational intervention easier, further elucidation of platelet, neutrophil, and pathogen interactions is still needed.

Transcription factor co-expression mediates lineage priming for embryonic and extra-embryonic differentiation

In early mammalian development, cleavage stage blastomeres and inner cell mass (ICM) cells co-express embryonic and extra-embryonic transcriptional determinants. Using a protein-based double reporter we identify an embryonic stem cell (ESC) population that co-expresses the extra-embryonic factor GATA6 alongside the embryonic factor SOX2. Based on single cell transcriptomics, we find this population resembles the unsegregated ICM, exhibiting enhanced differentiation potential for endoderm while maintaining epiblast competence. To relate transcription factor binding in these cells to future fate, we describe a complete enhancer set in both ESCs and naive extra-embryonic endoderm stem cells and assess SOX2 and GATA6 binding at these elements in the ICM-like ESC sub-population. Both factors support cooperative recognition in these lineages, with GATA6 bound alongside SOX2 on a fraction of pluripotency enhancers and SOX2 alongside GATA6 more extensively on endoderm enhancers, suggesting that cooperative binding between these antagonistic factors both supports self-renewal and prepares progenitor cells for later differentiation.

Chromatin priming elements direct tissue-specific gene activity before hematopoietic specification

Tissue-specific gene regulation during development involves the interplay between transcription factors and epigenetic regulators binding to enhancer and promoter elements. The pattern of active enhancers defines the cellular differentiation state. However, developmental gene activation involves a previous step called chromatin priming which is not fully understood. We recently developed a genome-wide functional assay that allowed us to functionally identify enhancer elements integrated in chromatin regulating five stages spanning the in vitro differentiation of embryonic stem cells to blood. We also measured global chromatin accessibility, histone modifications, and transcription factor binding. The integration of these data identified and characterised cis-regulatory elements which become activated before the onset of gene expression, some of which are primed in a signalling-dependent fashion. Deletion of such a priming element leads to a delay in the up-regulation of its associated gene in development. Our work uncovers the details of a complex network of regulatory interactions with the dynamics of early chromatin opening being at the heart of dynamic tissue-specific gene expression control.

From Genotype to Phenotype: How Enhancers Control Gene Expression and Cell Identity in Hematopoiesis

Blood comprises a wide array of specialized cells, all of which share the same genetic information and ultimately derive from the same precursor, the hematopoietic stem cell (HSC). This diversity of phenotypes is underpinned by unique transcriptional programs gradually acquired in the process known as hematopoiesis. Spatiotemporal regulation of gene expression depends on many factors, but critical among them are enhancers-sequences of DNA that bind transcription factors and increase transcription of genes under their control. Thus, hematopoiesis involves the activation of specific enhancer repertoires in HSCs and their progeny, driving the expression of sets of genes that collectively determine morphology and function. Disruption of this tightly regulated process can have catastrophic consequences: in hematopoietic malignancies, dysregulation of transcriptional control by enhancers leads to misexpression of oncogenes that ultimately drive transformation. This review attempts to provide a basic understanding of enhancers and their role in transcriptional regulation, with a focus on normal and malignant hematopoiesis. We present examples of enhancers controlling master regulators of hematopoiesis and discuss the main mechanisms leading to enhancer dysregulation in leukemia and lymphoma.

A bipartite function of ESRRB can integrate signaling over time to balance self-renewal and differentiation

Cooperative DNA binding of transcription factors (TFs) integrates the cellular context to support cell specification during development. Naive mouse embryonic stem cells are derived from early development and can sustain their pluripotent identity indefinitely. Here, we ask whether TFs associated with pluripotency evolved to directly support this state or if the state emerges from their combinatorial action. NANOG and ESRRB are key pluripotency factors that co-bind DNA. We find that when both factors are expressed, ESRRB supports pluripotency. However, when NANOG is absent, ESRRB supports a bistable culture of cells with an embryo-like primitive endoderm identity ancillary to pluripotency. The stoichiometry between NANOG and ESRRB allows quantitative titration of this differentiation, and in silico modeling of bipartite ESRRB activity suggests it safeguards plasticity in differentiation. Thus, the concerted activity of cooperative TFs can transform their effect to sustain intermediate cell identities and allow ex vivo expansion of immortal stem cells. A record of this paper's transparent peer review process is included in the supplemental information.

Shaping faces: genetic and epigenetic control of craniofacial morphogenesis

Major differences in facial morphology distinguish vertebrate species. Variation of facial traits underlies the uniqueness of human individuals, and abnormal craniofacial morphogenesis during development leads to birth defects that significantly affect quality of life. Studies during the past 40 years have advanced our understanding of the molecular mechanisms that establish facial form during development, highlighting the crucial roles in this process of a multipotent cell type known as the cranial neural crest cell. In this Review, we discuss recent advances in multi-omics and single-cell technologies that enable genes, transcriptional regulatory networks and epigenetic landscapes to be closely linked to the establishment of facial patterning and its variation, with an emphasis on normal and abnormal craniofacial morphogenesis. Advancing our knowledge of these processes will support important developments in tissue engineering, as well as the repair and reconstruction of the abnormal craniofacial complex.

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.

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.

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.

Dynamics of Chromatin Accessibility During Hematopoietic Stem Cell Differentiation Into Progressively Lineage-Committed Progeny

Epigenetic mechanisms regulate the multilineage differentiation capacity of hematopoietic stem cells (HSCs) into a variety of blood and immune cells. Mapping the chromatin dynamics of functionally defined cell populations will shed mechanistic insight into 2 major, unanswered questions in stem cell biology: how does epigenetic identity contribute to a cell type's lineage potential, and how do cascades of chromatin remodeling dictate ensuing fate decisions? Our recent work revealed evidence of multilineage gene priming in HSCs, where open cis-regulatory elements (CREs) exclusively shared between HSCs and unipotent lineage cells were enriched for DNA binding motifs of known lineage-specific transcription factors. Oligopotent progenitor populations operating between the HSCs and unipotent cells play essential roles in effecting hematopoietic homeostasis. To test the hypothesis that selective HSC-primed lineage-specific CREs remain accessible throughout differentiation, we used ATAC-seq to map the temporal dynamics of chromatin remodeling during progenitor differentiation. We observed epigenetic-driven clustering of oligopotent and unipotent progenitors into distinct erythromyeloid and lymphoid branches, with multipotent HSCs and MPPs associating with the erythromyeloid lineage. We mapped the dynamics of lineage-primed CREs throughout hematopoiesis and identified both unique and shared CREs as potential lineage reinforcement mechanisms at fate branch points. Additionally, quantification of genome-wide peak count and size revealed overall greater chromatin accessibility in HSCs, allowing us to identify HSC-unique peaks as putative regulators of self-renewal and multilineage potential. Finally, CRISPRi-mediated targeting of ATACseq-identified putative CREs in HSCs allowed us to demonstrate the functional role of selective CREs in lineage-specific gene expression. These findings provide insight into the regulation of stem cell multipotency and lineage commitment throughout hematopoiesis and serve as a resource to test functional drivers of hematopoietic lineage fate.

Functions and regulatory mechanisms of resting hematopoietic stem cells: a promising targeted therapeutic strategy

Hematopoietic stem cells (HSCs) are the common and essential precursors of all blood cells, including immune cells, and they are responsible for the lifelong maintenance and damage repair of blood tissue homeostasis. The vast majority (> 95%) of HSCs are in a resting state under physiological conditions and are only activated to play a functional role under stress conditions. This resting state affects their long-term survival and is also closely related to the lifelong maintenance of hematopoietic function; however, abnormal changes may also be an important factor leading to the decline of immune function in the body and the occurrence of diseases in various systems. While the importance of resting HSCs has attracted increasing research attention, our current understanding of this topic remains insufficient, and the direction of clinical targeted treatments is unclear. Here, we describe the functions of HSCs, analyze the regulatory mechanisms that affect their resting state, and discuss the relationship between resting HSCs and different diseases, with a view to providing guidance for the future clinical implementation of related targeted treatments.

Emergent Criticality in Coupled Boolean Networks

Early embryonic development involves forming all specialized cells from a fluid-like mass of identical stem cells. The differentiation process consists of a series of symmetry-breaking events, starting from a high-symmetry state (stem cells) to a low-symmetry state (specialized cells). This scenario closely resembles phase transitions in statistical mechanics. To theoretically study this hypothesis, we model embryonic stem cell (ESC) populations through a coupled Boolean network (BN) model. The interaction is applied using a multilayer Ising model that considers paracrine and autocrine signaling, along with external interventions. It is demonstrated that cell-to-cell variability can be interpreted as a mixture of steady-state probability distributions. Simulations have revealed that such models can undergo a series of first- and second-order phase transitions as a function of the system parameters that describe gene expression noise and interaction strengths. These phase transitions result in spontaneous symmetry-breaking events that generate new types of cells characterized by various steady-state distributions. Coupled BNs have also been shown to self-organize in states that allow spontaneous cell differentiation.

Molecular versatility during pluripotency progression

A challenge during development is to ensure lineage segregation while preserving plasticity. Using pluripotency progression as a paradigm, we review how developmental transitions are coordinated by redeployments, rather than global resettings, of cellular components. We highlight how changes in response to extrinsic cues (FGF, WNT, Activin/Nodal, Netrin-1), context- and stoichiometry-dependent action of transcription factors (Oct4, Nanog) and reconfigurations of epigenetic regulators (enhancers, promoters, TrxG, PRC) may confer robustness to naïve to primed pluripotency transition. We propose the notion of Molecular Versatility to regroup mechanisms by which molecules are repurposed to exert different, sometimes opposite, functions in close stem cell configurations.

Emergence of neuron types

Neuron types are the building blocks of the nervous system, and therefore, of functional circuits. Understanding the origin of neuronal diversity has always been an essential question in neuroscience and developmental biology. While knowledge on the molecular control of their diversification has largely increased during the last decades, it is now possible to reveal the dynamic mechanisms and the actual stepwise molecular changes occurring at single-cell level with the advent of single-cell omics technologies and analysis with high temporal resolution. Here, we focus on recent advances in the field and in technical and analytical tools that enable detailed insights into the emergence of neuron types in the central and peripheral nervous systems.

Early enforcement of cell identity by a functional component of the terminally differentiated state

How progenitor cells can attain a distinct differentiated cell identity is a challenging problem given the fluctuating signaling environment in which cells exist and that critical transcription factors are often not unique to a differentiation process. Here, we test the hypothesis that a unique differentiated cell identity can result from a core component of the differentiated state doubling up as a signaling protein that also drives differentiation. Using live single-cell imaging in the adipocyte differentiation system, we show that progenitor fat cells (preadipocytes) can only commit to terminally differentiate after up-regulating FABP4, a lipid buffer that is highly enriched in mature adipocytes. Upon induction of adipogenesis in mouse preadipocyte cells, we show that after a long delay, cells first abruptly start to engage a positive feedback between CEBPA and PPARG before then engaging, after a second delay, a positive feedback between FABP4 and PPARG. These sequential positive feedbacks both need to engage in order to drive PPARG levels past the threshold for irreversible differentiation. In the last step before commitment, PPARG transcriptionally increases FABP4 expression while fatty acid-loaded FABP4 increases PPARG activity. Together, our study suggests a control principle for robust cell identity whereby a core component of the differentiated state also promotes differentiation from its own progenitor state.

EBF1 primes B-lymphoid enhancers and limits the myeloid bias in murine multipotent progenitors

Hematopoietic stem cells (HSCs) and multipotent progenitors (MPPs) generate all cells of the blood system. Despite their multipotency, MPPs display poorly understood lineage bias. Here, we examine whether lineage-specifying transcription factors, such as the B-lineage determinant EBF1, regulate lineage preference in early progenitors. We detect low-level EBF1 expression in myeloid-biased MPP3 and lymphoid-biased MPP4 cells, coinciding with expression of the myeloid determinant C/EBPα. Hematopoietic deletion of Ebf1 results in enhanced myelopoiesis and reduced HSC repopulation capacity. Ebf1-deficient MPP3 and MPP4 cells exhibit an augmented myeloid differentiation potential and a transcriptome with an enriched C/EBPα signature. Correspondingly, EBF1 binds the Cebpa enhancer, and the deficiency and overexpression of Ebf1 in MPP3 and MPP4 cells lead to an up- and downregulation of Cebpa expression, respectively. In addition, EBF1 primes the chromatin of B-lymphoid enhancers specifically in MPP3 cells. Thus, our study implicates EBF1 in regulating myeloid/lymphoid fate bias in MPPs by constraining C/EBPα-driven myelopoiesis and priming the B-lymphoid fate.

Hematopoietic and Chronic Myeloid Leukemia Stem Cells: Multi-Stability versus Lineage Restriction

There is compelling evidence to support the view that the cell-of-origin for chronic myeloid leukemia is a hematopoietic stem cell. Unlike normal hematopoietic stem cells, the progeny of the leukemia stem cells are predominantly neutrophils during the disease chronic phase and there is a mild anemia. The hallmark oncogene for chronic myeloid leukemia is the BCR-ABLp210 fusion gene. Various studies have excluded a role for BCR-ABLp210 expression in maintaining the population of leukemia stem cells. Studies of BCR-ABLp210 expression in embryonal stem cells that were differentiated into hematopoietic stem cells and of the expression in transgenic mice have revealed that BCR-ABLp210 is able to veer hematopoietic stem and progenitor cells towards a myeloid fate. For the transgenic mice, global changes to the epigenetic landscape were observed. In chronic myeloid leukemia, the ability of the leukemia stem cells to choose from the many fates that are available to normal hematopoietic stem cells appears to be deregulated by BCR-ABLp210 and changes to the epigenome are also important. Even so, we still do not have a precise picture as to why neutrophils are abundantly produced in chronic myeloid leukemia.

Global genome decompaction leads to stochastic activation of gene expression as a first step toward fate commitment in human hematopoietic cells

When human cord blood-derived CD34+ cells are induced to differentiate, they undergo rapid and dynamic morphological and molecular transformations that are critical for fate commitment. In particular, the cells pass through a transitory phase known as "multilineage-primed" state. These cells are characterized by a mixed gene expression profile, different in each cell, with the coexpression of many genes characteristic for concurrent cell lineages. The aim of our study is to understand the mechanisms of the establishment and the exit from this transitory state. We investigated this issue using single-cell RNA sequencing and ATAC-seq. Two phases were detected. The first phase is a rapid and global chromatin decompaction that makes most of the gene promoters in the genome accessible for transcription. It results 24 h later in enhanced and pervasive transcription of the genome leading to the concomitant increase in the cell-to-cell variability of transcriptional profiles. The second phase is the exit from the multilineage-primed phase marked by a slow chromatin closure and a subsequent overall down-regulation of gene transcription. This process is selective and results in the emergence of coherent expression profiles corresponding to distinct cell subpopulations. The typical time scale of these events spans 48 to 72 h. These observations suggest that the nonspecificity of genome decompaction is the condition for the generation of a highly variable multilineage expression profile. The nonspecific phase is followed by specific regulatory actions that stabilize and maintain the activity of key genes, while the rest of the genome becomes repressed again by the chromatin recompaction. Thus, the initiation of differentiation is reminiscent of a constrained optimization process that associates the spontaneous generation of gene expression diversity to subsequent regulatory actions that maintain the activity of some genes, while the rest of the genome sinks back to the repressive closed chromatin state.

Zebrafish neuromesodermal progenitors undergo a critical state transition in vivo

The transition state model of cell differentiation proposes that a transient window of gene expression stochasticity precedes entry into a differentiated state. Here, we assess this theoretical model in zebrafish neuromesodermal progenitors (NMps) in vivo during late somitogenesis stages. We observed an increase in gene expression variability at the 24 somite stage (24ss) before their differentiation into spinal cord and paraxial mesoderm. Analysis of a published 18ss scRNA-seq dataset showed that the NMp population is noisier than its derivatives. By building in silico composite gene expression maps from image data, we assigned an 'NM index' to in silico NMps based on the expression of neural and mesodermal markers and demonstrated that cell population heterogeneity peaked at 24ss. Further examination revealed cells with gene expression profiles incongruent with their prospective fate. Taken together, our work supports the transition state model within an endogenous cell fate decision making event.

From developmental to atavistic bet-hedging: How cancer cells pervert the exploitation of random single-cell phenotypic fluctuations

Stochastic gene expression plays a leading developmental role through its contribution to cell differentiation. It is also proposed to promote phenotypic diversification in malignant cells. However, it remains unclear if these two forms of cellular bet-hedging are identical or rather display distinct features. Here we argue that bet-hedging phenomena in cancer cells are more similar to those occurring in unicellular organisms than to those of normal metazoan cells. We further propose that the atavistic bet-hedging strategies in cancer originate from a hijacking of the normal developmental bet-hedging of metazoans. Finally, we discuss the constraints that may shape the atavistic bet-hedging strategies of cancer cells.

Folate Carrier Deficiency Drives Differential Methylation and Enhanced Cellular Potency in the Neural Plate Border

The neural plate border (NPB) of vertebrate embryos segregates from the neural and epidermal regions, and it is comprised of an intermingled group of multipotent progenitor cells. Folate is the precursor of S-adenosylmethionine, the main methyl donor for DNA methylation, and it is critical for embryonic development, including the specification of progenitors which reside in the NPB. Despite the fact that several intersecting signals involved in the specification and territorial restriction of NPB cells are known, the role of epigenetics, particularly DNA methylation, has been a matter of debate. Here, we examined the temporal and spatial distribution of the methyl source and analyzed the abundance of 5mC/5 hmC and their epigenetic writers throughout the segregation of the neural and NPB territories. Reduced representation bisulfite sequencing (RRBS) on Reduced Folate Carrier 1 (RFC1)-deficient embryos leads to the identification of differentially methylated regions (DMRs). In the RFC1-deficient embryos, we identified several DMRs in the Notch1 locus, and the spatiotemporal expression of Notch1 and its downstream target gene Bmp4 were expanded in the NPB. Cell fate analysis on folate deficient embryos revealed a significant increase in the number of cells coexpressing both neural (SOX2) and NPB (PAX7) markers, which may represent an enhancing effect in the cellular potential of those progenitors. Taken together, our findings propose a model where the RFC1 deficiency drives methylation changes in specific genomic regions that are correlated with a dysregulation of pathways involved in early development such as Notch1 and BMP4 signaling. These changes affect the potency of the progenitors residing in the juncture of the neural plate and NPB territories, thus driving them to a primed state.

Cell and chromatin transitions in intestinal stem cell regeneration

The progeny of intestinal stem cells (ISCs) dedifferentiate in response to ISC attrition. The precise cell sources, transitional states, and chromatin remodeling behind this activity remain unclear. In the skin, stem cell recovery after injury preserves an epigenetic memory of the damage response; whether similar memories arise and persist in regenerated ISCs is not known. We addressed these questions by examining gene activity and open chromatin at the resolution of single Neurog3-labeled mouse intestinal crypt cells, hence deconstructing forward and reverse differentiation of the intestinal secretory (Sec) lineage. We show that goblet, Paneth, and enteroendocrine cells arise by multilineage priming in common precursors, followed by selective access at thousands of cell-restricted cis-elements. Selective ablation of the ISC compartment elicits speedy reversal of chromatin and transcriptional features in large fractions of precursor and mature crypt Sec cells without obligate cell cycle re-entry. ISC programs decay and reappear along a cellular continuum lacking discernible discrete interim states. In the absence of gross tissue damage, Sec cells simply reverse their forward trajectories, without invoking developmental or other extrinsic programs, and starting chromatin identities are effectively erased. These findings identify strikingly plastic molecular frameworks in assembly and regeneration of a self-renewing tissue.

Embryonic Origins of the Hematopoietic System: Hierarchies and Heterogeneity

The hierarchical framework of the adult blood system as we know it from current medical and hematology textbooks, displays a linear branching network of dividing and differentiated cells essential for the growth and maintenance of the healthy organism. This view of the hierarchy has evolved over the last 75 years. An amazing increase in cellular complexity has been realized; however, innovative single-cell technologies continue to uncover essential cell types and functions in animal models and the human blood system. The most potent cell of the hematopoietic hierarchy is the hematopoietic stem cell. Stem cells for adult tissues are the long-lived self-renewing cellular component, which ensure that differentiated tissue-specific cells are maintained and replaced through the entire adult lifespan. Although much blood research is focused on hematopoietic tissue homeostasis, replacement and regeneration during adult life, embryological studies have widened and enriched our understanding of additional developmental hierarchies and interacting cells of this life-sustaining tissue. Here, we review the current state of knowledge of the hierarchical organization and the vast heterogeneity of the hematopoietic system from embryonic to adult stages.

Globin vector regulatory elements are active in early hematopoietic progenitor cells

The globin genes are archetypal tissue-specific genes that are silent in most tissues but for late-stage erythroblasts upon terminal erythroid differentiation. The transcriptional activation of the β-globin gene is under the control of proximal and distal regulatory elements located on chromosome 11p15.4, including the β-globin locus control region (LCR). The incorporation of selected LCR elements in lentiviral vectors encoding β and β-like globin genes has enabled successful genetic treatment of the β-thalassemias and sickle cell disease. However, recent occurrences of benign clonal expansions in thalassemic patients and myelodysplastic syndrome in patients with sickle cell disease call attention to the non-erythroid functions of these powerful vectors. Here we demonstrate that lentivirally encoded LCR elements, in particular HS1 and HS2, can be activated in early hematopoietic cells including hematopoietic stem cells and myeloid progenitors. This activity is position-dependent and results in the transcriptional activation of a nearby reporter gene in these progenitor cell populations. We further show that flanking a globin vector with an insulator can effectively restrain this non-erythroid activity without impairing therapeutic globin expression. Globin lentiviral vectors harboring powerful LCR HS elements may thus expose to the risk of trans-activating cancer-related genes, which can be mitigated by a suitable insulator.

Detecting critical transition signals from single-cell transcriptomes to infer lineage-determining transcription factors

Analyzing single-cell transcriptomes promises to decipher the plasticity, heterogeneity, and rapid switches in developmental cellular state transitions. Such analyses require the identification of gene markers for semi-stable transition states. However, there are nontrivial challenges such as unexplainable stochasticity, variable population sizes, and alternative trajectory constructions. By advancing current tipping-point theory-based models with feature selection, network decomposition, accurate estimation of correlations, and optimization, we developed BioTIP to overcome these challenges. BioTIP identifies a small group of genes, called critical transition signal (CTS), to characterize regulated stochasticity during semi-stable transitions. Although methods rooted in different theories converged at the same transition events in two benchmark datasets, BioTIP is unique in inferring lineage-determining transcription factors governing critical transition. Applying BioTIP to mouse gastrulation data, we identify multiple CTSs from one dataset and validated their significance in another independent dataset. We detect the established regulator Etv2 whose expression change drives the haemato-endothelial bifurcation, and its targets together in CTS across three datasets. After comparing to three current methods using six datasets, we show that BioTIP is accurate, user-friendly, independent of pseudo-temporal trajectory, and captures significantly interconnected and reproducible CTSs. We expect BioTIP to provide great insight into dynamic regulations of lineage-determining factors.

FOXO Dictates Initiation of B Cell Development and Myeloid Restriction in Common Lymphoid Progenitors

The development of B cells relies on an intricate network of transcription factors critical for developmental progression and lineage commitment. In the B cell developmental trajectory, a temporal switch from predominant Foxo3 to Foxo1 expression occurs at the CLP stage. Utilizing VAV-iCre mediated conditional deletion, we found that the loss of FOXO3 impaired B cell development from LMPP down to B cell precursors, while the loss of FOXO1 impaired B cell commitment and resulted in a complete developmental block at the CD25 negative proB cell stage. Strikingly, the combined loss of FOXO1 and FOXO3 resulted in the failure to restrict the myeloid potential of CLPs and the complete loss of the B cell lineage. This is underpinned by the failure to enforce the early B-lineage gene regulatory circuitry upon a predominantly pre-established open chromatin landscape. Altogether, this demonstrates that FOXO3 and FOXO1 cooperatively govern early lineage restriction and initiation of B-lineage commitment in CLPs.

Differentiation Epitopes Define Hematopoietic Stem Cells and Change with Cell Cycle Passage

Hematopoietic stem cells express differentiation markers B220 and Gr1 and are proliferative. We have shown that the expression of these entities changes with cell cycle passage. Overall, we conclude that primitive hematopoietic stem cells alter their differentiation potential with cell cycle progression. Murine derived long-term hematopoietic stem cells (LT-HSC) are cycling and thus always changing phenotype. Here we show that over one half of marrow LT-HSC are in the population expressing differentiation epitopes and that B220 and Gr-1 positive populations are replete with LT-HSC after a single FACS separation but if subjected to a second separation these cells no longer contain LT-HSC. However, with second separated cells there is a population appearing that is B220 negative and replete with cycling c-Kit, Sca-1 CD150 positive LT-HSC. There is a 3-4 h interval between the first and second B220 or GR-1 FACS separation during which the stem cells continue to cycle. Thus, the LT-HSC have lost B220 or GR-1 expression as the cells progress through cell cycle, although they have maintained the c-kit, Sca-1 and CD150 stem cells markers over this time interval. These data indicate that cycling stem cells express differentiation epitopes and alter their differentiation potential with cell cycle passage.

Single-cell Transcriptomic Analysis Reveals the Cellular Heterogeneity of Mesenchymal Stem Cells

Ex vivo-expanded mesenchymal stem cells (MSCs) have been demonstrated to be a heterogeneous mixture of cells exhibiting varying proliferative, multipotential, and immunomodulatory capacities. However, the exact characteristics of MSCs remain largely unknown. By single-cell RNA sequencing of 61,296 MSCs derived from bone marrow and Wharton’s jelly, we revealed five distinct subpopulations. The developmental trajectory of these five MSC subpopulations was mapped, revealing a differentiation path from stem-like active proliferative cells (APCs) to multipotent progenitor cells, followed by branching into two paths: 1) unipotent preadipocytes or 2) bipotent prechondro-osteoblasts that were subsequently differentiated into unipotent prechondrocytes. The stem-like APCs, expressing the perivascular mesodermal progenitor markers CSPG4/MCAM/NES, uniquely exhibited strong proliferation and stemness signatures. Remarkably, the prechondrocyte subpopulation specifically expressed immunomodulatory genes and was able to suppress activated CD3⁺ T cell proliferation in vitro, supporting the role of this population in immunoregulation. In summary, our analysis mapped the heterogeneous subpopulations of MSCs and identified two subpopulations with potential functions in self-renewal and immunoregulation. Our findings advance the definition of MSCs by identifying the specific functions of their heterogeneous cellular composition, allowing for more specific and effective MSC application through the purification of their functional subpopulations.

Data-driven modeling predicts gene regulatory network dynamics during the differentiation of multipotential hematopoietic progenitors

Cellular differentiation during hematopoiesis is guided by gene regulatory networks (GRNs) comprising transcription factors (TFs) and the effectors of cytokine signaling. Based largely on analyses conducted at steady state, these GRNs are thought to be organized as a hierarchy of bistable switches, with antagonism between Gata1 and PU.1 driving red- and white-blood cell differentiation. Here, we utilize transient gene expression patterns to infer the genetic architecture—the type and strength of regulatory interconnections—and dynamics of a twelve-gene GRN including key TFs and cytokine receptors. We trained gene circuits, dynamical models that learn genetic architecture, on high temporal-resolution gene-expression data from the differentiation of an inducible cell line into erythrocytes and neutrophils. The model is able to predict the consequences of gene knockout, knockdown, and overexpression experiments and the inferred interconnections are largely consistent with prior empirical evidence. The inferred genetic architecture is densely interconnected rather than hierarchical, featuring extensive cross-antagonism between genes from alternative lineages and positive feedback from cytokine receptors. The analysis of the dynamics of gene regulation in the model reveals that PU.1 is one of the last genes to be upregulated in neutrophil conditions and that the upregulation of PU.1 and other neutrophil genes is driven by Cebpa and Gfi1 instead. This model inference is confirmed in an independent single-cell RNA-Seq dataset from mouse bone marrow in which Cebpa and Gfi1 expression precedes the neutrophil-specific upregulation of PU.1 during differentiation. These results demonstrate that full PU.1 upregulation during neutrophil development involves regulatory influences extrinsic to the Gata1-PU.1 bistable switch. Furthermore, although there is extensive cross-antagonism between erythroid and neutrophil genes, it does not have a hierarchical structure. More generally, we show that the combination of high-resolution time series data and data-driven dynamical modeling can uncover the dynamics and causality of developmental events that might otherwise be obscured.

The supply of blood cells is replenished by the maturation of hematopoietic progenitor cells into different cell types. Which cell type a progenitor cell develops into is determined by a complex network of genes whose protein products directly or indirectly regulate each others’ expression and that of downstream genes characteristic of the cell type. We inferred the nature and causality of the regulatory connections in a 12-gene network known to affect the decision between erythrocyte and neutrophil cell fates using a predictive machine-learning approach. Our analysis showed that the overall architecture of the network is densely interconnected and not hierarchical. Furthermore, the model inferred that PU.1, considered a master regulator of all white-blood cell lineages, is upregulated during neutrophil development by two other proteins, Cebpa and Gfi1. We validated this prediction by showing that Cebpa and Gfi1 expression precedes that of PU.1 in single-cell gene expression data from mouse bone marrow. These results revise the architecture of the gene network and the causality of regulatory events guiding hematopoiesis. The results also show that combining machine learning approaches with time course data can help resolve causality during development.

A regulatory network of microRNAs confers lineage commitment during early developmental trajectories of B and T lymphocytes

The commitment of hematopoietic multipotent progenitors (MPPs) toward a particular lineage involves activation of cell type-specific genes and silencing of genes that promote alternate cell fates. Although the gene expression programs of early-B and early-T lymphocyte development are mutually exclusive, we show that these cell types exhibit significantly correlated microRNA (miRNA) profiles. However, their corresponding miRNA targetomes are distinct and predominated by transcripts associated with natural killer, dendritic cell, and myeloid lineages, suggesting that miRNAs function in a cell-autonomous manner. The combinatorial expression of miRNAs miR-186-5p, miR-128-3p, and miR-330-5p in MPPs significantly attenuates their myeloid differentiation potential due to repression of myeloid-associated transcripts. Depletion of these miRNAs caused a pronounced de-repression of myeloid lineage targets in differentiating early-B and early-T cells, resulting in a mixed-lineage gene expression pattern. De novo motif analysis combined with an assay of promoter activities indicates that B as well as T lineage determinants drive the expression of these miRNAs in lymphoid lineages. Collectively, we present a paradigm that miRNAs are conserved between developing B and T lymphocytes, yet they target distinct sets of promiscuously expressed lineage-inappropriate genes to suppress the alternate cell-fate options. Thus, our studies provide a comprehensive compendium of miRNAs with functional implications for B and T lymphocyte development.

Transcriptome Dynamics of Epidermal Reprogramming during Direct Shoot Regeneration in Torenia fournieri

Shoot regeneration involves reprogramming of somatic cells and de novo organization of shoot apical meristems (SAMs). In the best-studied model system of shoot regeneration using Arabidopsis, regeneration is mediated by the auxin-responsive pluripotent callus formation from pericycle or pericycle-like tissues according to the lateral root development pathway. In contrast, shoot regeneration can be induced directly from fully differentiated epidermal cells of stem explants of Torenia fournieri (Torenia), without intervening the callus mass formation in culture with cytokinin; yet, its molecular mechanisms remain unaddressed. Here, we characterized this direct shoot regeneration by cytological observation and transcriptome analyses. The results showed that the gene expression profile rapidly changes upon culture to acquire a mixed signature of multiple organs/tissues, possibly associated with epidermal reprogramming. Comparison of transcriptomes between three different callus-inducing cultures (callus induction by auxin, callus induction by wounding and protoplast culture) of Arabidopsis and the Torenia stem culture identified genes upregulated in all the four culture systems as candidates of common factors of cell reprogramming. These initial changes proceeded independently of cytokinin, followed by cytokinin-dependent, transcriptional activations of nucleolar development and cell cycle. Later, SAM regulatory genes became highly expressed, leading to SAM organization in the foci of proliferating cells in the epidermal layer. Our findings revealed three distinct phases with different transcriptomic and regulatory features during direct shoot regeneration from the epidermis in Torenia, which provides a basis for further investigation of shoot regeneration in this unique culture system.

The cellular and molecular basis of somatosensory neuron development

Primary somatosensory neurons convey salient information about our external environment and internal state to the CNS, allowing us to detect, perceive, and react to a wide range of innocuous and noxious stimuli. Pseudo-unipolar in shape, and among the largest (longest) cells of most mammals, dorsal root ganglia (DRG) somatosensory neurons have peripheral axons that extend into skin, muscle, viscera, or bone and central axons that innervate the spinal cord and brainstem, where they synaptically engage the central somatosensory circuitry. Here, we review the diversity of mammalian DRG neuron subtypes and the intrinsic and extrinsic mechanisms that control their development. We describe classical and contemporary advances that frame our understanding of DRG neurogenesis, transcriptional specification of DRG neurons, and the establishment of morphological, physiological, and synaptic diversification across somatosensory neuron subtypes.

Intra-heterogeneity in transcription and chemoresistant property of leukemia-initiating cells in murine Setd2-/- acute myeloid leukemia

BACKGROUND

Heterogeneity of leukemia-initiating cells (LICs) is a major obstacle in acute myeloid leukemia (AML) therapy. Accumulated evidence indicates that the coexistence of multiple types of LICs with different pathogenicity in the same individual is a common feature in AML. However, the functional heterogeneity including the drug response of coexistent LICs remains unclear. Therefore, this study aimed to clarify the intra-heterogeneity in LICs that can help predict leukemia behavior and develop more effective treatments.

METHODS

Spleen cells from the primary Setd2^-/- -AML mouse were transplanted into C57BL/6 recipient mice to generate a transplantable model. Flow cytometry was used to analyze the immunophenotype of the leukemic mice. Whole-genome sequencing was conducted to detect secondary hits responsible for leukemia transformation. A serial transplantation assay was used to determine the self-renewal potential of Setd2^-/- -AML cells. A limiting-dilution assay was performed to identify the LIC frequency in different subsets of leukemia cells. Bulk and single-cell RNA sequencing were performed to analyze the transcriptional heterogeneity of LICs. Small molecular inhibitor screening and in vivo drug treatment were employed to clarify the difference in drug response between the different subsets of LICs.

RESULTS

In this study, we observed an aged Setd2^-/- mouse developing AML with co-mutation of Nras^G12S and Braf^K520E . Further investigation identified two types of LICs residing in the c-Kit⁺ B220⁺ Mac-1^- and c-Kit⁺ B220⁺ Mac-1⁺ subsets, respectively. In vivo transplantation assay disclosed the heterogeneity in differentiation between the coexistent LICs. Besides, an intrinsic doxorubicin-resistant transcriptional signature was uncovered in c-Kit⁺ B220⁺ Mac-1⁺ cells. Indeed, doxorubicin plus cytarabine (DA), the standard chemotherapeutic regimen used in AML treatment, could specifically kill c-Kit⁺ B220⁺ Mac-1^- cells, but it hardly affected c-Kit⁺ B220⁺ Mac-1⁺ cells. Transcriptome analysis unveiled a higher activation of RAS downstream signaling pathways in c-Kit⁺ B220⁺ Mac-1⁺ cells than in c-Kit⁺ B220⁺ Mac-1^- cells. Combined treatment with DA and RAS pathway inhibitors killed both c-Kit⁺ B220⁺ Mac-1^- and c-Kit⁺ B220⁺ Mac-1⁺ cells and attenuated disease progression.

CONCLUSIONS

This study identified two cell subsets enriched for LICs in murine Setd2^-/- -AML and disclosed the transcriptional and functional heterogeneity of LICs, revealing that the coexistence of different types of LICs in this model brings about diverse drug response.

B Lymphocyte Specification Is Preceded by Extensive Epigenetic Priming in Multipotent Progenitors

B lymphocyte development is dependent on the interplay between the chromatin landscape and lineage-specific transcription factors. It has been suggested that B lineage commitment is associated with major changes in the nuclear chromatin environment, proposing a critical role for lineage-specific transcription factors in the formation of the epigenetic landscape. In this report, we have used chromosome conformation capture in combination with assay for transposase-accessible chromatin sequencing analysis to enable highly efficient annotation of both proximal and distal transcriptional control elements to genes activated in B lineage specification in mice. A large majority of these genes were annotated to at least one regulatory element with an accessible chromatin configuration in multipotent progenitors. Furthermore, the majority of binding sites for the key regulators of B lineage specification, EBF1 and PAX5, occurred in already accessible regions. EBF1 did, however, cause a dynamic change in assay for transposase-accessible chromatin accessibility and was critical for an increase in distal promoter-enhancer interactions. Our data unravel an extensive epigenetic priming at regulatory elements annotated to lineage-restricted genes and provide insight into the interplay between the epigenetic landscape and transcription factors in cell specification.

Hematopoietic Stem Cells: Nature and Niche Nurture

Like all cells, hematopoietic stem cells (HSCs) and their offspring, the hematopoietic progenitor cells (HPCs), are highly sociable. Their capacity to interact with bone marrow niche cells and respond to environmental cytokines orchestrates the generation of the different types of blood and immune cells. The starting point for engineering hematopoiesis ex vivo is the nature of HSCs, and a longstanding premise is that they are a homogeneous population of cells. However, recent findings have shown that adult bone marrow HSCs are really a mixture of cells, with many having lineage affiliations. A second key consideration is: Do HSCs "choose" a lineage in a random and cell-intrinsic manner, or are they instructed by cytokines? Since their discovery, the hematopoietic cytokines have been viewed as survival and proliferation factors for lineage committed HPCs. Some are now known to also instruct cell lineage choice. These fundamental changes to our understanding of hematopoiesis are important for placing niche support in the right context and for fabricating an ex vivo environment to support HSC development.

Integrative Single-Cell RNA-Seq and ATAC-Seq Analysis of Human Developmental Hematopoiesis

Regulation of hematopoiesis during human development remains poorly defined. Here we applied single-cell RNA sequencing (scRNA-seq) and single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq) to over 8,000 human immunophenotypic blood cells from fetal liver and bone marrow. We inferred their differentiation trajectory and identified three highly proliferative oligopotent progenitor populations downstream of hematopoietic stem cells (HSCs)/multipotent progenitors (MPPs). Along this trajectory, we observed opposing patterns of chromatin accessibility and differentiation that coincided with dynamic changes in the activity of distinct lineage-specific transcription factors. Integrative analysis of chromatin accessibility and gene expression revealed extensive epigenetic but not transcriptional priming of HSCs/MPPs prior to their lineage commitment. Finally, we refined and functionally validated the sorting strategy for the HSCs/MPPs and achieved around 90% enrichment. Our study provides a useful framework for future investigation of human developmental hematopoiesis in the context of blood pathologies and regenerative medicine.

Computational Stem Cell Biology: Open Questions and Guiding Principles

Computational biology is enabling an explosive growth in our understanding of stem cells and our ability to use them for disease modeling, regenerative medicine, and drug discovery. We discuss four topics that exemplify applications of computation to stem cell biology: cell typing, lineage tracing, trajectory inference, and regulatory networks. We use these examples to articulate principles that have guided computational biology broadly and call for renewed attention to these principles as computation becomes increasingly important in stem cell biology. We also discuss important challenges for this field with the hope that it will inspire more to join this exciting area.

Chromatin accessibility maps provide evidence of multilineage gene priming in hematopoietic stem cells

Hematopoietic stem cells (HSCs) have the capacity to differentiate into vastly different types of mature blood cells. The epigenetic mechanisms regulating the multilineage ability, or multipotency, of HSCs are not well understood. To test the hypothesis that cis-regulatory elements that control fate decisions for all lineages are primed in HSCs, we used ATAC-seq to compare chromatin accessibility of HSCs with five unipotent cell types. We observed the highest similarity in accessibility profiles between megakaryocyte progenitors and HSCs, whereas B cells had the greatest number of regions with de novo gain in accessibility during differentiation. Despite these differences, we identified cis-regulatory elements from all lineages that displayed epigenetic priming in HSCs. These findings provide new insights into the regulation of stem cell multipotency, as well as a resource to identify functional drivers of lineage fate.

Epigenetic Signatures and Plasticity of Intestinal and Other Stem Cells

The cardinal properties of adult tissue stem cells are self-renewal and the ability to generate diverse resident cell types. The daily losses of terminally differentiated intestinal, skin, and blood cells require "professional" stem cells to produce replacements. This occurs by continuous expansion of stem cells and their immediate progeny, followed by coordinated activation of divergent transcriptional programs to generate stable cells with diverse functions. Other tissues turn over slowly, if at all, and vary widely in strategies for facultative stem cell activity or interconversion among mature resident cell types (transdifferentiation). Cell fate potential is programmed in tissue-specific configurations of chromatin, which restrict the complement of available genes and cis-regulatory elements, hence allowing specific cell types to arise. Using as a model the transcriptional and chromatin basis of cell differentiation and dedifferentiation in intestinal crypts, we discuss here how self-renewing and other tissues execute homeostatic and injury-responsive stem cell activity.

Tp63-expressing adult epithelial stem cells cross lineages boundaries revealing latent hairy skin competence

The formation of hair follicles, a landmark of mammals, requires complex mesenchymal–epithelial interactions and it is commonly believed that embryonic epidermal cells are the only cells that can respond to hair follicle morphogenetic signals in vivo. Here, we demonstrate that epithelial stem cells of non-skin origin (e.g. that of cornea, oesophagus, vagina, bladder, prostate) that express the transcription factor Tp63, a master gene for the development of epidermis and its appendages, can respond to skin morphogenetic signals. When exposed to a newborn skin microenvironment, these cells express hair-follicle lineage markers and contribute to hair follicles, sebaceous glands and/or epidermis renewal. Our results demonstrate that lineage restriction is not immutable and support the notion that all Tp63-expressing epithelial stem cells, independently of their embryonic origin, have latent skin competence explaining why aberrant hair follicles or sebaceous glands are sometimes observed in non-skin tissues (e.g. in cornea, vagina or thymus).

Adult stem cells are thought to be fate restricted to lineages distinct to their tissue of origin. Here, the authors demonstrate that Tp63 expressing epithelial stem cells from several disparate tissues can respond to skin morphogenetic signals and contribute to hair follicles, sebaceous glands and/or epidermis.

Dental cell type atlas reveals stem and differentiated cell types in mouse and human teeth

Understanding cell types and mechanisms of dental growth is essential for reconstruction and engineering of teeth. Therefore, we investigated cellular composition of growing and non-growing mouse and human teeth. As a result, we report an unappreciated cellular complexity of the continuously-growing mouse incisor, which suggests a coherent model of cell dynamics enabling unarrested growth. This model relies on spatially-restricted stem, progenitor and differentiated populations in the epithelial and mesenchymal compartments underlying the coordinated expansion of two major branches of pulpal cells and diverse epithelial subtypes. Further comparisons of human and mouse teeth yield both parallelisms and differences in tissue heterogeneity and highlight the specifics behind growing and non-growing modes. Despite being similar at a coarse level, mouse and human teeth reveal molecular differences and species-specific cell subtypes suggesting possible evolutionary divergence. Overall, here we provide an atlas of human and mouse teeth with a focus on growth and differentiation.

The transcriptional corepressor CBFA2T3 inhibits all-trans-retinoic acid-induced myeloid gene expression and differentiation in acute myeloid leukemia

CBFA2/RUNX1 partner transcriptional co-repressor 3 (CBFA2T3, also known as MTG16 or ETO2) is a myeloid translocation gene family protein that functions as a master transcriptional corepressor in hematopoiesis. Recently, it has been shown that CBFA2T3 maintains leukemia stem cell gene expression and promotes relapse in acute myeloid leukemia (AML). However, a role for CBFA2T3 in myeloid differentiation of AML has not been reported. Here, we show that CBFA2T3 represses retinoic acid receptor (RAR) target gene expression and inhibits all-trans-retinoic acid (ATRA)-induced myeloid differentiation of AML cells. ChIP-Seq revealed that CBFA2T3 targets the RARα/RXRα cistrome in U937 AML cells, predominantly at myeloid-specific enhancers associated with terminal differentiation. CRISPR/Cas9-mediated abrogation of CBFA2T3 resulted in spontaneous and ATRA-induced activation of myeloid-specific genes in a manner correlated with myeloid differentiation. Importantly, these effects were reversed by CBFA2T3 re-expression. Mechanistic studies showed that CBFA2T3 inhibits RAR target gene transcription by acting at an early step to regulate histone acetyltransferase recruitment, histone acetylation, and chromatin accessibility at RARα target sites, independently of the downstream, RAR-mediated steps of transcription. Finally, we validated the inhibitory effect of CBFA2T3 on RAR in multiple AML subtypes and patient samples. To our knowledge, this is the first study to show that CBFA2T3 down-regulation is both necessary and sufficient for enhancing ATRA-induced myeloid gene expression and differentiation of AML cells. Our findings suggest that CBFA2T3 can serve as a potential target for enhancing AML responsiveness to ATRA differentiation therapies.