Genomic characteristics of the T cell receptor (TRB) locus in the rabbit (Oryctolagus cuniculus) revealed by comparative and phylogenetic analyses

Squamate reptiles may have compensated for the lack of γδTCR with a duplication of the TRB locus

Squamate reptiles are amongst the most successful terrestrial vertebrate lineages, with over 10,000 species across a broad range of ecosystems. Despite their success, squamates are also amongst the least studied lineages immunologically. Recently, a universal lack of γδ T cells in squamates due to deletions of the genes encoding the T cell receptor (TCR) γ and δ chains was discovered. Here, we begin to address how the loss of γδ T cells may have impacted the evolution of the squamate immune system. Using the skink Tiliqua rugosa, we found that squamates have not significantly increased the complexity of conventional T cell receptor beta (TCRβ or TRB) chain V regions compared to that of the nearest living squamate relative, the tuatara, Sphenodon punctatus or other amniotes. Our analyses include a putative new TCR locus. This novel locus contains V, D, and J gene segments that undergo V(D)J recombination, albeit with a limited number of gene segments in most squamate species. Based on conserved residues, the predicted protein chain would be expected to form a heterodimer with TCRα. This new TCR locus appears to be derived from an ancient duplication of the TRB locus and is homologous to the recently described T cell receptor epsilon (TRE). TRE is absent from the genomes of the tuatara and all Archosaurs examined and appears squamate specific.

A Comprehensive Analysis of the Genomic and Expressed Repertoire of the T-Cell Receptor Beta Chain in Equus caballus

In this paper, we report a comprehensive and consistent annotation of the locus encoding the β-chain of the equine T-cell receptor (TRB), as inferred from recent genome assembly using bioinformatics tools. The horse TRB locus spans approximately 1 Mb, making it the largest locus among the mammalian species studied to date, with a significantly higher number of genes related to extensive duplicative events. In the region, 136 TRBV (belonging to 29 subgroups), 2 TRBD, 13 TRBJ, and 2 TRBC genes, were identified. The general genomic organization resembles that of other mammals, with a V cluster of 135 TRBV genes located upstream of two in-tandem aligned TRBD-J-C clusters and an inverted TRBV gene at the 3' end of the last TRBC gene. However, the horse b-chain repertoire would be affected by a high number of non-functional TRBV genes. Thus, we queried a transcriptomic dataset derived from splenic tissue of a healthy adult horse, using each TRBJ gene as a probe to analyze clonotypes encompassing the V(D)J junction. This analysis provided insights into the usage of the TRBV, TRBD, and TRBJ genes and the variability of the non-germline-encoded CDR3. Our results clearly demonstrated that the horse β-chain constitutes a complex level of variability, broadly like that described in other mammalian species.

Profiling the T Cell Receptor Alpha/Delta Locus in Salmonids

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

The expansion of the TRB and TRG genes in domestic goats (Capra hircus) is characteristic of the ruminant species

BACKGROUND

Goats (Capra hircus), one of the first domesticated species, are economically important for milk and meat production, and their broad geographical distribution reflects their successful adaptation to diverse environmental conditions. Despite the relevance of this species, the genetic research on the goat traits is limited compared to other domestic species. Thanks to the latest goat reference genomic sequence (ARS1), which is considered to be one of the most continuous assemblies in livestock, we deduced the genomic structure of the T cell receptor beta (TRB) and gamma (TRG) loci in this ruminant species.

RESULTS

Our analyses revealed that although the organization of the goat TRB locus is broadly similar to that of the other artiodactyl species, with three in-tandem D-J-C clusters located at the 3' end, a complex and extensive series of duplications have occurred in the V genes at the 5' end, leading to a marked expansion in the number of the TRBV genes. This phenomenon appears to be a feature of the ruminant lineage since similar gene expansions have also occurred in sheep and cattle. Likewise, the general organization of the goat TRG genes is typical of ruminant species studied so far, with two paralogous TRG loci, TRG1 and TRG2, located in two distinct and distant positions on the same chromosome as result of a split in the ancestral locus. Each TRG locus consists of reiterated V-J-J-C cassettes, with the goat TRG2 containing an additional cassette relative to the corresponding sheep and cattle loci.

CONCLUSIONS

Taken together, these findings demonstrate that strong evolutionary pressures in the ruminant lineage have selected for the development of enlarged sets of TRB and TRG genes that contribute to a diverse T cell receptor repertoire. However, differences observed among the goat, sheep and cattle TRB and TRG genes indicate that distinct evolutionary histories, with independent expansions and/or contractions, have also affected each ruminant species.

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

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

IMGT® Biocuration and Comparative Study of the T Cell Receptor Beta Locus of Veterinary Species Based on Homo sapiens TRB

IMGT®, the international ImMunoGeneTics information system® is the global reference in immunogenetics and immunoinformatics. By its creation in 1989 by Marie-Paule Lefranc (Université de Montpellier and CNRS), IMGT® marked the advent of immunoinformatics, which emerged at the interface between immunogenetics and bioinformatics. IMGT® is specialized in the immunoglobulins (IG) or antibodies, T cell receptors (TR), major histocompatibility (MH), and proteins of the IgSF and MhSF superfamilies. T cell receptors are divided into two groups, αβ and γδ TR, which express distinct TR containing either α and β, or γ and δ chains, respectively. The TRβ locus (TRB) was recently described and annotated by IMGT® biocurators for several veterinary species, i.e., cat (Felis catus), dog (Canis lupus familiaris), ferret (Mustela putorius furo), pig (Sus scrofa), rabbit (Oryctolagus cuniculus), rhesus monkey (Macaca mulatta), and sheep (Ovis aries). The aim of the present study is to compare the genes of the TRB locus among these different veterinary species based on Homo sapiens. The results reveal that there are similarities but also differences including the number of genes by subgroup which may demonstrate duplications and/or deletions during evolution.

Topology and expressed repertoire of the Felis catus T cell receptor loci

Background

The domestic cat (Felis catus) is an important companion animal and is used as a large animal model for human disease. However, the comprehensive study of adaptive immunity in this species is hampered by the lack of data on lymphocyte antigen receptor genes and usage. The objectives of this study were to annotate the feline T cell receptor (TR) loci and to characterize the expressed repertoire in lymphoid organs of normal cats using high-throughput sequencing.

Results

The Felis catus TRG locus contains 30 genes: 12 TRGV, 12 TRGJ and 6 TRGC, the TRB locus contains 48 genes: 33 TRBV, 2 TRBD, 11 TRBJ, 2 TRBC, the TRD locus contains 19 genes: 11 TRDV, 2 TRDD, 5 TRDJ, 1 TRDC, and the TRA locus contains 127 genes: 62 TRAV, 64 TRAJ, 1 TRAC. Functional feline V genes form monophyletic clades with their orthologs, and clustering of multimember subgroups frequently occurs in V genes located at the 5′ end of TR loci. Recombination signal (RS) sequences of the heptamer and nonamer of functional V and J genes are highly conserved. Analysis of the TRG expressed repertoire showed preferential intra-cassette over inter-cassette rearrangements and dominant usage of the TRGV2–1 and TRGJ1–2 genes. The usage of TRBV genes showed minor bias but TRBJ genes of the second J-C-cluster were more commonly rearranged than TRBJ genes of the first cluster. The TRA/TRD V genes almost exclusively rearranged to J genes within their locus. The TRAV/TRAJ gene usage was relatively balanced while the TRD repertoire was dominated by TRDJ3.

Conclusions

This is the first description of all TR loci in the cat. The genomic organization of feline TR loci was similar to that of previously described jawed vertebrates (gnathostomata) and is compatible with the birth-and-death model of evolution. The large-scale characterization of feline TR genes provides comprehensive baseline data on immune repertoires in healthy cats and will facilitate the development of improved reagents for the diagnosis of lymphoproliferative diseases in cats. In addition, these data might benefit studies using cats as a large animal model for human disease.

Characterization of the ferret TRB locus guided by V, D, J, and C gene expression analysis

The domestic ferret, Mustela putorius furo, is an important mammalian animal model to study human respiratory infection. However, insufficient genomic annotation hampers detailed studies of ferret T cell responses. In this study, we analyzed the published T cell receptor beta (TRB) locus and performed high-throughput sequencing (HTS) of peripheral blood of four healthy adult ferrets to identify expressed V, D, J, and C genes. The HTS data is used as a guide to manually curate the expressed V, D, J, and C genes. The ferret locus appears to be most similar to that of the dog. Like other mammalian TRB loci, the ferret TRB locus contains a library of variable genes located upstream of two D-J-C gene clusters, followed by a (in the ferret non-functional) V gene with an inverted transcriptional orientation. All TRB genes (expressed or not) reported here have been approved by the IMGT/WHO-IUIS nomenclature committee.

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

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

The Camel Adaptive Immune Receptors Repertoire as a Singular Example of Structural and Functional Genomics

The adaptive immune receptors repertoire is highly plastic, with its ability to produce antigen-binding molecules and select those with high affinity for their antigen. Species have developed diverse genetic and structural strategies to create their respective repertoires required for their survival in the different environments. Camelids, until now, considered as a case of evolutionary innovation because of their only heavy-chain antibodies, represent a new mammalian model particularly useful for understanding the role of diversity in the immune system function. Here, we review the structural and functional characteristics and the current status of the genomic organization of camel immunoglobulins (IG) or antibodies, α/ß and γ/δ T cell receptors (TR), and major histocompatibility complex (MHC). In camelid humoral response, in addition to the conventional antibodies, there are IG with “only-heavy-chain” (no light chain, and two identical heavy gamma chains lacking CH1 and with a VH domain designated as VHH). The unique features of these VHH offer advantages in biotechnology and for clinical applications. The TRG and TRD rearranged variable domains of Camelus dromedarius (Arabian camel) display somatic hypermutation (SHM), increasing the intrinsic structural stability in the γ/δ heterodimer and influencing the affinity maturation to a given antigen similar to immunoglobulin genes. The SHM increases the dromedary γ/δ repertoire diversity. In Camelus genus, the general structural organization of the TRB locus is similar to that of the other artiodactyl species, with a pool of TRBV genes positioned at the 5’ end of three in tandem D-J-C clusters, followed by a single TRBV gene with an inverted transcriptional orientation located at the 3’ end. At the difference of TRG and TRD, the diversity of the TRB variable domains is not shaped by SHM and depends from the classical combinatorial and junctional diversity. The MHC locus is located on chromosome 20 in Camelus dromedarius. Cytogenetic and comparative whole genome analyses revealed the order of the three major regions “Centromere-ClassII-ClassIII-ClassI”. Unexpectedly low extent of polymorphisms and haplotypes was observed in all Old World camels despite different geographic origins.

The T cell receptor (TRA) locus in the rabbit (Oryctolagus cuniculus): Genomic features and consequences for invariant T cells

The rabbit has been widely used in immunology and infectiology. Rabbit immunoglobulins have been extensively studied, leading to the discovery of their idiotypes, allotypic diversity, and of the diversification of the primary repertoire by hyperconversion. Much less is known about rabbit T cell receptors (TR), especially TRA. This isotype is particularly important for innate-like T cells, which typically express invariant TRA (iTRA). The presence of such cells in the rabbit remains an enigma. Rabbit NKT cells seem to be very rare, and lagomorphs lack MAIT cells. TRAV1, the variable gene expressed in the iTRA of these cells across most mammals, and MR1, the MH1-like receptor that present riboflavin derivatives to MAIT cells, are missing in rabbit. An alternative iTRA has been identified, that may be expressed by new innate-like T cells. To facilitate TRA repertoire analyses in rabbit, we report here a full description of TRA and TRD loci and a subgroup definition based on IMGT® classification. Rabbit TRA rearrangements follow the same temporal pattern that is observed in mouse and human. Rare transcripts expressing TRDV/TRDD/TRDJ rearrangements spliced to TRAC were detected. TRA and TRD genes have been made available in IMGT and IMGT/HighV-QUEST, allowing easy analysis of TRA/TRD RepSeq.

Comparative Analysis of the TRB Locus in the Camelus Genus

T cells can be separated into two major subsets based on the heterodimer that forms their T cell receptors. αβ T cells have receptors consisting of α and β chains, while γδ T cells are composed of γ and δ chains. αβ T cells play an essential role within the adaptive immune responses against pathogens. The recent genomic characterization of the Camelus dromedarius T cell receptor β (TRB) locus has allowed us to infer the structure of this locus from the draft genome sequences of its wild and domestic Bactrian congeners, Camelus ferus and Camelus bactrianus. The general structural organization of the wild and domestic Bactrian TRB locus is similar to that of the dromedary, with a pool of TRBV genes positioned at the 5′ end of D-J-C clusters, followed by a single TRBV gene located at the 3′ end with an inverted transcriptional orientation. Despite the fragmented nature of the assemblies, comparative genomics reveals the existence of a perfect co-linearity between the three Old World camel TRB genomic sequences, which enables the transfer of information from one sequence to another and the filling of gaps in the genomic sequences. A virtual camelid TRB locus is hypothesized with the presence of 33 TRBV genes distributed in 26 subgroups. Likewise, in the artiodactyl species, three in-tandem D-J-C clusters, each composed of one TRBD gene, six or seven TRBJ genes, and one TRBC gene, are placed at the 3′ end of the locus. As reported in the ruminant species, a group of four functional TRY genes at the 5′ end and only one gene at the 3′ end, complete the camelid TRB locus. Although the gene content is similar, differences are observed in the TRBV functional repertoire, and genes that are functional in one species are pseudogenes in the other species. Hence, variations in the functional repertoire between dromedary, wild and domestic Bactrian camels, rather than differences in the gene content, may represent the molecular basis explaining the disparity in the TRB repertoire between the Camelus species. Finally, our data contribute to the knowledge about the evolutionary history of Old World camelids.

Overview of the Germline and Expressed Repertoires of the TRB Genes in Sus scrofa

The α/β T cell receptor (TR) is a complex heterodimer that recognizes antigenic peptides and binds to major histocompatibility complex (MH) molecules. Both α and β chains are encoded by different genes localized on two distinct chromosomal loci: TRA and TRB. The present study employed the recent release of the swine genome assembly to define the genomic organization of the TRB locus. According to the sequencing data, the pig TRB locus spans approximately 400 kb of genomic DNA and consists of 38 TRBV genes belonging to 24 subgroups located upstream of three in tandem TRBD-J-C clusters, which are followed by a TRBV gene in an inverted transcriptional orientation. Comparative analysis confirms that the general organization of the TRB locus is similar among mammalian species, but the number of germline TRBV genes varies greatly even between species belonging to the same order, determining the diversity and specificity of the immune response. However, sequence analysis of the TRB locus also suggests the presence of blocks of conserved homology in the genomic region across mammals. Furthermore, by analysing a public cDNA collection, we identified the usage pattern of the TRBV, TRBD, and TRBJ genes in the adult pig TRB repertoire, and we noted that the expressed TRBV repertoire seems to be broader and more diverse than the germline repertoire, in line with the presence of a high level of TRBV gene polymorphisms. Because the nucleotide differences seems to be principally concentrated in the CDR2 region, it is reasonable to presume that most T cell β-chain diversity can be related to polymorphisms in pig MH molecules. Domestic pigs represent a valuable animal model as they are even more anatomically, genetically and physiologically similar to humans than are mice. Therefore, present knowledge on the genomic organization of the pig TRB locus allows the collection of increased information on the basic aspects of the porcine immune system and contributes to filling the gaps left by rodent models.

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

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

The occurrence of three D-J-C clusters within the dromedary TRB locus highlights a shared evolution in Tylopoda, Ruminantia and Suina

The αβ T cells are important components of the adaptive immune system and can recognize a vast array of peptides presented by MHC molecules. The ability of these T cells to recognize the complex depends on the diversity of the αβ TR, which is generated by a recombination of specific Variable, Diversity and Joining genes for the β chain, and Variable and Joining genes for the α chain. In this study, we analysed the genomic structure and the gene content of the TRB locus in Camelus dromedarius, which is a species belonging to the Tylopoda suborder. The most noteworthy result is the presence of three in tandem TRBD-J-C clusters in the dromedary TRB locus, which is similar to clusters found in sheep, cattle and pigs and suggests a common duplication event occurred prior to the Tylopoda/Ruminantia/Suina divergence. Conversely, a significant contraction of the dromedary TRBV genes, which was previously found in the TRG and TRD loci, was observed with respect to the other artiodactyl species.

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

Recently, many immune-related genes have been extensively studied in ducks, but relatively little is known about their TCR genes. Here, we determined the germline and expressed repertoire of TCR genes in White Peking duck. The genomic organization of the duck TCRα/δ, TCRγ and unconventional TCRδ2 loci are highly conserved with their counterparts in mammals or chickens. By contrast, the duck TCRβ locus is organized in an unusual pattern, (Vβ)_n-Dβ-(Jβ)₂-Cβ1-(Jβ)₄-Cβ2, which differs from the tandem-aligned clusters in mammals or the translocon organization in some teleosts. Excluding the first exon encoding the immunoglobulin domain, the subsequent exons of the two Cβ show significant diversity in nucleotide sequence and exon structure. Based on the nucleotide sequence identity, 49 Vα, 30 Vδ, 13 Vβ and 15 Vγ unique gene segments are classified into 3 Vα, 5 Vδ, 4 Vβ and 6 Vγ subgroups, respectively. Phylogenetic analyses revealed that most duck V subgroups, excluding Vβ1, Vγ5 and Vγ6, have closely related orthologues in chicken. The coding joints of all cDNA clones demonstrate conserved mechanisms that are used to increase junctional diversity. Collectively, these data provide insight into the evolution of TCRs in vertebrates and improve our understanding of the avian immune system.

An overview of the lagomorph immune system and its genetic diversity

Our knowledge of the lagomorph immune system remains largely based upon studies of the European rabbit (Oryctolagus cuniculus), a major model for studies of immunology. Two important and devastating viral diseases, rabbit hemorrhagic disease and myxomatosis, are affecting European rabbit populations. In this context, we discuss the genetic diversity of the European rabbit immune system and extend to available information about other lagomorphs. Regarding innate immunity, we review the most recent advances in identifying interleukins, chemokines and chemokine receptors, Toll-like receptors, antiviral proteins (RIG-I and Trim5), and the genes encoding fucosyltransferases that are utilized by rabbit hemorrhagic disease virus as a portal for invading host respiratory and gut epithelial cells. Evolutionary studies showed that several genes of innate immunity are evolving by strong natural selection. Studies of the leporid CCR5 gene revealed a very dramatic change unique in mammals at the second extracellular loop of CCR5 resulting from a gene conversion event with the paralogous CCR2. For the adaptive immune system, we review genetic diversity at the loci encoding antibody variable and constant regions, the major histocompatibility complex (RLA) and T cells. Studies of IGHV and IGKC genes expressed in leporids are two of the few examples of trans-species polymorphism observed outside of the major histocompatibility complex. In addition, we review some endogenous viruses of lagomorph genomes, the importance of the European rabbit as a model for human disease studies, and the anticipated role of next-generation sequencing in extending knowledge of lagomorph immune systems and their evolution.