Identification and characterization of a polyomavirus in the thornback skate (Raja clavata).

Members of the family Polyomaviridae have a circular double-stranded DNA genome that have been identified in various hosts ranging from mammals to arachnids. Here we report the identification and analysis of a complete genome sequence of a novel polyomavirus, Raja clavata polyomavirus (RcPyV1), from a cartilaginous fish, the thornback skate (Raja clavata). The genome sequence was determined using a metagenomics approach with an aim to provide baseline viral data in cartilaginous fish in different ecosystems. The RcPyV1 genome (4,195 nucleotides) had typical organization of polyomavirus, including early antigens (small T; Large T) encoded on one strand and late viral proteins (VP1; VP2) on the complementary strand. Maximum-likelihood phylogenetic analysis of the large T-antigen revealed that RcPyV1 clusters with a polyomavirus obtained from another cartilaginous fish, the guitarfish polyomavirus 1 (GfPyV1). These two share ~ 56% pairwise identity in LT and VP1 protein sequences. These analyses support the hypothesis that cartilaginous fishes have a specific lineage of polyomaviruses.

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.

Functional cross-species conservation of guanylate-binding proteins in innate immunity.

Guanylate binding proteins (GBPs) represent an evolutionary ancient protein family widely distributed among eukaryotes. They are interferon (IFN)-inducible guanosine triphosphatases that belong to the dynamin superfamily. GBPs are known to have a major role in the cell-autonomous innate immune response against bacterial, parasitic and viral infections and are also involved in inflammasome activation. Evolutionary studies depicted that GBPs present a pattern of gain and loss of genes in each family with several genes pseudogenized and some genes more divergent, indicative for the birth-and-death evolution process. Most species harbor large GBP gene clusters encoding multiple paralogs. Previous functional studies mainly focused on mouse and human GBPs, but more data are becoming available, broadening the understanding of this multifunctional protein family. In this review, we will provide new insights and give a broad overview about GBP evolution, conservation and their roles in all studied species, including plants, invertebrates and vertebrates, revealing how far the described features of GBPs can be transferred to other species.

TLR7 and TLR8 evolution in lagomorphs: different patterns in the different lineages.

Toll-like receptors (TLRs) are one of the most ancient and widely studied innate immune receptors responsible for host defense against invading pathogens. Among the known TLRs, TLR7 and TLR8 sense and recognize single-stranded (ss) RNAs with a dynamic evolutionary history. While TLR8 was lost in birds and duplicated in turtles and crocodiles, TLR7 is duplicated in some birds, but in other tetrapods, there is only one copy. In mammals, with the exception of lagomorphs, TLR7 and TLR8 are highly conserved. Here, we aim to study the evolution of TLR7 and TLR8 in mammals, with a special focus in the order Lagomorpha. By searching public sequence databases, conducting evolutionary analysis, and evaluating gene expression, we were able to confirm that TLR8 is absent in hares but widely expressed in the European rabbit. In contrast, TLR7 is absent in the European rabbit and quite divergent in hares. Our results suggest that, in lagomorphs, more in particular in leporids, TLR7 and TLR8 genes have evolved faster than in any other mammalian group. The long history of interaction with viruses and their location in highly dynamic telomeric regions might explain the pattern observed.

Real-world intensive safety monitoring of biosimilars rituximab and trastuzumab in a Portuguese oncology hospital.

PURPOSE: The aim of this study was to assess the safety profiles of two biosimilar medicines (rituximab and trastuzumab) in the treatment of cancer patients within a Portuguese oncology hospital. METHODS: This hospital-based prospective observational study followed a cohort event monitoring approach focused on signalling suspected adverse drug reactions (ADRs). Patients undergoing treatment with rituximab biosimilar CT-P10 (Truxima®) or trastuzumab biosimilar CT-P6 (Herzuma®) were recruited over an 11-month and a 6-month period, respectively. Clinicians identified eligible patients and used paper-based forms to report all ADRs associated with biosimilar medicines. ADR case reports were assessed for seriousness, expectedness and causality in the Pharmacovigilance Unit of Coimbra. RESULTS: Ninety-four patients received biosimilar medicines (rituximab, n = 35; trastuzumab, n = 59). Of those, 4 patients (11.4%) experienced 16 ADRs with rituximab and 1 patient (1.7%) experienced 5 ADRs with trastuzumab. All case reports contained serious and expected ADRs that were at least probably related with biosimilar medicines under study. Based on the MedDRA PT coding, the most reported ADR for rituximab CT-P10 was chest discomfort (n = 4; 19.1%), followed by odynophagia (n = 2; 9.5%). Trastuzumab CT-P6 was associated with back pain, headache, pain in extremity, tachypnoea and tremor (each, n = 1; 4.8%). CONCLUSION: The results of this study suggest that using biosimilar rituximab and biosimilar trastuzumab to treat cancer patients in the real-world clinical setting is associated with acceptable safety profiles. No new safety problems were identified.

Analysis of substitution rates showed that TLR5 is evolving at different rates among mammalian groups.

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.

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

Strong selection of the TLR2 coding region among the Lagomorpha suggests an evolutionary history that differs from other mammals.

Toll-like receptors (TLRs) are one of the first lines of defense against pathogens and are crucial for triggering an appropriate immune response. Among TLRs, TLR2 is functional in all vertebrates and has high ability in detecting bacterial and viral pathogen ligands. The mammals' phylogenetic tree of TLR2 showed longer branches for the Lagomorpha clade, raising the hypothesis that lagomorphs experienced an acceleration of the mutation rate. This hypothesis was confirmed by (i) Tajima's test of neutrality that revealed different evolutionary rates between lagomorphs and the remaining mammals with lagomorphs presenting higher nucleotide diversity; (ii) genetic distances were similar among lagomorphs and between lagomorphs and other mammals; and (iii) branch models reinforced the existence of an acceleration of the mutation rate in lagomorphs. These results suggest that the lagomorph TLR2 has been strongly involved in pathogen recognition, which probably caused a host-pathogen arms race that led to the observed acceleration of the mutation rate.

Host-Specific Glycans Are Correlated with Susceptibility to Infection by Lagoviruses, but Not with Their Virulence.

Rabbit hemorrhagic disease virus (RHDV) and European brown hare syndrome virus (EBHSV) are two lagoviruses from the family Caliciviridae that cause fatal diseases in two leporid genera, Oryctolagus and Lepus, respectively. In the last few years, several examples of host jumps of lagoviruses among leporids were recorded. In addition, a new pathogenic genotype of RHDV emerged, and many nonpathogenic strains of lagoviruses have been described. The molecular mechanisms behind host shifts and the emergence of virulence are unknown. Since RHDV uses glycans of the histo-blood group antigen type as attachment factors to initiate infection, we studied if glycan specificities of the new pathogenic RHDV genotype, nonpathogenic lagoviruses, and EBHSV potentially play a role in determining the host range and virulence of lagoviruses. We observed binding to A, B, or H antigens of the histo-blood group family for all strains known to primarily infect European rabbits (Oryctolagus cuniculus), which have recently been classified as GI strains. However, we could not explain the emergence of virulence, since similar glycan specificities were found in several pathogenic and nonpathogenic strains. In contrast, EBHSV, recently classified as GII.1, bound to terminal ß-linked N-acetylglucosamine residues of O-glycans. Expression of these attachment factors in the upper respiratory and digestive tracts in three lagomorph species (Oryctolagus cuniculus, Lepuseuropaeus, and Sylvilagus floridanus) showed species-specific patterns regarding susceptibility to infection by these viruses, indicating that species-specific glycan expression is likely a major contributor to lagovirus host specificity and range.IMPORTANCE Lagoviruses constitute a genus of the family Caliciviridae comprising highly pathogenic viruses, RHDV and EBHSV, that infect rabbits and hares, respectively. Recently, nonpathogenic strains were discovered and new pathogenic strains have emerged. In addition, host jumps between lagomorphs have been observed. The mechanisms responsible for the emergence of pathogenicity and host species range are unknown. Previous studies showed that RHDV strains attach to glycans expressed in the upper respiratory and digestive tracts of rabbits, the likely portals of virus entry. Here, we studied the glycan-binding properties of novel pathogenic and nonpathogenic strains looking for a link between glycan binding and virulence or between glycan specificity and host range. We found that glycan binding did not correlate with virulence. However, expression of glycan motifs in the upper respiratory and digestive tracts of lagomorphs revealed species-specific patterns associated with the host ranges of the virus strains, suggesting that glycan diversity contributes to lagovirus host ranges.

GI.1b/GI.1b/GI.2 recombinant rabbit hemorrhagic disease virus 2 (Lagovirus europaeus/GI.2) in Morocco, Africa.

Rabbit hemorrhagic disease virus (RHDV) is highly lethal to the European rabbit (Oryctolagus cuniculus). It was first reported in 1984 in China, but in 2010, a new variant of the virus was detected (GI.2) in France. Several recombination events with pathogenic and non-pathogenic strains have been described. Here, we report the first sequences of RHDV in Africa, isolated from Moroccan rabbits, and these resemble GI.1b/GI.1b/GI.2 recombinants found in the Iberian Peninsula. Monitoring and characterization of strains from future outbreaks are advised to guarantee the success of current programs on small-rabbit production for poverty alleviation in African countries.

Full genome sequences are key to disclose RHDV2 emergence in the Macaronesian islands.

A recent publication by Carvalho et al. in "Virus Genes" (June 2017) reported the presence of the new variant of rabbit hemorrhagic disease virus (RHDV2) in the two larger islands of the archipelago of Madeira. Based on the capsid protein sequence, the authors suggested that the high sequence identity, along with the short time spanning between outbreaks, points to dissemination from Porto Santo to Madeira. By including information of the full RHDV2 genome of strains from Azores, Madeira, and the Canary Islands, we confirm the results obtained by Carvalho et al., but further show that several subtypes of RHDV2 circulate in these islands: non-recombinant RHDV2 in the Canary Islands, G1/RHDV2 in Azores, Porto Santo and Madeira, and NP/RHDV2 also in Madeira. Here we conclude that RHDV2 has been independently introduced in these archipelagos, and that in Madeira at least two independent introductions must have occurred. We provide additional information on the dynamics of RHDV2 in the Macaronesian archipelagos of Azores, Madeira, and the Canary Islands and highlight the importance of analyzing RHDV2 complete genome.

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.

Proposal for a unified classification system and nomenclature of lagoviruses.

Lagoviruses belong to the Caliciviridae family. They were first recognized as highly pathogenic viruses of the European rabbit (Oryctolagus cuniculus) and European brown hare (Lepus europaeus) that emerged in the 1970-1980s, namely, rabbit haemorrhagic disease virus (RHDV) and European brown hare syndrome virus (EBHSV), according to the host species from which they had been first detected. However, the diversity of lagoviruses has recently expanded to include new related viruses with varying pathogenicity, geographic distribution and host ranges. Together with the frequent recombination observed amongst circulating viruses, there is a clear need to establish precise guidelines for classifying and naming lagovirus strains. Therefore, here we propose a new nomenclature based on phylogenetic relationships. In this new nomenclature, a single species of lagovirus would be recognized and called Lagovirus europaeus. The species would be divided into two genogroups that correspond to RHDV- and EBHSV-related viruses, respectively. Genogroups could be subdivided into genotypes, which could themselves be subdivided into phylogenetically well-supported variants. Based on available sequences, pairwise distance cutoffs have been defined, but with the accumulation of new sequences these cutoffs may need to be revised. We propose that an international working group could coordinate the nomenclature of lagoviruses and any proposals for revision.

Neofunctionalization of the Sec1 α1,2fucosyltransferase paralogue in leporids contributes to glycan polymorphism and resistance to rabbit hemorrhagic disease virus.

RHDV (rabbit hemorrhagic disease virus), a virulent calicivirus, causes high mortalities in European rabbit populations (Oryctolagus cuniculus). It uses α1,2fucosylated glycans, histo-blood group antigens (HBGAs), as attachment factors, with their absence or low expression generating resistance to the disease. Synthesis of these glycans requires an α1,2fucosyltransferase. In mammals, there are three closely located α1,2fucosyltransferase genes rSec1, rFut2 and rFut1 that arose through two rounds of duplications. In most mammalian species, Sec1 has clearly become a pseudogene. Yet, in leporids, it does not suffer gross alterations, although we previously observed that rabbit Sec1 variants present either low or no activity. Still, a low activity rSec1 allele correlated with survival to an RHDV outbreak. We now confirm the association between the α1,2fucosyltransferase loci and survival. In addition, we show that rabbits express homogenous rFut1 and rFut2 levels in the small intestine. Comparison of rFut1 and rFut2 activity showed that type 2 A, B and H antigens recognized by RHDV strains were mainly synthesized by rFut1, and all rFut1 variants detected in wild animals were equally active. Interestingly, rSec1 RNA levels were highly variable between individuals and high expression was associated with low binding of RHDV strains to the mucosa. Co-transfection of rFut1 and rSec1 caused a decrease in rFut1-generated RHDV binding sites, indicating that in rabbits, the catalytically inactive rSec1 protein acts as a dominant-negative of rFut1. Consistent with neofunctionalization of Sec1 in leporids, gene conversion analysis showed extensive homogenization between Sec1 and Fut2 in leporids, at variance with its limited degree in other mammals. Gene conversion additionally involving Fut1 was also observed at the C-terminus. Thus, in leporids, unlike in most other mammals where it became extinct, Sec1 evolved a new function with a dominant-negative effect on rFut1, contributing to fucosylated glycan diversity, and allowing herd protection from pathogens such as RHDV.

An update on the rabbit hemorrhagic disease virus (RHDV) strains circulating in Portugal in the 1990s: earliest detection of G3-G5 and G6.

Rabbit hemorrhagic disease virus (RHDV) causes devastating effects on European rabbit (Oryctolagus cuniculus) populations in the Iberian Peninsula. According to the information available, only genogroup 1 strains were circulating in Iberian wild rabbits until 2011; the antigenic variant G6 has been sporadically detected in rabbitries since 2007. Here, we show for the first time that G3-G5 strains were already present in mainland Portugal in 1998 and that G6 has been circulating since at least 1999. Moreover, we report a G3-G5 strain from the Azores collected in 1998, which is the likely ancestor of Azorean G3-G5like strains. These observations improve the current knowledge on RHDV epidemiology in the Iberian Peninsula and the Azores.

Recombination between G2 and G6 strains of rabbit hemorrhagic disease virus (RHDV) in China.

Rabbit hemorrhagic disease (RHD) is an acute fatal disease caused by the lagovirus rabbit hemorrhagic disease virus (RHDV), which was first reported in 1984 in China. Genetic characterization of RHDV has demonstrated that two different genogroups (G2 and G6) are present in China. To gain a better understanding of the molecular evolution of RHDV, we searched for recombination events by analyzing all full-length RHDV capsid VP60 sequences of Chinese isolates belonging to the genogroups 2 and 6. Our results revealed a recombinant origin for the NanBu/China/2011 isolate. This recombination event occurred between G2 and G6 strains with two breakpoints located at nucleotide positions 393 and 1079 of the VP60 sequence. Phylogenetically, the NanBu/China/2011 strain clustered with genogroup G6 in the entire capsid gene sequence except in the fragment between nucleotides 394 and 1078, where it clustered with genogroup G2. As the consequences of the presence of a G2/G6 recombinant strain in China are unpredictable, the circulation of RHDV in the populations should be carefully monitored.

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.

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.

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.