Abstract
Background:
Amphibians are valuable models for comparative immunology. In the caudate Ambystoma mexicanum, the architecture of immunoglobulin loci resembles that of the anuran Xenopus tropicalis, although some antibody gene features are absent. Evidence supports the presence of T lymphocytes in axolotl, the expression of T cell receptor alpha, beta, and delta chains, and a restricted diversity in the delta chain. Here, we describe the T cell receptor loci in the A. mexicanum genome and compare them with X. tropicalis and other tetrapods.
Methods:
T cell receptor loci were mapped and annotated in the A. mexicanum genome (UKY_AMEXF1_1) using reference sequences from axolotl, X. tropicalis, human, and mouse. Gene models were refined with RNA sequencing data from spleen, lung, and liver.
Results:
The T cell receptor alpha and delta locus in axolotl shows an overall conserved structure compared with other tetrapods. The alpha locus contained a higher number of variable genes than the beta and delta loci, with a predominance of functional genes (ratio 3.06). No gene encoding the pre-T cell receptor chain alpha was identified. The delta locus harbored two conventional variable genes, but no expression was detected in RNA sequencing data, suggesting pseudogenization. Neither delta chain diversity genes nor gamma chain elements were found in the genome or spleen transcriptome. The beta locus displayed structural similarity to that of other tetrapods and included five translocons with diversity, joining, and constant segments. One constant gene consisted of two exons encoding two constant domains. Functional variable genes predominated in the beta locus (ratio 3.6).
Conclusion:
Our study reveals conserved but distinctive features of axolotl T cell receptor loci, including restricted delta-chain diversity, absence of gamma chain and pre-T cell receptor alpha, and structural novelty in the beta locus. These findings provide new insights into the evolution of T cell receptors in amphibians and offer a genomic framework to explore the links between adaptive immunity and tissue regeneration in A. mexicanum.
1 Introduction
T-cell receptors (TR) recognize peptide antigens and other pathogen-derived molecules presented by antigen-presenting cells (APCs), a process essential for initiating adaptive cellular immunity. These receptors are expressed on the surface of T lymphocytes and enable specific antigen recognition through their variable extracellular domains. TR-mediated recognition is coupled to signaling via the CD3 complex (conformed by CD3γ, CD3δ, CD3ϵ chains), which transmits activation signals that drive T-cell activation, effector and memory differentiation, and clonal expansion (, ).
The general organization of TR genes has remained remarkably conserved throughout 400 million years of gnathostome evolution. Unlike the high variability of immunoglobulin loci, TR loci exhibit structural stability across vertebrates (). A conventional TR is a disulfide-linked heterodimer composed of α and β chains, or γ and δ chains. Each of these four types of TR chains comprises two immunoglobulin superfamily domains: a membrane-proximal constant domain (C) and the antigen-binding variable domain (V). The V domains of TRβ and TRδ are assembled via somatic recombination of variable, diversity (D), and joining (J) genes, whereas the V domains of TRα and TRγ are assembled only by V and J genes (). This recombination is mediated by recombination signal sequences (RSS) flanking each gene, consisting of a conserved heptamer and nonamer motif separated by either a 12- or 23-base pair spacer. According to the 12/23 rule, recombination typically occurs between one RSS with a 12-bp spacer and another with a 23-bp spacer, ensuring proper assembly of V(D)J junctions ().
The TR α, β, γ, and δ chains are found across all jawed vertebrates, exhibiting significant conservation in both sequence and genomic arrangement. A distinctive feature is that the T cell receptor alpha (TRA) locus is embedded within the T cell receptor delta (TRD) locus (TRA-TRD locus), an organization conserved in all jawed vertebrates studied, including fish, amphibians, reptiles, birds, and mammals (–). However, the availability of non-model vertebrate genome sequences provides valuable insights into the distant origins of rearranging gene systems and their links to both adaptive and innate recognition processes (). This approach has led to the identification of additional TR chains, such as the New Antigen Receptor (NAR-TCR) in sharks (–) and the T cell receptor μ (TRμ) chain in marsupials and monotremes, which originated from TRδ gene duplication during early mammalian evolution (). Furthermore, in Xenopus tropicalis, the TRD locus contains canonical variable δ (Vδ) genes and VH-like genes termed VHδ, which are VH domains related to the variable domain of the immunoglobulin heavy chain, adapted as V-domains for TRδ chains (). VHδ is also found in fish, birds, and monotremes. Overall, this highlights the remarkable evolutionary plasticity of TR evolution, likely due to selective pressure imposed by pathogen recognition (, , ).
Amphibians are well-suited models for comparative immune system analysis (, ), due to their key evolutionary relations as the first tetrapods, bridging aquatic vertebrates (e.g., fishes) and terrestrial vertebrates (, ). Their immune system comprises all major components of adaptive immunity, including T and B lymphocytes, immunoglobulins, and Major histocompatibility complex (MHC) molecules, enabling direct comparisons across both ancestral and derived vertebrate lineages. In addition, exhibit unique immunological features, such as unconventional TR gene arrangements or limited receptor diversity (), offering insights into the evolutionary plasticity of the immune system. The study of TR loci in amphibians like Ambystoma mexicanum is particularly relevant because the immune system is increasingly recognized as a critical player in tissue repair and regeneration. Characterizing the genomic architecture and diversity of these loci not only informs our understanding of adaptive immunity in urodele amphibians but also provides a framework for investigating how immune components modulate regenerative processes. Such knowledge could facilitate the development of species-specific immunological tools, enhancing both biomedical research and conservation strategies.
A. mexicanum, a neotenic urodele amphibian endemic to the Mexico City valley, is an endangered species () and has one of the largest genome (32 Gb) among vertebrates sequenced to date (). Previously, we characterized the immunoglobulin heavy (IGH) and lambda (IGL) loci in the Ambystoma mexicanum, finding that it shares the same general syntenic architecture with X. tropicalis, but lacks the kappa locus (IGK) and other antibody features described in X. tropicalis (). Pre-genomic studies in A. mexicanum revealed the presence of T lymphocytes found in the spleen and thymus, as well as the presence of T cells expressing α, β, and δ chains (–). Of note, the junctional diversity of the TRδ chain is minimal (, ) and so far, no description of TRγ chains has been provided. Two genome assemblies are currently available for A. mexicanum: The AmbMex60DD genome assembly, based on a highly inbred laboratory strain (d/d), a two-year-old leucistic male (, ) has 27,157 unmapped scaffolds and revealed several positional and orientation inconsistencies in the IGH locus, likely reflecting assembly errors (). Recently, a new assembly, UKY_AMEXF1_1, generated from an F1 hybrid between A. mexicanum and A. tigrinum, both of wild origin has been publicly released. This assembly, which is currently the reference genome in GenBank, presents an improved chromosomal organization with 21 chromosomes and only 220 unmapped scaffolds (BioProject: PRJNA1165261), suggesting a more accurate annotation of complex loci.
To further investigate the germ-line structure of T cell receptors in A. mexicanum, we present here a genomic characterization and annotation of TR loci in the axolotl compared with X. tropicalis. This is one of the few amphibian species whose adaptive immune system has been extensively characterized at the genomic level (, ). In our previous analysis of the IGH and IGL loci in A. mexicanum, we reported the absence of certain components, such as the IGK locus and the pseudogenization of the IgF isotype (). This feature had also been described in X. tropicalis. Building on these findings, one of the main objectives of the present study is to determine whether, as observed in the case of immunoglobulins, TR loci in A. mexicanum also exhibit missing or divergent components compared with other tetrapods.
2 Results
2.1 Ambystoma mexicanum TRA-TRD locus
A phylogenetically conserved feature of the TRA and TRD locus organization in vertebrates is that both loci are closely embedded (, , , ) near the centromere of chr13. Accordingly, in A. mexicanum, the TRA and TRD locus are closely linked, with some TRD genes nested within the TRA locus. The TRA-TRD locus in A. mexicanum is located in chromosome 13p: 264.6-285.3 Mbp (size 20.7 Mbp) (Figures 1A, B) (Supplementary File 1; Supplementary Figure S1; Supplementary File 2; Table 1 and Glossary) and is not interrupted by gaps. In X. tropicalis, the TRA-TRD locus was mapped to chromosome 1 (0.72-1.24 Mbp; size 0.52 Mbp) (). We identified 61 T cell receptor alpha variable genes (TRAV), 46 of which are functional, flanked by canonic Recombination Signal Sequences of Variable genes (V-RSS) with 23-bp spacers, corresponding to Functional/Pseudogenes (F/P) ratio of 3.06. We identified 36 T cell receptor alpha join genes (TRAJ), of which 33 are functional. All TRAJ genes encode the canonical FGXG motif and have a 12-bp spacer and conserved heptamer and nonamer in their Recombination Signal Sequences of Join genes (J-RSS) (Supplementary File 1; Supplementary Figure S2). We found three TRAJ pseudogenes, one of which (TRAJ_036) is intercalated within the intron of the T cell receptor alpha constant gene (TRAC). Its J-RSS lacks the conserved 5´-CAC-3´ motif in the RSS heptamer, which is required for Recombination activator gene (RAG) recognition during V(D)J recombination (). The pseudogenes TRAJ_015 and TRAJ_023 contain a frameshift (Supplementary File 2; Table 1, Glossary and GFF file). A single TRAC gene with canonical exons, including cytoplasmic (M2), transmembrane (M1), and C-Ig domain exons (Figures 1D-F; Table 1), has three glycosylation sites in 42-45 (NDTE), 72-75 (NDTQ), and 106-109 (NESF). In the transmembrane region (TM), residues Cys225, Glu237, Arg251, Lys256, Asn261, Tyr265, and Trp269, which interact with CD3, are mostly conserved in A. mexicanum except for Arg251, which is replaced with Lys (, ). The connecting peptide motif (FETDXXLN), another important site in the TM region for the transduction of activation upon antigen recognition (), is well conserved in the axolotl. Furthermore, Cα sequence alignment across human, mice, opossums, frogs, and axolotls reveal limited conservation of the AB loop in X. tropicalis and A. mexicanum regarding mammals (Supplementary File 1; Supplementary Figure S3).
Figure 1
Table 1
| Gene type | Number of functional genes | Percentage of functional genes (%) | Number of pseudogenes | Percentage of pseudogenes (%) | Total |
|---|---|---|---|---|---|
| TRAV | 46 | 75.4 | 15 | 24.59 | 61 |
| TRAJ | 33 | 91.6 | 3 | 8.3 | 36 |
| TRAC | 1 | 100 | 0 | 0 | 1 |
| TRDV | 1 | 50 | 1 | 50 | 2 |
| TRDD | There is no evidence of the presence of TRDD genes | ||||
| TRDJ | 2 | 100 | 0 | 0 | 2 |
| TRDC | 1 | 100 | 0 | 0 | 1 |
Summary of total number and percentage of functional and pseudogenes of V, D, J and C genes found in the TRA-TRD locus mapped in Chr 13.
The TRD locus harbors two conventional T cell receptor delta variable genes (TRDV) flanked by a canonical V-RSS with a 23-bp spacer; however, one of them is not transcribed based on spleen, lung, and liver RNA-seq data, suggesting it is a pseudogene. Only two T cell receptor delta join genes (TRDJ) were found (Figures 1C-E) (Supplementary File 1; Supplementary Figure S1). Both TRDJ genes appear to be functional; however, the consensus J motif (FGXG) encoding the di-glycine bulge is not present in Jδ1 (FKKG), whereas Jδ2 retains the canonical sequence. The J-RSS is conserved in both genes, with a 12-bp spacer (Supplementary File 1; Supplementary Figure S4). No Dδ genes or their corresponding recombination signal sequences of diversity (D-RSS) were identified. The exon organization of the single T cell receptor delta constant gene (TRDC) was found with canonical exons, including cytoplasmic (M2), transmembrane (M1), and C-Ig domain exons (Figure 1C; Table 1). This TRDC exon encodes an N-glycosylation site (NSSS, pos 36-39).
In all studied vertebrates so far, the TRA and the TRD locus are genetically linked (, –) in different vertebrates such as frog, human, mouse, and opossum. In all of these species, the TRA-TRD locus is flanked by the METTL3, SALL2, DAD1, and ABHD4 genes (, , , , –). However, the A. mexicanum locus is flanked by NUMP and LPCAT4 in the centromere direction and the KLHL33 gene in the telomeric direction (Figure 2; Supplementary File 3: Table 2).
Figure 2
Table 2
| Gene type | Number of functional genes | Percentage of functional genes (%) | Number of pseudogenes | Percentage of pseudogenes (%) | Total |
|---|---|---|---|---|---|
| TRBV | 18 | 78.26 | 5 | 21.7 | 23 |
| TRBD | 5 | 100 | 0 | 0 | 5 |
| TRBJ | 17 | 100 | 0 | 17 | |
| TRBC | 6 | 100 | 0 | 0 | 6 |
Summary of total number and percentage of functional and pseudogenesof V, D, J and C genes found in the TRB locus mapped in the Chr 3p.
2.2 Ambystoma mexicanum TRB locus
The T cell receptor beta (TRB) locus in X. tropicalis has not been previously described. We mapped the X. tropicalis TRB locus to chr7p (3–8 Mbp), flanked by the DPH-like and EPHRIN genes towards the centromere, and NOBOX and CNCL in the telomeric direction (Figure 3). A trypsin gene cluster (PRSS) was found between T cell receptor beta join genes (TRBJ) and T cell receptor beta variables clusters (TRBV) (Figures 4A-C) (Supplementary File 1; Supplementary Figure S5; Supplementary File 3: Table 2). A similar organization is observed in the human and mouse TRB locus; however, in X. tropicalis, the locus is inverted regarding the EPHRIN gene.
Figure 3
Figure 4
In A. mexicanum, the TRB locus was identified on chromosome 3p (30.03-34.85 Mbp, size 4.82) (Figures 4A, B; Supplementary File 2; Table 1 and GFF file) and contains a single gap in position 32.62 Mbp. The trypsin (PRSS) gene cluster divides the TRB locus into a canonical locus towards the centromere (referred to hereafter as cluster A), and a TRBV gene cluster (cluster B) towards the telomere (Figure 4B). Cluster A compromises five tandem TRBC-TRBJ-TRBD translocons and a TRBV gene cluster comprising 13 TRBV genes (Figure 4C). Cluster B contains 10 TRBV genes. Overall, there are 18 functional TRBV genes (9 in cluster A and 9 in cluster B) to F/P ratio of 3.6. Additionally, there are 5 functional T cell receptor beta diversity genes (TRBD). All TRBD genes exhibit 12-pb and 23-pb spacers. Recombination signal sequences of diversity genes (D-RSS) at their respective flanks (Supplementary File 1; Supplementary Figure S6). Furthermore, we found 17 functional TRBJ genes encoding the FGXG motif, with canonical 12-bp spacer J-RSS (Supplementary File 1; Supplementary Figure S7).
We identified five functional T cell receptor beta constant genes (TRBC), one associated with each translocon (Figures 4D, E). Spleen, lung, and liver RNA-seq data revealed that all functional genes are actively transcribed (Figures 4F, G). The functional genes feature the conserved Trp at position 41, the Leu at position 86, and the two characteristic Cys residues of the constant genes at positions 23 and 104 (Based on IMGT numbering of human TRBC1). Genes TRBC_001, 002, 004 and 005 share the typical TRB gene structure composed by a single Cβ exon with a glycosylation site (TRBC_001 100-103 (NITV), TRBC_002 100-103 (NITV), TRBC_004 7-10 (NVTQ), TRBC_005 51-54 (NRTK)), a M1 and M2 exons encoding a linker peptide and the transmembrane region, respectively (–). The TRBC_003 is unusual because it comprises two Cβ-domain exons, (TRBC_003_1 and TRBC_003_2) (Figures 4E-G; Table 2). The TRBC_003_1 exon encodes for a glycosylation site (4–7 NITQ), whereas the TRBC_003_2 exon lacks predicted N-glycosylation sites.
In mammal Cβ-domain, the FG loop is implicated in the interaction with the CD3 complex (). Sequence alignment of the Cβ domains from human, mouse, opossum, frog, and axolotl revealed the absence of the FG loop in all X. tropicalis and A. mexicanum Cβ. In all compared amphibian Cβ domains, the proline residue at position 232 is conserved, except in TRBC_003_2, where it is replaced by an alanine (Ala 228) (Supplementary File 1; Supplementary Figure S8) (, ).
2.3 No evidence of an Ambystoma mexicanum TRG locus
We use BLAST and HMMER alignment-based search tools, either in the genome or in spleen RNA-seq transcriptome data, we found no evidence of the existence of the T cell receptor gamma locus (TRG) in A. mexicanum. In X. tropicalis, the TRG locus is located on chromosome 6 (Chr6: 63.4-63.6 Mb), flanked proximally by the STARD3NL, EPDR1, SFP4, GPR141, ELMO1, AOAH, ANLN, and MATCAP2 genes and the AMPH, POU6F2, NPSR1, BMPER, BBS9, NT5C3A, RP9, and VELO1 genes at the distal flank. In humans, the TRG locus is located on chromosome 7 (Chr7: 38.24–38.36 Mb), with a genomic neighborhood like that of X. tropicalis, although the orientation of the flanking genes is inverted. In this case, the proximal flanking genes are AMPH, POU6F2, NPSR1, BMPER, BBS9, NT5C3A, RP9, and VELO1, while the distal flanking genes are STARD3NL, EPDR1, SFP4, GPR141, ELMO1, AOAH, ANLN, and MATCAP2. In contrast, A. mexicanum exhibits a genomic architecture similar to that observed in humans, although located on chromosome 5p (Chr5p: 834.5–863.3 Mb), with no evidence of the TRG locus in this region (Figure 5; Supplementary File 3; Table 2).
Figure 5
2.4 V-intron length
V genes are composed of two exons: Exon 1 encodes the L1 part of the leader peptide, whereas exon 2 encodes the L2-part of the leader peptide and the V-region. Both exons are separated by the V-intron that typically ranges in size from 80 to 250 bp (). In axolotl, the average V-intron length was 758 bp for TRA, 142 bp for TRD (considering two genes), and 237 bp for TRB (Supplementary File 2; Table 1, columns L-O). We analyzed the distribution of V-intron lengths in functional TRAV (P=0.0002) and TRBV (P=0.51) genes. No significant differences were observed between functional and non-functional TRBV genes. In contrast, a significant difference was detected in TRAV, suggesting that most TRAV genes with long V-introns are functional (Supplementary File 1; Supplementary Figure S9).
To investigate the relationship between V-intron length and V gene functionality in the TRA and TRB loci of A. mexicanum, we assessed whether the presence of long introns (>650 bp) was more frequent in non-functional TRAV and TRBV genes, as previously reported for the IGH and IGL loci (). Fisher’s exact test was applied to the TRAV and TRBV genes. For TRAV, the analysis revealed that the odds of finding a long intron in a non-functional gene were 0, resulting in an odds ratio of 0.0 (p=0.01035; 95% CI: 0.0–0.66). This suggests a significant depletion of long introns among non-functional TRAV genes. In the case of TRBV, no significant association was found (p=1; 95% CI: 0.0–233.15), likely due to the limited number of observations, making statistical comparison uninformative.
2.5 Ambystoma mexicanum PTCRA locus
In humans and mice, the pre-TCRα participates in αβ T cell early development in association with the TRβ chain at the surface of thymocytes. It is known that this receptor is absent in non-mammalian species such as Xenopus spp. and Danio rerio (Zebrafish) (). We performed BLAST and HMMER searches in the A. mexicanum genome and transcriptome and found no evidence of a PTCRA ortholog. Additionally, synteny analysis of the PTCRA locus across multiple species revealed conserved synteny between X. laevis, X. tropicalis, and A. mexicanum, confirming the absence of PTCRA in A. mexicanum in contrast to mammals and sauropsids (Figure 6).
Figure 6
3 Discussion
The UKY_AMEXF1_1 genome assembly of A. mexicanum enabled a comprehensive characterization of its TR loci, supported by an increased N50 = 1.5 compared with the AmbMex60DD genome version N50 = 1.2. The overall structure of the TRA-TRD locus is conserved, but the TRD locus exhibits strikingly low combinatorial diversity, and the TRG locus is absent, suggesting that bona fide TRγδ T cells may be lacking in axolotls. The TRB locus displays a conserved structure with tandem duplications of the TRBD-TRBJ-TRBC translocon, including a particular TRBC gene with two Cβ domains. As in X. tropicalis, the PTCRA gene is absent.
These findings provide important insights into the genomic organization and evolutionary constraints of TR loci in urodele amphibians. The limited diversity of TRD genes and absence of the TRG locus highlight unique features of the axolotl adaptive immune system, which may have implications for understanding T-cell function in regeneration and immune response. Overall, this work establishes a foundation for comparative immunogenomic studies across amphibians and other vertebrates.
Here we confirm such findings. We found similarities in the overall structure of the TRA-TRD locus, as well as a strikingly low combinatorial diversity of the TRD locus and the absence of the TRG locus, implying the absence of bona fide TRγδ T cells in axolotl. As for the TRB locus, we describe a conserved structure, with tandem duplications of the TRBD-TRBJ-TRBC translocon and a particular TRBC gene composed of two Cβ domains. As in X. tropicalis, the axolotl genome also lacks the PTCRA gene.
A. mexicanum’s exceptionally large genome (32 Gb) posed significant challenges for its assembly. The AmbMex60DD (white strain, d/d) version was released in 2021, based on 30× genomic coverage of 28 chromosomes and 27,157 unmapped scaffolds. This assembly presented several positional and orientation inconsistencies in the IGH locus, likely reflecting assembly errors (). The current genome assembly (UKY_AMEXF1_1) was released in 2024 and has an increased genomic coverage (48×), 21 chromosomes, and only 220 unmapped scaffolds, suggesting that a more accurate complex loci annotation can be achieved. However, as with AmbMex60DD, no publicly available data on local coverage is currently provided for UKY_AMEXF1_1. Therefore, it is not possible to assess coverage-based metrics for individual TR genes.
In agreement with our previous analyses using the AmbMex60DD assembly, the overall expression patterns of TR genes did not substantially change in the present study (Supplementary File 4; Supplementary Table 3). Remarkably, we identified two genes, TRAV_060 and TRDV_001, that exhibit a structurally complete configuration, including the SP, canonical donor and acceptor splice sites, conserved methionine’s at positions 23 and 104, and canonical RSS. However, neither of these genes showed detectable expression in the analyzed tissues. Conversely, two additional genes, TRAV_025 and TRBV_012, also retained an intact genomic organization and displayed transcriptional evidence of the V gene, yet lacked detectable expression of the corresponding SP. According to our classification criteria, these cases were therefore categorized as pseudogenes, since evidence of both V gene and SP expression was required to consider a gene as functional.
Despite the mentioned limitations, the remarkably stable of the TRA-TRD locus across species, maintaining a consistent genomic architecture for over 340 million years of evolutionary history, highlights strong evolutionary constraints on its organization (). The fact that TRA and TRD remain genetically linked in all examined species (, –) reinforces the functional importance of their physical association. This overall structure, conserved in X. tropicalis, alligators, birds, and mammals. In all of these species, the TRA-TRD locus is flanked by the METTL3, SALL2, DAD1, and ABHD4 genes (, , , , , , ) suggesting that this syntenic arrangement may be critical for maintaining locus integrity. In contrast, the distinct configuration in A. mexicanum; located on chromosome 13p, near the centromere and flanked by KLHL33, likely reflects a lineage-specific chromosomal rearrangement that preserved the internal gene order, indicating that positional changes do not necessarily disrupt locus function (Figure 2) (Supplementary File 1; Supplementary Figure S1).
In A. mexicanum, the residue equivalent to mammalian TRα Arg251 is replaced by Lys, a substitution that is unlikely to affect its functional contact with CD3δ, given the chemical similarity of both residues (, ). The connecting peptide motif in the Cα transmembrane region is conserved, maintaining its role in signal transduction from the αβ heterodimer to the CD3/ζ complex ().
Interestingly, as in X. tropicalis, A. mexicanum lacks the AB-loop, a structural feature essential for CD3ϵδ contact and T cell activation in mammals (, ). The absence of the AB-loop in amphibians may reflect a distinct co-evolutionary trajectory of TR and CD3 complexes compared to mammals ().
We confirm the restricted diversity of germline genes previously reported for the TRD locus. As described by André, et al. () the cause of such restricted diversity is determined by a single functional TRDV and two TRDJ genes, but significantly, we confirm the absence of TRDD genes in the germline. As in other tetrapod TRD loci, the V and J genes are flanked by 23-bp and 12-bp spaced RSS, respectively. Hence, direct TRDV-TRDJ junctions do not violate the 12/23 rule (). TRDV-TRDJ junctions have been described in a subset of human acute lymphoblastic leukemias (). To our knowledge, A. mexicanum is the first vertebrate capable of non-pathological direct TRDV-TRDJ recombination.
A notable feature of the TRA-TRD locus in the axolotl is that the only TRDV gene is of the conventional Vδ type and not of the VHδ type, in striking contrast to X. tropicalis, which contains 14 VHδ and 2 conventional Vδ genes (). The VHδ type is widely distributed among non-placental mammals and other vertebrates (, , , , , 53, 54), and it remains to be determined if this type of element was lost in all or some caudates or if it was never present.
Compared with the TRα and β chains, the presence of TRγ and TRδ chains shows more heterogeneity across different taxa (). In scaled reptiles, there is an absence of the TRD and TRG loci (55). Although the syntenic blocks flanking the TRG locus in Xenopus were identified in A. mexicanum, no genomic and transcriptomic evidence of the presence of the TRG locus was found. These results indicate that axolotl lacks true γδ T cells. It remains to be determined the functional role of the single TRδ chain, and if it pairs at all with itself or with another chain, but due to its invariant structural nature, it may function essentially as a Pattern Recognition Receptor (PRR), similarly to BTNL/Btnl family of innate γδ TR (56).
In this study, we updated and expanded the genomic annotation of the TRB locus located on chr3p previously described using cDNA libraries (57–59). The TRB locus in tetrapods is generally organized into TRBD-TRBJ-TRBC units, resembling the organization of the lambda chain locus (60). For instance, sheep possess three TRBD-TRBJ-TRBC tandem units, while rabbits, mice, and humans have two (, 61, 62). In comparison, X. tropicalis has a simpler configuration with a single unit. Remarkably, A. mexicanum displays a more complex arrangement, consisting of five tandem TRBD-TRBJ-TRBC units.
In conventional TRβ chains, the Cβ domain interacts with CD3 through the FG-loop, contributing to signal transduction upon MHC-peptide recognition. In mammals, mutations in the FG-loop alter the CD8+ and CD4+ T cell proportions, and it is associated with a poor antigen response (63). Our study in axolotl is in agreement with a previous report by Kim, et al., reporting the absence of the FG-loop in non-mammal vertebrates (Supplementary File 1; Supplementary Figure S8) (, ). We identified a novel TRBC gene (TRBC_003) that features two Cβ domains, along with transmembrane and cytoplasmic domains, each encoded by a separate exon. The TRBC_003 gene is actively transcribed in the spleen and associated with putatively functional Dβ and Jβ genes. Moreover, in the distal membrane exon (TRBC_003_1), there is a conserved Pro232, which is a relevant position of the FG-loop in mammals with a single N-glycosylation site. In contrast, in the proximal membrane exon (TRBC_003_2), there is an Ala in the 228 position and no N-glycosylation sites. This structural evidence suggests that the proximal exon may have arisen from a duplication of the distal exon. Whether a TRβ product of the TRBC_003 gene pairs with a TRα chain, how it interacts with CD3, and how MHC-peptide interaction takes place remain as open questions.
Genome sequencing of F1 crosses between A. mexicanum and A. tigrinum have revealed genomic regions of high polymorphism (64, 65), however the TR loci are outside these regions. We consider that spurious contributions of allelic variation to gene count are likely minimal. Moreover, genome information derives from a single individual, and our findings may thus fail to capture the extent of intra-species variability, particularly copy number variation in Adaptive Immune Receptor Repertoire (AIRR) loci, which are well documented in mice, macaques, and humans (66–68), that may explain discrepancies between the UKY_AMEXF1_1 and the AmbMex60DD assemblies.
We observed that many TRAV genes in A. mexicanum possess notably long V-introns, some exceeding 1,000 bp; and are considerably larger than IGH and IGL V-introns in A. mexicanum (). In the Gallus gallus genome, the average TRB V-intron length is 400 bp (69). The biological implication of longer TR V-introns is uncertain. Long introns have been reported to impose evolutionary costs by increasing the energetic demand due to the greater nucleotide investment and extended transcription time required. Additionally, they may compromise the fidelity of mature mRNA and expand the sequence space available for allelic variation and aberrant splicing (70). Intron length has also been suggested to play a functional role in evolutionary dynamics (71). Comeron and collaborators (72) proposed that extensively long introns may enhance the efficiency of natural selection by alleviating Hill-Robertson (HR) interference, a phenomenon where selection acting on linked loci reduces selective efficacy. In this context, long V-introns may act as spacer regions that decouple selective pressures acting on neighboring functional elements.
Moreover, HR interference might prevent their elimination via purifying selection if these introns are linked to active V genes. Consequently, such introns could facilitate the emergence of alternative splicing events or the production of non-functional transcripts, potentially affecting the expression and functionality of the antigen receptor repertoire (70, 73, 74). Collectively, these findings reinforce the idea that introns do not merely represent structural and energetic burdens but may also play key roles in gene generation, conservation, and diversification; particularly in immune-related loci such as those encoding T cell receptors.
It is noteworthy that the functional-to-pseudogene ratio in the A. mexicanum TRB is 3.6 and 3.06 in the TRA locus, compared with 0.9 and 2.7 in the IGH locus (AmbMex60DD) () and UKY_AMEXF1_11 assemblies, respectively. Within the framework of the “birth-and-death” model of T and B cell receptor gene evolution (75), higher ratios may indicate a slower accumulation of pseudogenes, potentially reflecting more recent functional gene birth events, stronger purifying selection, or lower rates of gene inactivation in T cell receptors compared with B cell receptors. While such differences could be stochastic, they may also reflect distinct selective pressures acting on these repertoires. The comparison between genome versions further shows that improvements in assembly quality can refine gene counts and alter calculated ratios, underscoring the importance of high-quality chromosome-level genome assemblies and thorough manual curation of AIRR loci for robust evolutionary inferences (76).
Its absence in amphibians implies that early αβ T cell development proceeds via alternative mechanisms, possibly involving different surrogate chains or signaling pathways. This reinforces the hypothesis that PTCRA originated as an amniote-specific innovation, rather than an ancestral gene lost independently in teleosts and amphibians. Identifying how amphibians compensate for the absence of PTCRA could provide insights into the evolution and diversification of T cell developmental programs in vertebrates.
The review and characterization of immune components in A. mexicanum are essential for understanding the cellular processes in which cells interact dynamically and persistently. These components have been shown to play a key role in modulating such interactions in other vertebrates with comparable regenerative capacity, influencing both the persistence of immune responses and the regulation of mechanisms underlying the regeneration of complex structures such as limbs, tail, heart, retina, and spinal cord. However, further studies are needed to clarify the specific role of T cells in this process in A. mexicanum (77–80).
In conclusion, the A. mexicanum TRB locus exhibits greater diversity than the TRA and TRD loci. Notably, no evidence of bona fide γδ T cells was found. This study leaves open questions regarding the composition of T-cell subpopulations and the pairing of TR chains, particularly the δ chain in the absence of the γ chain in axolotl. The presence of two constant domains in TRB_003 warrants further investigation to clarify their role in antigen recognition, functionality, and pairing with the α chain. Such insights are crucial to understanding the impact on this endangered species, which already presents a marked deficit of heterozygosity, reflecting substantial inbreeding and increasing vulnerability to infectious diseases ().
This work provides valuable insights for comparative evolutionary analyses in tetrapods and advances our knowledge of immune response in caudate amphibians. Moreover, it may aid in solving specific questions regarding the role of acquired immunity in the regulation of the immune response implicated in tissue regeneration (81).
4 Methods
4.1 Ambystoma mexicanum genome and transcriptome data
The published sequence of the A. mexicanum haploid genome (UKY_AmexF1_1; GenBank assembly: GCF_040938575.1) was generated from an adult female F1 hybrid (isolate Amex_F1_6; BioSample SAMN41071122) derived from a cross between a female A. mexicanum (isolate Mex_15411; BioSample SAMN43142723) and a male A. tigrinum (isolate Tig_M23; BioSample SAMN43142724). The phased haploid assembly has 48x coverage was sequenced using PacBio and Illumina HiSeq and Hi-C data (82). The UKY_AmexF1_1 assembly consists of 21 chromosomes and 220 unplaced scaffolds, with a scaffold N50 of 1.5 Gb.
The corresponding genome annotation file GCF_040938575.1_UKY_AmexF1_1_genomic.gff.gz; and gap positions GCF_040938575.1_UKY_AmexF1_1_genomic_gaps.txt.gz, were retrieved from https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/040/938/575/.
To validate gene models, we used previously published RNA-seq and transcriptome data available in the NCBI SRA database. Specifically, RNA-seq coverage bigWig (BW) files from spleen, liver, and lung (SRR15610271, SRR15610267, SRR15610267), obtained from the NCBI FTP repository (https://ftp.ncbi.nlm.nih.gov/genomes/all/GCF/040/938/575/GCF_040938575.1_UKY_AmexF1_1/RNASeq_coverage_graphs/), were used as visual support for manual curation without further processing. Additionally, transcriptome datasets (BioProject PRJNA378970) from spleen, liver, and lung (SRR5341570, SRR5341572, SRR5341571) () were aligned to the genome using STAR (83) with default parameters and a maximum intron size of 3000 bp. Gene-level quantification was performed using the –quantMode GeneCounts option.
Mapping statistics for each dataset and assembly including input reads, uniquely mapped reads, spliced alignments, mismatch rate per base, and reads discarded for being too short. RNA-seq datasets were obtained from adult tissues (NCBI SRA: SRR5341570: spleen; SRR5341572: liver; SRR15610267: lung) and aligned to both AMEX_F1_1 and AmbMex60DD assemblies. These statistics provide a benchmark for expression analysis and support the reliability of TR gene annotation across assemblies (Supplementary File 4; Supplementary Table 3).
4.2 TR loci mapping
Reference sequences (cDNA) for TRA, TRD, TRB, and TRG loci from X. tropicalis, X. laevis, and A. mexicanum were obtained from NCBI (https://www.ncbi.nlm.nih.gov/). These sequences were used to map the TRA, TRD, and TRB loci using TBLASTX and Exonerate (EST2genome alignment model). Hits with significant alignment scores (e-value < 1.0E-05 for BLASTX, score > 100 for Exonerate) were exported as GFF3 files. These annotations were visualized and manually curated using the Integrative Genomics Viewer (IGV) (84). To complement homology-based annotation and minimize the risk of missing novel or lineage-specific V(D)J genes, we developed a custom pipeline to detect RSSs according to the canonical 12/23 rule. The workflow comprised: (i) BLAST alignments with bitscore filtering and conversion to GFF to identify scaffolds or chromosomes of interest; (ii) Exonerate-based exon and gene detection; (iii) RSS search using HMMER with tbl-to-GFF conversion; (iv) redundancy reduction across gene and exon annotations; (v) overlap analyses to confirm V genes and their signal peptides; (vi) refinement of V gene and RSS-J coordinates with Miniprot (protein-to-genome aligner); and (vii) identification of candidate D genes based on flanking 5′ and 3′ RSSs. This pipeline was applied to the TRA, TRD, and TRB loci in the AmbMex60DD assembly to uncover additional putative V, D, and J genes. The search database was built from TR gene models previously described in X. tropicalis and X. laevis. Genes identified through RSS were integrated with homology-based results and manually curated. To refine annotation, the AmbMex60DD-derived sequences were subsequently aligned to the latest reference genome (AMEX_F1_1), yielding the final TR loci annotation. All TR genes were named by a provisional numeric identifier. All our annotations are compliant with the IUIS T-cell Receptor and Immunoglobulin Nomenclature Sub-Committee, except for the fact that individual gene coverage and loci coverage are not publicly available (85).
4.3 Definition of V, D, and J functionality
Functionality assessments were performed based on IMGT (86) criteria. To classify a V, D, or J gene as functional (F), each coding region was required to have an open reading frame, proper splicing sites, and recombination signals with 12/23 spacers. For V-exons, the presence of conserved residues Cys23, Trp41, Trp52, and Cys104 was mandatory. Genes were classified as pseudogenes (P) if they contained stop codons, lacked leader peptide exon and/or RSS and were frame-shifted. For J-exons, the di-glycine bulge (FGXG) was required.
4.4 Search for TRDD and TRBD genes
We constructed a Hidden Markov Model (HMM) profile with the HMMER3 (-hmmbuild option) (87) to represent the sequence structure of the genes of TRBD and their associated RSS’s. This profile was generated with multiple sequence alignments based on TRBD genes from A. mexicanum previosuly published by Fellah (59). The genes were flanked at the 5’ by a 12 bp-spaced RSS, and at the 3’ end by a 23 bp-spaced RSS. The same HMM profile was subsequently used to search for the TRA-TRD locus.
4.5 PTCRA gene search
A multiple sequence alignment of PTCRA orthologs described by Smelthy et al. () was used to build an HMM profile with the HMMER3 (-hmmbuild option). This probabilistic model captures evolutionary changes in conserved amino acids across related sequences (87). The resulting HMM profile was applied to the A. mexicanum proteome using -hmmsearch. Additional searches were performed with the PFAM model PF15028 for PTCRA.
For synteny analysis, we first identified PTCRA-flanking genes in human, opossum, birds, reptiles, and frogs. Orthologous regions were then located in A. mexicanum, followed by manual curation of the surrounding genes to ensure the annotation was correct.
4.6 Intron length analysis
V-intron length was calculated from the exon coordinates of the respective locus annotation file (Supplementary 2; Table 1, P column). Due to the presence of abnormally long intron, we used a non-parametric Wilcoxon rank-sum text to compare intron lengths between functional and pseudogene V genes. Enrichment of long introns in functional genes was further evaluated using Fisher’s exact test in R.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Author contributions
DP-O: Methodology, Data curation, Investigation, Formal Analysis, Writing – original draft, Writing – review & editing. SSR-H: Methodology, Writing – review & editing, Data curation, Investigation, Writing – original draft, Formal Analysis. EG-L: Methodology, Data curation, Formal Analysis, Writing – review & editing. JT-S: Data curation, Methodology, Writing – review & editing. HV-T: Writing – review & editing. EC-Q: Writing – review & editing. CL-M: Writing – review & editing, Conceptualization. JM-B: Data curation, Methodology, Writing – review & editing, Supervision, Investigation, Writing – original draft, Formal Analysis, Conceptualization, Project administration.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. DL-PO CVU 890789 and SSR-H CVU-921827 receive grants from the Secretariat of Science, Humanities, Technology, and Innovation (Secihti).
Acknowledgments
Robert D. Miller, Jordan Sampson, Horacio Mena González, and Luis Zambrano for helpful discussions.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Correction note
This article has been corrected with minor changes. These changes do not impact the scientific content of the article.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2025.1656386/full#supplementary-material
Footnotes
1.^Saint Remy-Hernández S, Pacheco-Olvera DL, Godoy-Lozano EE, Miguel-Ruiz JA, Téllez-Sosa J, Valdovinos-Torres H, et al. An update of immunoglobulin loci in Ambystoma mexicanum. (2025).
References
1
LitmanGWRastJPFugmannSD. The origins of vertebrate adaptive immunity. Nat Rev Immunol. (2010) 10:543–53. doi: 10.1038/nri2807
2
AltFWZhangYMengFLGuoCSchwerB. Mechanisms of programmed DNA lesions and genomic instability in the immune system. Cell. (2013) 152:417–29. doi: 10.1016/j.cell.2013.01.007
3
HansenVLMillerRD. The Evolution and Structure of Atypical T Cell Receptors. In: HsuEDu PasquierL, editors. Pathogen-Host Interactions: Antigenic Variation v Somatic Adaptations. Springer International Publishing, Cham (2015). p. 265–78.
4
KuhnsMSDavisMMGarciaKC. Deconstructing the form and function of the TCR/CD3 complex. Immunity. (2006) 24:133–9. doi: 10.1016/j.immuni.2006.01.006
5
CalisJJARosenbergBR. Characterizing immune repertoires by high throughput sequencing: strategies and applications. Trends Immunol. (2014) 35:581–90. doi: 10.1016/j.it.2014.09.004
6
ParraZEOhtaYCriscitielloMFFlajnikMFMillerRD. The dynamic TCRδ: TCRδ chains in the amphibian Xenopus tropicalis utilize antibody-like V genes. Eur J Immunol. (2010) 40:2319–29. doi: 10.1002/eji.201040515
7
MassariSLinguitiGGiannicoFD’AddabboPCiccareseSAntonacciR. The genomic organisation of the TRA/TRD locus validates the peculiar characteristics of dromedary δ-chain expression. Genes (Basel). (2021) 12:544. doi: 10.3390/genes12040544
8
LefrancMPLefrancG. IMGT®Homo sapiens IG and TR loci, gene order, CNV and haplotypes: New concepts as a paradigm for jawed vertebrates genome assemblies. Biomolecules. (2022) 12:381. doi: 10.3390/biom12030381
9
WangXHuangJWangPWangRWangCYuDet al. Analysis of the Chinese alligator TCRα/δ Loci reveals the evolutionary pattern of atypical TCRδ/TCRμ in tetrapods. J Immunol. (2020) 205:637–47. doi: 10.4049/jimmunol.2000257
10
ParraZEMillerRD. Comparative analysis of the chicken TCRα/δ locus. Immunogenetics. (2012) 64:641–5. doi: 10.1007/s00251-012-0621-5
11
SeelyeSLChenPLDeissTCCriscitielloMF. Genomic organization of the zebrafish (Danio rerio) T cell receptor alpha/delta locus and analysis of expressed products. Immunogenetics. (2016) 68:365–79. doi: 10.1007/s00251-016-0904-3
12
LitmanGWAndersonMKRastJP. Evolution of antigen binding receptors. Annu Rev Immunol. (1999) 17:109–47. doi: 10.1146/annurev.immunol.17.1.109
13
CriscitielloMFSaltisMFlajnikMF. An evolutionarily mobile antigen receptor variable region gene: doubly rearranging NAR-TcR genes in sharks. Proc Natl Acad Sci U S A. (2006) 103:5036–41. doi: 10.1073/pnas.0507074103
14
ParraZEBakerMLHathawayJLopezAMTrujilloJSharpAet al. Comparative genomic analysis and evolution of the T cell receptor loci in the opossum Monodelphis domestica. BMC Genomics. (2008) 9:111. doi: 10.1186/1471-2164-9-111
15
ParraZEMitchellKDalloulRAMillerRD. A second TCRδ locus in Galliformes uses antibody-like V domains: insight into the evolution of TCRδ and TCRμ genes in tetrapods. J Immunol. (2012) 188:3912–9. doi: 10.4049/jimmunol.1103521
16
ParraZELillieMMillerRD. A model for the evolution of the mammalian t-cell receptor α/δ and μ loci based on evidence from the duckbill Platypus. Mol Biol Evol. (2012) 29:3205–14. doi: 10.1093/molbev/mss128
17
MarchalonisJJBernsteinRMShenSXSchluterSF. Emergence of the immunoglobulin family: conservation in protein sequence and plasticity in gene organization. Glycobiology. (1996) 6:657–63. doi: 10.1093/glycob/6.7.657
18
CharlemagneJFellahJSDe GuerraAKerfournFPartulaS. T-cell receptors in ectothermic vertebrates. Immunol Rev. (1998) 166:87–102. doi: 10.1111/j.1600-065X.1998.tb01255.x
19
LitmanGWCannonJPDishawLJ. Reconstructing immune phylogeny: new perspectives. Nat Rev Immunol. (2005) 5:866–79. doi: 10.1038/nri1712
20
RobertJOhtaY. Comparative and developmental study of the immune system in Xenopus. Dev Dyn. (2009) 238:1249–70. doi: 10.1002/dvdy.21891
21
Du PasquierLSchwagerJFlajnikMF. The immune system of Xenopus. Annu Rev Immunol. (1989) 7:251–75. doi: 10.1146/annurev.iy.07.040189.001343
22
VossSRWoodcockMRZambranoL. A tale of two axolotls. Bioscience. (2015) 65:1134–40. doi: 10.1093/biosci/biv153
23
NowoshilowSSchloissnigSFeiJFDahlAPangAWCPippelMet al. The axolotl genome and the evolution of key tissue formation regulators. Nature. (2018) 554:50–5. doi: 10.1038/nature25458
24
Martinez-BarnetcheJGodoy-LozanoEESaint Remy-HernándezSPacheco-OlveraDLTéllez-SosaJValdovinos-TorresHet al. Characterization of immunoglobulin loci in the gigantic genome of Ambystoma mexicanum. Front Immunol. (2023) 14:1039274. doi: 10.3389/fimmu.2023.1039274
25
CharlemagneJ. Thymus independent anti-horse erythrocyte antibody response and suppressor T cells in the Mexican axolotl (Amphibia, Urodela, ambystoma mexicanum). Immunology. (1979) 36:643–8. https://pmc.ncbi.nlm.nih.gov/articles/PMC1457660/
26
FellahJSKerfournFDumayAMAubetGCharlemagneJ. Structure and diversity of the T-cell receptor alpha chain in the Mexican axolotl. Immunogenetics. (1997) 45:235–41. doi: 10.1007/s002510050198
27
FellahJSAndréSKerfournFGuerciADurandCAubetGet al. Structure, diversity and expression of the TCRdelta chains in the Mexican axolotl. Eur J Immunol. (2002) 32:1349–58. doi: 10.1002/1521-4141(200205)32:5<1349::AID-IMMU1349>3.0.CO;2-X
28
AndréSKerfournFAffaticatiPGuerciARavassardPFellahJS. Highly restricted diversity of TCR delta chains of the amphibian Mexican axolotl (Ambystoma mexicanum) in peripheral tissues. Eur J Immunol. (2007) 37:1621–33. doi: 10.1002/eji.200636375
29
SchloissnigSKawaguchiANowoshilowSFalconFOtsukiLTardivoPet al. The giant axolotl genome uncovers the evolution, scaling, and transcriptional control of complex gene loci. Proc Natl Acad Sci U S A. (2021) 118:e2017176118. doi: 10.1073/pnas.2017176118
30
HaireRNKitzan HaindfieldMKTurpenJBLitmanGW. Structure and diversity of T-lymphocyte antigen receptors alpha and gamma in Xenopus. Immunogenetics. (2002) 54:431–8. doi: 10.1007/s00251-002-0474-4
31
KubotaTWangJGöbelTWHockettRDCooperMDChenCH. Characterization of an avian (Gallus gallus domesticus) TCR alpha delta gene locus. J Immunol. (1999) 163:3858–66. doi: 10.4049/jimmunol.163.7.3858
32
AntonacciRMassariSLinguitiGCaputi JambrenghiAGiannicoFLefrancMPet al. Evolution of the T-cell receptor (TR) loci in the adaptive immune response: the tale of the TRG locus in mammals. Genes. (2020) 11:624. doi: 10.3390/genes11060624
33
LiuCZhangYLiuCCSchatzDG. Structural insights into the evolution of the RAG recombinase. Nat Rev Immunol. (2022) 22:353–70. doi: 10.1038/s41577-021-00628-6
34
BeddoeTChenZClementsCSElyLKBushellSRVivianJPet al. Antigen ligation triggers a conformational change within the constant domain of the alphabeta T cell receptor. Immunity. (2009) 30:777–88. doi: 10.1016/j.immuni.2009.03.018
35
BrazinKNMallisRJBoeszoermenyiAFengYYoshizawaARechePAet al. The T cell antigen receptor α transmembrane domain coordinates triggering through regulation of bilayer immersion and CD3 subunit associations. Immunity. (2018) 49:829–41.e6. doi: 10.1016/j.immuni.2018.09.007
36
BäckströmBTMiliaEPeterAJaureguiberryBBaldariCTPalmerE. A motif within the T cell receptor alpha chain constant region connecting peptide domain controls antigen responsiveness. Immunity. (1996) 5:437–47. doi: 10.1016/s1074-7613(00)80500-2
37
LefrancMPLefrancG. The T Cell Receptor FactsBook. Elsevier: London, UK (2001). p. 397.
38
PégorierPBertignacMChentliINguefack NgouneVFolchGJabado-MichaloudJet al. IMGT® biocuration and comparative study of the T cell receptor beta locus of veterinary species based on Homo sapiens TRB. Front Immunol. (2020) 11:821. doi: 10.3389/fimmu.2020.00821
39
BoudinotPNovasSJouneauLMondotSLefrancMPGrimholtUet al. Evolution of T cell receptor beta loci in salmonids. Front Immunol. (2023) 14:1238321. doi: 10.3389/fimmu.2023.1238321
40
KimSTToumaMTakeuchiKSunZYJDaveVPKappesDJet al. Distinctive CD3 heterodimeric ectodomain topologies maximize antigen-triggered activation of alpha beta T cell receptors. J Immunol. (2010) 185:2951–9. doi: 10.4049/jimmunol.1000732
41
PandeyPRRóżyckiBLipowskyRWeiklTR. Structural variability and concerted motions of the T cell receptor - CD3 complex. Elife. (2021) 10:e67195. doi: 10.7554/eLife.67195
42
GlusmanGRowenLLeeIBoysenCRoachJCSmitAFet al. Comparative genomics of the human and mouse T cell receptor loci. Immunity. (2001) 15:337–49. doi: 10.1016/S1074-7613(01)00200-X
43
SmeltyPMarchalCRenardRSinzelleLPolletNDunonDet al. Identification of the pre-T-cell receptor alpha chain in nonmammalian vertebrates challenges the structure-function of the molecule. Proc Natl Acad Sci U S A. (2010) 107:19991–6. doi: 10.1073/pnas.1010166107
44
ChienYHIwashimaMKaplanKBElliottJFDavisMM. A new T-cell receptor gene located within the alpha locus and expressed early in T-cell differentiation. Nature. (1987) 327:677–82. doi: 10.1038/327677a0
45
SatyanarayanaKHataSDevlinPRoncaroloMGDe VriesJESpitsHet al. Genomic organization of the human T-cell antigen-receptor alpha/delta locus. Proc Natl Acad Sci U S A. (1988) 85:8166–70. doi: 10.1073/pnas.85.21.8166
46
ConnelleyTKDegnanKLonghiCWMorrisonWI. Genomic analysis offers insights into the evolution of the bovine TRA/TRD locus. BMC Genomics. (2014) 15:994. doi: 10.1186/1471-2164-15-994
47
RastJPAndersonMKStrongSJLuerCLitmanRTLitmanGW. alpha, beta, gamma, and delta T cell antigen receptor genes arose early in vertebrate phylogeny. Immunity. (1997) 6:1–11. doi: 10.1016/S1074-7613(00)80237-X
48
SahaNROtaTLitmanGWHansenJParraZHsuEet al. Genome complexity in the coelacanth is reflected in its adaptive immune system. J Exp Zool B Mol Dev Evol. (2014) 322:438–63. doi: 10.1002/jez.b.22558
49
MariuzzaRAAgnihotriPOrbanJ. The structural basis of T-cell receptor (TCR) activation: An enduring enigma. J Biol Chem. (2020) 295:914–25. doi: 10.1016/S0021-9258(17)49904-2
50
BrazinKNMallisRJDasDKFengYHwangWWangJHet al. Structural features of the αβTCR mechanotransduction apparatus that promote pMHC discrimination. Front Immunol. (2015) 6:441. doi: 10.3389/fimmu.2015.00441
51
SchatzDGJiY. Recombination centres and the orchestration of V(D)J recombination. Nat Rev Immunol. (2011) 11:251–63. doi: 10.1038/nri2941
52
KelmMDarzentasFDarzentasNKotrovaMWesselsWBendigSet al. Dominant T-cell receptor delta rearrangements in B-cell precursor acute lymphoblastic leukemia: Leukemic markers or physiological γδ T repertoire? HemaSphere. (2023) 7:e948. doi: 10.1097/HS9.0000000000000948
53
CriscitielloMFOhtaYSaltisMMcKinneyECFlajnikMF. Evolutionarily conserved TCR binding sites, identification of T cells in primary lymphoid tissues, and surprising trans-rearrangements in nurse shark. J Immunol. (2010) 184:6950–60. doi: 10.4049/jimmunol.0902774
54
MorrisseyKAWegreckiMPraveenaTHansenVLBuLSivaramanKKet al. The molecular assembly of the marsupial γμ T cell receptor defines a third T cell lineage. Science. (2021) 371:1383–8. doi: 10.1126/science.abe7070
55
MorrisseyKASampsonJMRiveraMBuLHansenVLGemmellNJet al. Comparison of reptilian genomes reveals deletions associated with the natural loss of γδ T cells in squamates. J Immunol. (2022) 208:1960–7. doi: 10.4049/jimmunol.2101158
56
HaydayACVantouroutP. The innate biologies of adaptive antigen receptors. Annu Rev Immunol. (2020) 38:487–510. doi: 10.1146/annurev-immunol-102819-023144
57
FellahJSKerfournFCharlemagneJ. Evolution of T cell receptor genes. Extensive diversity of V beta families in the Mexican axolotl. J Immunol. (1994) 153:4539–45. doi: 10.4049/jimmunol.153.10.4539
58
KerfournFCharlemagneJFellahJS. The structure, rearrangement, and ontogenic expression of DB and JB gene segments of the Mexican axolotl T-cell antigen receptor beta chain (TCRB). Immunogenetics. (1996) 44:275–85. doi: 10.1007/s002510050124
59
FellahJSDurandCKerfournFCharlemagneJ. Complexity of the T cell receptor Cbeta isotypes in the Mexican axolotl: structure and diversity of the VDJCbeta3 and VDJCbeta4 chains. Eur J Immunol. (2001) 31:403–11. doi: 10.1002/1521-4141(200102)31:2<403::AID-IMMU403>3.0.CO;2-4
60
UdeyJABlombergB. Human lambda light chain locus: organization and DNA sequences of three genomic J regions. Immunogenetics. (1987) 25:63–70. doi: 10.1007/BF00768834
61
AntonacciRDi TommasoSLanaveCCribiuEPCiccareseSMassariS. Organization, structure and evolution of 41kb of genomic DNA spanning the D-J-C region of the sheep TRB locus. Mol Immunol. (2008) 45:493–509. doi: 10.1016/j.molimm.2007.05.023
62
AntonacciRGiannicoFCiccareseSMassariS. Genomic characteristics of the T cell receptor (TRB) locus in the rabbit (Oryctolagus cuniculus) revealed by comparative and phylogenetic analyses. Immunogenetics. (2014) 66:255–66. doi: 10.1007/s00251-013-0754-1
63
SasadaTToumaMChangHCClaytonLKWangJHReinherzEL. Involvement of the TCR Cbeta FG loop in thymic selection and T cell function. J Exp Med. (2002) 195:1419–31. doi: 10.1084/jem.20020119
64
SmithJJTimoshevskayaNTimoshevskiyVAKeinathMCHardyDVossSR. A chromosome-scale assembly of the axolotl genome. Genome Res. (2019) 29:317–24. doi: 10.1101/gr.241901.118
65
TimoshevskayaNVossSRLabiancaCNHighCRSmithJJ. Large-scale variation in single nucleotide polymorphism density within the laboratory axolotl (Ambystoma mexicanum). Dev Dyn. (2021) 250:822–37. doi: 10.1002/dvdy.257
66
WatsonCTGlanvilleJMarascoWA. The individual and population genetics of antibody immunity. Trends Immunol. (2017) 38:459–70. doi: 10.1016/j.it.2017.04.003
67
WatsonCTKosJTGibsonWSNewmanLDeikusGBusseCEet al. A comparison of immunoglobulin IGHV, IGHD and IGHJ genes in wild-derived and classical inbred mouse strains. Immunol Cell Biol. (2019) 97:888–901. doi: 10.1111/imcb.12288
68
Vázquez BernatNCorcoranMNowakIKadukMCastro DopicoXNarangSet al. Rhesus and cynomolgus macaque immunoglobulin heavy-chain genotyping yields comprehensive databases of germline VDJ alleles. Immunity. (2021) 54:355–66. doi: 10.1016/j.immuni.2020.12.018
69
LiuFLiJLinIYCYangXMaJChenYet al. The genome resequencing of TCR loci in Gallus gallus revealed their distinct evolutionary features in avians. ImmunoHorizons. (2020) 4:33–46. doi: 10.4049/immunohorizons.1900095
70
LynchMWalshB. The origins of genome architecture. New York, NY: Oxford University Press (2007). p. 389.
71
JoBSChoiSS. Introns: The functional benefits of introns in genomes. Genomics Inform. (2015) 13:112–8. doi: 10.5808/GI.2015.13.4.112
72
ComeronJMWillifordAKlimanRM. The Hill-Robertson effect: evolutionary consequences of weak selection and linkage in finite populations. Heredity (Edinb). (2008) 100:19–31. doi: 10.1038/sj.hdy.6801059
73
WangETSandbergRLuoSKhrebtukovaIZhangLMayrCet al. Alternative isoform regulation in human tissue transcriptomes. Nature. (2008) 456:470–6. doi: 10.1038/nature07509
74
BehlkeMLohEY. Alternative splicing of murine TCR beta transcripts. Nature. (1986) 322:379–82. doi: 10.1038/322379a0
75
NeiMGuXSitnikovaT. Evolution by the birth-and-death process in multigene families of the vertebrate immune system. Proc Natl Acad Sci U S A. (1997) 94:7799–806. doi: 10.1073/pnas.94.15.7799
76
von BoehmerH. Unique features of the pre-T-cell receptor alpha-chain: not just a surrogate. Nat Rev Immunol. (2005) 5:571–7. doi: 10.1038/nri1636
77
FahmyGHSicardRE. A role for effectors of cellular immunity in epimorphic regeneration of amphibian limbs. In Vivo. (2002) 16:179–84. https://europepmc.org/article/med/12182113
78
LeighNDDunlapGSJohnsonKMarianoROshiroRWongAYet al. Transcriptomic landscape of the blastema niche in regenerating adult axolotl limbs at single-cell resolution. Nat Commun. (2018) 9:5153. doi: 10.1038/s41467-018-07604-0
79
HuiSPShengDZSugimotoKGonzalez-RajalANakagawaSHesselsonDet al. Zebrafish regulatory T cells mediate organ-specific regenerative programs. Dev Cell. (2017) 43:659–72.e5. doi: 10.1016/j.devcel.2017.11.010
80
FranchiniABertolottiE. The thymus and tail regenerative capacity in Xenopus laevis tadpoles. Acta Histochem. (2012) 114:334–41. doi: 10.1016/j.acthis.2011.07.001
81
Bolaños-CastroLAWaltersHEGarcía VázquezROYunMH. Immunity in salamander regeneration: Where are we standing and where are we headed? Dev Dyn. (2021) 250:753–67. doi: 10.1002/dvdy.251
82
NCBI. Ambystoma mexicanum genome assembly UKY_AmexF1_1 (2025). Available online at: https://www.ncbi.nlm.nih.gov/datasets/genome/GCF_040938575.1/.
83
DobinADavisCASchlesingerFDrenkowJZaleskiCJhaSet al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. (2013) 29:15–21. doi: 10.1093/bioinformatics/bts635
84
ThorvaldsdottirHRobinsonJTMesirovJP. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Briefings Bioinf. (2013) 14:178–92. doi: 10.1093/bib/bbs017
85
IUIS Nomenclature Committee. T-cell receptor and immunoglobulin nomenclature sub-committee. Data Review and Naming of Human Immunoglobulin and T Cell Receptor Genes (2023). Available online at: https://wp-iuis.s3.eu-west-1.amazonaws.com/app/uploads/2024/10/04120932/Data.Review.and_.Naming.Dec2023.pdf (Accessed January 4, 2025 and June 20, 2025).
86
LefrancMP. IMGT, the international imMunoGeneTics information system. Cold Spring Harb Protoc. (2011) 2011:595–603. doi: 10.1101/pdb.top115
87
EddySR. Accelerated profile HMM searches. PloS Comput Biol. (2011) 7:e1002195. doi: 10.1371/journal.pcbi.1002195
Summary
Keywords
Ambystoma, T-cell receptor (TR), TRA/TRD locus, TRB locus, pseudogenes, amphibians
Citation
Pacheco-Olvera DL, Saint Remy-Hernández S, Godoy-Lozano EE, Téllez-Sosa J, Valdovinos-Torres H, Curiel-Quesada E, López-Macías C and Martínez-Barnetche J (2025) Genomic characterization of the T-cell receptor loci in Ambystoma mexicanum. Front. Immunol. 16:1656386. doi: 10.3389/fimmu.2025.1656386
Received
30 June 2025
Accepted
09 September 2025
Published
30 September 2025
Corrected
02 October 2025
Volume
16 - 2025
Edited by
Sofia Kossida, Université de Montpellier, France
Reviewed by
Anna E. Savage, University of Central Florida, United States
Long Ma, Zunyi Medical University, China
Updates
Copyright
© 2025 Pacheco-Olvera, Saint Remy-Hernández, Godoy-Lozano, Téllez-Sosa, Valdovinos-Torres, Curiel-Quesada, López-Macías and Martínez-Barnetche.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jesus Martínez-Barnetche, jmbarnet@insp.mx; Constantino López-Macías, constantino.lopez135@gmail.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.