Genetic variation at chemokine receptor CCR5 in leporids: alteration at the 2nd extracellular domain by gene conversion with CCR2 in Oryctolagus, but not in Sylvilagus and Lepus species

The rabbit as an animal model to study innate immunity genes: Is it better than mice?

The European rabbit (Oryctolagus cuniculus) was the first animal model used to understand human diseases like rabies and syphilis. Nowadays, the rabbit is still used to study several human infectious diseases like syphilis, HIV and papillomavirus. However, due to several mainly practical reasons, it has been replaced as an animal model by mice (Mus musculus). The rabbit and mouse share a recent common ancestor and are classified in the superorder Glires which arose at approximately 82 million years ago (mya). These species diverged from the Primates’ ancestor at around 92 million years ago and, as such, one expects the rabbit-human and mouse-human genetic distances to be very similar. To evaluate this hypothesis, we developed a set of tools for automatic data extraction, sequence alignment and similarity study, and a web application for visualization of the resulting data. We aligned and calculated the genetic distances for 2793 innate immune system genes from human, rabbit and mouse using sequences available in the NCBI database. The obtained results show that the rabbit-human genetic distance is lower than the mouse-human genetic distance for 88% of these genes. Furthermore, when we considered only genes with a difference in genetic distance higher than 0.05, this figure increase to 93%. These results can be explained by the increase of the mutation rates in the mouse lineage suggested by some authors and clearly show that, at least looking to the genetic distance to human genes, the European rabbit is a better model to study innate immune system genes than the mouse.

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.

Genome-wide identification and analysis of chemokine receptor superfamily in miiuy croaker, Miichthys miiuy

The chemokine receptor (ChemR) superfamily, which is divided into 4 subfamilies (CXCR, CCR, XCR, and CX3CR), is the main receptors of chemokines in innate immune responses. In the current study, we have identified 27 ChemRs in miiuy croaker: 13 CCR genes, 11 CXCR genes, and 3 XCR genes. Multiple characteristics of these genes, including phylogeny, gene structures, conserved motifs, chromosome locations, evolutionary mechanism, and expression levels upon the bacterial challenge were analyzed. Gene structure and location analysis showed that all ChemR genes contain fewer introns (≤4) and they are unevenly distributed on the 12 chromosomes. And the XCR subfamily of miiuy croaker don't have the DRY motif of ChemR. Phylogenetic and synteny analysis showed that these genes experienced tandem and segmental duplication event in several species, and tandem duplication might be the main expansion way in miiuy croaker. The major ChemRs of each orthologous group in vertebrates were selected for molecular evolution analysis, the results of which indicated that compared with vertebrates, ChemRs of teleost fishes may have a relatively high evolutionary dynamic. In addition, a total of 21 positively selected codons were detected in vertebrate ChemRs under Model 8. RNA-Seq analysis and qRT-PCR verification demonstrated that CXCR3.2, CXCR5, and XCR1 genes were up-regulated significantly upon the Vibrio harveyi infection. These results provide valuable information for investigating the evolutionary relationships of chemokine receptor superfamily in miiuy croaker and laid the basis for further functional analysis.

What Are the Keys to the Adaptive Success of European Wild Rabbit (Oryctolagus cuniculus) in the Iberian Peninsula?

The European wild rabbit (Oryctolagus cuniculus) plays an important ecological role in the ecosystems of the Iberian Peninsula. Recently, rabbit populations have drastically reduced, so the species is now considered endangered. However, in some places, this animal is considered a pest. This is the conservation paradox of the 21st century: the wild rabbit is both an invasive alien and an endangered native species. The authors of this review aimed to understand the keys to the adaptive success of European rabbits, addressing all aspects of their biology in order to provide the keys to the ecological management of this species. Aspects including nutrition, genetics, immunity interactions with the environment, behaviour, and conflict with human activities were reviewed. Ultimately, rabbits are resilient and adaptable. The main adaptations that explain the rabbit's adaptive success are its nutrition (wide adaptation to food and good nutritional use of caecotrophy), immune system (powerful and developed), and other aspects related to genetics and behaviour. Rabbits' relationship with humans has led them to colonise other places where they have become pests. Despite these adaptations, populations in native places have been drastically reduced in recent years. Since it serves as a bastion of the Mediterranean ecosystem, a specific conservation program for this species must be carried out. Therefore, a study of the rabbit's response to diseases and nutrition (especially protein), as well as the interaction between them, is of special interest.

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.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1390-7) contains supplementary material, which is available to authorized users.

Intragenus (Homo) variation in a chemokine receptor gene (CCR5)

Humans have a comparatively higher rate of more polymorphisms in regulatory regions of the primate CCR5 gene, an immune system gene with both general and specific functions. This has been interpreted as allowing flexibility and diversity of gene expression in response to varying disease loads. A broad expression repertoire is useful to humans-the only globally distributed primate-due to our unique adaptive pattern that increased pathogen exposure and disease loads (e.g., sedentism, subsistence practices). The main objective of the study was to determine if the previously observed human pattern of increased variation extended to other members of our genus, Homo. The data for this study are mined from the published genomes of extinct hominins (four Neandertals and two Denisovans), an ancient human (Ust'-Ishim), and modern humans (1000 Genomes). An average of 15 polymorphisms per individual were found in human populations (with a total of 262 polymorphisms). There were 94 polymorphisms identified across extinct Homo (an average of 13 per individual) with 41 previously observed in modern humans and 53 novel polymorphisms (32 in Denisova and 21 in Neandertal). Neither the frequency nor distribution of polymorphisms across gene regions exhibit significant differences within the genus Homo. Thus, humans are not unique with regards to the increased frequency of regulatory polymorphisms and the evolution of variation patterns across CCR5 gene appears to have originated within the genus. A broader evolutionary perspective on regulatory flexibility may be that it provided an advantage during the transition to confrontational foraging (and later hunting) that altered human-environment interaction as well as during migration to Eurasia and encounters with novel pathogens.

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.

Rabbits offer a powerful complement to rodents as a model for studying human immunology, disease pathology, and responses to infectious disease. A review from Pedro Esteves at the University of Porto, Portugal, Rose Mage of the National Institute of Allergy and Infectious Disease, Bethesda, USA and colleagues highlights some of the areas of research where rabbits offer an edge over rats and mice. Rabbits have a particularly sophisticated adaptive immune system, which could provide useful insights into human biology and produce valuable research and clinical reagents. They are also excellent models for studying - infectious diseases such as syphilis and tuberculosis, which produce pathology that closely resembles that of human patients. Rabbit-specific infections such as myxomatosis are giving researchers insights into how pathogens and hosts can shape each other’s evolution.

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

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

Mechanistic insights into molecular evolution of species-specific differential glycosaminoglycan binding surfaces in growth-related oncogene chemokines

Chemokines are chemotactic cytokines involved in leucocyte trafficking to infected tissue. Growth-related oncogene (GRO) chemokines namely CXCL1, CXCL2 and CXCL3 are neutrophil activating chemokines sharing a conserved three-dimensional structure, but encompassing functional diversity due to gene duplication and evolutionary events. However, the evolutionary mechanisms including selection pressures involved in diversification of GRO genes have not yet been characterized. Here, we performed comprehensive evolutionary analysis of GRO genes among different mammalian species. Phylogenetic analysis illustrated a species-specific evolution pattern. Selection analysis evidenced that these genes have undergone concerted evolution. Seventeen positively selected sites were obtained, although the majority of the protein is under purifying selection. Interestingly, these positively selected sites are more concentrated on the C-terminal/glycosaminoglycan (GAG) binding and dimerization segment compared to receptor binding domain. Substitution rate analysis confirmed the C-terminal domain of GRO genes as the highest substituted segment. Further, structural analysis established that the nucleotide alterations in the GAG binding domain are the source of surface charge modulation, thus generating the differential GAG binding surfaces and multiple binding sites as per evolutionary pressure, although the helical surface is primordial for GAG binding. Indeed, such variable electrostatic surfaces are crucial to regulate chemokine gradient formation during a host's defence against pathogens and also explain the significance of chemokine promiscuity.

Molecular characterization and analysis of TLR-1 in rabbit tissues

The rabbit has great commercial importance as a source of meat and fur, as well as its uses as a laboratory animal for the production of antibodies, used to detect the presence or absence of disease and for research in infectious diseases and immunology. One of the most critical problems in immunology is to understand how the immune system detects the presence of infectious agents and disposes the invader without destroying the self-tissues. Genetic characterization of Toll-like receptors has established that innate immunity is a skillful system that detects invasion of microbial pathogens. Our work aimed to identify, clone and express the Oryctolagus cuniculus (rabbit) TLR-1 mRNA and its encoding protein. We cloned the complete mRNA sequence of Oryctolagus cuniculus TLR-1 and deposit it in the GenBank under accession number (KC349941), which has 2388 base pair and it encodes encode an open reading frame (ORF) translated into 796 amino acids mRNA and consist of 20 types of amino acids. The analysis of amino acid sequence revealed that the rabbit TLR-1 has a typical protein components belonging to the TLR family. Rabbit TLR-1 was expressed in a wide variety of rabbit tissues, which indicate an important role in immune system in different organs.

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.

Evolution of CCL11: genetic characterization in lagomorphs and evidence of positive and purifying selection in mammals

The interactions between chemokines and their receptors are crucial for differentiation and activation of inflammatory cells. CC chemokine ligand 11 (CCL11) binds to CCR3 and to CCR5 that in leporids underwent gene conversion with CCR2. Here, we genetically characterized CCL11 in lagomorphs (leporids and pikas). All lagomorphs have a potentially functional CCL11, and the Pygmy rabbit has a mutation in the stop codon that leads to a longer protein. Other mammals also have mutations at the stop codon that result in proteins with different lengths. By employing maximum likelihood methods, we observed that, in mammals, CCL11 exhibits both signatures of purifying and positive selection. Signatures of purifying selection were detected in sites important for receptor binding and activation. Of the three sites detected as under positive selection, two were located close to the stop codon. Our results suggest that CCL11 is functional in all lagomorphs, and that the signatures of purifying and positive selection in mammalian CCL11 probably reflect the protein's biological roles.

Tissue-specific transcription profile of cytokine and chemokine genes associated with flavivirus control and non-lethal neuropathogenesis in rabbits

We previously showed that New Zealand White (NZWRs) and cottontail rabbits (CTRs) are a suitable model for studying immune mechanisms behind virus control and non-lethal neuropathogenesis associated with West Nile virus (WNV) and Murray Valley encephalitis virus (MVEV) infections. In the current study, we observed that MVEV infection induced high IFNα, TNFα, IL6, and CXCL10 transcript levels in the brains of weanling NZWRs, unlike infection with the less virulent WNVNSW2011. These transcript levels also correlated with encephalitis severity. Widespread STAT1 protein expression in brain with moderate neuropathology suggests that IFN-I signaling is crucial for limiting neural infection and mediating non-lethal neuropathogenesis. Unlike NZWRs, CTRs limit neuroinvasion without upregulation of many cytokine/chemokine transcripts, suggesting a species-dependent virus control mechanism. However, the common IFNγ, TNFα and IL6 transcript upregulation in specific lymphoid organs suggest some conserved elements in the response against flaviviruses, unique to all rabbits.

Chemokine receptor CXCR3 in turbot (Scophthalmus maximus): cloning, characterization and its responses to lipopolysaccharide

Chemokine (C-X-C motif) receptor 3, a member of the G protein-coupled receptors superfamily, regulates the responses of many immune responses. In this experiment, we cloned and characterized the cDNA of CXCR3 in Scophthalmus maximus (turbot). A 5'-UTR of 216-bp, a 259-bp 3'-UTR with a poly (A) tail and a 1089-bp CDS encoding 362 amino acids form the cDNA of CXCR3, which is 1564-bp long. Phylogenetic analyses indicated that turbot CXCR3 shared a high similarity with other CXCR3s and shared more similarity with CXCR5 than the other subfamilies of chemokines. The CXCR3 protein in turbot showed the highest similarity with the CXCR3b from rainbow trout (44.5%), which indicated that this CXCR3 gene/protein may be a CXCR3b isoform. Quantitative real-time PCR analysis showed that CXCR3 transcripts were constitutively expressed in all the tissues of the non-injected turbot used in this study, with the highest expression occurring in blood. Several immune-related tissues of fish, such as the spleen, head kidney, liver and blood, tissues, which were abundant of lymphocyte, were investigated in this study. CXCR3 gene was expressed at the highest level in blood than the other tested tissues. The injection experiment suggested that the CXCR3 expression level after LPS injection was significantly up-regulated in all immune-related tissues in turbot. These results improve our understanding of the functions of CXCR3 in the turbot immune response.

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.

Pseudogenization of CCL14 in the Ochotonidae (pika) family

The interaction between chemokines and their receptors is crucial for inflammatory cell trafficking. CCL14 binds with high affinity to CCR5. In leporids, CCR5 underwent gene conversion with CCR2. The study of CCR5 ligands in leporid species showed that CCL8 is pseudogenized, while CCL3, CCL4 and CCL5 are functional. Here, we study the evolution of CCL14 in mammals with emphasis in the order Lagomorpha. By employing maximum likelihood methods we detected six sites under positive selection. Some of these sites are located in regions crucial for CCL14 activation and binding to receptors. Sequencing of CCL14 in Ochotona species showed that O. princeps, O. pallasi, O. alpina and O. turuchanensis have a mutation at the start codon (Met > Thr), while O. hoffmanni, O. mantchurica, O. dauurica and O. rufescens present the mammalian conserved Met. Ochotona hyperborea has the two alleles. In O. pusilla, CCL14 is a pseudogene due to a seven base pair insertion. Like CCL3, CCL4 and CCL5, CCL14 is functional in all leporids but in the Ochotonidae family it underwent a pseudogenization process. This suggests that CCL14 has an important biological role in other mammals by evolving under positive selection that has been lost in Ochotonidae (subgenera Pika and Lagotona).

Field and experimental data indicate that the eastern cottontail (Sylvilagus floridanus) is susceptible to infection with European brown hare syndrome (EBHS) virus and not with rabbit haemorrhagic disease (RHD) virus

The eastern cottontail (Sylvilagus floridanus) is an American lagomorph. In 1966, it was introduced to Italy, where it is currently widespread. Its ecological niche is similar to those of native rabbits and hares and increasing overlap in distribution brings these species into ever closer contact. Therefore, cottontails are at risk of infection with the two lagoviruses endemically present in Italy: Rabbit Haemorrhagic Disease virus (RHDV) and European Brown Hare Syndrome Virus (EBHSV). To verify the susceptibility of Sylvilagus to these viruses, we analyzed 471 sera and 108 individuals from cottontail populations in 9 provinces of north-central Italy from 1999 to 2012. In total, 15–20% of the cottontails tested seropositive for EBHSV; most titres were low, but some were as high as 1/1280. All the cottontails virologically tested for RHDV and EBHSV were negative with the exception of one individual found dead with hares during a natural EBHS outbreak in December 2009. The cottontail and the hares showed typical EBHS lesions, and the EBHSV strain identified was the same in both species (99.9% identity). To experimentally confirm the diagnosis, we performed two trials in which we infected cottontails with both EBHSV and RHDV. One out of four cottontails infected with EBHSV died of an EBHS-like disease, and the three surviving animals developed high EBHSV antibody titres. In contrast, neither mortality nor seroconversion was detected after infection with RHDV. Taken together, these results suggest that Sylvilagus is susceptible to EBHSV infection, which occasionally evolves to EBHS-like disease; the eastern cottontail could therefore be considered a “spill over” or “dead end” host for EBHSV unless further evidence is found to confirm that it plays an active role in the epidemiology of EBHSV.

Detection of RHDV strains in the Iberian hare (Lepus granatensis): earliest evidence of rabbit lagovirus cross-species infection

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.

An intriguing shift occurs in the novel protein phosphatase 1 binding partner, TCTEX1D4: evidence of positive selection in a pika model

T-complex testis expressed protein 1 domain containing 4 (TCTEX1D4) contains the canonical phosphoprotein phosphatase 1 (PPP1) binding motif, composed by the amino acid sequence RVSF. We identified and validated the binding of TCTEX1D4 to PPP1 and demonstrated that indeed this protein is a novel PPP1 interacting protein. Analyses of twenty-one mammalian species available in public databases and seven Lagomorpha sequences obtained in this work showed that the PPP1 binding motif ₉₀RVSF₉₃ is present in all of them and is flanked by a palindromic sequence, PLGS, except in three species of pikas (Ochotona princeps, O. dauurica and O. pusilla). Furthermore, for the Ochotona species an extra glycosylation site, motif ₉₆NLS₉₈, and the loss of the palindromic sequence were observed. Comparison with other lagomorphs suggests that this event happened before the Ochotona radiation. The d_N/d_S for the sequence region comprising the PPP1 binding motif and the flanking palindrome highly supports the hypothesis that for Ochotona species this region has been evolving under positive selection. In addition, mutational screening shows that the ability of pikas TCTEX1D4 to bind to PPP1 is maintained, although the PPP1 binding motif is disrupted, and the N- and C-terminal surrounding residues are also abrogated. These observations suggest pika as an ideal model to study novel PPP1 complexes regulatory mechanisms.

Genetic characterization of CCL3, CCL4 and CCL5 in leporid genera Oryctolagus, Sylvilagus and Lepus

The genetic diversity of C-C motif chemokine receptor 5 (CCR5) ligands CCL3, CCL4 and CCL5 in the leporid genera Oryctolagus, Sylvilagus and Lepus was studied. Our results demonstrate that the three CCR5 chemokine ligands are under strong purifying selection as a result of possible functional binding constraints.

Insights into the European rabbit (Oryctolagus cuniculus) innate immune system: genetic diversity of the toll-like receptor 3 (TLR3) in wild populations and domestic breeds

Background

Toll-like receptors (TLRs) belong to the innate immune system and are a major class of pattern recognition receptors representing the first line of the innate immune response. The TLR molecule is structurally composed by an ectodomain that contains leucine rich repeats (LRRs) that interact with pathogen associated molecular patterns (PAMPs), a transmembrane domain and a conserved cytoplasmic domain designated TIR (Toll-IL1 receptor) that is responsible for the intracellular signaling. TLR3 has been associated with the direct recognition of double-stranded viral RNA resulting from viral replication, while TLR7 and TLR8 target single-stranded viral RNA. In the European rabbit (Oryctolagus cuniculus), TLR7 and TLR8 were reported to be absent and pseudogenised, respectively, making TLR3 the only available TLR for the recognition of viral RNA. Thus, the levels of diversity of TLR3 were evaluated in the European rabbit by analysing different genetic backgrounds and exposure to pathogen pressures.

Results

We detected 41 single nucleotide polymorphisms (SNPs) in the coding sequence of TLR3. The highest diversity was observed in the wild populations of Iberian Peninsula, between 22–33 polymorphic positions. In the French population, 18 SNPs were observed and only 4 polymorphic positions were detected in the domestic breeds. 14 non-synonymous substitutions were observed, most of them in the LRR molecules. The remaining were scattered across the transmembrane and TIR domains.

Conclusion

The study of TLR3 in European rabbit populations might be relevant to understand the interplay between RNA viruses and innate immunity. Wild rabbit populations presented more diversity than domestic breeds and other mammals previously studied. This might be linked to the absence of population bottlenecks during their evolution and to the almost inexistence of man-mediated selection. The observed variability might have also been potentiated by the contact of the wild populations with various pathogens. The study of these patterns of variability might reveal scenarios of host-pathogen interaction and identify TLR3 polymorphisms’ that arose due to viral pathogens affecting wild populations.

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.

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.

A family tree of vertebrate chemokine receptors for a unified nomenclature

Chemokines receptors are involved in the recruitment of various cell types in inflammatory and physiological conditions. There are 23 known chemokine receptor genes in the human genome. However, it is still unclear how many chemokine receptors exist in the genomes of various vertebrate species other than human and mouse. Moreover, the orthologous relationships are often obscure between the genes of higher and lower vertebrates. In order to provide a basis for a unified nomenclature system of the vertebrate chemokine receptor gene family, we have analysed the chemokine receptor genes from the genomes of 16 vertebrate species, and classify them into 29 orthologous groups using phylogenetic and comparative genomic analyses. The results reveal a continuous gene birth and death process during the vertebrate evolution and an interesting evolutionary history of the chemokine receptor genes after the emergence in agnathans.

Evidence of positive selection at codon sites localized in extracellular domains of mammalian CC motif chemokine receptor proteins

Background

CC chemokine receptor proteins (CCR1 through CCR10) are seven-transmembrane G-protein coupled receptors whose signaling pathways are known for their important roles coordinating immune system responses through targeted trafficking of white blood cells. In addition, some of these receptors have been identified as fusion proteins for viral pathogens: for example, HIV-1 strains utilize CCR5, CCR2 and CCR3 proteins to obtain cellular entry in humans. The extracellular domains of these receptor proteins are involved in ligand-binding specificity as well as pathogen recognition interactions.

In mammals, the majority of chemokine receptor genes are clustered together; in humans, seven of the ten genes are clustered in the 3p21-24 chromosome region. Gene conversion events, or exchange of DNA sequence between genes, have been reported in chemokine receptor paralogs in various mammalian lineages, especially between the cytogenetically closely located pairs CCR2/5 and CCR1/3. Datasets of mammalian orthologs for each gene were analyzed separately to minimize the potential confounding impact of analyzing highly similar sequences resulting from gene conversion events.

Molecular evolution approaches and the software package Phylogenetic Analyses by Maximum Likelihood (PAML) were utilized to investigate the signature of selection that has acted on the mammalian CC chemokine receptor (CCR) gene family. The results of neutral vs. adaptive evolution (positive selection) hypothesis testing using Site Models are reported. In general, positive selection is defined by a ratio of nonsynonymous/synonymous nucleotide changes (dN/dS, or ω) >1.

Results

Of the ten mammalian CC motif chemokine receptor sequence datasets analyzed, only CCR2 and CCR3 contain amino acid codon sites that exhibit evidence of positive selection using site based hypothesis testing in PAML. Nineteen of the twenty codon sites putatively indentified as likely to be under positive selection code for amino acid residues located in extracellular domains of the receptor protein products.

Conclusions

These results suggest that amino acid residues present in intracellular and membrane-bound domains are more selectively constrained for functional signal transduction and homo- or heterodimerization, whereas amino acid residues in extracellular domains of these receptor proteins evolve more quickly, perhaps due to heightened selective pressure resulting from ligand-binding and pathogen interactions of extracellular domains.

CCR5 polymorphism and plague resistance in natural populations of the black rat in Madagascar

Madagascar remains one of the world's largest plague foci. The black rat, Rattus rattus, is the main reservoir of plague in rural areas. This species is highly susceptible to plague in plague-free areas (low-altitude regions), whereas rats from the plague focus areas (central highlands) have evolved a disease-resistance polymorphism. We used the candidate gene CCR5 to investigate the genetic basis of plague resistance in R. rattus. We found a unique non-synonymous substitution (H184R) in a functionally important region of the gene. We then compared (i) CCR5 genotypes of dying and surviving plague-challenged rats and (ii) CCR5 allelic frequencies in plague focus and plague-free populations. Our results suggested a higher prevalence of the substitution in resistant animals compared to susceptible individuals, and a tendency for higher frequencies in plague focus areas compared to plague-free areas. Therefore, the CCR5 polymorphism may be involved in Malagasy black rat plague resistance. CCR5 and other undetermined plague resistance markers may provide useful biological information about host evolution and disease dynamics.

Genetic characterization of the chemokine receptor CXCR4 gene in lagomorphs: comparison between the families Ochotonidae and Leporidae

Chemokines receptors are transmembrane proteins that bind chemokines. Chemokines and their receptors are known to play a crucial role in the immune system and in pathogen entry. There is evidence that myxoma virus, the causative agent of myxomatosis, can use the chemokine receptor CXCR4 to infect cells. This virus causes a benign disease in its natural host, Sylvilagus, but in the European rabbit (Oryctolagus cuniculus) it causes a highly fatal and infectious disease known as myxomatosis. We have characterized the chemokine receptor CXCR4 gene in five genera of the order Lagomorpha, Ochotona (Ochotonidae), and Oryctolagus, Lepus, Bunolagus and Sylvilagus (Leporidae). In lagomorphs, the CXCR4 is highly conserved, with most of the protein diversity found at surface regions. Five amino acid replacements were observed, two in the intracellular loops, one in the transmembrane domain and two in the extracellular loops. Oryctolagus features unique amino acid changes at the intracellular domains, putting this genus apart of all other lagomorphs. Furthermore, in the 37 European rabbits analysed, which included healthy rabbits and rabbits with clinical symptoms of myxomatosis, 14 nucleotide substitutions were obtained but no amino acid differences were observed.

Extensive gene conversion between CCR2 and CCR5 in domestic cat (Felis catus)

Homogenization of the CC-motif chemokine receptors CCR2 and CCR5 of cat (Felis catus) is documented and shown to be the outcome of gene conversion within the feline lineage. All regions were concerned, except the three extracellular protein domains (N- and C-tails, and ECL2), suggesting that structural differentiation at these domains could be related to pathogen susceptibility.