Emergence of the erythroid lineage from multipotent hematopoiesis [preprint]

Red cell formation begins with the hematopoietic stem cell, but the manner by which it gives rise to erythroid progenitors, and their subsequent developmental path, remain unclear. Here we combined single-cell transcriptomics of murine hematopoietic tissues with fate potential assays to infer a continuous yet hierarchical structure for the hematopoietic network. We define the erythroid differentiation trajectory as it emerges from multipotency and diverges from 6 other blood lineages. With the aid of a new flow-cytometric sorting strategy, we validated predicted cell fate potentials at the single cell level, revealing a coupling between erythroid and basophil/mast cell fates. We uncovered novel growth factor receptor regulators of the erythroid trajectory, including the proinflammatory IL- 17RA, found to be a strong erythroid stimulator; and identified a global hematopoietic response to stress erythropoiesis. We further identified transcriptional and high-purity FACS gates for the complete isolation of all classically-defined erythroid burst-forming (BFU-e) and colony-forming progenitors (CFU-e), finding that they express a dedicated transcriptional program, distinct from that of terminally-differentiating erythroblasts. Intriguingly, profound remodeling of the cell cycle is intimately entwined with CFU-e developmental progression and with a sharp transcriptional switch that extinguishes the CFU-e stage and activates terminal differentiation. Underlying these results, our work showcases the utility of theoretic approaches linking transcriptomic data to predictive fate models, providing key insights into lineage development in vivo

A Large Fraction of Extragenic RNA Pol II Transcription Sites Overlap Enhancers

A substantial fraction of extragenic Pol II transcription sites coincides with transcriptional enhancers, which may be relevant for functional annotation of mammalian genomes

Fundamental limits on dynamic inference from single-cell snapshots

Single-cell expression profiling reveals the molecular states of individual cells with unprecedented detail. Because these methods destroy cells in the process of analysis, they cannot measure how gene expression changes over time. However, some information on dynamics is present in the data: the continuum of molecular states in the population can reflect the trajectory of a typical cell. Many methods for extracting single-cell dynamics from population data have been proposed. However, all such attempts face a common limitation: for any measured distribution of cell states, there are multiple dynamics that could give rise to it, and by extension, multiple possibilities for underlying mechanisms of gene regulation. Here, we describe the aspects of gene expression dynamics that cannot be inferred from a static snapshot alone and identify assumptions necessary to constrain a unique solution for cell dynamics from static snapshots. We translate these constraints into a practical algorithmic approach, population balance analysis (PBA), which makes use of a method from spectral graph theory to solve a class of high-dimensional differential equations. We use simulations to show the strengths and limitations of PBA, and then apply it to single-cell profiles of hematopoietic progenitor cells (HPCs). Cell state predictions from this analysis agree with HPC fate assays reported in several papers over the past two decades. By highlighting the fundamental limits on dynamic inference faced by any method, our framework provides a rigorous basis for dynamic interpretation of a gene expression continuum and clarifies best experimental designs for trajectory reconstruction from static snapshot measurements