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

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.

Evolution of Fc Receptor-Like Scavenger in Mammals

Fc receptor-like (FCRL) molecules comprise a large family of receptors, homologous to the receptors for the Fc portion of immunoglobulins (FCR). Within this family, an unusual gene known to exist in mice, rats and dogs, termed FCRLS, encodes a chimeric protein with both Ig-like FCRL and type B scavenger-receptor cysteine-rich (SRCR)-like domains. In mice, FCRLS is located next to the CD5L and KIRREL1 genes. Here, we show that the curious FCRLS gene is actually present across major mammalian groups, but its annotation is generally incorrect or absent. Anchored on mouse FCRLS and FCRL2 genomic sequence alignments, phylogenetic analyses demonstrated that many mammalian sequences currently annotated as FCRL2 cluster with FCRLS, supported by a conserved genetic synteny among organisms. This analysis shows that FCRLS is present in Rodentia, some Carnivora (Canidae and Ursidae), Chiroptera, Arctiodactyla, Proboscidae, and some Primata. Thus, the FCRLS most likely originated in a eutherian mammal ancestor since it is not present in Monotremata or Marsupialia. FCRLS has a peculiar distribution pattern across mammalian lineages, being present in some species, but absent in others from the same family, as in carnivores for example. The most parsimonious hypothesis to explain this FCRLS evolution is that it was convergently lost in several independent mammalian lineages. Analyses of branch-specific nucleotide evolutionary rates, show that FCRL2 and FCRLS have similar ranges of rates across mammals, suggesting that both genes have crucial, but separate functions in the immune system. Bayesian estimates of evolutionary rates for FCRLS in mammalian lineages revealed that carnivores display the highest mutation rate after rodents. Additionally, positive diversifying selection was detected for both FCRL2 and FCRLS. Our results show that the presence of the FCRLS gene is older and more widespread across mammals than previously thought and appears to be functional, being under positive selection. Its precise physiologic role should thus be investigated.

Not so unique to Primates: The independent adaptive evolution of TRIM5 in Lagomorpha lineage

The plethora of restriction factors with the ability to inhibit the replication of retroviruses have been widely studied and genetic hallmarks of evolutionary selective pressures in Primates have been well documented. One example is the tripartite motif-containing protein 5 alpha (TRIM5α), a cytoplasmic factor that restricts retroviral infection in a species-specific fashion. In Lagomorphs, similarly to what has been observed in Primates, the specificity of TRIM5 restriction has been assigned to the PRYSPRY domain. In this study, we present the first insight of an intra-genus variability within the Lagomorpha TRIM5 PRYSPRY domain. Remarkably, and considering just the 32 residue-long v1 region of this domain, the deduced amino acid sequences of Daurian pika (Ochotona dauurica) and steppe pika (O. pusilla) evidenced a high divergence when compared to the remaining Ochotona species, presenting values of 44% and 66% of amino acid differences, respectively. The same evolutionary pattern was also observed when comparing the v1 region of two Sylvilagus species members (47% divergence). However, and unexpectedly, the PRYSPRY domain of Lepus species exhibited a great conservation. Our results show a high level of variation in the PRYSPRY domain of Lagomorpha species that belong to the same genus. This suggests that, throughout evolution, the Lagomorpha TRIM5 should have been influenced by constant selective pressures, likely as a result of multiple different retroviral infections.

Genetic Diversity of IGHM and IGHE in the Leporids Revealed Different Patterns of Diversity in the Two European Rabbit Subspecies (O. cuniculus algirus and O. c. cuniculus)

Simple Summary

The study of European rabbit immunoglobulin genes has contributed decisively to the current knowledge on antibody structure and diversification. The European rabbit has also been increasingly used as an animal model for the study of many human diseases, such as syphilis, tuberculosis, and AIDS. As such, the study of its immune system genes is of crucial relevance, but the study of rabbit immunoglobulins has focused only on the IgG and IgA antibodies. In this study, we added to the knowledge of the rabbit immune system by investigating the genetic diversity of two antibodies, IgM and IgE, in wild and domestic rabbits as well as other rabbit close species. With the data obtained in this study, we showed a high similarity between the different rabbit close species studied and we pointed out important genetic differences in the wild and domestic rabbits. Our findings are a valuable tool for the management of rabbit wild populations and domestic breeds and may contribute to the identification of immunoglobulins genetic variants with greater efficiency against pathogens.

Abstract

The European rabbit (Oryctolagus cuniculus) has been an important model for immunological studies but the study of its immunoglobulins (Ig) has been restricted to its unique IgA and IgG. Here, we studied the genetic diversity of IgM and IgE in several species of leporids and performed population genetics studies on European rabbit wild populations and domestic breeds. The leporids sequencing showed that these Ig are well conserved (98% sequence similarity among leporids), For IgM the Cµ1 and Cµ4 were the most diverse and most conserved domains, respectively, while for IgE the Cε1 was the most diverse domain and Cε2 and Cε3 the most conserved domains. The differences in the pattern of most conserved and most diverse domain between the Ig isotypes are most likely related to each isotype function. The genetic population data showed contrasting results for IgM and IgE. For both Ig, as expected, a greater diversity was observed in the original species range, the Iberian Peninsula. However, unexpectedly the genetic diversity found for IgE in the domestic animals is higher than that for the French wild populations. These results will increase knowledge of the genetic diversity of leporids and wild and domestic rabbit populations and are important tools for the management of wild populations and rabbitries.

Overcoming species barriers: an outbreak of Lagovirus europaeus GI.2/RHDV2 in an isolated population of mountain hares (Lepus timidus)

Background

Prior to 2010, the lagoviruses that cause rabbit hemorrhagic disease (RHD) in European rabbits (Oryctolagus cuniculus) and European brown hare syndrome (EBHS) in hares (Lepus spp.) were generally genus-specific. However, in 2010, rabbit hemorrhagic disease virus 2 (RHDV2), also known as Lagovirus europaeus GI.2, emerged and had the distinguishing ability to cause disease in both rabbits and certain hare species. The mountain hare (Lepus timidus) is native to Sweden and is susceptible to European brown hare syndrome virus (EBHSV), also called Lagovirus europaeus GII.1. While most mountain hare populations are found on the mainland, isolated populations also exist on islands. Here we investigate a mortality event in mountain hares on the small island of Hallands Väderö where other leporid species, including rabbits, are absent.

Results

Post-mortem and microscopic examination of three mountain hare carcasses collected from early November 2016 to mid-March 2017 revealed acute hepatic necrosis consistent with pathogenic lagovirus infection. Using immunohistochemistry, lagoviral capsid antigen was visualized within lesions, both in hepatocytes and macrophages. Genotyping and immunotyping of the virus independently confirmed infection with L. europaeus GI.2, not GII.1. Phylogenetic analyses of the vp60 gene grouped mountain hare strains together with a rabbit strain from an outbreak of GI.2 in July 2016, collected approximately 50 km away on the mainland.

Conclusions

This is the first documented infection of GI.2 in mountain hares and further expands the host range of GI.2. Lesions and tissue distribution mimic those of GII.1 in mountain hares. The virus was most likely initially introduced from a concurrent, large-scale GI.2 outbreak in rabbits on the adjacent mainland, providing another example of how readily this virus can spread. The mortality event in mountain hares lasted for at least 4.5 months in the absence of rabbits, which would have required virus circulation among mountain hares, environmental persistence and/or multiple introductions. This marks the fourth Lepus species that can succumb to GI.2 infection, suggesting that susceptibility to GI.2 may be common in Lepus species. Measures to minimize the spread of GI.2 to vulnerable Lepus populations therefore are prudent.

Electronic supplementary material

The online version of this article (10.1186/s12917-018-1694-7) contains supplementary material, which is available to authorized users.

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.

Mechanistic insights into molecular evolution of species-specific differential glycosaminoglycan binding surfaces in growth-related oncogene chemokines

Chemokines are chemotactic cytokines involved in leucocyte trafficking to infected tissue. Growth-related oncogene (GRO) chemokines namely CXCL1, CXCL2 and CXCL3 are neutrophil activating chemokines sharing a conserved three-dimensional structure, but encompassing functional diversity due to gene duplication and evolutionary events. However, the evolutionary mechanisms including selection pressures involved in diversification of GRO genes have not yet been characterized. Here, we performed comprehensive evolutionary analysis of GRO genes among different mammalian species. Phylogenetic analysis illustrated a species-specific evolution pattern. Selection analysis evidenced that these genes have undergone concerted evolution. Seventeen positively selected sites were obtained, although the majority of the protein is under purifying selection. Interestingly, these positively selected sites are more concentrated on the C-terminal/glycosaminoglycan (GAG) binding and dimerization segment compared to receptor binding domain. Substitution rate analysis confirmed the C-terminal domain of GRO genes as the highest substituted segment. Further, structural analysis established that the nucleotide alterations in the GAG binding domain are the source of surface charge modulation, thus generating the differential GAG binding surfaces and multiple binding sites as per evolutionary pressure, although the helical surface is primordial for GAG binding. Indeed, such variable electrostatic surfaces are crucial to regulate chemokine gradient formation during a host's defence against pathogens and also explain the significance of chemokine promiscuity.

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.

Evolution of CCL11: genetic characterization in lagomorphs and evidence of positive and purifying selection in mammals

The interactions between chemokines and their receptors are crucial for differentiation and activation of inflammatory cells. CC chemokine ligand 11 (CCL11) binds to CCR3 and to CCR5 that in leporids underwent gene conversion with CCR2. Here, we genetically characterized CCL11 in lagomorphs (leporids and pikas). All lagomorphs have a potentially functional CCL11, and the Pygmy rabbit has a mutation in the stop codon that leads to a longer protein. Other mammals also have mutations at the stop codon that result in proteins with different lengths. By employing maximum likelihood methods, we observed that, in mammals, CCL11 exhibits both signatures of purifying and positive selection. Signatures of purifying selection were detected in sites important for receptor binding and activation. Of the three sites detected as under positive selection, two were located close to the stop codon. Our results suggest that CCL11 is functional in all lagomorphs, and that the signatures of purifying and positive selection in mammalian CCL11 probably reflect the protein's biological roles.

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 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.

Pseudogenization of CCL14 in the Ochotonidae (pika) family

The interaction between chemokines and their receptors is crucial for inflammatory cell trafficking. CCL14 binds with high affinity to CCR5. In leporids, CCR5 underwent gene conversion with CCR2. The study of CCR5 ligands in leporid species showed that CCL8 is pseudogenized, while CCL3, CCL4 and CCL5 are functional. Here, we study the evolution of CCL14 in mammals with emphasis in the order Lagomorpha. By employing maximum likelihood methods we detected six sites under positive selection. Some of these sites are located in regions crucial for CCL14 activation and binding to receptors. Sequencing of CCL14 in Ochotona species showed that O. princeps, O. pallasi, O. alpina and O. turuchanensis have a mutation at the start codon (Met > Thr), while O. hoffmanni, O. mantchurica, O. dauurica and O. rufescens present the mammalian conserved Met. Ochotona hyperborea has the two alleles. In O. pusilla, CCL14 is a pseudogene due to a seven base pair insertion. Like CCL3, CCL4 and CCL5, CCL14 is functional in all leporids but in the Ochotonidae family it underwent a pseudogenization process. This suggests that CCL14 has an important biological role in other mammals by evolving under positive selection that has been lost in Ochotonidae (subgenera Pika and Lagotona).