Contrasting patterns of diversity and population differentiation at the innate immunity gene toll-like receptor 2 (TLR2) in two sympatric rodent species

Comparing patterns of diversity and divergence between populations at immune genes and neutral markers can give insights into the nature and geographic scale of parasite-mediated selection. To date, studies investigating such patterns of selection in vertebrates have primarily focused on the acquired branch of the immune system, whereas it remains largely unknown how parasite-mediated selection shapes innate immune genes both within and across vertebrate populations. Here, we present a study on the diversity and population differentiation at the innate immune gene Toll-like receptor 2 (TLR2) across nine populations of yellow-necked mice (Apodemus flavicollis) and bank voles (Myodes glareolus) in southern Sweden. In yellow-necked mice, TLR2 diversity was very low, as was TLR2 population differentiation compared to neutral loci. In contrast, several TLR2 haplotypes co-occurred at intermediate frequencies within and across bank vole populations, and pronounced isolation by distance between populations was observed. The diversity and differentiation at neutral loci was similar in the two species. These results indicate that parasite-mediated selection has been acting in dramatically different ways on a given immune gene in ecologically similar and sympatric species. Furthermore, the finding of TLR2 population differentiation at a small geographical scale in bank voles highlights that vertebrate innate immune defense may be evolutionarily more dynamic than has previously been appreciated.     

Balancing selection is the main force shaping the evolution of innate immunity genes

The evolutionarily recent geographic expansion of humans, and the even more recent development of large, relatively dense human settlements, has exposed our species to new pathogenic environments. Potentially lethal pathogens are likely to have exerted important selective pressures on our genome, so immunity genes can be expected to show molecular signatures of the adaptation of human populations to these recent conditions. While genes related to the acquired immunity system have indeed been reported to show traces of local adaptation, little is known about the response of the innate immunity system. In this study, we analyze the variability patterns in different human populations of fifteen genes related to innate immunity. We have used both single nucleotide polymorphism and sequence data, and through the analysis of interpopulation differentiation, the linkage disequilibrium pattern, and intrapopulation diversity, we have discovered some signatures of positive and especially balancing selection in these genes, thus confirming the importance of the immune system genetic plasticity in the evolutionary adaptive process. Interestingly, the strongest evidence is found in three TLR genes and CD14. These innate immunity genes play a pivotal role, being involved in the primary recognition of pathogens. In general, more evidences of selection appear in the European populations, in some case possibly related to severe population specific pressures. However, we also describe evidence from African populations, which may reflect parallel or long-term selective forces acting in different geographic areas.     

Molecular evolution of the toll-like receptor multigene family in birds

Toll-like receptors (TLR) are membrane-bound sensors of the innate immune system that recognize invariant and distinctive molecular features of invading microbes and are also essential for initiating adaptive immunity in vertebrates. The genetic variation at TLR genes has been directly related to differential pathogen outcomes in humans and livestock. Nonetheless, new insights about the impact of TLRs polymorphism on the evolutionary ecology of infectious diseases can be gained through the investigation of additional vertebrate groups not yet investigated in detail. In this study, we have conducted the first characterization of the entire TLR multigene family (N = 10 genes) in non-model avian species. Using primers targeting conserved coding regions, we aimed at amplifying large segments of the extracellular domains (275-435 aa) involved in pathogen recognition across seven phylogenetically diverse bird species. Our analyses suggest avian TLRs are dominated by stabilizing selection, suggesting that slow rates of nonsynonymous substitution help preserve biological function. Overall, mean values of ω (= d(n)/d(s)) at each TLR locus ranged from 0.196 to 0.517. However, we also found patterns of positive selection acting on specific amino acid sites that could be linked to species-specific differences in pathogen-associated molecular pattern recognition. Only 39 of 2,875 (～1.35%) of the codons analyzed exhibited significant patterns of positive selection. At least one half of these positively selected codons can be mapped to putative ligand-binding regions, as suggested by crystallographic structures of TLRs and their ligands and mutagenic analyses. We also surveyed TLR polymorphism in wild populations of two bird species, the Lesser Kestrel Falco naumanni and the House Finch Carpodacus mexicanus. In general, avian TLRs displayed low to moderate single nucleotide polymorphism levels and an excess of silent nucleotide substitutions, but also conspicuous instances of positive directional selection. In particular, TLR5 and TLR15 exhibited the highest degree of genetic polymorphism and the highest occurrence of nonconservative amino acid substitutions. This study provides critical primers and a first look at the evolutionary patterns and implications of TLR polymorphism in non-model avian species and extends the list of candidate loci for avian eco-immunogenetics beyond the widely employed genes of the Major Histocompatibility Complex (MHC).     

Positive selection signatures in the TLR7 family

While adaptive immunity genes evolve rapidly under the influence of positive selection, innate immune system genes are known to evolve slowly due to strong purifying selection. Among the sensors of the innate immune system, Toll-like receptors (TLRs) are particularly important due to their ability to recognize and respond to pathogen-associated molecular patterns (PAMP), such as lipopolysaccharides, peptidoglycans, and nucleic acids from bacteria or viruses. In the present study, we examine the evolutionary process that has operated on the TLR7 family genes TLR7, TLR8, and TLR9. The results demonstrate that the average Ka/Ks (the ratio between nonsynonymous and synonymous substitution rates) of each TLR family gene is far lower than one regardless of estimating methods, supporting previous observations of strong purifying selection in this gene family. Interestingly, however, analysis of Ka/Ks ratios along the coding regions of TLR7 family genes by sliding-window analysis reveals a few narrow high peaks (Ka/Ks > 1). The most prominent peak corresponds to a specific region in the ectodomain, which exists only in the TLR7 family, suggesting that this unique structure of the TLR7 family might have been a target of positive selection in a variety of lineages. Furthermore, maximum likelihood model tests suggest that positive selection is the best explanation for a certain fraction of the amino acid substitutions in the TLR9.     

The immune gene repertoire encoded in the purple sea urchin genome

Echinoderms occupy a critical and largely unexplored phylogenetic vantage point from which to infer both the early evolution of bilaterian immunity and the underpinnings of the vertebrate adaptive immune system. Here we present an initial survey of the purple sea urchin genome for genes associated with immunity. An elaborate repertoire of potential immune receptors, regulators and effectors is present, including unprecedented expansions of innate pathogen recognition genes. These include a diverse array of 222 Toll-like receptor (TLR) genes and a coordinate expansion of directly associated signaling adaptors. Notably, a subset of sea urchin TLR genes encodes receptors with structural characteristics previously identified only in protostomes. A similarly expanded set of 203 NOD/NALP-like cytoplasmic recognition proteins is present. These genes have previously been identified only in vertebrates where they are represented in much lower numbers. Genes that mediate the alternative and lectin complement pathways are described, while gene homologues of the terminal pathway are not present. We have also identified several homologues of genes that function in jawed vertebrate adaptive immunity. The most striking of these is a gene cluster with similarity to the jawed vertebrate Recombination Activating Genes 1 and 2 (RAG1/2). Sea urchins are long-lived, complex organisms and these findings reveal an innate immune system of unprecedented complexity. Whether the presumably intense selective processes that molded these gene families also gave rise to novel immune mechanisms akin to adaptive systems remains to be seen. The genome sequence provides immediate opportunities to apply the advantages of the sea urchin model toward problems in developmental and evolutionary immunobiology.     

Polymorphisms at the innate immune receptor TLR2 are associated with Borrelia infection in a wild rodent population

The discovery of the key role of Toll-like receptors (TLRs) in initiating innate immune responses and modulating adaptive immunity has revolutionized our understanding of vertebrate defence against pathogens. Yet, despite their central role in pathogen recognition and defence initiation, there is little information on how variation in TLRs influences disease susceptibility in natural populations. Here, we assessed the extent of naturally occurring polymorphisms at TLR2 in wild bank voles (Myodes glareolus) and tested for associations between TLR2 variants and infection with Borrelia afzelii, a common tick-transmitted pathogen in rodents and one of the causative agents of human Lyme disease. Bank voles in our population had 15 different TLR2 haplotypes (10 different haplotypes at the amino acid level), which grouped in three well-separated clusters. In a large-scale capture–mark–recapture study, we show that voles carrying TLR2 haplotypes of one particular cluster (TLR2_c2) were almost three times less likely to be Borrelia infected than animals carrying other haplotypes. Moreover, neutrality tests suggested that TLR2 has been under positive selection. This is, to our knowledge, the first demonstration of an association between TLR polymorphism and parasitism in wildlife, and a striking example that genetic variation at innate immune receptors can have a large impact on host resistance.     

Small‐scale spatial variation in infection risk shapes the evolution of a Borrelia resistance gene in wild rodents

Spatial variation in pathogen‐mediated selection is predicted to influence the evolutionary trajectory of host populations and lead to spatial variation in their immunogenetic composition. However, to date few studies have been able to directly link small‐scale spatial variation in infection risk to host immune gene evolution in natural, nonhuman populations. Here, we use a natural rodent–Borrelia system to test for associations between landscape‐level spatial variation in Borrelia infection risk along replicated elevational gradients in the Swiss Alps and Toll‐like receptor 2 (TLR2) evolution, a candidate gene for Borrelia resistance, across bank vole (Myodes glareolus) populations. We found that Borrelia infection risk (i.e., the product of Borrelia prevalence in questing ticks and the average tick load of voles at a sampling site) was spatially variable and significantly negatively associated with elevation. Across sampling sites, Borrelia prevalence in bank voles was significantly positively associated with Borrelia infection risk along the elevational clines. We observed a significant association between naturally occurring TLR2 polymorphisms in hosts and their Borrelia infection status. The TLR2 variant associated with a reduced likelihood of Borrelia infection was most common in rodent populations at lower elevations that face a high Borrelia infection risk, and its frequency changed in accordance with the change in Borrelia infection risk along the elevational clines. These results suggest that small‐scale spatial variation in parasite‐mediated selection affects the immunogenetic composition of natural host populations, providing a striking example that the microbial environment shapes the evolution of the host's immune system in the wild.     

Distinct haplotype structure at the innate immune receptor Toll‐like receptor 2 across bank vole populations and lineages in Europe

Parasite‐mediated selection may contribute to the maintenance of genetic variation at host immune genes over long time scales. To date, the best evidence for the long‐term maintenance of immunogenetic variation in natural populations comes from studies on the major histocompatibility complex (MHC) genes, whereas evidence for such processes from other immune genes remains scarce. In the present study, we show that, despite pronounced population differentiation and the occurrence of numerous private alleles within populations, the innate immune gene Toll‐like receptor 2 (TLR2) displays a distinct haplotype structure in 21 bank vole (Myodes glareolus) populations across Europe. Haplotypes from all populations grouped in four clearly differentiated clusters, with the three main clusters co‐occurring in at least three previously described mitochondrial lineages. This pattern indicates that the distinct TLR2 haplotype structure may precede the split of the mitochondrial lineages 0.19–0.56 Mya and suggests that haplotype clusters at this innate immune receptor are maintained over prolonged time in wild bank vole populations.     

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.     

The evolutionary characteristics and structural biology of Gallus toll‐like receptor 21

Toll‐like receptors (TLRs) are an important part of the innate immune system, acting as a first line of defense against many invading pathogens. The ligand known to bind Gallus toll‐like receptor 21 (gTLR21) is the unmethylated cytosine phosphate guanine dideoxy nucleotide motif; however, the evolutionary characteristics and structural biology of gTLR21 are poorly elaborated. Our results suggest that gTLR21 is phylogenetically and evolutionarily related to the TLR11 family and is perhaps a close ortholog of the Mus TLR13. Structural biology of homology modeling of the gTLR21 ectodomain structure suggests that it has no Z‐loop like that seen in Mus TLR9. The cytosolic toll‐IL‐1 receptor region of gTLR21 contains a central 4‐stranded parallel β‐sheet (βA‐βD) surrounded by 5 α‐helices (αA‐αE) on both sides, a highly conserved structure also seen in other TLRs. Molecular docking analysis reveals that the gTLR21 ectodomain has the potential to distinguish between different ligands. Homodimer analysis results also suggest that Phe842 and Pro844 of the BB loop and Cys876 of the αC helix in gTLR21 are conserved in other cytosolic toll‐IL‐1 receptor domains of other TLRs and may contribute to the docking of homodimers. Our study on the evolutionary characteristics and structural biology of gTLR21 reveals that the molecule may have a broader role to play in innate immune system; however, further experimental validation is required to confirm our findings.     

Evolution of the bovine TLR gene family and member associations with Mycobacterium avium subspecies paratuberculosis infection

Members of the Toll-like receptor (TLR) gene family occupy key roles in the mammalian innate immune system by functioning as sentries for the detection of invading pathogens, thereafter provoking host innate immune responses. We utilized a custom next-generation sequencing approach and allele-specific genotyping assays to detect and validate 280 biallelic variants across all 10 bovine TLR genes, including 71 nonsynonymous single nucleotide polymorphisms (SNPs) and one putative nonsense SNP. Bayesian haplotype reconstructions and median joining networks revealed haplotype sharing between Bos taurus taurus and Bos taurus indicus breeds at every locus, and specialized beef and dairy breeds could not be differentiated despite an average polymorphism density of 1 marker/158 bp. Collectively, 160 tagSNPs and two tag insertion-deletion mutations (indels) were sufficient to predict 100% of the variation at 280 variable sites for both Bos subspecies and their hybrids, whereas 118 tagSNPs and 1 tagIndel predictively captured 100% of the variation at 235 variable sites for B. t. taurus. Polyphen and SIFT analyses of amino acid (AA) replacements encoded by bovine TLR SNPs indicated that up to 32% of the AA substitutions were expected to impact protein function. Classical and newly developed tests of diversity provide strong support for balancing selection operating on TLR3 and TLR8, and purifying selection acting on TLR10. An investigation of the persistence and continuity of linkage disequilibrium (r2≥0.50) between adjacent variable sites also supported the presence of selection acting on TLR3 and TLR8. A case-control study employing validated variants from bovine TLR genes recognizing bacterial ligands revealed six SNPs potentially eliciting small effects on susceptibility to Mycobacterium avium spp paratuberculosis infection in dairy cattle. The results of this study will broadly impact domestic cattle research by providing the necessary foundation to explore several avenues of bovine translational genomics, and the potential for marker-assisted vaccination.     

Tracing the evolution of novel features of human Toll‐like receptor 4

Toll‐like receptor 4 (TLR4) is a critical innate immune protein that activates inflammation in response to extracellular cues. Much of the work to understand how the protein works in humans has been done using mouse models. Although human and mouse TLR4 have many shared features, they have also diverged significantly since their last common ancestor, acquiring 277 sequence differences. Functional differences include the extent of ligand‐independent activation, whether lipid IVa acts as an antagonist or agonist, and the relative species cross‐compatibility of their MD‐2 cofactor. We set out to understand the evolutionary origins for these functional differences between human and mouse TLR4. Using a combination of phylogenetics, ancestral sequence reconstruction, and functional characterization, we found that evolutionary changes to the human TLR4, rather than changes to the mouse TLR4, were largely responsible for these functional changes. Human TLR4 repressed ancestral ligand‐independent activity and gained antagonism to lipid IVa. Additionally, mutations to the human TLR4 cofactor MD‐2 led to lineage‐specific incompatibility between human and opossum TLR4 complex members. These results were surprising, as mouse TLR4 has acquired many more mutations than human TLR4 since their last common ancestor. Our work has polarized this set of transitions and sets up work to study the mechanistic underpinnings for the evolution of new functions in TLR4.     

Genetic mechanisms of adaptive immunity emergence in vertebrates

The adaptive immune system in vertebrates emerged in a multistep process that can be reconstructed on the basis of the data concerning the structure of immune systems of modern cartilaginous and bony fishes, as well as of cyclostomes. The most probable evolutionary scenario is likely to be as follows: the T cell receptor loci emerged on the basis of NK cell-like receptor genes; the antibody loci evolved on the basis of T cell receptor loci; the MHC locus arose on the basis of the locus responsible for innate immunity of early chordates. The ancestral MHC molecules likely participated in the transplantation immunity before they acquired the ability of antigen peptide presentation.     

Toll-dependent and toll-independent innate antiviral immunity

Until recently, adaptive immunity and cytotoxic T cells were considered as the only essential components of the antiviral defence arsenal. Additional data that do not rule out the crucial role of these cells in the clearance of viral pathogens have, however, recently emerged. They indicate that innate immune cells such as macrophages, dendritic cells, gammadelta T cells as well as natural killer (NK) cells play a primordial role in this mechanism. It is now well established that innate immune cells can detect various pathogens (bacteria, viruses, fungi or parasites) very rapidly and respond to their presence through the activation of specific receptors. Once activated, these molecules trigger several signalling cascades that culminate in the establishment of very potent defence mechanisms. In addition, cytokines produced during this initial response are essential for the activation of the adaptive immune response which will add specificity and memory to the system. Among the innate immune receptors, attention has focused on the Toll-like receptors (TLR) and many reports indicate that some of the TLRs are clearly involved in defence against viral pathogens. However, new molecules, acting independently from any TLR, have recently been discovered. They define a second antiviral pathway which is presently the subject of intense research. In this article, we will review the role of the different molecules involved in each pathway within the framework of innate antiviral defence.     

Comparative fine maps of bovine toll-like receptor 4 and toll-like receptor 2 regions

Abstract Toll-like receptors are cell-surface receptors that activate innate and adaptive immune responses. We have used a 5000-rad, whole-genome radiation hybrid panel to map Toll-like receptor 4 (TLR4) to the distal end of bovine Chromosome (Chr) 8, and Toll-like receptor 2 (TLR2) to the proximal end of bovine Chr 17. To facilitate comparative mapping and contig construction, we have also used 5000- and 12,000-rad, whole-genome radiation hybrid panels to produce fine maps of the regions surrounding these genes in cattle. These fine maps triple the number of available markers in the TLR4 region and more than double the number of available markers in the TLR2 region. Comparative analyses show gene order conservation between the bovine Chr 8 region and human Chr 9, and between the bovine Chr 17 region and human Chr 4. In addition, the bovine Chr 8 region refines an evolutionary chromosomal breakpoint from a 10-megabase region to a 2.5-megabase region, and the bovine Chr 17 map suggests a new evolutionary chromosomal breakpoint.     

Drift,not selection,shapes toll‐like receptor variation among oceanic island populations

Understanding the relative role of different evolutionary forces in shaping the level and distribution of functional genetic diversity among natural populations is a key issue in evolutionary and conservation biology. To do so accurately genetic data must be analysed in conjunction with an unambiguous understanding of the historical processes that have acted upon the populations. Here, we focused on diversity at toll‐like receptor (TLR) loci, which play a key role in the vertebrate innate immune system and, therefore, are expected to be under pathogen‐mediated selection. We assessed TLR variation within and among 13 island populations (grouped into three archipelagos) of Berthelot's pipit, Anthus berthelotii, for which detailed population history has previously been ascertained. We also compared the variation observed with that found in its widespread sister species, the tawny pipit, Anthus campestris. We found strong evidence for positive selection at specific codons in TLR1LA, TLR3 and TLR4. Despite this, we found that at the allele frequency level, demographic history has played the major role in shaping patterns of TLR variation in Berthelot's pipit. Levels of diversity and differentiation within and across archipelagos at all TLR loci corresponded very closely with neutral microsatellite variation and with the severity of the bottlenecks that occurred during colonization. Our study shows that despite the importance of TLRs in combating pathogens, demography can be the main driver of immune gene variation within and across populations, resulting in patterns of functional variation that can persist over evolutionary timescales.     

Toll‐like receptor variation in the bottlenecked population of the Seychelles warbler: computer simulations see the ‘ghost of selection past’ and quantify the ‘drift debt’

Balancing selection can maintain immunogenetic variation within host populations, but detecting its signal in a postbottlenecked population is challenging due to the potentially overriding effects of drift. Toll‐like receptor genes (TLRs) play a fundamental role in vertebrate immune defence and are predicted to be under balancing selection. We previously characterized variation at TLR loci in the Seychelles warbler (Acrocephalus sechellensis), an endemic passerine that has undergone a historical bottleneck. Five of seven TLR loci were polymorphic, which is in sharp contrast to the low genomewide variation observed. However, standard population genetic statistical methods failed to detect a contemporary signature of selection at any TLR locus. We examined whether the observed TLR polymorphism could be explained by neutral evolution, simulating the population's demography in the software DIYABC. This showed that the posterior distributions of mutation rates had to be unrealistically high to explain the observed genetic variation. We then conducted simulations with an agent‐based model using typical values for the mutation rate, which indicated that weak balancing selection has acted on the three TLR genes. The model was able to detect evidence of past selection elevating TLR polymorphism in the prebottleneck populations, but was unable to discern any effects of balancing selection in the contemporary population. Our results show drift is the overriding evolutionary force that has shaped TLR variation in the contemporary Seychelles warbler population, and the observed TLR polymorphisms might be merely the ‘ghost of selection past’. Forecast models predict immunogenetic variation in this species will continue to be eroded in the absence of contemporary balancing selection. Such ‘drift debt’ occurs when a gene pool has not yet reached its new equilibrium level of polymorphism, and this loss could be an important threat to many recently bottlenecked populations.     

Phylogenetic and expression analysis of amphibian Xenopus Toll-like receptors

An anuran amphibian, South African clawed frog (Xenopus laevis), is used to study the immune system, as it possesses a set of acquired immune system represented by T and B lymphocytes and the immunoglobulins. The acquired immune system is impaired throughout the larva and the metamorphosis stage in the amphibians. On the other hand, the role of innate immune system in the tadpole remains unclear. Recently, insect Toll protein homologues, namely, Toll-like receptors (TLRs), have been identified as sensors recognizing microbe-pattern molecules in vertebrates. Whole-genome analysis of Xenopus tropicalis supported the existence of the tlr genes in the frog. In this study, we annotated 20 frog tlr gene nucleotide sequences from the latest genome assembly version 4.1 on the basis of homology and identified cDNAs of the predicted frog TLR proteins. Phylogenetic analysis showed that the repertoire of the frog TLRs consisted of both fish- and mammalian-type TLRs. We showed that the frog TLRs are constitutively expressed in the tadpole as well as in the adult frog. Our results suggest that tadpoles are protected from microbes by the innate system that includes TLRs, despite impaired acquired immune system in tadpoles. This is the first report on the properties of TLRs in the most primitive terrestrial animals like amphibia. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.