A comprehensive analysis of the germline and expressed TCR repertoire in White Peking duck

Macroevolution of avian T cell receptor C segments using genomic data

All jawed vertebrates have four T cell receptor (TCR) chains expressed by thymus-derived lymphocytes that play a significant role in animal immune defense. However, avian TCR studies have been limited to a few species, although their co-functional major histocompatibility complexes (MHCs) have been studied for decades, showing various copy numbers and polymorphisms. Here, using public genome data, we characterized the copy numbers, the phylogenic relationship and selection of T cell receptor complex (TCR-C) segments, and the genomic organization of TCR loci across birds. Various numbers of C segments were found in the TCRα/TCRδ, TCRβ, and TCRγ loci, and phylogenetic analysis reflected both ancient gene duplication events (two Cβ segments and Cδ segments divergent into CδI and CδII) and contemporary evolution (lineage-specific and species-specific characteristics). Most passerines lack CδII segments and a second TRD locus, except Meliphagidae and Maluridae. A relatively stable structure was verified in four TCR loci of birds, except for the arrangement of V segment groups. In this study, we explored the phylogenetic relationships of TCR-C segments across avians for the first time. We inferred gene duplication and loss events during the evolution process. The finding of diverse TCR germline repertoires provides a better understanding of the immune systems of birds.

Repertoire analysis of γδ T cells in the chicken enables functional annotation of the genomic region revealing highly variable pan-tissue TCR gamma V gene usage as well as identifying public and private repertoires

BACKGROUND

Despite increasing interest in γδ T cells and their non-classical behaviour, most studies focus on animals with low numbers of circulating γδ T cells, such as mice and humans. Arguably, γδ T cell functions might be more prominent in chickens where these cells form a higher proportion of the circulatory T cell compartment. The TCR repertoire defines different subsets of γδ T cells, and such analysis is facilitated by well-annotated TCR loci. γδ T cells are considered at the cusp of innate and adaptive immunity but most functions have been identified in γδ low species. A deeper understanding of TCR repertoire biology in γδ high and γδ low animals is critical for defining the evolution of the function of γδ T cells. Repertoire dynamics will reveal populations that can be classified as innate-like or adaptive-like as well as those that straddle this definition.

RESULTS

Here, a recent discrepancy in the structure of the chicken TCR gamma locus is resolved, demonstrating that tandem duplication events have shaped the evolution of this locus. Importantly, repertoire sequencing revealed large differences in the usage of individual TRGV genes, a pattern conserved across multiple tissues, including thymus, spleen and the gut. A single TRGV gene, TRGV3.3, with a highly diverse private CDR3 repertoire dominated every tissue in all birds. TRGV usage patterns were partly explained by the TRGV-associated recombination signal sequences. Public CDR3 clonotypes represented varying proportions of the repertoire of TCRs utilising different TRGVs, with one TRGV dominated by super-public clones present in all birds.

CONCLUSIONS

The application of repertoire analysis enabled functional annotation of the TCRG locus in a species with a high circulating γδ phenotype. This revealed variable usage of TCRGV genes across multiple tissues, a pattern quite different to that found in γδ low species (human and mouse). Defining the repertoire biology of avian γδ T cells will be key to understanding the evolution and functional diversity of these enigmatic lymphocytes in an animal that is numerically more reliant on them. Practically, this will reveal novel ways in which these cells can be exploited to improve health in medical and veterinary contexts.

Characterization of the chicken T cell receptor γ repertoire by high-throughput sequencing

BACKGROUND

As one of "γδ-high" species, chicken is an excellent model for the study of γδ T cells in non-mammalian animals. However, a comprehensive characterization of the TCRγδ repertoire is still missing in chicken. The objective of this study was to characterize the expressed TCRγ repertoire in chicken thymus using high-throughput sequencing.

METHODS

In this study, we first obtained the detailed genomic organization of the TCRγ locus of chicken based on the latest assembly of the red jungle fowl genome sequences (GRCg6a) and then characterized the TCRγ repertoire in the thymus of four chickens by using 5' Rapid Amplification of cDNA Ends (5' RACE) along with high-throughput sequencing (HTS).

RESULTS

The chicken TCRγ locus contains a single Cγ gene, three functional Jγ segments and 44 Vγ segments that could be classified into six subgroups, each containing six, nineteen, nine, four, three and three members. Dot-plot analysis of the chicken TCRγ locus against itself showed that almost all the entire zone containing Vγ segments had arisen through tandem duplication events, and the main homology unit, containing 9 or 10 Vγ gene segments, has tandemly duplicated for four times. For the analysis of chicken TCRγ repertoire, more than 100,000 unique Vγ-region nucleotide sequences were obtained from the thymus of each chicken. After alignment to the germline Vγ and Jγ segments identified above, we found that the four chickens had similar repertoire profile of TCRγ. In brief, four Vγ segments (including Vγ3.7, Vγ2.13, Vγ1.6 and Vγ1.3) and six Vγ-Jγ pairs (including Vγ3.7-Jγ3, Vγ2.13-Jγ1, Vγ2.13-Jγ3, Vγ1.6-Jγ3, Vγ3.7-Jγ1 and Vγ1.6-Jγ1) were preferentially utilized by all four individuals, and vast majority of the unique CDR3γ sequences encoded 4 to 22 amino acids with mean 12.90 amino acids, which exhibits a wider length distribution and/or a longer mean length than CDR3γ of human, mice and other animal species.

CONCLUSIONS

In this study, we present the first in-depth characterization of the TCRγ repertoire in chicken thymus. We believe that these data will facilitate the studies of adaptive immunology in birds.

Characterization of the domestic goat γδ T cell receptor gene loci and gene usage

Goats and cattle diverged 30 million years ago but retain similarities in immune system genes. Here, the caprine T cell receptor (TCR) gene loci and transcription of its genes were examined and compared to cattle. We annotated the TCR loci using an improved genome assembly (ARS1) of a highly homozygous San Clemente goat. This assembly has already proven useful for describing other immune system genes including antibody and leucocyte receptors. Both the TCRγ (TRG) and TCRδ (TRD) loci were similarly organized in goats as in cattle and the gene sequences were highly conserved. However, the number of genes varied slightly as a result of duplications and differences occurred in mutations resulting in pseudogenes. WC1⁺ γδ T cells in cattle have been shown to use TCRγ genes from only one of the six available cassettes. The structure of that Cγ gene product is unique and may be necessary to interact with WC1 for signal transduction following antigen ligation. Using RT-PCR and PacBio sequencing, we observed the same restriction for goat WC1⁺ γδ T cells. In contrast, caprine WC1⁺ and WC1^- γδ T cell populations had a diverse TCRδ gene usage although the propensity for particular gene usage differed between the two cell populations. Noncanonical recombination signal sequences (RSS) largely correlated with restricted expression of TCRγ and δ genes. Finally, caprine γδ T cells were found to incorporate multiple TRD diversity gene sequences in a single transcript, an unusual feature among mammals but also previously observed in cattle.

Analysis of the Chinese Alligator TCRα/δ Loci Reveals the Evolutionary Pattern of Atypical TCRδ/TCRμ in Tetrapods

Atypical TCRδ found in sharks, amphibians, birds, and monotremes and TCRμ found in monotremes and marsupials are TCR chains that use Ig or BCR-like variable domains (VHδ/Vμ) rather than conventional TCR V domains. These unconventional TCR are consistent with a scenario in which TCR and BCR, although having diverged from each other more than 400 million years ago, continue to exchange variable gene segments in generating diversity for Ag recognition. However, the process underlying this exchange and leading to the evolution of these atypical TCR receptor genes remains elusive. In this study, we identified two TCRα/δ gene loci in the Chinese alligator (Alligator sinensis). In total, there were 144 V, 154 Jα, nine Jδ, eight Dδ, two Cα, and five Cδ gene segments in the TCRα/δ loci of the Chinese alligator, representing the most complicated TCRα/δ gene system in both genomic structure and gene content in any tetrapod examined so far. A pool of 32 VHδ genes divided into 18 subfamilies was found to be scattered over the two loci. Phylogenetic analyses revealed that these VHδ genes could be related to bird VHδ genes, VHδ/Vμ genes in platypus or opossum, or alligator VH genes. Based on these findings, a model explaining the evolutionary pattern of atypical TCRδ/TCRμ genes in tetrapods is proposed. This study sheds new light on the evolution of TCR and BCR genes, two of the most essential components of adaptive immunity.

Evolution of the T-Cell Receptor (TR) Loci in the Adaptive Immune Response: The Tale of the TRG Locus in Mammals

T lymphocytes are the principal actors of vertebrates' cell-mediated immunity. Like B cells, they can recognize an unlimited number of foreign molecules through their antigen-specific heterodimer receptors (TRs), which consist of αβ or γδ chains. The diversity of the TRs is mainly due to the unique organization of the genes encoding the α, β, γ, and δ chains. For each chain, multi-gene families are arranged in a TR locus, and their expression is guaranteed by the somatic recombination process. A great plasticity of the gene organization within the TR loci exists among species. Marked structural differences affect the TR γ (TRG) locus. The recent sequencing of multiple whole genome provides an opportunity to examine the TR gene repertoire in a systematic and consistent fashion. In this review, we report the most recent findings on the genomic organization of TRG loci in mammalian species in order to show differences and similarities. The comparison revealed remarkable diversification of both the genomic organization and gene repertoire across species, but also unexpected evolutionary conservation, which highlights the important role of the T cells in the immune response.

Genomic organization of the chicken TCRβ locus originated by duplication of a Vβ segment combined with a trypsinogen gene

Based on the latest assembly of the red jungle fowl (Gallus gallus) genome sequence, we characterized the detailed genomic organization of the T cell receptor beta (TCRβ) locus of chicken. The chicken TCRβ locus spans approximately 210 kb, and is organized in a typical translocon organization as previously reported. Within this locus, a total of 16 germline Vβ gene segments were classified into three subgroups, containing 11, four, and one members, respectively. Phylogenetic analysis revealed that the chicken Vβ3.1 segment was homologous with the duck Vβ1 subgroup, and further clustered with Vβ  segments from reptiles but not amphibians. We also identified nine protease serine 1 (PRSS1) and three protease serine 2 (PRSS2) genes, which were interspersed within the chicken TCRβ locus. Dot-plot analysis of the chicken TCRβ locus against itself revealed that the 5' part of the locus had arisen through a series of tandem duplication events. The homology units were composed of one Vβ1 segment followed by a PRSS1 gene, or one Vβ2 segment followed by a PRSS2 gene. This duplication pattern, in which the Vβ segments and trypsinogen genes form a duplication unit, was unique to TCRβ loci of chicken and duck, but not observed in TCRβ loci of other tetrapods studied thus far. By analyzing the cloned TCRβ cDNA sequences, we found that the usage pattern of Vβ segments was consistent with the results of previous studies. These studies showed that members of the Vβ1 subgroup are preferentially utilized in V-D-J recombination. Furthermore, we found that the Vβ3.1 segment participated into V-D-J recombination, but at a very low frequency. The length distribution of the chicken complementarity-determining region 3β (CDR3β) showed a tendency similar to that observed for the duck.

Identification of Two Nonrearranging IgSF Genes in Chicken Reveals a Novel Family of Putative Remnants of an Antigen Receptor Precursor

In this study, we identified a pair of nonrearranging VJ-joined Ig superfamily genes, termed putative remnants of an Ag receptor precursor (PRARP) genes, in chicken. Both genes encode a single V-set Ig domain consisting of a canonical J-like segment and a potential immunoreceptor tyrosine-based inhibitory or switch motif in the cytoplasmic region. In vitro experiments showed that both genes were expressed at the cell surface as membrane proteins, and their recombinant products formed a monomer and a disulfide-linked homodimer or a heterodimer. These two genes were mainly expressed in B and T cells and were upregulated in response to stimulation with poly(I:C) in vitro and vaccination in vivo. Orthologs of PRARP have been identified in bony fish, amphibians, reptiles, and other birds, and a V-C1 structure similar to that of Ig or TCR chains was found in all these genes, with the exception of those in avian species, which appear to contain degenerated C1 domains or divergent Ig domains. Phylogenetic analyses suggested that the newly discovered genes do not belong to any known immune receptor family and appear to be a novel gene family. Further elucidation of the functions of PRARP and their origin might provide significant insights into the evolution of the immune system of jawed vertebrates.

The Florida manatee (Trichechus manatus latirostris) T cell receptor loci exhibit V subgroup synteny and chain-specific evolution

The Florida manatee (Trichechus manatus latirostris) has limited diversity in the immunoglobulin heavy chain. We therefore investigated the antigen receptor loci of the other arm of the adaptive immune system: the T cell receptor. Manatees are the first species from Afrotheria, a basal eutherian superorder, to have an in-depth characterization of all T cell receptor loci. By annotating the genome and expressed transcripts, we found that each chain has distinct features that correlates to their individual functions. The genomic organization also plays a role in modulating sequence conservation between species. There were extensive V subgroup synteny blocks in the TRA and TRB loci between T. m. latirostris and human. Increased genomic locus complexity correlated to increased locus synteny. We also identified evidence for a VHD pseudogene for the first time in a eutherian mammal. These findings emphasize the value of including species within this basal eutherian radiation in comparative studies.

Analysis of TCRβ and TCRγ genes in Chinese alligator provides insights into the evolution of TCR genes in jawed vertebrates

All jawed vertebrates have four T cell receptor (TCR) chains that are expressed by thymus-derived lymphocytes and play a major role in animal immune defence. However, few studies have investigated the TCR chains of crocodilians compared with those of birds and mammals, despite their key evolutionary position linking amphibians, reptiles, birds and mammals. Here, employing an Alligator sinensis genomic bacterial artificial chromosome (BAC) library and available genome data, we characterized the genomic organization, evolution and expression of TRB and TRG loci in Alligator sinensis. According to the sequencing data, the Alligator sinensis TRB locus spans approximately 500 Kb of genomic DNA containing two D-J-C clusters and 43 V gene segments and is organized as Vβ₍₃₉₎-pJβ1-pCβ1-pDβ1-Dβ2- Jβ2₍₁₂₎-Cβ2-Vβ₍₄₎, whereas the TRG locus spans 115 Kb of DNA genomic sequence consisting of 18 V gene segments, nine J gene segments and one C gene segment and is organized in a classical translocon pattern as Vγ₍₁₈₎-Jγ₍₉₎-Cγ. Moreover, syntenic analysis of TRB and TRG chain loci suggested a high degree of conserved synteny in the genomic regions across mammals, birds and Alligator sinensis. By analysing the cloned TRB/TRG cDNA, we identified the usage pattern of V families in the expressed TRB and TRG. An analysis of the junctions of the recombined VJ revealed the presence of N and P nucleotides in both expressed TRB and TRG sequences. Phylogenetic analysis revealed that TRB and TRG loci possess distinct evolutionary patterns. Most Alligator sinensis V subgroups have closely related orthologues in chicken and duck, and a small number of Alligator sinensis V subgroups have orthologues in mammals, which supports the hypothesis that crocodiles are the closest relatives of birds and mammals. Collectively, these data provide insights into TCR gene evolution in vertebrates and improve our understanding of the Alligator sinensis immune system.