Genomic and expression analyses of Tursiops truncatus T cell receptor gamma (TRG) and alpha/delta (TRA/TRD) loci reveal a similar basic public γδ repertoire in dolphin and human
- Giovanna Linguiti†1,
- Rachele Antonacci†1,
- Gianluca Tasco2,
- Francesco Grande3,
- Rita Casadio2,
- Serafina Massari4,
- Vito Castelli1,
- Arianna Consiglio5,
- Marie-Paule Lefranc6 and
- Salvatrice Ciccarese1Email authorView ORCID ID profile
© The Author(s). 2016
Received: 29 January 2016
Accepted: 15 June 2016
Published: 15 August 2016
The Erratum to this article has been published in BMC Genomics 2016 17:778
The bottlenose dolphin (Tursiops truncatus) is a mammal that belongs to the Cetartiodactyla and have lived in marine ecosystems for nearly 60 millions years. Despite its popularity, our knowledge about its adaptive immunity and evolution is very limited. Furthermore, nothing is known about the genomics and evolution of dolphin antigen receptor immunity.
Here we report a evolutionary and expression study of Tursiops truncatus T cell receptor gamma (TRG) and alpha/delta (TRA/TRD) genes. We have identified in silico the TRG and TRA/TRD genes and analyzed the relevant mature transcripts in blood and in skin from four subjects.
The dolphin TRG locus is the smallest and simplest of all mammalian loci as yet studied. It shows a genomic organization comprising two variable (V1 and V2), three joining (J1, J2 and J3) and a single constant (C), genes. Despite the fragmented nature of the genome assemblies, we deduced the TRA/TRD locus organization, with the recent TRDV1 subgroup genes duplications, as it is expected in artiodactyls.
Expression analysis from blood of a subject allowed us to assign unambiguously eight TRAV genes to those annotated in the genomic sequence and to twelve new genes, belonging to five different subgroups. All transcripts were productive and no relevant biases towards TRAV-J rearrangements are observed.
Blood and skin from four unrelated subjects expression data provide evidence for an unusual ratio of productive/unproductive transcripts which arise from the TRG V-J gene rearrangement and for a “public” gamma delta TR repertoire. The productive cDNA sequences, shared both in the same and in different individuals, include biases of the TRGV1 and TRGJ2 genes.
The high frequency of TRGV1-J2/TRDV1- D1-J4 productive rearrangements in dolphins may represent an interesting oligo-clonal population comparable to that found in human with the TRGV9- JP/TRDV2-D-J T cells and in primates.
Although the features of the TRG and TRA/TRD loci organization reflect those of the so far examined artiodactyls, genomic results highlight in dolphin an unusually simple TRG locus. The cDNA analysis reveal productive TRA/TRD transcripts and unusual ratios of productive/unproductive TRG transcripts. Comparing multiple different individuals, evidence is found for a “public” gamma delta TCR repertoire thus suggesting that in dolphins as in human the gamma delta TCR repertoire is accompanied by selection for public gamma chain.
KeywordsT cell receptor TRG locus TRGV TRGJ and TRGC genes TRA/TRD locus TRAV and TRDV genes Dolphin genome Expression analysis IMGT
Bottlenose dolphin (Tursiops truncatus) and the other cetaceans represent the most successful mammalian colonization of the aquatic environment and have undergone a radical transformation from the original mammalian bodyplan. The discovery of two archaic whales with morphological homology between Cetacea and Artiodactyla brought conclusive anatomical support to clade Cetartiodactyla [1, 2]. Whales and hippos shared a common semiaquatic ancestor that branched off from other artiodactyls around 60 million years ago [3–5]. One of the two branches would evolve into cetaceans, possibly beginning about 52 million years ago, with the protowhale Pakicetus, which underwent aquatic adaptation into the completely aquatic cetaceans . So far nothing is known about the genomic organization of dolphin immunoglobulins (IG) and T cell receptor (TR) loci. The only studies of antigen receptors immunity revealed that IgG are present in whales [6, 7] and IGHG and IGHA genes have been described in the Atlantic bottlenose dolphin . Within artiodactyls, the locus organization and expression of TRG and TRA/TRD genes have been characterized in ruminants; these species have been shown to possess a large TRG [9–11] and TRA/TRD [12–14] germline repertoire.
Here we present a evolutionary and expression analysis of Tursiops truncatus TRG and TRA/TRD genes. The surprising feature concerning TRG genes was, on the one hand, that the overall organization of the dolphin TRG locus resembles more the structure of a typical cassette of artiodactyls (IMGT®, the international ImMunoGeneTics information system®, http://www.imgt.org  > Locus representation: Sheep (Ovis aries) TRG1) than the structure typical of the human locus (IMGT® > Locus representation: Human (Homo sapiens) TRG). On the other hand, equally surprising was the finding of an unusual mechanism of biases in the V-J gene rearrangement usage, which is reminiscent of the most frequently used in the human peripheral γδ T cells repertoire of productively rearranged TRGV genes . Despite the fragmented and incomplete nature of the assembly, we have obtained important information on the genomics and the evolution of the TRA/TRD dolphin potential repertoire and its relationship with the expressed chains. Furthermore, the structural 3D visualization, computed by adopting a comparative procedure, using cDNA TRGV-J and TRDV-D-J rearranged amino acid sequences from a single individual, is consistent with the finding that the predicted γδ pairing, present both in the blood and in the skin, is shared among the organisms living in a controlled environment (kept under human care) as well as in those living in marine environment. This finding highlights in dolphin the existence of a basic “public” γδ repertoire of a given TR in a range of public T cell responses.
Genomic arrangement and evolution of the dolphin TRG locus
Genomic arrangement and evolution of the dolphin TRA/TRD genes
5′ RACE PCR and RT-PCR on blood and skin RNA identified the dolphin TRG, TRA and TRD chains repertoire
List of primers used in 5′ RACE, RT and genomic PCR
Location and sequence positions
TRGC EX3b 43364–43381
TRGC EX1b 37807–37825
TRGC EX1b 37635–37654
TRGC EX1b 37692–37712
nested, V1-V2 RT-PCR
TRGV1 L-Part1b 9989–10007
V1 RT-PCR, V1J2 genomic PCR
TRGV2 L-Part1b 21785–21805
V2 RT-PCR, V2J3 genomic PCR
V1J2 genomic PCR
V2J3 genomic PCR
TRAC EX1c 86762–86780
TRAC EX1c 86745–86763
TRAC EX1c 86574–86592
TRDC EX1d 85917–85934
TRDC EX1d 85938–85957
TRDC EX2d 86773–86791
TRDV4 EXd 91850–91871
V4 genomic PCR
TRDV4 EXd 92062–92081
V4 genomic PCR
Summary of the different 5′RACE and RT-PCR experiments and the obtained rearrangement types
FWD primer name
REV primer name
Total number of non-redundant clonotypes
Number of non-redundant out-of-frame clonotypes
Number of non-redundant in-frame clonotypes
Non-redundant in-frame clonotypes by rearrangement type
GenBank (GEDI) accession numbers
1 TRDV1*01–TRDJ4*01 c
To investigate the dolphin TRA chain repertoire, total RNA from the peripheral blood of a female dolphin (identified as L) was used as template in the single 5′ RACE experiment (Tables 1 and 2). A total of 41 different TRA clonotypes were obtained and sequenced (Fig. 3). All sequences were productive (in-frame junction and no stop codon), and the leader region was of 17 to 20 amino acids depending on the V subgroup. The CDR1- and CDR2-IMGT lengths of the transcripts [6.4.], [6.8.], [7.8.] corresponded to nine different TRAV subgroups and 29 different genes and were associated with diverse CDR3-IMGT of various length from 8 to 16 AA. In our cDNA collection, 8 TRAV genes (TRAV16, TRAV8-1, TRAV18-1, TRAV38-1, TRAV14, TRAV20-1D, TRAV9, TRAV17) were assigned unambiguously to genes annotated in the genomic sequence (Fig. 2), while 12 could be assigned to new genes, belonging to five different subgroups. One gene belongs to a new subgroup, TRAV13 (TRAV13S1), not yet identified in the dolphin genomic sequence. One gene belongs to subgroup TRDV1 (TRDV1S2), as shown by its CDR1- and CDR2-IMGT lengths [7.3.], and demonstrates that dolphin TRDV genes can, as in other species, participate to the synthesis of TRA chains by rearranging to a TRAJ gene (here, TRAJ11) . Among the other new genes, six belong to subgroup TRAV45 (TRAV18S2, TRAV18S3, TRAV18S4, TRAV18S5, TRAV18S6, TRAV18S7), three belong to subgroup TRAV20 (TRAV20S2, TRAV20S3, TRAV20S4) and one to subgroup 42 (TRAV8S2). As these last subgroups have several members, an IMGT approved provisional nomenclature was assigned (with the letter S), allowing these genes to be entered in IMGT/GENE-DB and IMGT® tools (IMGT/V-QUEST and IMGT/HighV-QUEST)  while waiting for the identification and location of these genes in the reference genomic sequence. Three 5′ RACE experiments on total RNA isolated from the peripheral blood of two unrelated adult animals (identified as M and L) (Tables 1 and 2) were carried out to investigate the dolphin TRD chain repertoire; only one of these three PCR amplifications produced 3 in-frame, 1 out-of-frame and 1 sterile germline, clonotypes from the animal identified as M (Fig. 3b).
Potential TRGV domain repertoire of productive and unproductive trancripts
Analyzing the TRG in-frame transcripts it is noteworthy that 5 TRG clonotypes were found identical in two or even three different individuals L and C (5RV1L2*/C2/C5), M, C and K (RTV1M1/C1/K2/K3), K and M (RTV1K7/5RV1M1), and M and C (RTV2M4/C7 and RTV2M5/C8) (Additional file 10). This observation was rather intriguing as they represented together 14/39 in-frame sequences whereas in contrast each out-of-frame clonotype was found in a single individual. These shared clonotypes result from V1-J2 rearrangements in L and C (CDR-IMGT lengths [8.7.13]) and in M, C and K (CDR-IMGT lengths [8.7.14]), from V1-J3 rearrangements in K and M (CDR-IMGT lengths [8.7.14]) and from V2-J3 rearrangements in M and C (CDR-IMGT lengths [8.6.14] (Fig. 4 and Additional file 10). This description of shared T cell clonoypes correspond to what is known in the literature as “public T cell response” in which T cells bearing identical TR may respond to the same antigenic epitope in different individuals . Although the number of the germline TRG genes is low, which implies a reduced potential in the V-J recombination, a sufficient diversity and variability of the TR gamma transcripts seems to be guaranteed in the dolphin by the classical process of CDR3 diversity formation during somatic rearrangement . Indeed, the creation of the CDR3 diversity results from the trimming of the 3′V-REGION (up to 12 nucleotides (nt) for the in-frame junctions, up to 17 for the out-of-frame junctions), from the trimming of the 5′J-REGION (up to 14 nt for the in-frame junctions, up to 22 for the out-of-frame junctions), and from the addition at random of the N nucleotides creating the N-REGION (up to 16 nt for the in-frame junctions, up to 23 for the out-of-frame junctions) (Fig. 4). This junction diversity is due to the activity of the terminal deoxynucleotidyl transferase (TdT) encoded by DNTT. The gene (NCBI ID: 101323636) has been identified in the bottle nosed dolphin genome and its amino acid sequence is 84 % identical to the human DNTT. The graphical representation of the number of in-frame versus out-of-frame sequences obtained for the 6 possible TRG rearrangements V1-J1, V1-J2, V1-J3, V2-J1, V2-J2 and V2-J3 display striking differences (Additional file 11). Both tests (Chi-squared p-value confirmed with Fisher’s p-value) reject the null hypothesis for V1-J2 and V2-J1 (Additional file 12). This result confirms what was noticed at first sight and it follows that V1-J2 gene rearrangements were dominant among the in-frame transcripts and were rare among the out-of-frame transcripts.
Computational analyses predict the pairing of the TRGV1-J2 and of the TRDV1-D1-J4 variable domains
In this study we report an extensive analysis of the genomic organization and expression of the TRG and TRA/TRD genes in dolphin. According to comparative analyses, dolphin TRG locus is the simplest and the smallest among the mammalian TRG loci identified to date [19, 20] and its organization is reminiscent of the structure of a typical single cassette of artiodactyls [9, 15] with a small number of genes, i.e. two TRGV, three TRGJ and one TRGC (Fig. 1).
The analysis of dolphin TRA/TRD locus confirmed that TRD genes are clustered within the TRA locus and that genes belonging to the TRDV1 subgroup are distributed among the TRAV genes as it is commonly expected in artiodactyls TRA/TRD locus [13, 14, 26, 32]. A total of 16 TRAV and 5 TRDV genes have been identified (Fig. 2). By the criterion that gene sequences having 75 % or greater nucleotide identity belong to the same subgroup, the TRAV and the TRDV genes belong to 13 and to three subgroups, respectively (Additional file 7). The sheep TRDV1 subgroup has been estimated to contain at least 40 genes , while only 25 TRDV1 genes have been identified in the genomic assembly . The phylogenetic analysis assigns the membership of the dolphin TRDV1 genes due to the monophyletic groupings marked by 25 sheep, 6 dromedary and 3 dolphin members in contrast with the single human one (Additional file 7B).
Dolphin TR alpha chain expression analysis allowed us to identify new TRAV genes, with respect to the available genomic sequence. Furthermore a bias towards rearrangements containing TRA genes belonging to the TRAV18 (12/40 cDNA) and TRAV20 (11/40 cDNA) gene subgroups, was observed (Fig. 3a). On the contrary, the usage of the 61 TRAJ genes is generally random with a slight increase in usage of TRAJ (Fig. 2) between 54 and 22 (31 of 50 functional rearrangements) (Additional file 9); this finding being consistent with the widely accepted view that TRAV-TRAJ recombination proceeds in a coordinated, sequential manner from proximal to progressively more distal TRAV and TRAJ genes [33, 34].
Dolphin TR gamma chain expression analysis demonstrated that the two TRGV and three TRGJ were used in every possible combination, although a bias towards some transcripts (TRGV1-TRGJ2 and TRGV2-TRGJ3) was noted. Furthermore, about half the transcripts using TRGV2 were unproductive due to the presence of stop codons in CDR3. The percentage values of the productive/unproductive rearrangements are similar for both cDNA (Fig. 4) and genomic clones (Fig. 5), in contrast with what is usually obseved (percentage of unproductive rearrangements lower in cDNA, due to nonsense-mediated decay of RNA).
In a previous work , it was reported that biased V-J gene rearrangement contributes to the regulation of the mature TRG repertoire. The biases in a given TR repertoire can stem from properties of the gene rearrangement process, as well as from thymic selection and the expansion of T cell clones. In the present work, we can make the following considerations: i) it seems to be a double preferentiality and that of the gene TRGV1 with respect to the gene TRGV2 as well as of the TRG V1-J2 rearrangement with respect to the five others (Fig. 4), the latter being supported given the comparison between the frequency of the in-frame and out-of-frame rearrangements both in cDNA and in genomic DNA (Additional file 11); ii) the fact that unrelated subjects show not only a biased usage of V-J genes, but also a biased number of nucleotides inserted/deleted at junction regions (Fig. 4 and Fig. 5c), could be explained by the presence of common antigens which can stimulate and expand T cells with a particular type of gamma chain, suggesting the existence of a basic “public” repertoire of a given TR in a range of public T cell responses; iii) finally we propose that the occurrence of clonotypes shared by different individuals who live both in marine and in artificial marine “habitat”, described as “convergent recombination” , could be strictly related to the biased V-J recombinational event.
The mechanisms that determine biases in genes use remain unclear. In a recent paper  a physical model of chromatin conformation at the TRB D-J genomic locus explains more than 80 % of the biases in TRBJ use that was measured in murine T cells. As a consequence of these structural and other biases, TR sequences are produced with different a priori frequencies, thus affecting their probability of becoming public TR that are shared among individuals. In dolphin, we could explain the abundance of TRGV1-J2 repertoire among individuals hypothizing that this combination could be produced by the rearrangement process with different a priori probabilities because an expanded role of chromatin conformation in TRGV-J rearrangement, which controls both the gene accessibility and the precise determination of gene use.
An evolutionary correlation between the dolphin TRGV1 and the human TRGV9 (Additional file 4A) genes and the dolphin TRGJ2 and the human TRGJP (Additional file 15) genes seems to exist, as in these two species the same mechanism pushes to an accurate determination of the J gene usage. In fact, dolphin TRGJ2 this work) and human TRGJP, are the most frequently used J genes in the peripheral γδ T cells  and occupy an intermediate position with respect to the other two J genes. At present we have knowledge of the position of the genes on the physical map for human (IMGT), dromedary , dolphin (this work) and sheep  (Additional file 1) and cattle  TRG loci.
It is admitted that the expressed γ/δ T cell repertoire partly depends upon preferentially rearranged TRGV-J gene combinations, indeed in human the gamma delta TCR repertoire is accompanied by selection for public gamma chain sequences such that many unrelated individuals overlap extensive in their circulating repertoire . As a conseguence, the high frequency of TRGV1-J2/TRDV1-D1-J4 productive rearrangements in dolphins may represent a situation of oligoclonality comparable to that found in human with TRGV9-JP/TRDV2-D-J T cells, and in primates.
The similarity in dolphin and human of a basic public γδ repertoire, seems to be correlated with other recent findings. McGowen discovered several genes, potentially under positive selection in the dolphin lineage, associated with the nervous system, including those related to human intellectual disabilities, synaptic plasticity and sleep . Moreover bottlenose dolphins are the only animals with man and apes, to be able to recognize themselves when confronted with a mirror , and have demonstrated the numerical skills . While here, in the present work, the functional convergence of γδ domains is suggested among mammals, recently it was proposed similarity of dual-function TRA and TRD genes in jawed vertebrates and in the VLRA and VLRC genes in jawless vertebrates and their differential expression in two major T cell lineages [41–43]. Therefore comparative immunobiology of different vertebrate lineages may reveal heretofore unrealized features.
The present study identifies the genomic organization and the gene content of the TRG and the TRA/TRD loci in the high quality draft sequence of the bottlenose dolphin (Tursiops truncatus) genome. The genomic structure of the smallest TRG locus thus described in mammals, includes two TRGV, three TRGV and only one TRGC genes. Through phylogenetic and expression analyses, 8 TRAV were assigned unambiguously to genes annotated in the TRA/TRD locus genomic sequence, while 12 TRAV could be assigned to new genes, belonging to five different subgroups. The presence of several variable genes belonging to the TRDV1 subgroup, makes the TRA/TRD dolphin locus more similar to the TRA/TRD locus of artiodactyls than to the human locus.
By comparing multiple different individuals, we provide evidence of an unusual ratio of productive/unproductive TRG transcripts and of a bias towards TRGV1-TRGJ2 rearrangements, which were dominant among the in-frame transcripts and were rare among the out-of-frame transcripts. Moreover, the cDNA analysis revealed sharing of in-frame TRG sequences within the same and in different individuals living in a controlled environment as well as in marine environment, suggesting expansion of “public” TCR by a common antigen. The selection for public gamma chain and the high frequency of TRGV1-J2/TRDV1-D1-J4 productive rearrangements in dolphins may represent a situation comparable to that found in human with TRGV9-JP/TRDV2-D-J T cells.
Genome and sequence analysis
The bottlenose dolphin (Tursiops truncatus) genome is being sequenced at ~2X coverage (BioProject: PRJNA20367) by the Human Genome Sequencing Center at the Baylor College of Medicine and the Broad Institute using a whole genome shotgun sequencing strategy . In 2008, Ensembl released the first low-coverage 2.59× assembly of the dolphin (turTru1). We employed these genome assemblies using BLAST algorithm to identify the TRG and TRA/TRD loci in this species.
For the TRG locus, two overlapping scaffolds were retrieved (GEDI ID: JH473572.1; BCM-HGSC ID: contig 425448–578749), respectively of 96017 and 284974 bp (gaps included). A sequence of 188414 bp was analysed. Amphiphysin (AMPH) and related to steroidogenic acute regulatory protein D3-N-terminal like (STARD3NL), flanking TRG locus at 5′ and 3′ ends, respectively, were included in the analysis. They were predicted to be functional in dolphin (GenBank ID: XM_004317564.1; Ensembl ID: ENSTTRT00000004099). The TRG genes were identified using both our dolphin cDNA collection (this work) and the corresponding human (GEDI ID: AF159056) and sheep (GEDI ID: DQ992075.1, DQ992074.1) genomic sequences. Locations of the TRG genes are provided in Fig. 1b.
For the TRA/TRD locus, we retrieved a sequence of 482052 bp from two GenBank scaffolds, JH484271.1 and JH481615.1, and five EMBL-EBI scaffolds, Ens_742, Ens_97, Ens_89, Ens_123 and Ens_112178. Scaffold Ens_97 and Ens_123 overlap for about 16,7 Kb, including TRA14/DV4, TRA9, TRA16 and TRA17 genes, while scaffold Ens_89 and JH484271.1 overlap for about 10 Kb, a region that includes two genes, TRAV1S1 and TRAV38.1. The TRA/TRD genes were identified using the corresponding human (GEDI ID: AE000521.1) genomic sequences. Sequences of all TRA/TRD genes are in Additional file 5. Computational analysis of the dolphin TR loci was conducted using the following programs: RepeatMasker for the identification of genome-wide repeats and low complexity regions  (RepeatMaskerathttp://www.repeatmasker.org) and Pipmaker  (http://www.pipmaker.bx.psu.edu/pipmaker/) for the alignment of the dolphin sequence with the human counterpart. ClustalW (http://www.ebi.ac.uk/Tools/msa/clustalw2/) and IMGT/V-QUEST (http://www.imgt.org/IMGT_vquest/share/textes/) tools allowed the identification and characterization of the TR genes.
The TRGV, TRGJ, TRAV and TRDV genes used for the phylogenetic analyses were retrieved from IMGT/LIGM-DB and GenBank databases with the following accession numbers: AF159056 (human TRG locus), DQ992075 (sheep TRG1 locus), DQ992074 (sheep TRG2 locus), JN165102 (dromedary TRGV1), JN172913 (dromedary TRGV1), AE000521.1 (human TRA/TRD locus); sheep TRA/TRD accession numbers  and FN298219- FN298227 (dromedary TRD genes) . Multiple alignments of the sequences under analysis were carried out with the MUSCLE program . Phylogenetic analyses were performed using MEGA version 6.06  and the bootstrap consensus tree inferred from 1000 replications using the Neighbor-Joining method [49, 50].
Animals (source of tissue)
Blood samples were provided by Zoomarine Italia S.p.A. (Rome, Italy) and were collected from three dolphins, two males (Marco and King) and one female (Leah). The three individuals were born and kept under human care and are unrelated. In particular, Marco was born in the dolphinarium in Bruges (Belgium) and Tex, Marco’s father, is from the United States (Texas, Gulf of Mexico). King was born in the dolphinarium in Albufeira (Portugal), and Sam, King’s father had Cuban origins. Leah was born in the dolphinarium in Benidorm (Spain) and Eduardo, Leah’s father has Cuban origins. The identifying letters are M, K and L, respectively. The Bank for the Tissues of Mediterranean Marine Mammals (Padua, Italy) provided us a sample of skin (epidermis plus dermis) belonging to a wild dolphin, that was found beached in the Northern Adriatic Sea; for this animal the identified letter is C.
5′ RACE and RT-PCR
Four types of 5′ RACE and three types of RT-PCR (total of six and six experiments, respectively) on total RNA from the peripheral blood of three unrelated adult animals (identified as M, K and L) and from the skin of animal (identified as C) (Table 1 and 2) were carried out to investigate the dolphin TRG, TRA and TRD chains repertoire. Two 5′RACE experiments from the peripheral blood of the animals (identified as M and L) and three types of RT-PCR, two from blood (K and M) and one from skin (C), were carried out to investigate the dolphin TRG chain repertoire. A single 5′RACE experiment from the peripheral blood of the animal identified as L was carried out to investigate the dolphin TRA chain repertoire. Three 5′RACE experiments from the peripheral blood of the animals (identified as M and L) were carried out to investigate the dolphin TRD chain repertoire.
Total RNA was isolated from peripheral blood leukocytes (PBL) or skin using the Trizol method according to the manufacturer’s protocol (Invitrogen, Carlsbad, CA), and integrity of RNA was verified on a 1 % agarose gel. About 5 μg of total RNA were reverse transcribed with Superscript II (Invitrogen, Carlsbad, CA) by using specific primers (Table 1), designed on the sequences of the first exon for each dolphin TR constant gene sequence (TC3L for gamma chain, TA1C1L for alpha chain and TD2CL for delta chain). After linking a poly-C tail at the 5′end of the cDNAss, the cDNAds was performed with Platinum Taq Polymerase (Invitrogen) by using specific primers as lower primers, TC1L1 for gamma chain, TA1C2L for alpha chain and TD1C2L for delta chain (Table 1) and an anchor oligonucleotide as upper primer (AAP) provided from the supplier (Invitrogen). PCR conditions were the following: one cycle at 94 °C for 1 min; 35 cycles at 94 °C for 30 s, 58 °C for 45 s, 72 °C for 1 min; a final cycle of 30 min at 72 °C. The products were then amplified in a subsequent nested PCR experiment by using specific lower primers, TC1L2 for gamma chain, TA1C3L for alpha chain and TD1CL1 for delta chain (Table 1) and AUAP oligonucleotide as upper primer, provided from the supplier (Invitrogen). Nested PCR conditions were the following: one cycle at 94 °C for 1 min; 30 cycles at 94 °C for 30 s, 58 °C for 35 s, 72 °C for 30 s; a final cycle of 30 min at 72 °C.RT-PCR experiments were carried out amplifing rearranged transcripts containing TRGV1 and TRGV2 genes. Upper primers containing TRGV1 (TV1LU) and TRGV2 (TV7LU) sequences, and lower primer containing the I exon of TRGC (TC1L2) sequence were used on sscDNA (Table 1 and 2). RT-PCR conditions were: one cycle at 94 °C for 2,30 min; 35 cycles at 94 °C for 30 s, 58 °C for 40 s, 72 °C for 40 s; a final cycle of 30 min at 72 °C. The RT-PCR and RACE products were then gel-purified and cloned using StrataClone PCR Cloning Kit (Statagene). Random selected positive clones for each cloning were sequenced by a commercial service. cDNA sequence data were processed and analyzed using the Blast program (http://www.blast.ncbi.nlm.nih.gov/Blast.cgi), Clustal W2 (http://www.ebi.ac.uk/Tools/msa/clustalw2/) and IMGT_ tools (IMGT/V-QUEST) [51, 52] with integrated IMGT/JunctionAnalysis tools [53, 54] and the IMGT unique numbering for V domain  (http://www.imgt.org/).
Genomic DNA isolation and PCR
Genomic DNA was extracted from whole blood of a female subject (animal identifiant letter L), with a salting-out method  with two modifications. First, whole blood was mixed with erythrocyte lysis buffer (155 mM NH4Cl, 10 mM KHCO3, 1 mM EDTA, pH 7.4) before the harvested white cell pellet was mixed with nucleus lysis buffer as described . Second, incubation with proteinase K was carried out for 2 h at 56 °C, instead of overnight at 37 °C. The quality of the genomic DNA was evaluated by agarose gel electrophoresis and concentration determined by 260 nm absorbance measurements. Genomic PCR was performed with 50 ng to 100 ng of genomic DNA as template using specific upper primers (TV1L1 and TV7LU) designed on the two TRGV (TRGV1 and TRGV2) gene sequences in combination with two lower primers (J2GL and J5BR) designed on the two TRGJ (TRGJ2 and TRGJ3) gene sequences (Table 1). Two genomic PCR were performed to amplify TRGV1-TRGJ2 and TRGV2-TRGJ3 rearrangement combinations, respectively. High-fidelity polymerase was used to minimize possible PCR errors. PCR were performed following the manifacture’s instruction for the DNA polymerase (Platinum®Taq DNA Polymerase, Life Technologies). V1-J2 genomic PCR conditions were the following: one cycle at 94 °C for 3 min; 35 cycles at 94 °C for 30 s, 62 °C for 30 s, 72 °C for 30 s; a final cycle of 30 min at 72 °C. V2-J3 genomic PCR conditions were the following: one cycle at 94 °C for 3 min; 35 cycles at 94 °C for 30 s, 60 °C for 30 s, 72 °C for 30 s; a final cycle of 30 min at 72 °C. The obtained fragments were agarose gel purified, cloned using StrataClone PCR Cloning Kit (Stratagene) and sequenced by a commercial service. Two genomic PCR were performed to amplify TRDV4 gene using a pair of primers designed based on the relative V-exon sequence (Table 1). PCR was performed following the manufacture’s instruction for the MyTaq™ HS DNA Polymerase, (Bioline). The following settings were used: 94 °C for 2 min, followed by 30 cycles each comprising a denaturation step at 94 °C for 30 s, an annealing step of 30 s at 55 °C (according to the melting temperature of the primers), an extension step at 72 °C for 30 s, and a final extension period of 7 min at 72 °C.
Statistical analyses were performed using 2 × 2 contingency tables. All the p-values shown in the Results were obtained using the Chi-squared test, considering as statistically significant a p- value <0.05. Fisher’s Exact test was used to confirm the significance of the Chi-squared test when the counts of observed samples had values <5. When performing multiple comparisons among in-frame and out-of-frame TRG cDNA (Additional file 12), the Chi-squared test p-values were adjusted using Benjamini–Hochberg false discovery rate . All the analyses were performed using the R software environment for statistical computing (https://www.r-project.org/).
Global alignments in protein secondary structure prediction and 3D visualization
Global alignment of the target and template sequences was performed with ClustalW (http://www.ebi.ac.uk/clustalw/index.html) . Furthermore, when necessary, alignment was manually adjusted after predicting the secondary structure of the target and aligned to that of the template as derived with the DSSP program . The secondary structure prediction was computed with SECPRED (http://gpcr.biocomp.unibo.it/cgi/predictors/s/pred_seccgi.cgi) and PSIPRED  and the target/template alignments were computed with YAP (http://gpcr.biocomp.unibo.it/cgi/predictors/alignss/alignss.cgi), that allows to align both primary and secondary structure at the same time. The template was selected from the Protein Data Bank (PDB) on the basis of sequence/function similarity with the target sequence and was the human γδ T cell receptor solved with an atomic solution of 3A°(PDB code and IMGT/3Dstructure-DB: 1hxm) [29, 30]. Target/ template alignments were then fed into Modeller version 9.8 . For a given alignment, 50 3D models were routinely built and, then, evaluated and validated with the PROCHECK  and PROSA2003  suites of programs. Models with the best stereochemical and energetics features were retained. 3D visualization (Additional file 14) of the RTV1M1/5R1D8 and of 5RV1M1/5R1D15 clones was computed, adopting as template the human γδ T cell receptor. The solvent accessibility was computed with DSSP program . The protein complex interface were computed by the online tool PDBePISA at the EBI server (http://www.ebi.ac.uk/msd-srv/prot_int/) and visualized by UCSF Chimera tool (http://www.cgl.ucsf.edu/chimera/). The IMGT Collier de Perles of RTV1M1, 5R1D8, 5RV1M1 and 5R1D15 cDNA clonotypes were obtained using the IMGT/Collier-de-Perles tool (http://www.imgt.org) , starting from amino acid sequences.
CDR, complementarity determining region; FR, framework region; IG, immunoglobulins; T cell receptor gamma locus; TR, T cell receptor; TRA/TRD locus, T cell receptor alpha/delta locus; TRG locus, TRGC, T cell receptor gamma constant; TRGJ, T cell receptor gamma joining; TRGV, T cell receptor gamma variable. All TR genes (functional, ORF, pseudogenes) reported here have been approved by the IMGT/WHO- IUIS nomenclature committee and their designations are in accord with the IMGT nomenclature for human (IMGT®, the international ImMunoGeneTics information system®, http://www.imgt.org)
The financial support of the University of Bari and of the University of Salento is gratefully acknowledged. A.C. is supported by Progetto MICROMAP PON01_02589". Thanks are due to COST Action TD1101 and COST Action BM1405 delivered to R.C. We thank Stefania Brandini and Angela Pala for technical assistance in cDNA cloning experiments and Prof. P. Barsanti for critically reading of the manuscript. We thank all the staff of Zoomarine Italy. We are especially thankful to all the trainers who dedicate themselves every day to dolphins. Through their work, they make possible the study of medical behaviours. We thank the members of the IMGT® team for their expertise and constant motivation.
Università degli Studi di Bari Aldo Moro (ex 60 %).
GL, RA, and SC designed research; GL, GT, VC, and AC performed research; FG contributed new reagents/analytic tools; RA, RC, SM, M-PL. and SC analyzed data; M-PL improved the manuscript; GL and SC wrote the paper. All authors have read and approved the final manuscript.
The author declares that he/she has no competing interests.
Consent for publication
Ethics approval and consent to participate
Zoomarine is a seaside water park. The availability of blood of this study was a byproduct of standard health checks. The blood samples were taken from the caudal vein on the ventral surface of the caudal fin. Training of medical behaviours brings the animals to collaborate completely, assuming and maintaining the positions that allow the veterinarian, constantly assisted by the trainers, to perform the necessary procedures to monitor their welfare. Zoomarine Italia SpA., via Casablanca 61, 00071 Pomezia (RM), Italy https://www.zoomarine.it/.
TR cDNA and genomic sequences were submitted to the GenBank database. TRG cDNAs are under accession numbers: JF755948-JF755968 and JN011996; TRG genomic sequences are under accession numbers: LN886662 - LN886693. TRA cDNAs are under accession numbers: LN610706 - LN610746; TRD cDNAs are under accession numbers: LN610747 - LN610749. Phylogenetic trees presented in the Additional file 4A and B, in the Additional file 7A and B and in the Additional file 15 were deposited in TreeBase (http://purl.org/phylo/treebase/phylows/study/TB2:S19396).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Gingerich PD, Haq M, Zalmout IS, Khan IH, Malkani MS. Origin of whales from early artiodactyls: hands and feet of Eocene Protocetidae from Pakistan. Science. 2001;293(5538):2239–42.View ArticlePubMedGoogle Scholar
- Price SA, Bininda-Emonds OR, Gittleman JL. A complete phylogeny of the whales, dolphins and even-toed hoofed mammals (Cetartiodactyla). Biol Rev Camb Philos Soc. 2005;80(3):445–73.View ArticlePubMedGoogle Scholar
- Boisserie JR, Lihoreau F, Brunet M. The position of Hippopotamidae within Cetartiodactyla. Proc Natl Acad Sci U S A. 2005;102(5):1537–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Bininda-Emonds OR, Cardillo M, Jones KE, MacPhee RD, Beck RM, Grenyer R, Price SA, Vos RA, Gittleman JL, Purvis A. The delayed rise of present-day mammals. Nature. 2007;446(7135):507–12.View ArticlePubMedGoogle Scholar
- Geisler JH, Theodor JM. Hippopotamus and whale phylogeny. Nature. 2009;458(7236):E1–4. discussion E5.View ArticlePubMedGoogle Scholar
- Andresdottir V, Magnadottir B, Andresson OS, Petursson G. Subclasses of IgG from whales. Dev Comp Immunol. 1987;11(4):801–6.View ArticlePubMedGoogle Scholar
- Mancia A, Romano TA, Gefroh HA, Chapman RW, Middleton DL, Warr GW, Lundqvist ML. The Immunoglobulin G Heavy Chain (IGHG) genes of the Atlantic bottlenose dolphin, Tursiops truncatus. Comp biochem physiol Part B, Biochem mol biol. 2006;144(1):38–46.View ArticleGoogle Scholar
- Mancia A, Romano TA, Gefroh HA, Chapman RW, Middleton DL, Warr GW, Lundqvist ML. Characterization of the immunoglobulin A heavy chain gene of the Atlantic bottlenose dolphin (Tursiops truncatus). Vet Immunol Immunopathol. 2007;118(3–4):304–9.View ArticlePubMedGoogle Scholar
- Vaccarelli G, Miccoli MC, Antonacci R, Pesole G, Ciccarese S. Genomic organization and recombinational unit duplication-driven evolution of ovine and bovine T cell receptor gamma loci. BMC Genomics. 2008;9:81.View ArticlePubMedPubMed CentralGoogle Scholar
- Conrad ML, Mawer MA, Lefranc MP, 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–56.View ArticlePubMedGoogle Scholar
- Ciccarese S, Vaccarelli G, Lefranc MP, Tasco G, Consiglio A, Casadio R, Linguiti G, Antonacci R. Characteristics of the somatic hypermutation in the Camelus dromedarius T cell receptor gamma (TRG) and delta (TRD) variable domains. Dev Comp Immunol. 2014;46(2):300–13.View ArticlePubMedGoogle Scholar
- 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–66.View ArticlePubMedGoogle Scholar
- Connelley TK, Degnan K, Longhi CW, Morrison WI. Genomic analysis offers insights into the evolution of the bovine TRA/TRD locus. BMC Genomics. 2014;15:994.View ArticlePubMedPubMed CentralGoogle Scholar
- Piccinni B, Massari S, Caputi Jambrenghi A, Giannico F, Lefranc MP, Ciccarese S, Antonacci R. Sheep (Ovis aries) T cell receptor alpha (TRA) and delta (TRD) genes and genomic organization of the TRA/TRD locus. BMC Genomics. 2015;16:709.View ArticlePubMedPubMed CentralGoogle Scholar
- Lefranc MP, Giudicelli V, Duroux P, Jabado-Michaloud J, Folch G, Aouinti S, Carillon E, Duvergey H, Houles A, Paysan-Lafosse T, et al. IMGT(R), the international ImMunoGeneTics information system(R) 25 years on. Nucleic Acids Res. 2015;43(Database issue):D413–22.View ArticlePubMedGoogle Scholar
- Lefranc MP, Lefranc G. The T cell receptor fact book. 2001.Google Scholar
- Lindblad-Toh K, Garber M, Zuk O, Lin MF, Parker BJ, Washietl S, Kheradpour P, Ernst J, Jordan G, Mauceli E, et al. A high-resolution map of human evolutionary constraint using 29 mammals. Nature. 2011;478(7370):476–82.View ArticlePubMedPubMed CentralGoogle Scholar
- 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–49.View ArticlePubMedGoogle Scholar
- Massari S, Bellahcene F, Vaccarelli G, Carelli G, Mineccia M, Lefranc MP, Antonacci R, Ciccarese S. The deduced structure of the T cell receptor gamma locus in Canis lupus familiaris. Mol Immunol. 2009;46(13):2728–36.View ArticlePubMedGoogle Scholar
- Massari S, Ciccarese S, Antonacci R. Structural and comparative analysis of the T cell receptor gamma (TRG) locus in Oryctolagus cuniculus. Immunogenetics. 2012;64(10):773–9.View ArticlePubMedGoogle Scholar
- Herzig C, Blumerman S, Lefranc MP, Baldwin C. Bovine T cell receptor gamma variable and constant genes: combinatorial usage by circulating gammadelta T cells. Immunogenetics. 2006;58(2–3):138–51.View ArticlePubMedGoogle Scholar
- Antonacci R, Vaccarelli G, Di Meo GP, Piccinni B, Miccoli MC, Cribiu EP, Perucatti A, Iannuzzi L, Ciccarese S. Molecular in situ hybridization analysis of sheep and goat BAC clones identifies the transcriptional orientation of T cell receptor gamma genes on chromosome 4 in bovids. Vet Res Commun. 2007;31(8):977–83.View ArticlePubMedGoogle Scholar
- Vaccarelli G, Antonacci R, Tasco G, Yang F, Giordano L, El Ashmaoui HM, Hassanane MS, Massari S, Casadio R, Ciccarese S. Generation of diversity by somatic mutation in the Camelus dromedarius T-cell receptor gamma variable domains. Eur J Immunol. 2012;42(12):3416–28.View ArticlePubMedGoogle Scholar
- Buresi C, Ghanem N, Huck S, Lefranc G, Lefranc MP. Exon Duplication and Triplication in the Human T-Cell Receptor-Gamma Constant Region Genes. Cytogenet Cell Genet. 1989;51(1–4):973.Google Scholar
- Lefranc MP, 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.View ArticlePubMedGoogle Scholar
- Herzig CTA, Lefranc MP, Baldwin CL. Annotation and classification of the bovine T cell receptor delta genes. BMC Genomics. 2010;11:100.View ArticlePubMedPubMed CentralGoogle Scholar
- Li S, Lefranc MP, Miles JJ, Alamyar E, Giudicelli V, Duroux P, Freeman JD, Corbin VD, Scheerlinck JP, Frohman MA, et al. IMGT/HighV QUEST paradigm for T cell receptor IMGT clonotype diversity and next generation repertoire immunoprofiling. Nat Commun. 2013;4:2333.PubMedPubMed CentralGoogle Scholar
- Venturi V, Kedzierska K, Tanaka MM, Turner SJ, Doherty PC, Davenport MP. Method for assessing the similarity between subsets of the T cell receptor repertoire. J Immunol Methods. 2008;329(1–2):67–80.View ArticlePubMedGoogle Scholar
- Xu B, Pizarro JC, Holmes MA, McBeth C, Groh V, Spies T, Strong RK. Crystal structure of a gammadelta T-cell receptor specific for the human MHC class I homolog MICA. Proc Natl Acad Sci U S A. 2011;108(6):2414–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Ehrenmann F, Kaas Q, Lefranc MP. IMGT/3Dstructure-DB and IMGT/DomainGapAlign: a database and a tool for immunoglobulins or antibodies, T cell receptors, MHC, IgSF and MhcSF. Nucleic Acids Res. 2010;38:D301–7.View ArticlePubMedGoogle Scholar
- Krissinel E, Henrick K. Inference of macromolecular assemblies from crystalline state. J Mol Biol. 2007;372(3):774–97.View ArticlePubMedGoogle Scholar
- Bovine Genome S, Analysis C, Elsik CG, Tellam RL, Worley KC, Gibbs RA, Muzny DM, Weinstock GM, Adelson DL, Eichler EE, et al. The genome sequence of taurine cattle: a window to ruminant biology and evolution. Science. 2009;324(5926):522–8.View ArticleGoogle Scholar
- Jouvin-Marche E, Fuschiotti P, Marche PN. Dynamic aspects of TCRalpha gene recombination: qualitative and quantitative assessments of the TCRalpha chain repertoire in man and mouse. Adv Exp Med Biol. 2009;650:82–92.View ArticlePubMedGoogle Scholar
- Krangel MS. Mechanics of T cell receptor gene rearrangement. Curr Opin Immunol. 2009;21(2):133–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Kohsaka H, Chen PP, Taniguchi A, Ollier WE, Carson DA. Regulation of the mature human T cell receptor gamma repertoire by biased V-J gene rearrangement. J Clin Invest. 1993;91(1):171–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Ndifon W, Gal H, Shifrut E, Aharoni R, Yissachar N, Waysbort N, Reich-Zeliger S, Arnon R, Friedman N. Chromatin conformation governs T-cell receptor Jbeta gene segment usage. Proc Natl Acad Sci U S A. 2012;109(39):15865–70.View ArticlePubMedPubMed CentralGoogle Scholar
- Pauza CD, Cairo C. Evolution and function of TCR Vgamma9 chain repertoire: It’s good to be public. Cell Immunol. 2015;296:22–30.View ArticlePubMedPubMed CentralGoogle Scholar
- McGowen MR, Grossman LI, Wildman DE. Dolphin genome provides evidence for adaptive evolution of nervous system genes and a molecular rate slowdown. Proc Biol sci / Royal Soc. 2012;279(1743):3643–51.View ArticleGoogle Scholar
- Reiss D, Marino L. Mirror self-recognition in the bottlenose dolphin: a case of cognitive convergence. Proc Natl Acad Sci U S A. 2001;98(10):5937–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Jaakkola K, Fellner W, Erb L, Rodriguez M, Guarino E. Understanding of the concept of numerically “less” by bottlenose dolphins (Tursiops truncatus). J Comp Psychol. 2005;119(3):296–303.View ArticlePubMedGoogle Scholar
- Das S, Li J, Hirano M, Sutoh Y, Herrin BR, Cooper MD. Evolution of two prototypic T cell lineages. Cell Immunol. 2015;296(1):87–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Hirano M, Das S, Guo P, Cooper MD. The evolution of adaptive immunity in vertebrates. Adv Immunol. 2011;109:125–57.View ArticlePubMedGoogle Scholar
- Krangel MS, Carabana J, Abbarategui I, Schlimgen R, Hawwari A. Enforcing order within a complex locus: current perspectives on the control of V(D)J recombination at the murine T-cell receptor alpha/delta locus. Immunol Rev. 2004;200:224–32.View ArticlePubMedGoogle Scholar
- Smit AFA, Hubley R & Green P: RepeatMasker open-4.0. 2013–2015 http://www.repeatmasker.org.
- Schwartz S, Zhang Z, Frazer KA, Smit A, Riemer C, Bouck J, Gibbs R, Hardison R, Miller W. PipMaker--a web server for aligning two genomic DNA sequences. Genome Res. 2000;10(4):577–86.View ArticlePubMedPubMed CentralGoogle Scholar
- Antonacci R, Mineccia M, Lefranc MP, Ashmaoui HM, Lanave C, Piccinni B, Pesole G, Hassanane MS, Massari S, Ciccarese S. Expression and genomic analyses of Camelus dromedarius T cell receptor delta (TRD) genes reveal a variable domain repertoire enlargement due to CDR3 diversification and somatic mutation. Mol Immunol. 2011;48(12–13):1384–96.View ArticlePubMedGoogle Scholar
- Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28(10):2731–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Felsenstein J. Confidence Limits on Phylogenies: An Approach Using the Bootstrap. Evolution. 1985;39(4):783–91.View ArticleGoogle Scholar
- Nei M, Kumar S. Molecular Evolution and Phylogenetics. 2000.Google Scholar
- Brochet X, Lefranc MP, 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 issue):W503–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Giudicelli V, Brochet X, Lefranc MP. IMGT/V-QUEST: IMGT standardized analysis of the immunoglobulin (IG) and T cell receptor (TR) nucleotide sequences. Cold Spring Harbor protocols. 2011;2011(6):695–715.PubMedGoogle Scholar
- Yousfi Monod M, Giudicelli V, Chaume D, Lefranc MP. IMGT/JunctionAnalysis: the first tool for the analysis of the immunoglobulin and T cell receptor complex V-J and V-D-J JUNCTIONs. Bioinformatics. 2004;20 Suppl 1:i379–85.View ArticlePubMedGoogle Scholar
- Giudicelli V, Lefranc MP. IMGT/junctionanalysis: IMGT standardized analysis of the V-J and V-D-J junctions of the rearranged immunoglobulins (IG) and T cell receptors (TR). Cold Spring Harbor protocols. 2011;2011(6):716–25.PubMedGoogle Scholar
- Miller SA, Dykes DD, Polesky HF. A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988;16(3):1215.View ArticlePubMedPubMed CentralGoogle Scholar
- Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J R Stat Soc Ser B Methodol. 1995;57(1):289–300.Google Scholar
- 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–8.View ArticlePubMedGoogle Scholar
- Kabsch W, Sander C. Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers. 1983;22(12):2577–637.View ArticlePubMedGoogle Scholar
- Jones DT. Protein secondary structure prediction based on position-specific scoring matrices. J Mol Biol. 1999;292(2):195–202.View ArticlePubMedGoogle Scholar
- Eswar N, Webb B, Marti-Renom MA, Madhusudhan MS, Eramian D, Shen MY, Pieper U, Sali A. Comparative protein structure modeling using Modeller. Current protocols in bioinformatics / editoral board, Andreas D Baxevanis [et al.] 2006, Chapter 5:Unit 5 6.Google Scholar
- Laskowski RA, MacArthur MW, Moss DS, Thornton JM. PROCHECK: a program to check the stereochemical quality of protein structures. J Appl Crystallogr. 1993;26(2):283–91.View ArticleGoogle Scholar
- Sippl MJ. Recognition of errors in three-dimensional structures of proteins. Proteins. 1993;17(4):355–62.View ArticlePubMedGoogle Scholar
- Lefranc MP, 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 issue):D1006–12.View ArticlePubMedGoogle Scholar
- 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–61.View ArticlePubMedGoogle Scholar
- Lefranc MP. IMGT Collier de Perles for the variable (V), constant (C), and groove (G) domains of IG, TR, MH, IgSF, and MhSF. Cold Spring Harbor protocols. 2011;2011(6):643–51.PubMedGoogle Scholar
- Lefranc MP. WHO-IUIS Nomenclature Subcommittee for immunoglobulins and T cell receptors report. Immunogenetics. 2007;59(12):899–902.View ArticlePubMedGoogle Scholar
- Lefranc MP. Antibody nomenclature: from IMGT-ONTOLOGY to INN definition. mAbs. 2011;3(1):1–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Lefranc MP, Pommie C, Kaas Q, Duprat E, Bosc N, Guiraudou D, Jean C, Ruiz M, Da Piedade I, Rouard M, et al. IMGT unique numbering for immunoglobulin and T cell receptor constant domains and Ig superfamily C-like domains. Dev Comp Immunol. 2005;29(3):185–203.View ArticlePubMedGoogle Scholar