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.

Alternated selection mechanisms maintain adaptive diversity in different demographic scenarios of a large carnivore.

BACKGROUND: Different population trajectories are expected to impact the signature of neutral and adaptive processes at multiple levels, challenging the assessment of the relative roles of different microevolutionary forces. Here, we integrate adaptive and neutral variability patterns to disentangle how adaptive diversity is driven under different demographic scenarios within the Iberian wolf (Canis lupus) range. We studied the persistent, the expanding and a small, isolated group within the Iberian wolf population, using 3 MHC class II genes (DRB1, DQA1, and DQB1), which diversity was compared with 39 microsatellite loci. RESULTS: Both the persistent and the expanding groups show evidence of balancing selection, revealed by a significant departure from neutrality at MHC loci, significant higher observed and expected heterozygosity and lower differentiation at MHC than at neutral loci, and signs of positive selection. However, despite exhibiting a significantly higher genetic diversity than the isolated group, the persistent group did not show significant excess of MHC heterozygotes. The expanding group, while showing a similar level of genetic diversity than the persistent group, displays by contrast a significant excess of MHC heterozygotes, which is compatible with the heterozygote advantage mechanism. Results are not clear regarding the role of drift and selection in the isolated group due to the small size of this population. Although diversity indices of MHC loci correspond to neutral expectations in the isolated group, accelerated MHC divergence, revealed by a higher differentiation at MHC than neutral loci, may indicate diversifying selection. CONCLUSION: Different selective pressures were observed in the three different demographic scenarios, which are possibly driven by different selection mechanisms to maintain adaptive diversity.

Evolution of CCL16 in Glires (Rodentia and Lagomorpha) shows an unusual random pseudogenization pattern.

BACKGROUND: The C-C motif chemokine ligand 16 (CCL16) is a potent pro-inflammatory chemokine and a chemoattractant for monocytes and lymphocytes. In normal plasma, it is present at high concentrations and elicits its effects on cells by interacting with cell surface chemokine receptors. In the European rabbit and in rodents such as mouse, rat and guinea pig, CCL16 was identified as a pseudogene, while in the thirteen-lined ground squirrel it appears to be potentially functional. To gain insight into the evolution of this gene in the superorder Glires (rodents and lagomorphs), we amplified the CCL16 gene from eleven Leporidae and seven Ochotonidae species. RESULTS: We compared our sequences with CCL16 sequences of twelve rodent species retrieved from public databases. The data show that for all leporid species studied CCL16 is a pseudogene. This is primarily due to mutations at the canonical Cys Cys motif, creating either premature stop codons, or disrupting amino acid replacements. In the Mexican cottontail, CCL16 is pseudogenized due to a frameshift deletion. Additionally, in the exon 1 (signal peptide), there are frameshift deletions present in all leporids studied. In contrast, in Ochotona species, CCL16 is potentially functional, except for an allele in Hoffmann's pika. In rodents, CCL16 is functional in a number of species, but patterns of pseudogenization similar to those observed in lagomorphs also exist. CONCLUSIONS: Our results suggest that while functional in the Glires ancestor, CCL16 underwent pseudogenization in some species. This process occurred stochastically or in specific lineages at different moments in the evolution of Glires. These observations suggest that the CCL16 had different evolutionary constrains in the Glires group that could be associated with the CCL16 biological function.

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.

Adaptive Gene Loss? Tracing Back the Pseudogenization of the Rabbit CCL8 Chemokine.

Studies of the process of pseudogenization have widened our understanding of adaptive evolutionary change. In Rabbit, an alteration at the second extra-cellular loop of the CCR5 chemokine receptor was found to be associated with the pseudogenization of one of its prime ligands, the chemokine CCL8. This relationship has raised questions about the existence of a causal link between both events, which would imply adaptive gene loss. This hypothesis is evaluated here by tracing back the history of the genetic modifications underlying the chemokine pseudogenization. The obtained data indicate that mutations at receptor and ligand genes occurred after the lineage split of New World Leporids versus Old World Leporids and prior to the generic split of the of Old World species studied, which occurred an estimated 8-9 million years ago. More important, they revealed the emergence, before this zoographical split, of a "slippery" nucleotide motif (CCCCGGG) at the 3' region of CCL8-exon2. Such motives are liable of generating +1G or -1G frameshifts, which could, however, be overcome by "translesion" synthesis or somatic reversion. The CCL8 pseudogenization in the Old World lineage was apparently initiated by three synapomorphic point mutations at the exon2-intron2 boundary which provide at short range premature terminating codons, independently of the reading frame imposed by the slippery motif. The presence of this motif in New World Leporids might allow verifying this scenario. The importance of CCL8-CCR5 signaling in parasite-host interaction would suggest that the CCL8 knock-out in Old World populations might be related to changes in pathogenic environment.

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.

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.

Pseudogenization of the MCP-2/CCL8 chemokine gene in European rabbit (genus Oryctolagus), but not in species of Cottontail rabbit (Sylvilagus) and Hare (Lepus).

BACKGROUND: Recent studies in human have highlighted the importance of the monocyte chemotactic proteins (MCP) in leukocyte trafficking and their effects in inflammatory processes, tumor progression, and HIV-1 infection. In European rabbit (Oryctolagus cuniculus) one of the prime MCP targets, the chemokine receptor CCR5 underwent a unique structural alteration. Until now, no homologue of MCP-2/CCL8a, MCP-3/CCL7 or MCP-4/CCL13 genes have been reported for this species. This is interesting, because at least the first two genes are expressed in most, if not all, mammals studied, and appear to be implicated in a variety of important chemokine ligand-receptor interactions. By assessing the Rabbit Whole Genome Sequence (WGS) data we have searched for orthologs of the mammalian genes of the MCP-Eotaxin cluster. RESULTS: We have localized the orthologs of these chemokine genes in the genome of European rabbit and compared them to those of leporid genera which do (i.e. Oryctolagus and Bunolagus) or do not share the CCR5 alteration with European rabbit (i.e. Lepus and Sylvilagus). Of the Rabbit orthologs of the CCL8, CCL7, and CCL13 genes only the last two were potentially functional, although showing some structural anomalies at the protein level. The ortholog of MCP-2/CCL8 appeared to be pseudogenized by deleterious nucleotide substitutions affecting exon1 and exon2. By analyzing both genomic and cDNA products, these studies were extended to wild specimens of four genera of the Leporidae family: Oryctolagus, Bunolagus, Lepus, and Sylvilagus. It appeared that the anomalies of the MCP-3/CCL7 and MCP-4/CCL13 proteins are shared among the different species of leporids. In contrast, whereas MCP-2/CCL8 was pseudogenized in every studied specimen of the Oryctolagus - Bunolagus lineage, this gene was intact in species of the Lepus - Sylvilagus lineage, and was, at least in Lepus, correctly transcribed. CONCLUSION: The biological function of a gene was often revealed in situations of dysfunction or gene loss. Infections with Myxoma virus (MYXV) tend to be fatal in European rabbit (genus Oryctolagus), while being harmless in Hares (genus Lepus) and benign in Cottontail rabbit (genus Sylvilagus), the natural hosts of the virus. This communication should stimulate research on a possible role of MCP-2/CCL8 in poxvirus related pathogenicity.

Rabbit haemorrhagic disease (RHD) and rabbit haemorrhagic disease virus (RHDV): a review.

Rabbit haemorrhagic disease virus (RHDV) is a calicivirus of the genus Lagovirus that causes rabbit haemorrhagic disease (RHD) in adult European rabbits (Oryctolagus cuniculus). First described in China in 1984, the virus rapidly spread worldwide and is nowadays considered as endemic in several countries. In Australia and New Zealand where rabbits are pests, RHDV was purposely introduced for rabbit biocontrol. Factors that may have precipitated RHD emergence remain unclear, but non-pathogenic strains seem to pre-date the appearance of the pathogenic strains suggesting a key role for the comprehension of the virus origins. All pathogenic strains are classified within one single serotype, but two subtypes are recognised, RHDV and RHDVa. RHD causes high mortality in both domestic and wild adult animals, with individuals succumbing between 48-72 h post-infection. No other species has been reported to be fatally susceptible to RHD. The disease is characterised by acute necrotising hepatitis, but haemorrhages may also be found in other organs, in particular the lungs, heart, and kidneys due to disseminated intravascular coagulation. Resistance to the disease might be explained in part by genetically determined absence or weak expression of attachment factors, but humoral immunity is also important. Disease control in rabbitries relies mainly on vaccination and biosecurity measures. Such measures are difficult to be implemented in wild populations. More recent research has indicated that RHDV might be used as a molecular tool for therapeutic applications. Although the study of RHDV and RHD has been hampered by the lack of an appropriate cell culture system for the virus, several aspects of the replication, epizootology, epidemiology and evolution have been disclosed. This review provides a broad coverage and description of the current knowledge on the disease and the virus.

Study of Sylvilagus rabbit TRIM5α species-specific domain: how ancient endoviruses could have shaped the antiviral repertoire in Lagomorpha.

BACKGROUND: Since the first report of the antiretroviral restriction factor TRIM5α in primates, several orthologs in other mammals have been described. Recent studies suggest that leporid retroviruses like RELIK, the first reported endogenous lentivirus ever, may have imposed positive selection in TRIM5α orthologs of the European rabbit and European brown hare. Considering that RELIK must already have been present in a common ancestor of the leporid genera Lepus, Sylvilagus and Oryctolagus, we extended the study of evolutionary patterns of TRIM5α to other members of the Leporidae family, particularly to the genus Sylvilagus. Therefore, we obtained the TRIM5α nucleotide sequences of additional subspecies and species of the three leporid genera. We also compared lagomorph TRIM5α deduced protein sequences and established TRIM5α gene and TRIM5α protein phylogenies. RESULTS: The deduced protein sequence of Iberian hare TRIM5α was 89% identical to European rabbit TRIM5α, although high divergence was observed at the PRYSPRY v1 region between rabbit and the identified alleles from this hare species (allele 1: 50% divergence; allele 2: 53% divergence). A high identity was expected between the Sylvilagus and Oryctolagus TRIM5α proteins and, in fact, the Sylvilagus TRIM5α was 91% identical to the Oryctolagus protein. Nevertheless, the PRYSPRY v1 region was only 50% similar between these genera. Selection analysis of Lagomorpha TRIM5α proteins identified 25 positively-selected codons, 11 of which are located in the PRYSPRY v1 region, responsible for species specific differences in viral capsid recognition. CONCLUSIONS: By extending Lagomorpha TRIM5α studies to an additional genus known to bear RELIK, we verified that the divergent species-specific pattern observed between the Oryctolagus and Lepus PRYSPRY-domains is also present in Sylvilagus TRIM5α. This work is one of the first known studies that compare the evolution of the antiretroviral restriction factor TRIM5α in different mammalian groups, Lagomorpha and Primates.

Molecular bases of genetic diversity and evolution of the immunoglobulin heavy chain variable region (IGHV) gene locus in leporids.

The rabbit has long been a model for studies of the immune system. Work using rabbits contributed both to the battle against infectious diseases such as rabies and syphilis, and to our knowledge, of antibodies' structure, function, and regulated expression. With the description of rabbit Ig allotypes, the discovery of different gene segments encoding immunoglobulins became possible. This challenged the "one gene-one protein" dogma. The observation that rabbit allotypic specificities of the variable regions were present on IgM and IgG molecules also led to the hypothesis of Ig class switching. Rabbit allotypes contributed to the documentation of phenomena such as allelic exclusion and imbalance in production of allelic gene products. During the last 30 years, the rabbit Ig allotypes revealed a number of unique features, setting them apart from mice, humans, and other mammals. Here, we review the most relevant findings concerning the rabbit IGHV. Among these are the preferential usage of one VH gene in VDJ rearrangements, the existence of trans-species polymorphism in the IGHV locus revealed by serology and confirmed by sequencing IGHV genes in Lepus, the unusually large genetic distances between allelic lineages and the fact that the antibody repertoire is diversified in this species only after birth. The whole genome sequence of a rabbit, plus re-sequencing of additional strains and related genera, will allow further evolutionary investigations of antibody variation. Continued research will help define the roles that genetic, allelic, and population diversity at antibody loci may play in host-parasite interactions.

Diversity and evolutionary history of the MHC DQA gene in leporids.

The European rabbit (Oryctolagus cuniculus) is used as a model for many human diseases, yet comparatively little is known of its genetics, particularly at important loci such as the major histocompatibility complex (MHC). This study investigated genetic diversity and evolutionary history of the DQA gene in a range of leporid species by analysing coding sequence diversity of exon 2 and intron 2 in 53 individuals of 16 different species. Fifty leporid DQA alleles were detected, including 13 novel European rabbit alleles. In the rabbit, the highest levels of diversity were observed in wild rabbits from Portugal, with wild rabbits from England and domestic rabbits showing less diversity. Within the sample, several recombination events were detected and trans-specific evolution of alleles was evidenced, both being general characteristics of mammalian MHC genes. Positive selection is implicated as operating on six codons within exon 2, which are also subject to positive selection in other mammals. Some of these positions are putative antigen recognition sites and underline the importance of pathogen-driven selection on these MHC genes.

Allelic variation at the VHa locus in natural populations of rabbit (Oryctolagus cuniculus, L.).

The large interallelic distances between the three rabbit Ig V(H)a lineages, a1, a2 and a3, suggest that the persistence time of the V(H)a polymorphism could amount to 50 million years, which is much longer than that of MHC polymorphisms. Rabbit originated in the Iberian Peninsula where two subspecies coexist, one of which is confined to Southwestern Iberia (Oryctolagus cuniculus algirus). We studied the V(H) loci in the original species range to obtain a better understanding of the evolutionary history of this unusual polymorphism. Serological surveys revealed that sera from the subspecies algirus, when tested with V(H)a locus-specific alloantisera, showed either cross-reactivity ("a-positive" variants) or no reaction at all ("a-blank"). Using RT-PCR, we determined 120 sequences of rearranged V(H) genes expressed in seven algirus rabbits that were typed as either a-positive or a-blank. The data show that the V(H) genes transcribed in a-positive rabbits are closely related to the V(H)1 alleles of domestic rabbits. In contrast, a-blank rabbits were found to preferentially use V(H) genes that, although clearly related to the known V(H)a genes, define a new major allotypic lineage, designated a4. The a4 sequences have hallmark rabbit V(H)a residues together with a number of unprecedented amino acid changes in framework region 2 and 3. The net protein distances between the V(H)a4 and the V(H)a1, a2, and a3 lineages were 20, 29, and 21% respectively. We conclude that at least four distantly related lineages of the rabbit V(H)a locus exist, one of which seems to be endemic in the Iberian range.

Heavy-chain antibodies in Camelidae; a case of evolutionary innovation.

The emergence in Camelidae species of functional antibodies devoid of light chains (referred to as heavy-chain antibodies or HCAbs) is an intriguing evolutionary event. Homodimeric HCAbs have also been documented in spotted ratfish ( Cos5-Abs) and nurse shark (NAR). To reveal the evolutionary history of HCAbs, we evaluated the phylogenetic and phenotypic relationships among HCAbs and conventional antibodies across taxa and confirmed the current viewpoint that different groups of HCAbs have evolved independently in the three lineages. At least, in the camelids, HCAbs are not the result of resuscitation of dormant genes. They are derived from the conventional antibodies within the Camelidae lineage, and are apparently the outcome of more recent adaptive changes occurring in the compartment of heteromeric antibodies. The shared structural properties of HCAbs across taxa are therefore explained by convergent evolution due to similar constraints related to the absence of pairing to the light chain. It appears that innovative evolutionary changes in Camelidae have led to a new level of antigen binding repertoire diversification and have allowed acquisition of novel antigen-receptor properties.