Clonal evolution in acute myeloid leukemia (AML) originates long before diagnosis and is a dynamic process that may affect survival. However, it remains uninvestigated during routine diagnostic workups. We hypothesized that the mutational status of bone marrow dysplastic cells and leukemic blasts, analyzed at the onset of AML using integrated multidimensional flow cytometry (MFC) immunophenotyping and fluorescence-activated cell sorting (FACS) with next-generation sequencing (NGS), could reconstruct leukemogenesis. Dysplastic cells were detected by MFC in 285 of 348 (82%) newly diagnosed patients with AML. Presence of dysplasia according to MFC and World Health Organization criteria had no prognostic value in older adults. NGS of dysplastic cells and blasts isolated at diagnosis identified 3 evolutionary patterns: stable (n= 12 of 21), branching (n= 4 of 21), and clonal evolution (n= 5 of 21). In patients achieving complete response (CR), integrated MFC and FACS with NGS showed persistent measurable residual disease (MRD) in phenotypically normal cell types, as well as the acquisition of genetic traits associated with treatment resistance. Furthermore, whole-exome sequencing of dysplastic and leukemic cells at diagnosis and of MRD uncovered different clonal involvement in dysplastic myelo-erythropoiesis, leukemic transformation, and chemoresistance. Altogether, we showed that it is possible to reconstruct leukemogenesis in 80% of patients with newly diagnosed AML, using techniques other than single-cell multiomics.

Early hematopoiesis is a continuous process in which hematopoietic stem and progenitor cells (HSPCs) gradually differentiate toward specific lineages. Aging and myeloid malignant transformation are characterized by changes in the composition and regulation of HSPCs. In this study, we used single-cell RNA sequencing (scRNA-seq) to characterize an enriched population of human HSPCs obtained from young and elderly healthy individuals. Based on their transcriptional profile, we identified changes in the proportions of progenitor compartments during aging, and differences in their functionality, as evidenced by gene set enrichment analysis. Trajectory inference revealed that altered gene expression dynamics accompanied cell differentiation, which could explain aging-associated changes in hematopoiesis. Next, we focused on key regulators of transcription by constructing gene regulatory networks (GRNs) and detected regulons that were specifically active in elderly individuals. Using previous findings in healthy cells as a reference, we analyzed scRNA-seq data obtained from patients with myelodysplastic syndrome (MDS) and detected specific alterations of the expression dynamics of genes involved in erythroid differentiation in all patients with MDS such as TRIB2. In addition, the comparison between transcriptional programs and GRNs regulating normal HSPCs and MDS HSPCs allowed identification of regulons that were specifically active in MDS cases such as SMAD1, HOXA6, POU2F2, and RUNX1 suggesting a role of these transcription factors (TFs) in the pathogenesis of the disease. In summary, we demonstrate that the combination of single-cell technologies with computational analysis tools enable the study of a variety of cellular mechanisms involved in complex biological systems such as early hematopoiesis and can be used to dissect perturbed differentiation trajectories associated with perturbations such as aging and malignant transformation. Furthermore, the identification of abnormal regulatory mechanisms associated with myeloid malignancies could be exploited for personalized therapeutic approaches in individual patients.

Identification of new markers associated with long-term efficacy in patients treated with CAR T cells is a current medical need, particularly in diseases such as multiple myeloma. In this study, we address the impact of CAR density on the functionality of BCMA CAR T cells. Functional and transcriptional studies demonstrate that CAR T cells with high expression of the CAR construct show an increased tonic signaling with up-regulation of exhaustion markers and increased in vitro cytotoxicity but a decrease in in vivo BM infiltration. Characterization of gene regulatory networks using scRNA-seq identified regulons associated to activation and exhaustion up-regulated in CARHigh T cells, providing mechanistic insights behind differential functionality of these cells. Last, we demonstrate that patients treated with CAR T cell products enriched in CARHigh T cells show a significantly worse clinical response in several hematological malignancies. In summary, our work demonstrates that CAR density plays an important role in CAR T activity with notable impact on clinical response.

Myeloid-derived suppressor cells (MDSCs) play a major role in cancer progression. In this study, we investigated the mechanisms by which complement C5a increases the capacity of polymorphonuclear MDSCs (PMN-MDSCs) to promote tumor growth and metastatic spread. Stimulation of PMN-MDSCs with C5a favored the invasion of cancer cells via a process dependent on the formation of neutrophil extracellular traps (NETs). NETosis was dependent on the production of high mobility group box 1 (HMGB1) by cancer cells. Moreover, C5a induced the surface expression of the HMGB1 receptors TLR4 and RAGE in PMN-MDSCs. In a mouse lung metastasis model, inhibition of C5a, C5a receptor-1 (C5aR1) or NETosis reduced the number of circulating-tumor cells (CTCs) and the metastatic burden. In support of the translational relevance of these findings, C5a was able to stimulate migration and NETosis in PMN-MDSCs obtained from lung cancer patients. Furthermore, myeloperoxidase (MPO)-DNA complexes, as markers of NETosis, were elevated in lung cancer patients and significantly correlated with C5a levels. In conclusion, C5a induces the formation of NETs from PMN-MDSCs in the presence of cancer cells, which may facilitate cancer cell dissemination and metastasis.

Understanding the regulation of normal and malignant human hematopoiesis requires comprehensive cell atlas of the hematopoietic stem cell (HSC) regulatory microenvironment. Here, we develop a tailored bioinformatic pipeline to integrate public and proprietary single-cell RNA sequencing (scRNA-seq) datasets. As a result, we robustly identify for the first time 14 intermediate cell states and 11 stages of differentiation in the endothelial and mesenchymal BM compartments, respectively. Our data provide the most comprehensive description to date of the murine HSC-regulatory microenvironment and suggest a higher level of specialization of the cellular circuits than previously anticipated. Furthermore, this deep characterization allows inferring conserved features in human, suggesting that the layers of microenvironmental regulation of hematopoiesis may also be shared between species. Our resource and methodology is a stepping-stone toward a comprehensive cell atlas of the BM microenvironment.

Large-scale immune monitoring is becoming routinely used in clinical trials to identify determinants of treatment responsiveness, particularly to immunotherapies. Flow cytometry remains one of the most versatile and high throughput approaches for single cell analysis; however, manual interpretation of multidimensional data poses a challenge when attempting to capture full cellular diversity and provide reproducible results. We present FlowCT, a semi-automated workspace empowered to analyze large data sets. It includes pre-processing, normalization, multiple dimensionality reduction techniques, automated clustering, and predictive modeling tools. As a proof of concept, we used FlowCT to compare the T-cell compartment in bone marrow (BM) with peripheral blood (PB) from patients with smoldering multiple myeloma (SMM), identify minimally invasive immune biomarkers of progression from smoldering to active MM, define prognostic T-cell subsets in the BM of patients with active MM after treatment intensification, and assess the longitudinal effect of maintenance therapy in BM T cells. A total of 354 samples were analyzed and immune signatures predictive of malignant transformation were identified in 150 patients with SMM (hazard ratio [HR], 1.7; P < .001). We also determined progression-free survival (HR, 4.09; P < .0001) and overall survival (HR, 3.12; P 5 .047) in 100 patients with active MM. New data also emerged about stem cell memory T cells, the concordance between immune profiles in BM and PB, and the immunomodulatory effect of maintenance therapy. FlowCT is a new open-source computational approach that can be readily implemented by research laboratories to perform quality control, analyze high-dimensional data, unveil cellular diversity, and objectively identify biomarkers in large immune monitoring studies. These trials were registered at www. clinicaltrials.gov as #NCT01916252 and #NCT02406144.

Vaccination using optimized strategies may increase response rates to immune checkpoint inhibitors (ICI) in some tumors. To enhance vaccine potency and improve thus responses to ICI, we analyzed the gene expression profile of an immunosuppressive dendritic cell (DC) population induced during vaccination, with the goal of identifying druggable inhibitory mechanisms. RNAseq studies revealed targetable genes, but their inhibition did not result in improved vaccines. However, we proved that immunosuppressive DC had a monocytic origin. Thus, monocyte depletion by gemcitabine administration reduced the generation of these DC and increased vaccine-induced immunity, which rejected about 20% of LLC-OVA and B16-OVA tumors, which are non-responders to anti-PD-1. This improved efficacy was associated with higher tumor T-cell infiltration and overexpression of PD-1/PD-L1. Therefore, the combination of vaccine + gemcitabine with anti-PD-1 was superior to anti-PD-1 monotherapy in both models. B16-OVA tumors benefited from a synergistic effect, reaching 75% of tumor rejection, but higher levels of exhausted T-cells in LLC-OVA tumors co-expressing PD-1, LAG3 and TIM3 precluded similar levels of efficacy. Our results indicate that gemcitabine is a suitable combination therapy with vaccines aimed at enhancing PD-1 therapies by targeting vaccine-induced immunosuppressive DC.

Normal cell counterparts of solid and myeloid tumors accumulate mutations years before disease onset; whether this occurs in B lymphocytes before lymphoma remains uncertain. We sequenced multiple stages of the B lineage in elderly individuals and patients with lymphoplasmacytic lymphoma, a singular disease for studying lymphomagenesis because of the high prevalence of mutated MYD88. We observed similar accumulation of random mutations in B lineages from both cohorts and unexpectedly found MYD88(L265P) in normal precursor and mature B lymphocytes from patients with lymphoma. We uncovered genetic and transcriptional pathways driving malignant transformation and leveraged these to model lymphoplasmacytic lymphoma in mice, based on mutated MYD88 in B cell precursors and BCL2 overexpression. Thus, MYD88(L265P) is a preneoplastic event, which challenges the current understanding of lymphomagenesis and may have implications for early detection of B cell lymphomas.

Adoptive cell transfer therapy using CD8+ T lymphocytes showed promising results eradicating metastatic malignancies. However, several regulatory mechanisms limit its efficacy. We studied the role of the expression of the transcription factor FOXP3 on CD8+ T cell function and anti-tumor immunity. Here we show that suboptimal T cell receptor stimulation of CD8+ T cells upregulates FOXP3 in vitro. Similarly, CD8 T cells transferred into tumor-bearing mice upregulate FOXP3 in vivo. Cell-intrinsic loss of FOXP3 by CD8+ T cells resulted in improved functionality after TCR stimulation and better antitumor responses in vivo. Inhibition of the FOXP3/NFAT interaction likewise improved CD8+ T cell functionality. Transcriptomic analysis of cells after TCR stimulation revealed an enrichment of genes implicated in the response to IFN-gamma, IFN-alpha, inflammatory response, IL-6/JAK/STAT, G2M checkpoint and IL-2/STAT signaling in FOXP3-deficient CD8+ T cells with respect to FOXP3-wt CD8+ T cells. Our results suggest that transient expression of FOXP3 by CD8+ T cells in the tumor microenvironment restrains their anti-tumor activity, with clear implications for improving T cell responses during immunotherapy.