Diversity of intratumoral regulatory T cells in B-cell non-Hodgkin lymphoma

Tumor-infiltrating regulatory T cells (Tregs) contribute to an immunosuppressive tumor microenvironment. Despite extensive studies, the prognostic impact of tumor-infiltrating Tregs in B-cell non-Hodgkin lymphomas (B-NHLs) remains unclear. Emerging studies suggest substantial heterogeneity in the phenotypes and suppressive capacities of Tregs, emphasizing the importance of understanding Treg diversity and the need for additional markers to identify highly suppressive Tregs. Here, we applied single-cell RNA sequencing and T-cell receptor sequencing combined with high-dimensional cytometry to decipher the heterogeneity of intratumoral Tregs in diffuse large B-cell lymphoma and follicular lymphoma (FL), compared with that in nonmalignant tonsillar tissue. We identified 3 distinct transcriptional states of Tregs: resting, activated, and unconventional LAG3+FOXP3- Tregs. Activated Tregs were enriched in B-NHL tumors, coexpressed several checkpoint receptors, and had stronger immunosuppressive activity compared with resting Tregs. In FL, activated Tregs were found in closer proximity to CD4+ and CD8+ T cells than other cell types. Furthermore, we used a computational approach to develop unique gene signature matrices, which were used to enumerate each Treg subset in cohorts with bulk gene expression data. In 2 independent FL cohorts, activated Tregs was the major subset, and high abundance was associated with adverse outcome. This study demonstrates that Tregs infiltrating B-NHL tumors are transcriptionally and functionally diverse. Highly immunosuppressive activated Tregs were enriched in tumor tissue but absent in the peripheral blood. Our data suggest that a deeper understanding of Treg heterogeneity in B-NHL could open new paths for rational drug design, facilitating selective targeting to improve antitumor immunity.

High-resolution alignment of single-cell and spatial transcriptomes with CytoSPACE

Recent studies have emphasized the importance of single-cell spatial biology, yet available assays for spatial transcriptomics have limited gene recovery or low spatial resolution. Here we introduce CytoSPACE, an optimization method for mapping individual cells from a single-cell RNA sequencing atlas to spatial expression profiles. Across diverse platforms and tissue types, we show that CytoSPACE outperforms previous methods with respect to noise tolerance and accuracy, enabling tissue cartography at single-cell resolution.

Profiling Cellular Ecosystems at Single-Cell Resolution and at Scale with EcoTyper

Tissues are composed of diverse cell types and cellular states that organize into distinct ecosystems with specialized functions. EcoTyper is a collection of machine learning tools for the large-scale delineation of cellular ecosystems and their constituent cell states from bulk, single-cell, and spatially resolved gene expression data. In this chapter, we provide a primer on EcoTyper and demonstrate its use for the discovery and recovery of cell states and ecosystems from healthy and diseased tissue specimens.

Identification of a minority population of LMO2+ breast cancer cells that integrate into the vasculature and initiate metastasis

Metastasis is responsible for most breast cancer-related deaths; however, identifying the cellular determinants of metastasis has remained challenging. Here, we identified a minority population of immature THY1+/VEGFA+ tumor epithelial cells in human breast tumor biopsies that display angiogenic features and are marked by the expression of the oncogene, LMO2. Higher abundance of LMO2+ basal cells correlated with tumor endothelial content and predicted poor distant recurrence-free survival in patients. Using MMTV-PyMT/Lmo2CreERT2 mice, we demonstrated that Lmo2 lineage-traced cells integrate into the vasculature and have a higher propensity to metastasize. LMO2 knockdown in human breast tumors reduced lung metastasis by impairing intravasation, leading to a reduced frequency of circulating tumor cells. Mechanistically, we find that LMO2 binds to STAT3 and is required for STAT3 activation by tumor necrosis factor-α and interleukin-6. Collectively, our study identifies a population of metastasis-initiating cells with angiogenic features and establishes the LMO2-STAT3 signaling axis as a therapeutic target in breast cancer metastasis.

Inferring gene expression from cell-free DNA fragmentation profiles

Profiling of circulating tumor DNA (ctDNA) in the bloodstream shows promise for noninvasive cancer detection. Chromatin fragmentation features have previously been explored to infer gene expression profiles from cell-free DNA (cfDNA), but current fragmentomic methods require high concentrations of tumor-derived DNA and provide limited resolution. Here we describe promoter fragmentation entropy as an epigenomic cfDNA feature that predicts RNA expression levels at individual genes. We developed 'epigenetic expression inference from cell-free DNA-sequencing' (EPIC-seq), a method that uses targeted sequencing of promoters of genes of interest. Profiling 329 blood samples from 201 patients with cancer and 87 healthy adults, we demonstrate classification of subtypes of lung carcinoma and diffuse large B cell lymphoma. Applying EPIC-seq to serial blood samples from patients treated with PD-(L)1 immune-checkpoint inhibitors, we show that gene expression profiles inferred by EPIC-seq are correlated with clinical response. Our results indicate that EPIC-seq could enable noninvasive, high-throughput tissue-of-origin characterization with diagnostic, prognostic and therapeutic potential.

T cell characteristics associated with toxicity to immune checkpoint blockade in patients with melanoma

Severe immune-related adverse events (irAEs) occur in up to 60% of patients with melanoma treated with immune checkpoint inhibitors (ICIs). However, it is unknown whether a common baseline immunological state precedes irAE development. Here we applied mass cytometry by time of flight, single-cell RNA sequencing, single-cell V(D)J sequencing, bulk RNA sequencing and bulk T cell receptor (TCR) sequencing to study peripheral blood samples from patients with melanoma treated with anti-PD-1 monotherapy or anti-PD-1 and anti-CTLA-4 combination ICIs. By analyzing 93 pre- and early on-ICI blood samples and 3 patient cohorts (n = 27, 26 and 18), we found that 2 pretreatment factors in circulation-activated CD4 memory T cell abundance and TCR diversity-are associated with severe irAE development regardless of organ system involvement. We also explored on-treatment changes in TCR clonality among patients receiving combination therapy and linked our findings to the severity and timing of irAE onset. These results demonstrate circulating T cell characteristics associated with ICI-induced toxicity, with implications for improved diagnostics and clinical management.

Atlas of clinically distinct cell states and ecosystems across human solid tumors

Determining how cells vary with their local signaling environment and organize into distinct cellular communities is critical for understanding processes as diverse as development, aging, and cancer. Here we introduce EcoTyper, a machine learning framework for large-scale identification and validation of cell states and multicellular communities from bulk, single-cell, and spatially resolved gene expression data. When applied to 12 major cell lineages across 16 types of human carcinoma, EcoTyper identified 69 transcriptionally defined cell states. Most states were specific to neoplastic tissue, ubiquitous across tumor types, and significantly prognostic. By analyzing cell-state co-occurrence patterns, we discovered ten clinically distinct multicellular communities with unexpectedly strong conservation, including three with myeloid and stromal elements linked to adverse survival, one enriched in normal tissue, and two associated with early cancer development. This study elucidates fundamental units of cellular organization in human carcinoma and provides a framework for large-scale profiling of cellular ecosystems in any tissue.

The landscape of tumor cell states and ecosystems in diffuse large B cell lymphoma

Biological heterogeneity in diffuse large B cell lymphoma (DLBCL) is partly driven by cell-of-origin subtypes and associated genomic lesions, but also by diverse cell types and cell states in the tumor microenvironment (TME). However, dissecting these cell states and their clinical relevance at scale remains challenging. Here, we implemented EcoTyper, a machine-learning framework integrating transcriptome deconvolution and single-cell RNA sequencing, to characterize clinically relevant DLBCL cell states and ecosystems. Using this approach, we identified five cell states of malignant B cells that vary in prognostic associations and differentiation status. We also identified striking variation in cell states for 12 other lineages comprising the TME and forming cell state interactions in stereotyped ecosystems. While cell-of-origin subtypes have distinct TME composition, DLBCL ecosystems capture clinical heterogeneity within existing subtypes and extend beyond cell-of-origin and genotypic classes. These results resolve the DLBCL microenvironment at systems-level resolution and identify opportunities for therapeutic targeting (https://ecotyper.stanford.edu/lymphoma).

Abstract 3443: Atlas of clinically-distinct cell states and cellular ecosystems across human solid tumors

Tumors are complex ecosystems consisting of malignant, immune, and stromal elements whose dynamic interactions drive patient survival and response to therapy. A comprehensive understanding of the diversity of cellular states within the tumor microenvironment (TME), and their patterns of co-occurrence, could provide new diagnostic tools for improved disease management and novel targets for therapeutic intervention. To address this challenge, we developed EcoTyper, a novel machine learning framework for large-scale identification of TME cell states and their co-association patterns from bulk, single-cell, and spatially resolved tumor expression data. EcoTyper starts by “purifying” cell type-specific gene expression profiles of epithelial cells, immune, and stromal cell types from bulk tissue transcriptomes using CIBERSORTx (Newman et al., Nat Biotechnol 2019). It then identifies transcriptional states for each cell type, validates them in scRNA-seq data, and uncovers co-occurrence patterns between cell states in order to define tumor cellular ecosystems. Applied to 6,475 tumor and adjacent normal samples from solid tumor types profiled by The Cancer Genome Atlas (TCGA), EcoTyper identified robust transcriptional states across 12 major cell types, including epithelial, fibroblast, endothelial, and 9 immune subsets. These states included both known and novel cellular phenotypes, nearly all of which could be validated in a compendium of scRNA-seq tumor atlases spanning ~140,000 cells. Most cell states were specific to neoplastic tissue, ubiquitous across tumor types, and significantly associated with overall survival, both in TCGA and in 9,062 held-out tumor specimens (Gentles/Newman et al., Nat Medicine 2015). We found that specific cell states co-occur in distinct cellular communities with characteristic patterns of ligand-receptor interactions, genomic features, clinical outcomes, and spatial organization. One such ecosystem defined a normal-like state that was strongly enriched in non-malignant samples. Others delineated novel pro- and anti-tumor inflammatory environments involving specific fibroblast, endothelial, and immune cell transcriptional programs. In summary, large-scale deconvolution of cell type-specific transcriptomes across thousands of solid tumors revealed a comprehensive atlas of TME cell states and cellular ecosystems. Our results provide a high-resolution portrait of cellular heterogeneity in the TME across multiple solid tumor types, with implications for novel diagnostics and immunotherapeutic targets.

References:

1. Newman, A.M., et al., Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat Biotechnol, 2019. 37(7): p. 773-782.

2. Gentles, A.J., et al., The prognostic landscape of genes and infiltrating immune cells across human cancers. Nature medicine, 2015. 21(8): p. 938.

Citation Format: Bogdan A. Luca, Chloé B. Steen, Armon Azizi, Magdalena Matusiak, Joanna Przybyl, Nastaran Neishaboori, Almudena Espín Pérez, Maximilian Diehn, Ash A. Alizadeh, Matt van de Rijn, Andrew J. Gentles, Aaron M. Newman. Atlas of clinically-distinct cell states and cellular ecosystems across human solid tumors [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 3443.

Profiling Cell Type Abundance and Expression in Bulk Tissues with CIBERSORTx

CIBERSORTx is a suite of machine learning tools for the assessment of cellular abundance and cell type-specific gene expression patterns from bulk tissue transcriptome profiles. With this framework, single-cell or bulk-sorted RNA sequencing data can be used to learn molecular signatures of distinct cell types from a small collection of biospecimens. These signatures can then be repeatedly applied to characterize cellular heterogeneity from bulk tissue transcriptomes without physical cell isolation. In this chapter, we provide a detailed primer on CIBERSORTx and demonstrate its capabilities for high-throughput profiling of cell types and cellular states in normal and neoplastic tissues.

Computational approaches for characterizing the tumor immune microenvironment

Recent advances in high-throughput molecular profiling technologies and multiplexed imaging platforms have revolutionized our ability to characterize the tumor immune microenvironment. As a result, studies of tumor-associated immune cells increasingly involve complex data sets that require sophisticated methods of computational analysis. In this review, we present an overview of key assays and related bioinformatics tools for analyzing the tumor-associated immune system in bulk tissues and at the single-cell level. In parallel, we describe how data science strategies and novel technologies have advanced tumor immunology and opened the door for new opportunities to exploit host immunity to improve cancer clinical outcomes.

Determining cell type abundance and expression from bulk tissues with digital cytometry

Single-cell RNA-sequencing has emerged as a powerful technique for characterizing cellular heterogeneity, but it is currently impractical on large sample cohorts and cannot be applied to fixed specimens collected as part of routine clinical care. We previously developed an approach for digital cytometry, called CIBERSORT, that enables estimation of cell type abundances from bulk tissue transcriptomes. We now introduce CIBERSORTx, a machine learning method that extends this framework to infer cell-type-specific gene expression profiles without physical cell isolation. By minimizing platform-specific variation, CIBERSORTx also allows the use of single-cell RNA-sequencing data for large-scale tissue dissection. We evaluated the utility of CIBERSORTx in multiple tumor types, including melanoma, where single-cell reference profiles were used to dissect bulk clinical specimens, revealing cell-type-specific phenotypic states linked to distinct driver mutations and response to immune checkpoint blockade. We anticipate that digital cytometry will augment single-cell profiling efforts, enabling cost-effective, high-throughput tissue characterization without the need for antibodies, disaggregation or viable cells.

T Cells Expressing Checkpoint Receptor TIGIT Are Enriched in Follicular Lymphoma Tumors and Characterized by Reversible Suppression of T-cell Receptor Signaling

Purpose: T cells infiltrating follicular lymphoma (FL) tumors are considered dysfunctional, yet the optimal target for immune checkpoint blockade is unknown. Characterizing coinhibitory receptor expression patterns and signaling responses in FL T-cell subsets might reveal new therapeutic targets.Experimental Design: Surface expression of 9 coinhibitory receptors governing T-cell function was characterized in T-cell subsets from FL lymph node tumors and from healthy donor tonsils and peripheral blood samples, using high-dimensional flow cytometry. The results were integrated with T-cell receptor (TCR)-induced signaling and cytokine production. Expression of T-cell immunoglobulin and ITIM domain (TIGIT) ligands was detected by immunohistochemistry.Results: TIGIT was a frequently expressed coinhibitory receptor in FL, expressed by the majority of CD8 T effector memory cells, which commonly coexpressed exhaustion markers such as PD-1 and CD244. CD8 FL T cells demonstrated highly reduced TCR-induced phosphorylation (p) of ERK and reduced production of IFNγ, while TCR proximal signaling (p-CD3ζ, p-SLP76) was not affected. The TIGIT ligands CD112 and CD155 were expressed by follicular dendritic cells in the tumor microenvironment. Dysfunctional TCR signaling correlated with TIGIT expression in FL CD8 T cells and could be fully restored upon in vitro culture. The costimulatory receptor CD226 was downregulated in TIGIT⁺ compared with TIGIT^- CD8 FL T cells, further skewing the balance toward immunosuppression.Conclusions: TIGIT blockade is a relevant strategy for improved immunotherapy in FL. A deeper understanding of the interplay between coinhibitory receptors and key T-cell signaling events can further assist in engineering immunotherapeutic regimens to improve clinical outcomes of cancer patients. Clin Cancer Res; 24(4); 870-81. ©2017 AACR.

ALIX and ESCRT-III coordinately control cytokinetic abscission during germline stem cell division in vivo

Abscission is the final step of cytokinesis that involves the cleavage of the intercellular bridge connecting the two daughter cells. Recent studies have given novel insight into the spatiotemporal regulation and molecular mechanisms controlling abscission in cultured yeast and human cells. The mechanisms of abscission in living metazoan tissues are however not well understood. Here we show that ALIX and the ESCRT-III component Shrub are required for completion of abscission during Drosophila female germline stem cell (fGSC) division. Loss of ALIX or Shrub function in fGSCs leads to delayed abscission and the consequent formation of stem cysts in which chains of daughter cells remain interconnected to the fGSC via midbody rings and fusome. We demonstrate that ALIX and Shrub interact and that they co-localize at midbody rings and midbodies during cytokinetic abscission in fGSCs. Mechanistically, we show that the direct interaction between ALIX and Shrub is required to ensure cytokinesis completion with normal kinetics in fGSCs. We conclude that ALIX and ESCRT-III coordinately control abscission in Drosophila fGSCs and that their complex formation is required for accurate abscission timing in GSCs in vivo.

Cytokinesis, the final step of cell division, concludes with a process termed abscission, during which the two daughter cells physically separate. In spite of their importance, the molecular machineries controlling abscission are poorly characterized especially in the context of living metazoan tissues. Here we provide molecular insight into the mechanism of abscission using the fruit fly Drosophila melanogaster as a model organism. We show that the scaffold protein ALIX and the ESCRT-III component Shrub are required for completion of abscission in Drosophila female germline stem cells (fGSCs). ESCRT-III has been implicated in topologically similar membrane scission events as abscission, namely intraluminal vesicle formation at endosomes and virus budding. Here we demonstrate that ALIX and Shrub co-localize and interact to promote abscission with correct timing in Drosophila fGSCs. We thus show that ALIX and ESCRT-III coordinately control abscission in Drosophila fGSCs cells and report an evolutionarily conserved function of the ALIX/ESCRT-III pathway during cytokinesis in a multi-cellular organism.