A series of products were amplified using a PCR strategy based on short minimally degenerate primers and R. eglanteria (clearnose skate) spleen cDNA as template. These products were used as probes to select corresponding cDNAs from a spleen cDNA library. The cDNA sequences exhibit significant identity with prototypic (alpha, beta, gamma, and delta T cell antigen receptor (TCR) genes. Characterization of cDNAs reveals extensive variable region diversity, putative diversity segments, and varying degrees of junctional diversification. This demonstrates expression of both alpha/beta and gamma/delta TCR genes at an early level of vertebrate phylogeny and indicates that the three major known classes of rearranging antigen receptors were present in the common ancestor of the present-day jawed vertebrates.

Alphabeta T cells and gammadelta T cells perform nonoverlapping immune functions. In mammalian species with a high percentage of very diverse gammadelta T cells, like ruminants and pigs, it is often assumed that alphabeta T cells are less diverse than gammadelta T cells. Based on the bovine genome, we have created a map of the bovine TRA/TRD locus and show that, in cattle, in addition to the anticipated >100 TRDV genes, there are also >300 TRAV or TRAV/DV genes. Among the V genes in the TRA/TRD locus, there are several genes that lack a CDR2 and are functionally rearranged and transcribed and, in some cases, have an extended CDR1. The number of bovine V genes is a multiple of the number in mice and humans and may encode T cell receptors that use a novel way of interacting with antigen.

Previously, somatic hypermutation (SHM) was considered to be exclusively associated with affinity maturation of antibodies, although it also occurred in T cells under certain conditions. More recently, it has been shown that SHM generates diversity in the variable domain of T cell receptor (TCR) in camel and shark. Here, we report somatic mutations in TCR alpha chain genes of the teleost fish, Ballan wrasse (Labrus bergylta), and show that this mechanism adds extra diversity to the polymorphic constant (C) region as well. The organization of the TCR alpha/delta locus in Ballan wrasse was obtained from a scaffold covering a single copy C alpha gene, 65 putative J alpha segments, a single copy C delta gene, 1 J delta segment, and 2 D delta segments. Analysis of 37 fish revealed 6 allotypes of the C alpha gene, each with 1-3 replacement substitutions. Somatic mutations were analyzed by molecular cloning of TCR alpha chain cDNA. Initially, 79 unique clones comprising four families of variable (V) alpha genes were characterized. Subsequently, a more restricted PCR was performed to focus on a specific V gene. Comparison of 48 clones indicated that the frequency of somatic mutations in the VJ region was 4.5/1,000 base pairs (bps), and most prevalent in complementary determining region 2 (CDR2). In total, 45 different J segments were identified among the 127 cDNA clones, counting for most of the CDR3 diversity. The number of mutations in the C alpha chain gene was 1.76 mutations/1,000 bps and A nucleotides were most frequently targeted, in contrast to the VJ region, where G nucleotides appeared to be mutational hotspots. The replacement/synonymous ratios in the VJ and C regions were 2.5 and 1.85, respectively. Only 7% of the mutations were found to be linked to the activation-induced cytidine deaminase hotspot motif (RGYW/WRCY).

In jawed vertebrates, two major T cell populations have been characterized. They are defined as α/β or γ/δ T cells, based on the expressed T cell receptor. Salmonids (family Salmonidae) include two key teleost species for aquaculture, rainbow trout (Oncorhynchus mykiss) and Atlantic salmon (Salmo salar) which constitute important models for fish immunology and important targets for vaccine development. The growing interest to decipher the dynamics of adaptive immune responses against pathogens or vaccines has resulted in recent efforts to sequence the immunoglobulin (IG) or antibodies and T cell receptor (TR) repertoire in these species. In this context, establishing a comprehensive and coherent locus annotation is the fundamental basis for the analysis of high-throughput repertoire sequencing data. We therefore decided to revisit the description and annotation of TRA/TRD locus in Atlantic salmon and two strains of rainbow trout (Swanson and Arlee) using the now available high-quality genome assemblies. Phylogenetic analysis of functional TRA/TRD V genes from these three genomes led to the definition of 25 subgroups shared by both species, some with particular feature. A total of 128 TRAJ genes were identified in Salmo, the majority with a close counterpart in Oncorhynchus. Analysis of expressed TRA repertoire indicates that most TRAV gene subgroups are expressed at mucosal and systemic level. The present work on TRA/TRD locus annotation along with the analysis of TRA repertoire sequencing data show the feasibility and advantages of a common salmonid TRA/TRD nomenclature that allows an accurate annotation and analysis of high-throughput sequencing results, across salmonid T cell subsets.

Predicting the binding specificity of T Cell Receptors (TCR) to MHC-peptide complexes (pMHCs) is essential for the development of repertoire-based biomarkers. This affinity may be affected by different components of the TCR, the peptide, and the MHC allele. Historically, the main element used in TCR-peptide binding prediction was the Complementarity Determining Region 3 (CDR3) of the beta chain. However, recently the contribution of other components, such as the alpha chain and the other V gene CDRs has been suggested. We use a highly accurate novel deep learning-based TCR-peptide binding predictor to assess the contribution of each component to the binding.

In mammals, T cells develop along two discrete pathways characterized by expression of either the αβ or the γδ T cell receptors. Human and mouse display a low peripheral blood γδ T cell percentage ("γδ low species") while sheep, bovine and pig accounts for a high proportion of γδ T lymphocytes ("γδ high species"). While the T cell receptor alpha (TRA) and delta (TRD) genes and the genomic organization of the TRA/TRD locus has been determined in human and mouse, this information is still poorly known in artiodactyl species, such as sheep.

Organ-specific autoimmune disease can be caused by alphabeta T cells that have escaped self-tolerance induction. Here we show that one of the causes of escape from self-tolerance is the coexpression of two different T cell receptors by the same cell, which can occur in up to 30% of all T cells in normal mice and can lead to low-level surface expression of an autospecific TCR. We found that double receptor-expressing T cells can escape tolerance even to ubiquitously expressed antigens but can nevertheless induce autoimmune diabetes when the relevant protein is expressed in pancreatic tissue. Such diabetogenic T cells are absent, however, among T cells expressing the autospecific TCR as the sole receptor.

Murine intestinal intraepithelial lymphocytes (iIELs) are made up of a heterogeneous mix of T cells with unique phenotypes. Whereas CD8(+) T cells in peripheral lymphoid organs use CD8alpha/beta and are selected on MHC class Ia molecules, a majority of iIELs use CD8alpha/alpha. Here, we report that the presence of CD8alpha/alpha TCR-alpha/beta cells in iIELs is independent of classical MHC class I molecules K(b) and D(b), as illustrated by their presence in K(b)/D(b) double-knockout mice and in mice lacking a nonclassical MHC class I molecule, CD1d. Most strikingly, their presence is decreased by approximately 70% in mice lacking transporter associated with antigen processing (TAP). The TAP-dependent nonclassical MHC class I molecule Qa-2 is strongly implicated in the presence of these cells, as inferred from the low numbers of CD8alpha/alpha TCR-alpha/beta T cells in mice deficient in Qa-2 genes. Second, a Qa-2-transgenic mouse made in a Qa-2(-) strain showed an increase in the numbers of CD8alpha/alpha cells among its iIELs. Thus, the presence of CD8alpha/alpha TCR-alpha/beta cells in iIELs is mainly dependent on the nonclassical MHC class I molecule Qa-2.

The complete TCR alpha/delta locus of Atlantic salmon (Salmo salar) has been characterized and annotated. In the 900 kb TCR alpha/delta locus, 292 Valpha/delta segments and 123 Jalpha/delta segments were identified. Of these, 128 Valpha/delta, 113 Jalpha, and a Jdelta segment appeared to be functional as they lacked frame shifts or stop codons. This represents the largest repertoire of Valpha/delta and Jalpha segments of any organism to date. The 128 functional Valpha/delta segments could be grouped into 29 subgroups based upon 70% nucleotide similarity. Expression data confirmed the usage of the diverse repertoire found at the genomic level. At least 99 Valpha, 13 Vdelta 86 Jalpha, 1 Jdelta, and 2 Ddelta segments were used in TCR alpha or delta transcription, and 652 unique genes were identified from a sample of 759 TCRalpha cDNA clones. Cumulatively, the genomic and expression data suggest that the Atlantic salmon T-cell receptor has enormous capacity to recognize a wide diversity of antigens.

Metastatic medullary thyroid cancer (MTC) is a rare but often aggressive thyroid malignancy with a 5-year survival rate of less than 40% and few effective therapeutic options. Adoptive T cell immunotherapy using chimeric antigen receptor (CAR)-modified T cells (CAR Ts) is showing encouraging results in the treatment of cancer, but development is challenged by the availability of suitable target antigens. We identified glial-derived neurotrophic factor (GDNF) family receptor alpha 4 (GFRα4) as a putative antigen target for CAR-based therapy of MTC. We show that GFRα4 is highly expressed in MTC, in parafollicular cells within the thyroid from which MTC originates, and in normal thymus. We isolated two single-chain variable fragments (scFvs) targeting GFRα4 isoforms a and b by antibody phage display. CARs bearing the CD3ζ and the CD137 costimulatory domains were constructed using these GFRα4-specific scFvs. GFRα4-specific CAR Ts trigger antigen-dependent cytotoxicity and cytokine production in vitro, and they are able to eliminate tumors derived from the MTC TT cell line in an immunodeficient mouse xenograft model of MTC. These data demonstrate the feasibility of targeting GFRα4 by CAR T and support this antigen as a promising target for adoptive T cell immunotherapy and other antibody-based therapies for MTC.

The generation of a productive "in-frame" T cell receptor beta (TCR beta), immunoglobulin (Ig) heavy (H) or Ig light (L) chain variable region gene can result in the cessation of rearrangement of the alternate allele, a process referred to as allelic exclusion. This process ensures that most alphabeta T cells express a single TCR beta chain and most B cells express single IgH and IgL chains. Assembly of TCR alpha and TCR gamma chain variable region genes exhibit allelic inclusion and alphabeta and gammadelta T cells can express two TCR alpha or TCR gamma chains, respectively. However, it was not known whether assembly of TCR delta variable regions genes is regulated in the context of allelic exclusion. To address this issue, we have analyzed TCR delta rearrangements in a panel of mouse splenic gammadelta T cell hybridomas. We find that, similar to TCR alpha and gamma variable region genes, assembly of TCR delta variable region genes exhibits properties of allelic inclusion. These findings are discussed in the context of gammadelta T cell development and regulation of rearrangement of TCR delta genes.

T cells specific for nucleosomal autoepitopes are selectively expanded in lupus mice and these Th cells drive autoimmune B cells to produce pathogenic antinuclear antibodies. We transfected the TCR-alpha and -beta chain genes of a representative, pathogenic autoantibody-inducing Th clone specific for the nucleosomal core histone peptide H471-94 into TCR-negative recipient cells. Although the autoimmune TCRs were originally derived from SNF1 (I-Ad/q) mice, the transfectants could recognize the nucleosomal autoepitope presented by APC-bearing I-A molecules of all haplotypes tested, as well as human DR molecules. Competition assays indicated that the autoepitopes bound to the MHC class II groove. Most remarkably, MHC-unrestricted recognition of the nucleosomal peptide epitope was conferred by the lupus TCR-alpha chain even when it paired with a TCR-beta chain of irrelevant specificity. Several other disease-relevant Th clones and splenic T cells of lupus mice had similar properties. The TCR-alpha chains of these murine lupus Th clones shared related motifs and charged residues in their CDRs, and similar motifs were apparent even in TCR-alpha chains of human lupus Th clones. The lupus TCR-alpha chains probably contact the nucleosomal peptide complexed with MHC with relatively high affinity/avidity to sustain TCR signaling, because CD4 coreceptor was not required for promiscuous recognition. Indeed, pathogenic autoantibody-inducing, CD4-negative, TCR-alphabeta+ Th cells are expanded in systemic lupus erythematosus. These results have implications regarding thymic selection and peripheral expansion of nucleosome-specific T cells in lupus. They also suggest that universally tolerogenic epitopes could be designed for therapy of lupus patients with diverse HLA alleles. We propose to designate nucleosomes and other antigens bearing universal epitopes "Pantigens" (for promiscuous antigens).

Knowledge of the complete nucleotide sequence of the mouse TCRAD locus allows an accurate determination V-J rearrangement status. Using multiplex genomic PCR assays and real time PCR analysis, we report a comprehensive and systematic analysis of the V-J recombination of TCR alpha chain in normal mouse thymocytes during development. These respective qualitative and quantitative approaches give rise to four major points describing the control of gene rearrangements. (a) The V-J recombination pattern is not random during ontogeny and generates a limited TCR alpha repertoire; (b) V-J rearrangement control is intrinsic to the thymus; (c) each V gene rearranges to a set of contiguous J segments with a gaussian-like frequency; (d) there are more rearrangements involving V genes at the 3' side than 5' end of V region. Taken together, this reflects a preferential association of V and J gene segments according to their respective positions in the locus, indicating that accessibility of both V and J regions is coordinately regulated, but in different ways. These results provide a new insight into TCR alpha repertoire size and suggest a scenario for V usage during differentiation.

Peripheral T-cell lymphomas (PTCLs) are aggressive non-Hodgkin lymphomas with generally poor outcomes following standard therapy. Few candidate therapeutic targets have been identified to date. Retinoic acid receptor alpha (RARA) is a transcription factor that modulates cell growth and differentiation in response to retinoids. While retinoids have been used to treat some cutaneous T-cell lymphomas (CTCLs), their mechanism of action and the role of RARA in CTCL and other mature T-cell lymphomas remain poorly understood. After identifying a PTCL with a RARAR394Q mutation, we sought to characterize the role of RARA in T-cell lymphoma cells. Overexpressing wild-type RARA or RARAR394Q significantly increased cell growth in RARAlow cell lines, while RARA knockdown induced G1 arrest and decreased expression of cyclin-dependent kinases CDK2/4/6 in RARAhigh cells. The retinoids, AM80 (tamibarotene) and all-trans retinoic acid, caused dose-dependent growth inhibition, G1 arrest, and CDK2/4/6 down-regulation. Genes down-regulated in transcriptome data were enriched for cell cycle and G1-S transition. Finally, RARA overexpression augmented chemosensitivity to retinoids. In conclusion, RARA drives cyclin-dependent kinase expression, G1-S transition, and cell growth in T-cell lymphoma. Synthetic retinoids inhibit these functions in a dose-dependent fashion and are most effective in cells with high RARA expression, indicating RARA may represent a therapeutic target in some PTCLs.

The complete germline repertoires of the channel catfish, Ictalurus punctatus, T cell receptor (TR) loci, TRAD, TRB, and TRG were obtained by analyzing genomic data from PacBio sequencing. The catfish TRB locus spans 214 kb, and contains 112 TRBV genes, a single TRBD gene, 31 TRBJ genes and two TRBC genes. In contrast, the TRAD locus is very large, at 1,285 kb. It consists of four TRDD genes, one TRDJ gene followed by the exons for TRDC, 125 TRAJ genes and the exons encoding the TRAC. Downstream of the TRAC, are 140 TRADV genes, and all of them are in the opposite transcriptional orientation. The catfish TRGC locus spans 151 kb and consists of four diverse V-J-C cassettes. Altogether, this locus contains 15 TRGV genes and 10 TRGJ genes. To place our data into context, we also analyzed the zebrafish TR germline gene repertoires. Overall, our findings demonstrated that catfish possesses a more restricted repertoire compared to the zebrafish. For example, the 140 TRADV genes in catfish form eight subgroups based on members sharing 75% nucleotide identity. However, the 149 TRAD genes in zebrafish form 53 subgroups. This difference in subgroup numbers between catfish and zebrafish is best explained by expansions of catfish TRADV subgroups, which likely occurred through multiple, relatively recent gene duplications. Similarly, 112 catfish TRBV genes form 30 subgroups, while the 51 zebrafish TRBV genes are placed into 36 subgroups. Notably, several catfish and zebrafish TRB subgroups share ancestor nodes. In addition, the complete catfish TR gene annotation was used to compile a TR gene segment database, which was applied in clonotype analysis of an available gynogenetic channel catfish transcriptome. Combined, the TR annotation and clonotype analysis suggested that the expressed TRA, TRB, and TRD repertoires were generated by different mechanisms. The diversity of the TRB repertoire depends on the number of TRBV subgroups and TRBJ genes, while TRA diversity relies on the many different TRAJ genes, which appear to be only minimally trimmed. In contrast, TRD diversity relies on nucleotide additions and the utilization of up to four TRDD segments.

Mice lacking the interleukin 7 receptor (IL-7R) generate alpha/beta T cells at a detectable but greatly reduced rate, but gamma/delta T cells are completely absent. The special role of IL-7R signaling in gamma/delta T cell development has remained unclear. IL-7Ralpha(-/-) mice exhibit a paucity of gamma gene rearrangements. This striking observation can be explained by a defect in T cell receptor (TCR)-gamma gene rearrangement, a defect in TCR-gamma gene transcription leading to death of gamma/delta lineage cells, and/or a requirement for IL-7R in commitment of cells to the gamma/delta lineage. To determine the role of IL-7R signaling in gamma/delta T cell development, we examined transcription of a prerearranged TCR-gamma transgene in IL-7Ralpha(-/-) mice, as well as the effects of IL-7 on transcription of endogenous, rearranged TCR-gamma genes in alpha/beta lineage cells. The results demonstrate that IL-7R-mediated signals are necessary for the normal expression of rearranged TCR-gamma genes. Equally significant, the results show that the poor expression of TCR-gamma genes in IL-7Ralpha(-/-) mice is responsible for the selective deficit in gamma/delta cells in these mice, since a high copy TCR-gamma transgene exhibited sufficient residual expression in IL-7Ralpha(-/-) mice to drive gamma/delta cell development. The results indicate that the absence of gamma/delta T cells in IL-7Ralpha(-/-) mice is due to insufficient TCR-gamma gene expression.

Human T cell receptors (TCRs) are critical for mediating immune responses to pathogens and tumors and regulating self-antigen recognition. Yet, variations in the genes encoding TCRs remain insufficiently defined. Detailed analysis of expressed TCR alpha, beta, gamma, and delta genes in 45 donors from four human populations-African, East Asian, South Asian, and European-revealed 175 additional TCR variable and junctional alleles. Most of these contained coding changes and were present at widely differing frequencies in the populations, a finding confirmed using DNA samples from the 1000 Genomes Project. Importantly, we identified three Neanderthal-derived, introgressed TCR regions including a highly divergent TRGV4 variant, which mediated altered butyrophilin-like molecule 3 (BTNL3) ligand reactivity and was frequent in all modern Eurasian population groups. Our results demonstrate remarkable variation in TCR genes in both individuals and populations, providing a strong incentive for including allelic variation in studies of TCR function in human biology.

Transcriptomics can be combined with TRA and TRB clonotype analysis at the single cell level. The aim of this study was to validate this approach on the ICELL8 Single-Cell system and to evaluate its usefulness to analyse clinical paucicellular samples. For this purpose, we carefully selected T cell lines with defined TRA/TRB clonotypes as well as clinical samples enriched for CD3+ T cells that possess a complex TCR repertoire. Low cell numbers of the different samples were dispensed in a chip on the ICELL8 Single-Cell System. Two sequencing libraries were generated from each single cell cDNA preparation, one for the TRA/TRB repertoire and one for the 5' ends of transcripts, and subsequently sequenced. Transcriptome analysis revealed that the cell lines on average express 2,268 unique genes/cell and T cells of clinical samples 770 unique genes/cell. The expected combined TRA/TRB clonotype was determined for on average 71% of the cells of the cell lines. In the clinical samples the TRA/TRB repertoire was more complex than those of the cell lines. Furthermore, the TRB clonotype distribution of the clinical samples was positively correlated to frequencies of TCRVβ families in CD3+ T cells obtained by a flow cytometry-based approach (Spearman's Rho correlation coefficient 0.81, P = 6.49 * 10-7). Combined analyses showed that transcriptome-based cell type-specific clusters in clinical samples corresponded to clinical features such as CMV status. In conclusion, we showed that the ICELL8 Single-Cell System enabled combined interrogation of both TRA/TRB repertoire and transcriptome of paucicellular clinical samples. This opens the way to study the response of single T cells within heterogeneous samples for both their transcriptome and TRA/TRB clonotypes in disease or upon treatment.

Accurate characterization of the repertoire of the T-cell receptor (TCR) alpha and beta chains is critical to understanding adaptive immunity. Such characterization has many applications across such fields as vaccine development and response, clone-tracking in cancer, and immunotherapy. Here we present a new methodology called single-cell TCRseq (scTCRseq) for the identification and assembly of full-length rearranged V(D)J T-cell receptor sequences from paired-end single-cell RNA sequencing reads. The method allows accurate identification of the V(D)J rearrangements for each individual T-cell and has the novel ability to recover paired alpha and beta segments. Source code is available at https://github.com/ElementoLab/scTCRseq .

Expression of a TCRgamma transgene in RAG-1-/- mice resulted in the development of a limited number of CD4+CD8+ (DP) thymocytes. In vivo treatments with anti-TCRgamma antibody enhanced the number of DP thymocytes, demonstrating that TCRgamma chains were expressed on the cell surface in the absence of delta, alpha, or beta chains. Mutations in pTalpha or CD3epsilon genes abolished transgene-induced DP cell development, indicating that TCRgamma can associate with pTalpha and CD3 to form a novel pre-TCR. With a transgene containing additional regulatory sequences, TCRgamma expression was down-regulated in DP cells, and little DP cell development occurred. Thus, the function of the endogenous TCRgamma/pTalpha is limited by the transcriptional down-regulation of TCRgamma genes that normally accompanies DP cell development.