Building and analyzing metacells in single-cell genomics data.

The advent of high-throughput single-cell genomics technologies has fundamentally transformed biological sciences. Currently, millions of cells from complex biological tissues can be phenotypically profiled across multiple modalities. The scaling of computational methods to analyze and visualize such data is a constant challenge, and tools need to be regularly updated, if not redesigned, to cope with ever-growing numbers of cells. Over the last few years, metacells have been introduced to reduce the size and complexity of single-cell genomics data while preserving biologically relevant information and improving interpretability. Here, we review recent studies that capitalize on the concept of metacells-and the many variants in nomenclature that have been used. We further outline how and when metacells should (or should not) be used to analyze single-cell genomics data and what should be considered when analyzing such data at the metacell level. To facilitate the exploration of metacells, we provide a comprehensive tutorial on the construction and analysis of metacells from single-cell RNA-seq data ( https://github.com/GfellerLab/MetacellAnalysisTutorial ) as well as a fully integrated pipeline to rapidly build, visualize and evaluate metacells with different methods ( https://github.com/GfellerLab/MetacellAnalysisToolkit ).

Metacells untangle large and complex single-cell transcriptome networks.

BACKGROUND: Single-cell RNA sequencing (scRNA-seq) technologies offer unique opportunities for exploring heterogeneous cell populations. However, in-depth single-cell transcriptomic characterization of complex tissues often requires profiling tens to hundreds of thousands of cells. Such large numbers of cells represent an important hurdle for downstream analyses, interpretation and visualization. RESULTS: We develop a framework called SuperCell to merge highly similar cells into metacells and perform standard scRNA-seq data analyses at the metacell level. Our systematic benchmarking demonstrates that metacells not only preserve but often improve the results of downstream analyses including visualization, clustering, differential expression, cell type annotation, gene correlation, imputation, RNA velocity and data integration. By capitalizing on the redundancy inherent to scRNA-seq data, metacells significantly facilitate and accelerate the construction and interpretation of single-cell atlases, as demonstrated by the integration of 1.46 million cells from COVID-19 patients in less than two hours on a standard desktop. CONCLUSIONS: SuperCell is a framework to build and analyze metacells in a way that efficiently preserves the results of scRNA-seq data analyses while significantly accelerating and facilitating them.

Tumor-specific cytolytic CD4 T cells mediate immunity against human cancer.

CD4 T cells have been implicated in cancer immunity for their helper functions. Moreover, their direct cytotoxic potential has been shown in some patients with cancer. Here, by mining single-cell RNA-seq datasets, we identified CD4 T cell clusters displaying cytotoxic phenotypes in different human cancers, resembling CD8 T cell profiles. Using the peptide-MHCII-multimer technology, we confirmed ex vivo the presence of cytolytic tumor-specific CD4 T cells. We performed an integrated phenotypic and functional characterization of these cells, down to the single-cell level, through a high-throughput nanobiochip consisting of massive arrays of picowells and machine learning. We demonstrated a direct, contact-, and granzyme-dependent cytotoxic activity against tumors, with delayed kinetics compared to classical cytotoxic lymphocytes. Last, we found that this cytotoxic activity was in part dependent on SLAMF7. Agonistic engagement of SLAMF7 enhanced cytotoxicity of tumor-specific CD4 T cells, suggesting that targeting these cells might prove synergistic with other cancer immunotherapies.

Deciphering the transcriptomic landscape of tumor-infiltrating CD8 lymphocytes in B16 melanoma tumors with single-cell RNA-Seq.

Recent studies have proposed that tumor-specific tumor-infiltrating CD8+ T lymphocytes (CD8 TIL) can be classified into two main groups: "exhausted" TILs, characterized by high expression of the inhibitory receptors PD-1 and TIM-3 and lack of transcription factor 1 (Tcf1); and "memory-like" TILs, with self-renewal capacity and co-expressing Tcf1 and PD-1. However, a comprehensive definition of the heterogeneity existing within CD8 TILs has yet to be clearly established. To investigate this heterogeneity at the transcriptomic level, we performed paired single-cell RNA and TCR sequencing of CD8 T cells infiltrating B16 murine melanoma tumors, including cells of known tumor specificity. Unsupervised clustering and gene-signature analysis revealed four distinct CD8 TIL states - exhausted, memory-like, naïve and effector memory-like (EM-like) - and predicted novel markers, including Ly6C for the EM-like cells, that were validated by flow cytometry. Tumor-specific PMEL T cells were predominantly found within the exhausted and memory-like states but also within the EM-like state. Further, T cell receptor sequencing revealed a large clonal expansion of exhausted, memory-like and EM-like cells with partial clonal relatedness between them. Finally, meta-analyses of public bulk and single-cell RNA-seq data suggested that anti-PD-1 treatment induces the expansion of EM-like cells. Our reference map of the transcriptomic landscape of murine CD8 TILs will help interpreting future bulk and single-cell transcriptomic studies and may guide the analysis of CD8IL subpopulations in response to therapeutic interventions.