Hitchhiking effect of a beneficial mutation spreading in a subdivided population

A central problem in population genetics is to detect and analyze positive natural selection by which beneficial mutations are driven to fixation. The hitchhiking effect of a rapidly spreading beneficial mutation, which results in local removal of standing genetic variation, allows such an analysis using DNA sequence polymorphism. However, the current mathematical theory that predicts the pattern of genetic hitchhiking relies on the assumption that a beneficial mutation increases to a high frequency in a single random-mating population, which is certainly violated in reality. Individuals in natural populations are distributed over a geographic space. The spread of a beneficial allele can be delayed by limited migration of individuals over the space and its hitchhiking effect can also be affected. To study this effect of geographic structure on genetic hitchhiking, we analyze a simple model of directional selection in a subdivided population. In contrast to previous studies on hitchhiking in subdivided populations, we mainly investigate the range of sufficiently high migration rates that would homogenize genetic variation at neutral loci. We provide a heuristic mathematical analysis that describes how the genealogical structure at a neutral locus linked to the locus under selection is expected to change in a population divided into two demes. Our results indicate that the overall strength of genetic hitchhiking--the degree to which expected heterozygosity decreases--is diminished by population subdivision, mainly because opportunity for the breakdown of hitchhiking by recombination increases as the spread of the beneficial mutation across demes is delayed when migration rate is much smaller than the strength of selection. Furthermore, the amount of genetic variation after a selective sweep is expected to be unequal over demes: a greater reduction in expected heterozygosity occurs in the subpopulation from which the beneficial mutation originates than in its neighboring subpopulations. This raises a possibility of detecting a "hidden" geographic structure of population by carefully analyzing the pattern of a selective sweep.     

The effect of hitch-hiking on genes linked to a balanced polymorphism in a subdivided population

The effect of multi-allelic balancing selection on nucleotide diversity at linked neutral sites was investigated by simulations of subdivided populations. The motivation is to understand the behaviour of self-recognition systems such as the MHC and plant self-incompatibility. For neutral sites, two types of subdivision are present: (1) into demes (connected by migration), and (2) into classes defined by different functional alleles at the selected locus (connected by recombination). Previous theoretical studies of each type of subdivision separately have shown that each increases diversity, and decreases the relative frequencies of low-frequency variants, at neutral sites or loci. We show here that the two types of subdivision act non-additively when sampling is at the whole population level, and that subdivision produces some non-intuitive results. For instance, in highly subdivided populations, genetic diversity at neutral sites may decrease with tighter linkage to a selected locus or site. Another conclusion is that, if there is population subdivision, balancing selection leads to decreased expected FST values for neutral sites linked to the selected locus. Finally, we show that the ability to detect balancing selection by its effects on linked variation, using tests such as Tajima's D, is reduced when genes in a subdivided population are sampled from the total population, rather than within demes.     

Effects of environmental heterogeneity on victim-exploiter coevolution

We study victim-exploiter coevolution in a spatially heterogeneous island model. In each species, fitness consequences of between-species interactions are controlled by a single haploid diallelic locus. Our emphasis is on the conditions for the maintenance of genetic variation, the dynamic patterns observed, the extent of local adaptation and genetic differentiation between different demes, and on how different parameters (such as the strength and heterogeneity in selection, migration rates, and the number of sites) affect the dynamic and static behavior of the system. We show that under spatially homogeneous selection the maintenance of genetic variation is possible through asynchronous nonlinear dynamics where the allele frequencies in a majority of demes quickly synchronize but the rest do not. Spatially heterogeneous selection can maintain genetic variation even if migration rates are maximal. This happens in an oscillatory way. Genetic variation is most likely to be maintained at high levels if the heterogeneity in selection is large. If there are some restrictions on migration, genetic variation can be maintained at a stable equilibrium. This behavior is most likely at intermediate migration rates. In this case, the system can exhibit high spatial subdivision as measured by F(ST) values but relatively low local adaptation.     

When can ecological speciation be detected with neutral loci?

It is not yet clear under what conditions empirical studies can reliably detect progress toward ecological speciation through the analysis of allelic variation at neutral loci. We use a simulation approach to investigate the range of parameter space under which such detection is, and is not, likely. We specifically test for the conditions under which divergent natural selection can cause a ‘generalized barrier to gene flow’ that is present across the genome. Our individual‐based numerical simulations focus on how population divergence at neutral loci varies in relation to recombination rate with a selected locus, divergent selection on that locus, migration rate and population size. We specifically test whether genetic differences at neutral markers are greater between populations in different environments than between populations in similar environments. We find that this expected signature of ecological speciation can be detected under part of the parameter space, most consistently when divergent selection is strong and migration is intermediate. By contrast, the expected signature of ecological speciation is not reliably detected when divergent selection is weak or migration is low or high. These findings provide insights into the strengths and weaknesses of using neutral markers to infer ecological speciation in natural systems.     

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.     

Chaotic genetic patchiness in the pelagic teleost fish Sardina pilchardus across the Siculo-Tunisian Strait

ABSTRACT

This investigation aims at assessing patterns of spatial genetic structure of the teleost fish Sardina pilchardus across the Siculo-Tunisian Strait (a well-known discontinuous biogeographic area) and delineating putative genetic stocks within the species. For this purpose, a total of 180 specimens, collected from 11 locations stretching across the western and eastern Mediterranean coasts of Tunisia, were analysed genetically by means of 18 nuclear allozyme loci. The outcome of this study revealed strong genetic differentiation among populations, with the marked genetic distinctiveness of the central Tunisian population at Mahdia. Despite the delineation of seven well-defined genetic groups, no significant correlation was found between genetic and geographic distances. Besides, the recorded population subdivision did not align with biogeographic boundaries, suggesting the presence of chaotic genetic patchiness. Recent genetic bottlenecks were evidenced in S. pilchardus populations. Patchy migration patterns were recorded among the examined pairs of sardine populations. Among the recorded 16 polymorphic loci, GPI-2 and SOD appeared to be subject to natural selection. Patterns of population genetic differentiation and structuring were not found to be driven by outlier loci that appeared to be under selection. Furthermore, the detected neutral GPI-1 locus was found to be responsible for most of the genetic variation among identified genetic clusters. Hence, natural selection cannot cause the detected genetic heterogeneity among sardine samples. Different explanations to the origin of chaotic genetic patterns, observed within S. pilchardus, were discussed.     

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.     

The maintenance of genetic diversity under host–parasite coevolution in finite,structured populations

As a corollary to the Red Queen hypothesis, host–parasite coevolution has been hypothesized to maintain genetic variation in both species. Recent theoretical work, however, suggests that reciprocal natural selection alone is insufficient to maintain variation at individual loci. As highlighted by our brief review of the theoretical literature, models of host–parasite coevolution often vary along multiple axes (e.g. inclusion of ecological feedbacks or abiotic selection mosaics), complicating a comprehensive understanding of the effects of interacting evolutionary processes on diversity. Here we develop a series of comparable models to explore the effect of interactions between spatial structures and antagonistic coevolution on genetic diversity. Using a matching alleles model in finite populations connected by migration, we find that, in contrast to panmictic populations, coevolution in a spatially structured environment can maintain genetic variation relative to neutral expectations with migration alone. These results demonstrate that geographic structure is essential for understanding the effect of coevolution on biological diversity.     

Perfect simulation from nonneutral population genetic models: variable population size and population subdivision

We show how the idea of monotone coupling from the past can produce simple algorithms for simulating samples at a nonneutral locus under a range of demographic models. We specifically consider a biallelic locus and either a general variable population size mode or a general migration model for population subdivision. We investigate the effect of demography on the efficacy of selection and the effect of selection on genetic divergence between populations.     

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.     

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.     

Population amalgamation and genetic variation: observations on artificially agglomerated tribal populations of Central and South America.

The interpretation of data on genetic variation with regard to the relative roles of different evolutionary factors that produce and maintain genetic variation depends critically on our assumptions concerning effective population size and the level of migration between neighboring populations. In humans, recent population growth and movements of specific ethnic groups across wide geographic areas mean that any theory based on assumptions of constant population size and absence of substructure is generally untenable. We examine the effects of population subdivision on the pattern of protein genetic variation in a total sample drawn from an artificial agglomerate of 12 tribal populations of Central and South America, analyzing the pooled sample as though it were a single population. Several striking findings emerge. (1) Mean heterozygosity is not sensitive to agglomeration, but the number of different alleles (allele count) is inflated, relative to neutral mutation/drift/equilibrium expectation. (2) The inflation is most serious for rare alleles, especially those which originally occurred as tribally restricted "private" polymorphisms. (3) The degree of inflation is an increasing function of both the number of populations encompassed by the sample and of the genetic divergence among them. (4) Treating an agglomerated population as though it were a panmictic unit of long standing can lead to serious biases in estimates of mutation rates, selection pressures, and effective population sizes. Current DNA studies indicate the presence of numerous genetic variants in human populations. The findings and conclusions of this paper are all fully applicable to the study of genetic variation at the DNA level as well.     

Alu insertion polymorphisms for the study of human genomic diversity

We analyze genetic variation at fused1, a locus that is close to the centromere of the X chromosome-autosome (X/4) fusion in Drosophila americana. In contrast to other X-linked and autosomal genes, for which a lack of population subdivision in D. americana has been observed at the DNA level, we find strong haplotype structure associated with the alternative chromosomal arrangements. There are several derived fixed differences at fused1 (including one amino acid replacement) between two haplotype classes of this locus. From these results, we obtain an estimate of an age of approximately 0.61 million years for the origin of the two haplotypes of the fused1 gene. Haplotypes associated with the X/4 fusion have less DNA sequence variation at fused1 than haplotypes associated with the ancestral chromosome arrangement. The X/4 haplotypes also exhibit clinal variation for the allele frequencies of the three most common amino acid replacement polymorphisms, but not for adjacent silent polymorphisms. These patterns of variation are best explained as a result of selection acting on amino acid substitutions, with geographic variation in selection pressures.     

Evolutionary Genomics of a Subdivided Species

The ways in which genetic variation is distributed within and among populations is a key determinant of the evolutionary features of a species. However, most comprehensive studies of these features have been restricted to studies of subdivision in settings known to have been driven by local adaptation, leaving our understanding of the natural dispersion of allelic variation less than ideal. Here, we present a geographic population-genomic analysis of 10 populations of the freshwater microcrustacean Daphnia pulex, an emerging model system in evolutionary genomics. These populations exhibit a pattern of moderate isolation-by-distance, with an average migration rate of 0.6 individuals per generation, and average effective population sizes of ∼650,000 individuals. Most populations contain numerous private alleles, and genomic scans highlight the presence of islands of excessively high population subdivision for more common alleles. A large fraction of such islands of population divergence likely reflect historical neutral changes, including rare stochastic migration and hybridization events. The data do point to local adaptive divergence, although the precise nature of the relevant variation is diffuse and cannot be associated with particular loci, despite the very large sample sizes involved in this study. In contrast, an analysis of between-species divergence highlights positive selection operating on a large set of genes with functions nearly nonoverlapping with those involved in local adaptation, in particular ribosome structure, mitochondrial bioenergetics, light reception and response, detoxification, and gene regulation. These results set the stage for using D. pulex as a model for understanding the relationship between molecular and cellular evolution in the context of natural environments.     

Variation after a selective sweep in a subdivided population

The effect of genetic hitchhiking on neutral variation is analyzed in subdivided populations with differentiated demes. After fixation of a favorable mutation, the consequences on particular subpopulations can be radically different. In the subpopulation where the mutation first appeared by mutation, variation at linked neutral loci is expected to be reduced, as predicted by the classical theory. However, the effect in the other subpopulations, where the mutation is introduced by migration, can be the opposite. This effect depends on the level of genetic differentiation of the subpopulations, the selective advantage of the mutation, the recombination frequency, and the population size, as stated by analytical derivations and computer simulations. The characteristic outcomes of the effect are three. First, the genomic region of reduced variation around the selected locus is smaller than that predicted in a panmictic population. Second, for more distant neutral loci, the amount of variation increases over the level they had before the hitchhiking event. Third, for these loci, the spectrum of gene frequencies is dominated by an excess of alleles at intermediate frequencies when compared with the neutral theory. At these loci, hitchhiking works like a system that takes variation from the between-subpopulation component and introduces it into the subpopulations. The mechanism can also operate in other systems in which the genetic variation is distributed in clusters with limited exchange of variation, such as chromosome arrangements or genomic regions closely linked to targets of balancing selection.     

Variation on islands: major histocompatibility complex (Mhc) polymorphism in populations of the Australian bush rat

Loss of genetic variation in small, isolated populations is commonly observed at neutral or nearly neutral loci. In this study, the loss of genetic variation was assessed in island populations for a locus of major histocompatibility complex (Mhc), a locus shown to be under the influence of balancing selection. A total of 36 alleles was found at the second exon of RT1.Ba in 14 island and two mainland populations of Rattus fuscipes greyii. Despite this high overall diversity, a substantial lack of variation was observed in the small island populations, with 13 islands supporting only one to two alleles. Two populations, Waldegrave and Williams Islands, showed moderately high levels of heterozygosity (52-56%) which were greater than expected under neutrality, suggesting the action of balancing selection. However, congruence between the level of variation at this Mhc locus and in previous allozyme electrophoresis and mitochondrial DNA studies highlights the dominant influence of genetic drift and population factors, such as bottlenecks and structuring in the founding population, in the loss of genetic variation in these small, isolated populations.     

The effects of the mating system on the evolution of migration in a spatially heterogeneous population

Summary Verbal explanations for the evolution of migration and dispersal often invoke inbreeding depression as an important force. Experimental work on plant populations indicates that while inbreeding depression may favor increased migration rates, adaptation to local environments may reduce the advantage to migrants. We formalize and test this hypothesis using a two-locus genetic model that incorporates lowered fitness in offspring produced by self-fertilization, and habitat differentiation. We also use the model to address questions about the general theory of genetic modifiers and the modifier reduction principle. We find that even under conditions when migration would increase the mean fitness of a population, migration may not be favored. This result is due to the associations that develop between genotypes at a locus subject to overdominant selection and at a neutral locus controlling the migration rate. Thus, it appears that, in this model, the forces of local adaptation, which favor a reduction in the migration rate, overwhelm those of inbreeding depression, which may favor dispersal.     

The Effect of Linkage on Establishment and Survival of Locally Beneficial Mutations

We study invasion and survival of weakly beneficial mutations arising in linkage to an established migration–selection polymorphism. Our focus is on a continent–island model of migration, with selection at two biallelic loci for adaptation to the island environment. Combining branching and diffusion processes, we provide the theoretical basis for understanding the evolution of islands of divergence, the genetic architecture of locally adaptive traits, and the importance of so-called “divergence hitchhiking” relative to other mechanisms, such as “genomic hitchhiking”, chromosomal inversions, or translocations. We derive approximations to the invasion probability and the extinction time of a de novo mutation. Interestingly, the invasion probability is maximized at a nonzero recombination rate if the focal mutation is sufficiently beneficial. If a proportion of migrants carries a beneficial background allele, the mutation is less likely to become established. Linked selection may increase the survival time by several orders of magnitude. By altering the timescale of stochastic loss, it can therefore affect the dynamics at the focal site to an extent that is of evolutionary importance, especially in small populations. We derive an effective migration rate experienced by the weakly beneficial mutation, which accounts for the reduction in gene flow imposed by linked selection. Using the concept of the effective migration rate, we also quantify the long-term effects on neutral variation embedded in a genome with arbitrarily many sites under selection. Patterns of neutral diversity change qualitatively and quantitatively as the position of the neutral locus is moved along the chromosome. This will be useful for population-genomic inference. Our results strengthen the emerging view that physically linked selection is biologically relevant if linkage is tight or if selection at the background locus is strong.     

Spatio-temporal variation in the strength and mode of selection acting on major histocompatibility complex diversity in water vole (Arvicola terrestris) metapopulations

Patterns of spatio-temporal genetic variation at a class II major histocompatibility complex (MHC) locus and multiple microsatellite loci were analysed within and between three water vole metapopulations in Scotland, UK. Comparisons of MHC and microsatellite spatial genetic differentiation, based on standardised tests between two demographically asynchronous zones within a metapopulation, suggested that spatial MHC variation was affected by balancing selection, directional selection and random genetic drift, but that the relative effects of these microevolutionary forces vary temporally. At the metapopulation level, between-year differentiation for MHC loci was significantly correlated with that of microsatellites, signifying that neutral factors such as migration and drift were primarily responsible for overall temporal genetic change at the metapopulation scale. Between metapopulations, patterns of genetic differentiation implied that, at large spatial scales, MHC variation was primarily affected by directional selection and drift. Levels of MHC heterozygosity in excess of Hardy–Weinberg expectations were consistent with overdominant balancing selection operating on MHC variation within metapopulations. However, this effect was not constant among all samples, indicating temporal variation in the strength of selection relative to other factors. The results highlight the benefit of contrasting variation at MHC with neutral markers to separate the effects of stochastic and deterministic microevolutionary forces, and add to a growing body of evidence showing that the mode and relative strength of selection acting on MHC diversity varies both spatially and temporally.     

Genetic profile of cosmopolitan populations: effects of hidden subdivision

Natural populations of many organisms exhibit excess of rare alleles in comparison with the predictions of the neutral mutation hypothesis. It has been shown before that either a population bottleneck or the presence of slightly deleterious mutations can explain this phenomenon. A third explanation is presented in this work, showing that hidden subdivision within a population can also lead to an excess of rare alleles in the total population when the expectations of the neutral model are based on the allele frequency profile of the entire population data. With two examples (mitochondrial DNA-morph distribution and isozyme allele frequency distributions), it is shown that most cosmopolitan human populations exhibit excess of rare as well as total allele counts, when these are compared with the expectations of the neutral mutation hypothesis. The mitochondrial data demonstrate that such excesses can be detected from genetic variation at a single locus as well, and this is not due to stochastic error of allele frequency distributions. Contrast of the present observations with the allele frequency profiles in agglomerated tribal populations from South and Central America shows that even when the neutral expectations hold for individual subpopulations, if all subpopulations are grouped into a single population, the pooled data exhibit an excess of total number of alleles that is mainly due to the excess of rare alleles. Therefore, a primary cause of the excess number of rare alleles could be the hidden subdivision, and the magnitude of the excess indicates the extent of substructuring. The two components of hidden subdivision are: 1) Number of subpopulations, and 2) the average genetic distance among them. The implications of this observation in estimating mutation rate are discussed indicating the difficulties of comparing mutation rates from different population surveys.