Contrasting patterns of body shape and neutral genetic divergence in marine and lake populations of threespine sticklebacks

Comparisons of neutral marker and quantitative trait divergence can provide important insights into the relative roles of natural selection and neutral genetic drift in population differentiation. We investigated phenotypic and genetic differentiation among Fennoscandian threespine stickleback (Gasterosteus aculeatus) populations, and found that the highest degree of differentiation occurred between sea and freshwater habitats. Within habitats, morphological divergence was highest among the different freshwater populations. Pairwise phenotypic and neutral genetic distances among populations were positively correlated, suggesting that genetic drift may have contributed to the morphological differentiation among habitats. On the other hand, the degree of phenotypic differentiation (PST) clearly surpassed the neutral expectation set by FST, suggesting a predominant role for natural selection over genetic drift as an explanation for the observed differentiation. However, separate PST/FST comparisons by habitats revealed that body shape divergence between lake and marine populations, and even among marine populations, can be strongly influenced by natural selection. On the other hand, genetic drift can play an important role in the differentiation among lake populations.     

Patterns of male sterility in a grasshopper hybrid zone imply accumulation of hybrid incompatibilities without selection

It is now widely accepted that post-zygotic reproductive isolation is the result of negative epistatic interactions between derived alleles fixed independently at different loci in diverging populations (the Dobzhansky-Muller model). What is less clear is the nature of the loci involved and whether the derived alleles increase in frequency through genetic drift, or as a result of natural or sexual selection. If incompatible alleles are fixed by selection, transient polymorphisms will be rare and clines for these alleles will be steep where divergent populations meet. If they evolve by drift, populations are expected to harbour substantial genetic variation in compatibility and alleles will introgress across hybrid zones once they recombine onto a genetic background with which they are compatible. Here we show that variation in male sterility in a naturally occurring Chorthippus parallelus grasshopper hybrid zone conforms to the neutral expectations. Asymmetrical clines for male sterility have long tails of introgression and populations distant from the zone centre show significant genetic variation for compatibility. Our data contrast with recent observations on 'speciation genes' that have diverged as a result of strong natural selection.     

Comparing the effects of genetic drift and fluctuating selection on genotype frequency changes in the scarlet tiger moth

One of the recurring arguments in evolutionary biology is whether evolution occurs principally through natural selection or through neutral processes such as genetic drift. A 60-year-long time series of changes in the genotype frequency of a colour polymorphism of the scarlet tiger moth, Callimorpha dominula, was used to compare the relative effects of genetic drift and variable natural selection. The analysis showed that most of the variation in frequency was the result of genetic drift. In addition, although selection was acting, mean fitness barely increased. This supports the 'Red Queen's hypothesis' that long-term improvements in fitness may not occur, because populations have to keep pace with changes in the environment.     

Microevolution in island forms: the roles of drift and directional selection in morphological divergence of a passerine bird

Abstract.— Theory predicts that in small isolated populations random genetic drift can lead to phenotypic divergence; however this prediction has rarely been tested quantitatively in natural populations. Here we utilize natural repeated island colonization events by members of the avian species complex, Zosterops lateralis , to assess whether or not genetic drift alone is an adequate explanation for the observed patterns of microevolutionary divergence in morphology. Morphological and molecular genetic characteristics of island and mainland populations are compared to test three predictions of drift theory: (1) that the pattern of morphological change is idiosyncratic to each island; (2) that there is concordance between morphological and neutral genetic shifts across island populations; and (3) for populations whose time of colonization is known, that the rate of morphological change is sufficiently slow to be accounted for solely by genetic drift. Our results are not consistent with these predictions. First, the direction of size shifts was consistently towards larger size, suggesting the action of a nonrandom process. Second, patterns of morphological divergence among recently colonized populations showed little concordance with divergence in neutral genetic characters. Third, rate tests of morphological change showed that effective population sizes were not small enough for random processes alone to account for the magnitude of microevolutionary change. Altogether, these three lines of evidence suggest that drift alone is not an adequate explanation of morphological differentiation in recently colonized island Zosterops and therefore we suggest that the observed microevolutionary changes are largely a result of directional natural selection.     

Genetic Drift of HIV Populations in Culture

Populations of Human Immunodeficiency Virus type 1 (HIV-1) undergo a surprisingly large amount of genetic drift in infected patients despite very large population sizes, which are predicted to be mostly deterministic. Several models have been proposed to explain this phenomenon, but all of them implicitly assume that the process of virus replication itself does not contribute to genetic drift. We developed an assay to measure the amount of genetic drift for HIV populations replicating in cell culture. The assay relies on creation of HIV populations of known size and measurements of variation in frequency of a neutral allele. Using this assay, we show that HIV undergoes approximately ten times more genetic drift than would be expected from its population size, which we defined as the number of infected cells in the culture. We showed that a large portion of the increase in genetic drift is due to non-synchronous infection of target cells. When infections are synchronized, genetic drift for the virus is only 3-fold higher than expected from its population size. Thus, the stochastic nature of biological processes involved in viral replication contributes to increased genetic drift in HIV populations. We propose that appreciation of these effects will allow better understanding of the evolutionary forces acting on HIV in infected patients.     

The Functional Gene Diversity in Natural Populations over Postglacial Areas: The Shaping Mechanisms Behind Genetic Composition of Longnose Dace (Rhinichthys cataractae) in Northeastern North America

The diversity of functional genes and the related processes are important issues for conservation biology. This is especially relevant for populations that have suffered from demographic reduction as a consequence of the processes of postglacial colonization. In this perspective, the aims of the present study are (1) to quantify the genetic diversity of functional genes and (2) to disentangle the long- and short-term effects of natural selection that shapes genetic diversity from those of drift, mutation, and allopatric fragmentation. This research was conducted using an extensive genetic polymorphism analysis of populations of longnose dace (Rhinichthys cataractae) living over an area once covered by Pleistocene glaciations. The sequence and diversity of one exon of three genes (MHC IIβ, growth hormone, and trypsin) were jointly analyzed with non-coding nuclear loci from 27 populations; these populations were sampled over four major basins of northeastern North America. The survey revealed a surprisingly low allelic richness, especially for the MHC gene, considering the number of individuals and populations sampled. The results suggest that there is a complex mixture of different evolutionary processes shaping the level of polymorphism among longnose dace. While our study underlines the importance of the short-term effects of neutral processes and the major impact of post-glacial colonization on gene diversity, locally dependent balancing selection was detected on MHC. From this perspective, our results support an understanding of the importance of drift on functional gene diversity but also highlight the transient effects of natural selection on allelic composition, even in populations that show drastic reduction of genetic diversity.     

Latitudinal divergence in a widespread amphibian: Contrasting patterns of neutral and adaptive genomic variation

Stochastic effects from demographic processes and selection are expected to shape the distribution of genetic variation in spatially heterogeneous environments. As the amount of genetic variation is central for long‐term persistence of populations, understanding how these processes affect variation over large‐scale geographical gradients is pivotal. We investigated the distribution of neutral and putatively adaptive genetic variation, and reconstructed demographic history in the moor frog (Rana arvalis) using 136 individuals from 15 populations along a 1,700‐km latitudinal gradient from northern Germany to northern Sweden. Using double digest restriction‐site associated DNA sequencing we obtained 27,590 single nucleotide polymorphisms (SNPs), and identified differentiation outliers and SNPs associated with growing season length. The populations grouped into a southern and a northern cluster, representing two phylogeographical lineages from different post‐glacial colonization routes. Hybrid index estimation and demographic model selection showed strong support for a southern and northern lineage and evidence of gene flow between regions located on each side of a contact zone. However, patterns of past gene flow over the contact zone differed between neutral and putatively adaptive SNPs. While neutral nucleotide diversity was higher along the southern than the northern part of the gradient, nucleotide diversity in differentiation outliers showed the opposite pattern, suggesting differences in the relative strength of selection and drift along the gradient. Variation associated with growing season length decreased with latitude along the southern part of the gradient, but not along the northern part where variation was lower, suggesting stronger climate‐mediated selection in the north. Outlier SNPs included loci involved in immunity and developmental processes.     

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.     

Simulated climate change provokes rapid genetic change in the Mediterranean shrub Fumana thymifolia

Rapid climate change will impose strong directional selection pressures on natural plant populations. Climate-linked genetic variation in natural populations indicates that an evolutionary response is possible. We investigated such a response by comparing individuals subjected to elevated drought and warming treatments with individuals establishing in an unmanipulated climate within the same population. We report that reduction in seedling establishment in response to climate manipulations is nonrandom and results from the selection pressure imposed by artificially warmed and droughted conditions. When compared against control samples, high single-locus genetic divergence occurred in drought and warming treatment samples, with genetic differentiation up to 37 times higher than background (mean neutral locus) genetic differentiation. These loci violate assumptions of selective neutrality, indicating the signature of natural selection by drought. Our results demonstrate that rapid evolution in response to climate change may be widespread in natural populations, based on genetic variation already present within the population.     

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.     

Were neandertal and modern human cranial differences produced by natural selection or genetic drift?

Most evolutionary explanations for cranial differences between Neandertals and modern humans emphasize adaptation by natural selection. Features of the crania of Neandertals could be adaptations to the glacial climate of Pleistocene Europe or to the high mechanical strains produced by habitually using the front teeth as tools, while those of modern humans could be adaptations for articulate speech production. A few researchers have proposed non-adaptive explanations. These stress that isolation between Neandertal and modern human populations would have lead to cranial diversification by genetic drift (chance changes in the frequencies of alleles at genetic loci contributing to variation in cranial morphology). Here we use a variety of statistical tests founded on explicit predictions from quantitative- and population-genetic theory to show that genetic drift can explain cranial differences between Neandertals and modern humans. These tests are based on thirty-seven standard cranial measurements from a sample of 2524 modern humans from 30 populations and 20 Neandertal fossils. As a further test, we compare our results for modern human cranial measurements with those for a genetic dataset consisting of 377 microsatellites typed for a sample of 1056 modern humans from 52 populations. We conclude that rather than requiring special adaptive accounts, Neandertal and modern human crania may simply represent two outcomes from a vast space of random evolutionary possibilities.     

Morphological Clines and Weak Drift along an Urbanization Gradient in the Butterfly,Pieris rapae

Urban areas are increasing globally, providing opportunities for biodiversity researchers to study the process in which species become established in novel, highly disturbed habitats. This ecological process can be understood through analyses of morphological and genetic variation, which can shed light on patterns of neutral and adaptive evolution. Previous studies have shown that urban populations often diverge genetically from non-urban source populations. This could occur due to neutral genetic drift, but an alternative is that selection could lead to allele frequency changes in urban populations. The development of genome scan methods provides an opportunity to investigate these outcomes from samples of genetic variation taken along an urbanization gradient. Here we examine morphological variation in wing size and diversity at neutral amplified fragment length polymorphisms in the butterfly Pieris rapae L. (Lepidoptera, Pieridae) sampled from the center to the periphery of Marseille. We utilize established and novel environmental correlation approaches to scan genetic variation for evidence of selection. We find significant morphological differences in urban populations, as well as weak genetic structure and decreased genetic diversity in urban versus non-urban sites. However, environmental correlation tests provide little support for selection in our dataset. Our comparison of different methods and allele frequency clines suggests that loci identified as significant are false positives. Although there is some indication that selection may be acting on wing size in urban butterflies, genetic analyses suggest P. rapae are undergoing neutral drift.     

Evolution of adaptive diversity and genetic connectivity in Arctic charr (Salvelinus alpinus) in Iceland

The ecological theory of adaptive radiation predicts that the evolution of phenotypic diversity within species is generated by divergent natural selection arising from different environments and competition between species. Genetic connectivity among populations is likely also to have an important role in both the origin and maintenance of adaptive genetic diversity. Our goal was to evaluate the potential roles of genetic connectivity and natural selection in the maintenance of adaptive phenotypic differences among morphs of Arctic charr, Salvelinus alpinus, in Iceland. At a large spatial scale, we tested the predictive power of geographic structure and phenotypic variation for patterns of neutral genetic variation among populations throughout Iceland. At a smaller scale, we evaluated the genetic differentiation between two morphs in Lake Thingvallavatn relative to historically explicit, coalescent-based null models of the evolutionary history of these lineages. At the large spatial scale, populations are highly differentiated, but weakly structured, both geographically and with respect to patterns of phenotypic variation. At the intralacustrine scale, we observe modest genetic differentiation between two morphs, but this level of differentiation is nonetheless consistent with strong reproductive isolation throughout the Holocene. Rather than a result of the homogenizing effect of gene flow in a system at migration-drift equilibrium, the modest level of genetic differentiation could equally be a result of slow neutral divergence by drift in large populations. We conclude that contemporary and recent patterns of restricted gene flow have been highly conducive to the evolution and maintenance of adaptive genetic variation in Icelandic Arctic charr.     

Drift versus selection as drivers of phenotypic divergence at small spatial scales: The case of Belgjarskógur threespine stickleback

Divergence in phenotypic traits is facilitated by a combination of natural selection, phenotypic plasticity, gene flow, and genetic drift, whereby the role of drift is expected to be particularly important in small and isolated populations. Separating the components of phenotypic divergence is notoriously difficult, particularly for multivariate phenotypes. Here, we assessed phenotypic divergence of threespine stickleback (Gasterosteus aculeatus) across 19 semi‐interconnected ponds within a small geographic region (~7.5 km²) using comparisons of multivariate phenotypic divergence (PST), neutral genetic (FST), and environmental (EST) variation. We found phenotypic divergence across the ponds in a suite of functionally relevant phenotypic traits, including feeding, defense, and swimming traits, and body shape (geometric morphometric). Comparisons of PSTs with FSTs suggest that phenotypic divergence is predominantly driven by neutral processes or stabilizing selection, whereas phenotypic divergence in defensive traits is in accordance with divergent selection. Comparisons of population pairwise PSTs with ESTs suggest that phenotypic divergence in swimming traits is correlated with prey availability, whereas there were no clear associations between phenotypic divergence and environmental difference in the other phenotypic groups. Overall, our results suggest that phenotypic divergence of these small populations at small geographic scales is largely driven by neutral processes (gene flow, drift), although environmental determinants (natural selection or phenotypic plasticity) may play a role.     

Diversifying selection drives parallel evolution of gill raker number and body size along the speciation continuum of European whitefish

Adaptive radiation is the evolution of ecological and phenotypical diversity. It arises via ecological opportunity that promotes the exploration of underutilized or novel niches mediating specialization and reproductive isolation. The assumed precondition for rapid local adaptation is diversifying natural selection, but random genetic drift could also be a major driver of this process. We used 27 populations of European whitefish (Coregonus lavaretus) from nine lakes distributed in three neighboring subarctic watercourses in northern Fennoscandia as a model to test the importance of random drift versus diversifying natural selection for parallel evolution of adaptive phenotypic traits. We contrasted variation for two key adaptive phenotypic traits correlated with resource utilization of polymorphic fish; the number of gill rakers and the total length of fish, with the posterior distribution of neutral genetic differentiation from 13 microsatellite loci, to test whether the observed phenotypic divergence could be achieved by random genetic drift alone. Our results show that both traits have been under diversifying selection and that the evolution of these morphs has been driven by isolation through habitat adaptations. We conclude that diversifying selection acting on gill raker number and body size has played a significant role in the ongoing adaptive radiation of European whitefish morphs in this region.     

A new method to uncover signatures of divergent and stabilizing selection in quantitative traits

While it is well understood that the pace of evolution depends on the interplay between natural selection, random genetic drift, mutation, and gene flow, it is not always easy to disentangle the relative roles of these factors with data from natural populations. One popular approach to infer whether the observed degree of population differentiation has been influenced by local adaptation is the comparison of neutral marker gene differentiation (as reflected in FST) and quantitative trait divergence (as reflected in QST). However, this method may lead to compromised statistical power, because FST and QST are summary statistics which neglect information on specific pairs of populations, and because current multivariate tests of neutrality involve an averaging procedure over the traits. Further, most FST-QST comparisons actually replace QST by its expectation over the evolutionary process and are thus theoretically flawed. To overcome these caveats, we derived the statistical distribution of population means generated by random genetic drift and used the probability density of this distribution to test whether the observed pattern could be generated by drift alone. We show that our method can differentiate between genetic drift and selection as a cause of population differentiation even in cases with FST=QST and demonstrate with simulated data that it disentangles drift from selection more accurately than conventional FST-QST tests especially when data sets are small.     

Inbreeding depression and low between-population heterosis in recently diverged experimental populations of a selfing species

In fragmented populations, genetic drift and selection reduce genetic diversity, which in turn results in a loss of fitness or in a loss of evolvability. Genetic rescue, that is, controlled input of diversity from distant populations, may restore evolutionary potential, whereas outbreeding depression might counteract the positive effect of this strategy. We carried out self-pollination and crosses within and between populations in an experimental subdivided population of a selfing species, Triticum aestivum L., to estimate the magnitude of these two phenomena. Surprisingly, for a self-fertilizing species, we found significant inbreeding depression within each population for four of the six traits studied, indicating that mildly deleterious mutations were still segregating in these populations. The progeny of within- and between-population crosses was very similar, indicating low between-population heterosis and little outbreeding depression. We conclude that relatively large population effective sizes prevented fixation of a high genetic load and that local adaptation was limited in these recently diverged populations. The kinship coefficient estimated between the parents using 20 neutral markers was a poor predictor of the progeny phenotypic values, indicating that there was a weak link between neutral diversity and genes controlling fitness-related traits. These results show that when assessing the viability of natural populations and the need for genetic rescue, the use of neutral markers should be complemented with information about the presence of local adaptation in the subdivided population.     

The influence of stochastic and selective forces in the population divergence of female colour polymorphism in damselflies of the genus Ischnura

Disentangling the relative importance and potential interactions of selection and genetic drift in driving phenotypic divergence of species is a classical research topic in population genetics and evolutionary biology. Here, we evaluate the role of stochastic and selective forces on population divergence of a colour polymorphism in seven damselfly species of the genus Ischnura, with a particular focus on I. elegans and I. graellsii. Colour-morph frequencies in Spanish I. elegans populations varied greatly, even at a local scale, whereas more similar frequencies were found among populations in eastern Europe. In contrast, I. graellsii and the other five Ischnura species showed little variation in colour-morph frequencies between populations. F(ST)-outlier analyses revealed that the colour locus deviated strongly from neutral expectations in Spanish populations of I. elegans, contrasting the pattern found in eastern European populations, and in I. graellsii, where no such discrepancy between morph divergence and neutral divergence could be detected. This suggests that divergent selection has been operating on the colour locus in Spanish populations of I. elegans, whereas processes such as genetic drift, possibly in combination with other forms of selection (such as negative frequency-dependent selection), appear to have been present in other regions, such as eastern Europe. Overall, the results indicate that both selective and stochastic processes operate on these colour polymorphisms, and suggest that the relative importance of factors varies between geographical regions.     

Novelty-seeking DRD4 polymorphisms are associated with human migration distance out-of-Africa after controlling for neutral population gene structure

Numerous lines of evidence suggest that Homo sapiens evolved as a distinct species in Africa by 150,000 years before the present (BP) and began major migrations out-of-Africa ～50,000 BP. By 20,000 BP, our species had effectively colonized the entire Old World, and by 12,000 BP H. sapiens had a global distribution. We propose that this rapid migration into new habitats selected for individuals with low reactivity to novel stressors. Certain dopamine receptor D4 (DRD4) polymorphisms are associated with low neuronal reactivity and increased exploratory behavior, novelty seeking, and risk taking, collectively considered novelty-seeking trait (NS). One previous report (Chen et al.: Evol Hum Behav 20 (1999) 309-324) demonstrated a correlation between migratory distance and the seven-repeat (7R) VNTR DRD4 allele at exon 3 for human populations. This study, however, failed to account for neutral genetic processes (drift and admixture) that might create such a correlation in the absence of natural selection. Furthermore, additional loci surrounding DRD4 are now recognized to influence NS. Herein we account for neutral genetic structure by modeling the nonindependence of neutral allele frequencies between human populations. We retest the DRD4 exon 3 alleles, and also test two other loci near DRD4 that are associated with NS. We conclude there is an association between migratory distance and DRD4 exon 3 2R and 7R alleles that cannot be accounted for by neutral genetic processes alone.     

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.