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.

Rabbit IgA Hinges That Resist IgA1 Protease Action Provide Options for Improved IgA-Based Therapeutic Agents.

Immunoglobulin A provides a major line of defence against pathogens and plays a key role in the maintenance of the commensal microbiota in the intestinal tract. Having been shown to be more effective at tumour cell killing than IgG and strongly active against pathogens present in the mucosae, IgA antibodies have been attracting significant attention in recent years for use as therapeutic antibodies. To improve their therapeutic potential, bioengineered IgA forms with increased serum half-life and neutralizing abilities have been developed but the IgA hinge, which impacts susceptibility to bacterial proteases and ability to bridge between target and effector cells, has not yet been explored. The European rabbit has 15 IgA subclasses with exclusive hinge region motifs and varying lengths, constituting a unique model to evaluate the functional capabilities offered by incorporation of longer IgA hinges into immunoglobulins. Hinge regions from rabbit IgAs, featuring different lengths and sequences, were inserted into human IgA1 heavy chain to substitute the IgA1 hinge. These hinges did not appear to affect antigen binding nor the ability of the engineered chimeric IgA1 to bind and trigger FcαRI, as detected by IgA-mediated cell agglutination and release of superoxide by neutrophils. All rabbit hinge-human IgA1 hybrids were resistant to Clostridrum ramosum IgA protease enzyme digestion, as predicted by the lack of the cleavage site in the rabbit hinges. Some IgA1s featuring long rabbit hinges were cleaved by Neisseria meningitidis IgA1 protease cleavage type 1 or 2 enzymes, despite the lack of the predicted cleavage sites. More interestingly, the hybrid featuring the rabbit IgA15 hinge was not affected by any of the IgA proteases. The IgA15 hinge is longer than that found in human IgA1 and is composed by a unique motif with a stretch of nine consecutive Ser residues. These characteristics allow the preservation of a long hinge, with associated ability to bridge distantly spaced antigens and provide higher avidity binding, while remaining resistant to IgA protease degradation. The data suggest that the rabbit Cα15 hinge represents an interesting alternative hinge sequence for therapeutic human IgA antibodies that remains resistant to proteolytic cleavage.

Feline Leukemia Virus-B Envelope Together With its GlycoGag and Human Immunodeficiency Virus-1 Nef Mediate Resistance to Feline SERINC5.

Human SERINC5 (SER5) protein is a recently described restriction factor against human immunodeficiency virus-1 (HIV-1), which is antagonized by HIV-1 Nef protein. Other retroviral accessory proteins such as the glycosylated Gag (glycoGag) from the murine leukemia virus (MLV) can also antagonize SER5. In addition, some viruses escape SER5 restriction by expressing a SER5-insensitive envelope (Env) glycoprotein. Here, we studied the activity of human and feline SER5 on HIV-1 and on the two pathogenic retroviruses in cats, feline immunodeficiency virus (FIV) and feline leukemia virus (FeLV). HIV-1 in absence of Nef is restricted by SER5 from domestic cats and protected by its Nef protein. The sensitivity of feline retroviruses FIV and FeLV to human and feline SER5 is considerably different: FIV is sensitive to feline and human SER5 and lacks an obvious mechanism to counteract SER5 activity, while FeLV is relatively resistant to SER5 inhibition. We speculated that similar to MLV, FeLV-A or FeLV-B express glycoGag proteins and investigated their function against human and feline SER5 in wild type and envelope deficient virus variants. We found that the endogenous FeLV recombinant virus, FeLV-B but not wild type exogenous FeLV-A envelope mediates a strong resistance against human and feline SER5. GlycoGag has an additional but moderate role to enhance viral infectivity in the presence of SER5 that seems to be dependent on the FeLV envelope. These findings may explain, why in vivo FeLV-B has a selective advantage and causes higher FeLV levels in infected cats compared to infections of FeLV-A only.

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.

IgA: Structure, Function, and Developability.

Immunoglobulin A (IgA) plays a key role in defending mucosal surfaces against attack by infectious microorganisms. Such sites present a major site of susceptibility due to their vast surface area and their constant exposure to ingested and inhaled material. The importance of IgA to effective immune defence is signalled by the fact that more IgA is produced than all the other immunoglobulin classes combined. Indeed, IgA is not just the most prevalent antibody class at mucosal sites, but is also present at significant concentrations in serum. The unique structural features of the IgA heavy chain allow IgA to polymerise, resulting in mainly dimeric forms, along with some higher polymers, in secretions. Both serum IgA, which is principally monomeric, and secretory forms of IgA are capable of neutralising and removing pathogens through a range of mechanisms, including triggering the IgA Fc receptor known as FcαRI or CD89 on phagocytes. The effectiveness of these elimination processes is highlighted by the fact that various pathogens have evolved mechanisms to thwart such IgA-mediated clearance. As the structure-function relationships governing the varied capabilities of this immunoglobulin class come into increasingly clear focus, and means to circumvent any inherent limitations are developed, IgA-based monoclonal antibodies are set to emerge as new and potent options in the therapeutic arena.

Author Correction: The wide utility of rabbits as models of human diseases.

This article was originally published under a CC BY-NC-SA License, but has now been made available under a CC BY 4.0 License.

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.

Evolutionary studies on the betaretrovirus RERV-H in the Leporidae family reveal an endogenization in the ancestor of Oryctolagus, Bunolagus and Pentalagus at 9 million years ago.

RERV-H was first identified in human tissues and mistaken for a human exogenous retrovirus. However, the integration sites carried by this virus showed that it was instead a European rabbit (Oryctolagus cuniculus) endogenous retrovirus. The first clones retrieved from European rabbit samples represented defective proviruses, although estimation of proviral copy numbers found in the European rabbit genome ranged from hundreds to thousands. Screening for the presence of RERV-H showed the absence of the virus in two other lagomorphs, pika (Ochotona) and hares (Lepus), which diverged from rabbits about 35 and 12 million years ago, respectively. Using a PCR-based approach, samples of seven different Lagomorph genera were tested for the presence of RERV-H. It was possible to amplify a proviral fragment corresponding to RNaseH from Oryctolagus, Bunolagus and Pentalagus genomic samples. The amplification of proviral DNA in species other than Oryctolagus revealed that this virus was endogenized in their common ancestor, roughly 9 million years ago. Using the European rabbit genome sequence OryCun2.0, it was possible to find multiple copies spread throughout the genome and several complete proviral genomes were retrieved. Some copies contained full open reading frames for all viral components. The lack of a complete genome in the other Lagomorph species did not allow further analyses of the provirus, although more deleterious mutations were found in Bunolagus and Pentalagus than in Oryctolagus RNaseH-amplified sequences. To what extent RERV-H and other endogenous viruses might have had an impact on the rabbit genome and its immune system remains elusive.

The antiviral activity of rodent and lagomorph SERINC3 and SERINC5 is counteracted by known viral antagonists.

A first step towards the development of a human immunodeficiency virus (HIV) animal model has been the identification and surmounting of species-specific barriers encountered by HIV along its replication cycle in cells from small animals. Serine incorporator proteins 3 (SERINC3) and 5 (SERINC5) were recently identified as restriction factors that reduce HIV-1 infectivity. Here, we compared the antiviral activity of SERINC3 and SERINC5 among mice, rats and rabbits, and their susceptibility to viral counteraction to their human counterparts. In the absence of viral antagonists, rodent and lagomorph SERINC3 and SERINC5 displayed anti-HIV activity in a similar range to human controls. Vesicular stomatitis virus G protein (VSV-G) pseudotyped virions were considerably less sensitive to restriction by all SERINC3/5 orthologs. Interestingly, HIV-1 Nef, murine leukemia virus (MLV) GlycoGag and equine infectious anemia virus (EIAV) S2 counteracted the antiviral activity of all SERINC3/5 orthologs with similar efficiency. Our results demonstrate that the antiviral activity of SERINC3/5 proteins is conserved in rodents and rabbits, and can be overcome by all three previously reported viral antagonists.

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.

The wide utility of rabbits as models of human diseases.

Studies using the European rabbit Oryctolagus cuniculus contributed to elucidating numerous fundamental aspects of antibody structure and diversification mechanisms and continue to be valuable for the development and testing of therapeutic humanized polyclonal and monoclonal antibodies. Additionally, during the last two decades, the use of the European rabbit as an animal model has been increasingly extended to many human diseases. This review documents the continuing wide utility of the rabbit as a reliable disease model for development of therapeutics and vaccines and studies of the cellular and molecular mechanisms underlying many human diseases. Examples include syphilis, tuberculosis, HIV-AIDS, acute hepatic failure and diseases caused by noroviruses, ocular herpes, and papillomaviruses. The use of rabbits for vaccine development studies, which began with Louis Pasteur's rabies vaccine in 1881, continues today with targets that include the potentially blinding HSV-1 virus infection and HIV-AIDS. Additionally, two highly fatal viral diseases, rabbit hemorrhagic disease and myxomatosis, affect the European rabbit and provide unique models to understand co-evolution between a vertebrate host and viral pathogens.

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.

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.

Cross-species comparison of mammalian saliva using an LC-MALDI based proteomic approach.

Despite the importance of saliva in the regulation of oral cavity homeostasis, few studies have been conducted to quantitatively compare the saliva of different mammal species. Aiming to define a proteome signature of mammals' saliva, an in-depth SDS-PAGE-LC coupled to MS/MS (GeLC-MS/MS) approach was used to characterize the saliva from primates (human), carnivores (dog), glires (rat and rabbit), and ungulates (sheep, cattle, horse). Despite the high variability in the number of distinct proteins identified per species, most protein families were shared by the mammals studied with the exception of cattle and horse. Alpha-amylase is an example that seems to reflect the natural selection related to digestion efficacy and food recognition. Casein protein family was identified in all species but human, suggesting an alternative to statherin in the protection of hard tissues. Overall, data suggest that different proteins might assure a similar role in the regulation of oral cavity homeostasis, potentially explaining the specific mammals' salivary proteome signature. Moreover, some protein families were identified for the first time in the saliva of some species, the presence of proline-rich proteins in rabbit's saliva being a good example.

Evolution of C, D and S-type cystatins in mammals: an extensive gene duplication in primates.

Cystatins are a family of inhibitors of cysteine peptidases that comprises the salivary cystatins (D and S-type cystatins) and cystatin C. These cystatins are encoded by a multigene family (CST3, CST5, CST4, CST1 and CST2) organized in tandem in the human genome. Their presence and functional importance in human saliva has been reported, however the distribution of these proteins in other mammals is still unclear. Here, we performed a proteomic analysis of the saliva of several mammals and studied the evolution of this multigene family. The proteomic analysis detected S-type cystatins (S, SA, and SN) in human saliva and cystatin D in rat saliva. The evolutionary analysis showed that the cystatin C encoding gene is present in species of the most representative mammalian groups, i.e. Artiodactyla, Rodentia, Lagomorpha, Carnivora and Primates. On the other hand, D and S-type cystatins are mainly retrieved from Primates, and especially the evolution of S-type cystatins seems to be a dynamic process as seen in Pongo abelii genome where several copies of CST1-like gene (cystatin SN) were found. In Rodents, a group of cystatins previously identified as D and S has also evolved. Despite the high divergence of the amino acid sequence, their position in the phylogenetic tree and their genome organization suggests a common origin with those of the Primates. These results suggest that the D and S type cystatins have emerged before the mammalian radiation and were retained only in Primates and Rodents. Although the mechanisms driving the evolution of cystatins are unknown, it seems to be a dynamic process with several gene duplications evolving according to the birth-and-death model of evolution. The factors that led to the appearance of a group of saliva-specific cystatins in Primates and its rapid evolution remain undetermined, but may be associated with an adaptive advantage.

An evolutionary perspective of mammal salivary peptide families: cystatins, histatins, statherin and PRPs.

Saliva's role in the oral cavity, such as lubrication, protection of tissues and antimicrobial action is a reflex of its composition, among which are several peptide families like statherin, histatins, proline rich proteins (PRPs) and some members of the cystatin family. These peptides present different evolutionary pathways being in the case of histatin, statherin and PRP families restricted to few millions and comprising few species when compared with cystatins, where duplication occurred at more than 650 mya. Though the recognized relevance of phylogenetic approaches to disclose relationships among different species, information on the salivary proteins that allow the association between peptide families-related structure and function in the oral cavity is scarce. In the present study, the four major salivary peptides classes are reviewed considering the few known phylogenetic studies focusing on their evolution among mammals. New perspectives and challenges for future and multidisciplinary experimental works are drawn.

Characterization of thymosin ß4 in mammals' saliva.

Thymosin ß4 (Tß4) is a low molecular weight peptide found in several mammalian tissues and is known mainly by its ability to bind cytoskeletal actin, influencing cell migration and differentiation, and promoting tissue repair. Considering the functional role of this peptide, the main goal of this work was to characterize Tß4 in mammals' saliva by using evolutionary and proteomic tools. For this, mammalian Tß4 sequences were retrieved from NCBI, SwissProt and Ensembl databases. The alignment of Tß4 amino acid sequences showed a high degree of conservation between species. The gene seems to be evolving under negative selection as indicated by a dN/dS ratio of 0.05. Whole saliva was collected from dog, human, rabbit, cow, horse and sheep and the salivary peptides were isolated through filtration and analyzed by LC-MS/MS. Spectra was processed against the database constructed with the retrieved Tß4 sequences. For the first time, the identification of this peptide was achieved in rat, dog, horse and bovine saliva. Detection in these mammal species and its amino acid conservation suggest an important role of Tß4 in the homeostasis of the mammalian oral cavity.