New frameworks for hematopoiesis derived from single-cell genomics.

ABSTRACT: Recent advancements in single-cell genomics have enriched our understanding of hematopoiesis, providing intricate details about hematopoietic stem cell biology, differentiation, and lineage commitment. Technological advancements have highlighted extensive heterogeneity of cell populations and continuity of differentiation routes. Nevertheless, intermediate "attractor" states signify structure in stem and progenitor populations that link state transition dynamics to fate potential. We discuss how innovative model systems quantify lineage bias and how stress accelerates differentiation, thereby reducing fate plasticity compared with native hematopoiesis. We conclude by offering our perspective on the current model of hematopoiesis and discuss how a more precise understanding can translate to strategies that extend healthy hematopoiesis and prevent disease.

Oncogene-induced TIM-3 ligand expression dictates susceptibility to anti-TIM-3 therapy in mice.

Leukemia relapse is a major cause of death after allogeneic hematopoietic cell transplantation (allo-HCT). We tested the potential of targeting T cell (Tc) immunoglobulin and mucin-containing molecule 3 (TIM-3) for improving graft-versus-leukemia (GVL) effects. We observed differential expression of TIM-3 ligands when hematopoietic stem cells overexpressed certain oncogenic-driver mutations. Anti-TIM-3 Ab treatment improved survival of mice bearing leukemia with oncogene-induced TIM-3 ligand expression. Conversely, leukemia cells with low ligand expression were anti-TIM-3 treatment resistant. In vitro, TIM-3 blockade or genetic deletion in CD8+ Tc enhanced Tc activation, proliferation, and IFN-γ production while enhancing GVL effects, preventing Tc exhaustion, and improving Tc cytotoxicity and glycolysis in vivo. Conversely, TIM-3 deletion in myeloid cells did not affect allogeneic Tc proliferation and activation in vitro, suggesting that anti-TIM-3 treatment-mediated GVL effects are Tc induced. In contrast to anti-programmed cell death protein 1 (anti-PD-1) and anti-cytotoxic T lymphocyte-associated protein 4 (anti-CTLA-4) treatment, anti-TIM-3-treatment did not enhance acute graft-versus-host disease (aGVHD). TIM-3 and its ligands were frequently expressed in acute myeloid leukemia (AML) cells of patients with post-allo-HCT relapse. We decipher the connections between oncogenic mutations found in AML and TIM-3 ligand expression and identify anti-TIM-3 treatment as a strategy for enhancing GVL effects via metabolic and transcriptional Tc reprogramming without exacerbation of aGVHD. Our findings support clinical testing of anti-TIM-3 Ab in patients with AML relapse after allo-HCT.

Targetable leukaemia dependency on noncanonical PI3Kγ signalling.

Phosphoinositide-3-kinase-γ (PI3Kγ) is implicated as a target to repolarize tumour-associated macrophages and promote antitumour immune responses in solid cancers1-4. However, cancer cell-intrinsic roles of PI3Kγ are unclear. Here, by integrating unbiased genome-wide CRISPR interference screening with functional analyses across acute leukaemias, we define a selective dependency on the PI3Kγ complex in a high-risk subset that includes myeloid, lymphoid and dendritic lineages. This dependency is characterized by innate inflammatory signalling and activation of phosphoinositide 3-kinase regulatory subunit 5 (PIK3R5), which encodes a regulatory subunit of PI3Kγ5 and stabilizes the active enzymatic complex. We identify p21 (RAC1)-activated kinase 1 (PAK1) as a noncanonical substrate of PI3Kγ that mediates this cell-intrinsic dependency and find that dephosphorylation of PAK1 by PI3Kγ inhibition impairs mitochondrial oxidative phosphorylation. Treatment with the selective PI3Kγ inhibitor eganelisib is effective in leukaemias with activated PIK3R5. In addition, the combination of eganelisib and cytarabine prolongs survival over either agent alone, even in patient-derived leukaemia xenografts with low baseline PIK3R5 expression, as residual leukaemia cells after cytarabine treatment have elevated G protein-coupled purinergic receptor activity and PAK1 phosphorylation. Together, our study reveals a targetable dependency on PI3Kγ-PAK1 signalling that is amenable to near-term evaluation in patients with acute leukaemia.

Multiomic profiling of human clonal hematopoiesis reveals genotype and cell-specific inflammatory pathway activation.

ABSTRACT: Clonal hematopoiesis (CH) is an age-associated phenomenon that increases the risk of hematologic malignancy and cardiovascular disease. CH is thought to enhance disease risk through inflammation in the peripheral blood.1 Here, we profile peripheral blood gene expression in 66 968 single cells from a cohort of 17 patients with CH and 7 controls. Using a novel mitochondrial DNA barcoding approach, we were able to identify and separately compare mutant Tet methylcytosine dioxygenase 2 (TET2) and DNA methyltransferase 3A (DNMT3A) cells with nonmutant counterparts. We discovered the vast majority of mutated cells were in the myeloid compartment. Additionally, patients harboring DNMT3A and TET2 CH mutations possessed a proinflammatory profile in CD14+ monocytes through previously unrecognized pathways such as galectin and macrophage inhibitory factor. We also found that T cells from patients with CH, although mostly unmutated, had decreased expression of GTPase of the immunity associated protein genes, which are critical to T-cell development, suggesting that CH impairs T-cell function.

Spectral Flow Cytometry Methods and Pipelines for Comprehensive Immunoprofiling of Human Peripheral Blood and Bone Marrow.

Profiling hematopoietic and immune cells provides important information about disease risk, disease status, and therapeutic responses. Spectral flow cytometry enables high-dimensional single-cell evaluation of large cohorts in a high-throughput manner. Here, we designed, optimized, and implemented new methods for deep immunophenotyping of human peripheral blood and bone marrow by spectral flow cytometry. Two blood antibody panels capture 48 cell-surface markers to assess more than 58 cell phenotypes, including subsets of T cells, B cells, monocytes, natural killer (NK) cells, and dendritic cells, and their respective markers of exhaustion, activation, and differentiation in less than 2 mL of blood. A bone marrow antibody panel captures 32 markers for 35 cell phenotypes, including stem/progenitor populations, T-cell subsets, dendritic cells, NK cells, and myeloid cells in a single tube. We adapted and developed innovative flow cytometric analysis algorithms, originally developed for single-cell genomics, to improve data integration and visualization. We also highlight technical considerations for users to ensure data fidelity. Our protocol and analysis pipeline accurately identifies rare cell types, discerns differences in cell abundance and phenotype across donors, and shows concordant immune landscape trends in patients with known hematologic malignancy. SIGNIFICANCE: This study introduces optimized methods and analysis algorithms that enhance capabilities in comprehensive immunophenotyping of human blood and bone marrow using spectral flow cytometry. This approach facilitates detection of rare cell types, enables measurement of cell variations across donors, and provides proof-of-concept in identifying known hematologic malignancies. By unlocking complexities of hematopoietic and immune landscapes at the single-cell level, this advancement holds potential for understanding disease states and therapeutic responses.

Deciphering cell states and genealogies of human haematopoiesis.

The human blood system is maintained through the differentiation and massive amplification of a limited number of long-lived haematopoietic stem cells (HSCs)1. Perturbations to this process underlie diverse diseases, but the clonal contributions to human haematopoiesis and how this changes with age remain incompletely understood. Although recent insights have emerged from barcoding studies in model systems2-5, simultaneous detection of cell states and phylogenies from natural barcodes in humans remains challenging. Here we introduce an improved, single-cell lineage-tracing system based on deep detection of naturally occurring mitochondrial DNA mutations with simultaneous readout of transcriptional states and chromatin accessibility. We use this system to define the clonal architecture of HSCs and map the physiological state and output of clones. We uncover functional heterogeneity in HSC clones, which is stable over months and manifests as both differences in total HSC output and biases towards the production of different mature cell types. We also find that the diversity of HSC clones decreases markedly with age, leading to an oligoclonal structure with multiple distinct clonal expansions. Our study thus provides a clonally resolved and cell-state-aware atlas of human haematopoiesis at single-cell resolution, showing an unappreciated functional diversity of human HSC clones and, more broadly, paving the way for refined studies of clonal dynamics across a range of tissues in human health and disease.

Ultraviolet radiation shapes dendritic cell leukaemia transformation in the skin.

Tumours most often arise from progression of precursor clones within a single anatomical niche. In the bone marrow, clonal progenitors can undergo malignant transformation to acute leukaemia, or differentiate into immune cells that contribute to disease pathology in peripheral tissues1-4. Outside the marrow, these clones are potentially exposed to a variety of tissue-specific mutational processes, although the consequences of this are unclear. Here we investigate the development of blastic plasmacytoid dendritic cell neoplasm (BPDCN)-an unusual form of acute leukaemia that often presents with malignant cells isolated to the skin5. Using tumour phylogenomics and single-cell transcriptomics with genotyping, we find that BPDCN arises from clonal (premalignant) haematopoietic precursors in the bone marrow. We observe that BPDCN skin tumours first develop at sun-exposed anatomical sites and are distinguished by clonally expanded mutations induced by ultraviolet (UV) radiation. A reconstruction of tumour phylogenies reveals that UV damage can precede the acquisition of alterations associated with malignant transformation, implicating sun exposure of plasmacytoid dendritic cells or committed precursors during BPDCN pathogenesis. Functionally, we find that loss-of-function mutations in Tet2, the most common premalignant alteration in BPDCN, confer resistance to UV-induced cell death in plasmacytoid, but not conventional, dendritic cells, suggesting a context-dependent tumour-suppressive role for TET2. These findings demonstrate how tissue-specific environmental exposures at distant anatomical sites can shape the evolution of premalignant clones to disseminated cancer.

"Candidatus Hydrogenisulfobacillus filiaventi" strain R50 gen. nov. sp. nov., a highly efficient producer of extracellular organic compounds from H2 and CO2.

Production of organic molecules is largely depending on fossil fuels. A sustainable alternative would be the synthesis of these compounds from CO2 and a cheap energy source, such as H2, CH4, NH3, CO, sulfur compounds or iron(II). Volcanic and geothermal areas are rich in CO2 and reduced inorganic gasses and therefore habitats where novel chemolithoautotrophic microorganisms for the synthesis of organic compounds could be discovered. Here we describe "Candidatus Hydrogenisulfobacillus filiaventi" R50 gen. nov., sp. nov., a thermoacidophilic, autotrophic H2-oxidizing microorganism, that fixed CO2 and excreted no less than 0.54 mol organic carbon per mole fixed CO2. Extensive metabolomics and NMR analyses revealed that Val, Ala and Ile are the most dominant form of excreted organic carbon while the aromatic amino acids Tyr and Phe, and Glu and Lys were present at much lower concentrations. In addition to these proteinogenic amino acids, the excreted carbon consisted of homoserine lactone, homoserine and an unidentified amino acid. The biological role of the excretion remains uncertain. In the laboratory, we noticed the production under high growth rates (0.034 h-1, doubling time of 20 h) in combination with O2-limitation, which will most likely not occur in the natural habitat of this strain. Nevertheless, this large production of extracellular organic molecules from CO2 may open possibilities to use chemolithoautotrophic microorganisms for the sustainable production of important biomolecules.

Chromatin complex dependencies reveal targeting opportunities in leukemia.

Chromatin regulators are frequently mutated in human cancer and are attractive drug targets. They include diverse proteins that share functional domains and assemble into related multi-subunit complexes. To investigate functional relationships among these regulators, here we apply combinatorial CRISPR knockouts (KOs) to test over 35,000 gene-gene pairings in leukemia cells, using a library of over 300,000 constructs. Top pairs that demonstrate either compensatory non-lethal interactions or synergistic lethality enrich for paralogs and targets that occupy the same protein complex. The screen highlights protein complex dependencies not apparent in single KO screens, for example MCM histone exchange, the nucleosome remodeling and deacetylase (NuRD) complex, and HBO1 (KAT7) complex. We explore two approaches to NuRD complex inactivation. Paralog and non-paralog combinations of the KAT7 complex emerge as synergistic lethal and specifically nominate the ING5 PHD domain as a potential therapeutic target when paired with other KAT7 complex member losses. These findings highlight the power of combinatorial screening to provide mechanistic insight and identify therapeutic targets within redundant networks.

CALR-mutated cells are vulnerable to combined inhibition of the proteasome and the endoplasmic reticulum stress response.

Cancer is driven by somatic mutations that provide a fitness advantage. While targeted therapies often focus on the mutated gene or its direct downstream effectors, imbalances brought on by cell-state alterations may also confer unique vulnerabilities. In myeloproliferative neoplasms (MPN), somatic mutations in the calreticulin (CALR) gene are disease-initiating through aberrant binding of mutant CALR to the thrombopoietin receptor MPL and ligand-independent activation of JAK-STAT signaling. Despite these mechanistic insights into the pathogenesis of CALR-mutant MPN, there are currently no mutant CALR-selective therapies available. Here, we identified differential upregulation of unfolded proteins, the proteasome and the ER stress response in CALR-mutant hematopoietic stem cells (HSCs) and megakaryocyte progenitors. We further found that combined pharmacological inhibition of the proteasome and IRE1-XBP1 axis of the ER stress response preferentially targets Calr-mutated HSCs and megakaryocytic-lineage cells over wild-type cells in vivo, resulting in an amelioration of the MPN phenotype. In serial transplantation assays following combined proteasome/IRE1 inhibition for six weeks, we did not find preferential depletion of Calr-mutant long-term HSCs. Together, these findings leverage altered proteostasis in Calr-mutant MPN to identify combinatorial dependencies that may be targeted for therapeutic benefit and suggest that eradicating disease-propagating Calr-mutant LT-HSCs may require more sustained treatment.

Targetable leukemia dependency on noncanonical PI3Kγ signaling.

Phosphoinositide 3-kinase gamma (PI3Kγ) is implicated as a target to repolarize tumor-associated macrophages and promote anti-tumor immune responses in solid cancers. However, cancer cell-intrinsic roles of PI3Kγ are unclear. Here, by integrating unbiased genome-wide CRISPR interference screening with functional analyses across acute leukemias, we define a selective dependency on the PI3Kγ complex in a high-risk subset that includes myeloid, lymphoid, and dendritic lineages. This dependency is characterized by innate inflammatory signaling and activation of phosphoinositide 3-kinase regulatory subunit 5 ( PIK3R5 ), which encodes a regulatory subunit of PI3Kγ and stabilizes the active enzymatic complex. Mechanistically, we identify p21 (RAC1) activated kinase 1 (PAK1) as a noncanonical substrate of PI3Kγ that mediates this cell-intrinsic dependency independently of Akt kinase. PI3Kγ inhibition dephosphorylates PAK1, activates a transcriptional network of NFκB-related tumor suppressor genes, and impairs mitochondrial oxidative phosphorylation. We find that treatment with the selective PI3Kγ inhibitor eganelisib is effective in leukemias with activated PIK3R5 , either at baseline or by exogenous inflammatory stimulation. Notably, the combination of eganelisib and cytarabine prolongs survival over either agent alone, even in patient-derived leukemia xenografts with low baseline PIK3R5 expression, as residual leukemia cells after cytarabine treatment have elevated G protein-coupled purinergic receptor activity and PAK1 phosphorylation. Taken together, our study reveals a targetable dependency on PI3Kγ/PAK1 signaling that is amenable to near-term evaluation in patients with acute leukemia.

Suppression of allograft rejection by regulatory B cells induced via TLR signaling.

B lymphocytes have long been recognized for their critical contributions to adaptive immunity, providing defense against pathogens through cognate antigen presentation to T cells and Ab production. More recently appreciated is that B cells are also integral in securing self-tolerance; this has led to interest in their therapeutic application to downregulate unwanted immune responses, such as transplant rejection. In this study, we found that PMA- and ionomycin-activated mouse B cells acquire regulatory properties following stimulation through TLR4/TLR9 receptors (Bregs-TLR). Bregs-TLR efficiently inhibited T cell proliferation in vitro and prevented allograft rejection. Unlike most reported Breg activities, the inhibition of alloimmune responses by Bregs-TLR relied on the expression of TGF-ß and not IL-10. In vivo, Bregs-TLR interrupted donor-specific T cell expansion and induced Tregs in a TGF-ß-dependent manner. RNA-Seq analyses corroborated the involvement of TGF-ß pathways in Breg-TLR function, identified potential gene pathways implicated in preventing graft rejection, and suggested targets to foster Breg regulation.

Clonal barcoding with qPCR detection enables live cell functional analyses for cancer research.

Single-cell analysis methods are valuable tools; however, current approaches do not easily enable live cell retrieval. That is a particular issue when further study of cells that were eliminated during experimentation could provide critical information. We report a clonal molecular barcoding method, called SunCatcher, that enables longitudinal tracking and live cell functional analysis. From complex cell populations, we generate single cell-derived clonal populations, infect each with a unique molecular barcode, and retain stocks of individual barcoded clones (BCs). We develop quantitative PCR-based and next-generation sequencing methods that we employ to identify and quantify BCs in vitro and in vivo. We apply SunCatcher to various breast cancer cell lines and combine respective BCs to create versions of the original cell lines. While the heterogeneous BC pools reproduce their original parental cell line proliferation and tumor progression rates, individual BCs are phenotypically and functionally diverse. Early spontaneous metastases can also be identified and quantified. SunCatcher thus provides a rapid and sensitive approach for studying live single-cell clones and clonal evolution, and performing functional analyses.

WAT3R: recovery of T-cell receptor variable regions from 3' single-cell RNA-sequencing.

SUMMARY: Diversity of the T-cell receptor (TCR) repertoire is central to adaptive immunity. The TCR is composed of α and ß chains, encoded by the TRA and TRB genes, of which the variable regions determine antigen specificity. To generate novel biological insights into the complex functioning of immune cells, combined capture of variable regions and single-cell transcriptomes provides a compelling approach. Recent developments enable the enrichment of TRA and TRB variable regions from widely used technologies for 3'-based single-cell RNA-sequencing (scRNA-seq). However, a comprehensive computational pipeline to process TCR-enriched data from 3' scRNA-seq is not available. Here, we present an analysis pipeline to process TCR variable regions enriched from 3' scRNA-seq cDNA. The tool reports TRA and TRB nucleotide and amino acid sequences linked to cell barcodes, enabling the reconstruction of T-cell clonotypes with associated transcriptomes. We demonstrate the software using peripheral blood mononuclear cells from a healthy donor and detect TCR sequences in a high proportion of single T cells. Detection of TCR sequences is low in non-T-cell populations, demonstrating specificity. Finally, we show that TCR clones are larger in CD8 Memory T cells than in other T-cell types, indicating an association between T-cell clonotypes and differentiation states. AVAILABILITY AND IMPLEMENTATION: The Workflow for Association of T-cell receptors from 3' single-cell RNA-seq (WAT3R), including test data, is available on GitHub (https://github.com/mainciburu/WAT3R), Docker Hub (https://hub.docker.com/r/mainciburu/wat3r) and a workflow on the Terra platform (https://app.terra.bio). The test dataset is available on GEO (accession number GSE195956). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Rise of the Clones: Myelodysplastic Syndrome to Secondary Acute Myeloid Leukemia.

Myelodysplastic syndrome (MDS) describes a family of blood disorders driven by the clonal expansion of mutated blood cells that can evolve into secondary acute myeloid leukemia (sAML). Two new studies use single-cell and deep sequencing to elucidate the progression of MDS to AML, revealing discrete clonal architectures and the driving role of signaling mutations.

A cellular hierarchy framework for understanding heterogeneity and predicting drug response in acute myeloid leukemia.

The treatment landscape of acute myeloid leukemia (AML) is evolving, with promising therapies entering clinical translation, yet patient responses remain heterogeneous, and biomarkers for tailoring treatment are lacking. To understand how disease heterogeneity links with therapy response, we determined the leukemia cell hierarchy makeup from bulk transcriptomes of more than 1,000 patients through deconvolution using single-cell reference profiles of leukemia stem, progenitor and mature cell types. Leukemia hierarchy composition was associated with functional, genomic and clinical properties and converged into four overall classes, spanning Primitive, Mature, GMP and Intermediate. Critically, variation in hierarchy composition along the Primitive versus GMP or Primitive versus Mature axes were associated with response to chemotherapy or drug sensitivity profiles of targeted therapies, respectively. A seven-gene biomarker derived from the Primitive versus Mature axis was associated with response to 105 investigational drugs. Cellular hierarchy composition constitutes a novel framework for understanding disease biology and advancing precision medicine in AML.

Single-Cell Multiomics Reveals Clonal T-Cell Expansions and Exhaustion in Blastic Plasmacytoid Dendritic Cell Neoplasm.

The immune system represents a major barrier to cancer progression, driving the evolution of immunoregulatory interactions between malignant cells and T-cells in the tumor environment. Blastic plasmacytoid dendritic cell neoplasm (BPDCN), a rare acute leukemia with plasmacytoid dendritic cell (pDC) differentiation, provides a unique opportunity to study these interactions. pDCs are key producers of interferon alpha (IFNA) that play an important role in T-cell activation at the interface between the innate and adaptive immune system. To assess how uncontrolled proliferation of malignant BPDCN cells affects the tumor environment, we catalog immune cell heterogeneity in the bone marrow (BM) of five healthy controls and five BPDCN patients by analyzing 52,803 single-cell transcriptomes, including 18,779 T-cells. We test computational techniques for robust cell type classification and find that T-cells in BPDCN patients consistently upregulate interferon alpha (IFNA) response and downregulate tumor necrosis factor alpha (TNFA) pathways. Integrating transcriptional data with T-cell receptor sequencing via shared barcodes reveals significant T-cell exhaustion in BPDCN that is positively correlated with T-cell clonotype expansion. By highlighting new mechanisms of T-cell exhaustion and immune evasion in BPDCN, our results demonstrate the value of single-cell multiomics to understand immune cell interactions in the tumor environment.