TMExplorer: A tumour microenvironment single-cell RNAseq database and search tool

Modeling T cell temporal response to cancer immunotherapy rationalizes development of combinatorial treatment protocols

Successful immunotherapy relies on triggering complex responses involving T cell dynamics in tumors and the periphery. Characterizing these responses remains challenging using static human single-cell atlases or mouse models. To address this, we developed a framework for in vivo tracking of tumor-specific CD8+ T cells over time and at single-cell resolution. Our tools facilitate the modeling of gene program dynamics in the tumor microenvironment (TME) and the tumor-draining lymph node (tdLN). Using this approach, we characterize two modes of anti-programmed cell death protein 1 (PD-1) activity, decoupling induced differentiation of tumor-specific activated precursor cells from conventional type 1 dendritic cell (cDC1)-dependent proliferation and recruitment to the TME. We demonstrate that combining anti-PD-1 therapy with anti-4-1BB agonist enhances the recruitment and proliferation of activated precursors, resulting in tumor control. These data suggest that effective response to anti-PD-1 therapy is dependent on sufficient influx of activated precursor CD8+ cells to the TME and highlight the importance of understanding system-level dynamics in optimizing immunotherapies.

Cytokine expression patterns: A single-cell RNA sequencing and machine learning based roadmap for cancer classification

Cytokines are small protein molecules that exhibit potent immunoregulatory properties, which are known as the essential components of the tumor immune microenvironment (TIME). While some cytokines are known to be universally upregulated in TIME, the unique cytokine expression patterns have not been fully resolved in specific types of cancers. To address this challenge, we develop a TIME single-cell RNA sequencing (scRNA-seq) dataset, which is designed to study cytokine expression patterns for precise cancer classification. The dataset, including 39 cancers, is constructed by integrating 684 tumor scRNA-seq samples from multiple public repositories. After screening and processing, the dataset retains only the expression data of immune cells. With a machine learning classification model, unique cytokine expression patterns are identified for various cancer categories and pioneering applied to cancer classification with an accuracy rate of 78.01%. Our method will not only boost the understanding of cancer-type-specific immune modulations in TIME but also serve as a crucial reference for future diagnostic and therapeutic research in cancer immunity.

SCAR: Single-cell and Spatially-resolved Cancer Resources

Advances in sequencing and imaging technologies offer a unique opportunity to unravel cell heterogeneity and develop new immunotherapy strategies for cancer research. There is an urgent need for a resource that effectively integrates a vast amount of transcriptomic profiling data to comprehensively explore cancer tissue heterogeneity and the tumor microenvironment. In this context, we developed the Single-cell and Spatially-resolved Cancer Resources (SCAR) database, a combined tumor spatial and single-cell transcriptomic platform, which is freely accessible at http://8.142.154.29/SCAR2023 or http://scaratlas.com. SCAR contains spatial transcriptomic data from 21 tumor tissues and single-cell transcriptomic data from 11 301 352 cells encompassing 395 cancer subtypes and covering a wide variety of tissues, organoids, and cell lines. This resource offers diverse functional modules to address key cancer research questions at multiple levels, including the screening of tumor cell types, metabolic features, cell communication and gene expression patterns within the tumor microenvironment. Moreover, SCAR enables the analysis of biomarker expression patterns and cell developmental trajectories. SCAR also provides a comprehensive analysis of multi-dimensional datasets based on 34 state-of-the-art omics techniques, serving as an essential tool for in-depth mining and understanding of cell heterogeneity and spatial location. The implications of this resource extend to both cancer biology research and cancer immunotherapy development.

signifinder enables the identification of tumor cell states and cancer expression signatures in bulk, single-cell and spatial transcriptomic data

Over the last decade, many studies and some clinical trials have proposed gene expression signatures as a valuable tool for understanding cancer mechanisms, defining subtypes, monitoring patient prognosis, and therapy efficacy. However, technical and biological concerns about reproducibility have been raised. Technical reproducibility is a major concern: we currently lack a computational implementation of the proposed signatures, which would provide detailed signature definition and assure reproducibility, dissemination, and usability of the classifier. Another concern regards intratumor heterogeneity, which has never been addressed when studying these types of biomarkers using bulk transcriptomics. With the aim of providing a tool able to improve the reproducibility and usability of gene expression signatures, we propose signifinder, an R package that provides the infrastructure to collect, implement, and compare expression-based signatures from cancer literature. The included signatures cover a wide range of biological processes from metabolism and programmed cell death, to morphological changes, such as quantification of epithelial or mesenchymal-like status. Collected signatures can score tumor cell characteristics, such as the predicted response to therapy or the survival association, and can quantify microenvironmental information, including hypoxia and immune response activity. signifinder has been used to characterize tumor samples and to investigate intra-tumor heterogeneity, extending its application to single-cell and spatial transcriptomic data. Through these higher-resolution technologies, it has become increasingly apparent that the single-sample score assessment obtained by transcriptional signatures is conditioned by the phenotypic and genetic intratumor heterogeneity of tumor masses. Since the characteristics of the most abundant cell type or clone might not necessarily predict the properties of mixed populations, signature prediction efficacy is lowered, thus impeding effective clinical diagnostics. Through signifinder, we offer general principles for interpreting and comparing transcriptional signatures, as well as suggestions for additional signatures that would allow for more complete and robust data inferences. We consider signifinder a useful tool to pave the way for reproducibility and comparison of transcriptional signatures in oncology.

Evaluation of single-cell RNAseq labelling algorithms using cancer datasets

Single-cell RNA sequencing (scRNA-seq) clustering and labelling methods are used to determine precise cellular composition of tissue samples. Automated labelling methods rely on either unsupervised, cluster-based approaches or supervised, cell-based approaches to identify cell types. The high complexity of cancer poses a unique challenge, as tumor microenvironments are often composed of diverse cell subpopulations with unique functional effects that may lead to disease progression, metastasis and treatment resistance. Here, we assess 17 cell-based and 9 cluster-based scRNA-seq labelling algorithms using 8 cancer datasets, providing a comprehensive large-scale assessment of such methods in a cancer-specific context. Using several performance metrics, we show that cell-based methods generally achieved higher performance and were faster compared to cluster-based methods. Cluster-based methods more successfully labelled non-malignant cell types, likely because of a lack of gene signatures for relevant malignant cell subpopulations. Larger cell numbers present in some cell types in training data positively impacted prediction scores for cell-based methods. Finally, we examined which methods performed favorably when trained and tested on separate patient cohorts in scenarios similar to clinical applications, and which were able to accurately label particularly small or under-represented cell populations in the given datasets. We conclude that scPred and SVM show the best overall performances with cancer-specific data and provide further suggestions for algorithm selection. Our analysis pipeline for assessing the performance of cell type labelling algorithms is available in https://github.com/shooshtarilab/scRNAseq-Automated-Cell-Type-Labelling.