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-nucleus RNA-seq dissection of choroid plexus tumor cell heterogeneity.

The genomic, genetic and cellular events regulating the onset, growth and survival of rare, choroid plexus neoplasms remain poorly understood. Here, we examine the heterogeneity of human choroid plexus tumors by single-nucleus transcriptome analysis of 23,906 cells from four disease-free choroid plexus and eleven choroid plexus tumors. The resulting expression atlas profiles cellular and transcriptional diversity, copy number alterations, and cell-cell interaction networks in normal and cancerous choroid plexus. In choroid plexus tumor epithelial cells, we observe transcriptional changes that correlate with genome-wide methylation profiles. We further characterize tumor type-specific stromal microenvironments that include altered macrophage and mesenchymal cell states, as well as changes in extracellular matrix components. This first single-cell dataset resource from such scarce samples should be valuable for divising therapies against these little-studied neoplasms.

Distinct myeloid progenitor-differentiation pathways identified through single-cell RNA sequencing.

According to current models of hematopoiesis, lymphoid-primed multi-potent progenitors (LMPPs) (Lin(-)Sca-1(+)c-Kit(+)CD34(+)Flt3(hi)) and common myeloid progenitors (CMPs) (Lin(-)Sca-1(+)c-Kit(+)CD34(+)CD41(hi)) establish an early branch point for separate lineage-commitment pathways from hematopoietic stem cells, with the notable exception that both pathways are proposed to generate all myeloid innate immune cell types through the same myeloid-restricted pre-granulocyte-macrophage progenitor (pre-GM) (Lin(-)Sca-1(-)c-Kit(+)CD41(-)FcγRII/III(-)CD150(-)CD105(-)). By single-cell transcriptome profiling of pre-GMs, we identified distinct myeloid differentiation pathways: a pathway expressing the gene encoding the transcription factor GATA-1 generated mast cells, eosinophils, megakaryocytes and erythroid cells, and a pathway lacking expression of that gene generated monocytes, neutrophils and lymphocytes. These results identify an early hematopoietic-lineage bifurcation that separates the myeloid lineages before their segregation from other hematopoietic-lineage potential.

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.

Single-Cell Analyses Reveal Megakaryocyte-Biased Hematopoiesis in Myelofibrosis and Identify Mutant Clone-Specific Targets.

Myelofibrosis is a severe myeloproliferative neoplasm characterized by increased numbers of abnormal bone marrow megakaryocytes that induce fibrosis, destroying the hematopoietic microenvironment. To determine the cellular and molecular basis for aberrant megakaryopoiesis in myelofibrosis, we performed single-cell transcriptome profiling of 135,929 CD34+ lineage- hematopoietic stem and progenitor cells (HSPCs), single-cell proteomics, genomics, and functional assays. We identified a bias toward megakaryocyte differentiation apparent from early multipotent stem cells in myelofibrosis and associated aberrant molecular signatures. A sub-fraction of myelofibrosis megakaryocyte progenitors (MkPs) are transcriptionally similar to healthy-donor MkPs, but the majority are disease specific, with distinct populations expressing fibrosis- and proliferation-associated genes. Mutant-clone HSPCs have increased expression of megakaryocyte-associated genes compared to wild-type HSPCs, and we provide early validation of G6B as a potential immunotherapy target. Our study paves the way for selective targeting of the myelofibrosis clone and illustrates the power of single-cell multi-omics to discover tumor-specific therapeutic targets and mediators of tissue fibrosis.

Unravelling Intratumoral Heterogeneity through High-Sensitivity Single-Cell Mutational Analysis and Parallel RNA Sequencing.

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for resolving transcriptional heterogeneity. However, its application to studying cancerous tissues is currently hampered by the lack of coverage across key mutation hotspots in the vast majority of cells; this lack of coverage prevents the correlation of genetic and transcriptional readouts from the same single cell. To overcome this, we developed TARGET-seq, a method for the high-sensitivity detection of multiple mutations within single cells from both genomic and coding DNA, in parallel with unbiased whole-transcriptome analysis. Applying TARGET-seq to 4,559 single cells, we demonstrate how this technique uniquely resolves transcriptional and genetic tumor heterogeneity in myeloproliferative neoplasms (MPN) stem and progenitor cells, providing insights into deregulated pathways of mutant and non-mutant cells. TARGET-seq is a powerful tool for resolving the molecular signatures of genetically distinct subclones of cancer cells.

Deep sequencing of short capped RNAs reveals novel families of noncoding RNAs.

In eukaryotes, capped RNAs include long transcripts such as messenger RNAs and long noncoding RNAs, as well as shorter transcripts such as spliceosomal RNAs, small nucleolar RNAs, and enhancer RNAs. Long capped transcripts can be profiled using cap analysis gene expression (CAGE) sequencing and other methods. Here, we describe a sequencing library preparation protocol for short capped RNAs, apply it to a differentiation time course of the human cell line THP-1, and systematically compare the landscape of short capped RNAs to that of long capped RNAs. Transcription initiation peaks associated with genes in the sense direction have a strong preference to produce either long or short capped RNAs, with one out of six peaks detected in the short capped RNA libraries only. Gene-associated short capped RNAs have highly specific 3' ends, typically overlapping splice sites. Enhancers also preferentially generate either short or long capped RNAs, with 10% of enhancers observed in the short capped RNA libraries only. Enhancers producing either short or long capped RNAs show enrichment for GWAS-associated disease SNPs. We conclude that deep sequencing of short capped RNAs reveals new families of noncoding RNAs and elucidates the diversity of transcripts generated at known and novel promoters and enhancers.

MARVEL: an integrated alternative splicing analysis platform for single-cell RNA sequencing data.

Alternative splicing is an important source of heterogeneity underlying gene expression between individual cells but remains an understudied area due to the paucity of computational tools to analyze splicing dynamics at single-cell resolution. Here, we present MARVEL, a comprehensive R package for single-cell splicing analysis applicable to RNA sequencing generated from the plate- and droplet-based methods. We performed extensive benchmarking of MARVEL against available tools and demonstrated its utility by analyzing multiple publicly available datasets in diverse cell types, including in disease. MARVEL enables systematic and integrated splicing and gene expression analysis of single cells to characterize the splicing landscape and reveal biological insights.

CD34+CD19-CD22+ B-cell progenitors may underlie phenotypic escape in patients treated with CD19-directed therapies.

CD19-directed immunotherapies have revolutionized the treatment of advanced B-cell acute lymphoblastic leukemia (B-ALL). Despite initial impressive rates of complete remission (CR) many patients ultimately relapse. Patients with B-ALL successfully treated with CD19-directed T cells eventually relapse, which, coupled with the early onset of CD22 expression during B-cell development, suggests that preexisting CD34+CD22+CD19- (pre)-leukemic cells represent an "early progenitor origin-related" mechanism underlying phenotypic escape to CD19-directed immunotherapies. We demonstrate that CD22 expression precedes CD19 expression during B-cell development. CD34+CD19-CD22+ cells are found in diagnostic and relapsed bone marrow samples of ∼70% of patients with B-ALL, and their frequency increases twofold in patients with B-ALL in CR after CD19 CAR T-cell therapy. The median of CD34+CD19-CD22+ cells before treatment was threefold higher in patients in whom B-ALL relapsed after CD19-directed immunotherapy (median follow-up, 24 months). Fluorescence in situ hybridization analysis in flow-sorted cell populations and xenograft modeling revealed that CD34+CD19-CD22+ cells harbor the genetic abnormalities present at diagnosis and initiate leukemogenesis in vivo. Our data suggest that preleukemic CD34+CD19-CD22+ progenitors underlie phenotypic escape after CD19-directed immunotherapies and reinforce ongoing clinical studies aimed at CD19/CD22 dual targeting as a strategy for reducing CD19- relapses. The implementation of CD34/CD19/CD22 immunophenotyping in clinical laboratories for initial diagnosis and subsequent monitoring of patients with B-ALL during CD19-targeted therapy is encouraged.

The quiescent fraction of chronic myeloid leukemic stem cells depends on BMPR1B, Stat3 and BMP4-niche signals to persist in patients in remission.

Chronic myelogenous leukemia arises from the transformation of hematopoietic stem cells by the BCR-ABL oncogene. Though transformed cells are predominantly BCR-ABL-dependent and sensitive to tyrosine kinase inhibitor treatment, some BMPR1B+ leukemic stem cells are treatment-insensitive and rely, among others, on the bone morphogenetic protein (BMP) pathway for their survival via a BMP4 autocrine loop. Here, we further studied the involvement of BMP signaling in favoring residual leukemic stem cell persistence in the bone marrow of patients having achieved remission under treatment. We demonstrate by single-cell RNA-Seq analysis that a sub-fraction of surviving BMPR1B+ leukemic stem cells are co-enriched in BMP signaling, quiescence and stem cell signatures, without modulation of the canonical BMP target genes, but enrichment in actors of the Jak2/Stat3 signaling pathway. Indeed, based on a new model of persisting CD34+CD38- leukemic stem cells, we show that BMPR1B+ cells display co-activated Smad1/5/8 and Stat3 pathways. Interestingly, we reveal that only the BMPR1B+ cells adhering to stromal cells display a quiescent status. Surprisingly, this quiescence is induced by treatment, while non-adherent BMPR1B+ cells treated with tyrosine kinase inhibitors continued to proliferate. The subsequent targeting of BMPR1B and Jak2 pathways decreased quiescent leukemic stem cells by promoting their cell cycle re-entry and differentiation. Moreover, while Jak2-inhibitors alone increased BMP4 production by mesenchymal cells, the addition of the newly described BMPR1B inhibitor (E6201) impaired BMP4-mediated production by stromal cells. Altogether, our data demonstrate that targeting both BMPR1B and Jak2/Stat3 efficiently impacts persisting and dormant leukemic stem cells hidden in their bone marrow microenvironment.

VALERIE: Visual-based inspection of alternative splicing events at single-cell resolution.

We present VALERIE (Visualising alternative splicing events from single-cell ribonucleic acid-sequencing experiments), an R package for visualising alternative splicing events at single-cell resolution. To explore any given specified genomic region, corresponding to an alternative splicing event, VALERIE generates an ensemble of informative plots to visualise cell-to-cell heterogeneity of alternative splicing profiles across single cells and performs statistical tests to compare percent spliced-in (PSI) values across the user-defined groups of cells. Among the features available, VALERIE displays PSI values, in lieu of read coverage, which is more suitable for representing alternative splicing profiles for a large number of samples typically generated by single-cell RNA-sequencing experiments. VALERIE is available on the Comprehensive R Archive Network (CRAN): https://cran.r-project.org/web/packages/VALERIE/index.html.

In vivo live imaging of RNA polymerase II transcription factories in primary cells.

Transcription steps are marked by different modifications of the C-terminal domain of RNA polymerase II (RNAPII). Phosphorylation of Ser5 and Ser7 by cyclin-dependent kinase 7 (CDK7) as part of TFIIH marks initiation, whereas phosphorylation of Ser2 by CDK9 marks elongation. These processes are thought to take place in localized transcription foci in the nucleus, known as "transcription factories," but it has been argued that the observed clusters/foci are mere fixation or labeling artifacts. We show that transcription factories exist in living cells as distinct foci by live-imaging fluorescently labeled CDK9, a kinase known to associate with active RNAPII. These foci were observed in different cell types derived from CDK9-mCherry knock-in mice. We show that these foci are very stable while highly dynamic in exchanging CDK9. Chromatin immunoprecipitation (ChIP) coupled with deep sequencing (ChIP-seq) data show that the genome-wide binding sites of CDK9 and initiating RNAPII overlap on transcribed genes. Immunostaining shows that CDK9-mCherry foci colocalize with RNAPII-Ser5P, much less with RNAPII-Ser2P, and not with CDK12 (a kinase reported to be involved in the Ser2 phosphorylation) or with splicing factor SC35. In conclusion, transcription factories exist in living cells, and initiation and elongation of transcripts takes place in different nuclear compartments.

The DNA-binding protein CTCF limits proximal Vκ recombination and restricts κ enhancer interactions to the immunoglobulin κ light chain locus.

Regulation of immunoglobulin (Ig) V(D)J gene rearrangement is dependent on higher-order chromatin organization. Here, we studied the in vivo function of the DNA-binding zinc-finger protein CTCF, which regulates interactions between enhancers and promoters. By conditional deletion of the Ctcf gene in the B cell lineage, we demonstrate that loss of CTCF allowed Ig heavy chain recombination, but pre-B cell proliferation and differentiation was severely impaired. In the absence of CTCF, the Igκ light chain locus showed increased proximal and reduced distal Vκ usage. This was associated with enhanced proximal Vκ and reduced Jκ germline transcription. Chromosome conformation capture experiments demonstrated that CTCF limits interactions of the Igκ enhancers with the proximal V(κ) gene region and prevents inappropriate interactions between these strong enhancers and elements outside the Igκ locus. Thus, although Ig gene recombination can occur in the absence of CTCF, it is a critical factor determining Vκ segment choice for recombination.

Cell-intrinsic depletion of Aml1-ETO-expressing pre-leukemic hematopoietic stem cells by K-Ras activating mutation.

Somatic mutations in acute myeloid leukemia are acquired sequentially and hierarchically. First, pre-leukemic mutations, such as t(8;21) that encodes AML1-ETO, are acquired within the hematopoietic stem cell (HSC) compartment, while signaling pathway mutations, including KRAS activating mutations, are late events acquired during transformation of leukemic progenitor cells and are rarely detectable in HSC. This raises the possibility that signaling pathway mutations are detrimental to clonal expansion of pre-leukemic HSC. To address this hypothesis, we used conditional genetics to introduce Aml1-ETO and K-RasG12D into murine HSC, either individually or in combination. In the absence of activated Ras, Aml1-ETO-expressing HSC conferred a competitive advantage. However, activated K-Ras had a marked detrimental effect on Aml1-ETO-expressing HSC, leading to loss of both phenotypic and functional HSC. Cell cycle analysis revealed a loss of quiescence in HSC co-expressing Aml1-ETO and K-RasG12D, accompanied by an enrichment in E2F and Myc target gene expression and depletion of HSC self-renewal-associated gene expression. These findings provide a mechanistic basis for the observed absence of KRAS signaling mutations in the pre-malignant HSC compartment.

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.

The genome-wide dynamics of the binding of Ldb1 complexes during erythroid differentiation.

One of the complexes formed by the hematopoietic transcription factor Gata1 is a complex with the Ldb1 (LIM domain-binding protein 1) and Tal1 proteins. It is known to be important for the development and differentiation of the erythroid cell lineage and is thought to be implicated in long-range interactions. Here, the dynamics of the composition of the complex-in particular, the binding of the negative regulators Eto2 and Mtgr1-are studied, in the context of their genome-wide targets. This shows that the complex acts almost exclusively as an activator, binding a very specific combination of sequences, with a positioning relative to transcription start site, depending on the type of the core promoter. The activation is accompanied by a net decrease in the relative binding of Eto2 and Mtgr1. A Chromosome Conformation Capture sequencing (3C-seq) assay also shows that the binding of the Ldb1 complex marks genomic interaction sites in vivo. This establishes the Ldb1 complex as a positive regulator of the final steps of erythroid differentiation that acts through the shedding of negative regulators and the active interaction between regulatory sequences.

Dynamic long-range chromatin interactions control Myb proto-oncogene transcription during erythroid development.

The key haematopoietic regulator Myb is essential for coordinating proliferation and differentiation. ChIP-Sequencing and Chromosome Conformation Capture (3C)-Sequencing were used to characterize the structural and protein-binding dynamics of the Myb locus during erythroid differentiation. In proliferating cells expressing Myb, enhancers within the Myb-Hbs1l intergenic region were shown to form an active chromatin hub (ACH) containing the Myb promoter and first intron. This first intron was found to harbour the transition site from transcription initiation to elongation, which takes place around a conserved CTCF site. Upon erythroid differentiation, Myb expression is downregulated and the ACH destabilized. We propose a model for Myb activation by distal enhancers dynamically bound by KLF1 and the GATA1/TAL1/LDB1 complex, which primarily function as a transcription elongation element through chromatin looping.

Genome-wide analysis shows that Ldb1 controls essential hematopoietic genes/pathways in mouse early development and reveals novel players in hematopoiesis.

The first site exhibiting hematopoietic activity in mammalian development is the yolk-sac blood island, which originates from the hemangioblast. Here we performed differentiation assays, as well as genome-wide molecular and functional studies in blast colony-forming cells to gain insight into the function of the essential Ldb1 factor in early primitive hematopoietic development. We show that the previously reported lack of yolk-sac hematopoiesis and vascular development in Ldb1(-/-) mouse result from a decreased number of hemangioblasts and a block in their ability to differentiate into erythroid and endothelial progenitor cells. Transcriptome analysis and correlation with the genome-wide binding pattern of Ldb1 in hemangioblasts revealed a number of direct-target genes and pathways misregulated in the absence of Ldb1. The regulation of essential developmental factors by Ldb1 defines it as an upstream transcriptional regulator of hematopoietic/endothelial development. We show the complex interplay that exists between transcription factors and signaling pathways during the very early stages of hematopoietic/endothelial development and the specific signaling occurring in hemangioblasts in contrast to more advanced hematopoietic developmental stages. Finally, by revealing novel genes and pathways not previously associated with early development, our study provides novel candidate targets to manipulate the differentiation of hematopoietic and/or endothelial cells.

r3Cseq: an R/Bioconductor package for the discovery of long-range genomic interactions from chromosome conformation capture and next-generation sequencing data.

The coupling of chromosome conformation capture (3C) with next-generation sequencing technologies enables the high-throughput detection of long-range genomic interactions, via the generation of ligation products between DNA sequences, which are closely juxtaposed in vivo. These interactions involve promoter regions, enhancers and other regulatory and structural elements of chromosomes and can reveal key details of the regulation of gene expression. 3C-seq is a variant of the method for the detection of interactions between one chosen genomic element (viewpoint) and the rest of the genome. We present r3Cseq, an R/Bioconductor package designed to perform 3C-seq data analysis in a number of different experimental designs. The package reads a common aligned read input format, provides data normalization, allows the visualization of candidate interaction regions and detects statistically significant chromatin interactions, thus greatly facilitating hypothesis generation and the interpretation of experimental results. We further demonstrate its use on a series of real-world applications.