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.

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.