VDJdb in 2019: database extension, new analysis infrastructure and a T-cell receptor motif compendium

TCR-ESM: Employing protein language embeddings to predict TCR-peptide-MHC binding

Cognate target identification for T-cell receptors (TCRs) is a significant barrier in T-cell therapy development, which may be overcome by accurately predicting TCR interaction with peptide-bound major histocompatibility complex (pMHC). In this study, we have employed peptide embeddings learned from a large protein language model- Evolutionary Scale Modeling (ESM), to predict TCR-pMHC binding. The TCR-ESM model presented outperforms existing predictors. The complementarity-determining region 3 (CDR3) of the hypervariable TCR is located at the center of the paratope and plays a crucial role in peptide recognition. TCR-ESM trained on paired TCR data with both CDR3α and CDR3β chain information performs significantly better than those trained on data with only CDR3β, suggesting that both TCR chains contribute to specificity, the relative importance however depends on the specific peptide-MHC targeted. The study illuminates the importance of MHC information in TCR-peptide binding which remained inconclusive so far and was thought dependent on the dataset characteristics. TCR-ESM outperforms existing approaches on external datasets, suggesting generalizability. Overall, the potential of deep learning for predicting TCR-pMHC interactions and improving the understanding of factors driving TCR specificity are highlighted. The prediction model is available at http://tcresm.dhanjal-lab.iiitd.edu.in/ as an online tool.

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.

VitTCR: A deep learning method for peptide recognition prediction

This study introduces VitTCR, a predictive model based on the vision transformer (ViT) architecture, aimed at identifying interactions between T cell receptors (TCRs) and peptides, crucial for developing cancer immunotherapies and vaccines. VitTCR converts TCR-peptide interactions into numerical AtchleyMaps using Atchley factors for prediction, achieving AUROC (0.6485) and AUPR (0.6295) values. Benchmark analysis indicates VitTCR's performance is comparable to other models, with further comparative studies suggested to understand its effectiveness in varied contexts. Additionally, integrating a positional bias weight matrix (PBWM), derived from amino acid contact probabilities in structurally resolved pMHC-TCR complexes, slightly improves VitTCR's accuracy. The model's predictions show weak yet statistically significant correlations with immunological factors like T cell clonal expansion and activation percentages, underscoring the biological relevance of VitTCR's predictive capabilities. VitTCR emerges as a valuable computational tool for predicting TCR-peptide interactions, offering insights for immunotherapy and vaccine development.

Tracking in situ checkpoint inhibitor-bound target T cells in patients with checkpoint-induced colitis

The success of checkpoint inhibitors (CPIs) for cancer has been tempered by immune-related adverse effects including colitis. CPI-induced colitis is hallmarked by expansion of resident mucosal IFNγ cytotoxic CD8+ T cells, but how these arise is unclear. Here, we track CPI-bound T cells in intestinal tissue using multimodal single-cell and subcellular spatial transcriptomics (ST). Target occupancy was increased in inflamed tissue, with drug-bound T cells located in distinct microdomains distinguished by specific intercellular signaling and transcriptional gradients. CPI-bound cells were largely CD4+ T cells, including enrichment in CPI-bound peripheral helper, follicular helper, and regulatory T cells. IFNγ CD8+ T cells emerged from both tissue-resident memory (TRM) and peripheral populations, displayed more restricted target occupancy profiles, and co-localized with damaged epithelial microdomains lacking effective regulatory cues. Our multimodal analysis identifies causal pathways and constitutes a resource to inform novel preventive strategies.

Unified cross-modality integration and analysis of T cell receptors and T cell transcriptomes by low-resource-aware representation learning

Single-cell RNA sequencing (scRNA-seq) and T cell receptor sequencing (TCR-seq) are pivotal for investigating T cell heterogeneity. Integrating these modalities, which is expected to uncover profound insights in immunology that might otherwise go unnoticed with a single modality, faces computational challenges due to the low-resource characteristics of the multimodal data. Herein, we present UniTCR, a novel low-resource-aware multimodal representation learning framework designed for the unified cross-modality integration, enabling comprehensive T cell analysis. By designing a dual-modality contrastive learning module and a single-modality preservation module to effectively embed each modality into a common latent space, UniTCR demonstrates versatility in connecting TCR sequences with T cell transcriptomes across various tasks, including single-modality analysis, modality gap analysis, epitope-TCR binding prediction, and TCR profile cross-modality generation, in a low-resource-aware way. Extensive evaluations conducted on multiple scRNA-seq/TCR-seq paired datasets showed the superior performance of UniTCR, exhibiting the ability of exploring the complexity of immune system.

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.

Novel insights into TCR-T cell therapy in solid neoplasms: optimizing adoptive immunotherapy

Adoptive immunotherapy in the T cell landscape exhibits efficacy in cancer treatment. Over the past few decades, genetically modified T cells, particularly chimeric antigen receptor T cells, have enabled remarkable strides in the treatment of hematological malignancies. Besides, extensive exploration of multiple antigens for the treatment of solid tumors has led to clinical interest in the potential of T cells expressing the engineered T cell receptor (TCR). TCR-T cells possess the capacity to recognize intracellular antigen families and maintain the intrinsic properties of TCRs in terms of affinity to target epitopes and signal transduction. Recent research has provided critical insight into their capability and therapeutic targets for multiple refractory solid tumors, but also exposes some challenges for durable efficacy. In this review, we describe the screening and identification of available tumor antigens, and the acquisition and optimization of TCRs for TCR-T cell therapy. Furthermore, we summarize the complete flow from  laboratory to clinical applications of TCR-T cells. Last, we emerge future prospects for improving therapeutic efficacy in cancer world with combination therapies or TCR-T derived products. In conclusion, this review depicts our current understanding of TCR-T cell therapy in solid neoplasms, and provides new perspectives for expanding its clinical applications and improving therapeutic efficacy.

Comparative single-cell RNA sequencing analysis of immune response to inactivated vaccine and natural SARS-CoV-2 infection

Uncovering the immune response to an inactivated SARS-CoV-2 vaccine (In-Vac) and natural infection is crucial for comprehending COVID-19 immunology. Here we conducted an integrated analysis of single-cell RNA sequencing (scRNA-seq) data from serial peripheral blood mononuclear cell (PBMC) samples derived from 12 individuals receiving In-Vac compared with those from COVID-19 patients. Our study reveals that In-Vac induces subtle immunological changes in PBMC, including cell proportions and transcriptomes, compared with profound changes for natural infection. In-Vac modestly upregulates IFN-α but downregulates NF-κB pathways, while natural infection triggers hyperactive IFN-α and NF-κB pathways. Both In-Vac and natural infection alter T/B cell receptor repertoires, but COVID-19 has more significant change in preferential VJ gene, indicating a vigorous immune response. Our study reveals distinct patterns of cellular communications, including a selective activation of IL-15RA/IL-15 receptor pathway after In-Vac boost, suggesting its potential role in enhancing In-Vac-induced immunity. Collectively, our study illuminates multifaceted immune responses to In-Vac and natural infection, providing insights for optimizing SARS-CoV-2 vaccine efficacy.

Transcriptional signature of durable effector T cells elicited by a replication defective HCMV vaccine

Human cytomegalovirus (HCMV) is a leading infectious cause of birth defects and the most common opportunistic infection that causes life-threatening diseases post-transplantation; however, an effective vaccine remains elusive. V160 is a live-attenuated replication defective HCMV vaccine that showed a 42.4% efficacy against primary HCMV infection among seronegative women in a phase 2b clinical trial. Here, we integrated the multicolor flow cytometry, longitudinal T cell receptor (TCR) sequencing, and single-cell RNA/TCR sequencing approaches to characterize the magnitude, phenotype, and functional quality of human T cell responses to V160. We demonstrated that V160 de novo induces IE-1 and pp65 specific durable polyfunctional effector CD8 T cells that are comparable to those induced by natural HCMV infection. We identified a variety of V160-responsive T cell clones which exhibit distinctive "transient" and "durable" expansion kinetics, and revealed a transcriptional signature that marks durable CD8 T cells post-vaccination. Our study enhances the understanding of human T-cell immune responses to V160 vaccination.

Local Enrichment with Convergence of Enriched T-Cell Clones Are Hallmarks of Effective Peptide Vaccination against B16 Melanoma

BACKGROUND

Some peptide anticancer vaccines elicit a strong T-cell memory response but fail to suppress tumor growth. To gain insight into tumor resistance, we compared two peptide vaccines, p20 and p30, against B16 melanoma, with both exhibiting good in vitro T-cell responses but different tumor suppression abilities.

METHODS

We compared activation markers and repertoires of T-lymphocytes from tumor-draining (dLN) and non-draining (ndLN) lymph nodes for the two peptide vaccines.

RESULTS

We showed that the p30 vaccine had better tumor control as opposed to p20. p20 vaccine induced better in vitro T-cell responsiveness but failed to suppress tumor growth. Efficient antitumor vaccination is associated with a higher clonality of cytotoxic T-cells (CTLs) in dLNs compared with ndLNs and the convergence of most of the enriched clones. With the inefficient p20 vaccine, the most expanded and converged were clones of the bystander T-cells without an LN preference.

CONCLUSIONS

Here, we show that the clonality and convergence of the T-cell response are the hallmarks of efficient antitumor vaccination. The high individual and methodological dependencies of these parameters can be avoided by comparing dLNs and ndLNs.

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.

Circulating Tumor Reactive KIR+CD8+ T cells Suppress Anti-Tumor Immunity in Patients with Melanoma

Effective anti-tumor immunity is largely driven by cytotoxic CD8+ T cells that can specifically recognize tumor antigens. However, the factors which ultimately dictate successful tumor rejection remain poorly understood. Here we identify a subpopulation of CD8+ T cells which are tumor antigen-specific in patients with melanoma but resemble KIR+CD8+ T cells with a regulatory function (Tregs). These tumor antigen-specific KIR+CD8+ T cells are detectable in both the tumor and the blood, and higher levels of this population are associated with worse overall survival. Our findings therefore suggest that KIR+CD8+ Tregs are tumor antigen-specific but uniquely suppress anti-tumor immunity in patients with melanoma.

Coordinated Immune Cell Networks in the Bone Marrow Microenvironment Define the Graft versus Leukemia Response with Adoptive Cellular Therapy

Understanding how intra-tumoral immune populations coordinate to generate anti-tumor responses following therapy can guide precise treatment prioritization. We performed systematic dissection of an established adoptive cellular therapy, donor lymphocyte infusion (DLI), by analyzing 348,905 single-cell transcriptomes from 74 longitudinal bone-marrow samples of 25 patients with relapsed myeloid leukemia; a subset was evaluated by protein-based spatial analysis. In acute myelogenous leukemia (AML) responders, diverse immune cell types within the bone-marrow microenvironment (BME) were predicted to interact with a clonally expanded population of ZNF683 + GZMB + CD8+ cytotoxic T lymphocytes (CTLs) which demonstrated in vitro specificity for autologous leukemia. This population, originating predominantly from the DLI product, expanded concurrently with NK and B cells. AML nonresponder BME revealed a paucity of crosstalk and elevated TIGIT expression in CD8+ CTLs. Our study highlights recipient BME differences as a key determinant of effective anti-leukemia response and opens new opportunities to modulate cell-based leukemia-directed therapy.

The signature of a T-cell response to KSHV persists across space and time in individuals with epidemic and endemic KS from Uganda

Inadequate T-cell control of Kaposi sarcoma-associated herpesvirus (KSHV) infection predisposes to development of Kaposi sarcoma (KS), but little is known about the T-cell response to KSHV. Postulating that KS tumors contain abundant KSHV-specific T-cells, we performed transcriptional profiling and T-cell receptor (TCR) repertoire analysis of tumor biopsies from 144 Ugandan adults with KS. We show that CD8+ T-cells and M2-polarized macrophages dominate the tumor micro-environment (TME). The TCR repertoire of KS tumor infiltrating lymphocytes (TIL) is shared across non-contiguous tumors and persists across time. Clusters of T-cells with predicted shared specificity for uncharacterized antigens, potentially encoded by KSHV, comprise ~25% of KS TIL, and are shared across tumors from different time points and individuals. Single-cell RNA-sequencing of blood identifies a non-proliferating effector memory phenotype and captured the TCRs in 14,698 putative KSHV-specific T-cells. These results suggest that a polyspecific KSHV-specific T-cell response inhibited by M2 macrophages exists within the KS TME, and provide a foundation for studies to define its specificity at a large scale.

BATMAN: Improved T cell receptor cross-reactivity prediction benchmarked on a comprehensive mutational scan database

Predicting T cell receptor (TCR) activation is challenging due to the lack of both unbiased benchmarking datasets and computational methods that are sensitive to small mutations to a peptide. To address these challenges, we curated a comprehensive database encompassing complete single amino acid mutational assays of 10,750 TCR-peptide pairs, centered around 14 immunogenic peptides against 66 TCRs. We then present an interpretable Bayesian model, called BATMAN, that can predict the set of peptides that activates a TCR. When validated on our database, BATMAN outperforms existing methods by 20% and reveals important biochemical predictors of TCR-peptide interactions.

Immune responses to SARS-CoV-2 mRNA vaccination in people with idiopathic CD4 lymphopenia

BACKGROUND

The immunogenicity of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) mRNA vaccines is variable in individuals with different inborn errors of immunity or acquired immune deficiencies and is yet unknown in people with idiopathic CD4 lymphopenia (ICL).

OBJECTIVE

We sought to determine the immunogenicity of mRNA vaccines in patients with ICL with a broad range of CD4 T-cell counts.

METHODS

Samples were collected from 25 patients with ICL and 23 age- and sex-matched healthy volunteers (HVs) after their second or third SARS-CoV-2 mRNA vaccine dose. Anti-spike and anti-receptor binding domain antibodies were measured. T-cell receptor sequencing and stimulation assays were performed to quantify SARS-CoV-2-specific T-cell responses.

RESULTS

The median age of ICL participants was 51 years, and their median CD4 count was 150 cells/μL; 11 participants had CD4 counts ≤100 cells/μL. Anti-spike IgG antibody levels were greater in HVs than in patients with ICL after 2 and 3 doses of mRNA vaccine. There was no detectable significant difference, however, in anti-S IgG between HVs and participants with ICL and CD4 counts >100 cells/μL. The depth of spike-specific T-cell responses by T-cell receptor sequencing was lower in individuals with ICL. Activation-induced markers and cytokine production of spike-specific CD4 T cells in participants with ICL did not differ significantly compared with HVs after 2 or 3 vaccine doses.

CONCLUSIONS

Patients with ICL and CD4 counts >100 cells/μL can mount vigorous humoral and cellular immune responses to SARS-CoV-2 vaccination; however, patients with more severe CD4 lymphopenia have blunted vaccine-induced immunity and may require additional vaccine doses and other risk mitigation strategies.

Characterization of T cell receptor repertoire in penile cancer

Tumor-infiltrating lymphocytes (TILs) play a key role in regulating the host immune response and shaping tumor microenvironment. It has been previously shown that T cell infiltration in penile tumors was associated with clinical outcomes. However, few studies have reported the T cell receptor (TCR) repertoire in patients with penile cancer. In the present study, we evaluated the TCR repertoires in tumor and adjacent normal tissues from 22 patients with penile squamous cell carcinoma (PSCC). Analysis of the T cell receptor beta-variable (TRBV) and joining (TRBJ) genes usage and analysis of complementarity determining region 3 (CDR3) length distribution did not show significant differences between tumor and matched normal tissues. Moreover, analysis of the median Jaccard index indicated a limited overlap of TCR repertoire between these groups. Compared with normal tissues, a significantly lower diversity and higher clonality of TCR repertoire was observed in tumor samples, which was associated with clinical characteristics. Further analysis of transcriptional profiles demonstrated that tumor samples with high clonality showed increased expression of genes associated with CD8 + T cells. In addition, we analyzed the TCR repertoire of CD4 + T cells and CD8 + T cells isolated from tumor tissues. We identified that expanded clonotypes were predominantly in the CD8 + T cell compartment, which presented with an exhausted phenotype. Overall, we comprehensively compared TCR repertoire between penile tumor and normal tissues and demonstrated the presence of distinct T cell immune microenvironments in patients with PSCC.

Neoantigen cancer vaccines: a new star on the horizon

Immunotherapy represents a promising strategy for cancer treatment that utilizes immune cells or drugs to activate the patient's own immune system and eliminate cancer cells. One of the most exciting advances within this field is the targeting of neoantigens, which are peptides derived from non-synonymous somatic mutations that are found exclusively within cancer cells and absent in normal cells. Although neoantigen-based therapeutic vaccines have not received approval for standard cancer treatment, early clinical trials have yielded encouraging outcomes as standalone monotherapy or when combined with checkpoint inhibitors. Progress made in high-throughput sequencing and bioinformatics have greatly facilitated the precise and efficient identification of neoantigens. Consequently, personalized neoantigen-based vaccines tailored to each patient have been developed that are capable of eliciting a robust and long-lasting immune response which effectively eliminates tumors and prevents recurrences. This review provides a concise overview consolidating the latest clinical advances in neoantigen-based therapeutic vaccines, and also discusses challenges and future perspectives for this innovative approach, particularly emphasizing the potential of neoantigen-based therapeutic vaccines to enhance clinical efficacy against advanced solid tumors.

Systems immunology spanning tumors, lymph nodes, and periphery

The immune system defines a complex network of tissues and cell types that orchestrate responses across the body in a dynamic manner. The local and systemic interactions between immune and cancer cells contribute to disease progression. Lymphocytes are activated in lymph nodes, traffic through the periphery, and impact cancer progression through their interactions with tumor cells. As a result, therapeutic response and resistance are mediated across tissues, and a comprehensive understanding of lymphocyte dynamics requires a systems-level approach. In this review, we highlight experimental and computational methods that can leverage the study of leukocyte trafficking through an immunomics lens and reveal how adaptive immunity shapes cancer.

GGNpTCR: A Generative Graph Structure Neural Network for Predicting Immunogenic Peptides for T-cell Immune Response

Identifying the interactions between T-cell receptor (TCRs) and human antigens is a crucial step in developing new vaccines, diagnostics, and immunotherapy. Current methods primarily focus on learning binding patterns from known TCR binding repertoires by using sequence information alone without considering the binding specificity of new antigens or exogenous peptides that have not appeared in the training set. Furthermore, the spatial structure of antigens plays a critical role in immune studies and immunotherapy, which should be addressed properly in the identification of interacting TCR-antigen pairs. In this study, we introduced a novel deep learning framework based on generative graph structures, GGNpTCR, for predicting interactions between TCR and peptides from sequence information. Results of real data analysis indicate that our model achieved excellent prediction for new antigens unseen in the training data set, making significant improvements compared to existing methods. We also applied the model to a large COVID-19 data set with no antigens in the training data set, and the improvement was also significant. Furthermore, through incorporation of additional supervised mechanisms, GGNpTCR demonstrated the ability to precisely forecast the locations of peptide-TCR interactions within 3D configurations. This enhancement substantially improved the model's interpretability. In summary, based on the performance on multiple data sets, GGNpTCR has made significant progress in terms of performance, universality, and interpretability.

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.

CAR+ and CAR- T cells share a differentiation trajectory into an NK-like subset after CD19 CAR T cell infusion in patients with B cell malignancies

Chimeric antigen receptor (CAR) T cell therapy is effective in treating B cell malignancies, but factors influencing the persistence of functional CAR+ T cells, such as product composition, patients' lymphodepletion, and immune reconstitution, are not well understood. To shed light on this issue, here we conduct a single-cell multi-omics analysis of transcriptional, clonal, and phenotypic profiles from pre- to 1-month post-infusion of CAR+ and CAR- T cells from patients from a CARTELL study (ACTRN12617001579381) who received a donor-derived 4-1BB CAR product targeting CD19. Following infusion, CAR+ T cells and CAR- T cells shows similar differentiation profiles with clonally expanded populations across heterogeneous phenotypes, demonstrating clonal lineages and phenotypic plasticity. We validate these findings in 31 patients with large B cell lymphoma treated with CD19 CAR T therapy. For these patients, we identify using longitudinal mass-cytometry data an association between NK-like subsets and clinical outcomes at 6 months with both CAR+ and CAR- T cells. These results suggest that non-CAR-derived signals can provide information about patients' immune recovery and be used as correlate of clinically relevant parameters.

Accurate TCR-pMHC interaction prediction using a BERT-based transfer learning method

Accurate prediction of TCR-pMHC binding is important for the development of cancer immunotherapies, especially TCR-based agents. Existing algorithms often experience diminished performance when dealing with unseen epitopes, primarily due to the complexity in TCR-pMHC recognition patterns and the scarcity of available data for training. We have developed a novel deep learning model, 'TCR Antigen Binding Recognition' based on BERT, named as TABR-BERT. Leveraging BERT's potent representation learning capabilities, TABR-BERT effectively captures essential information regarding TCR-pMHC interactions from TCR sequences, antigen epitope sequences and epitope-MHC binding. By transferring this knowledge to predict TCR-pMHC recognition, TABR-BERT demonstrated better results in benchmark tests than existing methods, particularly for unseen epitopes.

Globally shared TCR repertoires within the tumor-infiltrating lymphocytes of patients with metastatic gynecologic cancer

Gynecologic cancer, including ovarian cancer and endometrial cancer, is characterized by morphological and molecular heterogeneity. Germline and somatic testing are available for patients to screen for pathogenic variants in genes such as BRCA1/2. Tissue expression levels of immunogenomic markers such as PD-L1 are also being used in clinical research. The basic therapeutic approach to gynecologic cancer combines surgery with chemotherapy. Immunotherapy, while not yet a mainstream treatment for gynecologic cancers, is advancing, with Dostarlimab recently receiving approval as a treatment for endometrial cancer. The goal remains to harness stimulated immune cells in the bloodstream to eradicate multiple metastases, a feat currently deemed challenging in a typical clinical setting. For the discovery of novel immunotherapy-based tumor targets, tumor-infiltrating lymphocytes (TILs) give a key insight on tumor-related immune activities by providing T cell receptor (TCR) sequences. Understanding the TCR repertoires of TILs in metastatic tissues and the circulation is important from an immunotherapy standpoint, as a subset of T cells in the blood have the potential to help kill tumor cells. To explore the relationship between distant tissue biopsy regions and blood circulation, we investigated the TCR beta chain (TCRβ) in bulk tumor and matched blood samples from 39 patients with gynecologic cancer. We found that the TCR clones of TILs at different tumor sites were globally shared within patients and had high overlap with the TCR clones in peripheral blood.

Dominant CD4+ T cell receptors remain stable throughout antiretroviral therapy-mediated immune restoration in people with HIV

In people with HIV (PWH), the post-antiretroviral therapy (ART) window is critical for immune restoration and HIV reservoir stabilization. We employ deep immune profiling and T cell receptor (TCR) sequencing and examine proliferation to assess how ART impacts T cell homeostasis. In PWH on long-term ART, lymphocyte frequencies and phenotypes are mostly stable. By contrast, broad phenotypic changes in natural killer (NK) cells, γδ T cells, B cells, and CD4+ and CD8+ T cells are observed in the post-ART window. Whereas CD8+ T cells mostly restore, memory CD4+ T subsets and cytolytic NK cells show incomplete restoration 1.4 years post ART. Surprisingly, the hierarchies and frequencies of dominant CD4 TCR clonotypes (0.1%-11% of all CD4+ T cells) remain stable post ART, suggesting that clonal homeostasis can be independent of homeostatic processes regulating CD4+ T cell absolute number, phenotypes, and function. The slow restoration of host immunity post ART also has implications for the design of ART interruption studies.

Transcriptome-based identification of tumor-reactive and bystander CD8+ T cell receptor clonotypes in human pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) is generally refractory to immune checkpoint blockade, although patients with genetically unstable tumors can show modest therapeutic benefit. We previously demonstrated the presence of tumor-reactive CD8+ T cells in PDAC samples. Here, we charted the tumor-infiltrating T cell repertoire in PDAC by combining single-cell transcriptomics with functional testing of T cell receptors (TCRs) for reactivity against autologous tumor cells. On the basis of a comprehensive dataset including 93 tumor-reactive and 65 bystander TCR clonotypes, we delineated a gene signature that effectively distinguishes between these T cell subsets in PDAC, as well as in other tumor indications. This revealed a high frequency of tumor-reactive TCR clonotypes in three genetically unstable samples. In contrast, the T cell repertoire in six genetically stable PDAC tumors was largely dominated by bystander T cells. Nevertheless, multiple tumor-reactive TCRs were successfully identified in each of these samples, thereby providing a perspective for personalized immunotherapy in this treatment-resistant indication.

Specific T-cell receptor beta-rearrangements of gluten-triggered CD8+ T-cells are enriched in celiac disease patients' duodenal mucosa

Celiac disease (CeD) is an autoimmune disorder affecting the small intestine with gluten as disease trigger. Infections including Influenza A, increase the CeD risk. While gluten-specific CD4+ T-cells, recognizing HLA-DQ2/DQ8 presented gluten-peptides, initiate and sustain the celiac immune response, CD8+ α/β intraepithelial T-cells elicit mucosal damage. Here, we subjected TCRs from a cohort of 56 CeD patients and 22 controls to an analysis employing 749 published CeD-related TCRβ-rearrangements derived from gluten-specific CD4+ T-cells and gluten-triggered peripheral blood CD8+ T-cells. We show, that in addition to TCRs from gluten-specific CD4+ T-cells, TCRs of gluten-triggered CD8+ T-cells are significantly enriched in CeD duodenal tissue samples. TCRβ-rearrangements of gluten-triggered CD8+ T-cells were even more expanded in patients than TCRs from gluten-specific CD4+ T-cells (p < 0.0002) and highest in refractory CeD. Sequence alignments with TCR-antigen databases suggest that a subgroup of these most likely indirectly gluten-triggered TCRs recognize microbial, viral, and autoantigens.

An integrated approach to the characterization of immune repertoires using AIMS: An Automated Immune Molecule Separator

The adaptive immune system employs an array of receptors designed to respond with high specificity to pathogens or molecular aberrations faced by the host organism. Binding of these receptors to molecular fragments-collectively referred to as antigens-initiates immune responses. These antigenic targets are recognized in their native state on the surfaces of pathogens by antibodies, whereas T cell receptors (TCR) recognize processed antigens as short peptides, presented on major histocompatibility complex (MHC) molecules. Recent research has led to a wealth of immune repertoire data that are key to interrogating the nature of these molecular interactions. However, existing tools for the analysis of these large datasets typically focus on molecular sets of a single type, forcing researchers to separately analyze strongly coupled sequences of interacting molecules. Here, we introduce a software package for the integrated analysis of immune repertoire data, capable of identifying distinct biophysical differences in isolated TCR, MHC, peptide, antibody, and antigen sequence data. This integrated analytical approach allows for direct comparisons across immune repertoire subsets and provides a starting point for the identification of key interaction hotspots in complementary receptor-antigen pairs. The software (AIMS-Automated Immune Molecule Separator) is freely available as an open access package in GUI or command-line form.

Blood-based T cell receptor anti-viral CDR3s are associated with worse overall survival for neuroblastoma

With the advent of large collections of adaptive immune receptor recombination reads representing cancer, there is the opportunity to further investigate the adaptive immune response to viruses in the cancer setting. This is a particularly important goal due to longstanding but still not well-resolved questions about viral etiologies in cancer and viral infections as comorbidities. In this report, we assessed the T cell receptor complementarity determining region-3 (CDR3) amino acid (AA) sequences, for blood-sourced TCRs from neuroblastoma (NBL) cases, for exact AA sequence matches to previously identified anti-viral TCR CDR3 AA sequences. Results indicated the presence of anti-viral TCR CDR3 AA sequences in the NBL blood samples highly significantly correlated with worse overall survival. Furthermore, the TCR CDR3 AA sequences demonstrating chemical complementarity to many cytomegalovirus antigens represented cases with a worse outcome, including cases where such CDR3s were obtained from tumor samples. Overall, these results indicate a significant need for, and provide a novel strategy for assessing viral infection complications in NBL patients.

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.

Training of epitope-TCR prediction models with healthy donor-derived cancer-specific T cells

Discovery of epitope-specific T-cell receptors (TCRs) for cancer therapies is a time consuming and expensive procedure that usually requires a large amount of patient cells. To maximize information from and minimize the need of precious samples in cancer research, prediction models have been developed to identify in silico epitope-specific TCRs. In this chapter, we provide a step-by-step protocol to train a prediction model using the user-friendly TCRex webtool for the nearly universal tumor-associated antigen Wilms' tumor 1 (WT1)-specific TCR repertoire. WT1 is a self-antigen overexpressed in numerous solid and hematological malignancies with a high clinical relevance. Training of computational models starts from a list of known epitope-specific TCRs which is often not available for new cancer epitopes. Therefore, we describe a workflow to assemble a training data set consisting of TCR sequences obtained from WT137-45-reactive CD8 T cell clones expanded and sorted from healthy donor peripheral blood mononuclear cells.

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.

Diversification and expansion of the EBV-reactive cytotoxic T lymphocyte repertoire following autologous haematopoietic stem cell transplant for multiple sclerosis

Both genetic susceptibility and environmental exposures are thought to be involved in multiple sclerosis (MS) pathogenesis. Of all viruses potentially relevant to MS aetiology, Epstein-Barr virus (EBV) is the best-studied. EBV is a B cell lymphotropic virus which is able to evade the immune system by establishing latent infection in memory B cells, and EBV reactivation is restricted by CD8 cytotoxic T cell (CTL) responses in immune competent individuals. Autologous haematopoietic stem cell transplantation (AHSCT) is considered to be the most effective therapy in the treatment of relapsing MS even though chemotherapy-induced lymphopenia can associate with the re-emergence of latent viruses. Despite the increasing interest in EBV and MS pathogenesis the relationship between AHSCT, EBV and viral immunity in people with MS has not been investigated to date. This study analysed immune responses to EBV in a well characterised cohort of 13 individuals with MS by utilising pre-AHSCT, and 6-, 12- and 24-month post AHSCT bio-banked peripheral blood mononuclear cells and plasma samples. It is demonstrated that the infused stem cell product contains latently EBV-infected memory B cells, and that EBV viremia occurs in the immune-compromised recipient post-transplant. High throughput TCR analysis detected expansion and diversification of the CD8 CTL responses reactive with EBV lytic and latent antigens from 6 to 24 months following AHSCT. Increased levels of latent EBV infection found within the B cell pool following treatment, as measured by EBV genomic detection, did not associate with disease relapse. This is the first study of EBV immunity following application of AHSCT in the treatment of MS and not only raises important questions about the role of EBV infection in MS pathogenesis, but is of clinical importance given the expanding clinical trials of adoptive EBV-specific CTLs in MS.

Resolving the spatial architecture of myeloma and its microenvironment at the single-cell level

In multiple myeloma spatial differences in the subclonal architecture, molecular signatures and composition of the microenvironment remain poorly characterized. To address this shortcoming, we perform multi-region sequencing on paired random bone marrow and focal lesion samples from 17 newly diagnosed patients. Using single-cell RNA- and ATAC-seq we find a median of 6 tumor subclones per patient and unique subclones in focal lesions. Genetically identical subclones display different levels of spatial transcriptional plasticity, including nearly identical profiles and pronounced heterogeneity at different sites, which can include differential expression of immunotherapy targets, such as CD20 and CD38. Macrophages are significantly depleted in the microenvironment of focal lesions. We observe proportional changes in the T-cell repertoire but no site-specific expansion of T-cell clones in intramedullary lesions. In conclusion, our results demonstrate the relevance of considering spatial heterogeneity in multiple myeloma with potential implications for models of cell-cell interactions and disease progression.

Large-scale template-based structural modeling of T-cell receptors with known antigen specificity reveals complementarity features

Introduction

T-cell receptor (TCR) recognition of foreign peptides presented by the major histocompatibility complex (MHC) initiates the adaptive immune response against pathogens. While a large number of TCR sequences specific to different antigenic peptides are known to date, the structural data describing the conformation and contacting residues for TCR-peptide-MHC complexes is relatively limited. In the present study we aim to extend and analyze the set of available structures by performing highly accurate template-based modeling of these complexes using TCR sequences with known specificity.

Methods

Identification of CDR3 sequences and their further clustering, based on available spatial structures, V- and J-genes of corresponding T-cell receptors, and epitopes, was performed using the VDJdb database. Modeling of the selected CDR3 loops was conducted using a stepwise introduction of single amino acid substitutions to the template PDB structures, followed by optimization of the TCR-peptide-MHC contacting interface using the Rosetta package applications. Statistical analysis and recursive feature elimination procedures were carried out on computed energy values and properties of contacting amino acid residues between CDR3 loops and peptides, using R.

Results

Using the set of 29 complex templates (including a template with SARS-CoV-2 antigen) and 732 specificity records, we built a database of 1585 model structures carrying substitutions in either TCRα or TCRβ chains with some models representing the result of different mutation pathways for the same final structure. This database allowed us to analyze features of amino acid contacts in TCR - peptide interfaces that govern antigen recognition preferences and interpret these interactions in terms of physicochemical properties of interacting residues.

Conclusion

Our results provide a methodology for creating high-quality TCR-peptide-MHC models for antigens of interest that can be utilized to predict TCR specificity.

MIX-TPI: a flexible prediction framework for TCR-pMHC interactions based on multimodal representations

MOTIVATION

The interactions between T-cell receptors (TCR) and peptide-major histocompatibility complex (pMHC) are essential for the adaptive immune system. However, identifying these interactions can be challenging due to the limited availability of experimental data, sequence data heterogeneity, and high experimental validation costs.

RESULTS

To address this issue, we develop a novel computational framework, named MIX-TPI, to predict TCR-pMHC interactions using amino acid sequences and physicochemical properties. Based on convolutional neural networks, MIX-TPI incorporates sequence-based and physicochemical-based extractors to refine the representations of TCR-pMHC interactions. Each modality is projected into modality-invariant and modality-specific representations to capture the uniformity and diversities between different features. A self-attention fusion layer is then adopted to form the classification module. Experimental results demonstrate the effectiveness of MIX-TPI in comparison with other state-of-the-art methods. MIX-TPI also shows good generalization capability on mutual exclusive evaluation datasets and a paired TCR dataset.

AVAILABILITY AND IMPLEMENTATION

The source code of MIX-TPI and the test data are available at: https://github.com/Wolverinerine/MIX-TPI.

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.

Context-dependent T-cell Receptor Gene Repertoire Profiles in Proliferations of T Large Granular Lymphocytes

T cell large granular lymphocyte (T-LGL) lymphoproliferations constitute a disease spectrum ranging from poly/oligo to monoclonal. Boundaries within this spectrum of proliferations are not well established. T-LGL lymphoproliferations co-occur with a wide variety of other diseases ranging from autoimmune disorders, solid tumors, hematological malignancies, post solid organ, and hematopoietic stem cell transplantation, and can therefore arise as a consequence of a wide variety of antigenic triggers. Persistence of a dominant malignant T-LGL clone is established through continuous STAT3 activation. Using next-generation sequencing, we profiled a cohort of 27 well-established patients with T-LGL lymphoproliferations, aiming to identify the subclonal architecture of the T-cell receptor beta (TRB) chain gene repertoire. Moreover, we searched for associations between TRB gene repertoire patterns and clinical manifestations, with the ultimate objective of discriminating between T-LGL lymphoproliferations developing in different clinical contexts and/or displaying distinct clinical presentation. Altogether, our data demonstrates that the TRB gene repertoire of patients with T-LGL lymphoproliferations is context-dependent, displaying distinct clonal architectures in different settings. Our results also highlight that there are monoclonal T-LGL cells with or without STAT3 mutations that cause symptoms such as neutropenia on one end of a spectrum and reactive oligoclonal T-LGL lymphoproliferations on the other. Longitudinal analysis revealed temporal clonal dynamics and showed that T-LGL cells might arise as an epiphenomenon when co-occurring with other malignancies, possibly reactive toward tumor antigens.

Immunologic Characterization and T cell Receptor Repertoires of Expanded Tumor-infiltrating Lymphocytes in Patients with Renal Cell Carcinoma

The successful use of expanded tumor-infiltrating lymphocytes (TIL) in adoptive TIL therapies has been reported, but the effects of the TIL expansion, immunophenotype, function, and T cell receptor (TCR) repertoire of the infused products relative to the tumor microenvironment (TME) are not well understood. In this study, we analyzed the tumor samples (n = 58) from treatment-naïve patients with renal cell carcinoma (RCC), "pre-rapidly expanded" TILs (pre-REP TIL, n = 15) and "rapidly expanded" TILs (REP TIL, n = 25) according to a clinical-grade TIL production protocol, with single-cell RNA (scRNA)+TCRαβ-seq (TCRαβ sequencing), TCRβ-sequencing (TCRβ-seq), and flow cytometry. REP TILs encompassed a greater abundance of CD4⁺ than CD8⁺ T cells, with increased LAG-3 and low PD-1 expressions in both CD4⁺ and CD8⁺ T cell compartments compared with the pre-REP TIL and tumor T cells. The REP protocol preferentially expanded small clones of the CD4⁺ phenotype (CD4, IL7R, KLRB1) in the TME, indicating that the largest exhausted T cell clones in the tumor do not expand during the expansion protocol. In addition, by generating a catalog of RCC-associated TCR motifs from >1,000 scRNA+TCRαβ-seq and TCRβ-seq RCC, healthy and other cancer sample cohorts, we quantified the RCC-associated TCRs from the expansion protocol. Unlike the low-remaining amount of anti-viral TCRs throughout the expansion, the quantity of the RCC-associated TCRs was high in the tumors and pre-REP TILs but decreased in the REP TILs. Our results provide an in-depth understanding of the origin, phenotype, and TCR specificity of RCC TIL products, paving the way for a more rationalized production of TILs.

Significance

TILs are a heterogenous group of immune cells that recognize and attack the tumor, thus are utilized in various clinical trials. In our study, we explored the TILs in patients with kidney cancer by expanding the TILs using a clinical-grade protocol, as well as observed their characteristics and ability to recognize the tumor using in-depth experimental and computational tools.

Exploring the Potential of GANs in Biological Sequence Analysis

Biological sequence analysis is an essential step toward building a deeper understanding of the underlying functions, structures, and behaviors of the sequences. It can help in identifying the characteristics of the associated organisms, such as viruses, etc., and building prevention mechanisms to eradicate their spread and impact, as viruses are known to cause epidemics that can become global pandemics. New tools for biological sequence analysis are provided by machine learning (ML) technologies to effectively analyze the functions and structures of the sequences. However, these ML-based methods undergo challenges with data imbalance, generally associated with biological sequence datasets, which hinders their performance. Although various strategies are present to address this issue, such as the SMOTE algorithm, which creates synthetic data, however, they focus on local information rather than the overall class distribution. In this work, we explore a novel approach to handle the data imbalance issue based on generative adversarial networks (GANs), which use the overall data distribution. GANs are utilized to generate synthetic data that closely resembles real data, thus, these generated data can be employed to enhance the ML models' performance by eradicating the class imbalance problem for biological sequence analysis. We perform four distinct classification tasks by using four different sequence datasets (Influenza A Virus, PALMdb, VDjDB, Host) and our results illustrate that GANs can improve the overall classification performance.

Analysis of thymic generation of shared T-cell receptor α repertoire associated with recognition of tumor antigens shows no preference for neoantigens over wild-type antigens

BACKGROUND

The number of mutations in cancer cells is an important predictor of a positive response to cancer immunotherapy. It has been suggested that the neoantigens produced by these mutations are more immunogenic than nonmutated tumor antigens, which are likely to be protected by immunological tolerance. However, the mechanisms of tolerance as regards tumor antigens are incompletely understood.

METHODS

Here, we have analyzed the impact of thymic negative selection on shared T-cell receptor (TCR) repertoire associated with the recognition of either mutated or nonmutated tumor antigens by comparing previously known TCR-antigen-pairs to TCR repertoires of 21 immunologically healthy individuals.

RESULTS

Our results show that TCRα chains associated with either type of tumor antigens are readily generated in the thymus, at a frequency similar to TCRα chains associated with nonself. In the peripheral repertoire, the relative clone size of nonself-associated chains is higher than that of the tumor antigens, but importantly, there is no difference between TCRα chains associated with mutated or nonmutated tumor antigens.

CONCLUSION

This suggests that the tolerance mechanisms protecting nonmutated tumor antigens are non-deletional and therefore potentially reversible. As unmutated antigens are, unlike mutations, shared by a large number of patients, they may offer advantages in designing immunological approaches to cancer treatment.

Quantitative annotations of T-Cell repertoire specificity

The specificity of a T-cell receptor (TCR) repertoire determines personalized immune capacity. Existing methods have modeled the qualitative aspects of TCR specificity, while the quantitative aspects remained unaddressed. We developed a package, TCRanno, to quantify the specificity of TCR repertoires. We created deep-learning-based, epitope-aware vector embeddings to infer individual TCR specificity. Then we aggregated clonotype frequencies of TCRs to obtain a quantitative profile of repertoire specificity at epitope, antigen and organism levels. Applying TCRanno to 4195 TCR repertoires revealed quantitative changes in repertoire specificity upon infections, autoimmunity and cancers. Specifically, TCRanno found cytomegalovirus-specific TCRs in seronegative healthy individuals, supporting the possibility of abortive infections. TCRanno discovered age-accumulated fraction of severe acute respiratory syndrome coronavirus 2 specific TCRs in pre-pandemic samples, which may explain the aggressive symptoms and age-related severity of coronavirus disease 2019. TCRanno also identified the encounter of Hepatitis B antigens as a potential trigger of systemic lupus erythematosus. TCRanno annotations showed capability in distinguishing TCR repertoires of healthy and cancers including melanoma, lung and breast cancers. TCRanno also demonstrated usefulness to single-cell TCRseq+gene expression data analyses by isolating T-cells with the specificity of interest.

iTCep: a deep learning framework for identification of T cell epitopes by harnessing fusion features

Neoantigens recognized by cytotoxic T cells are effective targets for tumor-specific immune responses for personalized cancer immunotherapy. Quite a few neoantigen identification pipelines and computational strategies have been developed to improve the accuracy of the peptide selection process. However, these methods mainly consider the neoantigen end and ignore the interaction between peptide-TCR and the preference of each residue in TCRs, resulting in the filtered peptides often fail to truly elicit an immune response. Here, we propose a novel encoding approach for peptide-TCR representation. Subsequently, a deep learning framework, namely iTCep, was developed to predict the interactions between peptides and TCRs using fusion features derived from a feature-level fusion strategy. The iTCep achieved high predictive performance with AUC up to 0.96 on the testing dataset and above 0.86 on independent datasets, presenting better prediction performance compared with other predictors. Our results provided strong evidence that model iTCep can be a reliable and robust method for predicting TCR binding specificities of given antigen peptides. One can access the iTCep through a user-friendly web server at http://biostatistics.online/iTCep/, which supports prediction modes of peptide-TCR pairs and peptide-only. A stand-alone software program for T cell epitope prediction is also available for convenient installing at https://github.com/kbvstmd/iTCep/.

TCRmodel2: high-resolution modeling of T cell receptor recognition using deep learning

The cellular immune system, which is a critical component of human immunity, uses T cell receptors (TCRs) to recognize antigenic proteins in the form of peptides presented by major histocompatibility complex (MHC) proteins. Accurate definition of the structural basis of TCRs and their engagement of peptide-MHCs can provide major insights into normal and aberrant immunity, and can help guide the design of vaccines and immunotherapeutics. Given the limited amount of experimentally determined TCR-peptide-MHC structures and the vast amount of TCRs within each individual as well as antigenic targets, accurate computational modeling approaches are needed. Here, we report a major update to our web server, TCRmodel, which was originally developed to model unbound TCRs from sequence, to now model TCR-peptide-MHC complexes from sequence, utilizing several adaptations of AlphaFold. This method, named TCRmodel2, allows users to submit sequences through an easy-to-use interface and shows similar or greater accuracy than AlphaFold and other methods to model TCR-peptide-MHC complexes based on benchmarking. It can generate models of complexes in 15 minutes, and output models are provided with confidence scores and an integrated molecular viewer. TCRmodel2 is available at https://tcrmodel.ibbr.umd.edu.

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.

TPBTE: A model based on convolutional Transformer for predicting the binding of TCR to epitope

T cell receptors (TCRs) selectively bind to antigens to fight pathogens with specific immunity. Current tools focus on the nature of amino acids within sequences and take less into account the nature of amino acids far apart and the relationship between sequences, leading to significant differences in the results from different datasets. We propose TPBTE, a model based on convolutional Transformer for Predicting the Binding of TCR to Epitope. It takes epitope sequences and the complementary decision region 3 (CDR3) sequences of TCRβ chain as inputs. And it uses a convolutional attention mechanism to learn amino acid representations between different positions of the sequences based on learning local features of the sequences. At the same time, it uses cross attention to learn the interaction information between TCR sequences and epitope sequences. A comprehensive evaluation of the TCR-epitope data shows that the average area under the curve of TPBTE outperforms the baseline model, and demonstrate an intentional performance. In addition, TPBTE can give the probability of binding TCR to epitopes, which can be used as the first step of epitope screening, narrowing the scope of epitope search and reducing the time of epitope search.

Single-cell analyses and host genetics highlight the role of innate immune cells in COVID-19 severity

Mechanisms underpinning the dysfunctional immune response in severe acute respiratory syndrome coronavirus 2 infection are elusive. We analyzed single-cell transcriptomes and T and B cell receptors (BCR) of >895,000 peripheral blood mononuclear cells from 73 coronavirus disease 2019 (COVID-19) patients and 75 healthy controls of Japanese ancestry with host genetic data. COVID-19 patients showed a low fraction of nonclassical monocytes (ncMono). We report downregulated cell transitions from classical monocytes to ncMono in COVID-19 with reduced CXCL10 expression in ncMono in severe disease. Cell-cell communication analysis inferred decreased cellular interactions involving ncMono in severe COVID-19. Clonal expansions of BCR were evident in the plasmablasts of patients. Putative disease genes identified by COVID-19 genome-wide association study showed cell type-specific expressions in monocytes and dendritic cells. A COVID-19-associated risk variant at the IFNAR2 locus (rs13050728) had context-specific and monocyte-specific expression quantitative trait loci effects. Our study highlights biological and host genetic involvement of innate immune cells in COVID-19 severity.

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.

Single-cell immune profiling reveals immune responses in oral lichen planus

Introduction

Oral lichen planus (OLP) is a common chronic inflammatory disorder of the oral mucosa with an unclear etiology. Several types of immune cells are involved in the pathogenesis of OLP.

Methods

We used single-cell RNA sequencing and immune repertoire sequencing to characterize the mucosal immune microenvironment of OLP. The presence of tissue-resident memory CD8+ T cells are validated by multiplex immunofluorescence.

Results

We generated a transcriptome atlas from four OLP biopsy samples and their paired peripheral blood mononuclear cells (PBMCs), and compared them with two healthy tissues and three healthy PBMCs samples. Our analysis revealed activated tissue-resident memory CD8+ T cells in OLP tissues. T cell receptor repertoires displayed apperant clonal expansion and preferrential gene pairing in OLP patients. Additionally, obvious BCR clonal expansion was observed in OLP lesions. Plasmacytoid dendritic cells, a subtype that can promote dendritic cell maturation and enhance lymphocyte cytotoxicity, were identified in OLP. Conventional dendritic cells and macrophages are also found to exhibit pro-inflammatory activity in OLP. Cell-cell communication analysis reveals that fibroblasts might promote the recruitment and extravasation of immune cells into connective tissue.

Discussion

Our study provides insights into the immune ecosystem of OLP, serving as a valuable resource for precision diagnosis and therapy of OLP.