A single cell framework identifies functionally and molecularly distinct multipotent progenitors in adult human hematopoiesis.

Hematopoietic multipotent progenitors (MPPs) regulate blood cell production to appropriately meet the biological demands of the human body. Human MPPs remain ill-defined whereas mouse MPPs have been well characterized with distinct immunophenotypes and lineage potencies. Using multiomic single cell analyses and complementary functional assays, we identified new human MPPs and oligopotent progenitor populations within Lin-CD34+CD38dim/lo adult bone marrow with distinct biomolecular and functional properties. These populations were prospectively isolated based on expression of CD69, CLL1, and CD2 in addition to classical markers like CD90 and CD45RA. We show that within the canonical Lin-CD34+CD38dim/loCD90CD45RA-MPP population, there is a CD69+ MPP with long-term engraftment and multilineage differentiation potential, a CLL1+ myeloid-biased MPP, and a CLL1-CD69-erythroid-biased MPP. We also show that the canonical Lin-CD34+CD38dim/loCD90-CD45RA+ LMPP population can be separated into a CD2+ LMPP with lymphoid and myeloid potential, a CD2-LMPP with high lymphoid potential, and a CLL1+ GMP with minimal lymphoid potential. We used these new HSPC profiles to study human and mouse bone marrow cells and observe limited cell type specific homology between humans and mice and cell type specific changes associated with aging. By identifying and functionally characterizing new adult MPP sub-populations, we provide an updated reference and framework for future studies in human hematopoiesis.

AML/T cell interactomics uncover correlates of patient outcomes and the key role of ICAM1 in T cell killing of AML.

T cells are important for the control of acute myeloid leukemia (AML), a common and often deadly malignancy. We observed that some AML patient samples are resistant to killing by human-engineered cytotoxic CD4+ T cells. Single-cell RNA-seq of primary AML samples and CD4+ T cells before and after their interaction uncovered transcriptional programs that correlate with AML sensitivity or resistance to CD4+ T cell killing. Resistance-associated AML programs were enriched in AML patients with poor survival, and killing-resistant AML cells did not engage T cells in vitro. Killing-sensitive AML potently activated T cells before being killed, and upregulated ICAM1, a key component of the immune synapse with T cells. Without ICAM1, killing-sensitive AML became resistant to killing by primary ex vivo-isolated CD8+ T cells in vitro, and engineered CD4+ T cells in vitro and in vivo. While AML heterogeneity implies that multiple factors may determine their sensitivity to T cell killing, these data show that ICAM1 acts as an immune trigger, allowing T cell killing, and could play a role in AML patient survival in vivo.

Sarcoma microenvironment cell states and ecosystems are associated with prognosis and predict response to immunotherapy.

Characterization of the diverse malignant and stromal cell states that make up soft tissue sarcomas and their correlation with patient outcomes has proven difficult using fixed clinical specimens. Here, we employed EcoTyper, a machine-learning framework, to identify the fundamental cell states and cellular ecosystems that make up sarcomas on a large scale using bulk transcriptomes with clinical annotations. We identified and validated 23 sarcoma-specific, transcriptionally defined cell states, many of which were highly prognostic of patient outcomes across independent datasets. We discovered three conserved cellular communities or ecotypes associated with underlying genomic alterations and distinct clinical outcomes. We show that one ecotype defined by tumor-associated macrophages and epithelial-like malignant cells predicts response to immune-checkpoint inhibition but not chemotherapy and validate our findings in an independent cohort. Our results may enable identification of patients with soft tissue sarcomas who could benefit from immunotherapy and help develop new therapeutic strategies.

AML/T cell interactomics uncover correlates of patient outcomes and the key role of ICAM1 in T cell killing of AML.

T cells are important for the control of acute myeloid leukemia (AML), a common and often deadly malignancy. We observed that some AML patient samples are resistant to killing by human engineered cytotoxic CD4 + T cells. Single-cell RNA-seq of primary AML samples and CD4 + T cells before and after their interaction uncovered transcriptional programs that correlate with AML sensitivity or resistance to CD4 + T cell killing. Resistance-associated AML programs were enriched in AML patients with poor survival, and killing-resistant AML cells did not engage T cells in vitro . Killing-sensitive AML potently activated T cells before being killed, and upregulated ICAM1 , a key component of the immune synapse with T cells. Without ICAM1, killing-sensitive AML became resistant to killing to primary ex vivo -isolated CD8 + T cells in vitro , and engineered CD4 + T cells in vitro and in vivo . Thus, ICAM1 on AML acts as an immune trigger, allowing T cell killing, and could affect AML patient survival in vivo . SIGNIFICANCE: AML is a common leukemia with sub-optimal outcomes. We show that AML transcriptional programs correlate with susceptibility to T cell killing. Killing resistance-associated AML programs are enriched in patients with poor survival. Killing-sensitive, but not resistant AML activate T cells and upregulate ICAM1 that binds to LFA-1 on T cells, allowing immune synapse formation which is critical for AML elimination.

Profiling Cellular Ecosystems at Single-Cell Resolution and at Scale with EcoTyper.

Tissues are composed of diverse cell types and cellular states that organize into distinct ecosystems with specialized functions. EcoTyper is a collection of machine learning tools for the large-scale delineation of cellular ecosystems and their constituent cell states from bulk, single-cell, and spatially resolved gene expression data. In this chapter, we provide a primer on EcoTyper and demonstrate its use for the discovery and recovery of cell states and ecosystems from healthy and diseased tissue specimens.

T cell characteristics associated with toxicity to immune checkpoint blockade in patients with melanoma.

Severe immune-related adverse events (irAEs) occur in up to 60% of patients with melanoma treated with immune checkpoint inhibitors (ICIs). However, it is unknown whether a common baseline immunological state precedes irAE development. Here we applied mass cytometry by time of flight, single-cell RNA sequencing, single-cell V(D)J sequencing, bulk RNA sequencing and bulk T cell receptor (TCR) sequencing to study peripheral blood samples from patients with melanoma treated with anti-PD-1 monotherapy or anti-PD-1 and anti-CTLA-4 combination ICIs. By analyzing 93 pre- and early on-ICI blood samples and 3 patient cohorts (n = 27, 26 and 18), we found that 2 pretreatment factors in circulation-activated CD4 memory T cell abundance and TCR diversity-are associated with severe irAE development regardless of organ system involvement. We also explored on-treatment changes in TCR clonality among patients receiving combination therapy and linked our findings to the severity and timing of irAE onset. These results demonstrate circulating T cell characteristics associated with ICI-induced toxicity, with implications for improved diagnostics and clinical management.

Atlas of clinically distinct cell states and ecosystems across human solid tumors.

Determining how cells vary with their local signaling environment and organize into distinct cellular communities is critical for understanding processes as diverse as development, aging, and cancer. Here we introduce EcoTyper, a machine learning framework for large-scale identification and validation of cell states and multicellular communities from bulk, single-cell, and spatially resolved gene expression data. When applied to 12 major cell lineages across 16 types of human carcinoma, EcoTyper identified 69 transcriptionally defined cell states. Most states were specific to neoplastic tissue, ubiquitous across tumor types, and significantly prognostic. By analyzing cell-state co-occurrence patterns, we discovered ten clinically distinct multicellular communities with unexpectedly strong conservation, including three with myeloid and stromal elements linked to adverse survival, one enriched in normal tissue, and two associated with early cancer development. This study elucidates fundamental units of cellular organization in human carcinoma and provides a framework for large-scale profiling of cellular ecosystems in any tissue.

The landscape of tumor cell states and ecosystems in diffuse large B cell lymphoma.

Biological heterogeneity in diffuse large B cell lymphoma (DLBCL) is partly driven by cell-of-origin subtypes and associated genomic lesions, but also by diverse cell types and cell states in the tumor microenvironment (TME). However, dissecting these cell states and their clinical relevance at scale remains challenging. Here, we implemented EcoTyper, a machine-learning framework integrating transcriptome deconvolution and single-cell RNA sequencing, to characterize clinically relevant DLBCL cell states and ecosystems. Using this approach, we identified five cell states of malignant B cells that vary in prognostic associations and differentiation status. We also identified striking variation in cell states for 12 other lineages comprising the TME and forming cell state interactions in stereotyped ecosystems. While cell-of-origin subtypes have distinct TME composition, DLBCL ecosystems capture clinical heterogeneity within existing subtypes and extend beyond cell-of-origin and genotypic classes. These results resolve the DLBCL microenvironment at systems-level resolution and identify opportunities for therapeutic targeting (https://ecotyper.stanford.edu/lymphoma).

Determining cell type abundance and expression from bulk tissues with digital cytometry.

Single-cell RNA-sequencing has emerged as a powerful technique for characterizing cellular heterogeneity, but it is currently impractical on large sample cohorts and cannot be applied to fixed specimens collected as part of routine clinical care. We previously developed an approach for digital cytometry, called CIBERSORT, that enables estimation of cell type abundances from bulk tissue transcriptomes. We now introduce CIBERSORTx, a machine learning method that extends this framework to infer cell-type-specific gene expression profiles without physical cell isolation. By minimizing platform-specific variation, CIBERSORTx also allows the use of single-cell RNA-sequencing data for large-scale tissue dissection. We evaluated the utility of CIBERSORTx in multiple tumor types, including melanoma, where single-cell reference profiles were used to dissect bulk clinical specimens, revealing cell-type-specific phenotypic states linked to distinct driver mutations and response to immune checkpoint blockade. We anticipate that digital cytometry will augment single-cell profiling efforts, enabling cost-effective, high-throughput tissue characterization without the need for antibodies, disaggregation or viable cells.