New Insights into Rabbit Viral Diseases.

Viruses are responsible for many devastating rabbit diseases that impact their health and welfare and put their conservation and economic revenue at risk [...].

Comment on Shah et al. Genetic Characteristics and Phylogeographic Dynamics of Lagoviruses, 1988-2021. Viruses 2023, 15, 815.

Shah and colleagues [...].

Evolution of guanylate binding protein genes shows a remarkable variability within bats (Chiroptera).

Background: GBPs (guanylate binding proteins), an evolutionary ancient protein family, play a key role in the host's innate immune response against bacterial, parasitic and viral infections. In Humans, seven GBP genes have been described (GBP1-7). Despite the interest these proteins have received over the last years, evolutionary studies have only been performed in primates, Tupaia and rodents. These have shown a pattern of gene gain and loss in each family, indicative of the birth-and-death evolution process. Results: In this study, we analysed the evolution of this gene cluster in several bat species, belonging to the Yangochiroptera and Yinpterochiroptera sub-orders. Detailed analysis shows a conserved synteny and a gene expansion and loss history. Phylogenetic analysis showed that bats have GBPs 1,2 and 4-6. GBP2 has been lost in several bat families, being present only in Hipposideidae and Pteropodidae. GBPs1, 4 and 5 are present mostly as single-copy genes in all families but have suffered duplication events, particularly in Myotis myotis and Eptesicus fuscus. Most interestingly, we demonstrate that GBP6 duplicated in a Chiroptera ancestor species originating two genes, which we named GBP6a and GBP6b, with different subsequent evolutionary histories. GBP6a underwent several duplication events in all families while GBP6b is present as a single copy gene and has been lost in Pteropodidae, Miniopteridae and Desmodus rotundus, a Phyllostomidae. With 14 and 15 GBP genes, Myotis myotis and Eptesicus fuscus stand out as having far more copies than all other studied bat species. Antagonistically, Pteropodidae have the lowest number of GBP genes in bats. Conclusion: Bats are important reservoirs of viruses, many of which have become zoonotic diseases in the last decades. Further functional studies on bats GBPs will help elucidate their function, evolutionary history, and the role of bats as virus reservoirs.

Patterns of Evolution of TRIM Genes Highlight the Evolutionary Plasticity of Antiviral Effectors in Mammals.

The innate immune system of mammals is formed by a complex web of interacting proteins, which together constitute the first barrier of entry for infectious pathogens. Genes from the E3-ubiquitin ligase tripartite motif (TRIM) family have been shown to play an important role in the innate immune system by restricting the activity of different retrovirus species. For example, TRIM5 and TRIM22 have both been associated with HIV restriction and are regarded as crucial parts of the antiretroviral machinery of mammals. Our analyses of positive selection corroborate the great significance of these genes for some groups of mammals. However, we also show that many species lack TRIM5 and TRIM22 altogether. By analyzing a large number of mammalian genomes, here we provide the first comprehensive view of the evolution of these genes in eutherians, showcasing that the pattern of accumulation of TRIM genes has been dissimilar across mammalian orders. Our data suggest that these differences are caused by the evolutionary plasticity of the immune system of eutherians, which have adapted to use different strategies to combat retrovirus infections. Altogether, our results provide insights into the dissimilar evolution of a representative family of restriction factors, highlighting an example of adaptive and idiosyncratic evolution in the innate immune system.

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.

Corrigendum: Uncovering a 500 million year old history and evidence of pseudogenization for TLR15.

[This corrects the article DOI: 10.3389/fimmu.2022.1020601.].

The potential role of scavenging flies as mechanical vectors of Lagovirus europaeus/GI.2.

The European rabbit (Oryctolagus cuniculus) populations of the Iberian Peninsula have been severely affected by the emergence of the rabbit haemorrhagic disease virus (RHDV) Lagovirus europaeus/GI.2 (RHDV2/b). Bushflies and blowflies (Muscidae and Calliphoridae families, respectively) are important RHDV vectors in Oceania, but their epidemiological role is unknown in the native range of the European rabbit. In this study, scavenging flies were collected between June 2018 and February 2019 in baited traps at one site in southern Portugal, alongside a longitudinal capture-mark-recapture study of a wild European rabbit population, aiming to provide evidence of mechanical transmission of GI.2 by flies. Fly abundance, particularly from Calliphoridae and Muscidae families, peaked in October 2018 and in February 2019. By employing molecular tools, we were able to detect the presence of GI.2 in flies belonging to the families Calliphoridae, Muscidae, Fanniidae and Drosophilidae. The positive samples were detected during an RHD outbreak and absent in samples collected when no evidence of viral circulation in the local rabbit population was found. We were able to sequence a short viral genomic fragment, confirming its identity as RHDV GI.2. The results suggest that scavenging flies may act as mechanical vectors of GI.2 in the native range of the southwestern Iberian subspecies O. cuniculus algirus. Future studies should better assess their potential in the epidemiology of RHD and as a tool for monitoring viral circulation in the field.

Evidence for Extensive Duplication and Subfunctionalization of FCRL6 in Armadillo (Dasypus novemcinctus).

The control of infections by the vertebrate adaptive immune system requires careful modulation to optimize defense and minimize harm to the host. The Fc receptor-like (FCRL) genes encode immunoregulatory molecules homologous to the receptors for the Fc portion of immunoglobulin (FCR). To date, nine different genes (FCRL1-6, FCRLA, FCRLB and FCRLS) have been identified in mammalian organisms. FCRL6 is located at a separate chromosomal position from the FCRL1-5 locus, has conserved synteny in mammals and is situated between the SLAMF8 and DUSP23 genes. Here, we show that this three gene block underwent repeated duplication in Dasypus novemcinctus (nine-banded armadillo) resulting in six FCRL6 copies, of which five appear functional. Among 21 mammalian genomes analyzed, this expansion was unique to D. novemcinctus. Ig-like domains that derive from the five clustered FCRL6 functional gene copies show high structural conservation and sequence identity. However, the presence of multiple non-synonymous amino acid changes that would diversify individual receptor function has led to the hypothesis that FCRL6 endured subfunctionalization during evolution in D. novemcinctus. Interestingly, D. novemcinctus is noteworthy for its natural resistance to the Mycobacterium leprae pathogen that causes leprosy. Because FCRL6 is chiefly expressed by cytotoxic T and NK cells, which are important in cellular defense responses against M. leprae, we speculate that FCRL6 subfunctionalization could be relevant for the adaptation of D. novemcinctus to leprosy. These findings highlight the species-specific diversification of FCRL family members and the genetic complexity underlying evolving multigene families critical for modulating adaptive immune protection.

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.

Identification of a Novel Myxoma Virus C7-Like Host Range Factor That Enabled a Species Leap from Rabbits to Hares.

Myxoma virus (MYXV) is naturally found in rabbit Sylvilagus species and is known to cause lethal myxomatosis in European rabbits (Oryctolagus cuniculus). In 2019, an MYXV strain (MYXV strain Toledo [MYXV-Tol]) causing myxomatosis-like disease in Iberian hares (Lepus granatensis) was identified. MYXV-Tol acquired a recombinant region of ∼2.8 kb harboring several new genes, including a novel host range gene (M159) that we show to be an orthologous member of the vaccinia virus C7 host range family. Here, to test whether M159 alone has enabled MYXV to alter its host range to Iberian hares, several recombinant viruses were generated, including an MYXV-Tol ΔM159 (knockout) strain. While MYXV-Tol underwent fully productive infection in hare HN-R cells, neither the wild-type MYXV-Lau strain (lacking M159) nor vMyxTol-ΔM159 (deleted for M159) was able to infect and replicate, showing that the ability of MYXV-Tol to infect these cells and replicate depends on the presence of M159. Similar to other C7L family members, M159 was shown to be expressed as an early/late gene but was translocated into the nucleus at later time points, indicating that further studies are needed to elucidate its role in the nucleus. Finally, in rabbit cells, the M159 protein did not contribute to increased replication but was able to upregulate the replication levels of MYXV in nonpermissive and semipermissive human cancer cells, suggesting that the M159-targeted pathway is conserved across mammalian species. Altogether, these observations demonstrate that the M159 protein plays a critical role in determining the host specificity of MYXV-Tol in hare and human cells by imparting new host range functions. IMPORTANCE The coevolution of European rabbit populations and MYXV is a textbook example of an arms race between a pathogen and a host. Recently, a recombinant MYXV (MYXV-Tol) crossed the species barrier by jumping from leporid species to another species, causing lethal myxomatosis-like disease. Given the highly pathogenic nature of this new virus in hares and the incidences of other poxvirus cross-species spillovers into other animals, including humans, it is important to understand how and why MYXV-Tol was able to become virulent in a new host species. The results presented clearly demonstrate that M159 is the key factor allowing MYXV-Tol replication in hare cells by imparting new host range functions. These results have the potential to improve current knowledge about the virulence of poxviruses and provide a platform to better understand the new MYXV-Tol, rendering the virus capable of leaping into a new host species.

Evolution of CCR5 and CCR2 Genes in Bats Showed Multiple Independent Gene Conversion Events.

Chemokine receptors are an important determinant for the infectiousness of different pathogens, which are able to target the host cells by binding to the extracellular domains of these proteins. This is the mechanism of infection of HIV-1, among other concerning human diseases. Over the past years, it has been shown that two chemokine receptors, CCR2 and CCR5, have been shaped by events of gene conversion in different mammalian lineages, which has been linked to a possible selective advantage against pathogens. Here, by taking advantage of available bat genomes, we present the first insight of CCR2 and CCR5 evolution within the Chiroptera order. In total, four independent events of recombination between CCR2 and CCR5 were detected: two in a single species, Miniopterus natalensis; one in two species from the Rhinolophoidea superfamily; and one in four species from the Pteropodidae family. The regions affected by the gene conversions were generally extensive and always encompassed extracellular domains. Overall, we demonstrate that CCR2 and CCR5 have been subject to extensive gene conversion in multiple species of bats. Considering that bats are known to be large reservoirs of virus in nature, these results might indicate that chimeric CCR2-CCR5 genes might grant some bat species a selective advantage against viruses that rely in the extracellular portions of either CCR2 or CCR5 as gateways into the cell.

Evolution of TRIM5 and TRIM22 in Bats Reveals a Complex Duplication Process.

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.

Uncovering a 500 million year old history and evidence of pseudogenization for TLR15.

Introduction: Toll like receptors (TLRs) are at the front line of pathogen recognition and host immune response. Many TLR genes have been described to date with some being found across metazoans while others are restricted to specific lineages. A cryptic member of the TLR gene family, TLR15, has a unique phylogenetic distribution. Initially described in extant species of birds and reptiles, an ortholog has been reported for cartilaginous fish. Methods: Here, we significantly expanded the evolutionary analysis of TLR15 gene evolution, taking advantage of large genomic and transcriptomic resources available from different lineages of vertebrates. Additionally, we objectively search for TLR15 in lobe-finned and ray-finned fish, as well as in cartilaginous fish and jawless vertebrates. Results and discussion: We confirm the presence of TLR15 in early branching jawed vertebrates - the cartilaginous fish, as well as in basal Sarcopterygii - in lungfish. However, within cartilaginous fish, the gene is present in Holocephalans (all three families) but not in Elasmobranchs (its sister-lineage). Holocephalans have long TLR15 protein sequences that disrupt the typical TLR structure, and some species display a pseudogene sequence due to the presence of frameshift mutations and early stop codons. Additionally, TLR15 has low expression levels in holocephalans when compared with other TLR genes. In turn, lungfish also have long TLR15 protein sequences but the protein structure is not compromised. Finally, TLR15 presents several sites under negative selection. Overall, these results suggest that TLR15 is an ancient TLR gene and is experiencing ongoing pseudogenization in early-branching vertebrates.

Convergent Loss of the Necroptosis Pathway in Disparate Mammalian Lineages Shapes Viruses Countermeasures.

Programmed cell death is a vital process in the life cycle of organisms. Necroptosis, an evolutionary form of programmed necrosis, contributes to the innate immune response by killing pathogen-infected cells. This virus-host interaction pathway is organized around two components: the receptor-interacting protein kinase 3 (RIPK3), which recruits and phosphorylates the mixed lineage kinase-like protein (MLKL), inducing cellular plasma membrane rupture and cell death. Critically, the presence of necroptotic inhibitors in viral genomes validates necroptosis as an important host defense mechanism. Here, we show, counterintuitively, that in different mammalian lineages, central components of necroptosis, such as RIPK3 and MLKL, are deleted or display inactivating mutations. Frameshifts or premature stop codons are observed in all the studied species of cetaceans and leporids. In carnivores' genomes, the MLKL gene is deleted, while in a small number of species from afrotheria and rodentia premature stop codons are observed in RIPK3 and/or MLKL. Interestingly, we also found a strong correlation between the disruption of necroptosis in leporids and cetaceans and the absence of the N-terminal domain of E3-like homologs (responsible for necroptosis inhibition) in their naturally infecting poxviruses. Overall, our study provides the first comprehensive picture of the molecular evolution of necroptosis in mammals. The loss of necroptosis multiple times during mammalian evolution highlights the importance of gene/pathway loss for species adaptation and suggests that necroptosis is not required for normal mammalian development. Moreover, this study highlights a co-evolutionary relationship between poxviruses and their hosts, emphasizing the role of host adaptation in shaping virus evolution.

A Highly Complex, MHC-Linked, 350 Million-Year-Old Shark Nonclassical Class I Lineage.

Cartilaginous fish, or Chondrichthyes, are the oldest extant vertebrates to possess the MHC and the Ig superfamily-based Ag receptors, the defining genes of the gnathostome adaptive immune system. In this work, we have identified a novel MHC lineage, UEA, a complex multigene nonclassical class I family found in sharks (division Selachii) but not detected in chimaeras (subclass Holocephali) or rays (division Batoidea). This new lineage is distantly related to the previously reported nonclassical class I lineage UCA, which appears to be present only in dogfish sharks (order Squaliformes). UEA lacks conservation of the nine invariant residues in the peptide (ligand)-binding regions (PBR) that bind to the N and C termini of bound peptide in most vertebrate classical class I proteins, which are replaced by relatively hydrophobic residues compared with the classical UAA. In fact, UEA and UCA proteins have the most hydrophobic-predicted PBR of all identified chondrichthyan class I molecules. UEA genes detected in the whale shark and bamboo shark genome projects are MHC linked. Consistent with UEA comprising a very large gene family, we detected weak expression in different tissues of the nurse shark via Northern blotting and RNA sequencing. UEA genes fall into three sublineages with unique characteristics in the PBR. UEA shares structural and genetic features with certain nonclassical class I genes in other vertebrates, such as the highly complex XNC nonclassical class I genes in Xenopus, and we anticipate that each shark gene, or at least each sublineage, will have a unique function, perhaps in bacterial defense.

Spillover event of recombinant Lagovirus europaeus/GI.2 into the Iberian hare (Lepus granatensis) in Spain.

Viruses that affect lagomorphs have decades of reported history of spillover events. One of these viruses is the causative agent of the so-called rabbit or 'lagomorph' haemorrhagic disease (e.g. Lagovirus europaeus/GI.1 and L. europaeus/GI.2). In particular, L. europaeus/GI.2 has shown a great capacity to recombine with existing lagoviruses. In fact, it has replaced the former GI.1 genotype in the wild, and recently, an increase on spillover events has been detected among several lagomorph species including European and North American species of hares. In this study, we report for the first time the infection of a wild Iberian hare with GI.2 (RHDV2/b), potential shedding and associated histopathological alterations. We identify the recombinant GI.4P-GI.2 as causative of the infection and discuss plausible causes regarding the origin of the spillover event and its potential consequences for the Iberian hare wild populations, which is an endemic species of the Iberian Peninsula as well as an important game and prey species for many predators, including endangered species.

Brain and testis: more alike than previously thought?

Several strands of evidence indicate the presence of marked similarities between human brain and testis. Understanding these similarities and their implications has become a topic of interest among the scientific community. Indeed, an association of intelligence with some semen quality parameters has been reported and a relation between dysfunctions of the human brain and testis has also been evident. Numerous common molecular features are evident when these tissues are compared, which is reflected in the huge number of common proteins. At the functional level, human neurons and sperm share a number of characteristics, including the importance of the exocytotic process and the presence of similar receptors and signalling pathways. The common proteins are mainly involved in exocytosis, tissue development and neuron/brain-associated biological processes. With this analysis, we conclude that human brain and testis share several biochemical characteristics which, in addition to their involvement in the speciation process, could, at least in part, be responsible for the expression of a huge number of common proteins. Nonetheless, this is an underexplored topic, and the connection between these tissues needs to be clarified, which could help to understand the dysfunctions affecting brain and testis, as well as to develop improved therapeutic strategies.

A loss-of-function mutation in RORB disrupts saltatorial locomotion in rabbits.

Saltatorial locomotion is a type of hopping gait that in mammals can be found in rabbits, hares, kangaroos, and some species of rodents. The molecular mechanisms that control and fine-tune the formation of this type of gait are unknown. Here, we take advantage of one strain of domesticated rabbits, the sauteur d'Alfort, that exhibits an abnormal locomotion behavior defined by the loss of the typical jumping that characterizes wild-type rabbits. Strikingly, individuals from this strain frequently adopt a bipedal gait using their front legs. Using a combination of experimental crosses and whole genome sequencing, we show that a single locus containing the RAR related orphan receptor B gene (RORB) explains the atypical gait of these rabbits. We found that a splice-site mutation in an evolutionary conserved site of RORB results in several aberrant transcript isoforms incorporating intronic sequence. This mutation leads to a drastic reduction of RORB-positive neurons in the spinal cord, as well as defects in differentiation of populations of spinal cord interneurons. Our results show that RORB function is required for the performance of saltatorial locomotion in rabbits.