Search

Supplementary Figure 1. Pro-inflammatory cytokines are elevated in MF mice and reversed with JAK inhibitor treatment. Supplementary Figure 2. Principal component analysis of MF cells. Supplementary Figure 3. Single cell analysis of immunophenotypically-defined subpopulations from MF and control mice. Supplementary Figure 4. Elevated IL8 expression in MF patients. Supplementary Figure 5. Complete removal of Il6 from the mouse system shows only minor effects on MPLW515L disease. Supplementary Figure 6. Peripheral blood and cytometric flow analysis of MPLW515L transplanted mice with and without Stat3 deletion. Supplementary Figure 7. Constitutive activation of STAT3 in the BM of MPLW515L mice and primary MF patients. Supplementary Figure 8. Schematic illustration of bone marrow transplantation experiments using Stat3-deficient mice or littermate control mice as donors. Supplementary Figure 9. Excision of Stat3 post engraftment prolongs survival and reduces leukocytosis. Supplementary Figure 10. Mutant-restricted deletion of Stat3 does not affect disease severity in vivo. Supplementary Figure 11. Single cell analysis of populations from MF and control mice. Supplementary Figure 12. Mutant and non-mutant hematopoietic cells feature a higher percentage of cytokine secreting cells. Supplementary Figure 13. Single cell analysis of populations from Jak2V617F MF and control mice. Supplementary Figure 14. Mutant and non-mutant hematopoietic cells feature a higher percentage of cytokine secreting cells. Supplementary Figure 15. JAK1/2 inhibition reduces cytokine expression from CD45.1-positive wild-type cells and CD45.2-positive mutant cell population. Supplementary Table 1. Primary myelofibrosis patient characteristics. Supplementary Table 2. Genetic alterations of patient PMF#11.

Pooled CRISPR screens coupled with single-cell RNA-sequencing have enabled systematic interrogation of gene function and regulatory networks. Here, we introduce Cas13 RNA Perturb-seq (CaRPool-seq), which leverages the RNA-targeting CRISPR-Cas13d system and enables efficient combinatorial perturbations alongside multimodal single-cell profiling. CaRPool-seq encodes multiple perturbations on a cleavable CRISPR array that is associated with a detectable barcode sequence, allowing for the simultaneous targeting of multiple genes. We compared CaRPool-seq to existing Cas9-based methods, highlighting its unique strength to efficiently profile combinatorially perturbed cells. Finally, we apply CaRPool-seq to perform multiplexed combinatorial perturbations of myeloid differentiation regulators in an acute myeloid leukemia (AML) model system and identify extensive interactions between different chromatin regulators that can enhance or suppress AML differentiation phenotypes.

Single cell RNA sequencing (scRNA-seq) has become a core tool for researchers to understand biology. As scRNA-seq has become more ubiquitous, many applications demand higher scalability and sensitivity. Split-pool combinatorial barcoding makes it possible to scale projects to hundreds of samples and millions of cells, overcoming limitations of previous droplet based technologies. However, there is still a need for increased sensitivity for both droplet and combinatorial barcoding based scRNA-seq technologies. To meet this need, here we introduce an updated combinatorial barcoding method for scRNA-seq with dramatically improved sensitivity. To assess performance, we profile a variety of sample types, including cell lines, human peripheral blood mononuclear cells (PBMCs), mouse brain nuclei, and mouse liver nuclei. When compared to the previously best performing approach, we find up to a 2.6-fold increase in unique transcripts detected per cell and up to a 1.8-fold increase in genes detected per cell. These improvements to transcript and gene detection increase the resolution of the resulting data, making it easier to distinguish cell types and states in heterogeneous samples. Split-pool combinatorial barcoding already enables scaling to millions of cells, the ability to perform scRNA-seq on previously fixed and frozen samples, and access to scRNA-seq without the need to purchase specialized lab equipment. Our hope is that by combining these previous advantages with the dramatic improvements to sensitivity presented here, we will elevate the standards and capabilities of scRNA-seq for the broader community.

Recent technological advances have enabled massively parallel chromatin profiling with scATAC-seq (single-cell assay for transposase accessible chromatin by sequencing). Here we present ATAC with select antigen profiling by sequencing (ASAP-seq), a tool to simultaneously profile accessible chromatin and protein levels. Our approach pairs sparse scATAC-seq data with robust detection of hundreds of cell surface and intracellular protein markers and optional capture of mitochondrial DNA for clonal tracking, capturing three distinct modalities in single cells. ASAP-seq uses a bridging approach that repurposes antibody:oligonucleotide conjugates designed for existing technologies that pair protein measurements with single-cell RNA sequencing. Together with DOGMA-seq, an adaptation of CITE-seq (cellular indexing of transcriptomes and epitopes by sequencing) for measuring gene activity across the central dogma of gene regulation, we demonstrate the utility of systematic multi-omic profiling by revealing coordinated and distinct changes in chromatin, RNA and surface proteins during native hematopoietic differentiation and peripheral blood mononuclear cell stimulation and as a combinatorial decoder and reporter of multiplexed perturbations in primary T cells.

Continuous supply of immune cells throughout life relies on the delicate balance in the hematopoietic stem cell (HSC) pool between long-term maintenance and meeting the demands of both normal blood production and unexpected stress conditions. Here we identified distinct subsets of human long-term (LT)-HSCs that responded differently to regeneration-mediated stress: an immune checkpoint ligand CD112lo subset that exhibited a transient engraftment restraint (termed latency) before contributing to hematopoietic reconstitution and a primed CD112hi subset that responded rapidly. This functional heterogeneity and CD112 expression are regulated by INKA1 through direct interaction with PAK4 and SIRT1, inducing epigenetic changes and defining an alternative state of LT-HSC quiescence that serves to preserve self-renewal and regenerative capacity upon regeneration-mediated stress. Collectively, our data uncovered the molecular intricacies underlying HSC heterogeneity and self-renewal regulation and point to latency as an orchestrated physiological response that balances blood cell demands with preserving a stem cell reservoir. Dick and colleagues identify human LT-HSC subsets with distinct quiescent states. They link these differences to INKA1-mediated downregulation of the transmembrane protein CD112 and its interaction with the protein deacetylase SIRT1. INKA1 is inversely correlated with the histone H4K16Ac mark, which then distinguishes ‘latent’ CD112lo LT-HSCs from CD112hi LT-HSCs that are more readily activated in response to hematopoietic stress.

T cells are able to recognize and eliminate a wide variety of immunologic threats while maintaining self-tolerance. Pathogen recognition and clearance is ensured by a process called V(D)J recombination, during which a T cell obtains a unique set of V, D and J gene segments for all the chains (α and β or γ and δ) that make up its T cell receptor (TCR). Each recombined TCR detects a specific disease-associated antigen peptide, which triggers the appropriate adaptive immune response. Understanding the relationship between TCR sequences (clonotypes) and T cell activation during disease pathogenesis and progression can assist in the development of next-generation therapeutics with more favorable and sustainable outcomes. Recent advances in single-cell sequencing allowed for simultaneous profiling of TCRs and full transcriptomes leading to the characterization of key T cell populations with pathogen recognition and disease clearance capabilities. Despite their success, these methods rely on microfluidics devices or plate-based protocols with limited sensitivity and throughput (1,000s-10,000s of cells) making the study of disease-relevant T cells time consuming and costly. To overcome these limitations, we have extended our split pool combinatorial barcoding technology to simultaneously characterize the TCRs alongside the full transcriptomes of up to 1 million T cells. As a proof of concept experiment, we applied our method to primary T cells from 8 healthy human donors to characterize their TCR diversity. From the 920,000 T cells assayed we were able to recover TCR sequences from at least one chain in 88% of T cells (807,000 cells) and detected hundreds of thousands of unique alpha and beta chains across all donors. Furthermore we found that most of these chains represent rare clonotypes found in 1 or 2 cells. Lastly, we performed transcriptome-based clustering analyses to identify all major T cell subsets and discovered donor and subset-specific hyper-expanded clonotypes that could help us reconstruct each donor’s recent history of infection. In summary, we report on the extension of our highly flexible and scalable combinatorial barcoding technology to allow researchers to profile up to 1 million of T cells in a single experiment and investigate their functional responses during infection, cancer, autoimmunity, or therapeutic interventions. Citation Format: Efthymia Papalexi, Grace Kim, Bryan Hariadi, Sarah Schroeder, Peter Matulich, Vuong Tran, Daniel Diaz, Charles Roco, Alexander Rosenberg. Using combinatorial barcoding to simultaneously profile the transcriptome and immune repertoire of 1 million T cells [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 248.

In the thirteen years since its inception, single-cell RNA-sequencing (scRNA-seq) has rapidly spread across multiple fields of research, leading to many new discoveries. As technologies have matured, the number of cells that can be processed in a single experiment has seen exponential growth with workflows now assaying up to one million cells in an individual experiment. While high throughput sequencing methods have facilitated the discovery and characterization of various cell types, sequencing costs can be prohibitively high for routine use. Many applications of scRNA-seq are focused on cell type identification, gene regulatory networks, or biomarker discovery. These applications often do not require surveying the entire transcriptome, but rather require the interrogation of specific sets of well-characterized genes. In these cases, sequencing the entire transcriptome may be adding unnecessary project costs. To increase throughput and minimize sequencing costs, the development of a targeted gene enrichment method is required. Here we extend our whole transcriptome split-pool combinatorial barcoding technology, to enable enrichment of a subset of genes in the final single cell sequencing libraries. To illustrate the power of our technology, we enriched a whole transcriptome library of human bone marrow mononuclear cells (BMMCs) from four acute myeloid leukemia (AML) donors and four control donors. From a library of 27,000 cells, we used our immune gene panel to enrich 1,000 genes representing canonical immune cell markers and pathways. Our method increased the percent of reads on target from as low as 1.6% in the whole transcriptome libraries to 75% in the targeted libraries. Furthermore, despite a nearly ten-fold reduction in sequencing reads between unenriched and enriched libraries, the resulting clustering yielded very high concordance of cell type identities and preserved AML-specific signatures such as FLT3, MKI67, and CD19. Overall, we demonstrate our modular enrichment strategy preserves biological structure of the data and allows for deep characterization of gene signatures in health and disease. We envision our approach will enable researchers to simultaneously reduce sequencing costs while drastically scaling up the number of cells and samples across experiments. Citation Format: Grace Kim, Efthymia Papalexi, Peter Matulich, Ryan Koehler, Sarah Schroeder, Vuong Tran, Charles Roco, Alexander Rosenberg. Targeted transcriptome sequencing enables exponential scaling of combinatorial barcoding in AML samples [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 226.

The expression of inhibitory immune checkpoint molecules, such as programmed death-ligand (PD-L)1, is frequently observed in human cancers and can lead to the suppression of T cell-mediated immune responses. Here, we apply expanded CRISPR-compatible (EC)CITE-seq, a technology that combines pooled CRISPR screens with single-cell mRNA and surface protein measurements, to explore the molecular networks that regulate PD-L1 expression. We also develop a computational framework, mixscape, that substantially improves the signal-to-noise ratio in single-cell perturbation screens by identifying and removing confounding sources of variation. Applying these tools, we identify and validate regulators of PD-L1 and leverage our multimodal data to identify both transcriptional and post-transcriptional modes of regulation. Specifically, we discover that the Kelch-like protein KEAP1 and the transcriptional activator NRF2 mediate the upregulation of PD-L1 after interferon (IFN)-γ stimulation. Our results identify a new mechanism for the regulation of immune checkpoints and present a powerful analytical framework for the analysis of multimodal single-cell perturbation screens.

Summary The simultaneous measurement of multiple modalities represents an exciting frontier for single-cell genomics and necessitates computational methods that can define cellular states based on multimodal data. Here, we introduce “weighted-nearest neighbor” analysis, an unsupervised framework to learn the relative utility of each data type in each cell, enabling an integrative analysis of multiple modalities. We apply our procedure to a CITE-seq dataset of 211,000 human peripheral blood mononuclear cells (PBMCs) with panels extending to 228 antibodies to construct a multimodal reference atlas of the circulating immune system. Multimodal analysis substantially improves our ability to resolve cell states, allowing us to identify and validate previously unreported lymphoid subpopulations. Moreover, we demonstrate how to leverage this reference to rapidly map new datasets and to interpret immune responses to vaccination and coronavirus disease 2019 (COVID-19). Our approach represents a broadly applicable strategy to analyze single-cell multimodal datasets and to look beyond the transcriptome toward a unified and multimodal definition of cellular identity., Graphical abstract, Highlights • “Weighted nearest neighbor” analysis integrates multimodal single-cell data • A multimodal reference “atlas” of the circulating human immune system • Identification and validation of novel sources of lymphoid heterogeneity • “Reference-based” mapping of query datasets onto a multimodal atlas, A framework that allows for the integration of multiple data types using single cells is applied to understand distinct immune cell states, previously unidentified immune populations, and to interpret immune responses to vaccinations.

Recent technological advances have enabled massively parallel chromatin profiling with scATAC-seq (single-cell assay for transposase accessible chromatin by sequencing). Here we present ATAC with select antigen profiling by sequencing (ASAP-seq), a tool to simultaneously profile accessible chromatin and protein levels. Our approach pairs sparse scATAC-seq data with robust detection of hundreds of cell surface and intracellular protein markers and optional capture of mitochondrial DNA for clonal tracking, capturing three distinct modalities in single cells. ASAP-seq uses a bridging approach that repurposes antibody:oligonucleotide conjugates designed for existing technologies that pair protein measurements with single-cell RNA sequencing. Together with DOGMA-seq, an adaptation of CITE-seq (cellular indexing of transcriptomes and epitopes by sequencing) for measuring gene activity across the central dogma of gene regulation, we demonstrate the utility of systematic multi-omic profiling by revealing coordinated and distinct changes in chromatin, RNA and surface proteins during native hematopoietic differentiation and peripheral blood mononuclear cell stimulation and as a combinatorial decoder and reporter of multiplexed perturbations in primary T cells.

The expression of inhibitory immune checkpoint molecules such as PD-L1 is frequently observed in human cancers and can lead to the suppression of T cell-mediated immune responses. Here we apply ECCITE-seq, a technology which combines pooled CRISPR screens with single-cell mRNA and surface protein measurements, to explore the molecular networks that regulate PD-L1 expression. We also develop a computational framework, mixscape, that substantially improves the signal-to-noise ratio in single-cell perturbation screens by identifying and removing confounding sources of variation. Applying these tools, we identify and validate regulators of PD-L1, and leverage our multi-modal data to identify both transcriptional and post-transcriptional modes of regulation. In particular, we discover that the kelch-like protein KEAP1 and the transcriptional activator NRF2, mediate levels of PD-L1 upregulation after IFNγ stimulation. Our results identify a novel mechanism for the regulation of immune checkpoints and present a powerful analytical framework for the analysis of multi-modal single-cell perturbation screens.

The process whereby hematopoietic stem cells (HSCs) generate different blood cell types in the steady-state is poorly understood. Upadhaya et al. used inducible lineage tracing to characterize the earliest steps of adult HSC differentiation in vivo., Adult hematopoiesis has been studied in terms of progenitor differentiation potentials, whereas its kinetics in vivo is poorly understood. We combined inducible lineage tracing of endogenous adult hematopoietic stem cells (HSCs) with flow cytometry and single-cell RNA sequencing to characterize early steps of hematopoietic differentiation in the steady-state. Labeled cells, comprising primarily long-term HSCs and some short-term HSCs, produced megakaryocytic lineage progeny within 1 wk in a process that required only two to three cell divisions. Erythroid and myeloid progeny emerged simultaneously by 2 wk and included a progenitor population with expression features of both lineages. Myeloid progenitors at this stage showed diversification into granulocytic, monocytic, and dendritic cell types, and rare intermediate cell states could be detected. In contrast, lymphoid differentiation was virtually absent within the first 3 wk of tracing. These results show that continuous differentiation of HSCs rapidly produces major hematopoietic lineages and cell types and reveal fundamental kinetic differences between megakaryocytic, erythroid, myeloid, and lymphoid differentiation., Graphical Abstract

Patients with myeloproliferative neoplasms (MPNs) frequently progress to bone marrow failure or acute myeloid leukemia (AML), and mutations in epigenetic regulators such as the metabolic enzyme isocitrate dehydrogenase (IDH) are associated with poor outcomes. Here, we showed that combined expression of Jak2V617Fand mutant IDH1R132Hor Idh2R140Q induces MPN progression, alters stem/progenitor cell function, and impairs differentiation in mice. Jak2V617FIdh2R140Q–mutant MPNs were sensitive to small-molecule inhibition of IDH. Combined inhibition of JAK2 and IDH2 normalized the stem and progenitor cell compartments in the murine model and reduced disease burden to a greater extent than was seen with JAK inhibition alone. In addition, combined JAK2 and IDH2 inhibitor treatment also reversed aberrant gene expression in MPN stem cells and reversed the metabolite perturbations induced by concurrent JAK2 and IDH2 mutations. Combined JAK2 and IDH2 inhibitor therapy also showed cooperative efficacy in cells from MPN patients with both JAK2mutand IDH2mutmutations. Taken together, these data suggest that combined JAK and IDH inhibition May offer a therapeutic advantage in this high-risk MPN subtype., Damon Runyon Cancer Research Foundation (DRG-2241-15), Howard Hughes Medical Institute (Faculty Scholars Award), Stand Up To Cancer, National Cancer Institute (U.S.) (P50CA165962), National Cancer Institute (U.S.) (P30CA14051), Koch Institute for Integrative Cancer Research ( Dana-Farber Harvard Cancer Center Bridge Project), Leukemia & Lymphoma Society of America. Specialized Center of Research (SCOR) Program, National Institutes of Health (U.S.) (grant U54OD020355-01), National Institutes of Health (U.S.) (grant NCI R01CA172636), National Institutes of Health (U.S.) (grant R35CA197594), National Cancer Institute (U.S.) (Cancer Center Support Grant (P30 CA008747).

Advances in single-cell RNA sequencing (scRNA-seq) have allowed for comprehensive analysis of the immune system. In this Review, we briefly describe the available scRNA-seq technologies together with their corresponding strengths and weaknesses. We discuss in depth how scRNA-seq can be used to deconvolve immune system heterogeneity by identifying novel distinct immune cell subsets in health and disease, characterizing stochastic heterogeneity within a cell population and building developmental 'trajectories' for immune cells. Finally, we discuss future directions of the field and present integrated approaches to complement molecular information from a single cell with studies of the environment, epigenetic state and cell lineage.

Summary Genetic and functional studies underscore the central role of JAK/STAT signaling in myeloproliferative neoplasms (MPNs). However, the mechanisms that mediate transformation in MPNs are not fully delineated, and clinically utilized JAK inhibitors have limited ability to reduce disease burden or reverse myelofibrosis. Here we show that MPN progenitor cells are characterized by marked alterations in gene regulation through differential enhancer utilization, and identify nuclear factor κB (NF-κB) signaling as a key pathway activated in malignant and non-malignant cells in MPN. Inhibition of BET bromodomain proteins attenuated NF-κB signaling and reduced cytokine production in vivo . Most importantly, combined JAK/BET inhibition resulted in a marked reduction in the serum levels of inflammatory cytokines, reduced disease burden, and reversed bone marrow fibrosis in vivo .

The first insight into the complexity of post-transplant in vivo dynamics of hematopoietic stem cells (HSC) in humans was only recently showcased by examining longitudinal clonal contributions in gene therapy patients. Initial blood reconstitution was achieved by short-term (ST-)HSC, but in the longer-term hematopoiesis originated from long-term (LT-)HSC that only became recruited after a latency phase of 1-2 years (Scala et al 2018). Thus, deciphering the mechanisms governing how LT-HSC might resist transplant mediated activation and/or respond to the varying hematopoietic demands that occur during homeostasis would be an important goal for improving HSC cell therapies. We previously linked INKA1(=C3orf54/FAM212A) mediated PAK4 inhibition and reduced H4K16 acetylation (H4K16ac) to quiescence in human leukemia stem cells (Kaufmann, Blood 2019). Here, we interrogate the role of this signaling axis in normal human hematopoiesis. Immunostaining revealed distinct subsets defined by the dichotomy of INKA1 and H4K16ac within cord blood ST- and LT-HSC, mechanistically supported by BioID, chromatin fiber analysis and PLA data showing INKA1 interacting with the H4K16 deacetylase, SIRT1. Among quiescent LT-HSC, we found that the INKA1high fraction (~10 % of LT-HSC) had the lowest CDK6 levels and represented LT-HSC in an alternative state of quiescence. We then used protein interaction (BioID) data to show that a shared interactor of INKA1 and PAK4, CD112 could be used to sort for the subset of CD112low LT-HSC that was in this alternatively quiescent LT-HSC state (H4K16aclow, CDK6low, CellROXlow). Compared to primary cells, culture attenuates the strict dichotomy of INKA1 and H4K16ac but only in cells in which colocalization of PAK4 with both occurred, suggesting PAK4 interferes with SIRT1-INKA1 interaction thereby permitting H4K16 acetylation. In vitro time course analysis showed that H4K16ac and PAK4 levels preceded the upregulation of the G0-exit marker CDK6 implicating a role in cell cycle priming. Pseudo-time scRNAseq analysis for INKA1high versus PAK4high, CDK6high and CD112high expression in mobilized peripheral blood (modeling in vivo HSC activation) showed enrichment of early or late diffusion indicative of quiescent versus primed cell status, respectively. Strikingly, in xenograft assays CD112low LT-HSC exhibited delayed engraftment kinetics with higher secondary repopulation capacity than faster repopulating CD112high LT-HSC reflecting the subset of LT-HSC that resist early activation. Similarly, INKA1-OE or PAK4 knock-down in vivo resulted in an early restraint in engraftment levels (@4 w), whereas 20 w engraftment reached control levels. When measured in secondary transplant assays the HSC frequency was 4 to 8-fold higher in the groups that showed this early restraint. Thus, resisting early activation (latency) preserves the regenerative potential of such HSC. In vivo 5-Fluorouracil treatment at week 4 accelerated latency exit and further increased HSC frequency of INKA1-OE cells while abolishing serial transplantation potential of controls. Hence, the induction of proliferative stress via creating acute hematopoietic demand is able to lift the INKA1-OE imposed restraint resulting in increased hematopoietic output while also promoting HSC self-renewal. Collectively, our data decipher the molecular intricacies underlying HSC heterogeneity and self-renewal regulation and point to latency as an orchestrated physiological response that integrates quiescence control with HSC fate choices to preserve a stem cell reservoir. Disclosures Takayanagi: Kirin Holdings Company, Ltd: Current Employment. Dick:Bristol-Myers Squibb/Celgene: Research Funding.

The mechanisms governing how human hematopoietic stem cells (HSC) balance their capacity for maintaining long-term (LT) regeneration with meeting daily and/or stress demands are unknown. We previously linked INKA1-mediated PAK4 inhibition along with reduction of H4K16 acetylation (H4K16ac) to quiescence in leukemia stem cells. Here, we uncover mutually exclusive immunostaining for INKA1 and H4K16ac in cord blood HSC; supported by BioID and PLA data showing INKA1 interacts with the H4K16 deacetylase SIRT1. The proteomic data also pointed us to CD112, that enabled isolation of LT-HSC into transcriptionally distinct subsets that respond differently to regenerative stress: a CD112low subset (H4K16aclow, CDK6low, ROSlow, INKA1high) exhibited a transient restraint before contributing to hematopoietic reconstitution and a primed CD112high subset (H4K16achigh, CDK6high, ROShigh, INKA1low) that responded rapidly. Protein and scRNAseq pseudotime analysis of in vitro and in vivo (G-CSF-) activated HSC, respectively, confirmed upregulation of H4K16ac, PAK4, CDK6 and CD112 coinciding with cell-cycle priming. INKA1-overexpression (OE) or PAK4 knock-down (KD) induced transient restraint in reconstitution (@4w) with 20w engraftment reaching control levels; HSC frequency increased 4-8 fold as measures in secondary transplants. Thus, INKA1-OE or PAK4-KD enables HSC to resist early activation, thereby protracting their regenerative potential (a phenomenon we term latency). By elevating hematopoietic demand with 5-FU, the INKA1-OE imposed restraint was released engendering even further increased HSC frequency; in controls regenerative capacity was abolished. Collectively, our data uncover a new molecular axis for HSC self-renewal regulation and point to latency as an orchestrated physiological response that integrates quiescence control with HSC fate choices to maintain LT-HSC pool size.

Rapid technological progress in the recent years has allowed the high-throughput interrogation of different types of biomolecules from single cells. Combining several of these readouts into integrated multi-omic assays is essential to comprehensively understand and model cellular processes. Here, we report the development of Expanded CRISPR-compatible Cellular Indexing of Transcriptomes and Epitopes by sequencing (ECCITE-seq) for the high-throughput characterization of at least five modalities of information from each single cell: transcriptome, immune receptor clonotypes, surface markers, sample identity and sgRNAs. We demonstrate the use of ECCITE-seq to directly and efficiently capture sgRNA molecules and measure their effects on gene expression and protein levels, opening the possibility of performing high throughput single cell CRISPR screens with multimodal readout using existing libraries and commonly used vectors. Finally, by utilizing the combined phenotyping of clonotype and cell surface markers in immune cells, we apply ECCITE to study a lymphoma sample to discriminate cells and define molecular signatures of malignant cells within a heterogeneous population.

Single-cell transcriptomics has transformed our ability to characterize cell states, but deep biological understanding requires more than a taxonomic listing of clusters. As new methods arise to measure distinct cellular modalities, a key analytical challenge is to integrate these datasets to better understand cellular identity and function. Here, we develop a strategy to "anchor" diverse datasets together, enabling us to integrate single-cell measurements not only across scRNA-seq technologies, but also across different modalities. After demonstrating improvement over existing methods for integrating scRNA-seq data, we anchor scRNA-seq experiments with scATAC-seq to explore chromatin differences in closely related interneuron subsets and project protein expression measurements onto a bone marrow atlas to characterize lymphocyte populations. Lastly, we harmonize in situ gene expression and scRNA-seq datasets, allowing transcriptome-wide imputation of spatial gene expression patterns. Our work presents a strategy for the assembly of harmonized references and transfer of information across datasets.

Single cell transcriptomics (scRNA-seq) has transformed our ability to discover and annotate cell types and states, but deep biological understanding requires more than a taxonomic listing of clusters. As new methods arise to measure distinct cellular modalities, including high-dimensional immunophenotypes, chromatin accessibility, and spatial positioning, a key analytical challenge is to integrate these datasets into a harmonized atlas that can be used to better understand cellular identity and function. Here, we develop a computational strategy to “anchor” diverse datasets together, enabling us to integrate and compare single cell measurements not only across scRNA-seq technologies, but different modalities as well. After demonstrating substantial improvement over existing methods for data integration, we anchor scRNA-seq experiments with scATAC-seq datasets to explore chromatin differences in closely related interneuron subsets, and project single cell protein measurements onto a human bone marrow atlas to annotate and characterize lymphocyte populations. Lastly, we demonstrate how anchoring can harmonize in-situ gene expression and scRNA-seq datasets, allowing for the transcriptome-wide imputation of spatial gene expression patterns, and the identification of spatial relationships between mapped cell types in the visual cortex. Our work presents a strategy for comprehensive integration of single cell data, including the assembly of harmonized references, and the transfer of information across datasets.Availability: Installation instructions, documentation, and tutorials are available at: https://www.satijalab.org/seurat

Multimodal single-cell assays provide high-resolution snapshots of complex cell populations, but are mostly limited to transcriptome plus an additional modality. Here, we describe expanded CRISPR-compatible cellular indexing of transcriptomes and epitopes by sequencing (ECCITE-seq) for the high-throughput characterization of at least five modalities of information from each single cell. We demonstrate application of ECCITE-seq to multimodal CRISPR screens with robust direct single-guide RNA capture and to clonotype-aware multimodal phenotyping of cancer samples.

Computational single-cell RNA-seq (scRNA-seq) methods have been successfully applied to experiments representing a single condition, technology, or species to discover and define cellular phenotypes. However, identifying subpopulations of cells that are present across multiple data sets remains challenging. Here, we introduce an analytical strategy for integrating scRNA-seq data sets based on common sources of variation, enabling the identification of shared populations across data sets and downstream comparative analysis. We apply this approach, implemented in our R toolkit Seurat (http://satijalab.org/seurat/), to align scRNA-seq data sets of peripheral blood mononuclear cells under resting and stimulated conditions, hematopoietic progenitors sequenced using two profiling technologies, and pancreatic cell 'atlases' generated from human and mouse islets. In each case, we learn distinct or transitional cell states jointly across data sets, while boosting statistical power through integrated analysis. Our approach facilitates general comparisons of scRNA-seq data sets, potentially deepening our understanding of how distinct cell states respond to perturbation, disease, and evolution.

Hematopoietic stem cells (HSCs) give rise to diverse cell types in the blood system, yet our molecular understanding of the early trajectories that generate this enormous diversity in humans remains incomplete. Here, we leverage Drop‐seq, a massively parallel single‐cell RNA sequencing (scRNA‐seq) approach, to individually profile 20,000 progenitor cells from human cord blood, without prior enrichment or depletion for individual lineages based on surface markers. Our data reveal a transcriptional compendium of progenitor states in human cord blood, representing four committed lineages downstream from HSC, alongside the transcriptional dynamics underlying fate commitment. We identify intermediate stages that simultaneously co‐express “primed” programs for multiple downstream lineages, and also observe striking heterogeneity in the early molecular transitions between myeloid subsets. Integrating our data with a recently published scRNA‐seq dataset from human bone marrow, we illustrate the molecular similarity between these two commonly used systems and further explore the chromatin dynamics of “primed” transcriptional programs based on ATAC‐seq. Finally, we demonstrate that Drop‐seq data can be utilized to identify new heterogeneous surface markers of cell state that correlate with functional output.

Viny et al. show that depletion of the cohesin subunit Smc3 affects HSC function through disruption of chromatin structure and that cohesin haploinsufficiency cooperates with Flt3-ITD to induce transformation of stem/progenitor cells and leukemia development., Cohesin complex members have recently been identified as putative tumor suppressors in hematologic and epithelial malignancies. The cohesin complex guides chromosome segregation; however, cohesin mutant leukemias do not show genomic instability. We hypothesized that reduced cohesin function alters chromatin structure and disrupts cis-regulatory architecture of hematopoietic progenitors. We investigated the consequences of Smc3 deletion in normal and malignant hematopoiesis. Biallelic Smc3 loss induced bone marrow aplasia with premature sister chromatid separation and revealed an absolute requirement for cohesin in hematopoietic stem cell (HSC) function. In contrast, Smc3 haploinsufficiency increased self-renewal in vitro and in vivo, including competitive transplantation. Smc3 haploinsufficiency reduced coordinated transcriptional output, including reduced expression of transcription factors and other genes associated with lineage commitment. Smc3 haploinsufficiency cooperated with Flt3-ITD to induce acute leukemia in vivo, with potentiated Stat5 signaling and altered nucleolar topology. These data establish a dose dependency for cohesin in regulating chromatin structure and HSC function.

SummaryAlthough clinically tested JAK inhibitors reduce splenomegaly and systemic symptoms, molecular responses are not observed in most myeloproliferative neoplasm (MPN) patients. We previously demonstrated that MPN cells become persistent to type I JAK inhibitors that bind the active conformation of JAK2. We investigated whether CHZ868, a type II JAK inhibitor, would demonstrate activity in JAK inhibitor persistent cells, murine MPN models, and MPN patient samples. JAK2 and MPL mutant cell lines were sensitive to CHZ868, including type I JAK inhibitor persistent cells. CHZ868 showed significant activity in murine MPN models and induced reductions in mutant allele burden not observed with type I JAK inhibitors. These data demonstrate that type II JAK inhibition is a viable therapeutic approach for MPN patients.

The identification of JAK2/MPL mutations in patients with myeloproliferative neoplasms (MPN) has led to the clinical development of JAK kinase inhibitors, including ruxolitinib. Ruxolitinib reduces splenomegaly and systemic symptoms in myelofibrosis and improves overall survival; however, the mechanism by which JAK inhibitors achieve efficacy has not been delineated. Patients with MPN present with increased levels of circulating proinflammatory cytokines, which are mitigated by JAK inhibitor therapy. We sought to elucidate mechanisms by which JAK inhibitors attenuate cytokine-mediated pathophysiology. Single-cell profiling demonstrated that hematopoietic cells from myelofibrosis models and patient samples aberrantly secrete inflammatory cytokines. Pan-hematopoietic Stat3 deletion reduced disease severity and attenuated cytokine secretion, with similar efficacy as observed with ruxolitinib therapy. In contrast, Stat3 deletion restricted to MPN cells did not reduce disease severity or cytokine production. Consistent with these observations, we found that malignant and nonmalignant cells aberrantly secrete cytokines and JAK inhibition reduces cytokine production from both populations. Significance: Our results demonstrate that JAK–STAT3-mediated cytokine production from malignant and nonmalignant cells contributes to MPN pathogenesis and that JAK inhibition in both populations is required for therapeutic efficacy. These findings provide novel insight into the mechanisms by which JAK kinase inhibition achieves therapeutic efficacy in MPNs. Cancer Discov; 5(3); 316–31. ©2015 AACR. See related commentary by Belver and Ferrando, p. 234 This article is highlighted in the In This Issue feature, p. 213

JAK1 is a critical effector of pro-inflammatory cytokine signaling and plays important roles in immune function, while abnormal JAK1 activity has been linked to immunological and neoplastic diseases. Specific functions of JAK1 in the context of hematopoiesis, and specifically within hematopoietic stem cells (HSCs), have not clearly been delineated. Here, we show that conditional Jak1 loss in HSCs reduces their self-renewal and markedly alters lymphoid/myeloid differentiation invivo. Jak1-deficient HSCs exhibit decreased competitiveness invivo and are unable to rescue hematopoiesis in the setting of myelosuppression. They exhibit increased quiescence, an inability to enter the cell cycle in response to hematopoietic stress, and a marked reduction in cytokine sensing, including in response to type I interferons and IL-3. Moreover, Jak1 loss is not fully rescued byexpression of a constitutively active Jak2 allele. Together, these data highlight an essential role for Jak1 in HSC homeostasis and stress responses.

Reply to "Uveal melanoma cells are resistant to EZH2 inhibition regardless of BAP1 status"

The development of the dual Janus kinase 1/2 (JAK1/2) inhibitor ruxolitinib for the treatment of myeloproliferative neoplasms (MPNs) has led to studies of ruxolitinib in other clinical contexts, including JAK-mutated acute lymphoblastic leukemia (ALL). However, the limited ability of JAK inhibition to induce molecular or clinicopathological responses in MPNs suggests a need for development of better therapies for JAK kinase-dependent malignancies. Here, we demonstrate that heat shock protein 90 (HSP90) inhibition using a purine-scaffold HSP90 inhibitor in early clinical development is an effective therapeutic approach in JAK-dependent ALL and can overcome persistence to JAK-inhibitor therapy in ALL cells.

The tumor suppressors BAP1 and ASXL1 interact to form a polycomb deubiquitinase complex that removes monoubiquitin from histone H2A lysine 119 (H2AK119Ub). However, BAP1 and ASXL1 are mutated in distinct cancer types, consistent with independent roles in regulating epigenetic state and malignant transformation. Here we demonstrate that Bap1 loss in mice results in increased trimethylated histone H3 lysine 27 (H3K27me3), elevated enhancer of zeste 2 polycomb repressive complex 2 subunit (Ezh2) expression, and enhanced repression of polycomb repressive complex 2 (PRC2) targets. These findings contrast with the reduction in H3K27me3 levels seen with Asxl1 loss. Conditional deletion of Bap1 and Ezh2 in vivo abrogates the myeloid progenitor expansion induced by Bap1 loss alone. Loss of BAP1 results in a marked decrease in H4K20 monomethylation (H4K20me1). Consistent with a role for H4K20me1 in the transcriptional regulation of EZH2, expression of SETD8-the H4K20me1 methyltransferase-reduces EZH2 expression and abrogates the proliferation of BAP1-mutant cells. Furthermore, mesothelioma cells that lack BAP1 are sensitive to EZH2 pharmacologic inhibition, suggesting a novel therapeutic approach for BAP1-mutant malignancies.

Cytokine-mediated signal transduction is critical to hematopoiesis, immune responses, and other physiological processes. Aberrant production and secretion of pro-inflammatory cytokines disturbs homeostasis and proper immune function and if persistent results in symptoms of chronic inflammation. Previous studies have illustrated the importance of JAK1 as an effector of cytokine signaling, including in immunological and neoplastic diseases such that selective JAK1 inhibition is currently being investigated in clinical trials. However, the role of Jak1 in hematopoietic stem cell (HSC) function has not been delineated. This has led us to investigate the impact of loss of Jak1 signaling on HSC function by developing a novel conditional Jak1 knockout allele (Fig. 1a). Mice with conditional deletion of Jak1 in the hematopoietic system (hereafter referred to as Jak1 KO) are characterized by leukocytosis (Jak1 KO avg. 6.34K/ul, Jak1 WT avg. 10.76K/ul, P Consistent with the stem cell phenotype observed in JAK1 KO mice, we found that a significant larger proportion of Jak1-deficient stem cells lacks expression of the proliferation marker Ki67 and that Jak1-deficient stem cells fail to enter the cell cycle in response to hematopoietic stress. To begin to determine the mechanism by which Jak1 regulates normal stem cell function in vivo, we assessed the impact of loss of Jak1 on transcriptional output. Gene expression profiling of LT-HSCs from Jak1 KO and WT mice identified 259 significant genes, many of which were known to be Jak1 downstream targets. Gene set enrichment analysis (GSEA) revealed that the majority of genes that were altered following deletion of Jak1 corresponded to interferon signaling and inflammatory response pathways. Consistent with these findings, our functional in vitro and in vivo assays demonstrated that Jak1-deficient cells were insensitive to type I interferons as shown by lack of Stat1 and Stat5 activation (Fig. 1e), retained Sca1 surface expression, and an unchanged cell cycle status upon IFN stimulation. Moreover, the HSC defect observed in the setting of Jak1 loss was not fully rescued by expression of a constitutively active Jak2 allele, suggesting there is non-redundant signaling in HSCs within the JAK kinase family. Together, our data suggests that Jak1 functions as a central node for interferon signaling in HSCs and reveals an essential and nonredundant role of Jak1 in HSC homeostasis and stress response. Figure 1 a) Design of a conditional targeting vector and confirmation of gene deletion on protein level. b) Reduction of LSK cells in Jak1 KO mice. c) Competitive disadvantage of Jak1-deficient cells. d) Increased mortality of Jak1 KO mice when serially challenged with 5-FU. e) Jak1-deficient LSK cells are insensitive to type I interferon stimulation. Figure 1. a) Design of a conditional targeting vector and confirmation of gene deletion on protein level. b) Reduction of LSK cells in Jak1 KO mice. c) Competitive disadvantage of Jak1-deficient cells. d) Increased mortality of Jak1 KO mice when serially challenged with 5-FU. e) Jak1-deficient LSK cells are insensitive to type I interferon stimulation. Disclosures Koppikar: Amgen: Employment. Nolan:Fluidigm: Consultancy. Levine:Novartis: Consultancy; Qiagen: Membership on an entity's Board of Directors or advisory committees.

Cohesin complex members have recently been identified as putative tumor suppressors in hematologic and epithelial malignancies. The cohesin complex guides chromosome segregation, however cohesin-mutant leukemias do not show genomic instability suggesting an alternate role in malignant transformation. We hypothesized reduced cohesin function alters chromatin structure and disrupts cis-regulatory architecture of hematopoietic stem/progenitor cells. We therefore investigated the impact of both complete loss and haploinsufficiency of Smc3, an obligate member of the cohesin complex, in normal hematopoiesis and in myeloid transformation by developing a conditional Smc3 knockout allele. Somatic loss of Smc3 in hematopoietic cells induced lethal bone marrow aplasia (median survival 11 days; p To further explore the mechanism by which Smc3 loss cooperates with Flt3ITD to induce leukemia, we investigated chromatin cis-regulatory architecture with transposase hypersensitivity assays (ATAC-seq). We hypothesized that increased accessibility at cis-regulatory elements and the alterations in gene expression seen in cells with combined Smc3 haploinsufficiency and Flt3ITD may be in a large part driven by potentiated Stat signaling at chromatin. We analyzed 146 transcription factor recognition motifs within the THS differentially observed in Smc3Δ/+Flt3ITD and wild-type cells. Chromatin accessibility gained in Smc3Δ/+Flt3ITD cells are enriched in Stat family transcription factor binding sites, including Stat5. We also observed enrichment of the Stat5 gene expression signature in the Smc3Δ/+Flt3ITD cells compared to Smc3Δ/+, Flt3ITD and wild-type cells, suggesting the divergent mutations cooperate to potentiate oncogenic Stat5 signaling in HSPCs. Our results demonstrate a key dose-dependent role for the cohesin complex in hematopoiesis, and show that reduced cohesin functions to alter enhancer-mediated transcription and contribute to aberrant self-renewal and myeloid transformation. Figure 1. Figure 1. Disclosures Levine: Loxo Oncology: Membership on an entity's Board of Directors or advisory committees; CTI BioPharma: Membership on an entity's Board of Directors or advisory committees; Foundation Medicine: Consultancy.

Recurrent somatic loss-of-function mutations in ASXL1 (Addition of sex combs-like 1) are common genetic events in a spectrum of myeloid malignancies and these alterations demarcate patients with poor outcome. ASXL1 forms a chromatin regulatory complex with the ubiquitin hydrolase BAP1 (BRCA1 associated protein-1), a protein that has been found to be transcriptionally repressed in MDS patients. These data are consistent with BAP1 having tumor suppressive activity in MDS; however, the mechanism by which disruption of the ASXL1-BAP1 axis leads to transformation is not well understood. We conditionally deleted Bap1 in the murine hematopoietic system utilizing Mx1-Cre (hereafter referred to as Bap1 KO). One hundred percent of mice with confirmed Bap1 deletion developed a fully penetrant MDS-like disease characterized by leukocytosis, anemia, and splenomegaly. Bap1 KO mice have an expansion of the granulocyte macrophage progenitor compartment (GMP; Lin- c-Kit+ Sca1- CD34+ Fcϒ+). Given the role of BAP1 in epigenetic regulation, we investigated the effect of Bap1 loss on chromatin state and transcriptional output. We first assessed epigenetic changes in Bap1 KO mice by performing histone mass spectometry in control and Bap1 KO hematopoietic stem and progenitor cells (HSPCs, c-Kit+ enriched). Bap1 loss increased H3K27me2/3 at the expense of H3K27me0/1. We confirmed that H3K27me3 was increased in Bap1 KO bone marrow cells by completing H3K27me3 ChIP-Sequencing in HSPCs. Enumeration of H3K27me3 peaks in Bap1 KO versus control cells indicated an increase in H3K27me3 domains (Figure A). We next overlaid RNA-Sequencing from GMP sorted Bap1 KO bone marrow cells with genes marked by H3K27me3, as indicated by ChIP-Sequencing. We found that Bap1 loss resulted in a global decrease in gene expression (68% downregulated, 657/968 genes, p-adj Given the alterations in H3K27me3 in Bap1 KO mice, we investigated whether Bap1- deficient transformation could be rescued by abrogation of PRC2-mediated gene repression. We developed a genetic model with compound deletion of Bap1 and Ezh2, the catalytic component of the PRC2 complex. Co-deletion of Bap1 and Ezh2 resulted in a phenotypic rescue of Bap1 KO associated splenomegaly (spleen weights, Bap1 KO avg. 541.6 mg, Bap1/Ezh2 KO avg. 157.0 mg, p EZH2 small molecule inhibitors have proven effective in EZH2-dependent models of B cell lymphoma. To determine if Ezh2 inhibition was efficacious in the setting of Bap1 loss, we treated a cohort of Bap1 KO mice with either vehicle (NaCMC) or 500 mg/kg EPZ011989, an EZH2 inhibitor with in vivo activity. Treatment of Bap1 KO mice for 16 days resulted in significant reduction of splenomegaly (spleen weights, vehicle avg. 522.0, EPZ011989 treated avg. 216.2, p Figure 1. Figure 1. Disclosures Knutson: Epizyme, Inc: Employment. Campbell:Epizyme, Inc: Employment. Keilhack:Epizyme: Employment, Equity Ownership. Melnick:Janssen: Other: Research; ROCHE: Other: Research; Genentech: Speakers Bureau; Celgene: Consultancy; Eli Lilly: Consultancy; Epizyme: Consultancy. Armstrong:Epizyme, Inc: Consultancy. Levine:Foundation Medicine: Consultancy; CTI BioPharma: Membership on an entity's Board of Directors or advisory committees; Loxo Oncology: Membership on an entity's Board of Directors or advisory committees.