IL-10-expressing CAR T cells resist dysfunction and mediate durable clearance of solid tumors and metastases.

The success of chimeric antigen receptor (CAR) T cell therapy in treating several hematopoietic malignancies has been difficult to replicate in solid tumors, in part because of T cell exhaustion and eventually dysfunction. To counter T cell dysfunction in the tumor microenvironment, we metabolically armored CAR T cells by engineering them to secrete interleukin-10 (IL-10). We show that IL-10 CAR T cells preserve intact mitochondrial structure and function in the tumor microenvironment and increase oxidative phosphorylation in a mitochondrial pyruvate carrier-dependent manner. IL-10 secretion promoted proliferation and effector function of CAR T cells, leading to complete regression of established solid tumors and metastatic cancers across several cancer types in syngeneic and xenograft mouse models, including colon cancer, breast cancer, melanoma and pancreatic cancer. IL-10 CAR T cells also induced stem cell-like memory responses in lymphoid organs that imparted durable protection against tumor rechallenge. Our results establish a generalizable approach to counter CAR T cell dysfunction through metabolic armoring, leading to solid tumor eradication and long-lasting immune protection.

Semi-supervised integration of single-cell transcriptomics data.

Batch effects in single-cell RNA-seq data pose a significant challenge for comparative analyses across samples, individuals, and conditions. Although batch effect correction methods are routinely applied, data integration often leads to overcorrection and can result in the loss of biological variability. In this work we present STACAS, a batch correction method for scRNA-seq that leverages prior knowledge on cell types to preserve biological variability upon integration. Through an open-source benchmark, we show that semi-supervised STACAS outperforms state-of-the-art unsupervised methods, as well as supervised methods such as scANVI and scGen. STACAS scales well to large datasets and is robust to incomplete and imprecise input cell type labels, which are commonly encountered in real-life integration tasks. We argue that the incorporation of prior cell type information should be a common practice in single-cell data integration, and we provide a flexible framework for semi-supervised batch effect correction.

Activation of the transcription factor NFAT5 in the tumor microenvironment enforces CD8+ T cell exhaustion.

Persistent exposure to antigen during chronic infection or cancer renders T cells dysfunctional. The molecular mechanisms regulating this state of exhaustion are thought to be common in infection and cancer, despite obvious differences in their microenvironments. Here we found that NFAT5, an NFAT family transcription factor that lacks an AP-1 docking site, was highly expressed in exhausted CD8+ T cells in the context of chronic infections and tumors but was selectively required in tumor-induced CD8+ T cell exhaustion. Overexpression of NFAT5 in CD8+ T cells reduced tumor control, while deletion of NFAT5 improved tumor control by promoting the accumulation of tumor-specific CD8+ T cells that had reduced expression of the exhaustion-associated proteins TOX and PD-1 and produced more cytokines, such as IFNÉ£ and TNF, than cells with wild-type levels of NFAT5, specifically in the precursor exhausted PD-1+TCF1+TIM-3-CD8+ T cell population. NFAT5 did not promote T cell exhaustion during chronic infection with clone 13 of lymphocytic choriomeningitis virus. Expression of NFAT5 was induced by TCR triggering, but its transcriptional activity was specific to the tumor microenvironment and required hyperosmolarity. Thus, NFAT5 promoted the exhaustion of CD8+ T cells in a tumor-selective fashion.

T Cell Clonal Analysis Using Single-cell RNA Sequencing and Reference Maps.

T cells are endowed with T-cell antigen receptors (TCR) that give them the capacity to recognize specific antigens and mount antigen-specific adaptive immune responses. Because TCR sequences are distinct in each naïve T cell, they serve as molecular barcodes to track T cells with clonal relatedness and shared antigen specificity through proliferation, differentiation, and migration. Single-cell RNA sequencing provides coupled information of TCR sequence and transcriptional state in individual cells, enabling T-cell clonotype-specific analyses. In this protocol, we outline a computational workflow to perform T-cell states and clonal analysis from scRNA-seq data based on the R packages Seurat, ProjecTILs, and scRepertoire. Given a scRNA-seq T-cell dataset with TCR sequence information, cell states are automatically annotated by reference projection using the ProjecTILs method. TCR information is used to track individual clonotypes, assess their clonal expansion, proliferation rates, bias towards specific differentiation states, and the clonal overlap between T-cell subtypes. We provide fully reproducible R code to conduct these analyses and generate useful visualizations that can be adapted for the needs of the protocol user. Key features Computational analysis of paired scRNA-seq and scTCR-seq data Characterizing T-cell functional state by reference-based analysis using ProjecTILs Exploring T-cell clonal structure using scRepertoire Linking T-cell clonality to transcriptomic state to study relationships between clonal expansion and functional phenotype Graphical overview.

A CSF-1R-blocking antibody/IL-10 fusion protein increases anti-tumor immunity by effectuating tumor-resident CD8+ T cells.

Strategies to increase intratumoral concentrations of an anticancer agent are desirable to optimize its therapeutic potential when said agent is efficacious primarily within a tumor but also have significant systemic side effects. Here, we generate a bifunctional protein by fusing interleukin-10 (IL-10) to a colony-stimulating factor-1 receptor (CSF-1R)-blocking antibody. The fusion protein demonstrates significant antitumor activity in multiple cancer models, especially head and neck cancer. Moreover, this bifunctional protein not only leads to the anticipated reduction in tumor-associated macrophages but also triggers proliferation, activation, and metabolic reprogramming of CD8+ T cells. Furthermore, it extends the clonotype diversity of tumor-infiltrated T cells and shifts the tumor microenvironment (TME) to an immune-active state. This study suggests an efficient strategy for designing immunotherapeutic agents by fusing a potent immunostimulatory molecule to an antibody targeting TME-enriched factors.

Orthogonal cytokine engineering enables novel synthetic effector states escaping canonical exhaustion in tumor-rejecting CD8+ T cells.

To date, no immunotherapy approaches have managed to fully overcome T-cell exhaustion, which remains a mandatory fate for chronically activated effector cells and a major therapeutic challenge. Understanding how to reprogram CD8+ tumor-infiltrating lymphocytes away from exhausted effector states remains an elusive goal. Our work provides evidence that orthogonal gene engineering of T cells to secrete an interleukin (IL)-2 variant binding the IL-2Rßγ receptor and the alarmin IL-33 reprogrammed adoptively transferred T cells to acquire a novel, synthetic effector state, which deviated from canonical exhaustion and displayed superior effector functions. These cells successfully overcame homeostatic barriers in the host and led-in the absence of lymphodepletion or exogenous cytokine support-to high levels of engraftment and tumor regression. Our work unlocks a new opportunity of rationally engineering synthetic CD8+ T-cell states endowed with the ability to avoid exhaustion and control advanced solid tumors.

Dissecting the treatment-naive ecosystem of human melanoma brain metastasis.

Melanoma brain metastasis (MBM) frequently occurs in patients with advanced melanoma; yet, our understanding of the underlying salient biology is rudimentary. Here, we performed single-cell/nucleus RNA-seq in 22 treatment-naive MBMs and 10 extracranial melanoma metastases (ECMs) and matched spatial single-cell transcriptomics and T cell receptor (TCR)-seq. Cancer cells from MBM were more chromosomally unstable, adopted a neuronal-like cell state, and enriched for spatially variably expressed metabolic pathways. Key observations were validated in independent patient cohorts, patient-derived MBM/ECM xenograft models, RNA/ATAC-seq, proteomics, and multiplexed imaging. Integrated spatial analyses revealed distinct geography of putative cancer immune evasion and evidence for more abundant intra-tumoral B to plasma cell differentiation in lymphoid aggregates in MBM. MBM harbored larger fractions of monocyte-derived macrophages and dysfunctional TOX+CD8+ T cells with distinct expression of immune checkpoints. This work provides comprehensive insights into MBM biology and serves as a foundational resource for further discovery and therapeutic exploration.

A CD4+ T cell reference map delineates subtype-specific adaptation during acute and chronic viral infections.

CD4+ T cells are critical orchestrators of immune responses against a large variety of pathogens, including viruses. While multiple CD4+ T cell subtypes and their key transcriptional regulators have been identified, there is a lack of consistent definition for CD4+ T cell transcriptional states. In addition, the progressive changes affecting CD4+ T cell subtypes during and after immune responses remain poorly defined. Using single-cell transcriptomics, we characterized the diversity of CD4+ T cells responding to self-resolving and chronic viral infections in mice. We built a comprehensive map of virus-specific CD4+ T cells and their evolution over time, and identified six major cell states consistently observed in acute and chronic infections. During the course of acute infections, T cell composition progressively changed from effector to memory states, with subtype-specific gene modules and kinetics. Conversely, in persistent infections T cells acquired distinct, chronicity-associated programs. By single-cell T cell receptor (TCR) analysis, we characterized the clonal structure of virus-specific CD4+ T cells across individuals. Virus-specific CD4+ T cell responses were essentially private across individuals and most T cells differentiated into both Tfh and Th1 subtypes irrespective of their TCR. Finally, we showed that our CD4+ T cell map can be used as a reference to accurately interpret cell states in external single-cell datasets across tissues and disease models. Overall, this study describes a previously unappreciated level of adaptation of the transcriptional states of CD4+ T cells responding to viruses and provides a new computational resource for CD4+ T cell analysis.

scGate: marker-based purification of cell types from heterogeneous single-cell RNA-seq datasets.

SUMMARY: A common bioinformatics task in single-cell data analysis is to purify a cell type or cell population of interest from heterogeneous datasets. Here, we present scGate, an algorithm that automatizes marker-based purification of specific cell populations, without requiring training data or reference gene expression profiles. scGate purifies a cell population of interest using a set of markers organized in a hierarchical structure, akin to gating strategies employed in flow cytometry. scGate outperforms state-of-the-art single-cell classifiers and it can be applied to multiple modalities of single-cell data (e.g. RNA-seq, ATAC-seq, CITE-seq). scGate is implemented as an R package and integrated with the Seurat framework, providing an intuitive tool to isolate cell populations of interest from heterogeneous single-cell datasets. AVAILABILITY AND IMPLEMENTATION: scGate is available as an R package at https://github.com/carmonalab/scGate (https://doi.org/10.5281/zenodo.6202614). Several reproducible workflows describing the main functions and usage of the package on different single-cell modalities, as well as the code to reproduce the benchmark, can be found at https://github.com/carmonalab/scGate.demo (https://doi.org/10.5281/zenodo.6202585) and https://github.com/carmonalab/scGate.benchmark. Test data are available at https://doi.org/10.6084/m9.figshare.16826071. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Low-Dose Radiotherapy Reverses Tumor Immune Desertification and Resistance to Immunotherapy.

Developing strategies to inflame tumors is critical for increasing response to immunotherapy. Here, we report that low-dose radiotherapy (LDRT) of murine tumors promotes T-cell infiltration and enables responsiveness to combinatorial immunotherapy in an IFN-dependent manner. Treatment efficacy relied upon mobilizing both adaptive and innate immunity and depended on both cytotoxic CD4+ and CD8+ T cells. LDRT elicited predominantly CD4+ cells with features of exhausted effector cytotoxic cells, with a subset expressing NKG2D and exhibiting proliferative capacity, as well as a unique subset of activated dendritic cells expressing the NKG2D ligand RAE1. We translated these findings to a phase I clinical trial administering LDRT, low-dose cyclophosphamide, and immune checkpoint blockade to patients with immune-desert tumors. In responsive patients, the combinatorial treatment triggered T-cell infiltration, predominantly of CD4+ cells with Th1 signatures. Our data support the rational combination of LDRT with immunotherapy for effectively treating low T cell-infiltrated tumors. SIGNIFICANCE: Low-dose radiation reprogrammed the tumor microenvironment of tumors with scarce immune infiltration and together with immunotherapy induced simultaneous mobilization of innate and adaptive immunity, predominantly CD4+ effector T cells, to achieve tumor control dependent on NKG2D. The combination induced important responses in patients with metastatic immune-cold tumors.This article is highlighted in the In This Issue feature, p. 1.

SPICA: Swiss portal for immune cell analysis.

Single-cell transcriptomics allows the study of immune cell heterogeneity at an unprecedented level of resolution. The Swiss portal for immune cell analysis (SPICA) is a web resource dedicated to the exploration and analysis of single-cell RNA-seq data of immune cells. In contrast to other single-cell databases, SPICA hosts curated, cell type-specific reference atlases that describe immune cell states at high resolution, and published single-cell datasets analysed in the context of these atlases. Additionally, users can privately analyse their own data in the context of existing atlases and contribute to the SPICA database. SPICA is available at https://spica.unil.ch.

UCell: Robust and scalable single-cell gene signature scoring.

UCell is an R package for evaluating gene signatures in single-cell datasets. UCell signature scores, based on the Mann-Whitney U statistic, are robust to dataset size and heterogeneity, and their calculation demands less computing time and memory than other available methods, enabling the processing of large datasets in a few minutes even on machines with limited computing power. UCell can be applied to any single-cell data matrix, and includes functions to directly interact with Seurat objects. The UCell package and documentation are available on GitHub at https://github.com/carmonalab/UCell.

Interpretation of T cell states from single-cell transcriptomics data using reference atlases.

Single-cell RNA sequencing (scRNA-seq) has revealed an unprecedented degree of immune cell diversity. However, consistent definition of cell subtypes and cell states across studies and diseases remains a major challenge. Here we generate reference T cell atlases for cancer and viral infection by multi-study integration, and develop ProjecTILs, an algorithm for reference atlas projection. In contrast to other methods, ProjecTILs allows not only accurate embedding of new scRNA-seq data into a reference without altering its structure, but also characterizing previously unknown cell states that "deviate" from the reference. ProjecTILs accurately predicts the effects of cell perturbations and identifies gene programs that are altered in different conditions and tissues. A meta-analysis of tumor-infiltrating T cells from several cohorts reveals a strong conservation of T cell subtypes between human and mouse, providing a consistent basis to describe T cell heterogeneity across studies, diseases, and species.

STACAS: Sub-Type Anchor Correction for Alignment in Seurat to integrate single-cell RNA-seq data.

SUMMARY: STACAS is a computational method for the identification of integration anchors in the Seurat environment, optimized for the integration of single-cell (sc) RNA-seq datasets that share only a subset of cell types. We demonstrate that by (i) correcting batch effects while preserving relevant biological variability across datasets, (ii) filtering aberrant integration anchors with a quantitative distance measure and (iii) constructing optimal guide trees for integration, STACAS can accurately align scRNA-seq datasets composed of only partially overlapping cell populations. AVAILABILITY AND IMPLEMENTATION: Source code and R package available at https://github.com/carmonalab/STACAS; Docker image available at https://hub.docker.com/repository/docker/mandrea1/stacas_demo.

Immunoinformatics: Predicting Peptide-MHC Binding.

Immunoinformatics is a discipline that applies methods of computer science to study and model the immune system. A fundamental question addressed by immunoinformatics is how to understand the rules of antigen presentation by MHC molecules to T cells, a process that is central to adaptive immune responses to infections and cancer. In the modern era of personalized medicine, the ability to model and predict which antigens can be presented by MHC is key to manipulating the immune system and designing strategies for therapeutic intervention. Since the MHC is both polygenic and extremely polymorphic, each individual possesses a personalized set of MHC molecules with different peptide-binding specificities, and collectively they present a unique individualized peptide imprint of the ongoing protein metabolism. Mapping all MHC allotypes is an enormous undertaking that cannot be achieved without a strong bioinformatics component. Computational tools for the prediction of peptide-MHC binding have thus become essential in most pipelines for T cell epitope discovery and an inescapable component of vaccine and cancer research. Here, we describe the development of several such tools, from pioneering efforts to the current state-of-the-art methods, that have allowed for accurate predictions of peptide binding of all MHC molecules, even including those that have not yet been characterized experimentally.

NNAlign_MA; MHC Peptidome Deconvolution for Accurate MHC Binding Motif Characterization and Improved T-cell Epitope Predictions.

The set of peptides presented on a cell's surface by MHC molecules is known as the immunopeptidome. Current mass spectrometry technologies allow for identification of large peptidomes, and studies have proven these data to be a rich source of information for learning the rules of MHC-mediated antigen presentation. Immunopeptidomes are usually poly-specific, containing multiple sequence motifs matching the MHC molecules expressed in the system under investigation. Motif deconvolution -the process of associating each ligand to its presenting MHC molecule(s)- is therefore a critical and challenging step in the analysis of MS-eluted MHC ligand data. Here, we describe NNAlign_MA, a computational method designed to address this challenge and fully benefit from large, poly-specific data sets of MS-eluted ligands. NNAlign_MA simultaneously performs the tasks of (1) clustering peptides into individual specificities; (2) automatic annotation of each cluster to an MHC molecule; and (3) training of a prediction model covering all MHCs present in the training set. NNAlign_MA was benchmarked on large and diverse data sets, covering class I and class II data. In all cases, the method was demonstrated to outperform state-of-the-art methods, effectively expanding the coverage of alleles for which accurate predictions can be made, resulting in improved identification of both eluted ligands and T-cell epitopes. Given its high flexibility and ease of use, we expect NNAlign_MA to serve as an effective tool to increase our understanding of the rules of MHC antigen presentation and guide the development of novel T-cell-based therapeutics.

IEDB-AR: immune epitope database-analysis resource in 2019.

The Immune Epitope Database Analysis Resource (IEDB-AR, http://tools.iedb.org/) is a companion website to the IEDB that provides computational tools focused on the prediction and analysis of B and T cell epitopes. All of the tools are freely available through the public website and many are also available through a REST API and/or a downloadable command-line tool. A virtual machine image of the entire site is also freely available for non-commercial use and contains most of the tools on the public site. Here, we describe the tools and functionalities that are available in the IEDB-AR, focusing on the 10 new tools that have been added since the last report in the 2012 NAR webserver edition. In addition, many of the tools that were already hosted on the site in 2012 have received updates to newest versions, including NetMHC, NetMHCpan, BepiPred and DiscoTope. Overall, this IEDB-AR update provides a substantial set of updated and novel features for epitope prediction and analysis.

MS-Rescue: A Computational Pipeline to Increase the Quality and Yield of Immunopeptidomics Experiments.

LC-MS/MS has become the standard platform for the characterization of immunopeptidomes, the collection of peptides naturally presented by major histocompatibility complex molecules to the cell surface. The protocols and algorithms used for immunopeptidomics data analysis are based on tools developed for traditional bottom-up proteomics that address the identification of peptides generated by tryptic digestion. Such algorithms are generally not tailored to the specific requirements of MHC ligand identification and, as a consequence, immunopeptidomics datasets suffer from dismissal of informative spectral information and high false discovery rates. Here, a new pipeline for the refinement of peptide-spectrum matches (PSM) is proposed, based on the assumption that immunopeptidomes contain a limited number of recurring peptide motifs, corresponding to MHC specificities. Sequence motifs are learned directly from the individual peptidome by training a prediction model on high-confidence PSMs. The model is then applied to PSM candidates with lower confidence, and sequences that score significantly higher than random peptides are rescued as likely true ligands. The pipeline is applied to MHC class I immunopeptidomes from three different species, and it is shown that it can increase the number of identified ligands by up to 20-30%, while effectively removing false positives and products of co-precipitation. Spectral validation using synthetic peptides confirms the identity of a large proportion of rescued ligands in the experimental peptidome.

Footprints of antigen processing boost MHC class II natural ligand predictions.

BACKGROUND: Major histocompatibility complex class II (MHC-II) molecules present peptide fragments to T cells for immune recognition. Current predictors for peptide to MHC-II binding are trained on binding affinity data, generated in vitro and therefore lacking information about antigen processing. METHODS: We generate prediction models of peptide to MHC-II binding trained with naturally eluted ligands derived from mass spectrometry in addition to peptide binding affinity data sets. RESULTS: We show that integrated prediction models incorporate identifiable rules of antigen processing. In fact, we observed detectable signals of protease cleavage at defined positions of the ligands. We also hypothesize a role of the length of the terminal ligand protrusions for trimming the peptide to the MHC presented ligand. CONCLUSIONS: The results of integrating binding affinity and eluted ligand data in a combined model demonstrate improved performance for the prediction of MHC-II ligands and T cell epitopes and foreshadow a new generation of improved peptide to MHC-II prediction tools accounting for the plurality of factors that determine natural presentation of antigens.

Predicting HLA CD4 Immunogenicity in Human Populations.

BACKGROUND: Prediction of T cell immunogenicity is a topic of considerable interest, both in terms of basic understanding of the mechanisms of T cells responses and in terms of practical applications. HLA binding affinity is often used to predict T cell epitopes, since HLA binding affinity is a key requisite for human T cell immunogenicity. However, immunogenicity at the population it is complicated by the high level of variability of HLA molecules, potential other factors beyond HLA as well as the frequent lack of HLA typing data. To overcome those issues, we explored an alternative approach to identify the common characteristics able to distinguish immunogenic peptides from non-recognized peptides. METHODS: Sets of dominant epitopes derived from peer-reviewed published papers were used in conjunction with negative peptides from the same experiments/donors to train neural networks and generate an "immunogenicity score." We also compared the performance of the immunogenicity score with previously described method for immunogenicity prediction based on HLA class II binding at the population level. RESULTS: The immunogenicity score was validated on a series of independent datasets derived from the published literature, representing 57 independent studies where immunogenicity in human populations was assessed by testing overlapping peptides spanning different antigens. Overall, these testing datasets corresponded to over 2,000 peptides and tested in over 1,600 different human donors. The 7-allele method prediction and the immunogenicity score were associated with similar performance [average area under the ROC curve (AUC) values of 0.703 and 0.702, respectively] while the combined methods reached an average AUC of 0.725. This increase in average AUC value is significant compared with the immunogenicity score (p = 0.0135) and a strong trend toward significance is observed when compared to the 7-allele method (p = 0.0938). The new immunogenicity score method is now freely available using CD4 T cell immunogenicity prediction tool on the Immune Epitope Database website (http://tools.iedb.org/CD4episcore). CONCLUSION: The new immunogenicity score predicts CD4 T cell immunogenicity at the population level starting from protein sequences and with no need for HLA typing. Its efficacy has been validated in the context of different antigen sources, ethnicities, and disparate techniques for epitope identification.