Signatures of selection acting on the innate immunity gene Toll-like receptor 2 (TLR2) during the evolutionary history of rodents

Patterns of selection acting on immune defence genes have recently been the focus of considerable interest. Yet, when it comes to vertebrates, studies have mainly focused on the acquired branch of the immune system. Consequently, the direction and strength of selection acting on genes of the vertebrate innate immune defence remain poorly understood. Here, we present a molecular analysis of selection on an important receptor of the innate immune system of vertebrates, the Toll-like receptor 2 (TLR2), across 17 rodent species. Although purifying selection was the prevalent evolutionary force acting on most parts of the rodent TLR2, we found that codons in close proximity to pathogen-binding and TLR2-TLR1 heterodimerization sites have been subject to positive selection. This indicates that parasite-mediated selection is not restricted to acquired immune system genes like the major histocompatibility complex, but also affects innate defence genes. To obtain a comprehensive understanding of evolutionary processes in host-parasite systems, both innate and acquired immunity thus need to be considered.     

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.     

Borrelia burgdorferi sensu lato infection pressure shapes innate immune gene evolution in natural rodent populations across Europe

Although parasite-mediated selection is assumed to be the main driver of immune gene evolution, empirical evidence that parasites induce allele frequency changes at host immune genes in time and/or space remains scarce. Here, I show that the frequency of a protective gene variant of the innate immune receptor Toll-like receptor 2 in natural bank vole (Myodes glareolus) populations is positively associated with the strength of Borrelia burgdorferi sensu lato infection risk across the European continent. Thereby, this study provides rare evidence for the role of spatially variable infection pressures in moulding the vertebrate immune system.     

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.     

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.     

Sexual selection explains more functional variation in the mammalian major histocompatibility complex than parasitism

Understanding drivers of genetic diversity at the major histocompatibility complex (MHC) is vitally important for predicting how vertebrate immune defence might respond to future selection pressures and for preserving immunogenetic diversity in declining populations. Parasite-mediated selection is believed to be the major selective force generating MHC polymorphism, and while MHC-based mating preferences also exist for multiple species including humans, the general importance of mate choice is debated. To investigate the contributions of parasitism and sexual selection in explaining among-species variation in MHC diversity, we applied comparative methods and meta-analysis across 112 mammal species, including carnivores, bats, primates, rodents and ungulates. We tested whether MHC diversity increased with parasite richness and relative testes size (as an indicator of the potential for mate choice), while controlling for phylogenetic autocorrelation, neutral mutation rate and confounding ecological variables. We found that MHC nucleotide diversity increased with parasite richness for bats and ungulates but decreased with parasite richness for carnivores. By contrast, nucleotide diversity increased with relative testes size for all taxa. This study provides support for both parasite-mediated and sexual selection in shaping functional MHC polymorphism across mammals, and importantly, suggests that sexual selection could have a more general role than previously thought.     

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.     

TEMPORAL VARIATION AT THE MHC CLASS IIB IN WILD POPULATIONS OF THE GUPPY (POECILIA RETICULATA)

Understanding genetic diversity in natural populations is a fundamental objective of evolutionary biology. The immune genes of the major histocompatibility complex (MHC) are excellent candidates to study such diversity because they are highly polymorphic in populations. Although balancing selection may be responsible for maintaining diversity at these functionally important loci, temporal variation in selection pressure has rarely been examined. We examine temporal variation in MHC class IIB diversity in nine guppy (Poecilia reticulata) populations over two years. We found that five of the populations changed significantly more at the MHC than at neutral (microsatellite) loci as measured by F_ST, which suggests that the change at the MHC was due to selection and not neutral processes. Additionally, pairwise population differentiation measures at the MHC were higher in 2007 than in 2006, with the signature of selection changing from homogenizing to diversifying selection or neutral evolution. Interestingly, within the populations the magnitude of the change at the MHC between years was related to the change in the proportion of individuals infected by a common parasite, indicating a link between genetic structure and the parasite. Our data thereby implicate temporal variation in selective pressure as an important mechanism maintaining diversity at the MHC in wild populations.     

Patterns of selection and allele diversity of class I and class II major histocompatibility loci across the species range of sockeye salmon (Oncorhynchus nerka)

The major histocompatibility complex (MHC), an important component of the vertebrate immune system, provides an important suite of genes to examine the role of genetic diversity at non‐neutral loci for population persistence. We contrasted patterns of diversity at the two classical MHC loci in sockeye salmon (Oncorhynchus nerka), MHC class I (UBA) and MHC class II (DAB), and neutral microsatellite loci across 70 populations spanning the species range from Washington State to Japan. There was no correlation in allelic richness or heterozygosity between MHC loci or between MHC loci and microsatellites. The two unlinked MHC loci may be responding to different selective pressures; the distribution of F_ST values for the two loci was uncorrelated, and evidence for both balancing and directional selection on alleles and lineages of DAB and UBA was observed in populations throughout the species range but rarely on both loci within a population. These results suggest that fluctuating selection has resulted in the divergence of MHC loci in contemporary populations.     

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.Subject terms: Tropical ecology, Genetic variation     

Wild cyclic voles maintain high neutral and MHC diversity without strong evidence for parasite-mediated selection

The major histocompatibility complex (MHC) is an important component of vertebrate immune defense involved with self/nonself recognition and disease susceptibility. The high variability of genes of the MHC is thought to arise from both parasite-mediated and sexual selection. An outstanding question involves the degree to which balancing selection can oppose genetic drift to maintain high MHC diversity in the face of population bottlenecks. To address this question we examined genetic diversity and population structure at neutral (microsatellite) and MHC genes in montane voles [Microtus montanus (Peale, 1848)] subject to high amplitude population fluctuations, and compared these to measures of infection by common gastrointestinal parasites. We found high neutral and MHC allelic variability, indicating low impacts of genetic drift despite large fluctuations in population size. Greater MHC diversity did not predict lower parasite richness or infection by the two most common endoparasites (cestodes and coccidian protozoa), as might be expected if genotypic composition confers resistance to infection. One specific MHC allele predicted lower cestode intensity, but we found no other associations between MHC and infection measures. Neutral heterozygosity was positively associated with total parasite richness, possibly owing to greater parasite tolerance among heterozygous relative to more inbred hosts. Overall, these results suggest that factors beyond the parasites examined here, such as high inter-patch migration, mate choice, gene conversion or other infectious agents, are likely maintaining the high levels of MHC diversity observed in wild montane voles.     

Geographic variation in adaptation at the molecular level: a case study of plant immunity genes

Natural selection imposed by interacting species frequently varies among geographic locations and can lead to local adaptation, where alternative phenotypes are found in different populations. Little is known, however, about whether geographically variable selection acting on traits that mediate species interactions is consistent or strong enough to influence patterns of nucleotide variation at individual loci. To investigate this question, we examined patterns of nucleotide diversity and population structure at 16 plant innate immunity genes, with putative functions in defending plants against pathogens or herbivores, from six populations of teosinte (Zea mays ssp. parviglumis). Specifically, we tested whether patterns of population structure and within-population diversity at immunity genes differed from patterns found at nonimmunity (reference) loci and from neutral expectations derived from coalescent simulations of structured populations. For the majority of genes, we detected no strong evidence of geographically variable selection. However, in the wound-induced serine protease inhibitor (wip1), which inhibits the hydrolysis of dietary proteins in insect herbivores, one population showed unusually high levels of genetic differentiation, very low levels of nucleotide polymorphism, and was fixed for a novel replacement substitution in the active site of the protein. Taken together, these data suggest that wip1 experienced a recent selective sweep in one geographic region; this pattern may reflect local adaptation or an ongoing species-wide sweep. Overall, our results indicate that a signature of local adaptation at the molecular level may be uncommon-particularly for traits that are under complex genetic control.     

Unraveling the Effects of Selection and Demography on Immune Gene Variation in Free-Ranging Plains Zebra (Equus quagga) Populations

Demography, migration and natural selection are predominant processes affecting the distribution of genetic variation among natural populations. Many studies use neutral genetic markers to make inferences about population history. However, the investigation of functional coding loci, which directly reflect fitness, is critical to our understanding of species'' ecology and evolution. Immune genes, such as those of the Major Histocompatibility Complex (MHC), play an important role in pathogen recognition and provide a potent model system for studying selection. We contrasted diversity patterns of neutral data with MHC loci, ELA-DRA and -DQA, in two southern African plains zebra (Equus quagga) populations: Etosha National Park, Namibia, and Kruger National Park, South Africa. Results from neutrality tests, along with observations of elevated diversity and low differentiation across populations, supported previous genus-level evidence for balancing selection at these loci. Despite being low, MHC divergence across populations was significant and may be attributed to drift effects typical of geographically separated populations experiencing little to no gene flow, or alternatively to shifting allele frequency distributions driven by spatially variable and fluctuating pathogen communities. At the DRA, zebra exhibited geographic differentiation concordant with microsatellites and reduced levels of diversity in Etosha due to highly skewed allele frequencies that could not be explained by demography, suggestive of spatially heterogeneous selection and local adaptation. This study highlights the complexity in which selection affects immune gene diversity and warrants the need for further research on the ecological mechanisms shaping patterns of adaptive variation among natural populations.     

Geographic heterogeneity in natural selection on an MHC locus in sockeye salmon

Balancing selection maintains high levels of polymorphism and heterozygosity in genes of the MHC (major histocompatibility complex) of vertebrate organisms, and promotes long evolutionary persistence of individual alleles and strongly differentiated allelic lineages. In this study, genetic variation at the MHC class II DAB-beta1 locus was examined in 31 populations of sockeye salmon (Oncorhynchus nerka) inhabiting the Fraser River drainage of British Columbia, Canada. Twenty-five percent of variation at the locus was partitioned among sockeye populations, as compared with 5% at neutral genetic markers. Geographic heterogeneity of balancing selection was detected among four regions in the Fraser River drainage and among lake systems within regions. High levels of beta1 allelic diversity and heterozygosity, as well as distributions of alleles and allelic lineages that were more even than expected for a neutral locus, indicated the presence of balancing selection in populations throughout much of the interior Fraser drainage. However, proximate populations in the upper Fraser region, and four of six populations from the lower Fraser drainage, exhibited much lower levels of genetic diversity and had beta1 allele frequency distributions in conformance with those expected for a neutral locus, or a locus under directional selection. Pair-wise FST values for beta1 averaged 0.19 and tended to exceed the corresponding values estimated for neutral loci at all levels of population structure, although they were lower among populations experiencing balancing selection than among other populations. The apparent heterogeneity in selection resulted in strong genetic differentiation between geographically proximate populations with and without detectable levels of balancing selection, in stark contrast to observations at neutral loci. The strong partitioning and complex structure of beta1 diversity within and among sockeye populations on a small geographic scale illustrates the value of incorporating adaptive variation into conservation planning for the species.     

Genetic structure in insular and mainland populations of house sparrows (Passer domesticus) and their hemosporidian parasites

Small and isolated populations usually exhibit low levels of genetic variability, and thus, they are expected to have a lower capacity to adapt to changes in environmental conditions, such as exposure to pathogens and parasites. Comparing the genetic variability of selectively neutral versus functional loci allows one to assess the evolutionary history of populations and their future evolutionary potential. The genes of the major histocompatibility complex (MHC) control immune recognition of parasites, and their unusually high diversity is genes which is likely driven by parasite‐mediated balancing selection. Here, we examined diversity and differentiation of neutral microsatellite loci and functional MHC class I genes in house sparrows (Passer domesticus), living in six insular and six mainland populations, and we aimed to determine whether their diversity or differentiation correlates with the diversity and the prevalence of infection of hemosporidian parasites. We found that island bird populations tended to have lower neutral genetic variability, whereas MHC variability gene was similar between island and mainland populations. Similarly, island populations tended to show greater genetic differentiation than mainland populations, especially at microsatellite markers. The maintenance of MHC genetic diversity and its less marked structure in the island populations could be attributed to balancing‐selection. The greater MHC differentiation among populations was negatively correlated with similarity in blood parasites (prevalence and diversity of parasite strains) between populations. Even at low prevalence and small geographical scale, haemosporidian parasites might contribute to structure the variability of immune genes among populations of hosts.     

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.     

MHC-mediated local adaptation in reciprocally translocated Chinook salmon

Most Pacific salmonid populations have faced significant population declines over the past 30 years. In order to effectively conserve and manage these populations, knowledge of the evolutionary adaptive state of individuals and the scale of adaptation across populations is needed. The vertebrate major histocompatibility complex (MHC) represents an important adaptation to parasites, and genes encoding for the MHC are widely held to be undergoing balancing selection. However, the generality of balancing selection across populations at MHC loci is not well documented. Using Chinook salmon (Oncorhynchus tshawytscha) from two populations, we follow the survival of full-sib family replicates reared in their natal river and reciprocally transplanted to a foreign river to examine selection and local adaptation at the MHC class I and II loci. In both populations, we found evidence of a survivorship advantage associated with nucleotide diversity at the MHC class I locus. In contrast, we found evidence that MHC class II diversity was disadvantageous in one population. There was no evidence that these effects occurred in translocated families, suggesting some degree of local adaptation at the MHC loci. Thus, our results implicate balancing selection at the MHC class I but potentially differing selection across populations at the class II locus.     

Genome-wide polymorphisms show unexpected targets of natural selection

Natural selection can act on all the expressed genes of an individual, leaving signatures of genetic differentiation or diversity at many loci across the genome. New power to assay these genome-wide effects of selection comes from associating multi-locus patterns of polymorphism with gene expression and function. Here, we performed one of the first genome-wide surveys in a marine species, comparing purple sea urchins, Strongylocentrotus purpuratus, from two distant locations along the species' wide latitudinal range. We examined 9112 polymorphic loci from upstream non-coding and coding regions of genes for signatures of selection with respect to gene function and tissue- and ontogenetic gene expression. We found that genetic differentiation (F(ST)) varied significantly across functional gene classes. The strongest enrichment occurred in the upstream regions of E3 ligase genes, enzymes known to regulate protein abundance during development and environmental stress. We found enrichment for high heterozygosity in genes directly involved in immune response, particularly NALP genes, which mediate pro-inflammatory signals during bacterial infection. We also found higher heterozygosity in immune genes in the southern population, where disease incidence and pathogen diversity are greater. Similar to the major histocompatibility complex in mammals, balancing selection may enhance genetic diversity in the innate immune system genes of this invertebrate. Overall, our results show that how genome-wide polymorphism data coupled with growing databases on gene function and expression can combine to detect otherwise hidden signals of selection in natural populations.     

Signals of heterogeneous selection at an MHC locus in geographically proximate ecotypes of sockeye salmon

The genes of the major histocompatibility complex (MHC) are an important component of the vertebrate immune system and can provide insights into the role of pathogen‐mediated selection in wild populations. Here, we examined variation at the MHC class II peptide‐binding region in 27 populations of sockeye salmon (Oncorhynchus nerka), distributed among three distinct spawning ecotypes, from a complex of interconnected rivers and lakes in south‐western Alaska. We also obtained genotypes from 90 putatively neutral single nucleotide polymorphisms for each population to compare the relative roles of demography and selection in shaping the observed MHC variation. We found that MHC divergence was generally partitioned by spawning ecotype (lake beaches, rivers and streams) and was 30 times greater than variation at neutral markers. Additionally, we observed substantial differences in modes of selection and diversity among ecotypes, with beach populations displaying higher levels of directional selection and lower MHC diversity than the other two ecotypes. Finally, the level of MHC differentiation in our study system was comparable to that observed over much larger geographic ranges, suggesting that MHC variation does not necessarily increase with increasing spatial scale and may instead be driven by fine‐scale differences in pathogen communities or pathogen virulence. The low levels of neutral structure and spatial proximity of populations in our study system indicate that MHC differentiation can be maintained through strong selective pressure even when ample opportunities for gene flow exist.     

A temporal analysis shows major histocompatibility complex loci in the Scandinavian wolf population are consistent with neutral evolution

The major histocompatibility complex (MHC) has an integral role in the immune system, and hence diversity at its genes may be of particular importance for the health of populations. In large populations, balancing selection maintains diversity in MHC genes, but theoretical expectations indicate that this form of selection is absent or inefficient in small populations. We examine the level of diversity at three MHC class II loci in the wolf population of Scandinavia, a population naturally recolonized with a genetic contribution from as few as three founders, and in four neighbouring wolf populations. In the Scandinavian wolf population, two alleles were found for each locus and the distribution of alleles is compatible with their linkage into two haplotypes. Changes in the level of heterozygosity over time since recolonization demonstrate the effects of the proposed arrival of an immigrant wolf. The maintenance of diversity is shown to be compatible with a neutral, random allocation of alleles, in conjunction with crossing between packs. A total of 15 DRB1, seven DQA and 10 DQB1 alleles are found in four neighbouring wolf populations, with substantial sharing across populations. Even in these larger populations, bottlenecks and fragmentation with consequent genetic drift are likely to have resulted in few indicators for balancing selection and significant differentiation of populations.