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

Genomic landscapes of divergence among island bird populations: Evidence of parallel adaptation but at different loci?

When populations colonise new environments, they may be exposed to novel selection pressures but also suffer from extensive genetic drift due to founder effects, small population sizes and limited interpopulation gene flow. Genomic approaches enable us to study how these factors drive divergence, and disentangle neutral effects from differentiation at specific loci due to selection. Here, we investigate patterns of genetic diversity and divergence using whole-genome resequencing (>22× coverage) in Berthelot's pipit (Anthus berthelotii), a passerine endemic to the islands of three north Atlantic archipelagos. Strong environmental gradients, including in pathogen pressure, across populations in the species range, make it an excellent system in which to explore traits important in adaptation and/or incipient speciation. First, we quantify how genomic divergence accumulates across the speciation continuum, that is, among Berthelot's pipit populations, between sub species across archipelagos, and between Berthelot's pipit and its mainland ancestor, the tawny pipit (Anthus campestris). Across these colonisation timeframes (2.1 million-ca. 8000 years ago), we identify highly differentiated loci within genomic islands of divergence and conclude that the observed distributions align with expectations for non-neutral divergence. Characteristic signatures of selection are identified in loci associated with craniofacial/bone and eye development, metabolism and immune response between population comparisons. Interestingly, we find limited evidence for repeated divergence of the same loci across the colonisation range but do identify different loci putatively associated with the same biological traits in different populations, likely due to parallel adaptation. Incipient speciation across these island populations, in which founder effects and selective pressures are strong, may therefore be repeatedly associated with morphology, metabolism and immune defence.

Genotype-environment associations reveal genes potentially linked to avian malaria infection in populations of an endemic island bird

Patterns of pathogen prevalence are, at least partially, the result of coevolutionary host-pathogen interactions. Thus, exploring the distribution of host genetic variation in relation to infection by a pathogen within and across populations can provide important insights into mechanisms of host defence and adaptation. Here, we use a landscape genomics approach (Bayenv) in conjunction with genome-wide data (ddRADseq) to test for associations between avian malaria (Plasmodium) prevalence and host genetic variation across 13 populations of the island endemic Berthelot's pipit (Anthus berthelotii). Considerable and consistent spatial heterogeneity in malaria prevalence was observed among populations over a period of 15 years. The prevalence of malaria infection was also strongly positively correlated with pox (Avipoxvirus) prevalence. Multiple host loci showed significant associations with malaria prevalence after controlling for genome-wide neutral genetic structure. These sites were located near to or within genes linked to metabolism, stress response, transcriptional regulation, complement activity and the inflammatory response, many previously implicated in vertebrate responses to malarial infection. Our findings identify diverse genes - not just limited to the immune system - that may be involved in host protection against malaria and suggest that spatially variable pathogen pressure may be an important evolutionary driver of genetic divergence among wild animal populations, such as Berthelot's pipit. Furthermore, our data indicate that spatio-temporal variation in multiple different pathogens (e.g. malaria and pox in this case) may have to be studied together to develop a more holistic understanding of host pathogen-mediated evolution.

Non-invasive sampling reveals low mitochondrial genetic diversity for an island endemic species: The critically endangered Grenada Dove Leptotila wellsi

As an island endemic with a decreasing population, the critically endangered Grenada Dove Leptotila wellsi is threatened by accelerated loss of genetic diversity resulting from ongoing habitat fragmentation. Small, threatened populations are difficult to sample directly but advances in molecular methods mean that non-invasive samples can be used. We performed the first assessment of genetic diversity of populations of Grenada Dove by (a) assessing mtDNA genetic diversity in the only two areas of occupancy on Grenada, (b) defining the number of haplotypes present at each site and (c) evaluating evidence of isolation between sites. We used non-invasively collected samples from two locations: Mt Hartman (n = 18) and Perseverance (n = 12). DNA extraction and PCR were used to amplify 1751 bps of mtDNA from two mitochondrial markers: NADH dehydrogenase 2 (ND2) and Cytochrome b (Cyt b). Haplotype diversity (h) of 0.4, a nucleotide diversity (π) of 0.00023 and two unique haplotypes were identified within the ND2 sequences; a single haplotype was identified within the Cyt b sequences. Of the two haplotypes identified, the most common haplotype (haplotype A = 73.9%) was observed at both sites and the other (haplotype B = 26.1%) was unique to Perseverance. Our results show low mitochondrial genetic diversity and clear evidence for genetically isolated populations. The Grenada Dove needs urgent conservation action, including habitat protection and potentially augmentation of gene flow by translocation in order to increase genetic resilience and diversity with the ultimate aim of securing the long-term survival of this critically endangered species.

Genetic variation at innate and adaptive immune genes - contrasting patterns of differentiation and local adaptation in a wild gull

Immunogenetic variation in natural vertebrate populations is expected to respond to spatial and temporal fluctuations in pathogen assemblages. While spatial heterogeneity in pathogen-driven selection enhances local immunogenetic adaptations and population divergence, different immune genes may yield contrasting responses to the environment. Here, we investigated population differentiation at the key pathogen recognition genes of the innate and adaptive immune system in a colonial bird species, the black-headed gull Chroicocephalus ridibundus. We assessed genetic variation at three toll-like receptor (TLR) genes (innate immunity) and the major histocompatibility complex (MHC) class I and II genes (adaptive immunity) in gulls from seven colonies scattered across Poland. As expected, we found much greater polymorphism at the MHC than TLRs. Population differentiation at the MHC class II, but not MHC-I, was significantly stronger than at neutral microsatellite loci, suggesting local adaptation. This could reflect spatial variation in the composition of extracellular parasite communities (e.g., helminths), possibly driven by sharp differences in habitat structure between colonies. Despite contrasting patterns of population differentiation, both MHC classes showed similar regimes of diversifying selection. Some significant population differentiation was also observed at TLRs, suggesting that innate immune receptors may respond to fine-scale spatial variation in pathogen pressure, although this pattern could have been enhanced by drift. Our results suggested that local adaptation at the pathogen recognition immune genes can be maintained at relatively small or moderate spatial scales in species with high dispersal potential and they highlighted the complexity of immunogenetic responses of animals to heterogeneous environments.

Divergent evolution drives high diversity of toll-like receptors (TLRs) in passerine birds: Buntings and finches

Toll-like receptors (TLRs) form a key component of animal innate immunity, being responsible for recognition of conserved microbial structures. As such, TLRs may be subject to diversifying and balancing selection, which maintains allelic variation both within and between populations. However, most research on TLRs in non-model avian species is focused on bottlenecked populations with depleted genetic variation. Here, we assessed variation at the extracellular domains of three TLR genes (TLR1LA, TLR3, TLR4) across eleven species from two passerine families of buntings (Emberizidae) and finches (Fringillidae), all having large breeding population sizes (millions of individuals). We found extraordinary TLR polymorphism in our study taxa, with >100 alleles detected at TLR1LA and TLR4 across species and high haplotype diversity (>0.75) in several species. Despite recent species divergence, no nucleotide allelic variants were shared between species, suggesting rapid TLR evolution. Higher variation at TLR1LA and TLR4 than TLR3 was associated with a stronger signal of diversifying selection, as measured with nucleotide substitutions rates and the number of positively selected sites (PSS). Structural protein modelling of TLRs showed that some PSS detected within TLR1LA and TLR4 were previously recognized as functionally important sites or were located in their proximity, possibly affecting ligand recognition. Furthermore, we identified PSS responsible for major surface electrostatic charge clustering, which may indicate their adaptive importance. Our study provides compelling evidence for the divergent evolution of TLR genes in buntings and finches and indicates that high TLR variation may be adaptively maintained via diversifying selection acting on functional ligand binding sites.

Understanding the evolution of immune genes in jawed vertebrates

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

Runs of homozygosity reveal past bottlenecks and contemporary inbreeding across diverging populations of an island-colonizing bird

Genomes retain evidence of the demographic history and evolutionary forces that have shaped populations and drive speciation. Across island systems, contemporary patterns of genetic diversity reflect population demography, including colonization events, bottlenecks, gene flow and genetic drift. Here, we investigate genome-wide diversity and the distribution of runs of homozygosity (ROH) using whole-genome resequencing of individuals (>22× coverage) from six populations across three archipelagos of Berthelot's pipit (Anthus berthelotii)-a passerine that has recently undergone island speciation. We show the most dramatic reduction in diversity occurs between the mainland sister species (the tawny pipit) and Berthelot's pipit and is lowest in the populations that have experienced sequential bottlenecks (i.e., the Madeiran and Selvagens populations). Pairwise sequential Markovian coalescent (PSMC) analyses estimated that Berthelot's pipit diverged from its sister species ~2 million years ago, with the Madeiran archipelago founded 50,000 years ago, and the Selvagens colonized 8000 years ago. We identify many long ROH (>1 Mb) in these most recently colonized populations. Population expansion within the last 100 years may have eroded long ROH in the Madeiran archipelago, resulting in a prevalence of short ROH (<1 Mb). However, the extensive long and short ROH detected in the Selvagens suggest strong recent inbreeding and bottleneck effects, with as much as 38% of the autosomes consisting of ROH >250 kb. These findings highlight the importance of demographic history, as well as selection and genetic drift, in shaping contemporary patterns of genomic diversity across diverging populations.

Effect of population size and selection on Toll-like receptor diversity in populations of Galápagos mockingbirds

The interactions of evolutionary forces are difficult to analyse in free-living populations. However, when properly understood, they provide valuable insights into evolutionary biology and conservation genetics. This is particularly important for the interplay of genetic drift and natural selection in immune genes that confer resistance to disease. The Galápagos Islands are inhabited by four closely related species of mockingbirds (Mimus spp.). We used 12 different-sized populations of Galápagos mockingbirds and one population of their continental relative northern mockingbird (Mimus polyglottos) to study the effects of genetic drift on the molecular evolution of immune genes, the Toll-like receptors (TLRs: TLR1B, TLR4 and TLR15). We found that neutral genetic diversity was positively correlated with island size, indicating an important effect of genetic drift. However, for TLR1B and TLR4, there was little correlation between functional (e.g., protein) diversity and island size, and protein structural properties were largely conserved, indicating only a limited effect of genetic drift on molecular phenotype. By contrast, TLR15 was less conserved and even its putative functional polymorphism correlated with island size. The patterns observed for the three genes suggest that genetic drift does not necessarily dominate selection even in relatively small populations, but that the final outcome depends on the degree of selection constraint that is specific for each TLR locus.

Leveraging an existing whole genome resequencing population dataset to characterize toll‐like receptor gene diversity in a threatened bird

Species recovery programs are increasingly using genomic data to measure neutral genetic diversity and calculate metrics like relatedness. While these measures can inform conservation management, determining the mechanisms underlying inbreeding depression requires information about functional genes associated with adaptive or maladaptive traits. Toll‐like receptors (TLRs) are one family of functional genes, which play a crucial role in recognition of pathogens and activation of the immune system. Previously, these genes have been analysed using species‐specific primers and PCR. Here, we leverage an existing short‐read reference genome, whole‐genome resequencing population data set, and bioinformatic tools to characterize TLR gene diversity in captive and wild tchūriwat’/tūturuatu/shore plover (Thinornis novaeseelandiae), a threatened bird endemic to Aotearoa New Zealand. Our results show that TLR gene diversity in tchūriwat’/tūturuatu is low, and forms two distinct captive and wild genetic clusters. The bioinformatic approach presented here has broad applicability to other threatened species with existing genomic resources in Aotearoa New Zealand and beyond.

Genomic associations with poxvirus across divergent island populations in Berthelot's pipit

Understanding the mechanisms and genes that enable animal populations to adapt to pathogens is important from an evolutionary, health and conservation perspective. Berthelot's pipit (Anthus berthelotii) experiences extensive and consistent spatial heterogeneity in avian pox infection pressure across its range of island populations, thus providing an excellent system with which to examine how pathogen-mediated selection drives spatial variation in immunogenetic diversity. Here we test for evidence of genetic variation associated with avian pox at both an individual and population-level. At the individual level, we find no evidence that variation in MHC class I and TLR4 (both known to be important in recognising viral infection) was associated with pox infection within two separate populations. However, using genotype-environment association (Bayenv) in conjunction with genome-wide (ddRAD-seq) data, we detected strong associations between population-level avian pox prevalence and allele frequencies of single nucleotide polymorphisms (SNPs) at a number of sites across the genome. These sites were located within genes involved in cellular stress signalling and immune responses, many of which have previously been associated with responses to viral infection in humans and other animals. Consequently, our analyses indicates that pathogen-mediated selection may play a role in shaping genomic variation among relatively recently colonised island bird populations and highlights the utility of genotype-environment associations for identifying candidate genes potentially involved in host-pathogen interactions.

Stronger population differentiation at infection-sensing than infection-clearing innate immune loci in songbirds: Different selective regimes for different defenses

Parasite-mediated selection is widespread at loci involved in immune defense, but different defenses may experience different selective regimes. For defenses involved in clearing infections, purifying selection favoring a single most efficacious allele likely predominates. However, for defenses involved in sensing and recognizing infections, evolutionary arms races may make positive selection particularly important. This could manifest primarily within populations (e.g., balancing selection maintaining variation) or among them (e.g., spatially varying selection enhancing population differences in allele frequencies). We genotyped three toll-like receptors (TLR; involved in sensing infections) and three avian beta-defensins (involved in clearing infections) in 96 song sparrows (Melospiza melodia) from three breeding populations that differ in disease resistance. Variation-based indicators of selection (proportion of variable sites, proportion of nonsynonymous SNPs, proportion of sites bearing signatures of positive or purifying selection, rare allele frequencies) did not differ appreciably between the two locus types. However, differentiation was generally higher at infection-sensing than infection-clearing loci. Allele frequencies differed markedly at TLR3, driven by a variant predicted to alter protein function. Geographically structured variants at infection-sensing loci may reflect local adaptation to spatially heterogeneous parasite communities. Selective regimes experienced by infection-sensing versus infection-clearing loci may differ primarily due to parasite-mediated population differentiation.

Evolution of innate and adaptive immune genes in a non-model waterbird, the common tern

Toll-like receptors (TLRs) and the Major Histocompatibility Complex (MHC) are the key pathogen-recognition genes of vertebrate immune system and they have a crucial role in the initiation of innate and adaptive immune response, respectively. Recent advancements in sequencing technology sparked research on highly duplicated MHC genes in non-model species, but TLR variation in natural vertebrate populations has remained little studied and comparisons of polymorphism across both TLRs and MHC are scarce. Here, we aimed to compare variation across innate (four TLR loci) and adaptive (MHC class I and class II) immune genes in a non-model avian species, the common tern Sterna hirundo. We detected relatively high allelic richness at TLR genes (9-48 alleles per locus), which was similar to or even higher than the estimated per locus allelic richness at the MHC (24-30 alleles at class I and 13-16 alleles at class II under uniform sample sizes). Despite this, the total number of MHC alleles across all duplicated loci (four class I and three class II) was much higher and MHC alleles showed greater sequence divergence than TLRs. Positive selection targeted relatively more sites at the MHC than TLRs, but the strength of selection (dN/dS ratios) at TLRs was higher when compared to MHC class I. There were also differences in the signature of positive selection and recombination (gene conversion) between MHC class I and II (stronger signature at class II), suggesting that mechanisms maintaining variation at the MHC may vary between both classes. Our study indicates that allelic richness of both innate and adaptive immune receptors may be maintained at relatively high levels in viable avian populations and we recommend a transition from the traditional gene-specific to multi-gene approach in studying molecular evolution of vertebrate immune system.

Contemporary evolution of the innate immune receptor gene TLR3 in an isolated vertebrate population

Understanding where genetic variation exists, and how it influences fitness within populations is important from an evolutionary and conservation perspective. Signatures of past selection suggest that pathogen-mediated balancing selection is a key driver of immunogenetic variation, but studies tracking contemporary evolution are needed to help resolve the evolutionary forces and mechanism at play. Previous work in a bottlenecked population of Seychelles warblers (Acrocephalus sechellensis) show that functional variation has been maintained at the viral-sensing Toll-like receptor 3 (TLR3) gene, including one nonsynonymous SNP, resulting in two alleles. Here, we characterise evolution at this TLR3 locus over a 25-year period within the original remnant population of the Seychelles warbler, and in four other derived, populations. Results show a significant and consistent temporal decline in the frequency of the TLR3^C allele in the original population, and that similar declines in the TLR3^C allele frequency occurred in all the derived populations. Individuals (of both sexes) with the TLR3^CC genotype had lower survival, and males - but not females - that carry the TLR3^C allele had significantly lower lifetime reproductive success than those with only the TLR3^A allele. These results indicate that positive selection on the TLR3^A allele, caused by an as yet unknown agent, is driving TLR3 evolution in the Seychelles warbler. No evidence of heterozygote advantage was detected. However, whether the positive selection observed is part of a longer-term pattern of balancing selection (through fluctuating selection or rare-allele advantage) cannot be resolved without tracking the TLR3^C allele over an extended time period.

Habitat fragmentation differentially shapes neutral and immune gene variation in a tropical bird species

Habitat fragmentation is a major cause of biodiversity loss, responsible for an alteration of intraspecific patterns of neutral genetic diversity and structure. Although neutral genetic variation can be informative for demographic inferences, it may be a poor predictor of adaptive genetic diversity and thus of the consequences of habitat fragmentation on selective evolutionary processes. In this context, we contrasted patterns of genetic diversity and structure of neutral loci (microsatellites) and immune genes (i.e., toll-like receptors) in an understorey bird species, the wedge-billed woodcreeper Glyphorynchus spirurus. The objectives were (1) to investigate forest fragmentation effects on population genetic diversity, (2) to disentangle the relative role of demography (genetic drift and migration) and selection, and (3) to assess whether immunogenetic patterns could be associated with variation of ectoparasite (i.e., ticks) pressures. Our results revealed an erosion of neutral genetic diversity and a substantial genetic differentiation among fragmented populations, resulting from a decrease in landscape connectivity and leading to the divergence of distinct genetic pools at a small spatial scale. Patterns of genetic diversity observed for TLR4 and TLR5 were concordant with neutral genetic patterns, whereas those observed for TLR3 and TLR21 were discordant. This result underlines that the dominant evolutionary force shaping immunogenetic diversity (genetic drift vs. selection) may be different depending on loci considered. Finally, tick prevalence was higher in fragmented environments. We discussed the hypothesis that pathogen selective pressures may contribute to maintain adaptive genetic diversity despite the negative demographic effect of habitat fragmentation on neutral genetic diversity.

Genetic drift and bottleneck do not influence diversity in Toll-like receptor genes at a small spatial scale in a Himalayan passerine

Genetic diversity is important for long-term viability of a population. Low genetic diversity reduces persistence and survival of populations and increases susceptibility to diseases. Comparisons of the neutral markers with functional loci such as immune genes [Toll-like receptors; TLR] can provide useful insights into evolutionary potential of a species and how the diversity of pathogens and selection pressures on their hosts are directly linked to their environment. In this study, we compare genetic diversity in neutral (eleven microsatellite loci) and adaptive (seven TLR loci) loci to determine genetic variation in a nonmigratory western Himalayan passerine, the black-throated tit (Aegithalos concinnus), distributed across an elevation gradient with varying degree of pathogen-mediated selection pressure. We further compare the diversity in TLR loci with a high-elevation sister species, the white-throated tit (Aegithalos niveogularis). Our results indicate a lack of population genetic structure in the black-throated tit and signatures of a past bottleneck. In contrast, we found high diversity in TLR loci and locus-specific (TLR7) signatures of pathogen-mediated selection, which was comparable to diversity in the white-throated tit. Levels of diversity at TLR5 locus corresponded very closely with neutral microsatellite variation. We found evidence of positive selection in TLR1LA, TLR5, and TLR7 loci highlighting the importance in pathogen recognition. Our finding demonstrates that reduction in neutral variation does not necessarily lead to reduction in functional genetic diversity and probably helps in revival of population in a widespread species.

Evidence of Pathogen-Induced Immunogenetic Selection across the Large Geographic Range of a Wild Seabird

Over evolutionary time, pathogen challenge shapes the immune phenotype of the host to better respond to an incipient threat. The extent and direction of this selection pressure depend on the local pathogen composition, which is in turn determined by biotic and abiotic features of the environment. However, little is known about adaptation to local pathogen threats in wild animals. The Gentoo penguin (Pygoscelis papua) is a species complex that lends itself to the study of immune adaptation because of its circumpolar distribution over a large latitudinal range, with little or no admixture between different clades. In this study, we examine the diversity in a key family of innate immune genes-the Toll-like receptors (TLRs)-across the range of the Gentoo penguin. The three TLRs that we investigated present varying levels of diversity, with TLR4 and TLR5 greatly exceeding the diversity of TLR7. We present evidence of positive selection in TLR4 and TLR5, which points to pathogen-driven adaptation to the local pathogen milieu. Finally, we demonstrate that two positively selected cosegregating sites in TLR5 are sufficient to alter the responsiveness of the receptor to its bacterial ligand, flagellin. Taken together, these results suggest that Gentoo penguins have experienced distinct pathogen-driven selection pressures in different environments, which may be important given the role of the Gentoo penguin as a sentinel species in some of the world's most rapidly changing environments.

Genetic diversity of toll-like receptor genes in the vulnerable Chinese egret (Egretta eulophotes)

Toll-like receptor (TLR) genes have recently been employed to assess genetic diversity, as they can be used to infer both demographic history and adaptation to environments with different pathogen pressure. Here, we sampled 120 individuals of the Chinese egret (Egretta eulophotes), a globally vulnerable species, from four breeding populations across China. We assessed the levels of genetic diversity, selection pressure, and population differentiation at seven TLR loci (TLR1LB, TLR2A, TLR3, TLR4, TLR5, TLR7, and TLR15). Using a variety of metrics (SNPs, heterozygosity, nucleotides, haplotypes), our analyses showed that genetic diversity was lower at 4 of the 7 TLR loci in the vulnerable Chinese egret compared to the more common little egret (Egretta garzetta). The selection test indicated TLRs, except for TLR5, were under purifying selection in TLR evolution, suggesting that low TLR genetic diversity in the Chinese egret may be caused by purifying selection. Moreover, analysis of molecular variance indicated low but significant population differentiation among four populations at all of the TLR loci in this egret. However, some comparisons based on fixation index analyses did not show significant population differentiation, and Bayesian clustering showed admixture. Our finding suggested that these four populations of the Chinese egret in China may be considered a single unit for conservation planning. These results, the new report of TLR genetic diversity in a long-distance migratory vulnerable Ardeid species, will provide fundamental TLR information for further studies on the conservation genetics of the Chinese egret and other Ardeids.

Adaptive genetic diversity and evidence of population genetic structure in the endangered Sierra Madre Sparrow (Xenospiza baileyi)

The magnitude and distribution of genetic diversity through space and time can provide useful information relating to evolutionary potential and conservation status in threatened species. In assessing genetic diversity in species that are of conservation concern, several studies have focused on the use of Toll-like receptors (TLRs). TLRs are innate immune genes related to pathogen resistance, and polymorphisms may reflect not only levels of functional diversity, but may also be used to assess genetic diversity within and among populations. Here, we combined four potentially adaptive markers (TLRs) with one mitochondrial (COI) marker to evaluate genetic variation in the endangered Sierra Madre Sparrow (Xenospiza baileyi). This species offers an ideal model to investigate population and evolutionary genetic processes that may be occurring in a habitat restricted endangered species with disjunct populations (Mexico City and Durango), the census sizes of which differ by an order of magnitude. TLRs diversity in the Sierra Madre Sparrow was relatively high, which was not expected given its two small, geographically isolated populations. Genetic diversity was different (but not significantly so) between the two populations, with less diversity seen in the smaller Durango population. Population genetic structure between populations was due to isolation and different selective forces acting on different TLRs; population structure was also evident in COI. Reduction of genetic diversity in COI was observed over 20 years in the Durango population, a result likely caused by habitat loss, a factor which may be the main cause of diversity decline generally. Our results provide information related to the ways in which adaptive variation can be altered by demographic changes due to human-mediated habitat alterations. Furthermore, our findings may help to guide conservation schemes for both populations and their restricted habitat.

Post-glacial phylogeography and variation in innate immunity loci in a sylvatic rodent, bank vole Myodes glareolus

In the northern hemisphere, the spatial structure of many taxa has been shaped by migration patterns after the last glaciation, and phylogeography based on mtDNA variation may reflect the post-glacial demography. The mtDNA lineages are expected to differ in their adaptations to local conditions but little is known about the impact of these conditions on functional genetic variation. Here, we answer this question through an analysis of geographic variation and selection patterns in seven innate immunity genes in free-living bank voles Myodes glareolus from 10 localities across species range assigned to different lineages based on mtDNA. We found clear discrepancies between population structure in mtDNA and each of the studied innate immunity genes. There was no uniform pattern of spatial variation at immunity loci, they differed in the levels of polymorphism, and the results of neutrality tests were not consistent over loci. Each locus comprised a few common haplotypes shared between mitochondrial lineages and studied locations, plus numerous haplotypes unique for each studied site. Our results suggest that the diversity of innate immunity genes cannot be explained solely in terms of demographic processes, and that the observed polymorphism may be attributed to local selection. The strength and direction of selection differed between loci, even within the same gene family, which underlines how crucial it is to take a complex approach while studying the selection patterns acting on immune-related genes.

Toll-like receptor 4: A promising crossroads in the diagnosis and treatment of several pathologies

Toll-like receptor 4 (TLR4) is expressed in a wide variety of cells and is the central component of the mammalian innate immune system. Since its discovery in 1997, TLR4 has been assigned an ever-increasing number of functions that extend from pathogen recognition to tissue damage identification and promotion of the intrinsic "damage repair response" in pain, intestinal, respiratory and vascular disorders. Precisely, the finding of conserved sequence homology among species along with the molecular and functional characterisation of the TLR4 gene enabled researchers to envisage a common operating system in the activation of innate immunity and the initiation of plastic changes at the onset of chronic pain. Malfunctioning in other conditions was conceived in parallel. In this respect, "pivot" proteins and pathway redundancy are not just evolutionary leftovers but essential for normal functioning or cell survival. Indeed, at present, TLR4 single nucleotide polymorphisms (SNP) and their association with certain dysfunctions and diseases are being confirmed in different pools of patients. However, despite its ability to trigger pathogen infection or alternatively tissue injury communications to immune system, TLR4 targeting might not be considered a panacea. This review article represents a compilation of what we know about TLR4 from clinics and basic research on the 20th anniversary of its discovery. Understanding how to fine-tune the interaction between TLR4 and its specific ligands may lead in the next decades to the development of promising new treatments, reducing polypharmacy and probably having an impact on drug use in numerous pathologies.

Relaxation of selective constraints shapes variation of toll-like receptors in a colonial waterbird, the black-headed gull

Nonspecific innate immune response is activated by toll-like receptors (TLRs), which recognize conserved molecular motifs characteristic for a broad spectrum of pathogens. In this study, we examined nucleotide substitution patterns and allelic diversity at five TLR genes in a wild nonpasserine bird, the black-headed gull Chroicocephalus ridibundus. We hypothesized that balancing selection can maintain high allelic diversity of TLR genes in the black-headed gull because of its ecological characteristics, coloniality, and migratoriness, which are associated with increased exposure and transmission of pathogens. Although we found moderately high levels of sequence polymorphism (8-49 haplotypes retrieved per locus within a sample of 60 individuals), most of these haplotypes were recorded at low frequencies within our study population. At the same time, we found no convincing evidence for the role of balancing selection in the maintenance of this variation (Tajima's D < 0.5), and sites with a significant excess of nonsynonymous mutations (dN/dS > 1) were recorded only at two loci (TLR5 and TLR7). This pattern is consistent with relaxation of selective constraints, where most mutations are slightly deleterious and usually removed by purifying selection. No differences in the diversity and nucleotide substitution rates were found between endosomal loci responsible for viral RNA sensing and loci responsible for the recognition of extracellular pathogens. Our study provides the first information on evolutionary mechanisms shaping polymorphism of TLRs in a species from Lari suborder (gulls and allies) and suggests that TLR genes may be poorly responsive to ecological and life-history characteristics of hosts.

Adaptive landscape genetics and malaria across divergent island bird populations

Environmental conditions play a major role in shaping the spatial distributions of pathogens, which in turn can drive local adaptation and divergence in host genetic diversity. Haemosporidians, such as Plasmodium (malaria), are a strong selective force, impacting survival and fitness of hosts, with geographic distributions largely determined by habitat suitability for their insect vectors. Here, we have tested whether patterns of fine-scale local adaptation to malaria are replicated across discrete, ecologically differing island populations of Berthelot's pipits Anthus berthelotii. We sequenced TLR4, an innate immunity gene that is potentially under positive selection in Berthelot's pipits, and two SNPs previously identified as being associated with malaria infection in a genome-wide association study (GWAS) in Berthelot's pipits in the Canary Islands. We determined the environmental predictors of malaria infection, using these to estimate variation in malaria risk on Porto Santo, and found some congruence with previously identified environmental risk factors on Tenerife. We also found a negative association between malaria infection and a TLR4 variant in Tenerife. In contrast, one of the GWAS SNPs showed an association with malaria risk in Porto Santo, but in the opposite direction to that found in the Canary Islands GWAS. Together, these findings suggest that disease-driven local adaptation may be an important factor in shaping variation among island populations.

Purifying Selection in the Toll-Like Receptors of Song Sparrows Melospiza melodia

Variation in immune gene sequences is known to influence resistance to infectious diseases and parasites, and hence survival and mate choice, across animal taxa. Toll-like receptors (TLRs) comprise one essential gene family in the vertebrate innate immune system and recognize evolutionarily conserved structures from all major microorganism classes. However, the causes and consequences of TLR variation in passerine birds remain largely unexplored. We examined 7 TLR genes in song sparrows (Melospiza melodia), a species that is studied across North America. We then examined sequences from 4 unduplicated TLRs (TLR1LB, TLR3, TLR4, and TLR15) from birds in 2 parts of the species' range (N = 27, N = 6), tested for evidence of selection, and conducted pilot analyses of the role of TLR heterozygosity in survival. We identified 45 SNPs: 19 caused changes in amino acid sequences and 2 of these were likely deleterious. We found no evidence of codon-level episodic positive selection but detected purifying selection at codons in all TLRs. Contrary to expectations we found no strong correlation between heterozygosity at TLRs and inbreeding coefficient f (estimate ± standard error [SE] = -0.68 ± 0.37, Radj2 = 0.08, F1,25 = 3.38, P = 0.08). In addition, pilot analyses revealed no relationship between TLR heterozygosity and survival (β ± SE: 0.09 ± 2.00, P = 0.96), possibly due to small sample size. Further analyses of genetic diversity in TLRs are likely to advance understanding of the effects of innate immune gene diversity on the fitness and persistence of wild populations.

Global genetic diversity, lineage distribution, and Wolbachia infection of the alfalfa weevil Hypera postica (Coleoptera: Curculionidae)

The alfalfa weevil (Hypera postica) is a well-known example of a worldwide-distributed pest with high genetic variation. Based on the mitochondrial genes, the alfalfa weevil clusters into two main mitochondrial lineages. However, there is no clear picture of the global diversity and distribution of these lineages; neither the drivers of its diversification are known. However, it appears likely that historic demographic events including founder effects played a role. In addition, Wolbachia, a widespread intracellular parasite/symbiont, likely played an important role in the evolution of the species. Wolbachia infection so far was only detected in the Western lineage of H. postica with no information on the infecting strain, its frequency, and its consequences on the genetic diversity of the host. We here used a combination of mitochondrial and nuclear sequences of the host and sequence information on Wolbachia to document the distribution of strains and the degree of infection. The Eastern lineage has a higher genetic diversity and is found in the Mediterranean, the Middle East, Eastern Europe, and eastern America, whereas the less diverse Western lineage is found in Central Europe and the western America. Both lineages are infected with the same common strain of Wolbachia belonging to Supergroup B. Based on neutrality tests, selection tests, and the current distribution and diversification of Wolbachia in H. postica, we suggested the Wolbachia infection did not shape genetic diversity of the host. The introduced populations in the United States are generally genetically less diverse, which is in line with founder effects.

Adaptive genetic variation at three loci in South African vervet monkeys (Chlorocebus pygerythrus) and the role of selection within primates

Vervet monkeys (Chlorocebus pygerythrus) are one of the most widely distributed non-human primate species found in South Africa. They occur across all the South African provinces, inhabiting a large variety of habitats. These habitats vary sufficiently that it can be assumed that various factors such as pathogen diversity could influence populations in different ways. In turn, these factors could lead to varied levels of selection at specific fitness linked loci. The Toll-like receptor (TLR) gene family, which play an integral role in vertebrate innate immunity, is a group of fitness linked loci which has been the focus of much research. In this study, we assessed the level of genetic variation at partial sequences of two TLR loci (TLR4 and 7) and a reproductively linked gene, acrosin (ACR), across the different habitat types within the vervet monkey distribution range. Gene variation and selection estimates were also made among 11-21 primate species. Low levels of genetic variation for all three gene regions were observed within vervet monkeys, with only two polymorphic sites identified for TLR4, three sites for TLR7 and one site for ACR. TLR7 variation was positively correlated with high mean annual rainfall, which was linked to increased pathogen abundance. The observed genetic variation at TLR4 might have been influenced by numerous factors including pathogens and climatic conditions. The ACR exonic regions showed no variation in vervet monkeys, which could point to the occurrence of a selective sweep. The TLR4 and TLR7 results for the among primate analyses was mostly in line with previous studies, indicating a higher rate of evolution for TLR4. Within primates, ACR coding regions also showed signs of positive selection, which was congruent with previous reports on mammals. Important additional information to the already existing vervet monkey knowledge base was gained from this study, which can guide future research projects on this highly researched taxon as well as help conservation agencies with future management planning involving possible translocations of this species.

Genomic associations with bill length and disease reveal drift and selection across island bird populations

Island species provide excellent models for investigating how selection and drift operate in wild populations, and for determining how these processes act to influence local adaptation and speciation. Here, we examine the role of selection and drift in shaping genomic and phenotypic variation across recently separated populations of Berthelot's pipit (Anthus berthelotii), a passerine bird endemic to three archipelagos in the Atlantic. We first characterized genetic diversity and population structuring that supported previous inferences of a history of recent colonizations and bottlenecks. We then tested for regions of the genome associated with the ecologically important traits of bill length and malaria infection, both of which vary substantially across populations in this species. We identified a SNP associated with variation in bill length among individuals, islands, and archipelagos; patterns of variation at this SNP suggest that both phenotypic and genotypic variation in bill length is largely shaped by founder effects. Malaria was associated with SNPs near/within genes involved in the immune response, but this relationship was not consistent among archipelagos, supporting the view that disease resistance is complex and rapidly evolving. Although we found little evidence for divergent selection at candidate loci for bill length and malaria resistance, genome scan analyses pointed to several genes related to immunity and metabolism as having important roles in divergence and adaptation. Our findings highlight the utility and challenges involved with combining association mapping and population genetic analysis in nonequilibrium populations, to disentangle the effects of drift and selection on shaping genotypes and phenotypes.

MHC class II DRB diversity predicts antigen recognition and is associated with disease severity in California sea lions naturally infected with Leptospira interrogans

We examined the associations between California sea lion MHC class II DRB (Zaca-DRB) configuration and diversity, and leptospirosis. As Zaca-DRB gene sequences are involved with antigen presentation of bacteria and other extracellular pathogens, we predicted that they would play a role in determining responses to these pathogenic spirochaetes. Specifically, we investigated whether Zaca-DRB diversity (number of genes) and configuration (presence of specific genes) explained differences in disease severity, and whether higher levels of Zaca-DRB diversity predicted the number of specific Leptospira interrogans serovars that a sea lion's serum would react against. We found that serum from diseased sea lions with more Zaca-DRB loci reacted against a wider array of serovars. Specific Zaca-DRB loci were linked to reactions with particular serovars. Interestingly, sea lions with clinical manifestation of leptospirosis that had higher numbers of Zaca-DRB loci were less likely to recover from disease than those with lower diversity, and those that harboured Zaca-DRB.C or -G were 4.5 to 5.3 times more likely to die from leptospirosis, regardless of the infective serovars. We propose that for leptospirosis, a disadvantage of having a wider range of antigen presentation might be increased disease severity due to immunopathology. Ours is the first study to examine the importance of Zaca-DRB diversity for antigen detection and disease severity following natural exposure to infective leptospires.

Purifying selection and concerted evolution of RNA-sensing toll-like receptors in migratory waders

Migratory birds encounter a broad range of pathogens during their journeys, making them ideal models for studying immune gene evolution. Despite the potential value of these species to immunoecology and disease epidemiology, previous studies have typically focused on their adaptive immune gene repertoires. In this study, we examined the evolution of innate immune genes in three long-distance migratory waders (order Charadriiformes). We analysed two parts of the extracellular domains of two Toll-like receptors (TLR3 and TLR7) involved in virus recognition in the Sanderling (Calidris alba), Red-necked Stint (Calidris ruficollis), and Ruddy Turnstone (Arenaria interpres). Our analysis was extended to 50 avian species for which whole-genome sequences were available, including two additional waders. We found that the inferred relationships among avian TLR3 and TLR7 do not match the whole-genome phylogeny of birds. Further analyses showed that although both loci are predominantly under purifying selection, the evolution of the extracellular domain of avian TLR3 has also been driven by episodic diversifying selection. TLR7 was found to be duplicated in all five wader species and in two other orders of birds, Cuculiformes and Passeriformes. The duplication is likely to have occurred in the ancestor of each order, and the duplicated copies appear to be undergoing concerted evolution. The phylogenetic relationships of wader TLR7 matched those of the five wader species, but that of TLR3 did not. Instead, the tree inferred from TLR3 showed potential associations with the species' ecology, including migratory behaviour and exposure to pathogens. Our study demonstrates the importance of combining immunological and ecological knowledge to understand the impact of immune gene polymorphism on the evolutionary ecology of infectious diseases.

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.

Correlated patterns of genetic diversity and differentiation across an avian family

Comparative studies of closely related taxa can provide insights into the evolutionary forces that shape genome evolution and the prevalence of convergent molecular evolution. We investigated patterns of genetic diversity and differentiation in stonechats (genus Saxicola), a widely distributed avian species complex with phenotypic variation in plumage, morphology and migratory behaviour, to ask whether similar genomic regions have become differentiated in independent, but closely related, taxa. We used whole-genome pooled sequencing of 262 individuals from five taxa and found that levels of genetic diversity and divergence are strongly correlated among different stonechat taxa. We then asked whether these patterns remain correlated at deeper evolutionary scales and found that homologous genomic regions have become differentiated in stonechats and the closely related Ficedula flycatchers. Such correlation across a range of evolutionary divergence and among phylogenetically independent comparisons suggests that similar processes may be driving the differentiation of these independently evolving lineages, which in turn may be the result of intrinsic properties of particular genomic regions (e.g. areas of low recombination). Consequently, studies employing genome scans to search for areas important for reproductive isolation or adaptation should account for corresponding regions of differentiation, as these regions may not necessarily represent speciation islands or evidence of local adaptation.

Balancing selection and genetic drift create unusual patterns of MHCIIβ variation in Galápagos mockingbirds

The extracellular subunit of the major histocompatibility complex MHCIIβ plays an important role in the recognition of pathogens and the initiation of the adaptive immune response of vertebrates. It is widely accepted that pathogen-mediated selection in combination with neutral micro-evolutionary forces (e.g. genetic drift) shape the diversity of MHCIIβ, but it has proved difficult to determine the relative effects of these forces. We evaluated the effect of genetic drift and balancing selection on MHCIIβ diversity in 12 small populations of Galápagos mockingbirds belonging to four different species, and one larger population of the Northern mockingbird from the continental USA. After genotyping MHCIIβ loci by high-throughput sequencing, we applied a correlational approach to explore the relationships between MHCIIβ diversity and population size by proxy of island size. As expected when drift predominates, we found a positive effect of population size on the number of MHCIIβ alleles present in a population. However, the number of MHCIIβ alleles per individual and number of supertypes were not correlated with population size. This discrepancy points to an interesting feature of MHCIIβ diversity dynamics: some levels of diversity might be shaped by genetic drift while others are independent and possibly maintained by balancing selection.