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.

Characterization of old RHDV strains by complete genome sequencing identifies a novel genetic group.

Rabbit hemorrhagic disease (RHD) is a veterinary disease that affects the European rabbit and has a significant economic and ecological negative impact. In Portugal, rabbit hemorrhagic disease virus (RHDV) was reported in 1989 and still causes enzootic outbreaks. Several recombination events have been detected in RHDV strains, including in the first reported outbreak. Here we describe the occurrence of recombination in RHDV strains recovered from rabbit and Iberian hare samples collected in the mid-1990s in Portugal. Characterization of full genomic sequences revealed the existence of a single recombination breakpoint at the boundary of the non-structural and the structural encoding regions, further supporting the importance of this region as a recombination hotspot in lagoviruses. Phylogenetic analysis showed that in the structural region, the recombinant strains were similar to pathogenic G1 strains, but in the non-structural region they formed a new group that diverged ~13% from known strains. No further reports of such group exist, but this recombination event was also detected in an Iberian hare that was associated with the earliest species jump in RHDV. Our results highlight the importance of the characterization of full genomes to disclose RHDV evolution and show that lagoviruses' diversity has been significantly undersampled.

The phylogeny of pikas (Ochotona) inferred from a multilocus coalescent approach.

The clarification of the systematics of pikas (genus Ochotona) has been hindered by largely overlapping morphological characters among species and the lack of a comprehensive molecular phylogeny. Here we estimate the first multilocus phylogeny of the genus to date, by analysing 12 nuclear DNA markers (total of 7.5Kb) in 11 species of pikas from the four classified subgenera (Pika, Ochotona, Lagotona and Conothoa) using a multispecies coalescent-based framework. The species-tree confirmed the subgeneric classification by retrieving as monophyletic the subgenera represented here by more than one species. Contrary to previous phylogenies based on mtDNA alone, Lagotona was found to be sister to Pika. Also, support for the monophyly of the alpina group was not strong, thus caution should be used in future analyses of this group. A relaxed molecular clock calibrated using the Ochotonidae-Leporidae divergence resulted in more recent estimates of divergence times relative to previous studies. Strong concordance with inferences based on fossil records was found, suggesting that the initial diversification of the genus took place by the end of late Miocene. Finally, this work sets up methodologies and gathers molecular markers that can be used to extend the understanding of the evolutionary history of the genus.

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.

Signatures of positive selection in Toll-like receptor (TLR) genes in mammals.

BACKGROUND: Toll-like receptors (TLRs) are a major class of pattern recognition receptors (PRRs) expressed in the cell surface or membrane compartments of immune and non-immune cells. TLRs are encoded by a multigene family and represent the first line of defense against pathogens by detecting foreigner microbial molecular motifs, the pathogen-associated molecular patterns (PAMPs). TLRs are also important by triggering the adaptive immunity in vertebrates. They are characterized by the presence of leucine-rich repeats (LRRs) in the ectodomain, which are associated with the PAMPs recognition. The direct recognition of different pathogens by TLRs might result in different evolutionary adaptations important to understand the dynamics of the host-pathogen interplay. Ten mammal TLR genes, viral (TLR3, 7, 8, 9) and non-viral (TLR1-6, 10), were selected to identify signatures of positive selection that might have been imposed by interacting pathogens and to clarify if viral and non-viral TLRs might display different patterns of molecular evolution. RESULTS: By using Maximum Likelihood approaches, evidence of positive selection was found in all the TLRs studied. The number of positively selected codons (PSC) ranged between 2-26 codons (0.25%-2.65%) with the non-viral TLR4 as the receptor with higher percentage of positively selected codons (2.65%), followed by the viral TLR8 (2.50%). The results indicated that viral and non-viral TLRs are similarly under positive selection. Almost all TLRs have at least one PSC located in the LRR ectodomain which underlies the importance of the pathogen recognition by this region. CONCLUSIONS: Our results are not in line with previous studies on primates and birds that identified more codons under positive selection in non-viral TLRs. This might be explained by the fact that both primates and birds are homogeneous groups probably being affected by only a restricted number of related viruses with equivalent motifs to be recognized. The analyses performed in this work encompassed a large number of species covering some of the most representative mammalian groups - Artiodactyla, Rodents, Carnivores, Lagomorphs and Primates - that are affected by different families of viruses. This might explain the role of adaptive evolution in shaping viral TLR genes.

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.