Mapping the cellular biogeography of human bone marrow niches using single-cell transcriptomics and proteomic imaging.

Non-hematopoietic cells are essential contributors to hematopoiesis. However, heterogeneity and spatial organization of these cells in human bone marrow remain largely uncharacterized. We used single-cell RNA sequencing (scRNA-seq) to profile 29,325 non-hematopoietic cells and discovered nine transcriptionally distinct subtypes. We simultaneously profiled 53,417 hematopoietic cells and predicted their interactions with non-hematopoietic subsets. We employed co-detection by indexing (CODEX) to spatially profile over 1.2 million cells. We integrated scRNA-seq and CODEX data to link predicted cellular signaling with spatial proximity. Our analysis revealed a hyperoxygenated arterio-endosteal neighborhood for early myelopoiesis, and an adipocytic localization for early hematopoietic stem and progenitor cells (HSPCs). We used our CODEX atlas to annotate new images and uncovered mesenchymal stromal cell (MSC) expansion and spatial neighborhoods co-enriched for leukemic blasts and MSCs in acute myeloid leukemia (AML) patient samples. This spatially resolved, multiomic atlas of human bone marrow provides a reference for investigation of cellular interactions that drive hematopoiesis.

Decade-long leukaemia remissions with persistence of CD4+ CAR T cells.

The adoptive transfer of T lymphocytes reprogrammed to target tumour cells has demonstrated potential for treatment of various cancers1-7. However, little is known about the long-term potential and clonal stability of the infused cells. Here we studied long-lasting CD19-redirected chimeric antigen receptor (CAR) T cells in two patients with chronic lymphocytic leukaemia1-4 who achieved a complete remission in 2010. CAR T cells remained detectable more than ten years after infusion, with sustained remission in both patients. Notably, a highly activated CD4+ population emerged in both patients, dominating the CAR T cell population at the later time points. This transition was reflected in the stabilization of the clonal make-up of CAR T cells with a repertoire dominated by a small number of clones. Single-cell profiling demonstrated that these long-persisting CD4+ CAR T cells exhibited cytotoxic characteristics along with ongoing functional activation and proliferation. In addition, longitudinal profiling revealed a population of gamma delta CAR T cells that prominently expanded in one patient concomitant with CD8+ CAR T cells during the initial response phase. Our identification and characterization of these unexpected CAR T cell populations provide novel insight into the CAR T cell characteristics associated with anti-cancer response and long-term remission in leukaemia.

Bone marrow adipogenic lineage precursors are the major regulator of bone resorption in adult mice.

Bone resorption by osteoclasts is a critical step in bone remodeling, a process important for maintaining bone homeostasis and repairing injured bone. We previously identified a bone marrow mesenchymal subpopulation, marrow adipogenic lineage precursors (MALPs), and showed that its production of RANKL stimulates bone resorption in young mice using Adipoq-Cre. To exclude developmental defects and to investigate the role of MALPs-derived RANKL in adult bone, we generated inducible reporter mice (Adipoq-CreER Tomato) and RANKL deficient mice (Adipoq-CreER RANKLflox/flox, iCKO). Single cell-RNA sequencing data analysis, lineage tracing, and in situ hybridization revealed that Adipoq+ cells contain not only MALPs but also late mesenchymal progenitors capable of osteogenic differentiation. However, RANKLmRNA was only detected in MALPs, but not in osteogenic cells. RANKL deficiency in MALPs induced at 3 months of age rapidly increased trabecular bone mass in long bones as well as vertebrae within 1 month due to diminished bone resorption but had no effect on the cortical bone. Ovariectomy (OVX) induced trabecular bone loss at both sites. RANKL depletion either before OVX or at 6 weeks post OVX protected and restored trabecular bone mass. Furthermore, bone healing after drill-hole injury was delayed in iCKO mice. Together, our findings demonstrate that MALPs play a dominant role in controlling trabecular bone resorption and that RANKL from MALPs is essential for trabecular bone turnover in adult bone homeostasis, postmenopausal bone loss, and injury repair.

Mapping the Cellular Biogeography of Human Bone Marrow Niches Using Single-Cell Transcriptomics and Proteomic Imaging.

The bone marrow is the organ responsible for blood production. Diverse non-hematopoietic cells contribute essentially to hematopoiesis. However, these cells and their spatial organization remain largely uncharacterized as they have been technically challenging to study in humans. Here, we used fresh femoral head samples and performed single-cell RNA sequencing (scRNA-Seq) to profile 29,325 enriched non-hematopoietic bone marrow cells and discover nine transcriptionally distinct subtypes. We next employed CO-detection by inDEXing (CODEX) multiplexed imaging of 18 individuals, including both healthy and acute myeloid leukemia (AML) samples, to spatially profile over one million single cells with a novel 53-antibody panel. We discovered a relatively hyperoxygenated arterio-endosteal niche for early myelopoiesis, and an adipocytic, but not endosteal or perivascular, niche for early hematopoietic stem and progenitor cells. We used our atlas to predict cell type labels in new bone marrow images and used these predictions to uncover mesenchymal stromal cell (MSC) expansion and leukemic blast/MSC-enriched spatial neighborhoods in AML patient samples. Our work represents the first comprehensive, spatially-resolved multiomic atlas of human bone marrow and will serve as a reference for future investigation of cellular interactions that drive hematopoiesis.

Impact of Genetic Ancestry on T-cell Acute Lymphoblastic Leukemia Outcomes.

The influence of genetic ancestry on biology, survival outcomes, and risk stratification in T-cell Acute Lymphoblastic Leukemia (T-ALL) has not been explored. Genetic ancestry was genomically-derived from DNA-based single nucleotide polymorphisms in children and young adults with T-ALL treated on Children's Oncology Group trial AALL0434. We determined associations of genetic ancestry, leukemia genomics and survival outcomes; co-primary outcomes were genomic subtype, pathway alteration, overall survival (OS), and event-free survival (EFS). Among 1309 patients, T-ALL molecular subtypes varied significantly by genetic ancestry, including increased frequency of genomically defined ETP-like, MLLT10, and BCL11B-activated subtypes in patients of African ancestry. In multivariable Cox models adjusting for high-risk subtype and pathways, patients of Admixed American ancestry had superior 5-year EFS/OS compared with European; EFS/OS for patients of African and European ancestry were similar. The prognostic value of five commonly altered T-ALL genes varied by ancestry - including NOTCH1 , which was associated with superior OS for patients of European and Admixed American ancestry but non-prognostic among patients of African ancestry. Furthermore, a published five-gene risk classifier accurately risk stratified patients of European ancestry, but misclassified patients of African ancestry. We developed a penalized Cox model which successfully risk stratified patients across ancestries. Overall, 80% of patients had a genomic alteration in at least one gene with differential prognostic impact by genetic ancestry. T-ALL genomics and prognostic associations of genomic alterations vary by genetic ancestry. These data demonstrate the importance of incorporating genetic ancestry into analyses of tumor biology for risk classification algorithms.

A longitudinal single-cell and spatial multiomic atlas of pediatric high-grade glioma.

Pediatric high-grade glioma (pHGG) is an incurable central nervous system malignancy that is a leading cause of pediatric cancer death. While pHGG shares many similarities to adult glioma, it is increasingly recognized as a molecularly distinct, yet highly heterogeneous disease. In this study, we longitudinally profiled a molecularly diverse cohort of 16 pHGG patients before and after standard therapy through single-nucleus RNA and ATAC sequencing, whole-genome sequencing, and CODEX spatial proteomics to capture the evolution of the tumor microenvironment during progression following treatment. We found that the canonical neoplastic cell phenotypes of adult glioblastoma are insufficient to capture the range of tumor cell states in a pediatric cohort and observed differential tumor-myeloid interactions between malignant cell states. We identified key transcriptional regulators of pHGG cell states and did not observe the marked proneural to mesenchymal shift characteristic of adult glioblastoma. We showed that essential neuromodulators and the interferon response are upregulated post-therapy along with an increase in non-neoplastic oligodendrocytes. Through in vitro pharmacological perturbation, we demonstrated novel malignant cell-intrinsic targets. This multiomic atlas of longitudinal pHGG captures the key features of therapy response that support distinction from its adult counterpart and suggests therapeutic strategies which are targeted to pediatric gliomas.

Csf1 from marrow adipogenic precursors is required for osteoclast formation and hematopoiesis in bone.

Colony-stimulating factor 1 (Csf1) is an essential growth factor for osteoclast progenitors and an important regulator for bone resorption. It remains elusive which mesenchymal cells synthesize Csf1 to stimulate osteoclastogenesis. We recently identified a novel mesenchymal cell population, marrow adipogenic lineage precursors (MALPs), in bone. Compared to other mesenchymal subpopulations, MALPs expressed Csf1 at a much higher level and this expression was further increased during aging. To investigate its role, we constructed MALP-deficient Csf1 CKO mice using AdipoqCre. These mice had increased femoral trabecular bone mass, but their cortical bone appeared normal. In comparison, depletion of Csf1 in the entire mesenchymal lineage using Prrx1Cre led to a more striking high bone mass phenotype, suggesting that additional mesenchymal subpopulations secrete Csf1. TRAP staining revealed diminished osteoclasts in the femoral secondary spongiosa region of Csf1 CKOAdipoq mice, but not at the chondral-osseous junction nor at the endosteal surface of cortical bone. Moreover, Csf1 CKOAdipoq mice were resistant to LPS-induced calvarial osteolysis. Bone marrow cellularity, hematopoietic progenitors, and macrophages were also reduced in these mice. Taken together, our studies demonstrate that MALPs synthesize Csf1 to control bone remodeling and hematopoiesis.

Identification and targeting of treatment resistant progenitor populations in T-cell Acute Lymphoblastic Leukemia.

Refractoriness to initial chemotherapy and relapse after remission are the main obstacles to cure in T-cell Acute Lymphoblastic Leukemia (T-ALL). Biomarker guided risk stratification and targeted therapy have the potential to improve outcomes in high-risk T-ALL; however, cellular and genetic factors contributing to treatment resistance remain unknown. Previous bulk genomic studies in T-ALL have implicated tumor heterogeneity as an unexplored mechanism for treatment failure. To link tumor subpopulations with clinical outcome, we created an atlas of healthy pediatric hematopoiesis and applied single-cell multiomic (CITE-seq/snATAC-seq) analysis to a cohort of 40 cases of T-ALL treated on the Children's Oncology Group AALL0434 clinical trial. The cohort was carefully selected to capture the immunophenotypic diversity of T-ALL, with early T-cell precursor (ETP) and Near/Non-ETP subtypes represented, as well as enriched with both relapsed and treatment refractory cases. Integrated analyses of T-ALL blasts and normal T-cell precursors identified a bone-marrow progenitor-like (BMP-like) leukemia sub-population associated with treatment failure and poor overall survival. The single-cell-derived molecular signature of BMP-like blasts predicted poor outcome across multiple subtypes of T-ALL within two independent patient cohorts using bulk RNA-sequencing data from over 1300 patients. We defined the mutational landscape of BMP-like T-ALL, finding that NOTCH1 mutations additively drive T-ALL blasts away from the BMP-like state. We transcriptionally matched BMP-like blasts to early thymic seeding progenitors that have low NR3C1 expression and high stem cell gene expression, corresponding to a corticosteroid and conventional cytotoxic resistant phenotype we observed in ex vivo drug screening. To identify novel targets for BMP-like blasts, we performed in silico and in vitro drug screening against the BMP-like signature and prioritized BMP-like overexpressed cell-surface (CD44, ITGA4, LGALS1) and intracellular proteins (BCL-2, MCL-1, BTK, NF-κB) as candidates for precision targeted therapy. We established patient derived xenograft models of BMP-high and BMP-low leukemias, which revealed vulnerability of BMP-like blasts to apoptosis-inducing agents, TEC-kinase inhibitors, and proteasome inhibitors. Our study establishes the first multi-omic signatures for rapid risk-stratification and targeted treatment of high-risk T-ALL.

Integrative Bulk and Single-Cell Profiling of Premanufacture T-cell Populations Reveals Factors Mediating Long-Term Persistence of CAR T-cell Therapy.

The adoptive transfer of chimeric antigen receptor (CAR) T cells represents a breakthrough in clinical oncology, yet both between- and within-patient differences in autologously derived T cells are a major contributor to therapy failure. To interrogate the molecular determinants of clinical CAR T-cell persistence, we extensively characterized the premanufacture T cells of 71 patients with B-cell malignancies on trial to receive anti-CD19 CAR T-cell therapy. We performed RNA-sequencing analysis on sorted T-cell subsets from all 71 patients, followed by paired Cellular Indexing of Transcriptomes and Epitopes (CITE) sequencing and single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq) on T cells from six of these patients. We found that chronic IFN signaling regulated by IRF7 was associated with poor CAR T-cell persistence across T-cell subsets, and that the TCF7 regulon not only associates with the favorable naïve T-cell state, but is maintained in effector T cells among patients with long-term CAR T-cell persistence. These findings provide key insights into the underlying molecular determinants of clinical CAR T-cell function. SIGNIFICANCE: To improve clinical outcomes for CAR T-cell therapy, there is a need to understand the molecular determinants of CAR T-cell persistence. These data represent the largest clinically annotated molecular atlas in CAR T-cell therapy to date, and significantly advance our understanding of the mechanisms underlying therapeutic efficacy.This article is highlighted in the In This Issue feature, p. 2113.

Identification of functionally primitive and immunophenotypically distinct subpopulations in secondary acute myeloid leukemia by mass cytometry.

BACKGROUND: Background: Mass cytometry (CyTOF) is a powerful tool for analyzing cellular networks at the single cell level. Due to the high-dimensional nature of this approach, analysis algorithms have been developed to visualize and interpret mass cytometry data. In this study, we applied these approaches to a cohort of patients with secondary acute myeloid leukemia (sAML). METHODS: We utilized mass cytometry to interrogate localization and intensity of thrombopoietin-mediated intracellular signaling in sAML. Extracellular and intracellular phenotypes were dissected using SPADE, viSNE, and PhenoGraph. RESULTS: Healthy controls exhibited highly localized signaling responses largely restricted to the hematopoietic stem/progenitor cell (HSPC) compartment. In contrast, sAML samples contained subpopulations outside the HSPC compartment exhibiting thrombopoietin (TPO) sensitivity comparable to or greater than immunophenotypically defined HSPCs. We employed unsupervised clustering by PhenoGraph to elucidate distinct subpopulations within these heterogeneous samples. One metacluster composed almost exclusively of Lin- CD61+ CD34- CD38- CD45low cells was identified. This subpopulation was not readily identified by established manual gating approaches, and generally exhibited greater STAT phosphorylation in response to TPO stimulation than did Lin- CD61- CD34+ CD38- cells. Lin- CD61+ CD34- CD38- CD45low cells were identified in three additional sAML patients analyzed independently using a manual gating approach based upon PhenoGraph results. Each patient exhibited a similar TPO hypersensitivity to the PhenoGraph metacluster. CONCLUSIONS: The identification of this cellular subpopulation highlights the limitations of manual gating in sAML. Our study demonstrates the potential for mass cytometry to elucidate rare subpopulations in highly heterogeneous tumors by utilizing unsupervised high dimensional analysis. © 2018 International Clinical Cytometry Society.

Analysis of Signaling Networks at the Single-Cell Level Using Mass Cytometry.

Mass cytometry is a powerful technology that enables the measurement of >40 parameters at the single-cell level. The inherent spectral limitations of fluorescent flow cytometry are circumvented by the use of antibodies conjugated to metal isotope reporters, which are measured quantitatively using a CyTOF mass cytometer. The high dimensionality of mass cytometry is particularly useful for the analysis of cell signaling networks in complex biological samples. We describe here methods for cell preparation, antibody staining, data acquisition, and analysis of multidimensional data from a mass cytometry experiment.

Cholesterol esterification inhibition and imatinib treatment synergistically inhibit growth of BCR-ABL mutation-independent resistant chronic myelogenous leukemia.

Since the advent of tyrosine kinase inhibitors (TKIs) such as imatinib, nilotinib, and dasatinib, chronic myelogenous leukemia (CML) prognosis has improved greatly. However, ~30-40% of patients develop resistance to imatinib therapy. Although most resistance is caused by mutations in the BCR-ABL kinase domain, 50-85% of these patients develop resistance in the absence of new mutations. In these cases, targeting other pathways may be needed to regain clinical response. Using label-free Raman spectromicroscopy, we evaluated a number of leukemia cell lines and discovered an aberrant accumulation of cholesteryl ester (CE) in CML, which was found to be a result of BCR-ABL kinase activity. CE accumulation in CML was found to be a cancer-specific phenomenon as untransformed cells did not accumulate CE. Blocking cholesterol esterification with avasimibe, a potent inhibitor of acyl-CoA cholesterol acyltransferase 1 (ACAT-1), significantly suppressed CML cell proliferation in Ba/F3 cells with the BCR-ABLT315I mutation and in K562 cells rendered imatinib resistant without mutations in the BCR-ABL kinase domain (K562R cells). Furthermore, the combination of avasimibe and imatinib caused a profound synergistic inhibition of cell proliferation in K562R cells, but not in Ba/F3T315I. This synergistic effect was confirmed in a K562R xenograft mouse model. Analysis of primary cells from a BCR-ABL mutation-independent imatinib resistant patient by mass cytometry suggested that the synergy may be due to downregulation of the MAPK pathway by avasimibe, which sensitized the CML cells to imatinib treatment. Collectively, these data demonstrate a novel strategy for overcoming BCR-ABL mutation-independent TKI resistance in CML.