Detecting T cell receptors involved in immune responses from single repertoire snapshots

Statistical analysis of repertoire data demonstrates the influence of microhomology in V(D)J recombination

V(D)J recombination generates the diverse B and T cell receptors essential for recognizing a wide array of antigens. This diversity arises from the combinatorial assembly of V(D)J genes and the junctional deletion and insertion of nucleotides. While previous in vitro studies have shown that microhomology-short stretches of sequence homology between gene ends-can bias the recombination process, the extent of microhomology's impact in vivo, particularly in humans, remains unknown. In this paper, we assess how germline-encoded microhomology influences trimming and ligation during V(D)J recombination using statistical inference on previously published high-throughput TCRα repertoire sequencing data. We find that microhomology increases both trimming and ligation probabilities, making it an important predictor of recombination outcomes. These effects are consistent across other receptor loci and sequence types. Further, we demonstrate that accounting for germline microhomology effects significantly alters sequence annotation probabilities and rankings, highlighting its practical importance for accurately inferring the V(D)J recombination events that generated an observed sequence. Together, these results enhance our understanding of how germline-encoded microhomologous nucleotides shape the human V(D)J recombination process.

Simulation of adaptive immune receptors and repertoires with complex immune information to guide the development and benchmarking of AIRR machine learning

Machine learning (ML) has shown great potential in the adaptive immune receptor repertoire (AIRR) field. However, there is a lack of large-scale ground-truth experimental AIRR data suitable for AIRR-ML-based disease diagnostics and therapeutics discovery. Simulated ground-truth AIRR data are required to complement the development and benchmarking of robust and interpretable AIRR-ML methods where experimental data is currently inaccessible or insufficient. The challenge for simulated data to be useful is incorporating key features observed in experimental repertoires. These features, such as antigen or disease-associated immune information, cause AIRR-ML problems to be challenging. Here, we introduce LIgO, a software suite, which simulates AIRR data for the development and benchmarking of AIRR-ML methods. LIgO incorporates different types of immune information both on the receptor and the repertoire level and preserves native-like generation probability distribution. Additionally, LIgO assists users in determining the computational feasibility of their simulations. We show two examples where LIgO supports the development and validation of AIRR-ML methods: (i) how individuals carrying out-of-distribution immune information impacts receptor-level prediction performance and (ii) how immune information co-occurring in the same AIRs impacts the performance of conventional receptor-level encoding and repertoire-level classification approaches. LIgO guides the advancement and assessment of interpretable AIRR-ML methods.

Learning predictive signatures of HLA type from T-cell repertoires

T cells recognize a wide range of pathogens using surface receptors that interact directly with peptides presented on major histocompatibility complexes (MHC) encoded by the HLA loci in humans. Understanding the association between T cell receptors (TCR) and HLA alleles is an important step towards predicting TCR-antigen specificity from sequences. Here we analyze the TCR alpha and beta repertoires of large cohorts of HLA-typed donors to systematically infer such associations, by looking for overrepresentation of TCRs in individuals with a common allele.TCRs, associated with a specific HLA allele, exhibit sequence similarities that suggest prior antigen exposure. Immune repertoire sequencing has produced large numbers of datasets, however the HLA type of the corresponding donors is rarely available. Using our TCR-HLA associations, we trained a computational model to predict the HLA type of individuals from their TCR repertoire alone. We propose an iterative procedure to refine this model by using data from large cohorts of untyped individuals, by recursively typing them using the model itself. The resulting model shows good predictive performance, even for relatively rare HLA alleles.

Prediction of tumor-reactive T cell receptors from scRNA-seq data for personalized T cell therapy

The identification of patient-derived, tumor-reactive T cell receptors (TCRs) as a basis for personalized transgenic T cell therapies remains a time- and cost-intensive endeavor. Current approaches to identify tumor-reactive TCRs analyze tumor mutations to predict T cell activating (neo)antigens and use these to either enrich tumor infiltrating lymphocyte (TIL) cultures or validate individual TCRs for transgenic autologous therapies. Here we combined high-throughput TCR cloning and reactivity validation to train predicTCR, a machine learning classifier that identifies individual tumor-reactive TILs in an antigen-agnostic manner based on single-TIL RNA sequencing. PredicTCR identifies tumor-reactive TCRs in TILs from diverse cancers better than previous gene set enrichment-based approaches, increasing specificity and sensitivity (geometric mean) from 0.38 to 0.74. By predicting tumor-reactive TCRs in a matter of days, TCR clonotypes can be prioritized to accelerate the manufacture of personalized T cell therapies.

Analysis of Wilms' tumor protein 1 specific TCR repertoire in AML patients uncovers higher diversity in patients in remission than in relapsed

The Wilms' tumor protein 1 (WT1) is a well-known and prioritized tumor-associated antigen expressed in numerous solid and blood tumors. Its abundance and immunogenicity have led to the development of different WT1-specific immune therapies. The driving player in these therapies, the WT1-specific T-cell receptor (TCR) repertoire, has received much less attention. Importantly, T cells with high affinity against the WT1 self-antigen are normally eliminated after negative selection in the thymus and are thus rare in peripheral blood. Here, we developed computational models for the robust and fast identification of WT1-specific TCRs from TCR repertoire data. To this end, WT137-45 (WT1-37) and WT1126-134 (WT1-126)-specific T cells were isolated from WT1 peptide-stimulated blood of healthy individuals. The TCR repertoire from these WT1-specific T cells was sequenced and used to train a pattern recognition model for the identification of WT1-specific TCR patterns for the WT1-37 or WT1-126 epitopes. The resulting computational models were applied on an independent published dataset from acute myeloid leukemia (AML) patients, treated with hematopoietic stem cell transplantation, to track WT1-specific TCRs in silico. Several WT1-specific TCRs were found in AML patients. Subsequent clustering analysis of all repertoires indicated the presence of more diverse TCR patterns within the WT1-specific TCR repertoires of AML patients in complete remission in contrast to relapsing patients. We demonstrate the possibility of tracking WT1-37 and WT1-126-specific TCRs directly from TCR repertoire data using computational methods, eliminating the need for additional blood samples and experiments for the two studied WT1 epitopes.

Reading the repertoire: Progress in adaptive immune receptor analysis using machine learning

The adaptive immune system holds invaluable information on past and present immune responses in the form of B and T cell receptor sequences, but we are limited in our ability to decode this information. Machine learning approaches are under active investigation for a range of tasks relevant to understanding and manipulating the adaptive immune receptor repertoire, including matching receptors to the antigens they bind, generating antibodies or T cell receptors for use as therapeutics, and diagnosing disease based on patient repertoires. Progress on these tasks has the potential to substantially improve the development of vaccines, therapeutics, and diagnostics, as well as advance our understanding of fundamental immunological principles. We outline key challenges for the field, highlighting the need for software benchmarking, targeted large-scale data generation, and coordinated research efforts.

Investigating TCR-pMHC interactions for TCRs without identified epitopes by constructing a computational pipeline

The molecular mechanisms underlying epitope recognition by T cell receptors (TCRs) are critical for activating T cell immune responses and rationally designing TCR-based therapeutics. Single-cell sequencing techniques vastly boost the accumulation of TCR sequences, while the limitation of available TCR-pMHC structures hampers further investigations. In this study, we proposed a computational pipeline that incorporates structural information and single-cell sequencing data to investigate the epitope-recognition mechanisms for TCRs without identified epitopes. By antigen specificity clustering, we mapped the epitope sequences between epitope-known and epitope-unknown TCRs from COVID-19 patients. One reported SARS-CoV-2 epitope, NQKLIANQF (S919-927), was identified for a TCR expressed by 614 T cells (TCR-614). Epitope screening also identified a potential cross-reactive epitope, KLKTLVATA (NSP31790-1798), for a TCR expressed by 204 T cells (TCR-204). By molecular dynamics (MD) simulations, we revealed the detailed epitope-recognition mechanisms for both TCRs. The structural motifs responsible for epitope recognition revealed by the MD simulations are consistent with the sequential features recognized by the sequence-based clustering method. We hope that this strategy could facilitate the discovery and optimization of TCR-based therapeutics.

Bulk T cell repertoire sequencing (TCR-Seq) is a powerful technology for understanding inflammation-mediated diseases

Multiple alterations in the T cell repertoire were identified in many chronic inflammatory diseases such as inflammatory bowel disease, multiple sclerosis, and rheumatoid arthritis, suggesting that T cells might, directly or indirectly, be implicated in these pathologies. This has sparked a deep interest in characterizing disease-associated T cell clonotypes as well as in identifying and quantifying their contribution to the pathophysiology of different autoimmune and inflammation-mediated diseases. Bulk T cell repertoire sequencing (TCR-Seq) has emerged as a powerful method to profile the T cell repertoire of a sample in a high throughput fashion. Given the increasing utilization of TCR-Seq, we aimed here to provide a comprehensive, up-to-date review of the method, its extensions, and its ability to investigate chronic and autoimmune diseases. Specifically, we started by introducing the immunological basis of TCR repertoire generation and features, followed by discussing different experimental approach to perform TCR-Seq, then we describe different methods and frameworks for analyzing the generated datasets. Subsequently, different experimental techniques for investigating the antigenicity of T cell clonotypes are described. Lastly, we discuss recent studies that utilized TCR-Seq to understand different inflammation-mediated diseases, discuss fallbacks of the technology and potential future directions to overcome current limitations.

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.

Predicting Antigen-Specificities of Orphan T Cell Receptors from Cancer Patients with TCRpcDist

Approaches to analyze and cluster T-cell receptor (TCR) repertoires to reflect antigen specificity are critical for the diagnosis and prognosis of immune-related diseases and the development of personalized therapies. Sequence-based approaches showed success but remain restrictive, especially when the amount of experimental data used for the training is scarce. Structure-based approaches which represent powerful alternatives, notably to optimize TCRs affinity toward specific epitopes, show limitations for large-scale predictions. To handle these challenges, TCRpcDist is presented, a 3D-based approach that calculates similarities between TCRs using a metric related to the physico-chemical properties of the loop residues predicted to interact with the epitope. By exploiting private and public datasets and comparing TCRpcDist with competing approaches, it is demonstrated that TCRpcDist can accurately identify groups of TCRs that are likely to bind the same epitopes. Importantly, the ability of TCRpcDist is experimentally validated to determine antigen specificities (neoantigens and tumor-associated antigens) of orphan tumor-infiltrating lymphocytes (TILs) in cancer patients. TCRpcDist is thus a promising approach to support TCR repertoire analysis and TCR deorphanization for individualized treatments including cancer immunotherapies.

Predictability of antigen binding based on short motifs in the antibody CDRH3

Adaptive immune receptors, such as antibodies and T-cell receptors, recognize foreign threats with exquisite specificity. A major challenge in adaptive immunology is discovering the rules governing immune receptor-antigen binding in order to predict the antigen binding status of previously unseen immune receptors. Many studies assume that the antigen binding status of an immune receptor may be determined by the presence of a short motif in the complementarity determining region 3 (CDR3), disregarding other amino acids. To test this assumption, we present a method to discover short motifs which show high precision in predicting antigen binding and generalize well to unseen simulated and experimental data. Our analysis of a mutagenesis-based antibody dataset reveals 11 336 position-specific, mostly gapped motifs of 3-5 amino acids that retain high precision on independently generated experimental data. Using a subset of only 178 motifs, a simple classifier was made that on the independently generated dataset outperformed a deep learning model proposed specifically for such datasets. In conclusion, our findings support the notion that for some antibodies, antigen binding may be largely determined by a short CDR3 motif. As more experimental data emerge, our methodology could serve as a foundation for in-depth investigations into antigen binding signals.

Computational detection of antigen-specific B cell receptors following immunization

B cell receptors (BCRs) play a crucial role in recognizing and fighting foreign antigens. High-throughput sequencing enables in-depth sampling of the BCRs repertoire after immunization. However, only a minor fraction of BCRs actively participate in any given infection. To what extent can we accurately identify antigen-specific sequences directly from BCRs repertoires? We present a computational method grounded on sequence similarity, aimed at identifying statistically significant responsive BCRs. This method leverages well-known characteristics of affinity maturation and expected diversity. We validate its effectiveness using longitudinally sampled human immune repertoire data following influenza vaccination and SARS-CoV-2 infections. We show that different lineages converge to the same responding Complementarity Determining Region 3, demonstrating convergent selection within an individual. The outcomes of this method hold promise for application in vaccine development, personalized medicine, and antibody-derived therapeutics.

Microbiota dictate T cell clonal selection to augment graft-versus-host disease after stem cell transplantation

Allogeneic T cell expansion is the primary determinant of graft-versus-host disease (GVHD), and current dogma dictates that this is driven by histocompatibility antigen disparities between donor and recipient. This paradigm represents a closed genetic system within which donor T cells interact with peptide-major histocompatibility complexes (MHCs), though clonal interrogation remains challenging due to the sparseness of the T cell repertoire. We developed a Bayesian model using donor and recipient T cell receptor (TCR) frequencies in murine stem cell transplant systems to define limited common expansion of T cell clones across genetically identical donor-recipient pairs. A subset of donor CD4+ T cell clonotypes differentially expanded in identical recipients and were microbiota dependent. Microbiota-specific T cells augmented GVHD lethality and could target microbial antigens presented by gastrointestinal epithelium during an alloreactive response. The microbiota serves as a source of cognate antigens that contribute to clonotypic T cell expansion and the induction of GVHD independent of donor-recipient genetics.

Aberrant adaptive immune response underlies genetic susceptibility to tuberculosis

Mycobacterium tuberculosis (Mtb) remains a major threat worldwide, although only a fraction of infected individuals develops tuberculosis (TB). TB susceptibility is shaped by multiple genetic factors, and we performed comparative immunological analysis of two mouse strains to uncover relevant mechanisms underlying susceptibility and resistance. C57BL/6 mice are relatively TB-resistant, whereas I/St mice are prone to develop severe TB, partly due to the MHC-II allelic variant that shapes suboptimal CD4+ T cell receptor repertoire. We investigated the repertoires of lung-infiltrating helper T cells and B cells at the progressed stage in both strains. We found that lung CD4+ T cell repertoires of infected C57BL/6 but not I/St mice contained convergent TCR clusters with functionally confirmed Mtb specificity. Transcriptomic analysis revealed a more prominent Th1 signature in C57BL/6, and expression of pro-inflammatory IL-16 in I/St lung-infiltrating helper T cells. The two strains also showed distinct Th2 signatures. Furthermore, the humoral response of I/St mice was delayed, less focused, and dominated by IgG/IgM isotypes, whereas C57BL/6 mice generated more Mtb antigen-focused IgA response. We conclude that the inability of I/St mice to produce a timely and efficient anti-Mtb adaptive immune responses arises from a suboptimal helper T cell landscape that also impacts the humoral response, leading to diffuse inflammation and severe disease.

Single-cell immune repertoire analysis

Single-cell T cell and B cell antigen receptor-sequencing data analysis can potentially perform in-depth assessments of adaptive immune cells that inform on understanding immune cell development to tracking clonal expansion in disease and therapy. However, it has been extremely challenging to analyze and interpret T cells and B cells and their adaptive immune receptor repertoires at the single-cell level due to not only the complexity of the data but also the underlying biology. In this Review, we delve into the computational breakthroughs that have transformed the analysis of single-cell T cell and B cell antigen receptor-sequencing data.

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.

DNAJB1-PRKACA fusion neoantigens elicit rare endogenous T cell responses that potentiate cell therapy for fibrolamellar carcinoma

Fibrolamellar carcinoma (FLC) is a liver tumor with a high mortality burden and few treatment options. A promising therapeutic vulnerability in FLC is its driver mutation, a conserved DNAJB1-PRKACA gene fusion that could be an ideal target neoantigen for immunotherapy. In this study, we aim to define endogenous CD8 T cell responses to this fusion in FLC patients and evaluate fusion-specific T cell receptors (TCRs) for use in cellular immunotherapies. We observe that fusion-specific CD8 T cells are rare and that FLC patient TCR repertoires lack large clusters of related TCR sequences characteristic of potent antigen-specific responses, potentially explaining why endogenous immune responses are insufficient to clear FLC tumors. Nevertheless, we define two functional fusion-specific TCRs, one of which has strong anti-tumor activity in vivo. Together, our results provide insights into the fragmented nature of neoantigen-specific repertoires in humans and indicate routes for clinical development of successful immunotherapies for FLC.

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.

H3K27M neoepitope vaccination in diffuse midline glioma induces B and T cell responses across diverse HLA loci of a recovered patient

H3K27M, a driver mutation with T and B cell neoepitope characteristics, defines an aggressive subtype of diffuse glioma with poor survival. We functionally dissect the immune response of one patient treated with an H3K27M peptide vaccine who subsequently entered complete remission. The vaccine robustly expanded class II human leukocyte antigen (HLA)-restricted peripheral H3K27M-specific T cells. Using functional assays, we characterized 34 clonally unique H3K27M-reactive T cell receptors and identified critical, conserved motifs in their complementarity-determining region 3 regions. Using detailed HLA mapping, we further demonstrate that diverse HLA-DQ and HLA-DR alleles present immunogenic H3K27M epitopes. Furthermore, we identified and profiled H3K27M-reactive B cell receptors from activated B cells in the cerebrospinal fluid. Our results uncover the breadth of the adaptive immune response against a shared clonal neoantigen across multiple HLA allelotypes and support the use of class II-restricted peptide vaccines to stimulate tumor-specific T and B cells harboring receptors with therapeutic potential.

Dynamics of T-cell receptor repertoire in patients with ankylosing spondylitis after biologic therapy

INTRODUCTION

Ankylosing spondylitis (AS) is a chronic inflammatory autoimmune disease in which T-cell immune responses play important roles. AS has been characterized by altered T-cell receptor (TCR) repertoire profiles, which are thought to be caused by expansion of disease-related TCR clonotypes. However, how biological agents affect the TCR repertoire status and whether their therapeutic outcomes are associated with certain features or dynamic patterns of the TCR repertoire are still elusive.

MATERIAL AND METHODS

We collected clinical samples from AS patients pre- and post-treatment with biologics. TCR repertoire sequencing was conducted to investigate associations of TCRα and TCRβ repertoire characteristics with disease activity and inflammatory indicators/cytokines.

RESULTS

Our results showed that good responders were associated with an increase in the TCR repertoire diversity with higher proportions of contracted TCR clonotypes. Additionally, we further identified a positive correlation between TCR repertoire diversity and interleukin (IL)-23 levels in AS patients. A network analysis revealed that contracted AS-associated TCR clonotypes with the same complementary-determining region 3 (CDR3) motifs, which represented high probabilities of sharing TCR specificities to AS-related antigens, were dominant in good responders of AS after treatment with biologic therapies.

CONCLUSIONS

Our findings suggested an important connection between TCR repertoire changes and therapeutic outcomes in biologic-treated AS patients. The status and dynamics of TCR repertoire profiles are useful for assessing the prognosis of biologic treatments in AS patients.

Dandelion uses the single-cell adaptive immune receptor repertoire to explore lymphocyte developmental origins

Assessment of single-cell gene expression (single-cell RNA sequencing) and adaptive immune receptor (AIR) sequencing (scVDJ-seq) has been invaluable in studying lymphocyte biology. Here we introduce Dandelion, a computational pipeline for scVDJ-seq analysis. It enables the application of standard V(D)J analysis workflows to single-cell datasets, delivering improved V(D)J contig annotation and the identification of nonproductive and partially spliced contigs. We devised a strategy to create an AIR feature space that can be used for both differential V(D)J usage analysis and pseudotime trajectory inference. The application of Dandelion improved the alignment of human thymic development trajectories of double-positive T cells to mature single-positive CD4/CD8 T cells, generating predictions of factors regulating lineage commitment. Dandelion analysis of other cell compartments provided insights into the origins of human B1 cells and ILC/NK cell development, illustrating the power of our approach. Dandelion is available at https://www.github.com/zktuong/dandelion .

Local and Global Variability in Developing Human T-Cell Repertoires

The adaptive immune response relies on T cells that combine phenotypic specialization with diversity of T-cell receptors (TCRs) to recognize a wide range of pathogens. TCRs are acquired and selected during T-cell maturation in the thymus. Characterizing TCR repertoires across individuals and T-cell maturation stages is important for better understanding adaptive immune responses and for developing new diagnostics and therapies. Analyzing a dataset of human TCR repertoires from thymocyte subsets, we find that the variability between individuals generated during the TCR V(D)J recombination is maintained through all stages of T-cell maturation and differentiation. The interindividual variability of repertoires of the same cell type is of comparable magnitude to the variability across cell types within the same individual. To zoom in on smaller scales than whole repertoires, we defined a distance measuring the relative overlap of locally similar sequences in repertoires. We find that the whole repertoire models correctly predict local similarity networks, suggesting a lack of forbidden T-cell receptor sequences. The local measure correlates well with distances calculated using whole repertoire traits and carries information about cell types.

Steering and controlling evolution - from bioengineering to fighting pathogens

Control interventions steer the evolution of molecules, viruses, microorganisms or other cells towards a desired outcome. Applications range from engineering biomolecules and synthetic organisms to drug, therapy and vaccine design against pathogens and cancer. In all these instances, a control system alters the eco-evolutionary trajectory of a target system, inducing new functions or suppressing escape evolution. Here, we synthesize the objectives, mechanisms and dynamics of eco-evolutionary control in different biological systems. We discuss how the control system learns and processes information about the target system by sensing or measuring, through adaptive evolution or computational prediction of future trajectories. This information flow distinguishes pre-emptive control strategies by humans from feedback control in biotic systems. We establish a cost-benefit calculus to gauge and optimize control protocols, highlighting the fundamental link between predictability of evolution and efficacy of pre-emptive control.

Determining epitope specificity of T-cell receptors with transformers

SUMMARY

T-cell receptors (TCRs) on T cells recognize and bind to epitopes presented by the major histocompatibility complex in case of an infection or cancer. However, the high diversity of TCRs, as well as their unique and complex binding mechanisms underlying epitope recognition, make it difficult to predict the binding between TCRs and epitopes. Here, we present the utility of transformers, a deep learning strategy that incorporates an attention mechanism that learns the informative features, and show that these models pre-trained on a large set of protein sequences outperform current strategies. We compared three pre-trained auto-encoder transformer models (ProtBERT, ProtAlbert, and ProtElectra) and one pre-trained auto-regressive transformer model (ProtXLNet) to predict the binding specificity of TCRs to 25 epitopes from the VDJdb database (human and murine). Two additional modifications were performed to incorporate gene usage of the TCRs in the four transformer models. Of all 12 transformer implementations (four models with three different modifications), a modified version of the ProtXLNet model could predict TCR-epitope pairs with the highest accuracy (weighted F1 score 0.55 simultaneously considering all 25 epitopes). The modification included additional features representing the gene names for the TCRs. We also showed that the basic implementation of transformers outperformed the previously available methods, i.e. TCRGP, TCRdist, and DeepTCR, developed for the same biological problem, especially for the hard-to-classify labels. We show that the proficiency of transformers in attention learning can be made operational in a complex biological setting like TCR binding prediction. Further ingenuity in utilizing the full potential of transformers, either through attention head visualization or introducing additional features, can extend T-cell research avenues.

AVAILABILITY AND IMPLEMENTATION

Data and code are available on https://github.com/InduKhatri/tcrformer.

Targeted depletion of TRBV9+ T cells as immunotherapy in a patient with ankylosing spondylitis

Autoimmunity is intrinsically driven by memory T and B cell clones inappropriately targeted at self-antigens. Selective depletion or suppression of self-reactive T cells remains a holy grail of autoimmune therapy, but disease-associated T cell receptors (TCRs) and cognate antigenic epitopes remained elusive. A TRBV9-containing CD8+ TCR motif was recently associated with the pathogenesis of ankylosing spondylitis, psoriatic arthritis and acute anterior uveitis, and cognate HLA-B*27-presented epitopes were identified. Following successful testing in nonhuman primate models, here we report human TRBV9+ T cell elimination in ankylosing spondylitis. The patient achieved remission within 3 months and ceased anti-TNF therapy after 5 years of continuous use. Complete remission has now persisted for 4 years, with three doses of anti-TRBV9 administered per year. We also observed a profound improvement in spinal mobility metrics and the Bath Ankylosing Spondylitis Metrology Index (BASMI). This represents a possibly curative therapy of an autoimmune disease via selective depletion of a TRBV-defined group of T cells. The anti-TRBV9 therapy could potentially be applicable to other HLA-B*27-associated spondyloarthropathies. Such targeted elimination of the underlying cause of the disease without systemic immunosuppression could offer a new generation of safe and efficient therapies for autoimmunity.

T Cell Receptor Chain Centricity: The Phenomenon and Potential Applications in Cancer Immunotherapy

T cells are crucial players in adaptive anti-cancer immunity. The gene modification of T cells with tumor antigen-specific T cell receptors (TCRs) was a milestone in personalized cancer immunotherapy. TCR is a heterodimer (either α/β or γ/δ) able to recognize a peptide antigen in a complex with self-MHC molecules. Although traditional concepts assume that an α- and β-chain contribute equally to antigen recognition, mounting data reveal that certain receptors possess chain centricity, i.e., one hemi-chain TCR dominates antigen recognition and dictates its specificity. Chain-centric TCRs are currently poorly understood in terms of their origin and the functional T cell subsets that express them. In addition, the ratio of α- and β-chain-centric TCRs, as well as the exact proportion of chain-centric TCRs in the native repertoire, is generally still unknown today. In this review, we provide a retrospective analysis of studies that evidence chain-centric TCRs, propose patterns of their generation, and discuss the potential applications of such receptors in T cell gene modification for adoptive cancer immunotherapy.

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.

Novel bimodal TRBD1-TRBD2 rearrangements with dual or absent D-region contribute to TRB V-(D)-J combinatorial diversity

T-cell receptor (TR) diversity of the variable domains is generated by recombination of both the alpha (TRA) and beta (TRB) chains. The textbook process of TRB chain production starts with TRBD and TRBJ gene rearrangement, followed by the rearrangement of a TRBV gene to the partially rearranged D-J gene. Unsuccessful V-D-J TRB rearrangements lead to apoptosis of the cell. Here, we performed deep sequencing of the poorly explored pool of partial TRBD1-TRBD2 rearrangements in T-cell genomic DNA. We reconstructed full repertoires of human partial TRBD1-TRBD2 rearrangements using novel sequencing and validated them by detecting V-D-J recombination-specific byproducts: excision circles containing the recombination signal (RS) joint 5'D2-RS - 3'D1-RS. Identified rearrangements were in compliance with the classical 12/23 rule, common for humans, rats, and mice and contained typical V-D-J recombination footprints. Interestingly, we detected a bimodal distribution of D-D junctions indicating two active recombination sites producing long and short D-D rearrangements. Long TRB D-D rearrangements with two D-regions are coding joints D1-D2 remaining classically on the chromosome. The short TRB D-D rearrangements with no D-region are signal joints, the coding joint D1-D2 being excised from the chromosome. They both contribute to the TRB V-(D)-J combinatorial diversity. Indeed, short D-D rearrangements may be followed by direct V-J2 recombination. Long D-D rearrangements may recombine further with J2 and V genes forming partial D1-D2-J2 and then complete V-D1-D2-J2 rearrangement. Productive TRB V-D1-D2-J2 chains are present and expressed in thousands of clones of human antigen-experienced memory T cells proving their capacity for antigen recognition and actual participation in the immune response.

The landscape of the immunoglobulin repertoire in endemic pemphigus foliaceus

Introduction

Primarily driven by autoreactive B cells, pemphigus foliaceus (PF) is an uncommon autoimmune blistering skin disease of sporadic occurrence worldwide. However, PF reaches a prevalence of 3% in the endemic areas of Brazil, the highest ever registered for any autoimmune disease, which indicates environmental factors influencing the immune response in susceptible individuals. We aimed to provide insights into the immune repertoire of patients with PF living in the endemic region of the disease, compared to healthy individuals from the endemic region and a non-endemic area.

Methods

We characterized the B-cell repertoire in i) nontreated patients (n=5); ii) patients under immunosuppressive treatment (n=5); iii) patients in remission without treatment (n=6); and two control groups iv) from the endemic (n=6) and v) non-endemic areas in Brazil (n=4). We used total RNA extracted from peripheral blood mononuclear cells and performed a comprehensive characterization of the variable region of immunoglobulin heavy chain (IGH) in IgG and IgM using next-generation sequencing.

Results

Compared to individuals from a different area, we observed remarkably lower clonotype diversity in the B-cell immune repertoire of patients and controls from the endemic area (p < 0.02), suggesting that the immune repertoire in the endemic area is under geographically specific and intense environmental pressure. Moreover, we observed longer CDR3 sequences in patients, and we identified differential disease-specific usage of IGHV segments, including increased IGHV3-30 and decreased IGHV3-23 in patients with active disease (p < 0.04). Finally, our robust network analysis discovered clusters of CDR3 sequences uniquely observed in patients with PF.

Discussion

Our results indicate that environmental factors, in addition to disease state, impact the characteristics of the repertoire. Our findings can be applied to further investigation of the environmental factors that trigger pemphigus and expand the knowledge for identifying new targeted and more effective therapies.

Inhibitory IL-10-producing CD4+ T cells are T-bet-dependent and facilitate cytomegalovirus persistence via coexpression of arginase-1

Inhibitory CD4+ T cells have been linked with suboptimal immune responses against cancer and pathogen chronicity. However, the mechanisms that underpin the development of these regulatory cells, especially in the context of ongoing antigen exposure, have remained obscure. To address this knowledge gap, we undertook a comprehensive functional, phenotypic, and transcriptomic analysis of interleukin (IL)-10-producing CD4+ T cells induced by chronic infection with murine cytomegalovirus (MCMV). We identified these cells as clonally expanded and highly differentiated TH1-like cells that developed in a T-bet-dependent manner and coexpressed arginase-1 (Arg1), which promotes the catalytic breakdown of L-arginine. Mice lacking Arg1-expressing CD4+ T cells exhibited more robust antiviral immunity and were better able to control MCMV. Conditional deletion of T-bet in the CD4+ lineage suppressed the development of these inhibitory cells and also enhanced immune control of MCMV. Collectively, these data elucidated the ontogeny of IL-10-producing CD4+ T cells and revealed a previously unappreciated mechanism of immune regulation, whereby viral persistence was facilitated by the site-specific delivery of Arg1.

Unique features of the TCR repertoire of reactivated memory T cells in the experimental mouse tumor model

T cell engineering with T cell receptors (TCR) specific to tumor antigens has become a breakthrough towards personalized cancer adoptive cell immunotherapy. However, the search for therapeutic TCRs is often challenging, and effective strategies are strongly required for the identification and enrichment of tumor-specific T cells that express TCRs with superior functional characteristics. Using an experimental mouse tumor model, we studied sequential changes in TCR repertoire features of T cells involved in the primary and secondary immune responses to allogeneic tumor antigens. In-depth bioinformatics analysis of TCR repertoires showed differences in reactivated memory T cells compared to primarily activated effectors. After cognate antigen re-encounter, memory cells were enriched with clonotypes that express α-chain TCR with high potential cross-reactivity and enhanced strength of interaction with both MHC and docked peptides. Our findings suggest that functionally true memory T cells could be a better source of therapeutic TCRs for adoptive cell therapy. No marked changes were observed in the physicochemical characteristics of TCRβ in reactivated memory clonotypes, indicative of the dominant role of TCRα in the secondary allogeneic immune response. The results of this study could further contribute to the development of TCR-modified T cell products based on the phenomenon of TCR chain centricity.

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.

TCR sequencing: applications in immuno-oncology research

•T-cell receptor (TCR) interaction with major histocompatibility complex-antigen complexes leads to antitumour responses.•TCR sequencing analysis allows characterisation of T cells that recognise tumour neoantigens.•T-cell clonal revival and clonal replacement potentially underpin immunotherapy responses.

TCR_Explore: A novel webtool for T cell receptor repertoire analysis

T cells expressing either alpha-beta or gamma-delta T cell receptors (TCR) are critical sentinels of the adaptive immune system, with receptor diversity being essential for protective immunity against a broad array of pathogens and agents. Programs available to profile TCR clonotypic signatures can be limiting for users with no coding expertise. Current analytical pipelines can be inefficient due to manual processing steps, open to data entry errors and have multiple analytical tools with unique inputs that require coding expertise. Here we present a bespoke webtool designed for users irrespective of coding expertise, coined 'TCR_Explore', enabling analysis either derived via Sanger sequencing or next generation sequencing (NGS) platforms. Further, TCR_Explore incorporates automated quality control steps for Sanger sequencing. The creation of flexible and publication ready figures are enabled for different sequencing platforms following universal conversion to the TCR_Explore file format. TCR_Explore will enhance a user's capacity to undertake in-depth TCR repertoire analysis of both new and pre-existing datasets for identification of T cell clonotypes associated with health and disease. The web application is located at https://tcr-explore.erc.monash.edu for users to interactively explore TCR repertoire datasets.

Modeling and predicting the overlap of B- and T-cell receptor repertoires in healthy and SARS-CoV-2 infected individuals

Adaptive immunity's success relies on the extraordinary diversity of protein receptors on B and T cell membranes. Despite this diversity, the existence of public receptors shared by many individuals gives hope for developing population-wide vaccines and therapeutics. Using probabilistic modeling, we show many of these public receptors are shared by chance in healthy individuals. This predictable overlap is driven not only by biases in the random generation process of receptors, as previously reported, but also by their common functional selection. However, the model underestimates sharing between repertoires of individuals infected with SARS-CoV-2, suggesting strong specific antigen-driven convergent selection. We exploit this discrepancy to identify COVID-associated receptors, which we validate against datasets of receptors with known viral specificity. We study their properties in terms of sequence features and network organization, and use them to design an accurate diagnostic tool for predicting SARS-CoV-2 status from repertoire data.

Compartmentalization and persistence of dominant (regulatory) T cell clones indicates antigen skewing in juvenile idiopathic arthritis

Autoimmune inflammation is characterized by tissue infiltration and expansion of antigen-specific T cells. Although this inflammation is often limited to specific target tissues, it remains yet to be explored whether distinct affected sites are infiltrated with the same, persistent T cell clones. Here, we performed CyTOF analysis and T cell receptor (TCR) sequencing to study immune cell composition and (hyper-)expansion of circulating and joint-derived Tregs and non-Tregs in juvenile idiopathic arthritis (JIA). We studied different joints affected at the same time, as well as over the course of relapsing-remitting disease. We found that the composition and functional characteristics of immune infiltrates are strikingly similar between joints within one patient, and observed a strong overlap between dominant T cell clones, especially Treg, of which some could also be detected in circulation and persisted over the course of relapsing-remitting disease. Moreover, these T cell clones were characterized by a high degree of sequence similarity, indicating the presence of TCR clusters responding to the same antigens. These data suggest that in localized autoimmune disease, there is autoantigen-driven expansion of both Teffector and Treg clones that are highly persistent and are (re)circulating. These dominant clones might represent interesting therapeutic targets.

T cell receptor repertoires associated with control and disease progression following Mycobacterium tuberculosis infection

Antigen-specific, MHC-restricted αβ T cells are necessary for protective immunity against Mycobacterium tuberculosis, but the ability to broadly study these responses has been limited. In the present study, we used single-cell and bulk T cell receptor (TCR) sequencing and the GLIPH2 algorithm to analyze M. tuberculosis-specific sequences in two longitudinal cohorts, comprising 166 individuals with M. tuberculosis infection who progressed to either tuberculosis (n = 48) or controlled infection (n = 118). We found 24 T cell groups with similar TCR-β sequences, predicted by GLIPH2 to have common TCR specificities, which were associated with control of infection (n = 17), and others that were associated with progression to disease (n = 7). Using a genome-wide M. tuberculosis antigen screen, we identified peptides targeted by T cell similarity groups enriched either in controllers or in progressors. We propose that antigens recognized by T cell similarity groups associated with control of infection can be considered as high-priority targets for future vaccine development.

simAIRR: simulation of adaptive immune repertoires with realistic receptor sequence sharing for benchmarking of immune state prediction methods

BACKGROUND

Machine learning (ML) has gained significant attention for classifying immune states in adaptive immune receptor repertoires (AIRRs) to support the advancement of immunodiagnostics and therapeutics. Simulated data are crucial for the rigorous benchmarking of AIRR-ML methods. Existing approaches to generating synthetic benchmarking datasets result in the generation of naive repertoires missing the key feature of many shared receptor sequences (selected for common antigens) found in antigen-experienced repertoires.

RESULTS

We demonstrate that a common approach to generating simulated AIRR benchmark datasets can introduce biases, which may be exploited for undesired shortcut learning by certain ML methods. To mitigate undesirable access to true signals in simulated AIRR datasets, we devised a simulation strategy (simAIRR) that constructs antigen-experienced-like repertoires with a realistic overlap of receptor sequences. simAIRR can be used for constructing AIRR-level benchmarks based on a range of assumptions (or experimental data sources) for what constitutes receptor-level immune signals. This includes the possibility of making or not making any prior assumptions regarding the similarity or commonality of immune state-associated sequences that will be used as true signals. We demonstrate the real-world realism of our proposed simulation approach by showing that basic ML strategies perform similarly on simAIRR-generated and real-world experimental AIRR datasets.

CONCLUSIONS

This study sheds light on the potential shortcut learning opportunities for ML methods that can arise with the state-of-the-art way of simulating AIRR datasets. simAIRR is available as a Python package: https://github.com/KanduriC/simAIRR.

Comparing T cell receptor repertoires using optimal transport

The complexity of entire T cell receptor (TCR) repertoires makes their comparison a difficult but important task. Current methods of TCR repertoire comparison can incur a high loss of distributional information by considering overly simplistic sequence- or repertoire-level characteristics. Optimal transport methods form a suitable approach for such comparison given some distance or metric between values in the sample space, with appealing theoretical and computational properties. In this paper we introduce a nonparametric approach to comparing empirical TCR repertoires that applies the Sinkhorn distance, a fast, contemporary optimal transport method, and a recently-created distance between TCRs called TCRdist. We show that our methods identify meaningful differences between samples from distinct TCR distributions for several case studies, and compete with more complicated methods despite minimal modeling assumptions and a simpler pipeline.

T cell receptor convergence is an indicator of antigen-specific T cell response in cancer immunotherapies

T cells are potent at eliminating pathogens and playing a crucial role in the adaptive immune response. T cell receptor (TCR) convergence describes T cells that share identical TCRs with the same amino acid sequences but have different DNA sequences due to codon degeneracy. We conducted a systematic investigation of TCR convergence using single-cell immune profiling and bulk TCRβ-sequence (TCR-seq) data obtained from both mouse and human samples and uncovered a strong link between antigen-specificity and convergence. This association was stronger than T cell expansion, a putative indicator of antigen-specific T cells. By using flow-sorted tetramer+ single T cell data, we discovered that convergent T cells were enriched for a neoantigen-specific CD8+ effector phenotype in the tumor microenvironment. Moreover, TCR convergence demonstrated better prediction accuracy for immunotherapy response than the existing TCR repertoire indexes. In conclusion, convergent T cells are likely to be antigen-specific and might be a novel prognostic biomarker for anti-cancer immunotherapy.

TCR repertoire profiling revealed antigen-driven CD8+ T cell clonal groups shared in synovial fluid of patients with spondyloarthritis

Spondyloarthritis (SpA) comprises a number of inflammatory rheumatic diseases with overlapping clinical manifestations. Strong association with several HLA-I alleles and T cell infiltration into an inflamed joint suggest involvement of T cells in SpA pathogenesis. In this study, we performed high-throughput T cell repertoire profiling of synovial fluid (SF) and peripheral blood (PB) samples collected from a large cohort of SpA patients. We showed that synovial fluid is enriched with expanded T cell clones that are shared between patients with similar HLA genotypes and persist during recurrent synovitis. Using an algorithm for identification of TCRs involved in immune response we discovered several antigen-driven CD8+ clonal groups associated with risk HLA-B*27 or HLA-B*38 alleles. We further show that these clonal groups were enriched in SF and had higher frequency in PB of SpA patients vs healthy donors, implying their relevance to SpA pathogenesis. Several of the groups were shared among patients with different SpAs that suggests a common immunopathological mechanism of the diseases. In summary, our results provide evidence for the role of specific CD8+ T cell clones in pathogenesis of SpA.

CompAIRR: ultra-fast comparison of adaptive immune receptor repertoires by exact and approximate sequence matching

MOTIVATION

Adaptive immune receptor (AIR) repertoires (AIRRs) record past immune encounters with exquisite specificity. Therefore, identifying identical or similar AIR sequences across individuals is a key step in AIRR analysis for revealing convergent immune response patterns that may be exploited for diagnostics and therapy. Existing methods for quantifying AIRR overlap scale poorly with increasing dataset numbers and sizes. To address this limitation, we developed CompAIRR, which enables ultra-fast computation of AIRR overlap, based on either exact or approximate sequence matching.

RESULTS

CompAIRR improves computational speed 1000-fold relative to the state of the art and uses only one-third of the memory: on the same machine, the exact pairwise AIRR overlap of 104 AIRRs with 105 sequences is found in ∼17 min, while the fastest alternative tool requires 10 days. CompAIRR has been integrated with the machine learning ecosystem immuneML to speed up commonly used AIRR-based machine learning applications.

AVAILABILITY AND IMPLEMENTATION

CompAIRR code and documentation are available at https://github.com/uio-bmi/compairr. Docker images are available at https://hub.docker.com/r/torognes/compairr. The code to replicate the synthetic datasets, scripts for benchmarking and creating figures, and all raw data underlying the figures are available at https://github.com/uio-bmi/compairr-benchmarking.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Inactivated tick-borne encephalitis vaccine elicits several overlapping waves of T cell response

The development and implementation of vaccines have been growing exponentially, remaining one of the major successes of healthcare over the last century. Nowadays, active regular immunizations prevent epidemics of many viral diseases, including tick-borne encephalitis (TBE). Along with the generation of virus-specific antibodies, a highly effective vaccine should induce T cell responses providing long-term immune defense. In this study, we performed longitudinal high-throughput T cell receptor (TCR) sequencing to characterize changes in individual T cell repertoires of 11 donors immunized with an inactivated TBE vaccine. After two-step immunization, we found significant clonal expansion of both CD4+ and CD8+ T cells, ranging from 302 to 1706 vaccine-associated TCRβ clonotypes in different donors. We detected several waves of T cell clonal expansion generated by distinct groups of vaccine-responding clones. Both CD4+ and CD8+ vaccine-responding T cell clones formed 17 motifs in TCRβ sequences shared by donors with identical HLA alleles. Our results indicate that TBE vaccination leads to a robust T cell response due to the production of a variety of T cell clones with a memory phenotype, which recognize a large set of epitopes.

Reference-based comparison of adaptive immune receptor repertoires

B and T cell receptor (immune) repertoires can represent an individual's immune history. While current repertoire analysis methods aim to discriminate between health and disease states, they are typically based on only a limited number of parameters. Here, we introduce immuneREF: a quantitative multidimensional measure of adaptive immune repertoire (and transcriptome) similarity that allows interpretation of immune repertoire variation by relying on both repertoire features and cross-referencing of simulated and experimental datasets. To quantify immune repertoire similarity landscapes across health and disease, we applied immuneREF to >2,400 datasets from individuals with varying immune states (healthy, [autoimmune] disease, and infection). We discovered, in contrast to the current paradigm, that blood-derived immune repertoires of healthy and diseased individuals are highly similar for certain immune states, suggesting that repertoire changes to immune perturbations are less pronounced than previously thought. In conclusion, immuneREF enables the population-wide study of adaptive immune response similarity across immune states.

Resolving SARS-CoV-2 CD4+ T cell specificity via reverse epitope discovery

The current strategy to detect immunodominant T cell responses focuses on the antigen, employing large peptide pools to screen for functional cell activation. However, these approaches are labor and sample intensive and scale poorly with increasing size of the pathogen peptidome. T cell receptors (TCRs) recognizing the same epitope frequently have highly similar sequences, and thus, the presence of large sequence similarity clusters in the TCR repertoire likely identify the most public and immunodominant responses. Here, we perform a meta-analysis of large, publicly available single-cell and bulk TCR datasets from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-infected individuals to identify public CD4+ responses. We report more than 1,200 αβTCRs forming six prominent similarity clusters and validate histocompatibility leukocyte antigen (HLA) restriction and epitope specificity predictions for five clusters using transgenic T cell lines. Collectively, these data provide information on immunodominant CD4+ T cell responses to SARS-CoV-2 and demonstrate the utility of the reverse epitope discovery approach.

Machine Learning Approaches to TCR Repertoire Analysis

Sparked by the development of genome sequencing technology, the quantity and quality of data handled in immunological research have been changing dramatically. Various data and database platforms are now driving the rapid progress of machine learning for immunological data analysis. Of various topics in immunology, T cell receptor repertoire analysis is one of the most important targets of machine learning for assessing the state and abnormalities of immune systems. In this paper, we review recent repertoire analysis methods based on machine learning and deep learning and discuss their prospects.

Epidemiology of Δ8THC-Related Carcinogenesis in USA: A Panel Regression and Causal Inferential Study

The use of Δ8THC is increasing at present across the USA in association with widespread cannabis legalization and the common notion that it is "legal weed". As genotoxic actions have been described for many cannabinoids, we studied the cancer epidemiology of Δ8THC. Data on 34 cancer types was from the Centers for Disease Control Atlanta Georgia, substance abuse data from the Substance Abuse and Mental Health Services Administration, ethnicity and income data from the U.S. Census Bureau, and cannabinoid concentration data from the Drug Enforcement Agency, were combined and processed in R. Eight cancers (corpus uteri, liver, gastric cardia, breast and post-menopausal breast, anorectum, pancreas, and thyroid) were related to Δ8THC exposure on bivariate testing, and 18 (additionally, stomach, Hodgkins, and Non-Hodgkins lymphomas, ovary, cervix uteri, gall bladder, oropharynx, bladder, lung, esophagus, colorectal cancer, and all cancers (excluding non-melanoma skin cancer)) demonstrated positive average marginal effects on fully adjusted inverse probability weighted interactive panel regression. Many minimum E-Values (mEVs) were infinite. p-values rose from 8.04 × 10-78. Marginal effect calculations revealed that 18 Δ8THC-related cancers are predicted to lead to a further 8.58 cases/100,000 compared to 7.93 for alcoholism and -8.48 for tobacco. Results indicate that between 8 and 20/34 cancer types were associated with Δ8THC exposure, with very high effect sizes (mEVs) and marginal effects after adjustment exceeding tobacco and alcohol, fulfilling the epidemiological criteria of causality and suggesting a cannabinoid class effect. The inclusion of pediatric leukemias and testicular cancer herein demonstrates heritable malignant teratogenesis.

Profiling the baseline performance and limits of machine learning models for adaptive immune receptor repertoire classification

BACKGROUND

Machine learning (ML) methodology development for the classification of immune states in adaptive immune receptor repertoires (AIRRs) has seen a recent surge of interest. However, so far, there does not exist a systematic evaluation of scenarios where classical ML methods (such as penalized logistic regression) already perform adequately for AIRR classification. This hinders investigative reorientation to those scenarios where method development of more sophisticated ML approaches may be required.

RESULTS

To identify those scenarios where a baseline ML method is able to perform well for AIRR classification, we generated a collection of synthetic AIRR benchmark data sets encompassing a wide range of data set architecture-associated and immune state-associated sequence patterns (signal) complexity. We trained ≈1,700 ML models with varying assumptions regarding immune signal on ≈1,000 data sets with a total of ≈250,000 AIRRs containing ≈46 billion TCRβ CDR3 amino acid sequences, thereby surpassing the sample sizes of current state-of-the-art AIRR-ML setups by two orders of magnitude. We found that L1-penalized logistic regression achieved high prediction accuracy even when the immune signal occurs only in 1 out of 50,000 AIR sequences.

CONCLUSIONS

We provide a reference benchmark to guide new AIRR-ML classification methodology by (i) identifying those scenarios characterized by immune signal and data set complexity, where baseline methods already achieve high prediction accuracy, and (ii) facilitating realistic expectations of the performance of AIRR-ML models given training data set properties and assumptions. Our study serves as a template for defining specialized AIRR benchmark data sets for comprehensive benchmarking of AIRR-ML methods.

Utility of Bulk T-Cell Receptor Repertoire Sequencing Analysis in Understanding Immune Responses to COVID-19

Measuring immunity to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of coronavirus disease 19 (COVID-19), can rely on antibodies, reactive T cells and other factors, with T-cell-mediated responses appearing to have greater sensitivity and longevity. Because each T cell carries an essentially unique nucleic acid sequence for its T-cell receptor (TCR), we can interrogate sequence data derived from DNA or RNA to assess aspects of the immune response. This review deals with the utility of bulk, rather than single-cell, sequencing of TCR repertoires, considering the importance of study design, in terms of cohort selection, laboratory methods and analysis. The advances in understanding SARS-CoV-2 immunity that have resulted from bulk TCR repertoire sequencing are also be discussed. The complexity of sequencing data obtained by bulk repertoire sequencing makes analysis challenging, but simple descriptive analyses, clonal analysis, searches for specific sequences associated with immune responses to SARS-CoV-2, motif-based analyses, and machine learning approaches have all been applied. TCR repertoire sequencing has demonstrated early expansion followed by contraction of SARS-CoV-2-specific clonotypes, during active infection. Maintenance of TCR repertoire diversity, including the maintenance of diversity of anti-SARS-CoV-2 response, predicts a favourable outcome. TCR repertoire narrowing in severe COVID-19 is most likely a consequence of COVID-19-associated lymphopenia. It has been possible to follow clonotypic sequences longitudinally, which has been particularly valuable for clonotypes known to be associated with SARS-CoV-2 peptide/MHC tetramer binding or with SARS-CoV-2 peptide-induced cytokine responses. Closely related clonotypes to these previously identified sequences have been shown to respond with similar kinetics during infection. A possible superantigen-like effect of the SARS-CoV-2 spike protein has been identified, by means of observing V-segment skewing in patients with severe COVID-19, together with structural modelling. Such a superantigen-like activity, which is apparently absent from other coronaviruses, may be the basis of multisystem inflammatory syndrome and cytokine storms in COVID-19. Bulk TCR repertoire sequencing has proven to be a useful and cost-effective approach to understanding interactions between SARS-CoV-2 and the human host, with the potential to inform the design of therapeutics and vaccines, as well as to provide invaluable pathogenetic and epidemiological insights.

Pinpointing the tumor-specific T-cells via TCR clusters

Adoptive cell transfer (ACT) is a promising approach to cancer immunotherapy, but its efficiency fundamentally depends on the extent of tumor-specific T cell enrichment within the graft. This can be estimated via activation with identifiable neoantigens, tumor-associated antigens (TAAs), or living or lysed tumor cells, but these approaches remain laborious, time-consuming, and functionally limited, hampering clinical development of ACT. Here, we demonstrate that homology cluster analysis of T cell receptor (TCR) repertoires efficiently identifies tumor-reactive TCRs allowing to: (1) detect their presence within the pool of tumor-infiltrating lymphocytes (TILs); (2) optimize TIL culturing conditions, with IL-2_low/IL-21/anti-PD-1 combination showing increased efficiency; (3) investigate surface marker-based enrichment for tumor-targeting T cells in freshly isolated TILs (enrichment confirmed for CD4⁺ and CD8⁺ PD-1⁺/CD39⁺ subsets), or re-stimulated TILs (informs on enrichment in 4-1BB-sorted cells). We believe that this approach to the rapid assessment of tumor-specific TCR enrichment should accelerate T cell therapy development.