Genetic analysis of differentiation among breeding ponds reveals a candidate gene for local adaptation in Rana arvalis

One of the main questions in evolutionary and conservation biology is how geographical and environmental features of the landscape shape neutral and adaptive genetic variation in natural populations. The identification of genomic polymorphisms that account for adaptive variation can aid in finding candidate loci for local adaptation. Consequently, a comparison of spatial patterns in neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection at the landscape scale. Many amphibians breed in wetlands, which differ in environmental conditions and in the degree of isolation, enhancing the potential for local adaptation. We used microsatellite markers to measure genetic differentiation among 17 local populations of Rana arvalis breeding in a network of wetlands. We found that locus RC08604 deviated from neutral expectations, suggesting that it is a good candidate for directional selection. We used a genetic network analysis to show that the allele distribution in this locus is correlated with habitat characteristics, whereas this was not the case at neutral markers that displayed a different allele distribution and population network in the study area. The graph approach illustrated the genomic heterogeneity (neutral loci vs. the candidate locus for directional selection) of gene exchange and genetic divergence among populations under directional selection. Limited gene flow between wetlands was only observed at the candidate genomic region under directional selection. RC08604 is partially located inside an up‐regulated thyroid‐hormone receptor (TRβ) gene coordinating the expression of other genes during metamorphosis and appears to be linked with variation in larval life‐history traits found among R. arvalis populations. We suggest that directional selection on genes coding larval life‐history traits is strong enough to maintain the divergence in these genomic regions, reducing the effective recombination of locally adapted alleles but not in other regions of the genome. Integrating this knowledge into conservation plans at the landscape scale will improve the design of management strategies to preserve adaptive genetic diversity in wetland networks.     

On the relative roles of selection and genetic drift in shaping MHC variation

Genes of the major histocompatibility complex (MHC) have provided some of the clearest examples of how natural selection generates discordances between adaptive and neutral variation in natural populations. The type and intensity of selection as well as the strength of genetic drift are believed to be important in shaping the resulting pattern of MHC diversity. However, evaluating the relative contribution of multiple microevolutionary forces is challenging, and empirical studies have reported contrasting results. For instance, balancing selection has been invoked to explain high levels of MHC diversity and low population differentiation in comparison with other nuclear markers. Other studies have shown that genetic drift can sometimes overcome selection and then patterns of genetic variation at adaptive loci cannot be discerned from those occurring at neutral markers. Both empirical and simulated data also indicate that loss of genetic diversity at adaptive loci can occur faster than at neutral loci when selection and population bottlenecks act simultaneously. Diversifying selection, on the other hand, explains accelerated MHC divergence as the result of spatial variation in pathogen‐mediated selective regimes. Because of all these possible scenarios and outcomes, collecting information from as many study systems as possible, is crucial to enhance our understanding about the evolutionary forces driving MHC polymorphism. In this issue, Miller and co‐workers present an illuminating contribution by combining neutral markers (microsatellites) and adaptive MHC class I loci during the investigation of genetic differentiation across island populations of tuatara Sphenodon punctatus. Their study of geographical variation reveals a major role of genetic drift in shaping MHC variation, yet they also discuss some support for diversifying selection.     

Bayesian estimation of genomic clines

We developed a Bayesian genomic cline model to study the genetic architecture of adaptive divergence and reproductive isolation between hybridizing lineages. This model quantifies locus‐specific patterns of introgression with two cline parameters that describe the probability of locus‐specific ancestry as a function of genome‐wide admixture. ‘Outlier’ loci with extreme patterns of introgression relative to most of the genome can be identified. These loci are potentially associated with adaptive divergence or reproductive isolation. We simulated genetic data for admixed populations that included neutral introgression, as well as loci that were subject to directional, epistatic or underdominant selection, and analysed these data using the Bayesian genomic cline model. Under many demographic conditions, underdominance or directional selection had detectable and predictable effects on cline parameters, and ‘outlier’ loci were greatly enriched for genetic regions affected by selection. We also analysed previously published genetic data from two transects through a hybrid zone between Mus domesticus and M. musculus. We found considerable variation in rates of introgression across the genome and particularly low rates of introgression for two X‐linked markers. There were similarities and differences in patterns of introgression between the two transects, which likely reflects a combination of stochastic variability because of genetic drift and geographic variation in the genetic architecture of reproductive isolation. By providing a robust framework to quantify and compare patterns of introgression among genetic regions and populations, the Bayesian genomic cline model will advance our understanding of the genetics of reproductive isolation and the speciation process.     

Selection and drift influence genetic differentiation of insular Canada lynx (Lynx canadensis) on Newfoundland and Cape Breton Island

Island populations have long been important for understanding the dynamics and mechanisms of evolution in natural systems. While genetic drift is often strong on islands due to founder events and population bottlenecks, the strength of selection can also be strong enough to counteract the effects of drift. Here, we used several analyses to identify the roles of genetic drift and selection on genetic differentiation and diversity of Canada lynx (Lynx canadensis) across eastern Canada, including the islands of Cape Breton and Newfoundland. Specifically, we assessed whether we could identify a genetic component to the observed morphological differentiation that has been reported across insular and mainland lynx. We used a dinucleotide repeat within the promoter region of a functional gene that has been linked to mammalian body size, insulin‐like growth factor‐1 (IGF‐1). We found high genetic differentiation at neutral molecular markers but convergence of allele frequencies at the IGF‐1 locus. Thus, we showed that while genetic drift has influenced the observed genetic structure of lynx at neutral molecular markers, natural selection has also played a role in the observed patterns of genetic diversity at the IGF‐1 locus of insular lynx.     

Parallel and nonparallel genome‐wide divergence among replicate population pairs of freshwater and anadromous Atlantic salmon

Little is known about the genetic basis differentiating resident and anadromous forms found in many salmonid species. Using a medium‐density SNP array, we documented genomic diversity and divergence at 2336 genetically mapped loci among three pairs of North American anadromous and freshwater Atlantic salmon populations. Our results show that across the genome, freshwater populations have lower diversity and a smaller proportion of private polymorphism relative to anadromous populations. Moreover, differentiation was more pronounced among freshwater than among anadromous populations at multiple spatial scales, suggesting a large effect of genetic drift in these isolated freshwater populations. Using nonhierarchical and hierarchical genome scans, we identified hundreds of markers spread across the genome that are potentially under divergent selection between anadromous and freshwater populations, but few outlier loci were repeatedly found in all three freshwater–anadromous comparisons. Similarly, a sliding window analysis revealed numerous regions of high divergence that were nonparallel among the three comparisons. These last results show little evidence for the parallel evolution of alleles selected for in freshwater populations, but suggest nonparallel adaptive divergence at many loci of small effects distributed through the genome. Overall, this study emphasizes the important role of genetic drift in driving genome‐wide reduction in diversity and divergence in freshwater Atlantic salmon populations and suggests a complex multigenic basis of adaptation to resident and anadromous strategies with little parallelism.     

Spatial and temporal patterns of neutral and adaptive genetic variation in the endangered African wild dog (Lycaon pictus)

Deciphering patterns of genetic variation within a species is essential for understanding population structure, local adaptation and differences in diversity between populations. Whilst neutrally evolving genetic markers can be used to elucidate demographic processes and genetic structure, they are not subject to selection and therefore are not informative about patterns of adaptive variation. As such, assessments of pertinent adaptive loci, such as the immunity genes of the major histocompatibility complex (MHC), are increasingly being incorporated into genetic studies. In this study, we combined neutral (microsatellite, mtDNA) and adaptive (MHC class II DLA‐DRB1 locus) markers to elucidate the factors influencing patterns of genetic variation in the African wild dog (Lycaon pictus); an endangered canid that has suffered extensive declines in distribution and abundance. Our genetic analyses found all extant wild dog populations to be relatively small (N_e < 30). Furthermore, through coalescent modelling, we detected a genetic signature of a recent and substantial demographic decline, which correlates with human expansion, but contrasts with findings in some other African mammals. We found strong structuring of wild dog populations, indicating the negative influence of extensive habitat fragmentation and loss of gene flow between habitat patches. Across populations, we found that the spatial and temporal structure of microsatellite diversity and MHC diversity were correlated and strongly influenced by demographic stability and population size, indicating the effects of genetic drift in these small populations. Despite this correlation, we detected signatures of selection at the MHC, implying that selection has not been completely overwhelmed by genetic drift.     

Drift,selection and adaptive variation in small populations of a threatened rattlesnake

An important goal of conservation genetics is to determine if the viability of small populations is reduced by a loss of adaptive variation due to genetic drift. Here, we assessed the impact of drift and selection on direct measures of adaptive variation (toxin loci encoding venom proteins) in the eastern massasauga rattlesnake (Sistrurus catenatus), a threatened reptile that exists in small isolated populations. We estimated levels of individual polymorphism in 46 toxin loci and 1,467 control loci across 12 populations of this species, and compared the results with patterns of selection on the same loci following speciation of S. catenatus and its closest relative, the western massasauga (S. tergeminus). Multiple lines of evidence suggest that both drift and selection have had observable impacts on standing adaptive variation. In support of drift effects, we found little evidence for selection on toxin variation within populations and a significant positive relationship between current levels of adaptive variation and long‐ and short‐term estimates of effective population size. However, we also observed levels of directional selection on toxin loci among populations that are broadly similar to patterns predicted from interspecific selection analyses that pre‐date the effects of recent drift, and that functional variation in these loci persists despite small short‐term effective sizes. This suggests that much of the adaptive variation present in populations may represent an example of “drift debt,” a nonequilibrium state where present‐day levels of variation overestimate the amount of functional genetic diversity present in future populations.     

Major Histocompatibility Complex (MHC) class II sequence polymorphism in long-finned pilot whale (Globicephala melas) from the North Atlantic

Determining how intra-specific genetic diversity is apportioned among natural populations is essential for detecting local adaptation and identifying populations with inherently low levels of extant diversity which may become a conservation concern. Sequence polymorphism at two adaptive loci (MHC DRA and DQB) was investigated in long-finned pilot whales (Globicephala melas) from four regions in the North Atlantic and compared with previous data from New Zealand (South Pacific). Three alleles were resolved at each locus, with trans-species allele sharing and higher levels of non-synonymous to synonymous substitution, especially in the DQB locus. Overall nucleotide diversities of 0.49?±?0.38% and 4.60?±?2.39% were identified for the DRA and DQB loci, respectively, which are relatively low for MHC loci in the North Atlantic, but comparable to levels previously described in New Zealand (South Pacific). There were significant differences in allele frequencies within the North Atlantic and between the North Atlantic and New Zealand. Patterns of diversity and divergence are consistent with the long-term effects of balancing selection operating on the MHC loci, potentially mediated through the effects of host-parasite coevolution. Differences in allele frequency may reflect variation in pathogen communities, coupled with the effects of differential drift and gene flow.     

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.     

LANDSCAPE GENOMICS IN ATLANTIC SALMON (SALMO SALAR): SEARCHING FOR GENE–ENVIRONMENT INTERACTIONS DRIVING LOCAL ADAPTATION

A growing number of studies are examining the factors driving historical and contemporary evolution in wild populations. By combining surveys of genomic variation with a comprehensive assessment of environmental parameters, such studies can increase our understanding of the genomic and geographical extent of local adaptation in wild populations. We used a large‐scale landscape genomics approach to examine adaptive and neutral differentiation across 54 North American populations of Atlantic salmon representing seven previously defined genetically distinct regional groups. Over 5500 genome‐wide single nucleotide polymorphisms were genotyped in 641 individuals and 28 bulk assays of 25 pooled individuals each. Genome scans, linkage map, and 49 environmental variables were combined to conduct an innovative landscape genomic analysis. Our results provide valuable insight into the links between environmental variation and both neutral and potentially adaptive genetic divergence. In particular, we identified markers potentially under divergent selection, as well as associated selective environmental factors and biological functions with the observed adaptive divergence. Multivariate landscape genetic analysis revealed strong associations of both genetic and environmental structures. We found an enrichment of growth‐related functions among outlier markers. Climate (temperature–precipitation) and geological characteristics were significantly associated with both potentially adaptive and neutral genetic divergence and should be considered as candidate loci involved in adaptation at the regional scale in Atlantic salmon. Hence, this study significantly contributes to the improvement of tools used in modern conservation and management schemes of Atlantic salmon wild populations.     

AFLPs and Mitochondrial Haplotypes Reveal Local Adaptation to Extreme Thermal Environments in a Freshwater Gastropod

The way environmental variation shapes neutral and adaptive genetic variation in natural populations is a key issue in evolutionary biology. Genome scans allow the identification of the genetic basis of local adaptation without previous knowledge of genetic variation or traits under selection. Candidate loci for divergent adaptation are expected to show higher F_ST than neutral loci influenced solely by random genetic drift, migration and mutation. The comparison of spatial patterns of neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection among populations living in contrasting environments. Using the gastropod Radix balthica as a system, we analyzed 376 AFLP markers and 25 mtDNA COI haplotypes for candidate loci and associations with local adaptation among contrasting thermal environments in Lake Mývatn, a volcanic lake in northern Iceland. We found that 2% of the analysed AFLP markers were under directional selection and 12% of the mitochondrial haplotypes correlated with differing thermal habitats. The genetic networks were concordant for AFLP markers and mitochondrial haplotypes, depicting distinct topologies at neutral and candidate loci. Neutral topologies were characterized by intense gene flow revealed by dense nets with edges connecting contrasting thermal habitats, whereas the connections at candidate loci were mostly restricted to populations within each thermal habitat and the number of edges decreased with temperature. Our results suggest microgeographic adaptation within Lake Mývatn and highlight the utility of genome scans in detecting adaptive divergence.     

Population genomics of the killer whale indicates ecotype evolution in sympatry involving both selection and drift

The evolution of diversity in the marine ecosystem is poorly understood, given the relatively high potential for connectivity, especially for highly mobile species such as whales and dolphins. The killer whale (Orcinus orca) has a worldwide distribution, and individual social groups travel over a wide geographic range. Even so, regional populations have been shown to be genetically differentiated, including among different foraging specialists (ecotypes) in sympatry. Given the strong matrifocal social structure of this species together with strong resource specializations, understanding the process of differentiation will require an understanding of the relative importance of both genetic drift and local adaptation. Here we provide a high‐resolution analysis based on nuclear single‐nucleotide polymorphic markers and inference about differentiation at both neutral loci and those potentially under selection. We find that all population comparisons, within or among foraging ecotypes, show significant differentiation, including populations in parapatry and sympatry. Loci putatively under selection show a different pattern of structure compared to neutral loci and are associated with gene ontology terms reflecting physiologically relevant functions (e.g. related to digestion). The pattern of differentiation for one ecotype in the North Pacific suggests local adaptation and shows some fixed differences among sympatric ecotypes. We suggest that differential habitat use and resource specializations have promoted sufficient isolation to allow differential evolution at neutral and functional loci, but that the process is recent and dependent on both selection and drift.     

Adaptive and neutral genetic differentiation among Scottish and endangered Irish red grouse (<Emphasis Type="Italic">Lagopus lagopus scotica</Emphasis>)

Studying patterns of intra-specific genetic variation among populations allows for a better understanding of population structure and local adaptation. However, those patterns may differ according to the genetic markers applied, as neutral genetic markers reflect demographic processes and random genetic drift, whereas adaptive markers also carry the footprint of selection. In combination, neutral and adaptive genetic markers permit to assess the relative roles of drift and selection in shaping population structure. Among the best understood adaptive genetic loci are the genes of the major histocompatibility complex (MHC). We here study variation and differentiation at neutral SNP markers and MHC class II genes in red grouse (Lagopus lagopus scotica) from Ireland and Scotland. Irish red grouse populations are fragmented and drastically declining, but red grouse are abundant in Scotland. We find evidence for positive selection acting on the MHC genes and variation in MHC gene copy numbers among Irish individuals. Furthermore, there was significant population differentiation among red grouse from Ireland and Scotland at the neutral SNP markers (F_ST = 0.084) and the MHC-BLB genes (F_ST: BLB1 = 0.116, BLB2 = 0.090, BLB3 = 0.104). Differentiation at the MHC-BLB1 was significantly higher than at the neutral SNP markers, suggesting that selection plays an important role in shaping MHC variation, in addition to genetic drift. We speculate that the observed differentiation pattern might be due to local adaptation to different parasite regimes. These findings have strong conservation implications and we advise against the introduction of Scottish red grouse to supplement Irish populations.     

Conserving adaptive genetic diversity in dynamic landscapes

Genetic variation supplies the raw material for adaptation, evolution and survival of populations and has therefore been a key focus of conservation biology ever since its foundation (Soulé 1985). In previous decades, the neutral component of genetic diversity (generated by mutation and shaped by drift) has been the subject of intense scientific research, fuelled by the increasing availability of molecular markers. On the other hand, the adaptive component of genetic diversity, which is shaped by the action of natural selection, has long remained elusive and difficult to assess, especially at small spatial or temporal scales (Ouborg et al. 2010). Fortunately, new technological and methodological developments now make it possible to identify loci in the genome that are influenced by selection, and thus to get a more complete view of genetic diversity. One article featured in this issue of Molecular Ecology is a good example of this recent breakthrough. Richter-Boix et al. (2011) examined a network of moor frog populations breeding in contrasting habitats in order to understand how landscape features influence patterns of genetic variation. They combined information from both neutral markers and loci putatively under selection to quantify the relative roles of selection and isolation in the evolution of fine-scale local adaptations in these populations. This study nicely illustrates how data on polymorphisms of neutral and adaptive loci can now be judiciously synthesized to help identify the best strategies for preserving adaptive variation, and more generally to enlighten conservation and population-management plans.     

Genetic isolation of fragmented populations is exacerbated by drift and selection

Reduced genetic variation at marker loci in small populations has been well documented, whereas the relationship between quantitative genetic variation and population size has attracted little empirical investigation. Here we demonstrate that both neutral and quantitative genetic variation are reduced in small populations of a fragmented plant metapopulation, and that both drift and selective change are enhanced in small populations. Measures of neutral genetic differentiation (F(ST)) and quantitative genetic differentiation (Q(ST)) in two traits were higher among small demes, and Q(ST) between small populations exceeded that expected from drift alone. This suggests that fragmented populations experience both enhanced genetic drift and divergent selection on phenotypic traits, and that drift affects variation in both neutral markers and quantitative traits. These results highlight the need to integrate natural selection into conservation genetic theory, and suggests that small populations may represent reservoirs of genetic variation adaptive within a wide range of environments.     

Outlier SNP markers reveal fine‐scale genetic structuring across European hake populations (Merluccius merluccius)

Shallow population structure is generally reported for most marine fish and explained as a consequence of high dispersal, connectivity and large population size. Targeted gene analyses and more recently genome‐wide studies have challenged such view, suggesting that adaptive divergence might occur even when neutral markers provide genetic homogeneity across populations. Here, 381 SNPs located in transcribed regions were used to assess large‐ and fine‐scale population structure in the European hake (Merluccius merluccius), a widely distributed demersal species of high priority for the European fishery. Analysis of 850 individuals from 19 locations across the entire distribution range showed evidence for several outlier loci, with significantly higher resolving power. While 299 putatively neutral SNPs confirmed the genetic break between basins (F_CT = 0.016) and weak differentiation within basins, outlier loci revealed a dramatic divergence between Atlantic and Mediterranean populations (F_CT range 0.275–0.705) and fine‐scale significant population structure. Outlier loci separated North Sea and Northern Portugal populations from all other Atlantic samples and revealed a strong differentiation among Western, Central and Eastern Mediterranean geographical samples. Significant correlation of allele frequencies at outlier loci with seawater surface temperature and salinity supported the hypothesis that populations might be adapted to local conditions. Such evidence highlights the importance of integrating information from neutral and adaptive evolutionary patterns towards a better assessment of genetic diversity. Accordingly, the generated outlier SNP data could be used for tackling illegal practices in hake fishing and commercialization as well as to develop explicit spatial models for defining management units and stock boundaries.     

Population genomic evidence for adaptive differentiation in the Baltic Sea herring

Detecting and estimating the degree of genetic differentiation among populations of highly mobile marine fish having pelagic larval stages is challenging because their effective population sizes can be large, and thus, little genetic drift and differentiation is expected in neutral genomic sites. However, genomic sites subject to directional selection stemming from variation in local environmental conditions can still show substantial genetic differentiation, yet these signatures can be hard to detect with low‐throughput approaches. Using a pooled RAD‐seq approach, we investigated genomewide patterns of genetic variability and differentiation within and among 20 populations of Atlantic herring in the Baltic Sea (and adjacent Atlantic sites), where previous low‐throughput studies and/or studies based on few populations have found limited evidence for genetic differentiation. Stringent quality control was applied in the filtering of 1 791 254 SNPs, resulting in a final data set of 68 182 polymorphic loci. Clear differentiation was identified between Atlantic and Baltic populations in many genomic sites, while differentiation within the Baltic Sea area was weaker and geographically less structured. However, outlier analyses – whether including all populations or only those within the Baltic Sea – uncovered hundreds of directionally selected loci in which variability was associated with either salinity, temperature or both. Hence, our results support the view that although the degree of genetic differentiation among Baltic Sea herring populations is low, there are many genomic regions showing elevated divergence, apparently as a response to temperature‐ and salinity‐related natural selection. As such, the results add to the increasing evidence of local adaptation in highly mobile marine organisms, and those in the young Baltic Sea in particular.     

Trait variation in response to varying winter temperatures,diversity patterns and signatures of selection along the latitudinal distribution of the widespread grassland plant Arrhenatherum elatius

Across Europe, genetic diversity can be expected to decline toward the North because of stochastic and selective effects which may imply diminished phenotypic variation and less potential for future genetic adaptations to environmental change. Understanding such latitudinal patterns can aid provenance selection for breeding or assisted migration approaches. In an experiment simulating different winter temperatures, we assessed quantitative trait variation, genetic diversity, and differentiation for natural populations of the grass Arrhenatherum elatius originating from a large latitudinal gradient. In general, populations from the North grew smaller and had a lower flowering probability. Toward the North, the absolute plastic response to the different winter conditions as well as heritability for biomass production significantly declined. Genetic differentiation in plant height and probability of flowering were very strong and significantly higher than under neutral expectations derived from SNP data, suggesting adaptive differentiation. Differentiation in biomass production did not exceed but mirrored patterns for neutral genetic differentiation, suggesting that migration‐related processes caused the observed clinal trait variation. Our results demonstrate that genetic diversity and trait differentiation patterns for A. elatius along a latitudinal gradient are likely shaped by both local selection and genetic drift.     

Isolation by adaptation in Neochlamisus leaf beetles: host-related selection promotes neutral genomic divergence

This study tests how divergent natural selection promotes genomic differentiation during ecological speciation. Specifically, we use adaptive ecological divergence (here, population divergence in host plant use and preference) as a proxy for selection strength and evaluate the correlation between levels of adaptive and genetic differentiation across pairwise population comparisons. Positive correlations would reveal the pattern predicted by our hypothesis, that of 'isolation by adaptation' (IBA). Notably, IBA is predicted not only for selected loci but also for neutral loci. This may reflect the effects of divergent selection on neutral loci that are 'loosely linked' to divergently selected loci or on geneflow restriction that facilitates genetic drift at all loci, including neutral loci that are completely unlinked to those evolving under divergent selection. Here, we evaluate IBA in maple- and willow-associated populations of Neochlamisus bebbianae leaf beetles. To do so, we collected host preference data to construct adaptive divergence indices and used AFLPs (amplified fragment length polymorphisms) and mitochondrial sequences to quantify genetic differentiation. Partial Mantel tests showed significant IBA in 'pooled' analyses of putatively neutral and of putatively selected ('outlier') AFLP loci. This pattern was also recovered in 12% of 'locus-specific' analyses that separately evaluated genetic differentiation at individual neutral loci. These results provided evidence for widespread effects of selection on neutral genomic divergence. Our collective findings indicate that host-related selection may play important roles in the population genomic differentiation of both neutral and selected gene regions in herbivorous insects.     

Distinct evolutionary strategies of human leucocyte antigen loci in pathogen-rich environments

Human leucocyte antigen (HLA) loci have a complex evolution where both stochastic (e.g. genetic drift) and deterministic (natural selection) forces are involved. Owing to their extraordinary level of polymorphism, HLA genes are useful markers for reconstructing human settlement history. However, HLA variation often deviates significantly from neutral expectations towards an excess of genetic diversity. Because HLA molecules play a crucial role in immunity, this observation is generally explained by pathogen-driven-balancing selection (PDBS). In this study, we investigate the PDBS model by analysing HLA allelic diversity on a large database of 535 populations in relation to pathogen richness. Our results confirm that geographical distances are excellent predictors of HLA genetic differentiation worldwide. We also find a significant positive correlation between genetic diversity and pathogen richness at two HLA class I loci (HLA-A and -B), as predicted by PDBS, and a significant negative correlation at one HLA class II locus (HLA-DQB1). Although these effects are weak, as shown by a loss of significance when populations submitted to rapid genetic drift are removed from the analysis, the inverse relationship between genetic diversity and pathogen richness at different loci indicates that HLA genes have adopted distinct evolutionary strategies to provide immune protection in pathogen-rich environments.