A framework for highly multiplexed dextramer mapping and prediction of T cell receptor sequences to antigen specificity

TCRCluster: a novel approach to T-cell receptor latent featurization and clustering using contrastive learning-guided two-stage variational autoencoders

T cells play a vital role in adaptive immunity by targeting pathogen-infected or cancerous cells, but predicting their specificity remains challenging. Encoding T-cell receptor (TCR) sequences into informative feature spaces is therefore crucial for advancing specificity prediction and downstream applications. For this, we developed a variational autoencoder (VAE)-based model trained on paired TCR α-β chain data, incorporating all six complementarity-determining regions. A semi-supervised 'two-stage VAE' framework, integrating cosine triplet loss and a classifier, was found to further refine peptide-specific latent representations, outperforming sequence-based methods in specificity prediction. Clustering analyses leveraging our VAE latent space were evaluated using K-means, agglomerative clustering, and a novel graph-based method. Agglomerative clustering achieved the most biologically relevant results, balancing cluster purity and retention despite noise in TCR specificity annotations. We extended these insights to evaluate TCR repertoire data. Across datasets, VAE-based models outperformed sequence-based methods, particularly in retention metrics, with notable improvements in the SARS-CoV-2 repertoire dataset. Moreover, the cancer repertoire analysis highlighted the generalizability of our approach, where the model displayed high performance despite minimal similarity between the training and test data. Collectively, these results demonstrate the potential of VAE-based latent representations to offer a robust framework for prediction, clustering, and repertoire analysis.

TRAP: a contrastive learning-enhanced framework for robust TCR-pMHC binding prediction with improved generalizability

The binding of T cell receptors (TCRs) to peptide-MHC I (pMHC) complexes is critical for triggering adaptive immune responses to potential health threats. Developing highly accurate machine learning (ML) models to predict TCR-pMHC binding could significantly accelerate immunotherapy advancements. However, existing ML models for TCR-pMHC binding prediction often underperform with unseen epitopes, severely limiting their applicability. We introduce TRAP, which leverages contrastive learning to enhance model performance by aligning structural and sequence features of pMHC with TCR sequences. TRAP outperforms previous state-of-the-art models in both random and unseen epitope scenarios, achieving an AUPR of 0.84 (a 22% improvement over the second-best model) and an AUC of 0.92 in the random scenario, and an AUC of 0.75 (almost 11% higher than the second-best model) in the unseen epitope scenario. Furthermore, TRAP demonstrates a noteworthy capability to diagnose potential issues of cross-reactivity between TCRs and similar epitopes. This highly robust performance makes it a suitable tool for large-scale predictions in real-world settings. A specific case study confirmed that TRAP can discover hit TCRs with binding free energies comparable to referenced experimental results. These findings highlight TRAP's potential for practical applications and its role as a powerful tool in developing TCR-based immunotherapies.

Feature selection enhances peptide binding predictions for TCR-specific interactions

Introduction

T-cell receptors (TCRs) play a critical role in the immune response by recognizing specific ligand peptides presented by major histocompatibility complex (MHC) molecules. Accurate prediction of peptide binding to TCRs is essential for advancing immunotherapy, vaccine design, and understanding mechanisms of autoimmune disorders.

Methods

This study presents a theoretical approach that explores the impact of feature selection techniques on enhancing the predictive accuracy of peptide binding models tailored for specific TCRs. To evaluate our approach across different TCR systems, we utilized a dataset that includes peptide libraries tested against three distinct murine TCRs. A broad range of physicochemical properties, including amino acid composition, dipeptide composition, and tripeptide features, were integrated into the machine learning-based feature selection framework to identify key properties contributing to binding affinity.

Results

Our analysis reveals that leveraging optimized feature subsets not only simplifies the model complexity but also enhances predictive performance, enabling more precise identification of TCR peptide interactions. The results of our feature selection method are consistent with findings from hybrid approaches that utilize both sequence and structural data as input as well as experimental data.

Discussion

Our theoretical approach highlights the role of feature selection in peptide-TCR interactions, providing a quantitative tool for uncovering the molecular mechanisms of the T-cell response and assisting in the design of more advanced targeted therapeutics.

The Evolving T Cell Receptor Recognition Code: The Rules Are More Like Guidelines

αβ T cell receptor (TCR) recognition of peptide-MHC complexes lies at the core of adaptive immunity, balancing specificity and cross-reactivity to facilitate effective antigen discrimination. Early structural studies established basic frameworks helpful for understanding and contextualizing TCR recognition and features such as peptide specificity and MHC restriction. However, the growing TCR structural database and studies launched from structural work continue to reveal exceptions to common assumptions and simplifications derived from earlier work. Here we explore our evolving understanding of TCR recognition, illustrating how structural and biophysical investigations regularly uncover complex phenomena that push against paradigms and expand our understanding of how TCRs bind to and discriminate between peptide/MHC complexes. We discuss the implications of these findings for basic, translational, and predictive immunology, including the challenges in accounting for the inherent adaptability, flexibility, and occasional biophysical sloppiness that characterize TCR recognition.

TPepRet: a deep learning model for characterizing T-cell receptors-antigen binding patterns

MOTIVATION

T-cell receptors (TCRs) elicit and mediate the adaptive immune response by recognizing antigenic peptides, a process pivotal for cancer immunotherapy, vaccine design, and autoimmune disease management. Understanding the intricate binding patterns between TCRs and peptides is critical for advancing these clinical applications. While several computational tools have been developed, they neglect the directional semantics inherent in sequence data, which are essential for accurately characterizing TCR-peptide interactions.

RESULTS

To address this gap, we develop TPepRet, an innovative model that integrates subsequence mining with semantic integration capabilities. TPepRet combines the strengths of the Bidirectional Gated Recurrent Unit (BiGRU) network for capturing bidirectional sequence dependencies with the Large Language Model framework to analyze subsequences and global sequences comprehensively, which enables TPepRet to accurately decipher the semantic binding relationship between TCRs and peptides. We have evaluated TPepRet to a range of challenging scenarios, including performance benchmarking against other tools using diverse datasets, analysis of peptide binding preferences, characterization of T cells clonal expansion, identification of true binder in complex environments, assessment of key binding sites through alanine scanning, validation against expression rates from large-scale datasets, and ability to screen SARS-CoV-2 TCRs. The comprehensive results suggest that TPepRet outperforms existing tools. We believe TPepRet will become an effective tool for understanding TCR-peptide binding in clinical treatment.

AVAILABILITY AND IMPLEMENTATION

The source code can be obtained from https://github.com/CSUBioGroup/TPepRet.git.

Unveiling cross-reactivity: implications for immune response modulation in cancer

Antigen recognition by CD8+ T-cell receptors (TCR) is crucial for immune responses to pathogens and tumors. TCRs are cross-reactive, a single TCR can recognize multiple peptide-Human Leukocyte Antigen (HLA) complexes. The study of cross-reactivity can support the development of therapies focusing on immune modulation, such as the expansion of pre-existing T-cell clones to fight pathogens and tumors. The peptide-HLA (pHLA) surface has previously been used to identify TCR cross-reactivities. In the present work, we sought to perform a comprehensive analysis of peptide-HLA by selecting thousands of human and viral epitopes. We profit from established docking models to identify features from different spatial perspectives of HLA-A*02:01, explore similarities between self and non-self epitopes, and list potential cross-reactive epitopes of therapeutic interest. A total of 2631 unique epitopes from representative viral proteins or human proteins were modeled. We were able to demonstrate that cross-reactive CDR3 sequences from public databases recognize epitopes with similar electrostatic potential, charge, and spatial location. Using data from published studies that measured T-cell reactivity to mutated epitopes, we observed a negative correlation between epitope dissimilarity and T-cell activation. Most analysed cancer epitopes were more similar to self epitopes, yet we identified features distinguishing those more similar to viral antigens. Finally, we enumerated potential cross-reactivities between tumoral and viral epitopes and highlighted some challenges in their identification for therapeutic use. Moreover, the thousands of peptide-HLA complexes generated in our work constitute a valuable resource to study T-cell cross-reactivity.

Feature Selection Enhances Peptide Binding Predictions for TCR-Specific Interactions

T-cell receptors (TCRs) play a critical role in the immune response by recognizing specific ligand peptides presented by major histocompatibility complex (MHC) molecules. Accurate prediction of peptide binding to TCRs is essential for advancing immunotherapy, vaccine design, and understanding mechanisms of autoimmune disorders. This study presents a novel theoretical method that explores the impact of feature selection techniques on enhancing the predictive accuracy of peptide binding models tailored for specific TCRs. To evaluate the universality of our approach across different TCR systems, we utilized a dataset that includes peptide libraries tested against three distinct murine TCRs. A broad range of physicochemical properties, including amino acid composition, dipeptide composition, and tripeptide features, were integrated into the machine learning-based feature selection framework to identify key features contributing to binding affinity. Our analysis reveals that leveraging optimized feature subsets not only simplifies the model complexity but also enhances predictive performance, enabling more precise identification of TCR-peptide interactions. The results of our feature selection method are consistent with findings from hybrid approaches that utilize both sequence and structural data as input as well as experimental data. Our theoretical approach highlights the role of feature selection in peptide-TCR interactions, providing a powerful tool for uncovering the molecular mechanisms of the T-cell response and assisting in the design of more advanced targeted therapeutics.

Computational Methods for Predicting Key Interactions in T Cell-Mediated Adaptive Immunity

The adaptive immune system recognizes pathogen- and cancer-specific features and is endowed with memory, enabling it to respond quickly and efficiently to repeated encounters with the same antigens. T cells play a central role in the adaptive immune system by directly targeting intracellular pathogens and helping to activate B cells to secrete antibodies. Several fundamental protein interactions-including those between major histocompatibility complex (MHC) proteins and antigen-derived peptides as well as between T cell receptors and peptide-MHC complexes-underlie the ability of T cells to recognize antigens with great precision. Computational approaches to predict these interactions are increasingly being used for medically relevant applications, including vaccine design and prediction of patient response to cancer immunotherapies. We provide computational researchers with an accessible introduction to the adaptive immune system, review computational approaches to predict the key protein interactions underlying T cell-mediated adaptive immunity, and highlight remaining challenges.

Single-cell genomics-based immune and disease monitoring in blood malignancies

Achieving long-term disease control using therapeutic immunomodulation is a long-standing concept with a strong tradition in blood malignancies. Besides allogeneic hematopoietic stem cell transplantation that continues to provide potentially curative treatment for otherwise challenging diagnoses, recent years have seen impressive progress in immunotherapies for leukemias and lymphomas with immune checkpoint blockade, bispecific monoclonal antibodies, and CAR T cell therapies. Despite their success, non-response, relapse, and immune toxicities remain frequent, thus prioritizing the elucidation of the underlying mechanisms and identifying predictive biomarkers. The increasing availability of single-cell genomic tools now provides a system's immunology view to resolve the molecular and cellular mechanisms of immunotherapies at unprecedented resolution. Here, we review recent studies that leverage these technological advancements for tracking immune responses, the emergence of immune resistance, and toxicities. As single-cell immune monitoring tools evolve and become more accessible, we expect their wide adoption for routine clinical applications to catalyze more precise therapeutic steering of personal immune responses.

Artificial intelligence and neoantigens: paving the path for precision cancer immunotherapy

Cancer immunotherapy has witnessed rapid advancement in recent years, with a particular focus on neoantigens as promising targets for personalized treatments. The convergence of immunogenomics, bioinformatics, and artificial intelligence (AI) has propelled the development of innovative neoantigen discovery tools and pipelines. These tools have revolutionized our ability to identify tumor-specific antigens, providing the foundation for precision cancer immunotherapy. AI-driven algorithms can process extensive amounts of data, identify patterns, and make predictions that were once challenging to achieve. However, the integration of AI comes with its own set of challenges, leaving space for further research. With particular focus on the computational approaches, in this article we have explored the current landscape of neoantigen prediction, the fundamental concepts behind, the challenges and their potential solutions providing a comprehensive overview of this rapidly evolving field.

Designing meaningful continuous representations of T cell receptor sequences with deep generative models

T Cell Receptor (TCR) antigen binding underlies a key mechanism of the adaptive immune response yet the vast diversity of TCRs and the complexity of protein interactions limits our ability to build useful low dimensional representations of TCRs. To address the current limitations in TCR analysis we develop a capacity-controlled disentangling variational autoencoder trained using a dataset of approximately 100 million TCR sequences, that we name TCR-VALID. We design TCR-VALID such that the model representations are low-dimensional, continuous, disentangled, and sufficiently informative to provide high-quality TCR sequence de novo generation. We thoroughly quantify these properties of the representations, providing a framework for future protein representation learning in low dimensions. The continuity of TCR-VALID representations allows fast and accurate TCR clustering and is benchmarked against other state-of-the-art TCR clustering tools and pre-trained language models.

RACER-m leverages structural features for sparse T cell specificity prediction

Reliable prediction of T cell specificity against antigenic signatures is a formidable task, complicated by the immense diversity of T cell receptor and antigen sequence space and the resulting limited availability of training sets for inferential models. Recent modeling efforts have demonstrated the advantage of incorporating structural information to overcome the need for extensive training sequence data, yet disentangling the heterogeneous TCR-antigen interface to accurately predict MHC-allele-restricted TCR-peptide interactions has remained challenging. Here, we present RACER-m, a coarse-grained structural model leveraging key biophysical information from the diversity of publicly available TCR-antigen crystal structures. Explicit inclusion of structural content substantially reduces the required number of training examples and maintains reliable predictions of TCR-recognition specificity and sensitivity across diverse biological contexts. Our model capably identifies biophysically meaningful point-mutant peptides that affect binding affinity, distinguishing its ability in predicting TCR specificity of point-mutants from alternative sequence-based methods. Its application is broadly applicable to studies involving both closely related and structurally diverse TCR-peptide pairs.

Deep learning predictions of TCR-epitope interactions reveal epitope-specific chains in dual alpha T cells

T cells have the ability to eliminate infected and cancer cells and play an essential role in cancer immunotherapy. T cell activation is elicited by the binding of the T cell receptor (TCR) to epitopes displayed on MHC molecules, and the TCR specificity is determined by the sequence of its α and β chains. Here, we collect and curate a dataset of 17,715 αβTCRs interacting with dozens of class I and class II epitopes. We use this curated data to develop MixTCRpred, an epitope-specific TCR-epitope interaction predictor. MixTCRpred accurately predicts TCRs recognizing several viral and cancer epitopes. MixTCRpred further provides a useful quality control tool for multiplexed single-cell TCR sequencing assays of epitope-specific T cells and pinpoints a substantial fraction of putative contaminants in public databases. Analysis of epitope-specific dual α T cells demonstrates that MixTCRpred can identify α chains mediating epitope recognition. Applying MixTCRpred to TCR repertoires from COVID-19 patients reveals enrichment of clonotypes predicted to bind an immunodominant SARS-CoV-2 epitope. Overall, MixTCRpred provides a robust tool to predict TCRs interacting with specific epitopes and interpret TCR-sequencing data from both bulk and epitope-specific T cells.

Enhancing TCR specificity predictions by combined pan- and peptide-specific training, loss-scaling, and sequence similarity integration

Predicting the interaction between Major Histocompatibility Complex (MHC) class I-presented peptides and T-cell receptors (TCR) holds significant implications for vaccine development, cancer treatment, and autoimmune disease therapies. However, limited paired-chain TCR data, skewed towards well-studied epitopes, hampers the development of pan-specific machine-learning (ML) models. Leveraging a larger peptide-TCR dataset, we explore various alterations to the ML architectures and training strategies to address data imbalance. This leads to an overall improved performance, particularly for peptides with scant TCR data. However, challenges persist for unseen peptides, especially those distant from training examples. We demonstrate that such ML models can be used to detect potential outliers, which when removed from training, leads to augmented performance. Integrating pan-specific and peptide-specific models alongside with similarity-based predictions, further improves the overall performance, especially when a low false positive rate is desirable. In the context of the IMMREP22 benchmark, this modeling framework attained state-of-the-art performance. Moreover, combining these strategies results in acceptable predictive accuracy for peptides characterized with as little as 15 positive TCRs. This observation places great promise on rapidly expanding the peptide covering of the current models for predicting TCR specificity. The NetTCR 2.2 model incorporating these advances is available on GitHub (https://github.com/mnielLab/NetTCR-2.2) and as a web server at https://services.healthtech.dtu.dk/services/NetTCR-2.2/.

A comparison of clustering models for inference of T cell receptor antigen specificity

The vast potential sequence diversity of TCRs and their ligands has presented an historic barrier to computational prediction of TCR epitope specificity, a holy grail of quantitative immunology. One common approach is to cluster sequences together, on the assumption that similar receptors bind similar epitopes. Here, we provide the first independent evaluation of widely used clustering algorithms for TCR specificity inference, observing some variability in predictive performance between models, and marked differences in scalability. Despite these differences, we find that different algorithms produce clusters with high degrees of similarity for receptors recognising the same epitope. Our analysis strengthens the case for use of clustering models to identify signals of common specificity from large repertoires, whilst highlighting scope for improvement of complex models over simple comparators.

Development and use of machine learning algorithms in vaccine target selection

Computer-aided discovery of vaccine targets has become a cornerstone of rational vaccine design. In this article, I discuss how Machine Learning (ML) can inform and guide key computational steps in rational vaccine design concerned with the identification of B and T cell epitopes and correlates of protection. I provide examples of ML models, as well as types of data and predictions for which they are built. I argue that interpretable ML has the potential to improve the identification of immunogens also as a tool for scientific discovery, by helping elucidate the molecular processes underlying vaccine-induced immune responses. I outline the limitations and challenges in terms of data availability and method development that need to be addressed to bridge the gap between advances in ML predictions and their translational application to vaccine design.

pan-MHC and cross-Species Prediction of T Cell Receptor-Antigen Binding

Profiling the binding of T cell receptors (TCRs) of T cells to antigenic peptides presented by MHC proteins is one of the most important unsolved problems in modern immunology. Experimental methods to probe TCR-antigen interactions are slow, labor-intensive, costly, and yield moderate throughput. To address this problem, we developed pMTnet-omni, an Artificial Intelligence (AI) system based on hybrid protein sequence and structure information, to predict the pairing of TCRs of αβ T cells with peptide-MHC complexes (pMHCs). pMTnet-omni is capable of handling peptides presented by both class I and II pMHCs, and capable of handling both human and mouse TCR-pMHC pairs, through information sharing enabled this hybrid design. pMTnet-omni achieves a high overall Area Under the Curve of Receiver Operator Characteristics (AUROC) of 0.888, which surpasses competing tools by a large margin. We showed that pMTnet-omni can distinguish binding affinity of TCRs with similar sequences. Across a range of datasets from various biological contexts, pMTnet-omni characterized the longitudinal evolution and spatial heterogeneity of TCR-pMHC interactions and their functional impact. We successfully developed a biomarker based on pMTnet-omni for predicting immune-related adverse events of immune checkpoint inhibitor (ICI) treatment in a cohort of 57 ICI-treated patients. pMTnet-omni represents a major advance towards developing a clinically usable AI system for TCR-pMHC pairing prediction that can aid the design and implementation of TCR-based immunotherapeutics.

EPIC-TRACE: predicting TCR binding to unseen epitopes using attention and contextualized embeddings

MOTIVATION

T cells play an essential role in adaptive immune system to fight pathogens and cancer but may also give rise to autoimmune diseases. The recognition of a peptide-MHC (pMHC) complex by a T cell receptor (TCR) is required to elicit an immune response. Many machine learning models have been developed to predict the binding, but generalizing predictions to pMHCs outside the training data remains challenging.

RESULTS

We have developed a new machine learning model that utilizes information about the TCR from both α and β chains, epitope sequence, and MHC. Our method uses ProtBERT embeddings for the amino acid sequences of both chains and the epitope, as well as convolution and multi-head attention architectures. We show the importance of each input feature as well as the benefit of including epitopes with only a few TCRs to the training data. We evaluate our model on existing databases and show that it compares favorably against other state-of-the-art models.

AVAILABILITY AND IMPLEMENTATION

https://github.com/DaniTheOrange/EPIC-TRACE.

Benchmarking data-driven filtering for denoising of TCRpMHC single-cell data

Pairing of the T cell receptor (TCR) with its cognate peptide-MHC (pMHC) is a cornerstone in T cell-mediated immunity. Recently, single-cell sequencing coupled with DNA-barcoded MHC multimer staining has enabled high-throughput studies of T cell specificities. However, the immense variability of TCR-pMHC interactions combined with the relatively low signal-to-noise ratio in the data generated using current technologies are complicating these studies. Several approaches have been proposed for denoising single-cell TCR-pMHC specificity data. Here, we present a benchmark evaluating two such denoising methods, ICON and ITRAP. We applied and evaluated the methods on publicly available immune profiling data provided by 10x Genomics. We find that both methods identified approximately 75% of the raw data as noise. We analyzed both internal metrics developed for the purpose and performance on independent data using machine learning methods trained on the raw and denoised 10x data. We find an increased signal-to-noise ratio comparing the denoised to the raw data for both methods, and demonstrate an overall superior performance of the ITRAP method in terms of both data consistency and performance. In conclusion, this study demonstrates that Improving the data quality from high throughput studies of TCRpMHC-specificity by denoising is paramount in increasing our understanding of T cell-mediated immunity.

TAPIR: a T-cell receptor language model for predicting rare and novel targets

T-cell receptors (TCRs) are involved in most human diseases, but linking their sequences with their targets remains an unsolved grand challenge in the field. In this study, we present TAPIR (T-cell receptor and Peptide Interaction Recognizer), a T-cell receptor (TCR) language model that predicts TCR-target interactions, with a focus on novel and rare targets. TAPIR employs deep convolutional neural network (CNN) encoders to process TCR and target sequences across flexible representations (e.g., beta-chain only, unknown MHC allele, etc.) and learns patterns of interactivity via several training tasks. This flexibility allows TAPIR to train on more than 50k either paired (alpha and beta chain) or unpaired TCRs (just alpha or beta chain) from public and proprietary databases against 1933 unique targets. TAPIR demonstrates state-of-the-art performance when predicting TCR interactivity against common benchmark targets and is the first method to demonstrate strong performance when predicting TCR interactivity against novel targets, where no examples are provided in training. TAPIR is also capable of predicting TCR interaction against MHC alleles in the absence of target information. Leveraging these capabilities, we apply TAPIR to cancer patient TCR repertoires and identify and validate a novel and potent anti-cancer T-cell receptor against a shared cancer neoantigen target (PIK3CA H1047L). We further show how TAPIR, when extended with a generative neural network, is capable of directly designing T-cell receptor sequences that interact with a target of interest.

A transfer-learning approach to predict antigen immunogenicity and T-cell receptor specificity

Antigen immunogenicity and the specificity of binding of T-cell receptors to antigens are key properties underlying effective immune responses. Here we propose diffRBM, an approach based on transfer learning and Restricted Boltzmann Machines, to build sequence-based predictive models of these properties. DiffRBM is designed to learn the distinctive patterns in amino-acid composition that, on the one hand, underlie the antigen's probability of triggering a response, and on the other hand the T-cell receptor's ability to bind to a given antigen. We show that the patterns learnt by diffRBM allow us to predict putative contact sites of the antigen-receptor complex. We also discriminate immunogenic and non-immunogenic antigens, antigen-specific and generic receptors, reaching performances that compare favorably to existing sequence-based predictors of antigen immunogenicity and T-cell receptor specificity.

Quantitative approaches for decoding the specificity of the human T cell repertoire

T cell receptor (TCR)-peptide-major histocompatibility complex (pMHC) interactions play a vital role in initiating immune responses against pathogens, and the specificity of TCRpMHC interactions is crucial for developing optimized therapeutic strategies. The advent of high-throughput immunological and structural evaluation of TCR and pMHC has provided an abundance of data for computational approaches that aim to predict favorable TCR-pMHC interactions. Current models are constructed using information on protein sequence, structures, or a combination of both, and utilize a variety of statistical learning-based approaches for identifying the rules governing specificity. This review examines the current theoretical, computational, and deep learning approaches for identifying TCR-pMHC recognition pairs, placing emphasis on each method's mathematical approach, predictive performance, and limitations.

DeepAIR: A deep learning framework for effective integration of sequence and 3D structure to enable adaptive immune receptor analysis

Structural docking between the adaptive immune receptors (AIRs), including T cell receptors (TCRs) and B cell receptors (BCRs), and their cognate antigens are one of the most fundamental processes in adaptive immunity. However, current methods for predicting AIR-antigen binding largely rely on sequence-derived features of AIRs, omitting the structure features that are essential for binding affinity. In this study, we present a deep learning framework, termed DeepAIR, for the accurate prediction of AIR-antigen binding by integrating both sequence and structure features of AIRs. DeepAIR achieves a Pearson's correlation of 0.813 in predicting the binding affinity of TCR, and a median area under the receiver-operating characteristic curve (AUC) of 0.904 and 0.942 in predicting the binding reactivity of TCR and BCR, respectively. Meanwhile, using TCR and BCR repertoire, DeepAIR correctly identifies every patient with nasopharyngeal carcinoma and inflammatory bowel disease in test data. Thus, DeepAIR improves the AIR-antigen binding prediction that facilitates the study of adaptive immunity.

Can we predict T cell specificity with digital biology and machine learning?

Recent advances in machine learning and experimental biology have offered breakthrough solutions to problems such as protein structure prediction that were long thought to be intractable. However, despite the pivotal role of the T cell receptor (TCR) in orchestrating cellular immunity in health and disease, computational reconstruction of a reliable map from a TCR to its cognate antigens remains a holy grail of systems immunology. Current data sets are limited to a negligible fraction of the universe of possible TCR-ligand pairs, and performance of state-of-the-art predictive models wanes when applied beyond these known binders. In this Perspective article, we make the case for renewed and coordinated interdisciplinary effort to tackle the problem of predicting TCR-antigen specificity. We set out the general requirements of predictive models of antigen binding, highlight critical challenges and discuss how recent advances in digital biology such as single-cell technology and machine learning may provide possible solutions. Finally, we describe how predicting TCR specificity might contribute to our understanding of the broader puzzle of antigen immunogenicity.

BERTrand-peptide:TCR binding prediction using Bidirectional Encoder Representations from Transformers augmented with random TCR pairing

MOTIVATION

The advent of T-cell receptor (TCR) sequencing experiments allowed for a significant increase in the amount of peptide:TCR binding data available and a number of machine-learning models appeared in recent years. High-quality prediction models for a fixed epitope sequence are feasible, provided enough known binding TCR sequences are available. However, their performance drops significantly for previously unseen peptides.

RESULTS

We prepare the dataset of known peptide:TCR binders and augment it with negative decoys created using healthy donors' T-cell repertoires. We employ deep learning methods commonly applied in Natural Language Processing to train part a peptide:TCR binding model with a degree of cross-peptide generalization (0.69 AUROC). We demonstrate that BERTrand outperforms the published methods when evaluated on peptide sequences not used during model training.

AVAILABILITY AND IMPLEMENTATION

The datasets and the code for model training are available at https://github.com/SFGLab/bertrand.

The T cell receptor sequence influences the likelihood of T cell memory formation

T cell differentiation depends on activation through the T cell receptor (TCR), whose amino acid sequence varies cell to cell. Particular TCR amino acid sequences nearly guarantee Mucosal-Associated Invariant T (MAIT) and Natural Killer T (NKT) cell fates. To comprehensively define how TCR amino acids affects all T cell fates, we analyze the paired αβTCR sequence and transcriptome of 819,772 single cells. We find that hydrophobic CDR3 residues promote regulatory T cell transcriptional states in both the CD8 and CD4 lineages. Most strikingly, we find a set of TCR sequence features, concentrated in CDR2α, that promotes positive selection in the thymus as well as transition from naïve to memory in the periphery. Even among T cells that recognize the same antigen, these TCR sequence features help to explain which T cells form immunological memory, which is essential for effective pathogen response.

High-Throughput Screening of Functional Neo-Antigens and Their Specific T-Cell Receptors via the Jurkat Reporter System Combined with Droplet Microfluidics

T-cell receptor (TCR)-engineered T cells can precisely recognize a broad repertoire of targets derived from both intracellular and surface proteins of tumor cells. TCR-T adoptive cell therapy has shown safety and promising efficacy in solid tumor immunotherapy. However, antigen-specific functional TCR screening is time-consuming and expensive, which limits its application clinically. Here, we developed a novel integrated antigen-TCR screening platform based on droplet microfluidic technology, enabling high-throughput peptide-major histocompatibility complex (pMHC)-to-TCR paired screening with a high sensitivity and low background signal. We introduced DNA barcoding technology to label peptide antigen candidate-loaded antigen-presenting cells and Jurkat reporter cells to check the specificity of pMHC-TCR candidates. Coupled with the next-generation sequencing pipeline, interpretation of the DNA barcodes and the gene expression level of the Jurkat T-cell activation pathway provided a clear peptide-MHC-TCR recognition relationship. Our proof-of-principle study demonstrates that the platform could achieve pMHC-TCR paired high-throughput screening, which is expected to be used in the cross-reactivity and off-target high-throughput paired testing of candidate pMHC-TCRs in clinical applications.

Rapid TCR:Epitope Ranker (RAPTER): a primary human T cell reactivity screening assay pairing epitope and TCR at single cell resolution

Identifying epitopes that T cells respond to is critical for understanding T cell-mediated immunity. Traditional multimer and other single cell assays often require large blood volumes and/or expensive HLA-specific reagents and provide limited phenotypic and functional information. Here, we present the Rapid TCR:Epitope Ranker (RAPTER) assay, a single cell RNA sequencing (scRNA-SEQ) method that uses primary human T cells and antigen presenting cells (APCs) to assess functional T cell reactivity. Using hash-tag oligonucleotide (HTO) coding and T cell activation-induced markers (AIM), RAPTER defines paired epitope specificity and TCR sequence and can include RNA- and protein-level T cell phenotype information. We demonstrate that RAPTER identified specific reactivities to viral and tumor antigens at sensitivities as low as 0.15% of total CD8⁺ T cells, and deconvoluted low-frequency circulating HPV16-specific T cell clones from a cervical cancer patient. The specificities of TCRs identified by RAPTER for MART1, EBV, and influenza epitopes were functionally confirmed in vitro. In summary, RAPTER identifies low-frequency T cell reactivities using primary cells from low blood volumes, and the resulting paired TCR:ligand information can directly enable immunogenic antigen selection from limited patient samples for vaccine epitope inclusion, antigen-specific TCR tracking, and TCR cloning for further therapeutic development.

epiTCR: a highly sensitive predictor for TCR-peptide binding

MOTIVATION

Predicting the binding between T-cell receptor (TCR) and peptide presented by human leucocyte antigen molecule is a highly challenging task and a key bottleneck in the development of immunotherapy. Existing prediction tools, despite exhibiting good performance on the datasets they were built with, suffer from low true positive rates when used to predict epitopes capable of eliciting T-cell responses in patients. Therefore, an improved tool for TCR-peptide prediction built upon a large dataset combining existing publicly available data is still needed.

RESULTS

We collected data from five public databases (IEDB, TBAdb, VDJdb, McPAS-TCR, and 10X) to form a dataset of >3 million TCR-peptide pairs, 3.27% of which were binding interactions. We proposed epiTCR, a Random Forest-based method dedicated to predicting the TCR-peptide interactions. epiTCR used simple input of TCR CDR3β sequences and antigen sequences, which are encoded by flattened BLOSUM62. epiTCR performed with area under the curve (0.98) and higher sensitivity (0.94) than other existing tools (NetTCR, Imrex, ATM-TCR, and pMTnet), while maintaining comparable prediction specificity (0.9). We identified seven epitopes that contributed to 98.67% of false positives predicted by epiTCR and exerted similar effects on other tools. We also demonstrated a considerable influence of peptide sequences on prediction, highlighting the need for more diverse peptides in a more balanced dataset. In conclusion, epiTCR is among the most well-performing tools, thanks to the use of combined data from public sources and its use will contribute to the quest in identifying neoantigens for precision cancer immunotherapy.

AVAILABILITY AND IMPLEMENTATION

epiTCR is available on GitHub (https://github.com/ddiem-ri-4D/epiTCR).

Improved T cell receptor antigen pairing through data-driven filtering of sequencing information from single cells

Novel single-cell-based technologies hold the promise of matching T cell receptor (TCR) sequences with their cognate peptide-MHC recognition motif in a high-throughput manner. Parallel capture of TCR transcripts and peptide-MHC is enabled through the use of reagents labeled with DNA barcodes. However, analysis and annotation of such single-cell sequencing (SCseq) data are challenged by dropout, random noise, and other technical artifacts that must be carefully handled in the downstream processing steps. We here propose a rational, data-driven method termed ITRAP (improved T cell Receptor Antigen Paring) to deal with these challenges, filtering away likely artifacts, and enable the generation of large sets of TCR-pMHC sequence data with a high degree of specificity and sensitivity, thus outputting the most likely pMHC target per T cell. We have validated this approach across 10 different virus-specific T cell responses in 16 healthy donors. Across these samples, we have identified up to 1494 high-confident TCR-pMHC pairs derived from 4135 single cells.

Entropic analysis of antigen-specific CDR3 domains identifies essential binding motifs shared by CDR3s with different antigen specificities

Antigen-specific T cell receptor (TCR) sequences can have prognostic, predictive, and therapeutic value, but decoding the specificity of TCR recognition remains challenging. Unlike DNA strands that base pair, TCRs bind to their targets with different orientations and different lengths, which complicates comparisons. We present scanning parametrized by normalized TCR length (SPAN-TCR) to analyze antigen-specific TCR CDR3 sequences and identify patterns driving TCR-pMHC specificity. Using entropic analysis, SPAN-TCR identifies 2-mer motifs that decrease the diversity (entropy) of CDR3s. These motifs are the most common patterns that can predict CDR3 composition, and we identify "essential" motifs that decrease entropy in the same CDR3 α or β chain containing the 2-mer, and "super-essential" motifs that decrease entropy in both chains. Molecular dynamics analysis further suggests that these motifs may play important roles in binding. We then employ SPAN-TCR to resolve similarities in TCR repertoires against different antigens using public databases of TCR sequences.

Antigen-specificity measurements are the key to understanding T cell responses

Antigen-specific T cells play a central role in the adaptive immune response and come in a wide range of phenotypes. T cell receptors (TCRs) mediate the antigen-specificities found in T cells. Importantly, high-throughput TCR sequencing provides a fingerprint which allows tracking of specific T cells and their clonal expansion in response to particular antigens. As a result, many studies have leveraged TCR sequencing in an attempt to elucidate the role of antigen-specific T cells in various contexts. Here, we discuss the published approaches to studying antigen-specific T cells and their specific TCR repertoire. Further, we discuss how these methods have been applied to study the TCR repertoire in various diseases in order to characterize the antigen-specific T cells involved in the immune control of disease.

Measures of epitope binding degeneracy from T cell receptor repertoires

Adaptive immunity is driven by specific binding of hypervariable receptors to diverse molecular targets. The sequence diversity of receptors and targets are both individually known but because multiple receptors can recognize the same target, a measure of the effective "functional" diversity of the human immune system has remained elusive. Here, we show that sequence near-coincidences within T cell receptors that bind specific epitopes provide a new window into this problem and allow the quantification of how binding probability covaries with sequence. We find that near-coincidence statistics within epitope-specific repertoires imply a measure of binding degeneracy to amino acid changes in receptor sequence that is consistent across disparate experiments. Paired data on both chains of the heterodimeric receptor are particularly revealing since simultaneous near-coincidences are rare and we show how they can be exploited to estimate the number of epitope responses that created the memory compartment. In addition, we find that paired-chain coincidences are strongly suppressed across donors with different human leukocyte antigens, evidence for a central role of antigen-driven selection in making paired chain receptors public. These results demonstrate the power of coincidence analysis to reveal the sequence determinants of epitope binding in receptor repertoires.

Inferring the T cell repertoire dynamics of healthy individuals

The adaptive immune system is a diverse ecosystem that responds to pathogens by selecting cells with specific receptors. While clonal expansion in response to particular immune challenges has been extensively studied, we do not know the neutral dynamics that drive the immune system in the absence of strong stimuli. Here, we learn the parameters that underlie the clonal dynamics of the T cell repertoire in healthy individuals of different ages, by applying Bayesian inference to longitudinal immune repertoire sequencing (RepSeq) data. Quantifying the experimental noise accurately for a given RepSeq technique allows us to disentangle real changes in clonal frequencies from noise. We find that the data are consistent with clone sizes following a geometric Brownian motion and show that its predicted steady state is in quantitative agreement with the observed power-law behavior of the clone-size distribution. The inferred turnover time scale of the repertoire increases with patient age and depends on the clone size in some individuals.

NetTCR-2.1: Lessons and guidance on how to develop models for TCR specificity predictions

T cell receptors (TCR) define the specificity of T cells and are responsible for their interaction with peptide antigen targets presented in complex with major histocompatibility complex (MHC) molecules. Understanding the rules underlying this interaction hence forms the foundation for our understanding of basic adaptive immunology. Over the last decade, efforts have been dedicated to developing assays for high throughput identification of peptide-specific TCRs. Based on such data, several computational methods have been proposed for predicting the TCR-pMHC interaction. The general conclusion from these studies is that the prediction of TCR interactions with MHC-peptide complexes remains highly challenging. Several reasons form the basis for this including scarcity and quality of data, and ill-defined modeling objectives imposed by the high redundancy of the available data. In this work, we propose a framework for dealing with this redundancy, allowing us to address essential questions related to the modeling of TCR specificity including the use of peptide- versus pan-specific models, how to best define negative data, and the performance impact of integrating of CDR1 and 2 loops. Further, we illustrate how and why it is strongly recommended to include simple similarity-based modeling approaches when validating an improved predictive power of machine learning models, and that such validation should include a performance evaluation as a function of "distance" to the training data, to quantify the potential for generalization of the proposed model. The conclusion of the work is that, given current data, TCR specificity is best modeled using peptide-specific approaches, integrating information from all 6 CDR loops, and with negative data constructed from a combination of true and mislabeled negatives. Comparing such machine learning models to similarity-based approaches demonstrated an increased performance gain of the former as the "distance" to the training data was increased; thus demonstrating an improved generalization ability of the machine learning-based approaches. We believe these results demonstrate that the outlined modeling framework and proposed evaluation strategy form a solid basis for investigating the modeling of TCR specificities and that adhering to such a framework will allow for faster progress within the field. The final devolved model, NetTCR-2.1, is available at https://services.healthtech.dtu.dk/service.php?NetTCR-2.1.

Rapid identification of CMV-specific TCRs via reverse TCR cloning system based on bulk TCR repertoire data

Advances in next-generation sequencing (NGS) have improved the resolution of T-cell receptor (TCR) repertoire analysis, and recent single-cell sequencing has made it possible to obtain information about TCR pairs. In our previous study, cytomegalovirus (CMV) pp65-specific T-cell response restricted by a single human leukocyte antigen (HLA) class I allotype was observed in an individual. Therefore, to effectively clone an antigen-specific TCR from these T cells, we developed a TCR cloning system that does not require a single cell level. First, we established the improved Jurkat reporter cell line, which was TCRαβ double knock-out and expressed CD8αβ molecules. Furthermore, functional TCRs were directly obtained by reverse TCR cloning using unique CDR3-specific PCR primers after bulk TCR sequencing of activation marker-positive CD8 T cells by NGS. A total of 15 TCRα and 14 TCRβ strands were successfully amplified by PCR from cDNA of 4-1BB-positive CD8 T cells restricted by HLA-A*02:01, HLA-A*02:06, HLA-B*07:02, and HLA-B*40:06. The panels with combinations of TCRα and TCRβ genes were investigated using Jurkat reporter cell line and artificial antigen-presenting cells (APCs). In two TCR pairs restricted by HLA-A*02:01, one TCR pair by HLA-A*02:06, four TCR pairs by HLA-B*07:02, and one TCR pair by HLA-B*40:06, their specificity and affinity were confirmed. The TCR pair of A*02:01/1-1 showed alloreactivity to HLA-A*02:06. The one TCR pair showed a higher response to the naturally processed antigen than that of the peptide pool. This reverse TCR cloning system will not only provide functional information to TCR repertoire analysis by NGS but also help in the development of TCR-T therapy.

A guide to systems-level immunomics

The immune system is highly complex and distributed throughout an organism, with hundreds to thousands of cell states existing in parallel with diverse molecular pathways interacting in a highly dynamic and coordinated fashion. Although the characterization of individual genes and molecules is of the utmost importance for understanding immune-system function, high-throughput, high-resolution omics technologies combined with sophisticated computational modeling and machine-learning approaches are creating opportunities to complement standard immunological methods with new insights into immune-system dynamics. Like systems immunology itself, immunology researchers must take advantage of these technologies and form their own diverse networks, connecting with researchers from other disciplines. This Review is an introduction and 'how-to guide' for immunologists with no particular experience in the field of omics but with the intention to learn about and apply these systems-level approaches, and for immunologists who want to make the most of interdisciplinary networks.

Unraveling the Complexity of Liver Disease One Cell at a Time

The human liver is a complex organ made up of multiple specialized cell types that carry out key physiological functions. An incomplete understanding of liver biology limits our ability to develop therapeutics to prevent chronic liver diseases, liver cancers, and death as a result of organ failure. Recently, single-cell modalities have expanded our understanding of the cellular phenotypic heterogeneity and intercellular cross-talk in liver health and disease. This review summarizes these findings and looks forward to highlighting new avenues for the application of single-cell genomics to unravel unknown pathogenic pathways and disease mechanisms for the development of new therapeutics targeting liver pathology. As these technologies mature, their integration into clinical data analysis will aid in patient stratification and in developing treatment plans for patients suffering from liver disease.

Designing Personalized Antigen-Specific Immunotherapies for Autoimmune Diseases-The Case for Using Ignored Target Cell Antigen Determinants

We have proposed that antigen-specific immunotherapies (ASIs) for autoimmune diseases could be enhanced by administering target cell antigen epitopes (determinants) that are immunogenic but ignored by autoreactive T cells because these determinants may have large pools of naïve cognate T cells available for priming towards regulatory responses. Here, we identified an immunogenic preproinsulin determinant (PPI_L4-20) that was ignored by autoimmune responses in type 1 diabetes (T1D)-prone NOD mice. The size of the PPI_L4-20-specific splenic naive T cell pool gradually increased from 2–12 weeks in age and remained stable thereafter, while that of the major target determinant insulin B-chain_9-23 decreased greatly after 12 weeks in age, presumably due to recruitment into the autoimmune response. In 15–16 week old mice, insulin B-chain_9-23/alum immunization induced modest-low level of splenic T cell IL-10 and IL-4 responses, little or no spreading of these responses, and boosted IFNγ responses to itself and other autoantigens. In contrast, PPI_L4-20/alum treatment induced robust IL-10 and IL-4 responses, which spread to other autoantigens and increased the frequency of splenic IL-10-secreting Treg and Tr-1-like cells, without boosting IFNγ responses to ß-cell autoantigens. In newly diabetic NOD mice, PPI_L4-20, but not insulin B-chain_9-23 administered intraperitoneally (with alum) or intradermally (as soluble antigen) supplemented with oral GABA induced long-term disease remission. We discuss the potential of personalized ASIs that are based on an individual’s naïve autoantigen-reactive T cell pools and the use of HLA-appropriate ignored autoantigen determinants to safely enhance the efficacy of ASIs.

Analysis of T-Cell Receptor Repertoire in Transplantation: Fingerprint of T Cell-mediated Alloresponse

T cells play a key role in determining allograft function by mediating allogeneic immune responses to cause rejection, and recent work pointed their role in mediating tolerance in transplantation. The unique T-cell receptor (TCR) expressed on the surface of each T cell determines the antigen specificity of the cell and can be the specific fingerprint for identifying and monitoring. Next-generation sequencing (NGS) techniques provide powerful tools for deep and high-throughput TCR profiling, and facilitate to depict the entire T cell repertoire profile and trace antigen-specific T cells in circulation and local tissues. Tailing T cell transcriptomes and TCR sequences at the single cell level provides a full landscape of alloreactive T-cell clones development and biofunction in alloresponse. Here, we review the recent advances in TCR sequencing techniques and computational tools, as well as the recent discovery in overall TCR profile and antigen-specific T cells tracking in transplantation. We further discuss the challenges and potential of using TCR sequencing-based assays to profile alloreactive TCR repertoire as the fingerprint for immune monitoring and prediction of rejection and tolerance.

Clonotype pattern in T-cell lymphomas map the cell of origin to immature lymphoid precursors

Mature T-cell lymphomas (TCLs) are rare, clinically heterogeneous hematologic cancers of high medical need. TCLs have inferior prognosis which is attributed to poor understanding of their pathogenesis. Based on phenotypic similarities between normal and neoplastic lymphocytes it has been assumed that TCLs develop in the periphery, directly from various subtypes of normal T-cells. To address the debated question of the cell of origin in TCLs we analyzed to identify the highly variable complementarity determining regions (CDR3) regions of T-cell receptor (TCR) to trace the clonal history of the T-cells. We have collected previously published whole genome -exome, and -transcriptome sequencing data from 574 TCL patients. TCR clonotypes were identified by de novo assembly of CDR3 regions of TCR γ, β and α. We have found that the vast majority of TCLs are clonotypically oligoclonal, although the pattern oligoclonality varied. Anaplastic large cell lymphoma was most diverse comprising multiple clonotypes of TCRγ, β and α whereas adult T-cell lymphoma/leukemia and peripheral T-cell lymphomas often showed monoclonality for TCRγ and β but had diverse TCRα clonotypes. These patterns of rearrangements indicated that TCLs are initiated at the level of the lymphoid precursor. In keeping with this hypothesis, TCR rearrangements in TCLs resembled the pattern seen in the human thymus showing biased usage of V and J segments of high combinatorial probability resulting in recurrent, "public" CDR3 sequences shared across unrelated patients and different clinical TCL entities. Clonotypically diverse initiating cells may seed target tissues being responsible for disease relapses after therapy.

Using the T Cell Receptor as a Biomarker in Type 1 Diabetes

T cell receptors (TCRs) are unique markers that define antigen specificity for a given T cell. With the evolution of sequencing and computational analysis technologies, TCRs are now prime candidates for the development of next-generation non-cell based T cell biomarkers, which provide a surrogate measure to assess the presence of antigen-specific T cells. Type 1 diabetes (T1D), the immune-mediated form of diabetes, is a prototypical organ specific autoimmune disease in which T cells play a pivotal role in targeting pancreatic insulin-producing beta cells. While the disease is now predictable by measuring autoantibodies in the peripheral blood directed to beta cell proteins, there is an urgent need to develop T cell markers that recapitulate T cell activity in the pancreas and can be a measure of disease activity. This review focuses on the potential and challenges of developing TCR biomarkers for T1D. We summarize current knowledge about TCR repertoires and clonotypes specific for T1D and discuss challenges that are unique for autoimmune diabetes. Ultimately, the integration of large TCR datasets produced from individuals with and without T1D along with computational 'big data' analysis will facilitate the development of TCRs as potentially powerful biomarkers in the development of T1D.