Natural Selection and Neutral Evolution Jointly Drive Population Divergence between Alpine and Lowland Ecotypes of the Allopolyploid Plant Anemone multifida (Ranunculaceae)

Population differentiation can be driven in large part by natural selection, but selectively neutral evolution can play a prominent role in shaping patters of population divergence. The decomposition of the evolutionary history of populations into the relative effects of natural selection and selectively neutral evolution enables an understanding of the causes of population divergence and adaptation. In this study, we examined heterogeneous genomic divergence between alpine and lowland ecotypes of the allopolyploid plant, Anemone multifida. Using peak height and dominant AFLP data, we quantified population differentiation at non-outlier (neutral) and outlier loci to determine the potential contribution of natural selection and selectively neutral evolution to population divergence. We found 13 candidate loci, corresponding to 2.7% of loci, with signatures of divergent natural selection between alpine and lowland populations and between alpine populations (Fst  = 0.074–0.445 at outlier loci), but neutral population differentiation was also evident between alpine populations (F_ST  = 0.041–0.095 at neutral loci). By examining population structure at both neutral and outlier loci, we determined that the combined effects of selection and neutral evolution are associated with the divergence of alpine populations, which may be linked to extreme abiotic conditions and isolation between alpine sites. The presence of outlier levels of genetic variation in structured populations underscores the importance of separately analyzing neutral and outlier loci to infer the relative role of divergent natural selection and neutral evolution in population divergence.     

Ethnic genetics: relation between adaptive and neutral genetic differentiation of an ethnos

It is very important to be able to distinguish between selectively significant genetic variability and selectively-neutral one for quantitative analysis of genetic differentiation in human (and any other) populations. The key to the problem is to determine a start-point for detection of neutral genetic variability, which will help to establish alignment of adaptive and neutral forces operating in genetic differentiation. The purpose of this work is to adduce proofs in favour of mean Fst value for a sample of gene loci as of the start-point for measuring neutrality level of genetic differentiation. These proofs came from various demographic and onomastic characteristics of ethnic groups as well as from genetic chronology of ethnic history which is in good concordance with actual historical chronology of ethnic groups. Once the start-point for testing neutrality is determined, it becomes possible to reveal the selective pressure to which various human genes are undergone and to elucidate adaptive structure of mankind's genetic pool. It was shown that only 15 alleles from 49 belonging to 20 polymorphic loci can be considered selectively neutral.     

Discrepancies in population differentiation at microsatellites, mitochondrial DNA and plumage colour in the pied flycatcher--inferring evolutionary processes

Genetic differentiation between three populations of the pied flycatcher Ficedula hypoleuca (Norway, Czech Republic and Spain, respectively) was investigated at microsatellite loci and mitochondrial DNA (mtDNA) sequences and compared with the pattern of differentiation of male plumage colour. The Czech population lives sympatrically with the closely related collared flycatcher (F. albicollis) whereas the other two are allopatric. Allopatric populations are on average more conspicuously coloured than sympatric ones, a pattern that has been explained by sexual selection for conspicuous colour in allopatry and a character displacement on breeding plumage colour in sympatry that reduces the rate of hybridization with the collared flycatcher. The Czech population was genetically indistinguishable from the Norwegian population at microsatellite loci and mtDNA sequences. Recent isolation and/or gene flow may explain the lack of genetic differentiation. Accordingly, different selection on plumage colour in the two populations is either sufficiently strong so that gene flow has little impact on the pattern of colour variation, or differentiation of plumage colour occurred so recently that the (presumably) neutral, fast evolving markers employed here are unable to reflect the differentiation. Genetically, the Spanish population was significantly differentiated from the other populations, but the divergence was much more pronounced at mtDNA compared to microsatellites. This may reflect increased rate of differentiation by genetic drift at the mitochondrial, compared with the nuclear genome, caused by the smaller effective population size of the former genome. In accordance with this interpretation, a genetic pattern consistent with effects of small population size in the Spanish population (genetic drift and inbreeding) were also apparent at the microsatellites, namely reduced allelic diversity and heterozygous deficiency.     

Differentiation in neutral genes and a candidate gene in the pied flycatcher: using biological archives to track global climate change

Global climate change is one of the major driving forces for adaptive shifts in migration and breeding phenology and possibly impacts demographic changes if a species fails to adapt sufficiently. In Western Europe, pied flycatchers (Ficedula hypoleuca) have insufficiently adapted their breeding phenology to the ongoing advance of food peaks within their breeding area and consequently suffered local population declines. We address the question whether this population decline led to a loss of genetic variation, using two neutral marker sets (mitochondrial control region and microsatellites), and one potentially selectively non‐neutral marker (avian Clock gene). We report temporal changes in genetic diversity in extant populations and biological archives over more than a century, using samples from sites differing in the extent of climate change. Comparing genetic differentiation over this period revealed that only the recent Dutch population, which underwent population declines, showed slightly lower genetic variation than the historic Dutch population. As that loss of variation was only moderate and not observed in all markers, current gene flow across Western and Central European populations might have compensated local loss of variation over the last decades. A comparison of genetic differentiation in neutral loci versus the Clock gene locus provided evidence for stabilizing selection. Furthermore, in all genetic markers, we found a greater genetic differentiation in space than in time. This pattern suggests that local adaptation or historic processes might have a stronger effect on the population structure and genetic variation in the pied flycatcher than recent global climate changes.     

Patterns of differentiation in a colour polymorphism and in neutral markers reveal rapid genetic changes in natural damselfly populations

The existence and mode of selection operating on heritable adaptive traits can be inferred by comparing population differentiation in neutral genetic variation between populations (often using F(ST) values) with the corresponding estimates for adaptive traits. Such comparisons indicate if selection acts in a diversifying way between populations, in which case differentiation in selected traits is expected to exceed differentiation in neutral markers [F(ST )(selected) > F(ST )(neutral)], or if negative frequency-dependent selection maintains genetic polymorphisms and pulls populations towards a common stable equilibrium [F(ST) (selected) < F(ST) (neutral)]. Here, we compared F(ST) values for putatively neutral data (obtained using amplified fragment length polymorphism) with estimates of differentiation in morph frequencies in the colour-polymorphic damselfly Ischnura elegans. We found that in the first year (2000), population differentiation in morph frequencies was significantly greater than differentiation in neutral loci, while in 2002 (only 2 years and 2 generations later), population differentiation in morph frequencies had decreased to a level significantly lower than differentiation in neutral loci. Genetic drift as an explanation for population differentiation in morph frequencies could thus be rejected in both years. These results indicate that the type and/or strength of selection on morph frequencies in this system can change substantially between years. We suggest that an approach to a common equilibrium morph frequency across all populations, driven by negative frequency-dependent selection, is the cause of these temporal changes. We conclude that inferences about selection obtained by comparing F(ST) values from neutral and adaptive genetic variation are most useful when spatial and temporal data are available from several populations and time points and when such information is combined with other ecological sources of data.     

A program to compare genetic differentiation statistics across loci using resampling of individuals and loci

Comparisons of genetic differentiation across populations based on different loci can provide insight into the evolutionary patterns acting on various regions of genomes. Here, we develop a program to statistically compare population genetic differentiation statistics (F(ST) or G'(ST) ) calculated from different loci. The program employs a routine that resamples either or both of individuals and loci and calculates a bootstrap confidence interval in the statistics. Resampling individuals is important when fewer than 25 individuals are sampled per population and when confidence intervals are required for individual loci. Resampling loci provides confidence intervals for sets of loci, such as a set presumed to be neutral, but can be anticonservative if fewer than 20 loci are analysed. We demonstrate the program using previously published data on the genetic differentiation at a major histocompatibility complex locus and at microsatellite loci across 10 populations of the guppy (Poecilia reticulata).     

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.     

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.     

Adaptive and neutral markers both show continent‐wide population structure of mountain pine beetle (Dendroctonus ponderosae)

Assessments of population genetic structure and demographic history have traditionally been based on neutral markers while explicitly excluding adaptive markers. In this study, we compared the utility of putatively adaptive and neutral single‐nucleotide polymorphisms (SNPs) for inferring mountain pine beetle population structure across its geographic range. Both adaptive and neutral SNPs, and their combination, allowed range‐wide structure to be distinguished and delimited a population that has recently undergone range expansion across northern British Columbia and Alberta. Using an equal number of both adaptive and neutral SNPs revealed that adaptive SNPs resulted in a stronger correlation between sampled populations and inferred clustering. Our results suggest that adaptive SNPs should not be excluded prior to analysis from neutral SNPs as a combination of both marker sets resulted in better resolution of genetic differentiation between populations than either marker set alone. These results demonstrate the utility of adaptive loci for resolving population genetic structure in a nonmodel organism.     

Assessment of population genetic structure in common wild rice Oryza rufipogon Griff. using microsatellite and allozyme markers

The genetic structure of five natural populations of common wild rice Oryza rufipogon Griff. from China, was investigated with 21 microsatellite loci and compared to estimates of genetic diversity and genetic differentiation detected by 22 allozyme loci. Microsatellite loci, as expected, have much higher levels of genetic diversity (mean values of A = 3.1, P = 73.3%, Ho = 0.358 and He = 0.345) than allozyme loci (mean values of A = 1.2, P = 12.7%, Ho = 0.020 and He = 0.030). Genetic differentiation detected by microsatellite loci ( FST = 0.468, mean I = 0.472) was higher than that for allozyme loci ( FST =0.388, mean I = 0.976). However, microsatellite markers showed less deviation from Hardy-Weinberg expectation (Wright's inbreeding coefficient FIS = -0.069) than do allozymes ( FIS = 0.337). These results suggest that microsatellite markers are powerful high-resolution tools for the accurate assessment of important parameters in population biology and conservation genetics of O. rufipogon, and offer advantages over allozyme markers.     

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.     

Heterogeneous genomic differentiation between walking-stick ecotypes: "isolation by adaptation" and multiple roles for divergent selection

Genetic differentiation can be highly variable across the genome. For example, loci under divergent selection and those tightly linked to them may exhibit elevated differentiation compared to neutral regions. These represent "outlier loci" whose differentiation exceeds neutral expectations. Adaptive divergence can also increase genome-wide differentiation by promoting general barriers to neutral gene flow, thereby facilitating genomic divergence via genetic drift. This latter process can yield a positive correlation between adaptive phenotypic divergence and neutral genetic differentiation (described here as "isolation-by-adaptation"). Here, we examine both these processes by combining an AFLP genome scan of two host plant ecotypes of Timema cristinae walking-sticks with existing data on adaptive phenotypic divergence and ecological speciation in these insects. We found that about 8% of loci are outliers in multiple population comparisons. Replicated comparisons between population-pairs using the same versus different host species revealed that 1-2% of loci are subject to host-related selection specifically. Locus-specific analyses revealed that up to 10% of putatively neutral (nonoutlier) AFLP loci exhibit significant isolation-by-adaptation. Our results suggest that selection may affect differentiation directly, via linkage, or by facilitating genetic drift. They thus illustrate the varied and sometimes nonintuitive contributions of selection to heterogeneous genomic differentiation.     

A standardized genetic differentiation measure

Interpretation of genetic differentiation values is often problematic because of their dependence on the level of genetic variation. For example, the maximum level of GST is less than the average within population homozygosity so that for highly variable loci, even when no alleles are shared between subpopulations, GST may be low. To remedy this difficulty, a standardized measure of genetic differentiation is introduced here, one which has the same range, 0-1, for all levels of genetic variation. With this measure, the magnitude is the proportion of the maximum differentiation possible for the level of subpopulation homozygosity observed. This is particularly important for situations in which the mutation rate is of the same magnitude or higher than the rate of gene flow. The standardized measure allows comparison between loci with different levels of genetic variation, such as allozymes and microsatellite loci, or mtDNA and Y-chromosome genes, and for genetic differentiation for organisms with different effective population sizes.     

The effect on neutral gene flow of selection at a linked locus

The effect on gene flow at a neutral locus of a selective cline at a linked locus is investigated. A diffusion approximation for a two-locus island model is derived in which only one locus is subject to selection. The moments of the stationary distribution are obtained and compared to the corresponding moments from a one-locus, neutral island model. This comparison yields an effective migration rate. The effective migration rate is always less than the actual migration rate, but this effect is seen to be small for weak selection and loose linkage in the case of adult migration. The importance of selection at linked loci to the question of genetic differentiation in a subdivided population is discussed.     

Genetic Diversity Within and Among Populations of Viola tenuicornis with Reference to Sampling Strategies

Genetic diversity and population differentiation of Viola tenuicornis W. Beck. were studied in nine populations in Beijing using horizontal starch gel electrophoresis. Allozyme data for 11 loci of 8 enzymes indicated high level of genetic variability in this region with P=72.7％, A=2.4 and H e=0.243. In contrast, the population differentiation was lower than the average of other perennial herb with F st=0.196. It was found that distances among populations had little effect on population differentiation. A sample of 30 individuals could well represent a population, but less than 20 individuals per population might lead to significant bias in the estimation of genetic parameters, in particular the values of A and H e.     

Explorative genome scan to detect candidate loci for adaptation along a gradient of altitude in the common frog (Rana temporaria)

Today, with the rapid development of population genomics, the genetic basis of adaptation can be unraveled directly at the genome level, without any prerequisites about the selectively advantageous genes or traits. For nonmodel species, it is now possible to screen many markers randomly scattered across the genome and to distinguish between the neutral genetic background and outlier loci displaying an atypical behavior (e.g., a higher differentiation between populations). This study investigated the genetic frame of adaptation to a gradient of altitude in the common frog (Rana temporaria) by means of a genome scan based on 392 amplified fragment length polymorphism markers. Using two outlier detection methods never applied to dominant data so far, we sought for loci with a genetic differentiation diverging from neutral expectations when comparing populations from different altitudes. All the detected loci were sorted out according to their most probable cause for outlier behavior and classified as false positives, outliers due to local effects, or outliers associated with altitude. Altogether, eight good candidate loci were identified as potentially involved in adaptation to altitude because they were picked out in several independent interaltitude comparisons. This result illustrated the potential of genome-wide surveys to reveal selection signatures along selection gradients, where the association between environmental variables and fitness-related traits may be complex and/or cryptic. In this article, we also underlined the need for confirmation of the selection footprints for the outlier loci. Finally, we provided some preliminary insights into the genetic basis of adaptation along an altitudinal cline in the common frog.     

A comparative analysis of population structure of a forest tree, Eucalyptus globulus (Myrtaceae), using microsatellite markers and quantitative traits

Eucalyptus globulus (Myrtaceae) is a forest tree native to southeastern Australia, but is grown globally for pulpwood and timber. Eight microsatellite loci were used to determine the degree of selectively neutral differentiation between native populations of the geographic races of E. globulus that are used in a national breeding programme. Spatial differentiation was detected among 340 samples from across the species range (F _ST=0.09±0.02). Analysis of molecular variance showed that there was significant variation between the races, and an unweighted pair group method with arithmetic mean analysis of Nei’s genetic distance between races showed that geographically proximal races tended to be more closely related than geographically distant races. This contrasted markedly with analyses based on quantitative genetic data, where some races appeared to be highly divergent from their geographically closest neighbours. Comparison of racial differentiation based on quantitative (Q _ST) and molecular (F _ST) data suggested that at least five of the quantitative traits used for defining races of E. globulus have been influenced by natural selection, resulting in cases of both phenotypic divergence of parapatric races and phenotypic convergence of allopatric races. We conclude that selectively neutral molecular markers are more useful than quantitative genetic data for identifying the evolutionary affinities and lineages within E. globulus. However, both sources of information should be used in defining evolutionarily important units for conservation. The population structure observed in E. globulus has important consequences for future association studies and may also affect breeding strategies if significant genome co-adaptation has occurred.     

Divergent selection and heterogeneous genomic divergence

Levels of genetic differentiation between populations can be highly variable across the genome, with divergent selection contributing to such heterogeneous genomic divergence. For example, loci under divergent selection and those tightly physically linked to them may exhibit stronger differentiation than neutral regions with weak or no linkage to such loci. Divergent selection can also increase genome‐wide neutral differentiation by reducing gene flow (e.g. by causing ecological speciation), thus promoting divergence via the stochastic effects of genetic drift. These consequences of divergent selection are being reported in recently accumulating studies that identify: (i) ‘outlier loci’ with higher levels of divergence than expected under neutrality, and (ii) a positive association between the degree of adaptive phenotypic divergence and levels of molecular genetic differentiation across population pairs [‘isolation by adaptation’ (IBA)]. The latter pattern arises because as adaptive divergence increases, gene flow is reduced (thereby promoting drift) and genetic hitchhiking increased. Here, we review and integrate these previously disconnected concepts and literatures. We find that studies generally report 5–10% of loci to be outliers. These selected regions were often dispersed across the genome, commonly exhibited replicated divergence across different population pairs, and could sometimes be associated with specific ecological variables. IBA was not infrequently observed, even at neutral loci putatively unlinked to those under divergent selection. Overall, we conclude that divergent selection makes diverse contributions to heterogeneous genomic divergence. Nonetheless, the number, size, and distribution of genomic regions affected by selection varied substantially among studies, leading us to discuss the potential role of divergent selection in the growth of regions of differentiation (i.e. genomic islands of divergence), a topic in need of future investigation.     

The Genetic Structure of Natural Populations of DROSOPHILA MELANOGASTER. Xviii. Clinal and Uniform Genetic Variation over Populations

It has been reported in the previous papers of this series that in the eastern United States and Japan there is a north-to-south cline of additive genetic variance of viability and that the amount of the additive genetic variance in the northern population can be explained by mutation-selection balance. To determine whether or not the difference in the genetic variation in northern and southern populations can be explained by the differences in mutation rate and/or effective population size, numerical calculations were made using population genetic parameters. In addition, the average heterozygosities of the northern and southern populations at ten of 19 polymorphic structural loci surveyed were estimated in relation to the cline of additive genetic variance of viability, and the following findings were obtained. (1) The changes in mutation rate and population size cannot simultaneously explain the difference in additive genetic variance and inbreeding decline between the northern and southern populations. Thus, the operation of some kind of balancing selection, most likely diversifying selection, was suggested to explain the observed excess of additive genetic variance. (2) Estimates of the average heterozygosities of the southern population were not significantly different from those of the northern population. Thus, it was strongly suggested that the excess of additive genetic variance in the southern population cannot be caused by structural loci, but by factors outside the structural loci, and that protein polymorphisms are selectively neutral or nearly neutral.     

Is MHC diversity a better marker for conservation than neutral genetic diversity? A case study of two contrasting dolphin populations

Genetic diversity is essential for populations to adapt to changing environments. Measures of genetic diversity are often based on selectively neutral markers, such as microsatellites. Genetic diversity to guide conservation management, however, is better reflected by adaptive markers, including genes of the major histocompatibility complex (MHC). Our aim was to assess MHC and neutral genetic diversity in two contrasting bottlenose dolphin (Tursiops aduncus) populations in Western Australia—one apparently viable population with high reproductive output (Shark Bay) and one with lower reproductive output that was forecast to decline (Bunbury). We assessed genetic variation in the two populations by sequencing the MHC class II DQB, which encompasses the functionally important peptide binding regions (PBR). Neutral genetic diversity was assessed by genotyping twenty‐three microsatellite loci. We confirmed that MHC is an adaptive marker in both populations. Overall, the Shark Bay population exhibited greater MHC diversity than the Bunbury population—for example, it displayed greater MHC nucleotide diversity. In contrast, the difference in microsatellite diversity between the two populations was comparatively low. Our findings are consistent with the hypothesis that viable populations typically display greater genetic diversity than less viable populations. The results also suggest that MHC variation is more closely associated with population viability than neutral genetic variation. Although the inferences from our findings are limited, because we only compared two populations, our results add to a growing number of studies that highlight the usefulness of MHC as a potentially suitable genetic marker for animal conservation. The Shark Bay population, which carries greater adaptive genetic diversity than the Bunbury population, is thus likely more robust to natural or human‐induced changes to the coastal ecosystem it inhabits.