Identification of a non-canonical ciliate nuclear genetic code where UAA and UAG code for different amino acids

The genetic code is one of the most highly conserved features across life. Only a few lineages have deviated from the "universal" genetic code. Amongst the few variants of the genetic code reported to date, the codons UAA and UAG virtually always have the same translation, suggesting that their evolution is coupled. Here, we report the genome and transcriptome sequencing of a novel uncultured ciliate, belonging to the Oligohymenophorea class, where the translation of the UAA and UAG stop codons have changed to specify different amino acids. Genomic and transcriptomic analyses revealed that UAA has been reassigned to encode lysine, while UAG has been reassigned to encode glutamic acid. We identified multiple suppressor tRNA genes with anticodons complementary to the reassigned codons. We show that the retained UGA stop codon is enriched in the 3'UTR immediately downstream of the coding region of genes, suggesting that there is functional drive to maintain tandem stop codons. Using a phylogenomics approach, we reconstructed the ciliate phylogeny and mapped genetic code changes, highlighting the remarkable number of independent genetic code changes within the Ciliophora group of protists. According to our knowledge, this is the first report of a genetic code variant where UAA and UAG encode different amino acids.

Perivascular niche cells sense thrombocytopenia and activate hematopoietic stem cells in an IL-1 dependent manner

Hematopoietic stem cells (HSCs) residing in specialized niches in the bone marrow are responsible for the balanced output of multiple short-lived blood cell lineages in steady-state and in response to different challenges. However, feedback mechanisms by which HSCs, through their niches, sense acute losses of specific blood cell lineages remain to be established. While all HSCs replenish platelets, previous studies have shown that a large fraction of HSCs are molecularly primed for the megakaryocyte-platelet lineage and are rapidly recruited into proliferation upon platelet depletion. Platelets normally turnover in an activation-dependent manner, herein mimicked by antibodies inducing platelet activation and depletion. Antibody-mediated platelet activation upregulates expression of Interleukin-1 (IL-1) in platelets, and in bone marrow extracellular fluid in vivo. Genetic experiments demonstrate that rather than IL-1 directly activating HSCs, activation of bone marrow Lepr+ perivascular niche cells expressing IL-1 receptor is critical for the optimal activation of quiescent HSCs upon platelet activation and depletion. These findings identify a feedback mechanism by which activation-induced depletion of a mature blood cell lineage leads to a niche-dependent activation of HSCs to reinstate its homeostasis.

Deep sequencing of Phox2a nuclei reveals five classes of anterolateral system neurons

The anterolateral system (ALS) is a major ascending pathway from the spinal cord that projects to multiple brain areas and underlies the perception of pain, itch and skin temperature. Despite its importance, our understanding of this system has been hampered by the considerable functional and molecular diversity of its constituent cells. Here we use fluorescence-activated cell sorting to isolate ALS neurons belonging to the Phox2a-lineage for single-nucleus RNA sequencing. We reveal five distinct clusters of ALS neurons (ALS1-5) and document their laminar distribution in the spinal cord using in situ hybridization. We identify 3 clusters of neurons located predominantly in laminae I-III of the dorsal horn (ALS1-3) and two clusters with cell bodies located in deeper laminae (ALS4 & ALS5). Our findings reveal the transcriptional logic that underlies ALS neuronal diversity in the adult mouse and uncover the molecular identity of two previously identified classes of projection neurons. We also show that these molecular signatures can be used to target groups of ALS neurons using retrograde viral tracing. Overall, our findings provide a valuable resource for studying somatosensory biology and targeting subclasses of ALS neurons.

Panhematopoietic RNA barcoding enables kinetic measurements of nucleate and anucleate lineages and the activation of myeloid clones following acute platelet depletion

BACKGROUND

Platelets and erythrocytes constitute over 95% of all hematopoietic stem cell output. However, the clonal dynamics of HSC contribution to these lineages remains largely unexplored.

RESULTS

We use lentiviral genetic labeling of mouse hematopoietic stem cells to quantify output from all lineages, nucleate, and anucleate, simultaneously linking these with stem and progenitor cell transcriptomic phenotypes using single-cell RNA-sequencing. We observe dynamic shifts of clonal behaviors through time in same-animal peripheral blood and demonstrate that acute platelet depletion shifts the output of multipotent hematopoietic stem cells to the exclusive production of platelets. Additionally, we observe the emergence of new myeloid-biased clones, which support short- and long-term production of blood cells.

CONCLUSIONS

Our approach enables kinetic studies of multi-lineage output in the peripheral blood and transcriptional heterogeneity of individual hematopoietic stem cells. Our results give a unique insight into hematopoietic stem cell reactivation upon platelet depletion and of clonal dynamics in both steady state and under stress.

Single-cell gene and isoform expression analysis reveals signatures of ageing in haematopoietic stem and progenitor cells

Single-cell approaches have revealed that the haematopoietic hierarchy is a continuum of differentiation, from stem cell to committed progenitor, marked by changes in gene expression. However, many of these approaches neglect isoform-level information and thus do not capture the extent of alternative splicing within the system. Here, we present an integrated short- and long-read single-cell RNA-seq analysis of haematopoietic stem and progenitor cells. We demonstrate that over half of genes detected in standard short-read single-cell analyses are expressed as multiple, often functionally distinct, isoforms, including many transcription factors and key cytokine receptors. We observe global and HSC-specific changes in gene expression with ageing but limited impact of ageing on isoform usage. Integrating single-cell and cell-type-specific isoform landscape in haematopoiesis thus provides a new reference for comprehensive molecular profiling of heterogeneous tissues, as well as novel insights into transcriptional complexity, cell-type-specific splicing events and consequences of ageing.

Single-cell genomics reveals population structures from in vitro evolutionary studies of Salmonella

Single-cell DNA sequencing has the potential to reveal detailed hierarchical structures in evolving populations of cells. Single cell approaches are increasingly used to study clonal evolution in human ageing and cancer but have not yet been deployed to study evolving clonal microbial populations. Here, we present an approach for single bacterial genomic analysis for in vitro evolution experiments using FACS isolation of individual bacteria followed by whole-genome amplification and sequencing. We apply this to the experimental evolution of a hypermutator strain of Salmonella in response to antibiotic stress (ciprofloxacin). By analysing sequence polymorphisms in individual cells from populations we identified the presence and prevalence of sub-populations which have acquired polymorphisms in genes previously demonstrated to be associated with ciprofloxacin susceptibility. We were also able to identify that the population exposed to antibiotic stress was able to develop resistance whilst maintaining diversity. This population structure could not be resolved from bulk sequence data, and our results show how high-throughput single-cell sequencing can enhance experimental studies of bacterial evolution.

Into the multiverse: advances in single-cell multiomic profiling

Single-cell transcriptomic approaches have revolutionised the study of complex biological systems, with the routine measurement of gene expression in thousands of cells enabling construction of whole-organism cell atlases. However, the transcriptome is just one layer amongst many that coordinate to define cell type and state and, ultimately, function. In parallel with the widespread uptake of single-cell RNA-seq (scRNA-seq), there has been a rapid emergence of methods that enable multiomic profiling of individual cells, enabling parallel measurement of intercellular heterogeneity in the genome, epigenome, transcriptome, and proteomes. Linking measurements from each of these layers has the potential to reveal regulatory and functional mechanisms underlying cell behaviour in healthy development and disease.

Transcriptional characterization of human megakaryocyte polyploidization and lineage commitment

BACKGROUND

Megakaryocytes (MKs) originate from cells immuno-phenotypically indistinguishable from hematopoietic stem cells (HSCs), bypassing intermediate progenitors. They mature within the adult bone marrow and release platelets into the circulation. Until now, there have been no transcriptional studies of primary human bone marrow MKs.

OBJECTIVES

To characterize MKs and HSCs from human bone marrow using single-cell RNA sequencing, to investigate MK lineage commitment, maturation steps, and thrombopoiesis.

RESULTS

We show that MKs at different levels of polyploidization exhibit distinct transcriptional states. Although high levels of platelet-specific gene expression occur in the lower ploidy classes, as polyploidization increases, gene expression is redirected toward translation and posttranslational processing transcriptional programs, in preparation for thrombopoiesis. Our findings are in keeping with studies of MK ultrastructure and supersede evidence generated using in vitro cultured MKs. Additionally, by analyzing transcriptional signatures of a single HSC, we identify two MK-biased HSC subpopulations exhibiting unique differentiation kinetics. We show that human bone marrow MKs originate from these HSC subpopulations, supporting the notion that they display priming for MK differentiation. Finally, to investigate transcriptional changes in MKs associated with stress thrombopoiesis, we analyzed bone marrow MKs from individuals with recent myocardial infarction and found a specific gene expression signature. Our data support the modulation of MK differentiation in this thrombotic state.

CONCLUSIONS

Here, we use single-cell sequencing for the first time to characterize the human bone marrow MK transcriptome at different levels of polyploidization and investigate their differentiation from the HSC.

Non-neuronal expression of SARS-CoV-2 entry genes in the olfactory system suggests mechanisms underlying COVID-19-associated anosmia

Altered olfactory function is a common symptom of COVID-19, but its etiology is unknown. A key question is whether SARS-CoV-2 (CoV-2) - the causal agent in COVID-19 - affects olfaction directly, by infecting olfactory sensory neurons or their targets in the olfactory bulb, or indirectly, through perturbation of supporting cells. Here we identify cell types in the olfactory epithelium and olfactory bulb that express SARS-CoV-2 cell entry molecules. Bulk sequencing demonstrated that mouse, non-human primate and human olfactory mucosa expresses two key genes involved in CoV-2 entry, ACE2 and TMPRSS2. However, single cell sequencing revealed that ACE2 is expressed in support cells, stem cells, and perivascular cells, rather than in neurons. Immunostaining confirmed these results and revealed pervasive expression of ACE2 protein in dorsally-located olfactory epithelial sustentacular cells and olfactory bulb pericytes in the mouse. These findings suggest that CoV-2 infection of non-neuronal cell types leads to anosmia and related disturbances in odor perception in COVID-19 patients.

Defining Cell Identity with Single-Cell Omics

Cells are a fundamental unit of life, and the ability to study the phenotypes and behaviors of individual cells is crucial to understanding the workings of complex biological systems. Cell phenotypes (epigenomic, transcriptomic, proteomic, and metabolomic) exhibit dramatic heterogeneity between and within the different cell types and states underlying cellular functional diversity. Cell genotypes can also display heterogeneity throughout an organism, in the form of somatic genetic variation-most notably in the emergence and evolution of tumors. Recent technical advances in single-cell isolation and the development of omics approaches sensitive enough to reveal these aspects of cell identity have enabled a revolution in the study of multicellular systems. In this review, we discuss the technologies available to resolve the genomes, epigenomes, transcriptomes, proteomes, and metabolomes of single cells from a wide variety of living systems.

Pluripotent state transitions coordinate morphogenesis in mouse and human embryos

The foundations of mammalian development lie in a cluster of embryonic epiblast stem cells. In response to extracellular matrix signalling, these cells undergo epithelialization and create an apical surface in contact with a cavity, a fundamental event for all subsequent development. Concomitantly, epiblast cells transit through distinct pluripotent states, before lineage commitment at gastrulation. These pluripotent states have been characterized at the molecular level, but their biological importance remains unclear. Here we show that exit from an unrestricted naive pluripotent state is required for epiblast epithelialization and generation of the pro-amniotic cavity in mouse embryos. Embryonic stem cells locked in the naive state are able to initiate polarization but fail to undergo lumenogenesis. Mechanistically, exit from naive pluripotency activates an Oct4-governed transcriptional program that results in expression of glycosylated sialomucin proteins and the vesicle tethering and fusion events of lumenogenesis. Similarly, exit of epiblasts from naive pluripotency in cultured human post-implantation embryos triggers amniotic cavity formation and developmental progression. Our results add tissue-level architecture as a new criterion for the characterization of different pluripotent states, and show the relevance of transitions between these states during development of the mammalian embryo.

Proliferation Drives Aging-Related Functional Decline in a Subpopulation of the Hematopoietic Stem Cell Compartment

Aging of the hematopoietic stem cell (HSC) compartment is characterized by lineage bias and reduced stem cell function, the molecular basis of which is largely unknown. Using single-cell transcriptomics, we identified a distinct subpopulation of old HSCs carrying a p53 signature indicative of stem cell decline alongside pro-proliferative JAK/STAT signaling. To investigate the relationship between JAK/STAT and p53 signaling, we challenged HSCs with a constitutively active form of JAK2 (V617F) and observed an expansion of the p53-positive subpopulation in old mice. Our results reveal cellular heterogeneity in the onset of HSC aging and implicate a role for JAK2V617F-driven proliferation in the p53-mediated functional decline of old HSCs.

Power analysis of single-cell RNA-sequencing experiments

Single-cell RNA sequencing (scRNA-seq) has become an established and powerful method to investigate transcriptomic cell-to-cell variation, thereby revealing new cell types and providing insights into developmental processes and transcriptional stochasticity. A key question is how the variety of available protocols compare in terms of their ability to detect and accurately quantify gene expression. Here, we assessed the protocol sensitivity and accuracy of many published data sets, on the basis of spike-in standards and uniform data processing. For our workflow, we developed a flexible tool for counting the number of unique molecular identifiers (https://github.com/vals/umis/). We compared 15 protocols computationally and 4 protocols experimentally for batch-matched cell populations, in addition to investigating the effects of spike-in molecular degradation. Our analysis provides an integrated framework for comparing scRNA-seq protocols.

Single-cell transcriptome analysis of fish immune cells provides insight into the evolution of vertebrate immune cell types

The immune system of vertebrate species consists of many different cell types that have distinct functional roles and are subject to different evolutionary pressures. Here, we first analyzed conservation of genes specific for all major immune cell types in human and mouse. Our results revealed higher gene turnover and faster evolution of trans-membrane proteins in NK cells compared with other immune cell types, and especially T cells, but similar conservation of nuclear and cytoplasmic protein coding genes. To validate these findings in a distant vertebrate species, we used single-cell RNA sequencing of lck:GFP cells in zebrafish and obtained the first transcriptome of specific immune cell types in a nonmammalian species. Unsupervised clustering and single-cell TCR locus reconstruction identified three cell populations, T cells, a novel type of NK-like cells, and a smaller population of myeloid-like cells. Differential expression analysis uncovered new immune-cell-specific genes, including novel immunoglobulin-like receptors, and neofunctionalization of recently duplicated paralogs. Evolutionary analyses confirmed the higher gene turnover of trans-membrane proteins in NK cells compared with T cells in fish species, suggesting that this is a general property of immune cell types across all vertebrates.

Single-Cell Multiomics: Multiple Measurements from Single Cells

Single-cell sequencing provides information that is not confounded by genotypic or phenotypic heterogeneity of bulk samples. Sequencing of one molecular type (RNA, methylated DNA or open chromatin) in a single cell, furthermore, provides insights into the cell's phenotype and links to its genotype. Nevertheless, only by taking measurements of these phenotypes and genotypes from the same single cells can such inferences be made unambiguously. In this review, we survey the first experimental approaches that assay, in parallel, multiple molecular types from the same single cell, before considering the challenges and opportunities afforded by these and future technologies.

Unambiguous inference that a cellular phenotype is caused by a genotype can only be achieved by their measurement from the same single cell.

Estimating RNA and DNA copy number abundance in single cells is now possible using a variety of experimental approaches.

Parallel measurement of single-cell epigenomes and transcriptomes provides further insight into the regulation of cellular identity and phenotypes.

Parallel measurement of single-cell transcriptomes and protein abundance (by FACS, proximity ligation assays/PEA or mass cytometry) allows insight into expression dynamics.

Our understanding of cancer progression and diversity is likely to be advanced greatly by the multiomics investigation of single cells, as is our understanding of normal developmental and other disease processes.

Ongoing technological advances will see improvements in the coverage, sensitivity of multiomics approaches, as well the number of analytes that can be surveyed in parallel.

Genome-wide in vivo screen identifies novel host regulators of metastatic colonization

Metastasis is the leading cause of death for cancer patients. This multi-stage process requires tumour cells to survive in the circulation, extravasate at distant sites, then proliferate; it involves contributions from both the tumour cell and tumour microenvironment ('host', which includes stromal cells and the immune system). Studies suggest the early steps of the metastatic process are relatively efficient, with the post-extravasation regulation of tumour growth ('colonization') being critical in determining metastatic outcome. Here we show the results of screening 810 mutant mouse lines using an in vivo assay to identify microenvironmental regulators of metastatic colonization. We identify 23 genes that, when disrupted in mouse, modify the ability of tumour cells to establish metastatic foci, with 19 of these genes not previously demonstrated to play a role in host control of metastasis. The largest reduction in pulmonary metastasis was observed in sphingosine-1-phosphate (S1P) transporter spinster homologue 2 (Spns2)-deficient mice. We demonstrate a novel outcome of S1P-mediated regulation of lymphocyte trafficking, whereby deletion of Spns2, either globally or in a lymphatic endothelial-specific manner, creates a circulating lymphopenia and a higher percentage of effector T cells and natural killer (NK) cells present in the lung. This allows for potent tumour cell killing, and an overall decreased metastatic burden.

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.

Resolving early mesoderm diversification through single-cell expression profiling

In mammals, specification of the three major germ layers occurs during gastrulation, when cells ingressing through the primitive streak differentiate into the precursor cells of major organ systems. However, the molecular mechanisms underlying this process remain unclear, as numbers of gastrulating cells are very limited. In the mouse embryo at embryonic day 6.5, cells located at the junction between the extra-embryonic region and the epiblast on the posterior side of the embryo undergo an epithelial-to-mesenchymal transition and ingress through the primitive streak. Subsequently, cells migrate, either surrounding the prospective ectoderm contributing to the embryo proper, or into the extra-embryonic region to form the yolk sac, umbilical cord and placenta. Fate mapping has shown that mature tissues such as blood and heart originate from specific regions of the pre-gastrula epiblast, but the plasticity of cells within the embryo and the function of key cell-type-specific transcription factors remain unclear. Here we analyse 1,205 cells from the epiblast and nascent Flk1(+) mesoderm of gastrulating mouse embryos using single-cell RNA sequencing, representing the first transcriptome-wide in vivo view of early mesoderm formation during mammalian gastrulation. Additionally, using knockout mice, we study the function of Tal1, a key haematopoietic transcription factor, and demonstrate, contrary to previous studies performed using retrospective assays, that Tal1 knockout does not immediately bias precursor cells towards a cardiac fate.

Heterogeneity in Oct4 and Sox2 Targets Biases Cell Fate in 4-Cell Mouse Embryos

The major and essential objective of pre-implantation development is to establish embryonic and extra-embryonic cell fates. To address when and how this fundamental process is initiated in mammals, we characterize transcriptomes of all individual cells throughout mouse pre-implantation development. This identifies targets of master pluripotency regulators Oct4 and Sox2 as being highly heterogeneously expressed between blastomeres of the 4-cell embryo, with Sox21 showing one of the most heterogeneous expression profiles. Live-cell tracking demonstrates that cells with decreased Sox21 yield more extra-embryonic than pluripotent progeny. Consistently, decreasing Sox21 results in premature upregulation of the differentiation regulator Cdx2, suggesting that Sox21 helps safeguard pluripotency. Furthermore, Sox21 is elevated following increased expression of the histone H3R26-methylase CARM1 and is lowered following CARM1 inhibition, indicating the importance of epigenetic regulation. Therefore, our results indicate that heterogeneous gene expression, as early as the 4-cell stage, initiates cell-fate decisions by modulating the balance of pluripotency and differentiation.

Single-Cell RNA-Sequencing Reveals a Continuous Spectrum of Differentiation in Hematopoietic Cells

The transcriptional programs that govern hematopoiesis have been investigated primarily by population-level analysis of hematopoietic stem and progenitor cells, which cannot reveal the continuous nature of the differentiation process. Here we applied single-cell RNA-sequencing to a population of hematopoietic cells in zebrafish as they undergo thrombocyte lineage commitment. By reconstructing their developmental chronology computationally, we were able to place each cell along a continuum from stem cell to mature cell, refining the traditional lineage tree. The progression of cells along this continuum is characterized by a highly coordinated transcriptional program, displaying simultaneous suppression of genes involved in cell proliferation and ribosomal biogenesis as the expression of lineage specific genes increases. Within this program, there is substantial heterogeneity in the expression of the key lineage regulators. Overall, the total number of genes expressed, as well as the total mRNA content of the cell, decreases as the cells undergo lineage commitment.

Dynamic spatio-temporal contribution of single β5t+ cortical epithelial precursors to the thymus medulla

Intrathymic T‐cell development is critically dependent on cortical and medullary thymic epithelial cells (TECs). Both epithelial subsets originate during early thymus organogenesis from progenitor cells that express the thymoproteasome subunit β5t, a typical feature of cortical TECs. Using in vivo lineage fate mapping, we demonstrate in mice that β5t⁺ TEC progenitors give rise to the medullary TEC compartment early in life but significantly limit their contribution once the medulla has completely formed. Lineage‐tracing studies at single cell resolution demonstrate for young mice that the postnatal medulla is expanded from individual β5t⁺ cortical progenitors located at the cortico‐medullary junction. These results therefore not only define a developmental window during which the expansion of medulla is efficiently enabled by progenitors resident in the thymic cortex, but also reveal the spatio‐temporal dynamics that control the growth of the thymic medulla.

Parallel single-cell sequencing links transcriptional and epigenetic heterogeneity

We report scM&T-seq, a method for parallel single-cell genome-wide methylome and transcriptome sequencing, allowing discovery of associations between transcriptional and epigenetic variation. Profiling of 61 mouse embryonic stem cells confirmed known links between DNA methylation and transcription. Notably, the method reveals novel associations between heterogeneously methylated distal regulatory elements and transcription of key pluripotency genes.

Integrin Based Isolation Enables Purification of Murine Lineage Committed Cardiomyocytes

In contrast to mature cardiomyocytes which have limited regenerative capacity, pluripotent stem cells represent a promising source for the generation of new cardiomyocytes. The tendency of pluripotent stem cells to form teratomas and the heterogeneity from various differentiation stages and cardiomyocyte cell sub-types, however, are major obstacles to overcome before this type of therapy could be applied in a clinical setting. Thus, the identification of extracellular markers for specific cardiomyocyte progenitors and mature subpopulations is of particular importance. The delineation of cardiomyocyte surface marker patterns not only serves as a means to derive homogeneous cell populations by FACS, but is also an essential tool to understand cardiac development. By using single-cell expression profiling in early mouse embryonic hearts, we found that a combination of integrin alpha-1, alpha-5, alpha-6 and N-cadherin enables isolation of lineage committed murine cardiomyocytes. Additionally, we were able to separate trabecular cardiomyocytes from solid ventricular myocardium and atrial murine cells. These cells exhibit expected subtype specific phenotype confirmed by electrophysiological analysis. We show that integrin expression can be used for the isolation of living, functional and lineage-specific murine cardiomyocytes.

Combined Single-Cell Functional and Gene Expression Analysis Resolves Heterogeneity within Stem Cell Populations

Heterogeneity within the self-renewal durability of adult hematopoietic stem cells (HSCs) challenges our understanding of the molecular framework underlying HSC function. Gene expression studies have been hampered by the presence of multiple HSC subtypes and contaminating non-HSCs in bulk HSC populations. To gain deeper insight into the gene expression program of murine HSCs, we combined single-cell functional assays with flow cytometric index sorting and single-cell gene expression assays. Through bioinformatic integration of these datasets, we designed an unbiased sorting strategy that separates non-HSCs away from HSCs, and single-cell transplantation experiments using the enriched population were combined with RNA-seq data to identify key molecules that associate with long-term durable self-renewal, producing a single-cell molecular dataset that is linked to functional stem cell activity. Finally, we demonstrated the broader applicability of this approach for linking key molecules with defined cellular functions in another stem cell system.

Decoding the regulatory network of early blood development from single-cell gene expression measurements

Reconstruction of the molecular pathways controlling organ development has been hampered by a lack of methods to resolve embryonic progenitor cells. Here we describe a strategy to address this problem that combines gene expression profiling of large numbers of single cells with data analysis based on diffusion maps for dimensionality reduction and network synthesis from state transition graphs. Applying the approach to hematopoietic development in the mouse embryo, we map the progression of mesoderm toward blood using single-cell gene expression analysis of 3,934 cells with blood-forming potential captured at four time points between E7.0 and E8.5. Transitions between individual cellular states are then used as input to develop a single-cell network synthesis toolkit to generate a computationally executable transcriptional regulatory network model of blood development. Several model predictions concerning the roles of Sox and Hox factors are validated experimentally. Our results demonstrate that single-cell analysis of a developing organ coupled with computational approaches can reveal the transcriptional programs that underpin organogenesis.

Population and single-cell genomics reveal the Aire dependency, relief from Polycomb silencing, and distribution of self-antigen expression in thymic epithelia

Promiscuous gene expression (PGE) by thymic epithelial cells (TEC) is essential for generating a diverse T cell antigen receptor repertoire tolerant to self-antigens, and thus for avoiding autoimmunity. Nevertheless, the extent and nature of this unusual expression program within TEC populations and single cells are unknown. Using deep transcriptome sequencing of carefully identified mouse TEC subpopulations, we discovered a program of PGE that is common between medullary (m) and cortical TEC, further elaborated in mTEC, and completed in mature mTEC expressing the autoimmune regulator gene (Aire). TEC populations are capable of expressing up to 19,293 protein-coding genes, the highest number of genes known to be expressed in any cell type. Remarkably, in mouse mTEC, Aire expression alone positively regulates 3980 tissue-restricted genes. Notably, the tissue specificities of these genes include known targets of autoimmunity in human AIRE deficiency. Led by the observation that genes induced by Aire expression are generally characterized by a repressive chromatin state in somatic tissues, we found these genes to be strongly associated with H3K27me3 marks in mTEC. Our findings are consistent with AIRE targeting and inducing the promiscuous expression of genes previously epigenetically silenced by Polycomb group proteins. Comparison of the transcriptomes of 174 single mTEC indicates that genes induced by Aire expression are transcribed stochastically at low cell frequency. Furthermore, when present, Aire expression-dependent transcript levels were 16-fold higher, on average, in individual TEC than in the mTEC population.

Transcriptional diversity during lineage commitment of human blood progenitors

Blood cells derive from hematopoietic stem cells through stepwise fating events. To characterize gene expression programs driving lineage choice, we sequenced RNA from eight primary human hematopoietic progenitor populations representing the major myeloid commitment stages and the main lymphoid stage. We identified extensive cell type-specific expression changes: 6711 genes and 10,724 transcripts, enriched in non-protein-coding elements at early stages of differentiation. In addition, we found 7881 novel splice junctions and 2301 differentially used alternative splicing events, enriched in genes involved in regulatory processes. We demonstrated experimentally cell-specific isoform usage, identifying nuclear factor I/B (NFIB) as a regulator of megakaryocyte maturation-the platelet precursor. Our data highlight the complexity of fating events in closely related progenitor populations, the understanding of which is essential for the advancement of transplantation and regenerative medicine.

Myelodysplastic syndromes are propagated by rare and distinct human cancer stem cells in vivo

Evidence for distinct human cancer stem cells (CSCs) remains contentious and the degree to which different cancer cells contribute to propagating malignancies in patients remains unexplored. In low- to intermediate-risk myelodysplastic syndromes (MDS), we establish the existence of rare multipotent MDS stem cells (MDS-SCs), and their hierarchical relationship to lineage-restricted MDS progenitors. All identified somatically acquired genetic lesions were backtracked to distinct MDS-SCs, establishing their distinct MDS-propagating function in vivo. In isolated del(5q)-MDS, acquisition of del(5q) preceded diverse recurrent driver mutations. Sequential analysis in del(5q)-MDS revealed genetic evolution in MDS-SCs and MDS-progenitors prior to leukemic transformation. These findings provide definitive evidence for rare human MDS-SCs in vivo, with extensive implications for the targeting of the cells required and sufficient for MDS-propagation.

Lymphomyeloid contribution of an immune-restricted progenitor emerging prior to definitive hematopoietic stem cells

In jawed vertebrates, development of an adaptive immune-system is essential for protection of the born organism against otherwise life-threatening pathogens. Myeloid cells of the innate immune system are formed early in development, whereas lymphopoiesis has been suggested to initiate much later, following emergence of definitive hematopoietic stem cells (HSCs). Herein, we demonstrate that the embryonic lymphoid commitment process initiates earlier than previously appreciated, prior to emergence of definitive HSCs, through establishment of a previously unrecognized entirely immune-restricted and lymphoid-primed progenitor. Notably, this immune-restricted progenitor appears to first emerge in the yolk sac and contributes physiologically to the establishment of lymphoid and some myeloid components of the immune-system, establishing the lymphomyeloid lineage restriction process as an early and physiologically important lineage-commitment step in mammalian hematopoiesis.

Platelet-biased stem cells reside at the apex of the haematopoietic stem-cell hierarchy

The blood system is maintained by a small pool of haematopoietic stem cells (HSCs), which are required and sufficient for replenishing all human blood cell lineages at millions of cells per second throughout life. Megakaryocytes in the bone marrow are responsible for the continuous production of platelets in the blood, crucial for preventing bleeding--a common and life-threatening side effect of many cancer therapies--and major efforts are focused at identifying the most suitable cellular and molecular targets to enhance platelet production after bone marrow transplantation or chemotherapy. Although it has become clear that distinct HSC subsets exist that are stably biased towards the generation of lymphoid or myeloid blood cells, we are yet to learn whether other types of lineage-biased HSC exist or understand their inter-relationships and how differently lineage-biased HSCs are generated and maintained. The functional relevance of notable phenotypic and molecular similarities between megakaryocytes and bone marrow cells with an HSC cell-surface phenotype remains unclear. Here we identify and prospectively isolate a molecularly and functionally distinct mouse HSC subset primed for platelet-specific gene expression, with enhanced propensity for short- and long-term reconstitution of platelets. Maintenance of platelet-biased HSCs crucially depends on thrombopoietin, the primary extrinsic regulator of platelet development. Platelet-primed HSCs also frequently have a long-term myeloid lineage bias, can self-renew and give rise to lymphoid-biased HSCs. These findings show that HSC subtypes can be organized into a cellular hierarchy, with platelet-primed HSCs at the apex. They also demonstrate that molecular and functional priming for platelet development initiates already in a distinct HSC population. The identification of a platelet-primed HSC population should enable the rational design of therapies enhancing platelet output.

Canonical Wnt signaling negatively regulates platelet function

Wnts regulate important intracellular signaling events, and dysregulation of the Wnt pathway has been linked to human disease. Here, we uncover numerous Wnt canonical effectors in human platelets where Wnts, their receptors, and downstream signaling components have not been previously described. We demonstrate that the Wnt3a ligand inhibits platelet adhesion, activation, dense granule secretion, and aggregation. Wnt3a also altered platelet shape change and inhibited the activation of the small GTPase RhoA. In addition, we found the Wnt-beta-catenin signaling pathway to be functional in platelets. Finally, disruption of the Wnt Frizzled 6 receptor in the mouse resulted in a hyperactivatory platelet phenotype and a reduced sensitivity to Wnt3a. Taken together our studies reveal a novel functional role for Wnt signaling in regulating anucleate platelet function and may provide a tractable target for future antiplatelet therapy.

Identification of variation in the platelet transcriptome associated with glycoprotein 6 haplotype

Platelet Glycoprotein VI (GPVI) is the activatory collagen signalling receptor that transmits an outside-in signal via the FcR gamma-chain. In Caucasians two GP6 haplotypes have been identified which encode GPVI isoforms that differ by five amino-acids. The minor haplotype is associated with a modest but statistically significant reduction in GPVI abundance and reduced downstream signalling events. As GPVI is also expressed on megakaryocytes, different GPVI isoforms may imprint on the platelet transcriptome. We investigated the association of GP6 haplotype with transcription by comparing the transcriptomes of platelets from individuals homozygous for the major ('a') and minor ('b') haplotypes to identify differentially expressed (DE) transcripts. Platelet RNA was isolated from apheresis concentrates from 16 'aa' donors and eight 'bb' donors. mRNA was amplified using a template-switching PCR based protocol and fluorescently labelled. Samples were randomly paired both within and between haplotypes and compared on a cDNA microarray. No consistently DE transcripts were identified within the 'aa' haplotype but 52 significantly DE transcripts were observed between haplotypes. Generally the fold differences were low (two to four-fold) but were confirmed by qRT-PCR for selected transcripts (TUBB1, P = 0.0004; VWF, P = 0.0126). The results of this study indicate that there are subtle differences between the platelet transcriptomes of individuals who differ by GP6 haplotype. The identification of DE genes may identify critical pathways and nodes not previously known to be involved in platelet development and function.

The novel inhibitory receptor G6B is expressed on the surface of platelets and attenuates platelet function in vitro

The G6B cell-surface receptor, which contains a single Ig-like domain, has been shown to bind to SHP-1 and SHP-2 after phosphorylation of 2 immunoreceptor tyrosine-based inhibitory motifs (ITIMs) in its cytoplasmic tail, classifying this protein as a new member of the family of inhibitory receptors. In this study, we demonstrate by real-time polymerase chain reaction (PCR) and Western-blot analysis that G6B is expressed on platelets. Cross-linking of G6B with polyclonal antisera has a significant inhibitory effect on platelet aggregation and activation by agonists such as ADP and collagen-related peptide (CRP). This inhibitory function of G6B appears to operate in a calcium-independent manner. Our results suggest that G6B represents a novel inhibitory receptor found on the surface of platelets and that it could be an antithrombotic drug target.

Comparative gene expression profiling of in vitro differentiated megakaryocytes and erythroblasts identifies novel activatory and inhibitory platelet membrane proteins

To identify previously unknown platelet receptors we compared the transcriptomes of in vitro differentiated megakaryocytes (MKs) and erythroblasts (EBs). RNA was obtained from purified, biologically paired MK and EB cultures and compared using cDNA microarrays. Bioinformatical analysis of MK-up-regulated genes identified 151 transcripts encoding transmembrane domain-containing proteins. Although many of these were known platelet genes, a number of previously unidentified or poorly characterized transcripts were also detected. Many of these transcripts, including G6b, G6f, LRRC32, LAT2, and the G protein-coupled receptor SUCNR1, encode proteins with structural features or functions that suggest they may be involved in the modulation of platelet function. Immunoblotting on platelets confirmed the presence of the encoded proteins, and flow cytometric analysis confirmed the expression of G6b, G6f, and LRRC32 on the surface of platelets. Through comparative analysis of expression in platelets and other blood cells we demonstrated that G6b, G6f, and LRRC32 are restricted to the platelet lineage, whereas LAT2 and SUCNR1 were also detected in other blood cells. The identification of the succinate receptor SUCNR1 in platelets is of particular interest, because physiologically relevant concentrations of succinate were shown to potentiate the effect of low doses of a variety of platelet agonists.

Platelet genomics and proteomics in human health and disease

Proteomic and genomic technologies provide powerful tools for characterizing the multitude of events that occur in the anucleate platelet. These technologies are beginning to define the complete platelet transcriptome and proteome as well as the protein-protein interactions critical for platelet function. The integration of these results provides the opportunity to identify those proteins involved in discrete facets of platelet function. Here we summarize the findings of platelet proteome and transcriptome studies and their application to diseases of platelet function.

Analysing the platelet transcriptome

Historical studies, using Northern blot hybridization, RT-PCR and cDNA library construction have demonstrated the presence of a variety of mRNA molecules in platelets. The development of microarray technology has allowed further characterization of the transcripts represented in the platelet transcriptome. In this review, these studies will be summarized and their findings in relation to the study of platelet function and the identification of disease risk genes discussed.