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

Natural selection directing molecular evolution in vertebrate viral sensors

Diseases caused by pathogens contribute to molecular adaptations in host immunity. Variety of viral pathogens challenging animal immunity can drive positive selection diversifying receptors recognising the infections. However, whether distinct virus sensing systems differ across animals in their evolutionary modes remains unclear. Our review provides a comparative overview of natural selection shaping molecular evolution in vertebrate viral-binding pattern recognition receptors (PRRs). Despite prevailing negative selection arising from the functional constraints, multiple lines of evidence now suggest diversifying selection in the Toll-like receptors (TLRs), NOD-like receptors (NLRs), RIG-I-like receptors (RLRs) and oligoadenylate synthetases (OASs). In several cases, location of the positively selected sites in the ligand-binding regions suggests effects on viral detection although experimental support is lacking. Unfortunately, in most other PRR families including the AIM2-like receptor family, C-type lectin receptors (CLRs), and cyclic GMP-AMP synthetase studies characterising their molecular evolution are rare, preventing comparative insight. We indicate shared characteristics of the viral sensor evolution and highlight priorities for future research.

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.

Immune genotypes, immune responses, and survival in a wild bird population

Individuals vary in their immune genotype, inbreeding coefficient f, immune responses, survival to adulthood, and adult longevity. However, whether immune genes predict survival or longevity, whether such relationships are mediated through immune responses, and how f affects immune genotype remain unclear. We use a wild song sparrow (Melospiza melodia) population in which survival to adulthood, adult longevity, and f were measured precisely, and in which immune responses have previously been assessed. We investigate four toll-like receptor (TLR) and the major histocompatibility complex (MHC) class IIB exon 2 genes. We test whether immune genes predict fitness (survival to adulthood or adult longevity); whether immune genes predict immune response; whether immune response predicts fitness and whether fitness, immune responses, or immune genotypes are correlated with f. We find that survival to adulthood is not associated with immune gene variation, but adult longevity is decreased by high MHC allele diversity (especially in birds that were relatively outbred), and by the presence of a specific MHC supertype. Immune responses were affected by specific immune genotypes. Survival to adulthood and adult longevity were not predicted by immune response, implying caution in the use of immune response as a predictor for fitness. We also found no relationship between f and immune genotype. This finding indicates that immune gene associations with longevity and immune response are not artefacts of f, and suggests that pathogen-mediated selection at functional loci can slow the loss of genetic variation arising from genetic drift and small population size.

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.

Adaptation to climate change through seasonal migration revealed by climatic versus demographic niche models

Predicting the geographic range of species and their response to climatic variation and change are entwined goals in conservation and evolutionary ecology. Species distribution models (SDMs) are foundational in this effort and used to visualize the geographic range of species as the spatial representation of its realized niche. SDMs are also used to forecast range shifts under climate change, but often in the absence of empirical evidence that climate limits population growth. We explored the influence of climate on demography, seasonal migration, and the extent of the geographic range in song sparrows (Melospiza melodia), a species thought to display marked local adaptation to regional climate. To do so, we developed SDMs to predict the demographic and climate niches of migratory and resident song sparrows across our study area in western North America from California to Alaska, using 48 years of demographic data from a focal population in British Columbia and 1.2 million continental-scale citizen science observations. Spatial agreement of our demographic and climate niche models in the region of our focal population was strong (76%), supporting the hypothesis that demographic performance and the occurrence of seasonal migration varied predictably with climatic conditions. In contrast, agreement at the northern (58%) and southern (40%) extents of our study area was lower, as expected if the factors limiting population growth vary regionally. Our results support the hypothesis that local climate drives spatial variation in the occurrence of seasonal migration in song sparrows by limiting the fitness of year-round residents, and suggest that climate warming has favored range expansions and facilitated an upward shift in elevational range song sparrows that forgo seasonal migration. Our work highlights the potential role of seasonal migration in climate adaptation and limits on the reliability of climate niche models not validated with demographic data.

Selection balancing at innate immune genes: adaptive polymorphism maintenance in Toll-like receptors

Balancing selection is a classic mechanism for maintaining variability in immune genes involved in host–pathogen interactions. However, it remains unclear how widespread the mechanism is across immune genes other than the major histocompatibility complex (MHC). Although occasional reports suggest that balancing selection (heterozygote advantage, negative frequency-dependent selection, and fluctuating selection) may act on other immune genes, the current understanding of the phenomenon in non-MHC immune genes is far from solid. In this review, we focus on Toll-like receptors (TLRs), innate immune genes directly involved in pathogen recognition and immune response activation, as there is a growing body of research testing the assumptions of balancing selection in these genes. After reviewing infection- and fitness-based evidence, along with evidence based on population allelic frequencies and heterozygosity levels, we conclude that balancing selection maintains variation in TLRs, though it tends to occur under specific conditions in certain evolutionary lineages rather than being universal and ubiquitous. Our review also identifies key gaps in current knowledge and proposes promising areas for future research. Improving our understanding of host–pathogen interactions and balancing selection in innate immune genes are increasingly important, particularly regarding threats from emerging zoonotic diseases.

Genomic signatures of inbreeding in a critically endangered parrot, the kākāpō

Events of inbreeding are inevitable in critically endangered species. Reduced population sizes and unique life-history traits can increase the severity of inbreeding, leading to declines in fitness and increased risk of extinction. Here, we investigate levels of inbreeding in a critically endangered flightless parrot, the kākāpō (Strigops habroptilus), wherein a highly inbred island population and one individual from the mainland of New Zealand founded the entire extant population. Genotyping-by-sequencing (GBS), and a genotype calling approach using a chromosome-level genome assembly, identified a filtered set of 12,241 single-nucleotide polymorphisms (SNPs) among 161 kākāpō, which together encompass the total genetic potential of the extant population. Multiple molecular-based estimates of inbreeding were compared, including genome-wide estimates of heterozygosity (FH), the diagonal elements of a genomic-relatedness matrix (FGRM), and runs of homozygosity (RoH, FRoH). In addition, we compared levels of inbreeding in chicks from a recent breeding season to examine if inbreeding is associated with offspring survival. The density of SNPs generated with GBS was sufficient to identify chromosomes that were largely homozygous with RoH distributed in similar patterns to other inbred species. Measures of inbreeding were largely correlated and differed significantly between descendants of the two founding populations. However, neither inbreeding nor ancestry was found to be associated with reduced survivorship in chicks, owing to unexpected mortality in chicks exhibiting low levels of inbreeding. Our study highlights important considerations for estimating inbreeding in critically endangered species, such as the impacts of small population sizes and admixture between diverse lineages.

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.

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.

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.

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.

Molecular evolution in immune genes across the avian tree of life

All organisms encounter pathogens, and birds are especially susceptible to infection by malaria parasites and other haemosporidians. It is important to understand how immune genes, primarily innate immune genes which are the first line of host defense, have evolved across birds, a highly diverse group of tetrapods. Here, we find that innate immune genes are highly conserved across the avian tree of life and that although most show evidence of positive or diversifying selection within specific lineages or clades, the number of sites is often proportionally low in this broader context of putative constraint. Rather, the evidence shows a much higher level of negative or purifying selection in these innate immune genes - rather than adaptive immune genes - which is consistent with birds' long coevolutionary history with pathogens and the need to maintain a rapid response to infection. We further explored avian responses to haemosporidians by comparing differential gene expression in wild birds (1) uninfected with haemosporidians, (2) infected with Plasmodium, and (3) infected with Haemoproteus (Parahaemoproteus). We found patterns of significant differential expression with some genes unique to infection with each genus and a few shared between "treatment" groups, but none that overlapped with the genes included in the phylogenetic study.