Rabbit hemorrhagic disease virus (RHDV) is a highly lethal Lagovirus, family Caliciviridae, that threatens European rabbits (Oryctolagus cuniculus). Although a related virus severely affects hares, cross-species infection was only recently described for new variant RHDV in Cape hares (Lepus capensis mediterraneus). We sequenced two strains from dead Iberian hares (Lepus granatensis) collected in the 1990s in Portugal. Clinical signs were compatible with a Lagovirus infection. Phylogenetic analysis of the complete capsid gene positioned them in the RHDV genogroup that circulated on the Iberian Peninsula at that time. This is the earliest evidence of RHDV affecting a species other than European rabbits.

The innate immunological response in mammals involves a diverse and complex network of many proteins. Over the last years, the tripartite motif-containing protein 5 (TRIM5) and 22 (TRIM22) have shown promise as restriction factors of a plethora of viruses that infect primates. Although there have been studies describing the evolution of these proteins in a wide range of mammals, no prior studies of the TRIM6/34/5/22 gene cluster have been performed in the Chiroptera order. Here, we provide a detailed analysis of the evolution of this gene cluster in several bat genomes. Examination of different yangochiroptera and yinpterochiroptera bat species revealed a dynamic history of gene expansion occurring in TRIM5 and TRIM22 genes. Multiple copies of TRIM5 were found in the genomes of several bats, demonstrating a very low degree of conservation in the synteny of this gene among species of the Chiroptera order. Our findings also reveal that TRIM22 is often found duplicated in yangochiroptera bat species, an evolutionary phenomenon not yet observed in any other lineages of mammals. In total, we identified 31 TRIM5 and 19 TRIM22 amino acids to be evolving under positive selection, with most of the residues being placed in the PRYSPRY domain, known to be responsible for binding to the viral capsid during restriction in the primate orthologous TRIM proteins. Altogether, our results help to shed light on the distinctive role of bats in nature as reservoirs of viruses, many of which have become threatening zoonotic diseases through virus spillover in the last decades.

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.

Background: Toll-like receptors (TLRs) are the most widely studied innate immunity receptors responsible for recognition of invading pathogens. Among the TLR family, TLR5 is the only that senses and recognizes flagellin, the major protein of bacterial flagella. TLR5 has been reported to be under overall purifying selection in mammals, with a small proportion of codons under positive selection. However, the variation of substitution rates among major mammalian groups has been neglected. Here, we studied the evolution of TLR5 in mammals, comparing the substitution rates among groups.Results: In this study we analysed the TLR5 substitution rates in Euungulata, Carnivora, Chiroptera, Primata, Rodentia and Lagomorpha, groups. For that, Tajima’s relative rate test, Bayesian inference of evolutionary rates and genetic distances were estimated with CODEML’s branch model and RELAX. The combined results showed that in the Lagomorpha, Rodentia, Carnivora and Chiroptera lineages TLR5 is evolving at a higher substitution rate. The RELAX analysis further suggested a significant relaxation of selective pressures for the Lagomorpha (K = 0.22,p < 0.01), Rodentia (K = 0.58,p < 0.01) and Chiroptera (K = 0.65,p < 0.01) lineages and for the Carnivora ancestral branches (K = 0.13,p < 0.01).Conclusions: Our results show that the TLR5 substitution rate is not uniform among mammals. In fact, among the different mammal groups studied, the Lagomorpha, Rodentia, Carnivora and Chiroptera are evolving faster. This evolutionary pattern could be explained by 1) the acquisition of new functions of TLR5 in the groups with higher substitution rate, i.e. TLR5 neofunctionalization, 2) by the beginning of a TLR5 pseudogenization in these groups due to some redundancy between the TLRs genes, or 3) an arms race between TLR5 and species-specific parasites.

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.

Rabbit hemorrhagic disease is caused by a calicivirus, rabbit hemorrhagic disease virus (RHDV), which is responsible for high mortality in domestic and wild European rabbits (Oryctolagus cuniculus). RHDV strains were sequenced from wild European rabbits (Oryctolagus cuniculus algirus) collected in the Azorean island of Pico, Portugal. Phylogenetic analyses showed that the Pico RHDV strains diverge from all of the others described so far, but cluster with the genogroups 1–5 (G1–G5). The genetic distance between the Pico RHDV sequences and each G1, G2 and G3–G5 genogroup (~0.08) is compatible with an RHDV introduction at least 17 years ago. Our results show that in Pico, RHDV is the outcome of an independent evolution from the original RHDV strain that appeared in its European rabbit population. These are the first sequences of RHDV obtained in the subspecies O. c. algirus, outside of its original region, the Iberian Peninsula. Furthermore, we discuss the risk of rabbit translocations from the Azores to the Iberian Peninsula, where the rabbit wild populations are suffering high mortalities.