Coordinated Immune Cell Networks in the Bone Marrow Microenvironment Define the Graft versus Leukemia Response with Adoptive Cellular Therapy.

Understanding how intra-tumoral immune populations coordinate to generate anti-tumor responses following therapy can guide precise treatment prioritization. We performed systematic dissection of an established adoptive cellular therapy, donor lymphocyte infusion (DLI), by analyzing 348,905 single-cell transcriptomes from 74 longitudinal bone-marrow samples of 25 patients with relapsed myeloid leukemia; a subset was evaluated by protein-based spatial analysis. In acute myelogenous leukemia (AML) responders, diverse immune cell types within the bone-marrow microenvironment (BME) were predicted to interact with a clonally expanded population of ZNF683 + GZMB + CD8+ cytotoxic T lymphocytes (CTLs) which demonstrated in vitro specificity for autologous leukemia. This population, originating predominantly from the DLI product, expanded concurrently with NK and B cells. AML nonresponder BME revealed a paucity of crosstalk and elevated TIGIT expression in CD8+ CTLs. Our study highlights recipient BME differences as a key determinant of effective anti-leukemia response and opens new opportunities to modulate cell-based leukemia-directed therapy.

Integrative genotyping of cancer and immune phenotypes by long-read sequencing.

Single-cell transcriptomics has become the definitive method for classifying cell types and states, and can be augmented with genotype information to improve cell lineage identification. Due to constraints of short-read sequencing, current methods to detect natural genetic barcodes often require cumbersome primer panels and early commitment to targets. Here we devise a flexible long-read sequencing workflow and analysis pipeline, termed nanoranger, that starts from intermediate single-cell cDNA libraries to detect cell lineage-defining features, including single-nucleotide variants, fusion genes, isoforms, sequences of chimeric antigen and TCRs. Through systematic analysis of these classes of natural 'barcodes', we define the optimal targets for nanoranger, namely those loci close to the 5' end of highly expressed genes with transcript lengths shorter than 4 kB. As proof-of-concept, we apply nanoranger to longitudinal tracking of subclones of acute myeloid leukemia (AML) and describe the heterogeneous isoform landscape of thousands of marrow-infiltrating immune cells. We propose that enhanced cellular genotyping using nanoranger can improve the tracking of single-cell tumor and immune cell co-evolution.

A multidimensional analysis reveals distinct immune phenotypes and tertiary lymphoid structure-like aggregates in the bone marrow of pediatric acute myeloid leukemia.

Because of the low mutational burden, children with acute myeloid leukemia (AML) are thought to have a 'cold' tumor microenvironment and consequently, a low likelihood of response to T cell-directed immunotherapies. Here, we provide a multidimensional overview of the tumor immune microenvironment in newly diagnosed pediatric AML. On a cohort level, we demonstrate wide variation in T cell infiltration with nearly one-third of cases harboring an immune-infiltrated bone marrow. These immune-infiltrated cases are characterized by a decreased abundance of M2-like macrophages, which we find to be associated with response to T cell-directed immunotherapy in adult AML. On an organizational level, we reveal the composition of spatially organized immune aggregates in pediatric AML, and show that in the adult setting such aggregates in post-treatment bone marrow and extramedullary sites associate with response to ipilimumab-based therapy. Altogether, our study provides immune correlates of response to T cell-directed immunotherapies and indicates starting points for further investigations into immunomodulatory mechanisms in AML.

Preexisting tumor-resident T cells with cytotoxic potential associate with response to neoadjuvant anti-PD-1 in head and neck cancer.

About 50% of patients with locally advanced head and neck squamous cell carcinoma (HNSCC) experience recurrences after definitive therapy. The presurgical administration of anti-programmed cell death protein 1 (PD-1) immunotherapy results in substantial pathologic tumor responses (pTR) within the tumor microenvironment (TME). However, the mechanisms underlying the dynamics of antitumor T cells upon neoadjuvant PD-1 blockade remain unresolved, and approaches to increase pathologic responses are lacking. In a phase 2 trial (NCT02296684), we observed that 45% of patients treated with two doses of neoadjuvant pembrolizumab experienced marked pTRs (≥50%). Single-cell analysis of 17,158 CD8+ T cells from 14 tumor biopsies, including 6 matched pre-post neoadjuvant treatment, revealed that responding tumors had clonally expanded putative tumor-specific exhausted CD8+ tumor-infiltrating lymphocytes (TILs) with a tissue-resident memory program, characterized by high cytotoxic potential (CTX+) and ZNF683 expression, within the baseline TME. Pathologic responses after 5 weeks of PD-1 blockade were consistent with activation of preexisting CTX+ZNF683+CD8+ TILs, paralleling loss of viable tumor and associated tumor antigens. Response was associated with high numbers of CD103+PD-1+CD8+ T cells infiltrating pretreatment lesions, whereas revival of nonexhausted persisting clones and clonal replacement were modest. By contrast, nonresponder baseline TME exhibited a relative absence of ZNF683+CTX+ TILs and subsequent accumulation of highly exhausted clones. In HNSCC, revival of preexisting ZNF683+CTX+ TILs is a major mechanism of response in the immediate postneoadjuvant setting.

Mechanisms of response and resistance to combined decitabine and ipilimumab for advanced myeloid disease.

The challenge of eradicating leukemia in patients with acute myelogenous leukemia (AML) after initial cytoreduction has motivated modern efforts to combine synergistic active modalities including immunotherapy. Recently, the ETCTN/CTEP 10026 study tested the combination of the DNA methyltransferase inhibitor decitabine together with the immune checkpoint inhibitor ipilimumab for AML/myelodysplastic syndrome (MDS) either after allogeneic hematopoietic stem cell transplantation (HSCT) or in the HSCT-naïve setting. Integrative transcriptome-based analysis of 304 961 individual marrow-infiltrating cells for 18 of 48 subjects treated on study revealed the strong association of response with a high baseline ratio of T to AML cells. Clinical responses were predominantly driven by decitabine-induced cytoreduction. Evidence of immune activation was only apparent after ipilimumab exposure, which altered CD4+ T-cell gene expression, in line with ongoing T-cell differentiation and increased frequency of marrow-infiltrating regulatory T cells. For post-HSCT samples, relapse could be attributed to insufficient clearing of malignant clones in progenitor cell populations. In contrast to AML/MDS bone marrow, the transcriptomes of leukemia cutis samples from patients with durable remission after ipilimumab monotherapy showed evidence of increased infiltration with antigen-experienced resident memory T cells and higher expression of CTLA-4 and FOXP3. Altogether, activity of combined decitabine and ipilimumab is impacted by cellular expression states within the microenvironmental niche of leukemic cells. The inadequate elimination of leukemic progenitors mandates urgent development of novel approaches for targeting these cell populations to generate long-lasting responses. This trial was registered at www.clinicaltrials.gov as #NCT02890329.

In Vivo Modeling of CLL Transformation to Richter Syndrome Reveals Convergent Evolutionary Paths and Therapeutic Vulnerabilities.

Transformation to aggressive disease histologies generates formidable clinical challenges across cancers, but biological insights remain few. We modeled the genetic heterogeneity of chronic lymphocytic leukemia (CLL) through multiplexed in vivo CRISPR-Cas9 B-cell editing of recurrent CLL loss-of-function drivers in mice and recapitulated the process of transformation from indolent CLL into large cell lymphoma [i.e., Richter syndrome (RS)]. Evolutionary trajectories of 64 mice carrying diverse combinatorial gene assortments revealed coselection of mutations in Trp53, Mga, and Chd2 and the dual impact of clonal Mga/Chd2 mutations on E2F/MYC and interferon signaling dysregulation. Comparative human and murine RS analyses demonstrated tonic PI3K signaling as a key feature of transformed disease, with constitutive activation of the AKT and S6 kinases, downmodulation of the PTEN phosphatase, and convergent activation of MYC/PI3K transcriptional programs underlying enhanced sensitivity to MYC/mTOR/PI3K inhibition. This robust experimental system presents a unique framework to study lymphoid biology and therapy. SIGNIFICANCE: Mouse models reflective of the genetic complexity and heterogeneity of human tumors remain few, including those able to recapitulate transformation to aggressive disease histologies. Herein, we model CLL transformation into RS through multiplexed in vivo gene editing, providing key insight into the pathophysiology and therapeutic vulnerabilities of transformed disease. This article is highlighted in the In This Issue feature, p. 101.

Mitochondrial DNA Mutations as Natural Barcodes for Lineage Tracing of Murine Tumor Models.

Murine models are indispensable tools for functional genomic studies and preclinical testing of novel therapeutic approaches. Mitochondrial single-cell assay for transposase-accessible chromatin using sequencing (mtscATAC-seq) enables the dissection of cellular heterogeneity and clonal dynamics by capturing chromatin accessibility, copy-number variations (CNV), and mitochondrial DNA (mtDNA) mutations, yet its applicability to murine studies remains unexplored. By leveraging mtscATAC-seq in novel chronic lymphocytic leukemia and Richter syndrome mouse models, we report the detection of mtDNA mutations, particularly in highly proliferative murine cells, alongside CNV and chromatin state changes indicative of clonal evolution upon secondary transplant. This study thus demonstrates the feasibility and utility of multi-modal single-cell and natural barcoding approaches to characterize murine cancer models. SIGNIFICANCE: mtDNA mutations can serve as natural barcodes to enable lineage tracing in murine cancer models, which can be used to provide new insights into disease biology and to identify therapeutic vulnerabilities.

Immune Phenotypes and Target Antigens of Clonally Expanded Bone Marrow T Cells in Treatment-Naïve Multiple Myeloma.

Multiple myeloma is a hematologic malignancy of monoclonal plasma cells that accumulate in the bone marrow. Despite their clinical and pathophysiologic relevance, the roles of bone marrow-infiltrating T cells in treatment-naïve patients are incompletely understood. We investigated whether clonally expanded T cells (i) were detectable in multiple myeloma bone marrow, (ii) showed characteristic immune phenotypes, and (iii) whether dominant clones recognized antigens selectively presented on multiple myeloma cells. Single-cell index sorting and T-cell receptor (TCR) αß sequencing of bone marrow T cells from 13 treatment-naïve patients showed dominant clonal expansion within CD8+ cytolytic effector compartments, and only a minority of expanded T-cell clones expressed the classic immune-checkpoint molecules PD-1, CTLA-4, or TIM-3. To identify their molecular targets, TCRs of 68 dominant bone marrow clones from five selected patients were reexpressed and incubated with multiple myeloma and non-multiple myeloma cells from corresponding patients. Only 1 of 68 TCRs recognized antigen presented on multiple myeloma cells. This TCR was HLA-C-restricted, self-peptide-specific and could be activated by multiple myeloma cells of multiple patients. The remaining dominant T-cell clones did not recognize multiple myeloma cells and were, in part, specific for antigens associated with chronic viral infections. In conclusion, we showed that dominant bone marrow T-cell clones in treatment-naïve patients rarely recognize antigens presented on multiple myeloma cells and exhibit low expression of classic immune-checkpoint molecules. Our data provide experimental context for experiences from clinical immune-checkpoint inhibition trials and will inform future T cell-dependent therapeutic strategies.

Rapid single-cell identification of Epstein-Barr virus-specific T-cell receptors for cellular therapy.

BACKGROUND AND AIMS: Epstein-Barr virus (EBV) is associated with solid and hematopoietic malignancies. After allogeneic stem cell transplantation, EBV infection or reactivation represents a potentially life-threatening condition with no specific treatment available in clinical routine. In vitro expansion of naturally occurring EBV-specific T cells for adoptive transfer is time-consuming and influenced by the donor's T-cell receptor (TCR) repertoire and requires a specific memory compartment that is non-existent in seronegative individuals. The authors present highly efficient identification of EBV-specific TCRs that can be expressed on human T cells and recognize EBV-infected cells. METHODS AND RESULTS: Mononuclear cells from six stem cell grafts were expanded in vitro with three HLA-B*35:01- or four HLA-A*02:01-presented peptides derived from six EBV proteins expressed during latent and lytic infection. Epitope-specific T cells expanded on average 42-fold and were single-cell-sorted and TCRαß-sequenced. To confirm specificity, 11 HLA-B*35:01- and six HLA-A*02:01-restricted dominant TCRs were expressed on reporter cell lines, and 16 of 17 TCRs recognized their presumed target peptides. To confirm recognition of virus-infected cells and assess their value for adoptive therapy, three selected HLA-B*35:01- and four HLA-A*02:01-restricted TCRs were expressed on human peripheral blood lymphocytes. All TCR-transduced cells recognized EBV-infected lymphoblastoid cell lines. CONCLUSIONS: The authors' approach provides sets of EBV epitope-specific TCRs in two different HLA contexts. Resulting cellular products do not require EBV-seropositive donors, can be adjusted to cell subsets of choice with exactly defined proportions of target-specific T cells, can be tracked in vivo and will help to overcome unmet clinical needs in the treatment and prophylaxis of EBV reactivation and associated malignancies.

Improved T-cell Immunity Following Neoadjuvant Chemotherapy in Ovarian Cancer.

PURPOSE: Although local tissue-based immune responses are critical for elucidating direct tumor-immune cell interactions, peripheral immune responses are increasingly recognized as occupying an important role in anticancer immunity. We evaluated serial blood samples from patients with advanced epithelial ovarian cancer (EOC) undergoing standard-of-care neoadjuvant carboplatin and paclitaxel chemotherapy (including dexamethasone for prophylaxis of paclitaxel-associated hypersensitivity reactions) to characterize the evolution of the peripheral immune cell function and composition across the course of therapy. EXPERIMENTAL DESIGN: Serial blood samples from 10 patients with advanced high-grade serous ovarian cancer treated with neoadjuvant chemotherapy (NACT) were collected before the initiation of chemotherapy, after the third and sixth cycles, and approximately 2 months after completion of chemotherapy. T-cell function was evaluated using ex vivo IFNγ ELISpot assays, and the dynamics of T-cell repertoire and immune cell composition were assessed using bulk and single-cell RNA sequencing (RNAseq). RESULTS: T cells exhibited an improved response to viral antigens after NACT, which paralleled the decrease in CA125 levels. Single-cell analysis revealed increased numbers of memory T-cell receptor (TCR) clonotypes and increased central memory CD8+ and regulatory T cells throughout chemotherapy. Finally, administration of NACT was associated with increased monocyte frequency and expression of HLA class II and antigen presentation genes; single-cell RNAseq analyses showed that although driven largely by classical monocytes, increased class II gene expression was a feature observed across monocyte subpopulations after chemotherapy. CONCLUSIONS: NACT may alleviate tumor-associated immunosuppression by reducing tumor burden and may enhance antigen processing and presentation. These findings have implications for the successful combinatorial applications of immune checkpoint blockade and therapeutic vaccine approaches in EOC.

Coevolving JAK2V617F+relapsed AML and donor T cells with PD-1 blockade after stem cell transplantation: an index case.

Relapse of myeloproliferative neoplasms (MPNs) after allogeneic hematopoietic stem cell transplantation (HSCT) is associated with poor outcomes, as therapeutic approaches to reinstate effective graft-versus-leukemia (GVL) responses remain suboptimal. Immune escape through overexpression of PD-L1 in JAK2V617F-mutated MPN provides a rationale for therapeutic PD-1 blockade, and indeed, clinical activity of nivolumab in relapsed MPN post-HSCT has been observed. Elucidation of the features of response following PD-1 blockade in such patients could inform novel therapeutic concepts that enhance GVL. Here, we report an integrated high-dimensional analysis using single-cell RNA sequencing, T-cell receptor sequencing, cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq), and assay for transposase-accessible chromatin using sequencing (scATAC-seq), together with mass cytometry, in peripheral blood mononuclear cells collected at 6 timepoints before, during, and after transient response to PD-1 blockade from an index case of relapsed MPN following HSCT. Before nivolumab infusion, acute myeloid leukemia (AML) blasts demonstrated high expression of chemokines, and T cells were characterized by expression of interferon-response genes. This baseline inflammatory signature disappeared after nivolumab infusion. Clinical response was characterized by transient expansion of a polyclonal CD4+ T-cell population and contraction of an AML subpopulation that exhibited megakaryocytic features and elevated PD-L1 expression. At relapse, the proportion of the AML subpopulation with progenitor-like features progressively increased, suggesting coevolution of AML blasts and donor-derived T cells. We thus demonstrate how single-cell technologies can provide complementary insight into cellular mechanisms underlying response to PD-1 blockade, motivating future longitudinal high-dimensional single-cell studies of GVL responses in relapsed myeloid disease.

Longitudinal Single-Cell Dynamics of Chromatin Accessibility and Mitochondrial Mutations in Chronic Lymphocytic Leukemia Mirror Disease History.

While cancers evolve during disease progression and in response to therapy, temporal dynamics remain difficult to study in humans due to the lack of consistent barcodes marking individual clones in vivo. We employ mitochondrial single-cell assay for transposase-accessible chromatin with sequencing to profile 163,279 cells from 9 patients with chronic lymphocytic leukemia (CLL) collected across disease course and utilize mitochondrial DNA (mtDNA) mutations as natural genetic markers of cancer clones. We observe stable propagation of mtDNA mutations over years in the absence of strong selective pressure, indicating clonal persistence, but dramatic changes following tight bottlenecks, including disease transformation and relapse posttherapy, paralleled by acquisition of copy-number variants and changes in chromatin accessibility and gene expression. Furthermore, we link CLL subclones to distinct chromatin states, providing insight into nongenetic sources of relapse. mtDNA mutations thus mirror disease history and provide naturally occurring genetic barcodes to enable patient-specific study of cancer subclonal dynamics. SIGNIFICANCE: Single-cell multi-omic profiling of CLL reveals the utility of somatic mtDNA mutations as in vivo barcodes, which mark subclones that can evolve over time along with changes in accessible chromatin and gene expression profiles to capture dynamics of disease evolution. See related commentary by Hilton and Scott, p. 2965. This article is highlighted in the In This Issue feature, p. 2945.

Molecular and cellular features of CTLA-4 blockade for relapsed myeloid malignancies after transplantation.

Relapsed myeloid disease after allogeneic stem cell transplantation (HSCT) remains largely incurable. We previously demonstrated the potent activity of immune checkpoint blockade in this clinical setting with ipilimumab or nivolumab. To define the molecular and cellular pathways by which CTLA-4 blockade with ipilimumab can reinvigorate an effective graft-versus-leukemia (GVL) response, we integrated transcriptomic analysis of leukemic biopsies with immunophenotypic profiling of matched peripheral blood samples collected from patients treated with ipilimumab following HSCT on the Experimental Therapeutics Clinical Trials Network 9204 trial. Response to ipilimumab was associated with transcriptomic evidence of increased local CD8+ T-cell infiltration and activation. Systemically, ipilimumab decreased naïve and increased memory T-cell populations and increased expression of markers of T-cell activation and costimulation such as PD-1, HLA-DR, and ICOS, irrespective of response. However, responding patients were characterized by higher turnover of T-cell receptor sequences in peripheral blood and showed increased expression of proinflammatory chemokines in plasma that was further amplified by ipilimumab. Altogether, these data highlight the compositional T-cell shifts and inflammatory pathways induced by ipilimumab both locally and systemically that associate with successful GVL outcomes. This trial was registered at www.clinicaltrials.gov as #NCT01822509.

Natural Barcodes for Longitudinal Single Cell Tracking of Leukemic and Immune Cell Dynamics.

Blood malignancies provide unique opportunities for longitudinal tracking of disease evolution following therapeutic bottlenecks and for the monitoring of changes in anti-tumor immunity. The expanding development of multi-modal single-cell sequencing technologies affords newer platforms to elucidate the mechanisms underlying these processes at unprecedented resolution. Furthermore, the identification of molecular events that can serve as in-vivo barcodes now facilitate the tracking of the trajectories of malignant and of immune cell populations over time within primary human samples, as these permit unambiguous identification of the clonal lineage of cell populations within heterogeneous phenotypes. Here, we provide an overview of the potential for chromosomal copy number changes, somatic nuclear and mitochondrial DNA mutations, single nucleotide polymorphisms, and T and B cell receptor sequences to serve as personal natural barcodes and review technical implementations in single-cell analysis workflows. Applications of these methodologies include the study of acquired therapeutic resistance and the dissection of donor- and host cellular interactions in the context of allogeneic hematopoietic stem cell transplantation.

Personal tumor antigens in blood malignancies: genomics-directed identification and targeting.

Hematological malignancies have long been at the forefront of the development of novel immune-based treatment strategies. The earliest successful efforts originated from the extensive body of work in the field of allogeneic hematopoietic stem cell transplantation. These efforts laid the foundation for the recent exciting era of cancer immunotherapy, which includes immune checkpoint blockade, personal neoantigen vaccines, and adoptive T cell transfer. At the heart of the specificity of these novel strategies is the recognition of target antigens presented by malignant cells to T cells. Here, we review the advances in systematic identification of minor histocompatibility antigens and neoantigens arising from personal somatic alterations or recurrent driver mutations. These exciting efforts pave the path for the implementation of personalized combinatorial cancer therapy.

Localization-associated immune phenotypes of clonally expanded tumor-infiltrating T cells and distribution of their target antigens in rectal cancer.

The degree and type of T cell infiltration influence rectal cancer prognosis regardless of classical tumor staging. We asked whether clonal expansion and tumor infiltration are restricted to selected-phenotype T cells; which clones are accessible in peripheral blood; and what the spatial distribution of their target antigens is. From five rectal cancer patients, we isolated paired tumor-infiltrating T cells (TILs) and T cells from unaffected rectum mucosa (TUM) using 13-parameter FACS single cell index sorting. TCRαß sequences, cytokine, and transcription factor expression were determined with single cell sequencing. TILs and TUM occupied distinct phenotype compartments and clonal expansion predominantly occurred within CD8+ T cells. Expanded TIL clones identified by paired TCRαß sequencing and exclusively detectable in the tumor showed characteristic PD-1 and TIM-3 expression. TCRß repertoire sequencing identified 49 out of 149 expanded TIL clones circulating in peripheral blood and 41 (84%) of these were PD-1- TIM-3-. To determine whether clonal expansion of predominantly tumor-infiltrating T cell clones was driven by antigens uniquely presented in tumor tissue, selected TCRs were reconstructed and incubated with cells isolated from corresponding tumor or unaffected mucosa. The majority of clones exclusively detected in the tumor recognized antigen at both sites. In summary, rectal cancer is infiltrated with expanded distinct-phenotype T cell clones that either i) predominantly infiltrate the tumor, ii) predominantly infiltrate the unaffected mucosa, or iii) overlap between tumor, unaffected mucosa, and peripheral blood. However, the target antigens of predominantly tumor-infiltrating TIL clones do not appear to be restricted to tumor tissue.

FACS single cell index sorting is highly reliable and determines immune phenotypes of clonally expanded T cells.

FACS index sorting allows the isolation of single cells with retrospective identification of each single cell's high-dimensional immune phenotype. We experimentally determine the error rate of index sorting and combine the technology with T cell receptor sequencing to identify clonal T cell expansion in aplastic anemia bone marrow as an example.

Clonal Expansion and Interrelatedness of Distinct B-Lineage Compartments in Multiple Myeloma Bone Marrow.

Multiple myeloma is characterized by the clonal expansion of malignant plasma cells in the bone marrow. But the phenotypic diversity and the contribution of less predominant B-lineage clones to the biology of this disease have been controversial. Here, we asked whether cells bearing the dominant multiple myeloma immunoglobulin rearrangement occupy phenotypic compartments other than that of plasma cells. To accomplish this, we combined 13-parameter FACS index sorting and t-Stochastic Neighbor Embedding (t-SNE) visualization with high-throughput single-cell immunoglobulin sequencing to track selected B-lineage clones across different stages of human B-cell development. As expected, the predominant clones preferentially mapped to aberrant plasma cell compartments, albeit phenotypically altered from wild type. Interestingly, up to 1.2% of cells of the predominant clones colocalized with B-lineage cells of a normal phenotype. In addition, minor clones with distinct immunoglobulin sequences were detected in up to 9% of sequenced cells, but only 2 out of 12 of these clones showed aberrant immune phenotypes. The majority of these minor clones showed intraclonal silent nucleotide differences within the CDR3s and varying frequencies of somatic mutations in the immunoglobulin genes. Therefore, the phenotypic range of multiple myeloma cells in the bone marrow is not confined to aberrant-phenotype plasma cells but extends to low frequencies of normal-phenotype B cells, in line with the recently reported success of B cell-targeting cellular therapies in some patients. The majority of minor clones result from parallel nonmalignant expansion. Cancer Immunol Res; 5(9); 744-54. ©2017 AACR.