The T cell receptor (TCR) repertoire is an essential component of the CD8 T-cell immune response. Here, we seek to investigate factors that drive selection of TCR repertoires specific to the HLA-A2-restricted immunodominant epitope BRLF1109-117 (YVLDHLIVV) over the course of primary Epstein Barr virus (EBV) infection. Using single-cell paired TCRαβ sequencing of tetramer sorted CD8 T cells ex vivo, we show at the clonal level that recognition of the HLA-A2-restricted BRLF1 (YVL-BR, BRLF-1109) epitope is mainly driven by the TCRα chain. For the first time, we identify a CDR3α (complementarity determining region 3 α) motif, KDTDKL, resulting from an obligate AV8.1-AJ34 pairing that was shared by all four individuals studied. This observation coupled with the fact that this public AV8.1-KDTDKL-AJ34 TCR pairs with multiple different TCRβ chains within the same donor (median 4; range: 1–9), suggests that there are some unique structural features of the interaction between the YVL-BR/MHC and the AV8.1-KDTDKL-AJ34 TCR that leads to this high level of selection. Newly developed TCR motif algorithms identified a lysine at position 1 of the CDR3α motif that is highly conserved and likely important for antigen recognition. Crystal structure analysis of the YVL-BR/HLA-A2 complex revealed that the MHC-bound peptide bulges at position 4, exposing a negatively charged aspartic acid that may interact with the positively charged lysine of CDR3α. TCR cloning and site-directed mutagenesis of the CDR3α lysine ablated YVL-BR-tetramer staining and substantially reduced CD69 upregulation on TCR mutant-transduced cells following antigen-specific stimulation. Reduced activation of T cells expressing this CDR3 motif was also observed following exposure to mutated (D4A) peptide. In summary, we show that a highly public TCR repertoire to an immunodominant epitope of a common human virus is almost completely selected on the basis of CDR3α and provide a likely structural basis for the selection. These studies emphasize the importance of examining TCRα, as well as TCRβ, in understanding the CD8 T cell receptor repertoire.

Recognition modes of individual T cell receptors (TCRs) are well studied, but factors driving the selection of TCR repertoires from primary through persistent human virus infections are less well understood. Using deep sequencing, we demonstrate a high degree of diversity of Epstein-Barr virus (EBV)-specific clonotypes in acute infectious mononucleosis (AIM). Only 9% of unique clonotypes detected in AIM persisted into convalescence; the majority (91%) of unique clonotypes detected in AIM were not detected in convalescence and were seeming replaced by equally diverse “de novo” clonotypes. The persistent clonotypes had a greater probability of being generated than nonpersistent clonotypes due to convergence recombination of multiple nucleotide sequences to encode the same amino acid sequence, as well as the use of shorter complementarity-determining regions 3 (CDR3s) with fewer nucleotide additions (i.e., sequences closer to germ line). Moreover, the two most immunodominant HLA-A2-restricted EBV epitopes, BRLF1109 and BMLF1280, show highly distinct antigen-specific public (i.e., shared between individuals) features. In fact, TCRα CDR3 motifs played a dominant role, while TCRβ played a minimal role, in the selection of TCR repertoire to an immunodominant EBV epitope, BRLF1. This contrasts with the majority of previously reported repertoires, which appear to be selected either on TCRβ CDR3 interactions with peptide/major histocompatibility complex (MHC) or in combination with TCRα CDR3. Understanding of how TCR-peptide-MHC complex interactions drive repertoire selection can be used to develop optimal strategies for vaccine design or generation of appropriate adoptive immunotherapies for viral infections in transplant settings or for cancer. IMPORTANCE Several lines of evidence suggest that TCRα and TCRβ repertoires play a role in disease outcomes and treatment strategies during viral infections in transplant patients and in cancer and autoimmune disease therapy. Our data suggest that it is essential that we understand the basic principles of how to drive optimum repertoires for both TCR chains, α and β. We address this important issue by characterizing the CD8 TCR repertoire to a common persistent human viral infection (EBV), which is controlled by appropriate CD8 T cell responses. The ultimate goal would be to determine if the individuals who are infected asymptomatically develop a different TCR repertoire than those that develop the immunopathology of AIM. Here, we begin by doing an in-depth characterization of both CD8 T cell TCRα and TCRβ repertoires to two immunodominant EBV epitopes over the course of AIM, identifying potential factors that may be driving their selection.

Virus specific CD8+ T cells expand dramatically during acute Epstein Barr virus (EBV) infection, and their persistence is important for lifelong control of EBV-related disease. To better define the generation and maintenance of these effective CD8+ T cell responses, we used microarrays to characterize gene expression in total and EBV-specific CD8+ T cells isolated from the peripheral blood of ten individuals followed from acute infectious mononucleosis (AIM) into convalescence (CONV). In total CD8+ T cells, differential expression of genes in AIM and CONV was most pronounced among those encoding proteins important in T cell activation/differentiation, cell division/metabolism, chemokines/cytokines and receptors, signaling and transcription factors (TF), immune effector functions, and negative regulators. Within these categories, we identified 28 genes that correlated with CD8+ T cell expansion in response to an acute EBV infection. In EBV-specific CD8+ T cells, we identified 33 genes that were differentially expressed in AIM and CONV. Two important TF, T-bet and Eomesodermin (Eomes), were upregulated and maintained at similar levels in both AIM and CONV; by contrast, protein expression declined from AIM to CONV. Expression of these TF varied among cells with different epitope specificities. Altogether, gene and protein expression patterns suggest that a large proportion, if not a majority of CD8+ T cells in AIM are virus-specific, activated, dividing, and primed to exert effector activities. High expression of T-bet and Eomes may help to maintain effector mechanisms in activated cells, and to enable proliferation and transition to earlier differentiation states in CONV.

Background With more T cell receptor sequence data becoming available, the need for bioinformatics approaches to predict T cell receptor specificity is even more pressing. Here we present SwarmTCR, a method that uses labeled sequence data to predict the specificity of T cell receptors using a nearest-neighbor approach. SwarmTCR works by optimizing the weights of the individual CDR regions to maximize classification performance. Results We compared the performance of SwarmTCR against another nearest-neighbor method and showed that SwarmTCR performs well both with bulk sequencing data and with single cell data. In addition, we show that the weights returned by SwarmTCR are biologically interpretable. Conclusions Computationally predicting the specificity of T cell receptors can be a powerful tool to shed light on the immune response against infectious diseases and cancers, autoimmunity, cancer immunotherapy, and immunopathology. SwarmTCR is distributed freely under the terms of the GPL-3 license. The source code and all sequencing data are available at GitHub (https://github.com/thecodingdoc/SwarmTCR).

Motivation: Computationally predicting the specificity of T cell receptors can be a powerful tool to shed light on the immune response against infectious diseases and cancers, autoimmunity, cancer immunotherapy, and immunopathology. With more T cell receptor sequence data becoming available, the need for bioinformatics approaches to tackle this problem is even more pressing. Here we present SwarmTCR, a method that uses labeled sequence data to predict the specificity of T cell receptors using a nearest-neighbor approach. SwarmTCR works by optimizing the weights of the individual CDR regions to maximize classification performance. Results: We compared the performance of SwarmTCR against a state-of-the-art method (TCRdist) and showed that SwarmTCR performed significantly better on epitopes EBV-BRLF1, NS4B with single cell data and epitopes EBV-BRLF1, IAV-M1 with bulk sequencing data (alpha and beta chains). In addition, we show that the weights returned by SwarmTCR are biologically interpretable. Availability: SwarmTCR is distributed freely under the terms of the GPL-3 license. The source code and all sequencing data are available at GitHub (https://github.com/thecodingdoc/SwarmTCR).