Alternative platelet differentiation pathways initiated by nonhierarchically related hematopoietic stem cells.

Rare multipotent stem cells replenish millions of blood cells per second through a time-consuming process, passing through multiple stages of increasingly lineage-restricted progenitors. Although insults to the blood-forming system highlight the need for more rapid blood replenishment from stem cells, established models of hematopoiesis implicate only one mandatory differentiation pathway for each blood cell lineage. Here, we establish a nonhierarchical relationship between distinct stem cells that replenish all blood cell lineages and stem cells that replenish almost exclusively platelets, a lineage essential for hemostasis and with important roles in both the innate and adaptive immune systems. These distinct stem cells use cellularly, molecularly and functionally separate pathways for the replenishment of molecularly distinct megakaryocyte-restricted progenitors: a slower steady-state multipotent pathway and a fast-track emergency-activated platelet-restricted pathway. These findings provide a framework for enhancing platelet replenishment in settings in which slow recovery of platelets remains a major clinical challenge.

Single-Cell RNA-Seq Reveals Lineage and X Chromosome Dynamics in Human Preimplantation Embryos.

Mouse studies have been instrumental in forming our current understanding of early cell-lineage decisions; however, similar insights into the early human development are severely limited. Here, we present a comprehensive transcriptional map of human embryo development, including the sequenced transcriptomes of 1,529 individual cells from 88 human preimplantation embryos. These data show that cells undergo an intermediate state of co-expression of lineage-specific genes, followed by a concurrent establishment of the trophectoderm, epiblast, and primitive endoderm lineages, which coincide with blastocyst formation. Female cells of all three lineages achieve dosage compensation of X chromosome RNA levels prior to implantation. However, in contrast to the mouse, XIST is transcribed from both alleles throughout the progression of this expression dampening, and X chromosome genes maintain biallelic expression while dosage compensation proceeds. We envision broad utility of this transcriptional atlas in future studies on human development as well as in stem cell research.

The heterogeneity of human CD127(+) innate lymphoid cells revealed by single-cell RNA sequencing.

Innate lymphoid cells (ILCs) are increasingly appreciated as important participants in homeostasis and inflammation. Substantial plasticity and heterogeneity among ILC populations have been reported. Here we have delineated the heterogeneity of human ILCs through single-cell RNA sequencing of several hundreds of individual tonsil CD127(+) ILCs and natural killer (NK) cells. Unbiased transcriptional clustering revealed four distinct populations, corresponding to ILC1 cells, ILC2 cells, ILC3 cells and NK cells, with their respective transcriptomes recapitulating known as well as unknown transcriptional profiles. The single-cell resolution additionally divulged three transcriptionally and functionally diverse subpopulations of ILC3 cells. Our systematic comparison of single-cell transcriptional variation within and between ILC populations provides new insight into ILC biology during homeostasis, with additional implications for dysregulation of the immune system.

Initial seeding of the embryonic thymus by immune-restricted lympho-myeloid progenitors.

The final stages of restriction to the T cell lineage occur in the thymus after the entry of thymus-seeding progenitors (TSPs). The identity and lineage potential of TSPs remains unclear. Because the first embryonic TSPs enter a non-vascularized thymic rudiment, we were able to directly image and establish the functional and molecular properties of embryonic thymopoiesis-initiating progenitors (T-IPs) before their entry into the thymus and activation of Notch signaling. T-IPs did not include multipotent stem cells or molecular evidence of T cell-restricted progenitors. Instead, single-cell molecular and functional analysis demonstrated that most fetal T-IPs expressed genes of and had the potential to develop into lymphoid as well as myeloid components of the immune system. Moreover, studies of embryos deficient in the transcriptional regulator RBPJ demonstrated that canonical Notch signaling was not involved in pre-thymic restriction to the T cell lineage or the migration of T-IPs.

Phenotypic variation of transcriptomic cell types in mouse motor cortex.

Cortical neurons exhibit extreme diversity in gene expression as well as in morphological and electrophysiological properties1,2. Most existing neural taxonomies are based on either transcriptomic3,4 or morpho-electric5,6 criteria, as it has been technically challenging to study both aspects of neuronal diversity in the same set of cells7. Here we used Patch-seq8 to combine patch-clamp recording, biocytin staining, and single-cell RNA sequencing of more than 1,300 neurons in adult mouse primary motor cortex, providing a morpho-electric annotation of almost all transcriptomically defined neural cell types. We found that, although broad families of transcriptomic types (those expressing Vip, Pvalb, Sst and so on) had distinct and essentially non-overlapping morpho-electric phenotypes, individual transcriptomic types within the same family were not well separated in the morpho-electric space. Instead, there was a continuum of variability in morphology and electrophysiology, with neighbouring transcriptomic cell types showing similar morpho-electric features, often without clear boundaries between them. Our results suggest that neuronal types in the neocortex do not always form discrete entities. Instead, neurons form a hierarchy that consists of distinct non-overlapping branches at the level of families, but can form continuous and correlated transcriptomic and morpho-electrical landscapes within families.

Molecular spikes: a gold standard for single-cell RNA counting.

Single-cell sequencing methods rely on molecule-counting strategies to account for amplification biases, yet no experimental strategy to evaluate counting performance exists. Here, we introduce molecular spikes-RNA spike-ins containing built-in unique molecular identifiers (UMIs) that we use to identify critical experimental and computational conditions for accurate RNA counting in single-cell RNA-sequencing (scRNA-seq). Using molecular spikes, we uncovered impaired RNA counting in methods that were not informative for cellular RNA abundances due to inflated UMI counts. We further leverage molecular spikes to improve estimates of total endogenous RNA amounts in cells, and introduce a strategy to correct experiments with impaired RNA counting. The molecular spikes and the accompanying R package UMIcountR ( https://github.com/cziegenhain/UMIcountR ) will improve the validation of new methods, better estimate and adjust for cellular mRNA amounts and enable more indepth characterization of RNA counting in scRNA-seq.

Genomic encoding of transcriptional burst kinetics.

Mammalian gene expression is inherently stochastic1,2, and results in discrete bursts of RNA molecules that are synthesized from each allele3-7. Although transcription is known to be regulated by promoters and enhancers, it is unclear how cis-regulatory sequences encode transcriptional burst kinetics. Characterization of transcriptional bursting, including the burst size and frequency, has mainly relied on live-cell4,6,8 or single-molecule RNA fluorescence in situ hybridization3,5,8,9 recordings of selected loci. Here we determine transcriptome-wide burst frequencies and sizes for endogenous mouse and human genes using allele-sensitive single-cell RNA sequencing. We show that core promoter elements affect burst size and uncover synergistic effects between TATA and initiator elements, which were masked at mean expression levels. Notably, we provide transcriptome-wide evidence that enhancers control burst frequencies, and demonstrate that cell-type-specific gene expression is primarily shaped by changes in burst frequencies. Together, our data show that burst frequency is primarily encoded in enhancers and burst size in core promoters, and that allelic single-cell RNA sequencing is a powerful model for investigating transcriptional kinetics.

LCM-seq reveals unique transcriptional adaptation mechanisms of resistant neurons and identifies protective pathways in spinal muscular atrophy.

Somatic motor neurons are selectively vulnerable in spinal muscular atrophy (SMA), which is caused by a deficiency of the ubiquitously expressed survival of motor neuron protein. However, some motor neuron groups, including oculomotor and trochlear (ocular), which innervate eye muscles, are for unknown reasons spared. To reveal mechanisms of vulnerability and resistance in SMA, we investigate the transcriptional dynamics in discrete neuronal populations using laser capture microdissection coupled with RNA sequencing (LCM-seq). Using gene correlation network analysis, we reveal a TRP53-mediated stress response that is intrinsic to all somatic motor neurons independent of their vulnerability, but absent in relatively resistant red nucleus and visceral motor neurons. However, the temporal and spatial expression analysis across neuron types shows that the majority of SMA-induced modulations are cell type-specific. Using Gene Ontology and protein network analyses, we show that ocular motor neurons present unique disease-adaptation mechanisms that could explain their resilience. Specifically, ocular motor neurons up-regulate (1) Syt1, Syt5, and Cplx2, which modulate neurotransmitter release; (2) the neuronal survival factors Gdf15, Chl1, and Lif; (3) Aldh4, that protects cells from oxidative stress; and (4) the caspase inhibitor Pak4. Finally, we show that GDF15 can rescue vulnerable human spinal motor neurons from degeneration. This confirms that adaptation mechanisms identified in resilient neurons can be used to reduce susceptibility of vulnerable neurons. In conclusion, this in-depth longitudinal transcriptomics analysis in SMA reveals novel cell type-specific changes that, alone and combined, present compelling targets, including Gdf15, for future gene therapy studies aimed toward preserving vulnerable motor neurons.

Transcriptional bursts explain autosomal random monoallelic expression and affect allelic imbalance.

Transcriptional bursts render substantial biological noise in cellular transcriptomes. Here, we investigated the theoretical extent of allelic expression resulting from transcriptional bursting and how it compared to the amount biallelic, monoallelic and allele-biased expression observed in single-cell RNA-sequencing (scRNA-seq) data. We found that transcriptional bursting can explain the allelic expression patterns observed in single cells, including the frequent observations of autosomal monoallelic gene expression. Importantly, we identified that the burst frequency largely determined the fraction of cells with monoallelic expression, whereas the burst size had little effect on monoallelic observations. The high consistency between the bursting model predictions and scRNA-seq observations made it possible to assess the heterogeneity of a group of cells as their deviation in allelic observations from the expected. Finally, both burst frequency and size contributed to allelic imbalance observations and reinforced that studies of allelic imbalance can be confounded from the inherent noise in transcriptional bursting. Altogether, we demonstrate that allele-level transcriptional bursting renders widespread, although predictable, amounts of monoallelic and biallelic expression in single cells and cell populations.

A TRPV1-to-secretagogin regulatory axis controls pancreatic ß-cell survival by modulating protein turnover.

Ca2+-sensor proteins are generally implicated in insulin release through SNARE interactions. Here, secretagogin, whose expression in human pancreatic islets correlates with their insulin content and the incidence of type 2 diabetes, is shown to orchestrate an unexpectedly distinct mechanism. Single-cell RNA-seq reveals retained expression of the TRP family members in ß-cells from diabetic donors. Amongst these, pharmacological probing identifies Ca2+-permeable transient receptor potential vanilloid type 1 channels (TRPV1) as potent inducers of secretagogin expression through recruitment of Sp1 transcription factors. Accordingly, agonist stimulation of TRPV1s fails to rescue insulin release from pancreatic islets of glucose intolerant secretagogin knock-out(-/-) mice. However, instead of merely impinging on the SNARE machinery, reduced insulin availability in secretagogin-/- mice is due to ß-cell loss, which is underpinned by the collapse of protein folding and deregulation of secretagogin-dependent USP9X deubiquitinase activity. Therefore, and considering the desensitization of TRPV1s in diabetic pancreata, a TRPV1-to-secretagogin regulatory axis seems critical to maintain the structural integrity and signal competence of ß-cells.

Identification of spatial expression trends in single-cell gene expression data.

As methods for measuring spatial gene expression at single-cell resolution become available, there is a need for computational analysis strategies. We present trendsceek, a method based on marked point processes that identifies genes with statistically significant spatial expression trends. trendsceek finds these genes in spatial transcriptomic and sequential fluorescence in situ hybridization data, and also reveals significant gene expression gradients and hot spots in low-dimensional projections of dissociated single-cell RNA-seq data.

Random monoallelic expression of autosomal genes: stochastic transcription and allele-level regulation.

Random monoallelic expression (RME) of genes represents a striking example of how stochastic molecular processes can result in cellular heterogeneity. Recent transcriptome-wide studies have revealed both mitotically stable and cell-to-cell dynamic forms of autosomal RME, with the latter presumably resulting from burst-like stochastic transcription. Here, we discuss the distinguishing features of these two forms of RME and revisit literature on their nature, pervasiveness and regulation. Finally, we explore how RME may contribute to phenotypic variation, including the incomplete penetrance and variable expressivity often seen in genetic disease.

Disruption of the Extracellular Matrix Progressively Impairs Central Nervous System Vascular Maturation Downstream of ß-Catenin Signaling.

Objective- The Wnt/ß-catenin pathway orchestrates development of the blood-brain barrier, but the downstream mechanisms involved at different developmental windows and in different central nervous system (CNS) tissues have remained elusive. Approach and Results- Here, we create a new mouse model allowing spatiotemporal investigations of Wnt/ß-catenin signaling by induced overexpression of Axin1, an inhibitor of ß-catenin signaling, specifically in endothelial cells ( Axin1 iEC- OE). AOE (Axin1 overexpression) in Axin1 iEC- OE mice at stages following the initial vascular invasion of the CNS did not impair angiogenesis but led to premature vascular regression followed by progressive dilation and inhibition of vascular maturation resulting in forebrain-specific hemorrhage 4 days post-AOE. Analysis of the temporal Wnt/ß-catenin driven CNS vascular development in zebrafish also suggested that Axin1 iEC- OE led to CNS vascular regression and impaired maturation but not inhibition of ongoing angiogenesis within the CNS. Transcriptomic profiling of isolated, ß-catenin signaling-deficient endothelial cells during early blood-brain barrier-development (E11.5) revealed ECM (extracellular matrix) proteins as one of the most severely deregulated clusters. Among the 20 genes constituting the forebrain endothelial cell-specific response signature, 8 ( Adamtsl2, Apod, Ctsw, Htra3, Pglyrp1, Spock2, Ttyh2, and Wfdc1) encoded bona fide ECM proteins. This specific ß-catenin-responsive ECM signature was also repressed in Axin1 iEC- OE and endothelial cell-specific ß-catenin-knockout mice ( Ctnnb1-KOiEC) during initial blood-brain barrier maturation (E14.5), consistent with an important role of Wnt/ß-catenin signaling in orchestrating the development of the forebrain vascular ECM. Conclusions- These results suggest a novel mechanism of establishing a CNS endothelium-specific ECM signature downstream of Wnt-ß-catenin that impact spatiotemporally on blood-brain barrier differentiation during forebrain vessel development. Visual Overview- An online visual overview is available for this article.

Intussusceptive Vascular Remodeling Precedes Pathological Neovascularization.

Objective- Pathological neovascularization is crucial for progression and morbidity of serious diseases such as cancer, diabetic retinopathy, and age-related macular degeneration. While mechanisms of ongoing pathological neovascularization have been extensively studied, the initiating pathological vascular remodeling (PVR) events, which precede neovascularization remains poorly understood. Here, we identify novel molecular and cellular mechanisms of preneovascular PVR, by using the adult choriocapillaris as a model. Approach and Results- Using hypoxia or forced overexpression of VEGF (vascular endothelial growth factor) in the subretinal space to induce PVR in zebrafish and rats respectively, and by analyzing choriocapillaris membranes adjacent to choroidal neovascular lesions from age-related macular degeneration patients, we show that the choriocapillaris undergo robust induction of vascular intussusception and permeability at preneovascular stages of PVR. This PVR response included endothelial cell proliferation, formation of endothelial luminal processes, extensive vesiculation and thickening of the endothelium, degradation of collagen fibers, and splitting of existing extravascular columns. RNA-sequencing established a role for endothelial tight junction disruption, cytoskeletal remodeling, vesicle- and cilium biogenesis in this process. Mechanistically, using genetic gain- and loss-of-function zebrafish models and analysis of primary human choriocapillaris endothelial cells, we determined that HIF (hypoxia-induced factor)-1α-VEGF-A-VEGFR2 signaling was important for hypoxia-induced PVR. Conclusions- Our findings reveal that PVR involving intussusception and splitting of extravascular columns, endothelial proliferation, vesiculation, fenestration, and thickening is induced before neovascularization, suggesting that identifying and targeting these processes may prevent development of advanced neovascular disease in the future. Visual Overview- An online visual overview is available for this article.

Mouse Model of Alagille Syndrome and Mechanisms of Jagged1 Missense Mutations.

BACKGROUND & AIMS: Alagille syndrome is a genetic disorder characterized by cholestasis, ocular abnormalities, characteristic facial features, heart defects, and vertebral malformations. Most cases are associated with mutations in JAGGED1 (JAG1), which encodes a Notch ligand, although it is not clear how these contribute to disease development. We aimed to develop a mouse model of Alagille syndrome to elucidate these mechanisms. METHODS: Mice with a missense mutation (H268Q) in Jag1 (Jag1+/Ndr mice) were outbred to a C3H/C57bl6 background to generate a mouse model for Alagille syndrome (Jag1Ndr/Ndr mice). Liver tissues were collected at different timepoints during development, analyzed by histology, and liver organoids were cultured and analyzed. We performed transcriptome analysis of Jag1Ndr/Ndr livers and livers from patients with Alagille syndrome, cross-referenced to the Human Protein Atlas, to identify commonly dysregulated pathways and biliary markers. We used species-specific transcriptome separation and ligand-receptor interaction assays to measure Notch signaling and the ability of JAG1Ndr to bind or activate Notch receptors. We studied signaling of JAG1 and JAG1Ndr via NOTCH 1, NOTCH2, and NOTCH3 and resulting gene expression patterns in parental and NOTCH1-expressing C2C12 cell lines. RESULTS: Jag1Ndr/Ndr mice had many features of Alagille syndrome, including eye, heart, and liver defects. Bile duct differentiation, morphogenesis, and function were dysregulated in newborn Jag1Ndr/Ndr mice, with aberrations in cholangiocyte polarity, but these defects improved in adult mice. Jag1Ndr/Ndr liver organoids collapsed in culture, indicating structural instability. Whole-transcriptome sequence analyses of liver tissues from mice and patients with Alagille syndrome identified dysregulated genes encoding proteins enriched at the apical side of cholangiocytes, including CFTR and SLC5A1, as well as reduced expression of IGF1. Exposure of Notch-expressing cells to JAG1Ndr, compared with JAG1, led to hypomorphic Notch signaling, based on transcriptome analysis. JAG1-expressing cells, but not JAG1Ndr-expressing cells, bound soluble Notch1 extracellular domain, quantified by flow cytometry. However, JAG1 and JAG1Ndr cells each bound NOTCH2, and signaling from NOTCH2 signaling was reduced but not completely inhibited, in response to JAG1Ndr compared with JAG1. CONCLUSIONS: In mice, expression of a missense mutant of Jag1 (Jag1Ndr) disrupts bile duct development and recapitulates Alagille syndrome phenotypes in heart, eye, and craniofacial dysmorphology. JAG1Ndr does not bind NOTCH1, but binds NOTCH2, and elicits hypomorphic signaling. This mouse model can be used to study other features of Alagille syndrome and organ development.

Single-cell analyses of X Chromosome inactivation dynamics and pluripotency during differentiation.

Pluripotency, differentiation, and X Chromosome inactivation (XCI) are key aspects of embryonic development. However, the underlying relationship and mechanisms among these processes remain unclear. Here, we systematically dissected these features along developmental progression using mouse embryonic stem cells (mESCs) and single-cell RNA sequencing with allelic resolution. We found that mESCs grown in a ground state 2i condition displayed transcriptomic profiles diffused from preimplantation mouse embryonic cells, whereas EpiStem cells closely resembled the post-implantation epiblast. Sex-related gene expression varied greatly across distinct developmental states. We also identified novel markers that were highly enriched in each developmental state. Moreover, we revealed that several novel pathways, including PluriNetWork and Focal Adhesion, were responsible for the delayed progression of female EpiStem cells. Importantly, we "digitalized" XCI progression using allelic expression of active and inactive X Chromosomes and surprisingly found that XCI states exhibited profound variability in each developmental state, including the 2i condition. XCI progression was not tightly synchronized with loss of pluripotency and increase of differentiation at the single-cell level, although these processes were globally correlated. In addition, highly expressed genes, including core pluripotency factors, were in general biallelically expressed. Taken together, our study sheds light on the dynamics of XCI progression and the asynchronicity between pluripotency, differentiation, and XCI.

Cell shape determines gene expression: cardiomyocyte morphotypic transcriptomes.

Cardiomyocytes undergo considerable changes in cell shape. These can be due to hemodynamic constraints, including changes in preload and afterload conditions, or to mutations in genes important for cardiac function. These changes instigate significant changes in cellular architecture and lead to the addition of sarcomeres, at the same time or at a later stage. However, it is currently unknown whether changes in cell shape on their own affect gene expression and the aim of this study was to fill that gap in our knowledge. We developed a single-cell morphotyping strategy, followed by single-cell RNA sequencing, to determine the effects of altered cell shape in gene expression. This enabled us to profile the transcriptomes of individual cardiomyocytes of defined geometrical morphotypes and characterize them as either normal or pathological conditions. We observed that deviations from normal cell shapes were associated with significant downregulation of gene expression and deactivation of specific pathways, like oxidative phosphorylation, protein kinase A, and cardiac beta-adrenergic signaling pathways. In addition, we observed that genes involved in apoptosis of cardiomyocytes and necrosis were upregulated in square-like pathological shapes. Mechano-sensory pathways, including integrin and Src kinase mediated signaling, appear to be involved in the regulation of shape-dependent gene expression. Our study demonstrates that cell shape per se affects the regulation of the transcriptome in cardiac myocytes, an effect with possible implications for cardiovascular disease.

Exploring parasite heterogeneity using single-cell RNA-seq reveals a gene signature among sexual stage Plasmodium falciparum parasites.

The malaria parasite has a complex lifecycle, including several events of differentiation and stage progression, while actively evading immunity in both its mosquito and human hosts. Important parasite gene expression and regulation during these events remain hidden in rare populations of cells. Here, we combine a capillary-based platform for cell isolation with single-cell RNA-sequencing to transcriptionally profile 165 single infected red blood cells (iRBCs) during the intra-erythrocytic developmental cycle (IDC). Unbiased analyses of single-cell data grouped the cells into eight transcriptional states during IDC. Interestingly, we uncovered a gene signature from the single iRBC analyses that can successfully discriminate between developing asexual and sexual stage parasites at cellular resolution, and we verify five, previously undefined, gametocyte stage specific genes. Moreover, we show the capacity of detecting expressed genes from the variable gene families in single parasites, despite the sparse nature of data. In total, the single parasite transcriptomics holds promise for molecular dissection of rare parasite phenotypes throughout the malaria lifecycle.

Sequentially acting Sox transcription factors in neural lineage development.

Pluripotent embryonic stem (ES) cells can generate all cell types, but how cell lineages are initially specified and maintained during development remains largely unknown. Different classes of Sox transcription factors are expressed during neurogenesis and have been assigned important roles from early lineage specification to neuronal differentiation. Here we characterize the genome-wide binding for Sox2, Sox3, and Sox11, which have vital functions in ES cells, neural precursor cells (NPCs), and maturing neurons, respectively. The data demonstrate that Sox factor binding depends on developmental stage-specific constraints and reveal a remarkable sequential binding of Sox proteins to a common set of neural genes. Interestingly, in ES cells, Sox2 preselects for neural lineage-specific genes destined to be bound and activated by Sox3 in NPCs. In NPCs, Sox3 binds genes that are later bound and activated by Sox11 in differentiating neurons. Genes prebound by Sox proteins are associated with a bivalent chromatin signature, which is resolved into a permissive monovalent state upon binding of activating Sox factors. These data indicate that a single key transcription factor family acts sequentially to coordinate neural gene expression from the early lineage specification in pluripotent cells to later stages of neuronal development.