# Annotation and classification of the bovine T cell receptor delta genes

C. Herzig, Marie-Paule Lefranc, C. Baldwin

<figure><img src="image_1.png" /></figure>

## To cite this version:

C. Herzig, Marie-Paule Lefranc, C. Baldwin. Annotation and classification of the bovine T cell receptor delta genes. BMC Genomics, 2010, 11 (1), pp.100. 10.1186/1471-2164-11-100. hal-00456810

HAL Id: hal-00456810

https://hal.science/hal-00456810v1

Submitted on 26 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Annotation and classification of the bovine T cell receptor delta genes

Carolyn TA Herzig<sup>1,4*</sup>, Marie-Paule Lefranc<sup>2</sup>, Cynthia L Baldwin<sup>1,3</sup>

## Abstract

**Background:** γδ T cells differ from αβ T cells with regard to the types of antigen with which their T cell receptors interact; γδ T cell antigens are not necessarily peptides nor are they presented on MHC. Cattle are considered a "γδ T cell high" species indicating they have an increased proportion of γδ T cells in circulation relative to that in "γδ T cell low" species such as humans and mice. Prior to the onset of the studies described here, there was limited information regarding the genes that code for the T cell receptor delta chains of this γδ T cell high species.

**Results:** By annotating the bovine (*Bos taurus*) genome Btau_3.1 assembly the presence of 56 distinct T cell receptor delta (TRD) variable (V) genes were found, 52 of which belong to the TRDV1 subgroup and were co-mingled with the T cell receptor alpha variable (TRAV) genes. In addition, two genes belonging to the TRDV2 subgroup and single TRDV3 and TRDV4 genes were found. We confirmed the presence of five diversity (D) genes, three junctional (J) genes and a single constant (C) gene and describe the organization of the TRD locus. The TRDV4 gene is found downstream of the C gene and in an inverted orientation of transcription, consistent with its orthologs in humans and mice. cDNA evidence was assessed to validate expression of the variable genes and showed that one to five D genes could be incorporated into a single transcript. Finally, we grouped the bovine and ovine TRDV1 genes into sets based on their relatedness.

**Conclusions:** The bovine genome contains a large and diverse repertoire of TRD genes when compared to the genomes of "γδ T cell low" species. This suggests that in cattle γδ T cells play a more important role in immune function since they would be predicted to bind a greater variety of antigens.

## Background

T lymphocytes can be subdivided into at least two types based on the expression of either the αβ or γδ T cell receptor. Although both perceive antigen, they differ in the types of antigens with which they react. αβ T cell receptors react with antigenic protein peptides in the context of self major histocompatibility complex (MHC) proteins while γδ T cell receptors may react with proteins but this does not involve MHC presentation. They also react with autologous molecules on cells [1-6] as well as nonproteinaceous molecules [7]. The gene repertoires that code for the γδ T cell receptor chains and the T cell receptor gamma (TRG) and delta (TRD) gene locus organizations have been extensively described for humans and mice but to a lesser extent for the artiodac-tyls which includes ruminants and swine. These latter

are "γδ T cell high" species because of their high levels of γδ T cells in circulation ("γδ T cell low" species exhibit much lower levels of γδ T cells in circulation). It is clear that while both αβ and γδ T cells have large and diverse T cell receptor gene repertoires [8-14], "γδ T cell high" species have a TRD gene repertoire that is much more extensive than that in the "γδ T cell low" species [15-21]. The bovine (*Bos taurus*) locus organization and TRD gene repertoire are the subject of this work.

T cell receptor delta and beta chains are encoded by the rearrangement of variable (V), joining (J) and diversity (D) genes making them more complex than the T cell receptor gamma and alpha chains which lack D gene products. In all mammals evaluated the genes encoding the T cell receptor beta and gamma chains are found at the T cell receptor beta (TRB) and TRG loci, respectively. The genes that encode the T cell receptor delta and alpha chains are found at a single chromosomal location with the TRD genes embedded within the

<footnote>* Correspondence: cth2115@columbia.edu</footnote>

<footnote n="1">¹Department of Veterinary and Animal Sciences, University of Massachusetts, Amherst, MA 01003, USA</footnote>
T cell receptor alpha (TRA) locus. For both humans and mice these are located on chromosome 14 and encompass over 1 megabase (Mb) for the combined TRA/TRD locus [13]. The TRD locus embedded within the TRA locus spans approximately 60 kb in humans and 275 kb in mice. The loci comprise TRDV genes (two in humans and five in mice), followed by TRDD genes (three in humans and two in mice), TRDJ genes (four in humans and two in mice), a single TRDC gene and an additional TRDV gene that is located 3' of the TRDC gene in an inverted transcriptional orientation [10,14]. In addition, five and ten functional TRAV/DV genes (that rearrange to either TRDD or TRAJ) have been identified for humans and mice, respectively, and these are found upstream of the embedded TRD locus [10,14].

Limited evidence derived from analysis of cDNA clones from cattle, sheep and swine, as well as limited germline information for sheep, suggests that the general organization of the bovine TRA/TRD locus does not differ greatly from that of humans and mice [17,20]. It also suggests that a much larger repertoire of TRDV genes exists for cattle, as well as for the other "γδ T cell high" species sheep and swine [15,17,20,21]. Indeed, recent mapping of the bovine TRA/TRD locus identified over 100 TRDV and over 300 TRAV or TRAV/DV genes in cattle [22]. For cattle, four genes belonging to three different small subgroups (TRDV2, TRDV3 and TRDV4) have been identified that are orthologous to those same subgroups in sheep [17,23]. In contrast, while the bovine TRDV1 genes have been found to be related to the single TRDV1 gene that occurs in humans, expansion of the TRDV1 subgroup accounts for the larger number of TRDV genes in the "γδ T cell high" species [15,17,20-22]. The swine and sheep TRDV1 subgroups have been estimated to contain at least 31 and 40 genes, respectively [15,21] while the bovine TRDV1 subgroup has been reported to contain at least 37 genes [20]. However, current knowledge about ruminant and swine TRD genes is predominantly based on cDNA evidence rather than on genomic DNA. Because of this, it has been unclear whether the large number of cDNA TRDV1 sequences is the result of multiple genes or polymorphisms among animals; thus, there has been no clear way to classify and name these transcripts. By annotating the bovine genome Btau_3.1 assembly, here we demonstrated the presence of 56 TRDV genes, 52 of which belong to the TRDV1 subgroup. We also proposed TRDV1 sets to classify the bovine and ovine TRDV1 genes based on their phylogenetic grouping.

The existence of multiple D and J genes also contributes to diversity of possible T cell receptor delta chain amino acid sequences since a particular TRDV gene is expected to be able to recombine with any D and J gene. The region of the protein derived from the V-D-J

gene rearrangement is known as the complementarity determining region 3 (CDR3), while the CDR1 and CDR2 loops are germline-encoded by the TRDV gene. The V-D-J gene rearrangement involves recognition of the recombination signal (RS) sequences for subsequent DNA cleavage by the enzymes RAG1 and RAG2 [24]. T cell receptor delta chain sequence diversity is augmented by junctional flexibility and the addition of N and P nucleotides. RS flank the V, D and J genes and are composed of highly conserved heptamer and nonamer sequences separated by either a 12 base pair (bp) spacer (located 5' of TRDD and TRDJ genes) or a 23 bp spacer (located 3' of TRDV and TRDD genes). Since the conserved heptamer and nonamer sequences and spacer lengths have been found to be important for efficient recombination [25] they were also evaluated here. Finally, recent evidence demonstrated that diversity of T cell receptor delta chain sequences is amplified by the occurrence of multiple TRDD genes within a single CDR3 [20,26]. Thus we also evaluated TRDD gene usage by analyzing cDNA sequences and report those findings here.

## Results

### Gene structure and genomic organization of bovine TRD genes

Here we set out to annotate and classify the bovine TRD genes. Prior to the initiation of the studies described here many cDNA sequences representing TRDV1 gene transcripts had been reported [16-20,27] but the actual number of functional bovine TRDV1 genes and their sequences was unknown. As a result it was not possible to differentiate between cDNA sequences derived from distinct genes and those that represented polymorphisms of the same gene among animals. Furthermore, the existence of TRDV2, TRDV3 and TRDV4 genes had been demonstrated [17], along with the TRDJ and TRDD genes [20], but their numbers and placement within the bovine genome had not been definitively described.

The exon-intron structures of TRDV genes (Figure 1A) and of the single TRDC gene (Figure 1B) were determined. A TRDV gene encompasses a coding region of approximately 560 bp and comprises two exons. The first exon, L-PART1, is 49, 55 and 37 bp in length for TRDV1/TRDV2, TRDV3 and TRDV4, respectively. The second exon, V-EXON, is between 271 and 335 bp in length for TRDV1, and is 299, 293 and 309 bp long for TRDV2, TRDV3 and TRDV4, respectively. The TRDC gene encompasses a coding region of 1282 bp comprising three exons. RS sequences were identified, as expected, adjacent to the 3' end of each V-EXON (Figure 1A) while the termination codon was at the 3' end of TRDC exon 3.
A

<figure><img src="image_3.png" /></figure>

B

<figure><img src="image_4.png" /></figure>

**Figure 1** Schematic representation of TRD gene exon-intron structures. (A) Representative structures of TRDV genes (based on genomic sequence for TRDV1αu, GLEAN_15708), which contain two exons, and (B) the single TRDC gene (GLEAN_19705), which contains three exons, are shown. Scale is shown in base pair increments beneath the schematic. V-RS, recombination signal sequence of a V gene.

Table 1 summarizes the proposed functional (or Open Reading Frame [ORF]) TRD genes that were identified within the bovine Btau_3.1 assembly. Most genes were identified within regions that were not placed on a chromosome (i.e. were found on ChrUn). However, it is expected that, like in primates and rodents, these genes are all located within the TRA/TRD locus which in cattle is on chromosome 10. The bovine consensus gene set (known as GLEAN) was used to identify gene prediction models of TRD genes. Where applicable, GLEAN numbers are listed in Table 1, along with the Bovine Genome Scaffold identification number, location and orientation. It should be noted that the orientation of each contig individually, as well as in relation to other contigs, is not definitive at this point for the bovine Btau_3.1 assembly. All putatively functional (or ORF) genes (as determined based on criteria described in Methods) are shown. Based on analysis of flanking and intronic genomic sequence it was determined that many of these gene models might represent the same genes although they were found more than once. It is possible that this is due to an anomaly in the Btau_3.1 assembly, although it cannot be ruled out that those gene models do indeed represent duplicated genes that are present in the genome. 67 TRDV1, three TRDV2, two TRDV3 and one TRDV4 genes were identified; after consideration of multiple gene models that might be representative of a single gene we propose the presence of 52 TRDV1, two TRDV2, one TRDV3 and one TRDV4 genes. The number of TRDV2, TRDV3 and

TRDV4 genes is fairly consistent with what has been previously reported [17,20], although evidence for two TRDV3 genes has been presented by others [20]. In addition to the TRDV genes, five TRDD genes, three TRDJ genes and a single TRDC gene were identified, also consistent with previous reports [17,20].

The number of functional TRD genes among mammalian species is compared in Table 2. All species evaluated have genes for at least three TRDV subgroups. It is notable that humans have only a single TRDV1 gene (the mouse genes most closely related to human TRDV1 are named TRDV2-1 and TRDV2-2) while the TRDV1 subgroup in cattle and other artiodactyls are large, multigene families. Because of the disparity in numbers of TRDV1 genes in the human and cattle genomes, organization of their TRA/TRD loci cannot be compared. Furthermore, additional assembly of the bovine genome is required in order to definitively determine the size of the TRA/TRD locus, to refine the gene organization, and to evaluate duplication events that resulted in the large TRDV1 subgroup. Schematic representations for regions in which three or more TRDV genes were identified are shown (Figure 2A). Schematics are shown to scale and gene orientation is as indicated, based on the Btau_3.1 assembly. As in other species, a single TRDC gene (see Figure 2B) as well as one TRDV gene located downstream of TRDC in an inverted transcriptional orientation (named TRDV4 in cattle) were identified [17]. In some cases multiple gene models
<table><caption>Table 1 Annotated TRD genes</caption><thead><tr><th>Gene designation<sup>1</sup></th><th>GLEAN number<sup>2</sup></th><th>Scaffold Identification Number</th><th>Start</th><th>End</th><th>Orientation within the scaffold</th></tr></thead><tbody><tr><td>TRDV1a</td><td>GLEAN_22158</td><td>ChrUn.139</td><td>314906</td><td>315465</td><td>+</td></tr><tr><td>TRDV1b</td><td>GLEAN_22157</td><td>ChrUn.139</td><td>303549</td><td>304112</td><td>+</td></tr><tr><td>TRDV1c<sup>55555</sup></td><td>GLEAN_22161</td><td>ChrUn.139</td><td>356577</td><td>357146</td><td>+</td></tr><tr><td>TRDV1d*</td><td>GLEAN_22150</td><td>ChrUn.139</td><td>75036</td><td>75601</td><td>-</td></tr><tr><td>TRDV1e*</td><td>GLEAN_22151</td><td>ChrUn.139</td><td>16460</td><td>17025</td><td>+</td></tr><tr><td>TRDV1f</td><td>GLEAN_22162</td><td>ChrUn.139</td><td>373301</td><td>373837</td><td>+</td></tr><tr><td>TRDV1g</td><td>GLEAN_17632</td><td>ChrUn.298</td><td>28570</td><td>29136</td><td>+</td></tr><tr><td>TRDV1h</td><td>GLEAN_17646</td><td>ChrUn.298</td><td>223460</td><td>224017</td><td>+</td></tr><tr><td>TRDV1i</td><td>GLEAN_17645</td><td>ChrUn.298</td><td>206764</td><td>207325</td><td>+</td></tr><tr><td>TRDV1j</td><td>GLEAN_17648</td><td>ChrUn.298</td><td>263421</td><td>263882</td><td>+</td></tr><tr><td>TRDV1k</td><td>GLEAN_02242</td><td>ChrUn.1758</td><td>21943</td><td>22497</td><td>-</td></tr><tr><td>TRDV1l</td><td>GLEAN_02239</td><td>ChrUn.1758</td><td>58028</td><td>58581</td><td>-</td></tr><tr><td>TRDV1m</td><td>GLEAN_18280</td><td>ChrUn.221</td><td>6077</td><td>6638</td><td>-</td></tr><tr><td>TRDV1n</td><td>GLEAN_18279</td><td>ChrUn.221</td><td>30283</td><td>30837</td><td>-</td></tr><tr><td>TRDV1o<sup>****</sup></td><td>GLEAN_18269</td><td>ChrUn.221</td><td>174860</td><td>175436</td><td>-</td></tr><tr><td>TRDV1p</td><td>NA</td><td>ChrUn.221</td><td>220755</td><td>221325</td><td>-</td></tr><tr><td>TRDV1q<sup>555</sup></td><td>GLEAN_18266</td><td>ChrUn.221</td><td>240734</td><td>241302</td><td>-</td></tr><tr><td>TRDV1r</td><td>GLEAN_18263</td><td>ChrUn.221</td><td>264437</td><td>265003</td><td>-</td></tr><tr><td>TRDV1s</td><td>GLEAN_05951</td><td>ChrUn.41</td><td>188111</td><td>188668</td><td>-</td></tr><tr><td>TRDV1t<sup>**</sup></td><td>GLEAN_05971</td><td>ChrUn.41</td><td>313360</td><td>313923</td><td>+</td></tr><tr><td>TRDV1u</td><td>GLEAN_05970</td><td>ChrUn.41</td><td>294987</td><td>295550</td><td>+</td></tr><tr><td>TRDV1v<sup>***</sup></td><td>GLEAN_05970</td><td>ChrUn.41</td><td>282485</td><td>282946</td><td>+</td></tr><tr><td>TRDV1w</td><td>GLEAN_05969</td><td>ChrUn.41</td><td>272858</td><td>273319</td><td>+</td></tr><tr><td>TRDV1x</td><td>GLEAN_05952</td><td>ChrUn.41</td><td>161524</td><td>162099</td><td>-</td></tr><tr><td>TRDV1y<sup>**</sup></td><td>GLEAN_05991</td><td>ChrUn.41</td><td>661894</td><td>662457</td><td>+</td></tr><tr><td>TRDV1z</td><td>GLEAN_05990</td><td>ChrUn.41</td><td>642343</td><td>642906</td><td>+</td></tr><tr><td>TRDV1aa</td><td>GLEAN_05984</td><td>ChrUn.41</td><td>480039</td><td>480605</td><td>+</td></tr><tr><td>TRDV1ab<sup>***</sup></td><td>GLEAN_05989</td><td>ChrUn.41</td><td>629990</td><td>630287</td><td>+</td></tr><tr><td>TRDV1ac</td><td>GLEAN_16627</td><td>ChrUn.158</td><td>178923</td><td>179480</td><td>-</td></tr><tr><td>TRDV1ad</td><td>GLEAN_16623</td><td>ChrUn.158</td><td>237111</td><td>237679</td><td>-</td></tr><tr><td>TRDV1ae</td><td>GLEAN_16625</td><td>ChrUn.158</td><td>209540</td><td>210077</td><td>-</td></tr><tr><td>TRDV1af</td><td>GLEAN_03708</td><td>ChrUn.129</td><td>414487</td><td>415056</td><td>+</td></tr><tr><td>TRDV1ag</td><td>GLEAN_03704</td><td>ChrUn.129</td><td>331405</td><td>331977</td><td>+</td></tr><tr><td>TRDV1ah</td><td>GLEAN_03700</td><td>ChrUn.129</td><td>300288</td><td>300855</td><td>+</td></tr><tr><td>TRDV1ai<sup>****</sup></td><td>GLEAN_03707</td><td>ChrUn.129</td><td>390952</td><td>391528</td><td>+</td></tr><tr><td>TRDV1aj</td><td>GLEAN_03683</td><td>ChrUn.129</td><td>5640</td><td>6395</td><td>+</td></tr><tr><td>TRDV1al<sup>°</sup></td><td>GLEAN_10351</td><td>ChrUn.857</td><td>4172</td><td>4726</td><td>+</td></tr><tr><td>TRDV1am</td><td>GLEAN_10353</td><td>ChrUn.857</td><td>19494</td><td>20055</td><td>+</td></tr><tr><td>TRDV1an<sup>§</sup></td><td>GLEAN_10354</td><td>ChrUn.857</td><td>28383</td><td>28939</td><td>+</td></tr><tr><td>TRDV1ao<sup>§§</sup></td><td>GLEAN_10356</td><td>ChrUn.857</td><td>50213</td><td>50768</td><td>+</td></tr><tr><td>TRDV1ap</td><td>GLEAN_10558</td><td>ChrUn.857</td><td>74024</td><td>74592</td><td>+</td></tr><tr><td>TRDV1aq</td><td>GLEAN_15702</td><td>Chr10.28</td><td>20379</td><td>20949</td><td>-</td></tr><tr><td>TRDV1ar<sup>****</sup></td><td>GLEAN_15705</td><td>Chr10.28</td><td>58187</td><td>58723</td><td>+</td></tr><tr><td>TRDV1as</td><td>GLEAN_15706</td><td>Chr10.28</td><td>84670</td><td>85131</td><td>+</td></tr><tr><td>TRDV1at<sup>°</sup></td><td>GLEAN_15706</td><td>Chr10.28</td><td>97187</td><td>97750</td><td>+</td></tr><tr><td>TRDV1au<sup>§</sup></td><td>GLEAN_15708</td><td>Chr10.28</td><td>115470</td><td>116026</td><td>+</td></tr><tr><td>TRDV1av<sup>§§</sup></td><td>GLEAN_15709</td><td>Chr10.28</td><td>132846</td><td>133401</td><td>+</td></tr><tr><td>TRDV1aw<sup>55555</sup></td><td>GLEAN_15699</td><td>Chr10.28</td><td>198285</td><td>198854</td><td>-</td></tr><tr><td>TRDV1ax</td><td>GLEAN_15696</td><td>Chr10.28</td><td>245707</td><td>246273</td><td>-</td></tr><tr><td>TRDV1ay<sup>§§§</sup></td><td>GLEAN_19588</td><td>ChrUn.2578</td><td>31691</td><td>32259</td><td>-</td></tr><tr><td>TRDV1az</td><td>GLEAN_06130</td><td>ChrUn.1227</td><td>79925</td><td>80480</td><td>-</td></tr></tbody></table>
<table><caption>**Table 1: Annotated TRD genes (Continued)**</caption><tr><td>TRDV1ba</td><td>GLEAN_06133</td><td>ChrUn.1227</td><td>39695</td><td>40263</td><td>-</td></tr><tr><td>TRDV1bb</td><td>GLEAN_06131</td><td>ChrUn.1227</td><td>62992</td><td>63563</td><td>-</td></tr><tr><td>TRDV1bc<sup>***</sup></td><td>GLEAN_26185</td><td>ChrUn.2188</td><td>10476</td><td>11031</td><td>+</td></tr><tr><td>TRDV1bd</td><td>GLEAN_26186</td><td>ChrUn.2188</td><td>33687</td><td>34255</td><td>+</td></tr><tr><td>TRDV1be<sup>5555</sup></td><td>GLEAN_21863</td><td>ChrUn.5529</td><td>9610</td><td>10170</td><td>-</td></tr><tr><td>TRDV1bf<sup>55555</sup></td><td>GLEAN_20556</td><td>ChrUn.4913</td><td>10675</td><td>11244</td><td>-</td></tr><tr><td>TRDV1bg</td><td>GLEAN_05967</td><td>ChrUn.41</td><td>250250</td><td>250849</td><td>+</td></tr><tr><td>TRDV1bh</td><td>GLEAN_04691</td><td>ChrUn.3566</td><td>5067</td><td>5528</td><td>+</td></tr><tr><td>TRDV1bj<sup>°</sup></td><td>GLEAN_04693</td><td>ChrUn.3566</td><td>19054</td><td>19618</td><td>+</td></tr><tr><td>TRDV1bj<sup>°°</sup></td><td>GLEAN_10738</td><td>ChrUn.3620</td><td>538</td><td>1129</td><td>+</td></tr><tr><td>TRDV1bk<sup>°°</sup></td><td>GLEAN_12353</td><td>ChrUn.3793</td><td>3915</td><td>4506</td><td>+</td></tr><tr><td>TRDV1bl<sup>°°</sup></td><td>GLEAN_22424</td><td>ChrUn.3970</td><td>9390</td><td>9945</td><td>+</td></tr><tr><td>TRDV1bm</td><td>GLEAN_20555</td><td>ChrUn.10149</td><td>1843</td><td>2409</td><td>-</td></tr><tr><td>TRDV1bn<sup>****</sup></td><td>GLEAN_08205</td><td>ChrUn.3458</td><td>16972</td><td>17538</td><td>-</td></tr><tr><td>TRDV1bo<sup>****</sup></td><td>GLEAN_19099</td><td>ChrUn.2982</td><td>15173</td><td>15739</td><td>+</td></tr><tr><td>TRDV1bp<sup>****</sup></td><td>GLEAN_05881</td><td>ChrUn.4123</td><td>14961</td><td>15497</td><td>+</td></tr><tr><td>TRDV2-1</td><td>NA</td><td>ChrUn.158</td><td>345140</td><td>345657</td><td>-</td></tr><tr><td>TRDV2-2<sup>°</sup></td><td>GLEAN_10313</td><td>ChrUn.1907</td><td>43992</td><td>44509</td><td>+</td></tr><tr><td>TRDV2-2b<sup>°</sup></td><td>GLEAN_25296</td><td>ChrUn.1948</td><td>11302</td><td>11819</td><td>+</td></tr><tr><td>TRDV3-1<sup>nn</sup></td><td>NA</td><td>ChrUn.7503</td><td>3142</td><td>3667</td><td>-</td></tr><tr><td>TRDV3-1b<sup>nn</sup></td><td>GLEAN_19095</td><td>ChrUn.1002</td><td>77577</td><td>78102</td><td>+</td></tr><tr><td>TRDV4</td><td>GLEAN_19724</td><td>Chr10.30</td><td>396663</td><td>397245</td><td>+</td></tr><tr><td>TRDD1</td><td>NA</td><td>Chr10.30</td><td>504600</td><td>504612</td><td>-</td></tr><tr><td>TRDD2</td><td>NA</td><td>Chr10.30</td><td>488743</td><td>488757</td><td>-</td></tr><tr><td>TRDD3</td><td>NA</td><td>Chr10.30</td><td>449660</td><td>449672</td><td>-</td></tr><tr><td>TRDD4</td><td>NA</td><td>Chr10.30</td><td>441241</td><td>441249</td><td>-</td></tr><tr><td>TRDD5</td><td>NA</td><td>Chr10.30</td><td>417750</td><td>417760</td><td>-</td></tr><tr><td>TRDJ1</td><td>NA</td><td>Chr10.30</td><td>416717</td><td>416770</td><td>-</td></tr><tr><td>TRDJ2</td><td>NA</td><td>Chr10.30</td><td>406656</td><td>406598</td><td>-</td></tr><tr><td>TRDJ3</td><td>NA</td><td>Chr10.30</td><td>410369</td><td>410417</td><td>-</td></tr><tr><td>TRDC</td><td>GLEAN_19705</td><td>Chr10.30</td><td>402706</td><td>403988</td><td>-</td></tr></table>

<table><caption>**Table 2 Number of TRD genes among species**</caption><thead><tr><th rowspan="3">Locus</th><th colspan="6">TRD</th></tr><tr><th colspan="2">TRDV</th><th rowspan="2">TRDD</th><th rowspan="2">TRDJ</th><th rowspan="2">TRDC</th><th rowspan="2">Reference</th></tr><tr><th>Group</th><th>TRDV1</th><th>TRDV2 to TRDV5</th></tr></thead><tbody><tr><td>Human</td><td>1</td><td>2</td><td>3</td><td>4</td><td>1</td><td>IMGT http://www.imgt.org and [14]</td></tr><tr><td>Mouse</td><td>1</td><td>5</td><td>2</td><td>2</td><td>1</td><td>IMGT http://www.imgt.org and [10]</td></tr><tr><td>Bovine</td><td>52</td><td>4</td><td>5</td><td>3</td><td>1</td><td>Described here</td></tr><tr><td>Ovine</td><td>&gt; 40</td><td>3</td><td>?</td><td>4</td><td>1</td><td>[15,23]</td></tr><tr><td>Swine</td><td>&gt; 31</td><td>6</td><td>?</td><td>4</td><td>1</td><td>[21]</td></tr></tbody></table>

were found to represent the same gene (see Table 1, as indicated above, and Figure 2) and therefore subsequent analyses were performed after removing redundant sequences. It is likely that the occurrence of redundant sequences is a result of assembly anomalies; however, it is also possible that these duplicated TRDV genes are in fact present within the bovine genome.

**TRDV gene analysis**

Deduced amino acid sequences of the 56 TRDV genes identified here were initially aligned using ClustalW2 [28]. Alignments were refined in BioEdit according to the IMGT unique numbering for the V-DOMAIN [29], visualized with Jalview [30] and are shown (Figure 3A) with IMGT numbering and framework regions (FR) and

<footnote n="1">¹Gene designations with the same symbols denote gene models that are thought to be representative of the same gene based on analysis of genomic sequence.</footnote>

<footnote n="2">²NA denotes identified genes lacking glean models.</footnote>
A

<figure><img src="image_5.png" /></figure>

B

<figure><img src="image_6.png" /></figure>

**Figure 2 Schematic representation of the bovine (Bos taurus) TRD locus organization.** (A) 52 unique TRDV1 genes were identified and are all predicted to be functional (or ORF) and found on chromosome 10, although many of them were unplaced in the Bttau_3.1 assembly. Multiple TRDV genes were identified in two or three gene prediction models, as indicated, and labels correspond to those used in Table 1. However, evidence suggests that in each case only a single gene is represented. Genomic organization is shown in those cases when three or more TRDV genes were found within a single scaffolded region. TRA genes are not shown although in many cases occur among the TRDV genes. In addition to some TRDV1 genes, the schematic includes TRDV2-1 on ChrUn.158 (no other TRDV genes are included in the schematic). (B) Genomic organization of the five TRDD genes, 3 TRDJ genes and the single TRDC gene is shown, along with TRDV4 and TRAC. Gene designations, orientations and Bovine Genome Scaffold identifications are as indicated. The determination of TRDV gene orientation was based on the assembly of the Bovine Genome Scaffolds, it is possible that some scaffolds were assembled in an incorrect orientation individually and/or in relation to other scaffolds. Diagrams are shown to scale with base pair increments beneath the schematics. With the exception of TRAC, TRA genes are not shown. TRDV2-2 and TRDV3 are not shown because those genes were identified on scaffolds lacking additional TRDV genes.
<figure><img src="image_7.png" /><figcaption>**Figure 3** TRDV gene sequences. (A) TRDV deduced amino acid sequences were aligned with ClustalW2 using the default parameters and visualized with JalView. Analysis included all non-redundant genomic sequences (see Table 1 for GLEAN identification numbers). IMGT unique numbering for V-DOMAIN [29,32] is indicated above the alignment and conserved cysteines are indicated below the alignment. (B) IMGT Collier de Perles [33] are shown for TRDV1 (based on TRDV1a, GLEAN_22158), TRDV2-1 (no glean identification), TRDV3-1 (no glean identification) and TRDV4 (GLEAN_19724) and were determined using IMGT/V-QUEST [32]. The CDR-IMGT lengths are indicated.</figcaption></figure>
CDR, as indicated. TRDV1 genes reported here are tem-
porarily designated TRDV1a - TRDV1bp until the entire
chromosomal sequence has been assembled. TRDV1
genes encode CDR1 that range in length from five to
ten amino acids (except for TRDV1o which encodes an
18 amino acid long CDR1) while TRDV2, TRDV3 and
TRDV4 genes all encode CDR1 that are seven amino
acids in length. In contrast, most CDR2 of TRDV1,
TRDV2 and TRDV3 genes were three amino acids in
length, except for four of the TRDV1 genes as noted
below, while that of TRDV4 was five amino acids long.
In all but two cases (TRDV1bc and TRDV1w), TRDV1
and TRDV2 genes encode a QXS motif in the CDR2.
Interestingly, TRDV1f, TRDV1ae, TRDV1ar and
TRDV1o (see Figure 3A) all contained a deletion which
spans the last FR2 amino acid to the fifth FR3 amino
acid, resulting in a lack of CDR2 in those genes. An
identical deletion had also been found in the sheep
TRDV1S29 gene [15]. Also, the RS of TRDV1bg and
TRDV1bj were found to occur farther downstream than
those of the other TRDV1 genes identified here, result-
ing in a putative transcript that extends 15 (for
TRDV1bg) or 12 (for TRDV1bj) amino acids beyond the
conserved 2nd-CYS 104, as opposed to extending only
three or four amino acids which is the case for the
other TRDV1 genes. Because of the unusual nature of
those six genes, until their expression by γδ T cells is
validated they will be considered ORF instead of func-
tional as defined by IMGT [31]. 2D structure graphical
representations of TRDV1, TRDV2, TRDV3 and TRDV4
genes (Figure 3B) were generated using IMGT/V-
QUEST [32] and the IMGT/Collier-de-Perles tool pro-
gram [33]. The CDR1 and CDR2 loops are indicated in
red and yellow, respectively. The amino acids 1st-CYS
23, CONSERVED-TRP 41, hydrophobic amino acid
(Leu) 89 and 2nd-CYS 104, which are conserved in all
four TRDV subgroups, are shown with letters in red.

As for sheep, genes in the bovine TRDV1 family share
many characteristics including a 20 amino acid long
highly conserved leader region and the YFC motif found
in the 3' end of FR3 (see Figure 3A). However, bovine
TRDV1 genes described here were further classified into
eleven sets by constructing a phylogenetic tree using the
neighbor-joining method allowing us to compare them
to equivalent human, mouse and ovine sequences
(Figure 4). The grouping of ovine TRDV1 genes by this
method corresponded with what had been previously
described [15]. The tree shows grouping of bovine
TRDV1 gene sequences into eleven sets, as indicated,
and the percentage interior test values (based on 1000
replicates) are shown for each set. Classification was
further supported based on CDR characteristics used to
classify ovine TRDV1 genes [15] (Table 3). Previously
described characteristics used to classify TRDV1 genes

included CDR1 length (five, seven, nine or ten amino
acids), the chemical characteristics of the amino acid at
position 57 (within CDR2) and the presence or absence
of Trp 107 in CDR3. The number of TRDV1 genes eval-
uated here that occur within each set is also shown
(Table 4) with set 1 being the largest and set 4 being
the smallest. While some sets can share the same fea-
tures reported in Table 3 (i.e. sets 5 and 8, 4 and 6) the
classification of TRDV1 genes into those separate sets is
supported based on the phylogenetic analysis (Figure 4)
and additional sequence features (data not shown).

When bovine TRDV1 RS (total of 52 sequences) and
ovine TRDV1 RS (total of 12 sequences) were evaluated,
a high level of sequence conservation was observed both
among and within species, as depicted in the sequence
logos shown in Figure 5A. When the RS of bovine
TRDV4 orthologs were compared (including those for
human TRDV3, mouse TRDV5 and ovine TRDV4) a
similarly high level of sequence conservation was
observed among species (Figure 5B). In contrast, com-
parison of the remaining TRDV RS (those for bovine
TRDV2 and TRDV3, ovine TRDV2 and TRDV3, mouse
TRDV1, TRDV2 and TRDV4, human TRDV1 and
TRDV2) revealed conservation in only the heptamer and
nonamer sequences (Figure 5C). This is consistent with
the lack of relatedness of these genes as determined
based on sequence analysis (refer to Figure 4). Overall,
sequence conservation was observed within the RS hep-
tamer and nonamer of all TRDV genes, regardless of
species (compare Figure 5A, Figure 5B and Figure 5C).

TRDD and TRDJ gene analysis

Genomic locations, sequences and RS of bovine TRDD
and TRDJ genes were also determined (Table 1 and
Additional File 1) and confirmed what had been
reported previously [20]. Evidence for five TRDD genes
was found and deduced amino acid sequences for all
three open reading frames are shown in Additional File
1. Additional putative TRDD genes were identified
based on potential RS (data not shown); however, subse-
quent cDNA sequence analyses provided little or no evi-
dence that those putative TRDD genes were valid and
thus they are not presented. Three TRDJ genes were
found and their nucleotide and deduced amino acid
sequences are also reported in Additional File 1. Align-
ments comparing bovine, ovine, swine, human and
mouse TRDJ gene nucleotide and deduced amino acid
sequences are shown (Figure 6A and Figure 6B). A phy-
logenetic tree (Figure 6C), constructed using the neigh-
bor-joining method to evaluate the relatedness of the
above sequences, reflected the alignment.

Bovine, human and murine 5'D-RS and 3'D-RS of TRDD
genes and J-RS of TRDJ genes were evaluated for sequence
conservation (data not shown). Despite a low number of
<figure><img src="image_8.png" /></figure>

sequences used in the comparisons (because few TRDD and TRDJ genes exist in the genomes of those species) there appeared to be conserved heptamer and nonamer sequences among the RS for all species evaluated. Although heptamers were found to vary, as expected, they all contained the 5' CAC (or 3' GTG for 5'D-RS and J'RS) consensus sequence which has been found to be critical

for recombination [25]. Nonamers generally contained the previously reported consensus sequence ACAAAAACC [34] (or GGTTTTTGT for 5'D-RS and J'RS) although this was found to be less conserved. This is consistent with findings that the nonamer sequence requirements are less rigid than those for the heptamer sequence in order to obtain efficient recombination [25].
<table><caption>**Table 3 Features of bovine TRDV1 genes**</caption><thead><tr><th>Phylogenetic set</th><th>Relationship to previously named ovine sets [15]<sup>1</sup></th><th>CDR1 length</th><th>CDR2 (amino acid at IMGT 57)</th><th>CDR3 (presence of W at IMGT 107)</th><th>Exceptions</th></tr></thead><tbody><tr><td>1</td><td>1C</td><td>7</td><td>G</td><td>-</td><td></td></tr><tr><td>2</td><td>1B</td><td>7</td><td>Y</td><td>+/-</td><td>TRDV1t (CDR1 length of 8)</td></tr><tr><td>3</td><td>1D/1C</td><td>9</td><td>E</td><td>+/-</td><td></td></tr><tr><td>4</td><td>1A</td><td>9</td><td>Y</td><td>-</td><td></td></tr><tr><td>5</td><td>1A</td><td>9</td><td>Y</td><td>+</td><td>TRDV1d (Q at IMGT-107)</td></tr><tr><td>6</td><td>1A</td><td>9</td><td>Y or N</td><td>-</td><td></td></tr><tr><td>7</td><td>NA</td><td>5 or 10</td><td>A or V</td><td>+</td><td></td></tr><tr><td>8</td><td>1A</td><td>9</td><td>Y</td><td>+</td><td></td></tr><tr><td>9</td><td>1E</td><td>7</td><td>lacks CDR2</td><td>+/-</td><td>TRDV1o (CDR1 length of 18)</td></tr><tr><td>10</td><td>1E</td><td>7</td><td>N or G</td><td>+</td><td></td></tr><tr><td>11</td><td>1E</td><td>7</td><td>G or V</td><td>+</td><td></td></tr></tbody></table>

<table><caption>**Table 4 TRDV phylogenetic sets**</caption><thead><tr><th>TRDV gene</th><th>TRDV1 set<sup>1</sup></th><th>Bovine genes in sets<sup>2</sup></th><th>CDR-IMGT lengths</th><th>Number of TRDV genes in each set<sup>2</sup></th></tr></thead><tbody><tr><td rowspan="10">TRDV1</td><td>1</td><td>1a, 1ac, 1af, 1b, 1h, 1k, 1l, 1m, 1n, 1s</td><td>[7.3.4]</td><td>10</td></tr><tr><td rowspan="3">2</td><td>1am, 1an</td><td>[7.3.4]</td><td rowspan="3">4</td></tr><tr><td>1az</td><td>[7.3.3]</td></tr><tr><td>1t</td><td>[8.3.4]</td></tr><tr><td>3</td><td>1aa, 1ah, 1ax, 1bm, 1bn, 1g, 1r</td><td>[9.3.4]</td><td>7</td></tr><tr><td>4</td><td>1c, 1aq</td><td>[9.3.4]</td><td>2</td></tr><tr><td>5</td><td>1aj, 1ap, 1ba, 1bd, 1d</td><td>[9.3.4]</td><td>5</td></tr><tr><td rowspan="3">6</td><td>1ag, 1aw, 1q</td><td>[9.3.4]</td><td rowspan="3">5</td></tr><tr><td>1bg</td><td>[9.3.15]</td></tr><tr><td>1bj</td><td>[9.3.12]</td></tr><tr><td rowspan="2">7</td><td>1bb</td><td>[10.3.4]</td><td rowspan="2">3</td></tr><tr><td>1bc, 1ao</td><td>[5.3.4]</td></tr><tr><td>8</td><td>1ad, 1p, 1x</td><td>[9.3.4]</td><td>3</td></tr><tr><td rowspan="2">9</td><td>1ae, 1ar, 1f</td><td>[7.0.4]</td><td rowspan="2">4</td></tr><tr><td>1o</td><td>[18.0.4]</td></tr><tr><td>10</td><td>1i, 1al, 1u, 1z</td><td>[7.3.4]</td><td>4</td></tr><tr><td>11</td><td>1w, 1j, 1bh, 1v, 1as</td><td>[7.3.4]</td><td>5</td></tr><tr><td>TRDV2</td><td>NA</td><td>2-1, 2-2</td><td>[7.3.5]</td><td>2</td></tr><tr><td>TRDV3</td><td>NA</td><td>3-1</td><td>[7.3.5]</td><td>1</td></tr><tr><td>TRDV4</td><td>NA</td><td>4</td><td>[7.5.4]</td><td>1</td></tr></tbody></table>

### cDNA analysis as evidence for TRDV gene usage

Prior to these studies the number of TRDV1 genes found in the bovine genome, and their sequences, had not been resolved. Therefore, it was impossible to classify cloned TRDV1 sequences based on their encoding genes. Furthermore, TRDV1 genes are highly polymorphic but the subgroup also contains genes that differ by only one or three nucleotides in their coding regions (for example TRDV1aa and TRDV1bn, TRDV1v and TRDV1as). Therefore, it has been very difficult to differentiate between expressed TRDV1 sequences with regard to whether they represent the same or different genes. Here we classified previously-identified TRDV1

cloned cDNA sequences based on several criteria which included their relative placement in trees based on nucleotide and deduced amino acid sequences, specific sequence characteristics and percent identities based on nucleotide sequence alignments (data not shown). Results are summarized in Table 5 with the provisional IMGT nomenclature, corresponding genes and TRDV1 sets (as defined in Figure 4 and Tables 3 and 4) indicated.

In some cases it was impossible to classify TRDV1 cDNA sequences definitively, because some TRDV1 genes share very high sequence identity; those cases are indicated in Table 5. Furthermore, in 20 cases (see

<footnote n="1">¹NA indicates a bovine set lacking a related ovine set.</footnote>

<footnote n="1">¹NA indicates TRDV subgroup for which no set is defined.</footnote>

<footnote n="2">²For TRDV1 subgroup genes or subgroups (for the others).</footnote>
<figure><img src="image_9.png" /></figure>

Table 5) no corresponding TRDV1 genes could be identified although cDNA sequences were found to have features characteristic of particular sets, as indicated. It is possible that no corresponding genes were found because additional TRDV1 genes are present in the bovine genome but were not identified in this study due to gaps in the analyzed regions of the Btau_3.1 assembly. Only cDNA sequences that differed from potential

corresponding gene sequences by six or fewer nucleotides were classified; however, it remains possible that those cDNA sequences are not in fact representative of the genes indicated. Overall, corresponding genes were identified for most cDNA sequences (35 out of 55 sequences), although 36 of the 52 predicted TRDV1 genes reported here lacked cDNA confirmation and those genes were found to be distributed among all sets.
A

<table><tr><td>huTRDJ1</td><td>--------ACACCGATAAACTCATCTTTGGAAAAGGAACCCGTGTGACTGTGGAACCAA</td><td>51</td></tr><tr><td>moTRDJ1</td><td>--------CT......C......G......C......AA..................</td><td>51</td></tr><tr><td>swTRDJ1</td><td>------ATG.T..A..................G..T.AGC..GT...........</td><td>54</td></tr><tr><td>boTRDJ1</td><td>------ATG.G..T..C..G...........G......CC.C.TC...........</td><td>54</td></tr><tr><td>ovTRDJ1</td><td>------G.G..T..C..G....C...........G......CC.C.TC...........</td><td>52</td></tr><tr><td>huTRDJ4</td><td>----------AGACCC..G...........C...........TA.C..GAG..AC..A.C</td><td>48</td></tr><tr><td>swTRDJ4</td><td>----------CATCT.CC...G.CA..C..CC...........TA.C..AC..T......C</td><td>49</td></tr><tr><td>boTRDJ3</td><td>----------CAGT..CC...A..A......C...........TA.C..AC...........G</td><td>49</td></tr><tr><td>ovTRDJ3</td><td>----------CAGA..CC...A..A......C...........TA.C..AC...........G</td><td>49</td></tr><tr><td>huTRDJ2</td><td>------CTTTG..A.CAC......T...........G....A.AAC.C.TC...........</td><td>54</td></tr><tr><td>swTRDJ2</td><td>--ATACACTG.TG.CAC.TT...........T..G..T.AAC.C.GG...........</td><td>57</td></tr><tr><td>huTRDJ3</td><td>CTCCTGGG.....CGAC.GA.GT.T..C....CT..C.T.AAAC.CTTC......G..CC</td><td>59</td></tr><tr><td>moTRDJ2</td><td>CTCCTGGG.....CGAC.GA.GT.T...........CT..C.TAGAGC.CTT......G..CC</td><td>59</td></tr><tr><td>swTRDJ3</td><td>CTCCTGGG.....CGAC.GA.GTAT...........GCT..C..AAAC.CTT......G..CC</td><td>59</td></tr><tr><td>ovTRDJ2</td><td>CTCCTGGG.....CGGC.GA.GT.T...........GCT..C..AAAC.CTTC......G..CC</td><td>59</td></tr><tr><td>boTRDJ2</td><td>CTCCTGGG.....CGAC.GA.AT.T...........GCT..C..AAAC.CTTC......G..CC</td><td>59</td></tr></table>

B

<table><tr><td>huTRDJ1</td><td>--TDKLIFGKGTRVTVEP</td><td>16</td></tr><tr><td>moTRDJ1</td><td>-------V..Q..Q......</td><td>16</td></tr><tr><td>swTRDJ1</td><td>--D............QLV...</td><td>17</td></tr><tr><td>boTRDJ1</td><td>--E............LI...</td><td>17</td></tr><tr><td>ovTRDJ1</td><td>------T......LI...</td><td>16</td></tr><tr><td>huTRDJ4</td><td>-----RP...........YLE.QQ</td><td>15</td></tr><tr><td>swTRDJ4</td><td>-----P.T..Q..YLN...</td><td>14</td></tr><tr><td>boTRDJ3</td><td>--QYP...........YLN...</td><td>16</td></tr><tr><td>ovTRDJ3</td><td>------P...........YLN...</td><td>14</td></tr><tr><td>huTRDJ2</td><td>--L.AQ.F......QLI...</td><td>17</td></tr><tr><td>swTRDJ2</td><td>--LMAHF..N..QLR...</td><td>17</td></tr><tr><td>huTRDJ3</td><td>SWD.RQMF..T.IKLF...</td><td>19</td></tr><tr><td>moTRDJ2</td><td>SWD.RQMF..T.IELF...</td><td>19</td></tr><tr><td>swTRDJ3</td><td>SWD.RQMY..A..KLF...</td><td>19</td></tr><tr><td>ovTRDJ2</td><td>SWD.RQMF..A..KLF...</td><td>19</td></tr><tr><td>boTRDJ2</td><td>SWD.RQIF..A..KLF...</td><td>19</td></tr></table>

C

<figure><img src="image_10.png" /><figcaption>**Figure 6 TRDJ gene sequence alignments and phylogenetic tree.** Alignments of TRDJ gene (A) nucleotide and (B) deduced amino acid sequences, using bovine (bo) genomic TRDJ sequences identified here and previously classified ovine (ov), swine (sw), human (hu) and murine (mo) TRDJ sequences (described in Methods) are shown with identities indicated by a dot (.) and gaps indicated by a dash (-). The Neighbor-Joining method [47] was used to infer evolutionary history of the above mentioned sequences. The optimal tree (C) with the sum of branch length = 2.43145392 is shown. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) are shown next to the branches [48]. The Maximum Composite Likelihood method [49] was used to compute evolutionary distances. Complete deletion to eliminate gaps was performed and the final dataset included a total of 42 positions.</figcaption></figure>

Last, for the sake of completeness we note that cDNA evidence for bovine TRDV2, TRDV3 and TRDV4 genes reported here has already been demonstrated [17] and has been confirmed [20].

Identification of the TRDD genes within the bovine genome allowed us to evaluate TRDD gene usage in cDNA sequences that had been previously sequenced in our laboratory. CDR3 of 73 TRDV1, 17 TRDV2, 17 TRDV3 and 36 TRDV4 containing sequences were analyzed (Table 6). All D genes were represented in transcripts containing genes of the four TRDV subgroups. P nucleotides were observed flanking untrimmed TRDD regions and extensive N-regions were frequently found between TRDD regions. Some cases were ambiguous because nucleotides could be attributed to either P- or N-region but also to the presence of more than one

TRDD gene. Therefore, we applied the criteria that at least six nucleotides of a particular TRDD gene must be present to claim usage of that gene with reasonable confidence. Examples of TRDD gene usage with nucleotide and deduced amino acid sequences are shown (Figure 7). These demonstrated the incorporation of one to five TRDD genes within a single CDR3. Only cases in which the order of TRDD genes in the CDR3 corresponded to that found in the genome were retained.

Discussion

In order to evaluate and characterize the bovine T cell receptor delta gene repertoire and their genomic organization we annotated the TRD genes in the Btau_3.1 assembly. Here we describe the existence of two TRDV2 genes, as well as a single gene each for
<table><caption>Table 5 Classification of bovine TRDV1 cDNA sequences</caption><thead><tr><th>cDNA clone name</th><th>Accession number<sup>1</sup></th><th>Provisional IMGT nomenclature<sup>2</sup></th><th>Closest germline gene found in this study<sup>3</sup></th><th>Number of nucleotide differences<sup>4</sup></th><th>TRDV set</th></tr></thead><tbody><tr><td>Vdelta1.36</td><td>EF175175</td><td>NN</td><td>TRDV1I</td><td>5</td><td>1</td></tr><tr><td>BVd1.26</td><td>U73393</td><td>TRDV1S11</td><td>TRDV1m</td><td>6</td><td>1</td></tr><tr><td>SB-C11</td><td>FJ907530</td><td>NN</td><td>TRDV1af</td><td>1</td><td>1</td></tr><tr><td>SB-C17</td><td>-</td><td>NN</td><td>TRDV1af</td><td>1</td><td>1</td></tr><tr><td>BTDV27</td><td>D16114</td><td>TRDV1S14</td><td>TRDV1af</td><td>4</td><td>1</td></tr><tr><td>BVd1.13</td><td>U73380</td><td>TRDV1S15</td><td>TRDV1a</td><td>2</td><td>1</td></tr><tr><td>BTDV2</td><td>D13656</td><td>TRDV1S18</td><td>TRDV1a</td><td>2</td><td>1</td></tr><tr><td>SB-B1</td><td>FJ907521</td><td>NN</td><td>TRDV1a</td><td>5</td><td>1</td></tr><tr><td>SB-C25</td><td>FJ907532</td><td>NN</td><td>TRDV1a</td><td>5</td><td>1</td></tr><tr><td>SB-C3</td><td>-</td><td>NN</td><td>TRDV1a</td><td>4</td><td>1</td></tr><tr><td>BTDV4</td><td>D13660</td><td>TRDV1S13</td><td>TRDV1h</td><td>1</td><td>1</td></tr><tr><td>SB-A5</td><td>FJ907519</td><td>NN</td><td>TRDV1h</td><td>6</td><td>1</td></tr><tr><td>SB-C1</td><td>FJ907526</td><td>NN</td><td>TRDV1h</td><td>1</td><td>1</td></tr><tr><td>SB-C15</td><td>-</td><td>NN</td><td>TRDV1h</td><td>1</td><td>1</td></tr><tr><td>SB-A11</td><td>FJ907520</td><td>NN</td><td>TRDV1h</td><td>2</td><td>1</td></tr><tr><td>SB-C23</td><td>-</td><td>NN</td><td>TRDV1h</td><td>1</td><td>1</td></tr><tr><td>SB-B7</td><td>FJ907524</td><td>NN</td><td>TRDV1b</td><td>5</td><td>1</td></tr><tr><td>BVd1.23</td><td>U73390</td><td>TRDV1S4</td><td>TRDV1b</td><td>0</td><td>1</td></tr><tr><td>SB-C7</td><td>FJ907528</td><td>NN</td><td>ND</td><td>NA</td><td>1</td></tr><tr><td>BTDV1</td><td>D16116</td><td>TRDV1S6</td><td>ND</td><td>NA</td><td>1</td></tr><tr><td>BTDV1</td><td>D13655</td><td>TRDV1S15</td><td>ND</td><td>NA</td><td>1</td></tr><tr><td>BVd1.14</td><td>U73381</td><td>TRDV1S7</td><td>ND</td><td>NA</td><td>1</td></tr><tr><td>SB-B6</td><td>FJ907523</td><td>NN</td><td>ND</td><td>NA</td><td>1</td></tr><tr><td>BTDV8</td><td>D16112</td><td>TRDV1S12</td><td>ND</td><td>NA</td><td>1</td></tr><tr><td>SB-A7</td><td>-</td><td>NN</td><td>ND</td><td>NA</td><td>1</td></tr><tr><td>SB-B11</td><td>FJ907525</td><td>NN</td><td>ND</td><td>NA</td><td>1</td></tr><tr><td>BTDV33</td><td>D13658</td><td>TRDV1S1</td><td>ND</td><td>NA</td><td>1</td></tr><tr><td>SB-C19</td><td>-</td><td>NN</td><td>ND</td><td>NA</td><td>1</td></tr><tr><td>BTDV3</td><td>D13657</td><td>TRDV1S16</td><td>ND</td><td>NA</td><td>1</td></tr><tr><td>SB-B3</td><td>FJ907522</td><td>NN</td><td>ND</td><td>NA</td><td>1</td></tr><tr><td>BTDV28</td><td>D16115</td><td>TRDV1S17</td><td>ND</td><td>NA</td><td>1</td></tr><tr><td>BVd1.24</td><td>U73391</td><td>TRDV1S2</td><td>ND</td><td>NA</td><td>1</td></tr><tr><td>SB-C9</td><td>FJ907529</td><td>NN</td><td>TRDV1an</td><td>2</td><td>2</td></tr><tr><td>BVd1.16</td><td>U73383</td><td>TRDV1S9</td><td>TRDV1an</td><td>1</td><td>2</td></tr><tr><td>BVd1.18</td><td>U73385</td><td>TRDV1S10</td><td>TRDV1an</td><td>5</td><td>2</td></tr><tr><td>BVd1.15</td><td>U73382</td><td>TRDV1S3</td><td>TRDV1am</td><td>0</td><td>2</td></tr><tr><td>SB-C5</td><td>FJ907527</td><td>NN</td><td>TRDV1az<sup>5</sup></td><td>4</td><td>2</td></tr><tr><td>SB-C13</td><td>-</td><td>NN</td><td>TRDV1az<sup>5</sup></td><td>4</td><td>2</td></tr><tr><td>SB-A3</td><td>-</td><td>NN</td><td>TRDV1az</td><td>1</td><td>2</td></tr><tr><td>BVd1.17</td><td>U73384</td><td>TRDV1S8</td><td>TRDV1az</td><td>1</td><td>2</td></tr><tr><td>SB-A9</td><td>-</td><td>NN</td><td>TRDV1az</td><td>2</td><td>2</td></tr><tr><td>SB-B9</td><td>-</td><td>NN</td><td>TRDV1az<sup>5</sup></td><td>4</td><td>2</td></tr><tr><td>BVd1.21</td><td>U73388</td><td>TRDV1S23</td><td>TRDV1ax</td><td>1</td><td>3</td></tr><tr><td>BVd1.20</td><td>U73387</td><td>TRDV1S22</td><td>TRDV1aa</td><td>2</td><td>3</td></tr><tr><td>BTDV6</td><td>D13661</td><td>TRDV1S24</td><td>TRDV1r</td><td>1</td><td>3</td></tr><tr><td>Vdelta1.37</td><td>EF175176</td><td>NN</td><td>TRDV1aq</td><td>4<sup>6</sup></td><td>4</td></tr><tr><td>Vdelta1.33</td><td>EF175171</td><td>NN</td><td>ND</td><td>NA</td><td>5</td></tr><tr><td>BTDV26</td><td>D16113</td><td>TRDV1S25</td><td>ND</td><td>NA</td><td>5</td></tr><tr><td>BVd1.22</td><td>U73389</td><td>TRDV1S21</td><td>ND</td><td>NA</td><td>5</td></tr><tr><td>BVd1.25</td><td>U73392</td><td>TRDV1S5</td><td>ND</td><td>NA</td><td>7</td></tr><tr><td>Vdelta1.34</td><td>EF175172</td><td>NN</td><td>ND</td><td>NA</td><td>7</td></tr></tbody></table>
<table><caption>Table 5: Classification of bovine TRDV1 cDNA sequences (Continued)</caption><tr><td>BVd1.19</td><td>U73386</td><td>TRDV1S20</td><td>TRDV1ad</td><td>1</td><td>8</td></tr><tr><td>Vdelta1.35</td><td>EF175173</td><td>NN</td><td>TRDV1ae</td><td>0</td><td>9</td></tr><tr><td>SB-C21</td><td>FJ907531</td><td>NN</td><td>TRDV1u</td><td>5</td><td>10</td></tr><tr><td>BTDV36</td><td>D13659</td><td>TRDV1S19</td><td>ND</td><td>NA</td><td>11</td></tr></table>

<table><caption>Table 6 TRDD gene usage in CDR3 of bovine rearranged sequences by TRDV genes of different subgroups</caption><thead><tr><th rowspan="2">TRDV subgroup</th><th rowspan="2">Number of sequences evaluated</th><th colspan="5">Number of times TRDD gene identified (Percentage of all TRDD genes used)</th></tr><tr><th>TRDD1</th><th>TRDD2</th><th>TRDD3</th><th>TRDD4</th><th>TRDD5</th></tr></thead><tbody><tr><td>TRDV1</td><td>73</td><td>33 (24.3)</td><td>33 (24.3)</td><td>18 (13.2)</td><td>23 (16.9)</td><td>29 (21.3)</td></tr><tr><td>TRDV2</td><td>17</td><td>6 (16.2)</td><td>9 (24.3)</td><td>7 (18.9)</td><td>9 (24.3)</td><td>6 (16.2)</td></tr><tr><td>TRDV3</td><td>17</td><td>9 (29.0)</td><td>3 (9.7)</td><td>10 (32.2)</td><td>4 (12.9)</td><td>5 (16.1)</td></tr><tr><td>TRDV4</td><td>36</td><td>9 (14.8)</td><td>16 (26.2)</td><td>11 (18.0)</td><td>16 (26.2)</td><td>9 (14.8)</td></tr></tbody></table>

<table><thead><tr><th>TRDV</th><th>P/N</th><th>TRDD1</th><th>P/N</th><th>TRDD2</th><th>P/N</th><th>TRDD3</th><th>P/N</th><th>TRDD4</th><th>P/N</th><th>TRDD5</th><th>P/N</th><th>TRDJ</th></tr></thead><tbody><tr><td>Germline</td><td></td><td>GTGGACTGGGGG</td><td></td><td>GGTGGGATTACGAG</td><td></td><td>GGTGGGGTACTAC</td><td></td><td>ACAACGTAC</td><td></td><td>GGTGGGATACG</td><td></td><td></td></tr><tr><td>SBA5</td><td>GCTCTCTGC</td><td></td><td></td><td>GGATTTCAC</td><td>G</td><td>GGTGGGGTACT</td><td></td><td></td><td></td><td>GGGACAGTAT</td><td>G Q Y</td><td>CCACTA P L</td></tr><tr><td>CE21</td><td>A L C</td><td></td><td></td><td>R I Y</td><td>G</td><td>G G V L</td><td></td><td></td><td></td><td>G G Y</td><td></td><td>AATCCACTA N P L</td></tr><tr><td>SBB7</td><td>GCTGGTACA GC</td><td></td><td></td><td></td><td></td><td>GCTGGGGTAC</td><td></td><td></td><td></td><td>GTGGGATAC</td><td>R G I R</td><td>CAAGCTC K L</td></tr><tr><td>CH442</td><td>GCTCTCCAT GAAGAAAGAGC</td><td>GTGGACTGGGGG</td><td>TAC</td><td></td><td></td><td>GGTGGGGT</td><td></td><td>ACGTAC</td><td></td><td>GGGATA</td><td></td><td>GAATAT E Y</td></tr><tr><td>CH335</td><td>A L H E E R A</td><td>L D W G</td><td>T</td><td></td><td></td><td>W W G</td><td></td><td>T Y</td><td></td><td>G I</td><td></td><td>CCACTA P L</td></tr><tr><td></td><td>GCTCTCTGC TCTAGTAGCGTACAA</td><td>GGACTGGG</td><td></td><td>ATTACGA</td><td>TAGCGGAGGTATC</td><td>GGTGTAC</td><td>ACCGTAC</td><td>TC</td><td>GGTGGGATAC</td><td>R W D T</td><td>CCCATATGAG ACTGACAAGCTC</td><td>P Y E T O K L</td></tr><tr><td></td><td>A L C S S S V Q</td><td>G L G</td><td></td><td>E T</td><td>I A E V Y</td><td>R V</td><td>T Y</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>SBVd27</td><td>GCTCTCTGG GACAGCGTAC</td><td></td><td></td><td>GGTGGGATTIA</td><td>V G F</td><td></td><td></td><td></td><td></td><td></td><td></td><td>ATCCACTA N P L</td></tr><tr><td>SBVd23</td><td>GCTCTCTGG GAAGGAACTTTA</td><td>GTGGACTGGG</td><td>AAGATC</td><td></td><td></td><td></td><td></td><td>CAACGT</td><td></td><td>GCAG</td><td>Q</td><td>AATCCACTA N P L</td></tr><tr><td>SBVd21</td><td>A L W E G T L</td><td>V G L G</td><td>R S</td><td></td><td></td><td>GGTGGGATTTCAC</td><td></td><td>CAACGTAC</td><td></td><td>GGTGGGAT</td><td>G G M</td><td>AATCCACTG N P L</td></tr><tr><td>CH285</td><td>GCTTCTCTGG GAAGTTTACA</td><td>V L W E V Y</td><td></td><td></td><td></td><td>R W D L</td><td></td><td>P T Y</td><td></td><td>GGTGGGATACG</td><td>R W D T</td><td>CCACTA P L</td></tr><tr><td></td><td>GCTCTCTGG GAAGCTATATGGGGAGCGGGCGTA</td><td>GACTGGGG</td><td></td><td></td><td></td><td>ATTACG</td><td></td><td>ACACGTAC</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>A L W E A I W G A G G</td><td>R L G</td><td></td><td></td><td></td><td>I Y</td><td></td><td>N V</td><td></td><td></td><td></td><td>ATGCAGTAT</td></tr><tr><td>CH03D3</td><td>GCTCTCAGA GAGTACAGCGTAC</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>CCACTA P L</td></tr><tr><td>CH289</td><td>A V R E Y S V</td><td>GGACTGGGGGG</td><td>R T G G</td><td></td><td></td><td></td><td></td><td></td><td></td><td>GGGATAC</td><td>G I</td><td>CCACTA P L</td></tr><tr><td>CH287</td><td>GCTCTCAGA GAGTTTG</td><td></td><td></td><td>GGTGGGATTTCAC</td><td>T A V</td><td></td><td></td><td>ACAACG</td><td></td><td>GGGATAC</td><td></td><td>CCACTA P L</td></tr><tr><td></td><td>A V R E L</td><td></td><td></td><td>G G I Y</td><td></td><td></td><td></td><td>T T</td><td></td><td></td><td></td><td>ATCAGTAT</td></tr><tr><td></td><td>GCTCTCAGA GAGTTTCC</td><td>GTGGACTGG</td><td>L D W</td><td></td><td></td><td></td><td></td><td></td><td></td><td>GGGATAC</td><td>G I</td><td>CCACTA P L</td></tr><tr><td>CH374</td><td>GCITTGAGT CTTAATC</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>ATCAGTAT</td></tr><tr><td>CH383</td><td>GCITTGAGG CTCCGGAGTACGG</td><td></td><td></td><td>GGGATTTCACGAG</td><td>G F T S</td><td>CCGCTTATTT</td><td></td><td></td><td></td><td>GGTGGGATAC</td><td>R W D T</td><td>AGCT A D K L</td></tr><tr><td>CH378</td><td>A L R L R E Y G</td><td>GTGGACTGGG</td><td></td><td>GGGATTTCACGAG</td><td>A L F</td><td>CTCT</td><td></td><td>ACACGTAC</td><td></td><td>GTGGGATAC</td><td>R G Y</td><td>CCACTA P L</td></tr><tr><td>CH367</td><td>GCITTGAGG GGCACAGCATCT</td><td>V G L G</td><td></td><td>G I Y E</td><td>L</td><td></td><td></td><td>N V</td><td></td><td>GTGGGATAC</td><td>G G M</td><td>CCACTA P L</td></tr><tr><td></td><td>A L R G T A S</td><td></td><td></td><td>GGGATTTCACGAG</td><td>G I Y E</td><td></td><td></td><td></td><td></td><td>GTGGGATAC</td><td>R W D T</td><td>CCACTA P L</td></tr><tr><td></td><td>GCITTGAGG CTCGGG</td><td></td><td></td><td>G G I Y E</td><td></td><td></td><td></td><td></td><td></td><td>GTGGGATAC</td><td>V G Y</td><td>CCACTA P L</td></tr></tbody></table>

**Figure 7 TRD rearranged CDR3 sequences and TRDD gene usage.** TRD rearranged CDR3 nucleotide sequences, and their corresponding deduced amino acid sequences, derived from cDNA from peripheral blood mononuclear cells, were aligned to demonstrate TRDD usage. Germline TRDD sequences are shown above representative TRDV1, TRDV2, TRDV3 and TRDV4 sequences containing between one and five TRDD genes. TRDD genes are shaded and gene usage was determined by the presence of at least six nucleotides of a particular TRDD gene. Sufficient cDNA sequence was analyzed to determine TRDV gene usage; however, the sequences shown here have been truncated after the 2nd-CYS 104 and at the beginning of the TRDJ gene.

TRDV3 and TRDV4. We confirmed the presence of three TRDJ and five TRDD genes [20] and evaluated TRD junctional diversity incorporating these genes. Furthermore, we described the existence of 52 TRDV1 genes (46 of them functional and six of them ORF), evaluated them for their sequence characteristics and further classified them into eleven sets based on

phylogenetic analysis. As for other mammals, the bovine TRD locus was found embedded within the TRA locus and thus TRDV1 genes were distributed amongst the TRAV genes. Indeed, there were a number of genes for which assignment to the TRAV or TRDV group was ambiguous. It is important to note that the presence of TRAV/DV genes is certain and

<footnote n="1">¹A dash indicates redundant sequence for which a GenBank accession number was not obtained.</footnote>

<footnote n="2">²NN indicates sequence for which no Provisional Nomenclature exists.</footnote>

<footnote n="3">³ND indicates cDNA sequences for which no corresponding gene could be identified in this study.</footnote>

<footnote n="4">⁴Number of nucleotides that are different between the Provisional IMGT Nomenclature and the corresponding gene sequences between amino acid positions 1 and 104 (2nd-CYS); NA indicates sequence that lacks corresponding gene.</footnote>

<footnote n="5">⁵Genes listed as corresponding are the closest matches; however, classification of these genes was unclear.</footnote>

<footnote n="6">⁶cDNA sequence used for comparison was partial.</footnote>
has been reported [17] but the ability of TRDV or
TRAV genes (as defined based on sequence character-
istics) to rearrange with either TRDD or TRAJ must
be determined experimentally and will be evaluated in
future studies. Evaluation of the genes contained
within the TRA/TRD locus was reported [22] following
the initial submission of this work with overall findings
similar to those discussed here.

The question remains as to whether or not the find-
ings reported here represent the complete TRD gene
repertoire because of the gaps in the Btau_3.1 assembly.
The presence of additional TRDV1 genes is almost cer-
tain because of extensive gaps in this region although it
is unlikely that additional TRDD or TRDJ genes are pre-
sent in the genome since that region is fully scaffolded
and has been carefully evaluated for the presence of
additional genes. That region also contains TRDV4 and
its analysis is likely definitive. In addition, previous
cDNA evidence supports the existence of only two
TRDV2 genes [17] and is consistent with our findings.
The existence of two TRDV3 genes has been reported
[20] but is not supported by our findings which may be
due to gaps in the genome assembly. Even with gaps,
locus organization can be assisted by comparison to
other species. Indeed, all mammals studied have genes
that are closely related to human TRDV1 with some
being distributed among the TRAV genes. The bovine
TRDV4 gene is located downstream of the single TRDC
gene in inverse orientation, as found for the human
TRDV3 and mouse TRDV5 genes. It is notable that
there is significant sequence conservation among the
orthologous bovine TRDV4, human TRDV3, mouse
TRDV5 and ovine TRDV4 genes as well as their entire
RS. This suggests that they derive from a common
ancestral gene that appeared before the separation of
these species. By comparison, there are no human or
mouse genes that are equivalent to bovine TRDV2 or
TRDV3. The enormous expansion of the TRDV1 reper-
toire that has been observed for cattle, sheep and pig
[15,21] is striking when compared with the single
TRDV1 gene found in human (and its two orthologs in
mouse) and thus does not assist us in predicting the
extent of the possible missing information.

TRDV1 gene subgroup expansion is also striking when
compared with other TRDV gene subgroups which, even
in artiodactyls, are found to contain only one or two
members. While the number of TRGV genes found in
each species does not vary considerably [31,35,36],
expansion of the TRDV1 genes in bovine, sheep and pig
is consistent with them being “γδ T cell high” species,
implying that this increase in antigen receptor diversity
permits them to be stimulated by and to respond to
more types of antigens. The presence of fewer TRDV
genes in “γδ T cell low” species, such as humans and

mice, suggests that these species evolved so that γδ
T cells play a less important role in their immune func-
tion. Sequence conservation among bovine TRDV1 RS
(including the 23-bp spacer sequences), supports the
concept that the large TRDV1 subgroup arose from
multiple duplication events. Many bovine TRDV1 genes
share high percent identities, indicating that in some
cases this subgroup is not as diverse as the number of
individual genes might suggest. For example, some
TRDV1 sequences differed by as few as four (TRDV1m
and TRDV1s), three (TRDV1v and TRDV1as) or even
one (TRDV1aa and TRDV1bn) nucleotide in their cod-
ing regions even though all were determined to be dis-
tinct genes based on analysis of the flanking genomic
sequences. However, it cannot be ruled out that some of
these gene models represent the same gene but do not
appear as such due to sequencing and/or assembly
anomalies. The question of why numerous and highly
similar TRDV1 genes have been maintained in the
bovine genome is compelling and studies to evaluate
their functions will be pursued.

Based on phylogenetic sequence analysis and gene features, we subdivided the bovine TRDV1 genes into eleven sets. Each set corresponds to intraspecies gene duplications of ancestral genes that were present before the separation between bovine and ovine species. We also assigned bovine cDNA sequences to germline genes identified in this work (see Table 5). Of 55 previously identified TRDV1 cDNA sequences evaluated, 20 could not be assigned to a germline gene, possibly because of gaps in that region of the genome assembly. Of the 52 TRDV1 genes identified in this study, cDNA evidence existed for only 16. It is possible that, because a majority of the previously identified TRDV1 cDNA sequences were derived from peripheral blood mononuclear cells [16-20], the remaining TRDV1 genes that lack cDNA verification might be expressed in other tissues that have not been evaluated. It is also possible that other factors, such as age or immunological experience, might impact the expression of particular genes and thus prevented verification of those genes. Future studies will focus on specifically amplifying individual TRDV1 transcripts in order to verify that they are functional and to determine whether they are expressed in an age and/or tissue dependent manner.

The potential importance of the role of the T cell
receptor delta rearranged CDR3 is evident instudies
evaluating the binding of the T10 and T22 antigens
[26,37,38] by murine γδ T cells since antigen recognition
by the CDR3 was found to be autonomous, that is, not
dependent on the TRDV gene associated with it [39]. It
could be reasoned that, because the genomes of humans
and mice lack the extensive TRDV repertoirethat has
been observed in artiodactyls, thiswould be compensated
for by diversity of the genes that rearrange to form the CDR3 (i.e. TRDD and TRDJ genes). Moreover, because many bovine TRDV1 amino acid sequences are very similar, structurally it might be expected that their ligand binding regions including CDR3 would also need to be as diverse as that in humans and mice. Humans have three TRDD and four TRDJ genes while mice have two TRDD and two TRDJ genes. With the exception of the human TRDJ genes, these numbers all fall short of the five TRDD and three TRDJ genes identified in cattle. We demonstrated here that in cattle the TRD CDR3 contained combinations that include between one and five TRDD genes and would range in length from nine to 37 amino acids. With the exception of TRDD4, the germline bovine TRDD genes predominantly encode glycine and the hydrophobic residues valine and tryptophan, with the neutral residues threonine and tyrosine being encoded to a lesser extent (see Additional File 1). This is reminiscent of findings for the equine VH CDR3 loops in which a higher proportion of glycine and lower proportion of cysteine content was identified when sequences were compared with those from humans, sheep and pigs [40]. As was suggested for equine VH CDR3, this indicates that bovine TRDV CDR3 loops have increased flexibility and are therefore better suited to recognize a large number of antigenic conformations. Considering this, it is unlikely that CDR3 diversity in mice and humans compensates for their less diverse TRDV repertoire and, in fact, it seems that cattle have more scope for antigen binding on all levels.

If the highly diverse CDR3 is indeed the main struc-
tural component involved in antigen recognition, then
the TRDV gene products might be most important for
interactions with co-receptors or antigen presenting
cells, whatever their nature may be. While the interac-
tions between γδ T cell receptor and the antigens that
they bind are still largely unknown, it has been demon-
strated that for αβ T cell receptor the CDR1 loops
(which are germline encoded) are involved with binding
peptide as well as MHC. The CDR2 loops (also germline
encoded) bind only MHC. The TRDV CDR1 in cattle
exhibited diverse lengths and were generally found to be
five to ten amino acids long as for mouse and human
TRDV, TRBV and TRAV CDR1. This diverse bovine
TRDV CDR1 repertoire may contribute to antigen bind-
ing in contrast to what has been reported for murine γδ
T cell recognition of T22 as described above. It was
found that bovine TRDV1 CDR2 were either only three
amino acids long or absent completely (e.g., TRDV1f,
TRDV1ae, TRDV1ar and TRDV1o) in contrast to TRBV
CDR2 lengths of five to seven amino acids in mammals.
It is yet to be determined whether the TRDV1 genes
lacking CDR2 are functional, however, if they are found
to be so this lack of, or abbreviation of, the CDR2

inbovine TRDV1 genes is logical since γδ T cells are not
MHC-restricted.

Conclusions

Based on annotations of the bovine genome we identi-
fied the TRD genes including 56 TRDV genes, five
TRDD genes, three TRDJ genes and the single TRDC
gene and described their organization within the TRD
locus. Furthermore, we report that the TRDV1 subgroup
contains 52 genes, indicating that this subgroup under-
went expansion as has been found for sheep and swine.
Also, because this large gene subgroup has been main-
tained in the bovine genome this indicates a unique role
for these genes in γδ T cell biology.

Methods

Genome Annotation

In conjunction with the Bovine Genome Sequencing
Consortium http://genomes.arc.georgetown.edu/
bovine/, manual annotation of the T cell receptor
delta genes was performed using the Apollo Genome
Annotation and Curation Tool, version 1.6.5 [41] and
the bovine genome assembly Btau_3.1 [42]. First, a
BLAST search of bovine T cell receptor delta cDNA
sequences against the Bovine Official Gene Set
(GLEAN) was performed in order to identify predicted
gene models. These were then analyzed using the
Apollo software and the following actions were per-
formed where applicable: (i) models were checked for
correct exon-intron structure, (ii) initiation and termi-
nation codons were identified where applicable and in
other cases 5' and 3' ends were set based on RS and
splice sites, (iii) exons were either added or deleted if
it was determined that the coding region in the pre-
dicted model was incorrect, (iv) predicted gene models
were split when a single model encompassed more
than one gene or merged when two models coded for
a single gene and (v) RS were identified for TRDV,
TRDD and TRDJ genes. Predicted gene models identi-
fied from the BLAST search were considered pseudo-
genes, and were not included in subsequent analyses,
when premature stop codons or frameshifts occurred
in areas where the sequence integrity was deemed ade-
quate. Furthermore, predicted gene models were con-
sidered Open Reading Frame (ORF), as defined by
IMGT [31] if the coding region had an open reading
frame but the sequence structure differed significantly
from known TRDV sequences; ORF gene models were
included in subsequent analyses. Following annotation,
predicted gene model identity and classification were
verified using BLAST searches and IMGT/V-QUEST
[32]; only gene models that were unambiguously iden-
tified as TRDV genes, and not TRAV/DV genes, based
on those results were included in subsequent analyses.
## Sequence analyses

Nucleotide sequences were aligned and analyzed using BioEdit version 7.0.5.3 [43]. Exon-intron structure schematics were based on alignments of cDNA and genomic DNA sequence using SIM4 [44] and visualized with LalnView http://pbil.univ-lyon1.fr/software/lalnview.html. GLEAN numbers of the annotated bovine gene sequences used in sequence analyses are reported in Table 1. Published T cell receptor delta sequences, derived from cloned RT-PCR products, were submitted to GenBank (http://www.ncbi.nlm.nih.gov/Genbank/index.html; see Table 5 for accession numbers) and assigned to TRDV genes identified here. CDR3 sequence data described in Table 6, from bovine TRD rearranged genes, are available upon request. Rearranged CDR3 lengths were determined from positions 105 to 117, that is, between the 2nd-CYS 104 and J-PHE 118 of the FGXG motif, according to the IMGT unique numbering for V-DOMAIN [29]. The accession numbers of additional sequences used in analyses are as follows: human (Homo sapiens): TRDV1 (M22198), TRDV2 (X15207), TRDV3 (M23326), TRDJ1 (M20289, AE000661), TRDJ2 (L36386, AE000661), TRDJ3 (M21508, AE000661), TRDJ4 (AJ249814, AE000661), TRDD1 (M23325), TRDD2 (M22153), TRDD3 (M22152); mouse (Mus musculus): TRDV1 (M94080), TRDV2-1 (AE008686), TRDV2-2 (AE008686), TRDV4 (AE008686), TRDV5 (AE008686), TRDJ1 (AF019412), TRDJ2 (X64903, NT_039614), TRDD1 (X64900), TRDD2 (X64901); pig (Sus scrofa): TRDJ1 (D49560, AB053451), TRDJ2 (D49561, AB053451), TRDJ3 (D49562, AB053451), TRDJ4 (D49563, AB053451); cattle (Bos taurus): TRAV (BC148926), TRGJ1-1; sheep (Ovis aries): TRDV1S1 (Z12970, AJ786827), TRDV1S2 (Z12981), TRDV1S3 (Z12988), TRDV1S4 (Z12989), TRDV1S5 (Z12990, AJ786828), TRDV1S7 (Z12992), TRDV1S8 (Z12993, AJ786830), TRDV1S9 (Z12994, AJ786831), TRDV1S10 (Z12971, AJ786832), TRDV1S11 (Z12972), TRDV1S12 (Z12973), TRDV1S13 (Z12974, AJ786833), TRDV1S14 (Z12975), TRDV1S15 (Z12976), TRDV1S16 (AJ786834), TRDV1S17 (Z12978), TRDV1S18 (Z12979, AJ786835), TRDV1S19 (Z12980), TRDV1S20 (Z12982), TRDV1S21 (Z12983), TRDV1S22 (Z12984), TRDV1S24 (Z12986), TRDV1S25 (Z12987), TRDV1S29 (AJ786836), TRDV1S30 (AJ786837), TRDV1S31 (AJ786838), TRDV1S32 (AJ786839), TRDV1S33 (AJ786840), TRDV1S34 (AJ786841), TRDV1S35 (AJ809502), TRDV1S36 (AJ809503), TRDV1S37 (AJ809504), TRDV1S38 (AJ809505), TRDV1S39 (AJ809506), TRDV1S40 (AJ809507), TRDV2 (Z12995), TRDV3 (Z12996), TRDV4 (Z12997, AJ810117), TRDJ1 (AJ277510), TRDJ2 (AJ277511), TRDJ3 (AJ277512).

Multiple sequence alignments were performed using ClustalW2 (http://www.ebi.ac.uk/Tools/clustalw2/index.html; [28]) and the default parameters and were visualized using JalView [30]. IMGT unique numbering was determined for the V-DOMAIN [29] and TRDV secondary structures were evaluated and displayed as IMGT Collier de Perles [33] as determined using IMGT/V-QUEST [32]. Sequence logos were generated using WebLogo3 [45].

## Phylogenetic analyses

Phylogenetic analyses were performed using nucleotide sequences, unless otherwise indicated, of functional (and ORF) TRDV and TRDJ genes. TRDV sequences were truncated to include the V-REGION (IMGT amino acid positions 1 through 108) and comprised those from bovine (identified here in Btau_3.1), human, mouse and ovine (see accession numbers above). TRDJ sequences were truncated to include sequence between the RS (5' end) and splice site (3' end) and comprised those from bovine (identified here in Btau_3.1), human, mouse, ovine and swine (see accession numbers above).

Phylogenies were constructed in MEGA4 [46] using the neighbor-joining method and the p-distance model for nucleotide sequences. Phylogenies were tested using an interior branch test or bootstrap analysis, as indicated, with 1000 replications.

## cDNA analysis

Bovine TRDV cDNA was obtained as described in previous studies [16,17]. Briefly, blood was collected via jugular venipuncture into heparin and peripheral blood mononuclear cells were isolated via density gradient centrifugation over Ficoll-Paque™ PLUS (GE Healthcare Bio-Sciences, Piscataway, NJ) according to the manufacturer's protocol. Animal use complied with federal guidelines and had IACUC approval. Pelleted cells were resuspended in TRIzol (Invitrogen, Carlsbad, CA) and RNA was isolated according to the manufacturer's protocol. Polymerase chain reaction (PCR) was performed following reverse transcription (RT) and using primers designed to specifically amplify TRDV transcripts as previously described [16,17]. In addition, amplification of TRD rearranged CDR3 was performed using a common reverse primer called TRDC-rev (5'-CTC CTT CAC CAA ACA AGC GAC G-3') and the following forward primers: TRDV1S1S2-for (5'-GCA GAT AAA TCC ATC AGC CTC ACC-3'), TRDV1S3-for (5'-CAC AGA ACT CCA TCA GCC TCA CC-3'), TRDV1S4-for (5'-TCA CGT AAA GCC ATC AGC CTC ATT-3'), TRDV2-for (5'-GGA AAA AAA ATC ATC AGC CTC ACC-3'), TRDV3-for (5'-TAC AAA CCC AAC CAA ATG CTG AAA-3'), TRDV4-for (5'-AGC ATG AGC CAA AAA ACC TTC CAC-3'). PCR products were analyzed on TAE agarose gels and cloned into the pCR2.1 vector (Invitrogen) according to the manufacturer's
protocol. cDNA clones were sequenced commercially (GeneWiz, South Plainfield, NJ) in order to verify the insert identity and for subsequent sequence analysis.

Additional file 1: Germline TRDD and TRDJ gene sequences

Click here for file

[http://www.biomedcentral.com/content/supplementary/1471-2164-11-100-S1.PDF]

Acknowledgements

This work was conducted as part of the Bovine Genome Annotation Sequencing Consortium. The work was supported by USDA-NRI grants nos. 2006-01691 and 2005-01812 and CSREES-USDA Massachusetts Agricultural Experiment Station under project no. MAS00913 and Multi-State NRSP-008 project no. MAS00955.

Author details

¹Department of Veterinary and Animal Sciences, University of Massachusetts, Amherst, MA 01003, USA. ²IMGT, Institut Universitaire de France, Laboratoire d'Immuno Génétique Moléculaire, Université Montpellier 2, UPR CNRS 1142, Montpellier, France. ³Program in Molecular and Cellular Biology, University of Massachusetts, Amherst, MA 01003, USA. ⁴Current address: Department of Epidemiology, Columbia University, New York, NY 10032, USA.

Authors' contributions

CH carried out the annotations, molecular studies and the sequence analyses and helped to draft the manuscript. MPL participated in interpretation of data and helped to draft the manuscript. CB participated in the design of the study, interpretation of data, securing funding and helped to draft the manuscript. All authors read and approved the final manuscript.

Received: 2 June 2009 Accepted: 9 February 2010

Published: 9 February 2010

References

1. Chien YH, Jores R, Crowley MP: Recognition by gamma/delta T cells. *Annu Rev Immunol* 1996, **14**:511-532.

2. Egan PJ, Carding SR: Downmodulation of the inflammatory response to bacterial infection by gammadelta T cells cytotoxic for activated macrophages. *J Exp Med* 2000, **191**(12):2145-2158.

3. Groh V, Steinle A, Bauer S, Spies T: Recognition of stress-induced MHC molecules by intestinal epithelial gammadelta T cells. *Science* 1998, **279**(5357):1737-1740.

4. Havran WL, Chien YH, Allison JP: Recognition of self antigens by skin-derived T cells with invariant gamma delta antigen receptors. *Science* 1991, **252**(5011):1430-1432.

5. Okragly AJ, Hanby-Flarida M, Mann D, Baldwin CL: Bovine gamma/delta T-cell proliferation is associated with self-derived molecules constitutively expressed in vivo on mononuclear phagocytes. *Immunology* 1996, **87**(1):71-79.

6. Sathiyaseelan T, Naiman B, Welte S, Machugh N, Black SJ, Baldwin CL: Immunological characterization of a gammadelta T-cell stimulatory ligand on autologous monocytes. *Immunology* 2002, **105**(2):181-189.

7. Morita CT, Parker CM, Brenner MB, Band H: TCR usage and functional capabilities of human gamma delta T cells at birth. *J Immunol* 1994, **153**(9):3979-3988.

8. Arden B, Clark SP, Kabelitz D, Mak TW: Mouse T-cell receptor variable gene segment families. *Immunogenetics* 1995, **42**(6):501-530.

9. Arden B, Clark SP, Kabelitz D, Mak TW: Human T-cell receptor variable gene segment families. *Immunogenetics* 1995, **42**(6):455-500.

10. Bosc N, Lefranc M-P: The mouse (*Mus musculus*) T cell receptor alpha (TRA) and delta (TRD) variable genes. *Dev Comp Immunol* 2003, **27**(6-7):465-497.

11. Clark SP, Arden B, Kabelitz D, Mak TW: Comparison of human and mouse T-cell receptor variable gene segment subfamilies. *Immunogenetics* 1995, **42**(6):531-540.

12. Giudicelli V, Chaume D, Lefranc M-P: IMGT/GENE-DB: a comprehensive database for human and mouse immunoglobulin and T cell receptor genes. *Nucleic Acids Res* 2005, **33** Database: D256-261.

13. Glusman G, Rowen L, Lee I, Boysen C, Roach JC, Smit AF, Wang K, Koop BF. Hood L: Comparative genomics of the human and mouse T cell receptor loci. *Immunity* 2001, **15**(3):337-349.

14. Lefranc M-P, Lefranc G: The T cell receptor FactsBook. London, UK: Academic Press 2001.

15. Antonacci R, Lanave C, Del Faro L, Vaccarelli G, Ciccarese S, Massari S: Artiodactyl emergence is accompanied by the birth of an extensive pool of diverse germline TRDV1 genes. *Immunogenetics* 2005, **57**(3-4):254-266.

16. Blumerman SL, Herzig CT, Rogers AN, Telfer JC, Baldwin CL: Differential TCR gene usage between WC1- and WC1+ ruminant gammadelta T cell subpopulations including those responding to bacterial antigen. *Immunogenetics* 2006, **58**(8):680-692.

17. Herzig CT, Blumerman SL, Baldwin CL: Identification of three new bovine T-cell receptor delta variable gene subgroups expressed by peripheral blood T cells. *Immunogenetics* 2006, **58**(9):746-757.

18. Ishiguro N, Aida Y, Shinagawa T, Shinagawa M: Molecular structures of cattle T-cell receptor gamma and delta chains predominantly expressed on peripheral blood lymphocytes. *Immunogenetics* 1993, **38**(6):437-443.

19. Takeuchi N, Ishiguro N, Shinagawa M: Molecular cloning and sequence analysis of bovine T-cell receptor gamma and delta chain genes. *Immunogenetics* 1992, **35**(2):89-96.

20. Van Rhijn I, Spiering R, Smits M, van Blokland MT, de Weger R, van Eden W, Rutten VP, Koets AP: Highly diverse TCR delta chain repertoire in bovine tissues due to the use of up to four D segments per delta chain. *Mol Immunol* 2007, **44**(12):3155-3161.

21. Yang YG, Ohta S, Yamada S, Shimizu M, Takagaki Y: Diversity of T cell receptor delta-chain cDNA in the thymus of a one-month-old pig. *J Immunol* 1995, **155**(4):1981-1993.

22. Reinink P, Van Rhijn I: The bovine T cell receptor alpha/delta locus contains over 400 V genes and encodes V genes without CDR2. *Immunogenetics* 2009, **61**(7):541-549.

23. Hein WR, Dudler L: Divergent evolution of T cell repertoires: extensive diversity and developmentally regulated expression of the sheep gamma delta T cell receptor. *EMBO J* 1993, **12**(2):715-724.

24. Oettinger MA, Schatz DG, Gorka C, Baltimore D: RAG-1 and RAG-2, adjacent genes that synergistically activate V(D)J recombination. *Science* 1990, **248**(4962):1517-1523.

25. Hesse JE, Lieber MR, Mizuuchi K, Gellert M: V(D)J recombination: a functional definition of the joining signals. *Genes Dev* 1989, **3**(7):1053-1061.

26. Chien YH, Bonneville M: Gamma delta T cell receptors. *Cell Mol Life Sci* 2006, **63**(18):2089-2094.

27. Hein WR, Dudler L: TCR gamma delta+ cells are prominent in normal bovine skin and express a diverse repertoire of antigen receptors. *Immunology* 1997, **91**(1):58-64.

28. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, et al: Clustal W and Clustal X version 2.0. *Bioinformatics* 2007, **23**(21):2947-2948.

29. Lefranc M-P, Pommie C, Ruiz M, Giudicelli V, Foulquier E, Truong L, Thouvenin-Contet V, Lefranc G: IMGT unique numbering for immunoglobulin and T cell receptor variable domains and Ig superfamily V-like domains. *Dev Comp Immunol* 2003, **27**(1):55-77.

30. Clamp M, Cuff J, Searle SM, Barton GJ: The Jalview Java alignment editor. *Bioinformatics* 2004, **20**(3):426-427.

31. Lefranc M-P, Giudicelli V, Ginestoux C, Jabado-Michaloud J, Folch G, Bellahcene F, Wu Y, Gemrot E, Brochet X, Lane J, et al: IMGT, the international ImMunoGeneTics information system. *Nucleic Acids Res* 2009, **37** Database: D1006-1012.

32. Brochet X, Lefranc M-P, Giudicelli V: IMGT/V-QUEST: the highly customized and integrated system for IG and TR standardized V-J and V-D-J sequence analysis. *Nucleic Acids Res* 2008, **36** Web Server: W5037-508.

33. Ruiz M, Lefranc M-P: IMGT gene identification and Colliers de Perles of human immunoglobulins with known 3D structures. *Immunogenetics* 2002, **53**(10-11):857-883.

34. Tonegawa S: Somatic generation of antibody diversity. *Nature* 1983, **302**(5090):575-581.

35. Conrad ML, Mawer MA, Lefranc M-P, McKinnell L, Whitehead J, Davis SK, Pettman R, Koop BF: The genomic sequence of the bovine T cell receptor
gamma TRG loci and localization of the TRGC5 cassette. *Vet Immunol Immunopathol* 2007, 115(3-4):346-356.

36. Herzig C, Blumerman S, Lefranc M-P, Baldwin C: Bovine T cell receptor gamma variable and constant genes: combinatorial usage by circulating gammadelta T cells. *Immunogenetics* 2006, 58(2-3):138-151.

37. Adams EJ, Chien YH, Garcia KC: Structure of a gammadelta T cell receptor in complex with the nonclassical MHC T22. *Science* 2005, 308(5719):227-231.

38. Shin S, El-Diwany R, Schaffert S, Adams EJ, Garcia KC, Pereira P, Chien YH: Antigen recognition determinants of gammadelta T cell receptors. *Science* 2005, 308(5719):252-255.

39. Adams EJ, Strop P, Shin S, Chien YH, Garcia KC: An autonomous CDR3delta is sufficient for recognition of the nonclassical MHC class I molecules T10 and T22 by gammadelta T cells. *Nat Immunol* 2008, 9(7):777-784.

40. Almagro JC, Martinez L, Smith SL, Alagon A, Estevez J, Paniagua J: Analysis of the horse V(H) repertoire and comparison with the human IGHV germline genes, and sheep, cattle and pig V(H) sequences. *Mol Immunol* 2006, 43(11):1836-1845.

41. Lewis SE, Searle SM, Harris N, Gibson M, Lyer V, Richter J, Wiel C, Bayraktaroglir L, Birney E, Crosby MA, et al: *Apollo: a sequence annotation editor*. *Genome Biol* 2002, 3(12):RESEARCH0082.

42. Elsik CG, Tellam RL, Worley KC, Gibbs RA, Muzny DM, Weinstock GM, Adelson DL, Eichler EE, Elnitski L, Guigo R, et al: The genome sequence of taurine cattle: a window to ruminant biology and evolution. *Science* 2009, 324(5926):522-528.

43. Hall T: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. *Nucleic Acids Symp Ser* 1999, 41:95-98.

44. Florea L, Hartzell G, Zhang Z, Rubin GM, Miller W: A computer program for aligning a cDNA sequence with a genomic DNA sequence. *Genome Res* 1998, 8(9):967-974.

45. Crooks GE, Hon G, Chandonia JM, Brenner SE: WebLogo: a sequence logo generator. *Genome Res* 2004, 14(6):1188-1190.

46. Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. *Mol Biol Evol* 2007, 24(8):1596-1599.

47. Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. *Mol Biol Evol* 1987, 4(4):406-425.

48. Felsenstein J: Confidence limits on phylogenies: An approach using the bootstrap. *Evolution* 1985, 39:783-791.

49. Tamura K, Nei M, Kumar S: Prospects for inferring very large phylogenies by using the neighbor-joining method. *Proc Natl Acad Sci USA* 2004, 101(30):11030-11035.

doi:10.1186/1471-2164-11-100

Cite this article as: Herzig et al.: Annotation and classification of the bovine T cell receptor delta genes. BMC Genomics 2010 11:100.

Submit your next manuscript to BioMed Central
and take full advantage of:

• Convenient online submission

• Thorough peer review

• No space constraints or color figure charges

• Immediate publication on acceptance

• Inclusion in PubMed, CAS, Scopus and Google Scholar

• Research which is freely available for redistribution

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central