Sequencing of modern Lepus VDJ genes shows that the usage of VHn genes has been retained in both Oryctolagus and Lepus that diverged 12 million years ago

The IgA of hares (Lepus sp.) and rabbit confirms that the leporids IgA explosion is old and reveals a new case of trans-species polymorphism

Background

Immunoglobulin A (IgA) is the mammalian mucosal antibody, providing an important line of defense against pathogens. With 15 IgA subclasses, the European rabbit has an extremely complex IgA system, strikingly more complex than most other mammals, which have only one IgA or, in the case of hominoids, two IgA subclasses. Similar to the two hominoid primate IGHA genes, the expansion of the rabbit IGHA genes appears to have begun in an ancestral lagomorph since multiple IgA copies were found by Southern blot analysis for the genera Sylvilagus, Lepus, and Ochotona.

Results

To gain a better insight into the extraordinary lagomorph IgA evolution, we sequenced, for the first time, expressed IgA genes for two Lepus species, L. europaeus and L. granatensis. These were aligned with the 15 rabbit IgA isotypes, and evolutionary analyses were conducted. The obtained phylogenetic tree shows that the Lepus IgA sequences cluster with and among the rabbit IgA isotypes, and the interspecies and intraspecies nucleotide genetic distances are similar. A comparison of the amino acid sequences of the Lepus and rabbit IgA confirms that there are two trans-species polymorphisms and that the rabbit and Lepus sequences share a common genetic pool. In fact, the main differences between the studied leporids IgAs reside in the characteristics of the hinge region.

Conclusion

The Lepus IgA sequences we have obtained strongly suggest that the great expansion of the leporid IGHA genes occurred in a common ancestral species and was then maintained in the descendants. A strong selective pressure caused the extraordinary expansion of the IGHA genes but then subsided, leading to the maintenance of the acquired polymorphisms in the descendants, with little subsequent divergence. This is a unique evolutionary pattern in which an ancient gene expansion has been maintained for approximately 18 million years.

Evolutionary analyses of polymeric immunoglobulin receptor (pIgR) in the mammals reveals an outstanding mutation rate in the lagomorphs

BACKGROUND

The transcytosis of polymeric immunoglobulins, IgA and IgM, across the epithelial barrier to the luminal side of mucosal tissues is mediated by the polymeric immunoglobulin receptor (pIgR). At the luminal side the extracellular ligand binding region of pIgR, the secretory component (SC), is cleaved and released bound to dimeric IgA (dIgA), protecting it from proteolytic degradation, or in free form, protecting the mucosa form pathogens attacks. The pIgR was first cloned for rabbit in early 1980's and since then has been described for all vertebrates, from fish to mammals. The existence of more than one functional pIgR alternative-spliced variant in the European rabbit, the complete pIgR as other mammals and a shorter pIgR lacking two SC exons, raised the question whether other lagomorphs share the same characteristics and how has the PIGR gene evolved in these mammals.

RESULTS

To investigate these questions, we sequenced expressed pIgR genes for other leporid genus, Lepus spp., and obtained and aligned pIgR sequences from representative species of all mammalian orders. The obtained mammalian phylogeny, as well as the Bayesian inference of evolutionary rates and genetic distances, show that Lagomorpha pIgR is evolving at a higher substitution rate. Codon-based analyses of positive selection show that mammalian pIgR is evolving under strong positive selection, with strong incidence in the domains excised from the rabbit short pIgR isoform. We further confirmed that the hares also express the two rabbit pIgR isoforms.

CONCLUSIONS

The Lagomorpha pIgR unique evolutionary pattern may reflect a group specific adaptation. The pIgR evolution may be linked to the unusual expansion of IgA genes observed in lagomorphs, or to neofunctionalization in this group. Further studies are necessary to clarify the driving forces behind the unique lagomorph pIgR evolution.

Sequencing of VDJ genes in Lepus americanus confirms a correlation between VHn expression and the leporid species continent of origin

Leporid VH genes used in the generation of their primary antibody repertoire exhibit highly divergent lineages. For the European rabbit (Oryctolagus cuniculus) four VHa lineages have been described, the a1, a2, a3 and a4. Hares (Lepus spp.) and cottontail (Sylvilagus floridanus) express one VHa lineage each, the a2L and the a5, respectively, along with a more ancient lineage, the Lepus spp. sL and S. floridanus sS. Both the European rabbit and the Lepus europaeus use a third lineage, VHn, in a low proportion of their VDJ rearrangements. The VHn genes are a conserved ancestral polymorphism that is being maintained in the leporid genome.Their usage in a low proportion of VDJ rearrangements by both European rabbit and L. europaeus but not S. floridanus has been argued to be a remnant of an ancient European leporid immunologic response to pathogens. To address this hypothesis, in this study we sequenced VDJ rearranged genes for another North American leporid, L. americanus. Our results show that L. americanus expressed these genes less frequently and in a highly modified fashion compared to the European Lepus species. Our results suggest that the American leporid species use a different VH repertoire than the European species which may be related with an immune adaptation to different environmental conditions, such as different pathogenic agents.

Identification of a new European rabbit IgA with a serine-rich hinge region

In mammals, the most striking IgA system belongs to Lagomorpha. Indeed, 14 IgA subclasses have been identified in European rabbits, 11 of which are expressed. In contrast, most other mammals have only one IgA, or in the case of hominoids, two IgA subclasses. Characteristic features of the mammalian IgA subclasses are the length and amino acid sequence of their hinge regions, which are often rich in Pro, Ser and Thr residues and may also carry Cys residues. Here, we describe a new IgA that was expressed in New Zealand White domestic rabbits of IGHVa1 allotype. This IgA has an extended hinge region containing an intriguing stretch of nine consecutive Ser residues and no Pro or Thr residues, a motif exclusive to this new rabbit IgA. Considering the amino acid properties, this hinge motif may present some advantage over the common IgA hinge by affording novel functional capabilities. We also sequenced for the first time the IgA14 CH2 and CH3 domains and showed that IgA14 and IgA3 are expressed.

Computational Strategies for Dissecting the High-Dimensional Complexity of Adaptive Immune Repertoires

The adaptive immune system recognizes antigens via an immense array of antigen-binding antibodies and T-cell receptors, the immune repertoire. The interrogation of immune repertoires is of high relevance for understanding the adaptive immune response in disease and infection (e.g., autoimmunity, cancer, HIV). Adaptive immune receptor repertoire sequencing (AIRR-seq) has driven the quantitative and molecular-level profiling of immune repertoires, thereby revealing the high-dimensional complexity of the immune receptor sequence landscape. Several methods for the computational and statistical analysis of large-scale AIRR-seq data have been developed to resolve immune repertoire complexity and to understand the dynamics of adaptive immunity. Here, we review the current research on (i) diversity, (ii) clustering and network, (iii) phylogenetic, and (iv) machine learning methods applied to dissect, quantify, and compare the architecture, evolution, and specificity of immune repertoires. We summarize outstanding questions in computational immunology and propose future directions for systems immunology toward coupling AIRR-seq with the computational discovery of immunotherapeutics, vaccines, and immunodiagnostics.

The remnant of the European rabbit (Oryctolagus cuniculus) IgD gene

Although IgD first appeared, along with IgM, in the cartilaginous fishes and has been retained throughout subsequent vertebrate evolution, it has been lost in a diverse group of vertebrate species. We previously showed that, unlike vertebrates that express IgD, the rabbit lacks an IgD (Cδ) gene within 13.5 kb downstream of the IgM gene. We report here that, by conducting BLAST searches of rabbit Ig heavy chain genomic DNA with known mammalian IgD exons, we identified the remnant of the rabbit Cδ gene approximately 21 kb downstream of the IgM gene. The remnant Cδ locus lacks the δCH1 and hinge exons, but contains truncated δCH2 and δCH3 exons, as well as largely intact, but non-functional, secretory and transmembrane exons. In addition, we report that the Cδ gene probably became non-functional in leporids at least prior to the divergence of rabbits and hares ~12 million years ago.

Evolutionary study of leporid CD4 reveals a hotspot of genetic variability within the D2 domain

CD4 is the major receptor on T helper cells involved in the uptake of human immunodeficiency virus (HIV) or simian immunodeficiency virus (SIV) into their host cells. Evolutionary studies of CD4 in primates revealed signatures of positive selection in the D1 domain that interacts with primate exogenous lentivirus gp120 proteins. Here, we studied the evolution of CD4 in lagomorphs by comparing sequences obtained for the genera Oryctolagus, Sylvilagus, Lepus, and Ochotona. Our results reveal an overall higher divergence in lagomorphs compared to primates with highest divergence in the D2 domain. A detailed analysis of a small fragment of 33 nucleotides coding for amino acids 169 to 179 in the D2 domain showed dramatic amino acid alterations with a dN/dS value of 3.2 for lagomorphs, suggesting that CD4 is under strong positive selection in this particular region. Within each leporid genus, no significant amino acid changes were observed for the D2 domain which indicates that the genetic differentiation occurred in the ancestor of each genus before the species radiation. The rabbit endogenous lentivirus type K (RELIK) found in leporids shares high structural similarity with HIV which suggests a possible interaction between RELIK and CD4. The presence of RELIK in the studied leporids, the high structural similarity to modern-day exogenous lentiviruses and the absence of exogenous lentiviruses in leporids, allows us to hypothesize that this endogenous retrovirus, that was most probably exogenous in the past, drove the divergent evolution of leporid CD4.

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.

Genetic characterization of interleukins (IL-1α, IL-1β, IL-2, IL-4, IL-8, IL-10, IL-12A, IL-12B, IL-15 and IL-18) with relevant biological roles in lagomorphs

ILs, as essential innate immune modulators, are involved in an array of biological processes. In the European rabbit (Oryctolagus cuniculus) IL-1α, IL-1β, IL-2, IL-4, IL-8, IL-10, IL-12A, IL-12B, IL-15 and IL-18 have been implicated in inflammatory processes and in the immune response against rabbit hemorrhagic disease virus and myxoma virus infections. In this study we characterized these ILs in six Lagomorpha species (European rabbit, pygmy rabbit, two cottontail rabbit species, European brown hare and American pika). Overall, these ILs are conserved between lagomorphs, including in their exon/intron structure. Most differences were observed between leporids and American pika. Indeed, when comparing both, some relevant differences were observed in American pika, such as the location of the stop codon in IL-1α and IL-2, the existence of a different transcript in IL8 and the number of cysteine residues in IL-1β. Changes at N-glycosylation motifs were also detected in IL-1, IL-10, IL-12B and IL-15. IL-1α is the protein that presents the highest evolutionary distances, which is in contrast to IL-12A where the distances between lagomorphs are the lowest. For all these ILs, sequences of human and European rabbit are more closely related than between human and mouse or European rabbit and mouse.

Genetic diversity comparison of the DQA gene in European rabbit (Oryctolagus cuniculus) populations

The European rabbit (Oryctolagus cuniculus) natural populations within the species native region, the Iberian Peninsula, are considered a reservoir of genetic diversity. Indeed, the Iberia was a Pleistocene refuge to the species and currently two subspecies are found in the peninsula (Oryctolagus cuniculus cuniculus and Oryctolagus cuniculus algirus). The genes of the major histocompatibility complex (MHC) have been substantially studied in wild populations due to their exceptional variability, believed to be pathogen driven. They play an important function as part of the adaptive immune system affecting the individual fitness and population viability. In this study, the MHC variability was assessed by analysing the exon 2 of the DQA gene in several European rabbit populations from Portugal, Spain and France and in domestic breeds. Twenty-eight DQA alleles were detected, among which 18 are described for the first time. The Iberian rabbit populations are well differentiated from the French population and domestic breeds. The Iberian populations retained the higher allelic diversity with the domestic breeds harbouring the lowest; in contrast, the DQA nucleotide diversity was higher in the French population. Signatures of positive selection were detected in four codons which are putative peptide-binding sites and have been previously detected in other mammals. The evolutionary relationships showed instances of trans-species polymorphism. Overall, our results suggest that the DQA in European rabbits is evolving under selection and genetic drift.

Endogenization of mouse mammary tumor virus (MMTV)-like elements in genomes of pikas (Ochotona sp.)

Despite the finding in European rabbit and other leporid genomes of the first ever described endogenous lentivirus and of a European rabbit exclusive endogenous gammaretrovirus, until now no exogenous retroviruses have been isolated in Lagomorpha species. Nevertheless, looking for the presence of endogenous retroviruses (ERVs) in the species genomes could lead to the discovery of retroviral lineages yet to be found in Lagomorpha. Different mammalian genomes harbor endogenous viral sequences phylogenetically close to the betaretrovirus mouse mammary tumor virus (MMTV), propelling us to look for such retroviral "fossil" in American pika (Ochotona princeps) and European rabbit (Oryctolagus cuniculus) genomes. By performing genomic mining using MMTV gag and LTR as query sequences, we found that such viral elements were absent from the European rabbit genome. Oppositely, significant matches were found in American pika, and more importantly, a nearly complete MMTV-like virus (Pika-BERV) was identified. Using Pika-BERV gag and LTR as templates, we found similar sequences endogenized in different pika (Ochotona sp.) species. The orthology of the LTR flanking region between some pika species supported shared ancestry of specific endogenous betaretroviruses, while in other pika species similar sequences, but not orthologous, should have resulted from independent insertions. Our study supports the possible existence of infecting exogenous betaretroviruses for a long term, after the divergence of Ochotonidae from Leporidae, but yet to be identified.

Leporid immunoglobulin G shows evidence of strong selective pressure on the hinge and CH3 domains

Immunoglobulin G (IgG) is the predominant serum immunoglobulin and has the longest serum half-life of all the antibody classes. The European rabbit IgG has been of significant importance in immunological research, and is therefore well characterized. However, the IgG of other leporids has been disregarded. To evaluate the evolution of this gene in leporids, we sequenced the complete IGHG for six other genera: Bunolagus, Brachylagus, Lepus, Pentalagus, Romerolagus and Sylvilagus. The newly sequenced leporid IGHG gene has an organization and structure similar to that of the European rabbit IgG. A gradient in leporid IgG constant domain diversity was observed, with the CH1 being the most conserved and the CH3 the most variable domain. Positive selection was found to be acting on all constant domains, but with a greater incidence in the CH3 domain, where a cluster of three positively selected sites was identified. In the hinge region, only three polymorphic positions were observed. The same hinge length was observed for all leporids. Unlike the variation observed for the European rabbit, all 11 Lepus species studied share exactly the same hinge motif, suggesting its maintenance as a result of an advantageous structure or conformation.

Field and experimental data indicate that the eastern cottontail (Sylvilagus floridanus) is susceptible to infection with European brown hare syndrome (EBHS) virus and not with rabbit haemorrhagic disease (RHD) virus

The eastern cottontail (Sylvilagus floridanus) is an American lagomorph. In 1966, it was introduced to Italy, where it is currently widespread. Its ecological niche is similar to those of native rabbits and hares and increasing overlap in distribution brings these species into ever closer contact. Therefore, cottontails are at risk of infection with the two lagoviruses endemically present in Italy: Rabbit Haemorrhagic Disease virus (RHDV) and European Brown Hare Syndrome Virus (EBHSV). To verify the susceptibility of Sylvilagus to these viruses, we analyzed 471 sera and 108 individuals from cottontail populations in 9 provinces of north-central Italy from 1999 to 2012. In total, 15–20% of the cottontails tested seropositive for EBHSV; most titres were low, but some were as high as 1/1280. All the cottontails virologically tested for RHDV and EBHSV were negative with the exception of one individual found dead with hares during a natural EBHS outbreak in December 2009. The cottontail and the hares showed typical EBHS lesions, and the EBHSV strain identified was the same in both species (99.9% identity). To experimentally confirm the diagnosis, we performed two trials in which we infected cottontails with both EBHSV and RHDV. One out of four cottontails infected with EBHSV died of an EBHS-like disease, and the three surviving animals developed high EBHSV antibody titres. In contrast, neither mortality nor seroconversion was detected after infection with RHDV. Taken together, these results suggest that Sylvilagus is susceptible to EBHSV infection, which occasionally evolves to EBHS-like disease; the eastern cottontail could therefore be considered a “spill over” or “dead end” host for EBHSV unless further evidence is found to confirm that it plays an active role in the epidemiology of EBHSV.

Sequencing of Sylvilagus VDJ genes reveals a new VHa allelic lineage and shows that ancient VH lineages were retained differently in leporids

Antigen recognition by immunoglobulins depends upon initial rearrangements of heavy chain V, D, and J genes. In leporids, a unique system exists for the VH genes usage that exhibit highly divergent lineages: the VHa allotypes, the Lepus sL lineage and the VHn genes. For the European rabbit (Oryctolagus cuniculus), four VHa lineages have been described, the a1, a2, a3 and a4. For hares (Lepus sp.), one VHa lineage was described, the a2L, as well as a more ancient sL lineage. Both genera use the VHn genes in a low frequency of their VDJ rearrangements. To address the hypothesis that the VH specificities could be associated with different environments, we sequenced VDJ genes from a third leporid genus, Sylvilagus. We found a fifth and equally divergent VHa lineage, the a5, and an ancient lineage, the sS, related to the hares' sL, but failed to obtain VHn genes. These results show that the studied leporids employ different VH lineages in the generation of the antibody repertoire, suggesting that the leporid VH genes are subject to strong selective pressure likely imposed by specific pathogens.

An intriguing shift occurs in the novel protein phosphatase 1 binding partner, TCTEX1D4: evidence of positive selection in a pika model

T-complex testis expressed protein 1 domain containing 4 (TCTEX1D4) contains the canonical phosphoprotein phosphatase 1 (PPP1) binding motif, composed by the amino acid sequence RVSF. We identified and validated the binding of TCTEX1D4 to PPP1 and demonstrated that indeed this protein is a novel PPP1 interacting protein. Analyses of twenty-one mammalian species available in public databases and seven Lagomorpha sequences obtained in this work showed that the PPP1 binding motif ₉₀RVSF₉₃ is present in all of them and is flanked by a palindromic sequence, PLGS, except in three species of pikas (Ochotona princeps, O. dauurica and O. pusilla). Furthermore, for the Ochotona species an extra glycosylation site, motif ₉₆NLS₉₈, and the loss of the palindromic sequence were observed. Comparison with other lagomorphs suggests that this event happened before the Ochotona radiation. The d_N/d_S for the sequence region comprising the PPP1 binding motif and the flanking palindrome highly supports the hypothesis that for Ochotona species this region has been evolving under positive selection. In addition, mutational screening shows that the ability of pikas TCTEX1D4 to bind to PPP1 is maintained, although the PPP1 binding motif is disrupted, and the N- and C-terminal surrounding residues are also abrogated. These observations suggest pika as an ideal model to study novel PPP1 complexes regulatory mechanisms.