Contrasting patterns of polymorphism and selection in bacterial-sensing toll-like receptor 4 in two house mouse subspecies

Phenogenomic resources immortalized in a panel of wild-derived strains of five species of house mice

The house mouse, Mus musculus, is a widely used animal model in biomedical research, with classical laboratory strains (CLS) being the most frequently employed. However, the limited genetic variability in CLS hinders their applicability in evolutionary studies. Wild-derived strains (WDS), on the other hand, provide a suitable resource for such investigations. This study quantifies genetic and phenotypic data of 101 WDS representing 5 species, 3 subspecies, and 8 natural Y consomic strains and compares them with CLS. Genetic variability was estimated using whole mtDNA sequences, the Prdm9 gene, and copy number variation at two sex chromosome-linked genes. WDS exhibit a large natural variation with up to 2173 polymorphic sites in mitogenomes, whereas CLS display 92 sites. Moreover, while CLS have two Prdm9 alleles, WDS harbour 46 different alleles. Although CLS resemble M. m. domesticus and M. m. musculus WDS, they differ from them in 10 and 14 out of 16 phenotypic traits, respectively. The results suggest that WDS can be a useful tool in evolutionary and biomedical studies with great potential for medical applications.

Temporal dynamics of mildly deleterious nonsynonymous substitutions in mitochondrial gene sequences in rodents and moles

We have previously estimated the evolutionary rate (number of substitutions/site/million years) of mitochondrial cytochrome b gene (Cytb) sequences in rodents and moles to be about 0.11 at more recent divergence times of tens of thousands of years, and to decrease rapidly to about 0.03 at more distant divergence times. Because this time dependency is thought to be caused by the removal of mildly deleterious substitutions in later generations, we focused in this study on the abundance of nonsynonymous substitutions. We collected 23 haplogroups of Cytb with signals of late Quaternary population expansion events from rodents and moles and categorized them into three groups for comparison based on predicted expansion start time: 5,000-15,000 years ago (Group I), ca. 53,000 years ago (Group II) and 130,000-230,000 years ago (Group III). We counted the numbers of nonsynonymous and synonymous substitutions in all haplogroups. The rates of nonsynonymous substitutions were lowest in Groups II and III (0.08-0.22), whereas those in Group I varied markedly. We further classified Group I into two subgroups based on high (0.29-0.43) and low (0.09-0.20) nonsynonymous substitution rates, which were likely to be associated with the start of the expansion within 10,000 years and at around 15,000 years ago, respectively. The Group II and III networks had two- or three-step star-shaped structures and tended to exhibit frequent and less frequent nonsynonymous substitutions on exterior and interior branches, respectively. Based on temporal dynamics, nonsynonymous mitochondrial DNA (mtDNA) substitutions in small mammals accounted for at most 40% of all substitutions during the early evolutionary stage and then rapidly declined, dropping to approximately 15%. The results of this study provide a good explanation of the time-dependent trend in the mtDNA evolution rate predicted in previous work.

The genetic diversity of the genus Mus (Linnaeus, 1758) in the eastern part of the North Caucasus

Background

There are two species of Mus in the Caucasus: M. musculus and M. macedonicus. M. musculus is widespread in the Caucasus, where the species is found everywhere from the Black to the Caspian Sea. M. macedonicus is ubiquitous Transcaucasia. The most north-astern border of its distribution in the Caucasus, according to the literature, is located in the Derbent region, near the border between Dagestan and Azerbaijan.

Results

Cytochrome b mt-DNA of genus Mus research in this study in the Eastern Caucasus. About 70% of M. musculus haplotypes from the lowlands of Dagestan were recorded for the first time. One of these haplotypes accounts for approximately 25% of the total species diversity of haplotypes. M. macedonicus was found in only one locality, the Sarykum barchans, where this species prevails in number and accounts for 70% of the total number mice of the genus Mus. The species is characterized by low values of genetic diversity and nucleotide variability, which may indicate that the population originated from a small number of founders and may explain its relative isolation from the main range. The dating of the appearance of the ancestors of M. musculus in the east of the Russian Caucasus corresponds to 99-66 thousand years ago (at a mutation rate of 3-10% per million years).

Conclusion

The results obtained suggest that the history of the appearance of M. musculus in the Eastern Caucasus is more ancient and is not associated with human agricultural activities.

We believe that possibly the ancestral range of M. musculus covered the eastern and western coasts of the Caspian Sea in the territory of southern Dagestan, Azerbaijan, and Iran. In this paper M. macedonicus, a Balkan-Asia Minor species, was registered for the first time in the North Caucasus. This species was registered in the center of Dagestan, where it inhabits sympatrically (on the territory) and syntopically (on the same biotope) with M. musculus. The low values of genetic diversity of M. macedonicus in the North Caucasus suggest that the population originated from a small group of founders.

Genetic characterizations of Toll-like receptors in the brown rat and their associations with pathogen infections

Toll-like receptors (TLRs) are important initiators of innate immune responses that target host-pathogen interactions. However, further research into the molecular characteristics of TLRs in wild populations is required, as well as how TLRs genetically influenced pathogen infections in the brown rat (Rattus norvegicus). Here, we explored the genetic characterization and evolution of 2 sensing nucleic acid TLRs (TLR7 and TLR8) and 2 sensing non-nucleic acid TLRs (TLR2 and TLR4) in the wild brown rat, and assessed their associations with 2 RNA viruses (Seoul hantavirus and rat hepatitis E virus (HEV)) and 2 bacteria (Leptospira and Bartonella). In these 4 TLRs, we discovered a total of 16 variants. Furthermore, TLR8 had high genetic diversity among 7 variants, while TLR2 had low genetic diversity with only 1 variant. According to selective pressure analyses, TLR4, TLR7, and TLR8 genes evolved under purifying selection. Interestingly, significant associations were found between 3 TLR8 variants and HEV infection, as well as 1 TLR2 variant and Bartonella infection. Overall, our findings provided a glimpse into the genetic characterization of TLRs in the brown rat, and further demonstrated that TLR2 and TLR8 genetic variations were related to Bartonella and HEV infection, respectively. Especially, TLR8 may be a good candidate immune gene for future research on molecular ecology and functional adaptation in wild populations.

Positive selection and convergent evolution shape molecular phenotypic traits of innate immunity receptors in tits (Paridae)

Despite widespread variability and redundancy abounding animal immunity, little is currently known about the rate of evolutionary convergence (functionally analogous traits not inherited from a common ancestor) in host molecular adaptations to parasite selective pressures. Toll-like receptors (TLRs) provide the molecular interface allowing hosts to recognize pathogenic structures and trigger early danger signals initiating an immune response. Using a novel combination of bioinformatic approaches, here we explore genetic variation in ligand-binding regions of bacteria-sensing TLR4 and TLR5 in 29 species belonging to the tit family of passerine birds (Aves: Paridae). Three out of the four consensual positively selected sites in TLR4 and six out of 14 positively selected positions in TLR5 were located on the receptor surface near the functionally important sites, and based on the phylogenetic pattern evolved in a convergent (parallel) manner. This type of evolution was also seen at one N-glycosylation site and two positively selected phosphorylation sites, providing the first evidence of convergence in post-translational modifications in evolutionary immunology. Finally, the overall mismatch between phylogeny and the clustering of surface charge distribution demonstrates that convergence is common in overall TLR4 and TLR5 molecular phenotypes involved in ligand binding. Our analysis did not reveal any broad ecological traits explaining the convergence observed in electrostatic potentials, suggesting that information on microbial symbionts may be needed to explain TLR evolution. Adopting state-of-the-art predictive structural bionformatics, we have outlined a new broadly applicable methodological approach to estimate the functional significance of positively selected variation and test for the adaptive molecular convergence in protein-coding polymorphisms.

Pathogen-associated selection on innate immunity genes (TLR4, TLR7) in a neotropical rodent in landscapes differing in anthropogenic disturbance

Toll-like receptors (TLRs) form part of the innate immune system and can recognize structurally conserved pathogen-associated molecular pattern (PAMP) molecules. Their functional importance in the resistance to pathogens has been documented in laboratory experimental settings and in humans. TLR diversity, however, has been rarely investigated in wildlife species. How the genetic diversity of TLRs is associated with various pathogens and how it is shaped by habitat disturbance are understudied. Therefore, we investigated the role of genetic diversity in the functionally important parts of TLR4 and TLR7 genes in resistance towards gastrointestinal nematodes and Hepacivirus infection. We chose a generalist study species, the rodent Proechimys semispinosus, because it is highly abundant in three Panamanian landscapes that differ in their degree of anthropogenic modification. We detected only two TLR7 haplotypes that differed by one synonymous single-nucleotide polymorphism (SNP) position. The TLR4 variability was higher, and we detected four TLR4 haplotypes that differed at one synonymous SNP and at three amino acid positions within the leucine-rich repeat region. Only TLR4 haplotypes had different frequencies in each landscape. Using generalized linear models, we found evidence that nematode loads and virus prevalence were influenced by both specific TLR4 haplotypes and landscape. Here, the variable “landscape” served as a surrogate for the important influential ecological factors distinguishing landscapes in our study, i.e. species diversity and host population density. Individuals carrying the common TLR4_Ht1 haplotype were less intensely infected by the most abundant strongyle nematode. Individuals carrying the rare TLR4_Ht3 haplotype were all Hepacivirus-positive, where those carrying the rare haplotype TLR4_Ht4 were less often infected by Hepacivirus than individuals with other haplotypes. Our study highlights the role of TLR diversity in pathogen resistance and the importance of considering immune genetic as well as ecological factors in order to understand the effects of anthropogenic changes on wildlife health.

Intensity of infection with intracellular Eimeria spp. and pinworms is reduced in hybrid mice compared to parental subspecies

Genetic diversity in animal immune systems is usually beneficial. In hybrid recombinants, this is less clear, as the immune system could also be impacted by genetic conflicts. In the European house mouse hybrid zone, the long-standing impression that hybrid mice are more highly parasitized and less fit than parentals persists despite the findings of recent studies. Working across a novel transect, we assessed infections by intracellular protozoans, Eimeria spp., and infections by extracellular macroparasites, pinworms. For Eimeria, we found lower intensities in hybrid hosts than in parental mice but no evidence of lowered probability of infection or increased mortality in the centre of the hybrid zone. This means ecological factors are very unlikely to be responsible for the reduced load of infected hybrids. Focusing on parasite intensity (load in infected hosts), we also corroborated reduced pinworm loads reported for hybrid mice in previous studies. We conclude that intensity of diverse parasites, including the previously unstudied Eimeria, is reduced in hybrid mice compared to parental subspecies. We suggest caution in extrapolating this to differences in hybrid host fitness in the absence of, for example, evidence for a link between parasitemia and health.

Avian Toll-like receptor allelic diversity far exceeds human polymorphism: an insight from domestic chicken breeds

Immune genes show remarkable levels of adaptive variation shaped by pathogen-mediated selection. Compared to humans, however, population polymorphism in animals has been understudied. To provide an insight into immunogenetic diversity in birds, we sequenced complete protein-coding regions of all Toll-like receptor (TLR) genes with direct orthology between mammals and birds (TLR3, TLR4, TLR5 and TLR7) in 110 domestic chickens from 25 breeds and compared their variability with a corresponding human dataset. Chicken TLRs (chTLRs) exhibit on average nine-times higher nucleotide diversity than human TLRs (hTLRs). Increased potentially functional non-synonymous variability is found in chTLR ligand-binding ectodomains. While we identified seven sites in chTLRs under positive selection and found evidence for convergence between alleles, no selection or convergence was detected in hTLRs. Up to six-times more alleles were identified in fowl (70 chTLR4 alleles vs. 11 hTLR4 alleles). In chTLRs, high numbers of alleles are shared between the breeds and the allelic frequencies are more equal than in hTLRs. These differences may have an important impact on infectious disease resistance and host-parasite co-evolution. Though adaptation through high genetic variation is typical for acquired immunity (e.g. MHC), our results show striking levels of intraspecific polymorphism also in poultry innate immune receptors.

Signatures of diversifying selection and convergence acting on passerine Toll-like receptor 4 in an evolutionary context

Positive selection acting on Toll-like receptors (TLRs) has been recently investigated to reveal evolutionary mechanisms of host-pathogen molecular co-adaptation. Much of this research, however, has focused mainly on the identification of sites predicted to be under positive selection, bringing little insight into the functional differences and similarities among species and a limited understanding of convergent evolution in the innate immune molecules. In this study, we provide evidence of phenotypic variability in the avian TLR4 ligand-binding region (LBR), the direct interface between host and pathogen molecular structures. We show that 55 passerine species vary substantially in the distribution of electrostatic potential on the surface of the receptor, and based on these distinct patterns, we identified four species clusters. Seven of the 34 evolutionarily nonconservative and positively selected residues correspond topologically to sites previously identified as being important for lipopolysaccharide, lipid IVa or MD-2 binding. Five of these positions codetermine the identity of the charge clusters. Groups of species that host-related communities of pathogens were predicted to cluster based on their TLR4 LBR charge. Despite some evidence for convergence among taxa, there were no clear associations between the TLR4 LBR charge distribution and any of the general ecological characteristics compared (migration, latitudinal distribution and diet). Closely related species, however, mostly belonged to the same surface charge cluster indicating that phylogenetic constraints are key determinants shaping TLR4 adaptive evolution. Our results suggest that host innate immune evolution is consistent with Fahrenholz's rule on the cospeciation of hosts and their parasites.

Characterization of the peripheral blood transcriptome and adaptive evolution of the MHC I and TLR gene families in the wolf (Canis lupus)

Background

The wolf (Canis lupus) is one of the most widely distributed terrestrial mammals, because it is well adapted to various ecological niches and their corresponding pathogen environments. Immunological competence is a crucial factor involved in adapting to a changing environment and fighting pathogen infection in animals. In this study, the peripheral blood transcriptome of wolves was generated via RNA-seq to advance understanding of the wolf immunome, with a special focus on the major histocompatibility complex class I (MHC I) and toll-like receptor (TLR) gene families, which are involved in pathogen recognition and defense.

Results

The blood transcriptomic libraries of eight wolves originating from Tibet and Inner Mongolia were sequenced, and approximately 383 million reads were generated. Using a genome-guided assembly strategy, we obtained 123,851 unigenes, with a mean length of 845 bp and an N50 length of 1121 bp. On the basis of BLAST searches against the NCBI non-redundant protein database (Nr), a total of 36,192 (29.22%) unigenes were annotated. For functional classification, 24,663 unigenes were assigned to 13,016 Gene Ontology (GO) terms belonging to 51 sub-categories of the three main GO categories. Additionally, 7682 unigenes were classified into 6 Kyoto Encyclopedia of Genes and Genomes (KEGG) categories, in which the most represented functional sub-categories were signal transduction and the immune system, and 16,238 unigenes were functionally classified into 25 Eukaryotic Orthologous Groups (KOG) categories. We observed an overall higher ω (d_N/d_S) value at antigen-binding sites (ABSs) than at non-ABS regions as well as clear evidence of intergenic/intragenic recombination events at wolf MHC I loci. Additionally, our analysis revealed that carnivorous TLRs were dominated by purifying selection, with mean ω values at each TLR locus ranging from 0.173 to 0.527. However, we also found significant instances of positive selection that acted on several codons in pathogen recognition domains and were linked to species-specific differences in pathogen recognition.

Conclusions

This study represents the first attempt to characterize the blood transcriptome of the wolf and to highlight the value of investigating the immune system. Balancing selection and recombination have contributed to the historical evolution of wolf MHC I genes. Moreover, TLRs in carnivores have undergone adaptive evolution against the background of purifying selection, and a high level of adaptive evolution was detected in the wolf TLR system.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3983-0) contains supplementary material, which is available to authorized users.

Microevolution of bank voles (Myodes glareolus) at neutral and immune-related genes during multiannual dynamic cycles: Consequences for Puumala hantavirus epidemiology

Understanding how host dynamics, including variations of population size and dispersal, may affect the epidemiology of infectious diseases through ecological and evolutionary processes is an active research area. Here we focus on a bank vole (Myodes glareolus) metapopulation surveyed in Finland between 2005 and 2009. Bank vole is the reservoir of Puumala hantavirus (PUUV), the agent of nephropathia epidemica (NE, a mild form of hemorrhagic fever with renal symptom) in humans. M. glareolus populations experience multiannual density fluctuations that may influence the level of genetic diversity maintained in bank voles, PUUV prevalence and NE occurrence. We examine bank vole metapopulation genetics at presumably neutral markers and immune-related genes involved in susceptibility to PUUV (Tnf-promoter, Tlr4, Tlr7 and Mx2 gene) to investigate the links between population dynamics, microevolutionary processes and PUUV epidemiology. We show that genetic drift slightly and transiently affects neutral and adaptive genetic variability within the metapopulation. Gene flow seems to counterbalance its effects during the multiannual density fluctuations. The low abundance phase may therefore be too short to impact genetic variation in the host, and consequently viral genetic diversity. Environmental heterogeneity does not seem to affect vole gene flow, which might explain the absence of spatial structure previously detected in PUUV in this area. Besides, our results suggest the role of vole dispersal on PUUV circulation through sex-specific and density-dependent movements. We find little evidence of selection acting on immune-related genes within this metapopulation. Footprint of positive selection is detected at Tlr-4 gene in 2008 only. We observe marginally significant associations between Mx2 genotype and PUUV genogroups. These results show that neutral processes seem to be the main factors affecting the evolution of these immune-related genes at a contemporary scale, although the relative effects of neutral and adaptive forces could vary temporally with density fluctuations. Immune related gene polymorphism may in turn partly influence PUUV epidemiology in this metapopulation.

Polymorphisms of the Toll-Like Receptor 2 of Goats (Capra hircus) may be Associated with Somatic Cell Count in Milk

Toll-like receptor 2 (TLR2) plays an important role in recognition by the innate immune system of Gram-positive bacteria. As Gram-positive bacteria cause mastitis, we examined variations in the region of the TLR2 gene that codes for the extracellular domain. Samples of forty goats from a single dairy herd were collected, half with low SCC (≤200,000 cells/mL), and half with higher SCC. Two synonymous single nucleotide polymorphisms (SNPs) were identified: 840G > A and 1083A > G. One nonsynonymous SNP 739G > A was identified. This coded for valine or isoleucine, which have similar physiochemical properties, and was not in a region coding for a known functional domain. Surprisingly, the least square mean SCC of the heterozygous goats (146,220) was significantly lower than the SCC of homozygous GG goats (537,700; p = 0.004), although these two groups were similar in days in milk (p = 0.984), and there was no significant difference by breed (p = 0.941). Because factors other than mastitis can affect SCC and our sample sizes were limited, additional studies are needed to corroborate an association between TLR2 genotype and SCC or mastitis in goats.

Trans-Species Polymorphism in Immune Genes: General Pattern or MHC-Restricted Phenomenon?

Immunity exhibits extraordinarily high levels of variation. Evolution of the immune system in response to host-pathogen interactions in particular ecological contexts appears to be frequently associated with diversifying selection increasing the genetic variability. Many studies have documented that immunologically relevant polymorphism observed today may be tens of millions years old and may predate the emergence of present species. This pattern can be explained by the concept of trans-species polymorphism (TSP) predicting the maintenance and sharing of favourable functionally important alleles of immune-related genes between species due to ongoing balancing selection. Despite the generality of this concept explaining the long-lasting adaptive variation inherited from ancestors, current research in TSP has vastly focused only on major histocompatibility complex (MHC). In this review we summarise the evidence available on TSP in human and animal immune genes to reveal that TSP is not a MHC-specific evolutionary pattern. Further research should clearly pay more attention to the investigation of TSP in innate immune genes and especially pattern recognition receptors which are promising candidates for this type of evolution. More effort should also be made to distinguish TSP from convergent evolution and adaptive introgression. Identification of balanced TSP variants may represent an accurate approach in evolutionary medicine to recognise disease-resistance alleles.

Adaptive Evolution of Toll-Like Receptors (TLRs) in the Family Suidae

Members of the family Suidae have diverged over extended evolutionary periods in diverse environments, suggesting that adaptation in response to endemic infectious agents may have occurred. Toll-like receptors (TLRs) comprise a multigene family that acts as the first line of defense against infectious microbes at the host-environment interface. We hypothesized that across the Suidae, positive selection mediated by infectious agents has contributed to the evolution of TLR diversity. Thus, we analyzed Sus scrofa, Sus barbatus, Sus verrucosus, Sus celebensis, Sus scebifrons, Babyrousa babyrussa, Potamochoerus larvatus, Potamochoerus porcus and Phacochoerus africanus genomes. Specifically, analyses were performed to identify evidence of positive selection using Maximum likelihood (ML) methods within a phylogenetic framework for bacterial and viral sensing Suidae TLR extracellular domains. Our analyses did not reveal evidence of positive selection for TLR3 and TLR7, suggesting strong functional conservation among these two genes for members of the Suidae. Positive selection was inferred for Suidae TLR1, TLR2, TLR6 and TLR8 evolution. ML methods identified amino acid sites of the bacterial sensing TLR1, TLR2, TLR6 and the viral sensing TLR8 to be under persistent positive selection. Some of these sites are in close proximity to functionally relevant sites, further strengthening the case for pathogen mediated selection for these sites. The branch leading to the genus Sus demonstrated evidence of episodic positive selection for TLR1, indicating selection mediated by infectious agents encountered within the specific geographic origin of the Sus. These results indicate that species of the Suidae have positively selected residues within functional domains of TLRs reflective of prior infections. Thus, TLR genes represent candidates for experimental validation to determine their functional role in antibacterial and antiviral activity within members of the Suidae.