Effects of metapopulation processes on measures of genetic diversity

Many species persist as a metapopulation under a balance between the local extinction of subpopulations or demes and their recolonization through dispersal from occupied patches. Here we review the growing body of literature dealing with the genetic consequences of such population turnover. We focus our attention principally on theoretical studies of a classical metapopulation with a 'finite-island' model of population structure, rather than on 'continent-island' models or 'source-sink' models. In particular, we concern ourselves with the subset of geographically subdivided population models in which it is assumed that all demes are liable to extinction from time to time and that all demes receive immigrants. Early studies of the genetic effects of population turnover focused on population differentiation, such as measured by F(ST). A key advantage of F(ST) over absolute measures of diversity is its relative independence of the mutation process, so that different genes in the same species may be compared. Another advantage is that F(ST) will usually equilibrate more quickly following perturbations than will absolute levels of diversity. However, because F(ST) is a ratio of between-population differentiation to total diversity, the genetic effects of metapopulation processes may be difficult to interpret in terms of F(ST) on its own, so that the analysis of absolute measures of diversity in addition is likely to be informative. While population turnover may either increase or decrease F(ST), depending on the mode of colonization, recurrent extinction and recolonization is expected always to reduce levels of both within-population and species-wide diversity (piS and piT, respectively). One corollary of this is that piS cannot be used as an unbiased estimate of the scaled mutation rate, theta, as it can, with some assumptions about the migration process, in species whose demes do not fluctuate in size. The reduction of piT in response to population turnover reflects shortened mean coalescent times, although the distribution of coalescence times under extinction colonization equilibrium is not yet known. Finally, we review current understanding of the effect of metapopulation dynamics on the effective population size.     

Genetic load,inbreeding depression and heterosis in an age‐structured metapopulation

In a metapopulation, the process of recurrent local extinction and recolonization gives rise to an age structure among demes. Recently established demes will tend to differ from older demes in terms of the levels of genetic diversity found within them and the way this diversity is distributed among demes in the same and different ages. The effects of population turnover on average levels of genetic diversity among demes in a metapopulation have been the focus of much attention, both for neutral and nonneutral loci, but much less is known about the distribution of nonneutral genetic diversity among demes of different ages. In this paper, we used computer simulations to study the distribution of genetic load, inbreeding depression and heterosis in an age‐structured metapopulation. We found that, for mildly deleterious mutations, within‐deme inbreeding depression increased, whereas heterosis and genetic load decreased with deme age following severe colonization bottlenecks. In contrast, recessive lethal alleles tended to be purged during colonization, with older populations showing higher genetic load and higher within‐deme inbreeding depression. Heterosis caused by recessive lethal alleles and resulting from gene flow among different demes tended to be greatest for young demes, because the mutations responsible tended to be purged in the first few generations after colonization, but its effects increased again as populations grow older as a result of immigration. Our results point to a need for estimates of genetic diversity, genetic load, within‐deme inbreeding depression and heterosis in demes of different age classes separately.     

Coalescence in a metapopulation with recurrent local extinction and recolonization

Abstract Many species exist as metapopulations in balance between local population extinction and recolonization. The effect of these processes on average population differentiation, within-deme diversity, and specieswide diversity has been considered previously. In this paper, coalescent simulations of Slatkin's propagule-pool and migrant-pool models are used to characterize the distribution of neutral genetic diversity within demes (π_s), diversity in the metapopulation a whole (TT_T), the ratio F _ST= (π_t–π_S)/π_T, Tajima's D statistic, and several ratios of gene-tree branch lengths. Using these distributions, power to detect differences in key metapopulation parameter values is determined under contrasting sampling regimes. The results indicate that it will be difficult to use sequence data from a single locus to detect a history of extinctions and recolonizations in a metapopulation because of high genealogical variance, the loss of diversity due to reductions in effective population size, and the fact that a genealogy of lineages from different demes under Slatkin's model differs from a neutral coalescent only in its time scale. Genetic indices of gene-tree shape that capture the effects of extinction/recolonization on both external branches and the length of the genealogy as a whole will provide the best indication of metapopulation dynamics if several lineages are sampled from several different demes.     

Linking extinction-colonization dynamics to genetic structure in a salamander metapopulation

Theory predicts that founder effects have a primary role in determining metapopulation genetic structure. However, ecological factors that affect extinction-colonization dynamics may also create spatial variation in the strength of genetic drift and migration. We tested the hypothesis that ecological factors underlying extinction-colonization dynamics influenced the genetic structure of a tiger salamander (Ambystoma tigrinum) metapopulation. We used empirical data on metapopulation dynamics to make a priori predictions about the effects of population age and ecological factors on genetic diversity and divergence among 41 populations. Metapopulation dynamics of A. tigrinum depended on wetland area, connectivity and presence of predatory fish. We found that newly colonized populations were more genetically differentiated than established populations, suggesting that founder effects influenced genetic structure. However, ecological drivers of metapopulation dynamics were more important than age in predicting genetic structure. Consistent with demographic predictions from metapopulation theory, genetic diversity and divergence depended on wetland area and connectivity. Divergence was greatest in small, isolated wetlands where genetic diversity was low. Our results show that ecological factors underlying metapopulation dynamics can be key determinants of spatial genetic structure, and that habitat area and isolation may mediate the contributions of drift and migration to divergence and evolution in local populations.     

Founder events as determinants of within-island and among-island genetic structure of Daphnia metapopulations

The genetic structure of metapopulations offers insights into the genetic consequences of local extinction and recolonization. We studied allozyme variation in rock pool metapopulations of two species of waterfleas (Daphnia) with the aim to understand how these dynamics influence genetic differentiation. We screened 138 populations of D. magna and 65 populations of D. longispina from an area in the archipelago of southern Finland. The pools from which they were sampled are separated by distances between 1.5 and 4710 m and located on a total of 38 islands. The genetic population structure of the two species was strikingly similar, consistent with their similar metapopulation ecology. The mean F(PT) value (differentiation among pools with respect to the total metapopulation) was 0.55 and a hierarchical analysis showed that genetic differentiation was strong (>0.25) among pools within islands as well as among whole islands. Within islands, pairwise genetic differentiation increased with geographic distance, indicating isolation by distance due to spatially limited dispersal. Previous studies have shown strong founder events occurring during colonization in our metapopulation. We suggest that the genetic population structure in the studied metapopulations is largely explained by three consequences of these founder events: (i) strong drift during colonization, (ii) local inbreeding, which results in hybrid vigour and increased effective migration rates after subsequent immigration, and (iii) effects of selection through hitchhiking of neutral genes with linked loci under selection.     

Migration and recovery of the genetic diversity during the increasing density phase in cyclic vole populations

In cyclic populations, high genetic diversity is currently reported despite the periodic low numbers experienced by the populations during the low phases. Here, we report spatio-temporal monitoring at a very fine scale of cyclic populations of the fossorial water vole (Arvicola terrestris) during the increasing density phase. This phase marks the transition from a patchy structure (demes) during low density to a continuous population in high density. We found that the genetic diversity was effectively high but also that it displayed a local increase within demes over the increasing phase. The genetic diversity remained relatively constant when considering all demes together. The increase in vole abundance was also correlated with a decrease of genetic differentiation among demes. Such results suggest that at the end of the low phase, demes are affected by genetic drift as the result of being small and geographically isolated. This leads to a loss of local genetic diversity and a spatial differentiation among demes. This situation is counterbalanced during the increasing phase by the spatial expansion of demes and the increase of the effective migration among differentiated demes. We provide evidences that in cyclic populations of the fossorial water voles, the relative influence of drift operating during low density populations and migration occurring principally while population size increases interacts closely to maintain high genetic diversity.     

NEUTRAL GENETIC DIVERSITY IN A METAPOPULATION WITH RECURRENT LOCAL EXTINCTION AND RECOLONIZATION

Many species exist as metapopulations in balance between local population extinction and recolonization, processes that may strongly affect the distribution of neutral genetic diversity within demes and in the metapopulation as a whole. In this paper we use both the infinite-alleles and the infinite-sites models to reframe Slatkin's propagulepool and migrant-pool models in terms of mean within-deme and among-deme genetic diversity; the infinite-sites model is particularly relevant to DNA sequence data. Population turnover causes a major reduction in neutral genetic diversity within demes, π_S, and in the metapopulation as a whole, π_t. This effect is particularly strong for propagulepool colonization, in which colonists are drawn from a single extant deme. Because metapopulation dynamics affect both within-deme and total metapopulation diversity similarly, comparisons between species with different ecologies on the basis of ratios such as F_ST are difficult to interpret and absolute measures of divergence between populations should be used as well. Although the value of F_ST in a metapopulation with local extinction depends strongly on the mode of colonization, this has almost no effect on the numerator of the F_ST ratio, π_t – π_S, so that F_ST is influenced mainly by the effect of the colonization mode on the denominator (π_t). Our results also indicate that it is inappropriate to use measures of average within-deme diversity in species with population turnover to estimate the scaled mutation rate, θ, because extinction can greatly reduce π_S. Finally, we discuss the effect of population turnover on the effective size of a metapopulation.     

Determinants of Genetic Structure in a Nonequilibrium Metapopulation of the Plant Silene latifolia

Population genetic differentiation will be influenced by the demographic history of populations, opportunities for migration among neighboring demes and founder effects associated with repeated extinction and recolonization. In natural populations, these factors are expected to interact with each other and their magnitudes will vary depending on the spatial distribution and age structure of local demes. Although each of these effects has been individually identified as important in structuring genetic variance, their relative magnitude is seldom estimated in nature. We conducted a population genetic analysis in a metapopulation of the angiosperm, Silene latifolia, from which we had more than 20 years of data on the spatial distribution, demographic history, and extinction and colonization of demes. We used hierarchical Bayesian methods to disentangle which features of the populations contributed to among population variation in allele frequencies, including the magnitude and direction of their effects. We show that population age, long-term size and degree of connectivity all combine to affect the distribution of genetic variance; small, recently-founded, isolated populations contributed most to increase F _ST in the metapopulation. However, the effects of population size and population age are best understood as being modulated through the effects of connectivity to other extant populations, i.e. F _ST diminishes as populations age, but at a rate that depends how isolated the population is. These spatial and temporal correlates of population structure give insight into how migration, founder effect and within-deme genetic drift have combined to enhance and restrict genetic divergence in a natural metapopulation.     

Low genetic diversity and strong but shallow population differentiation suggests genetic homogenization by metapopulation dynamics in a social spider

Mating systems and population dynamics influence genetic diversity and structure. Species that experience inbreeding and limited gene flow are expected to evolve isolated, divergent genetic lineages. Metapopulation dynamics with frequent extinctions and colonizations may, on the other hand, deplete and homogenize genetic variation, if extinction rate is sufficiently high compared to the effect of drift in local demes. We investigated these theoretical predictions empirically in social spiders that are highly inbred. Social spiders show intranest mating, female‐biased sex ratio, and frequent extinction and colonization events, factors that deplete genetic diversity within nests and populations and limit gene flow. We characterized population genetic structure in Stegodyphus sarasinorum, a social spider distributed across the Indian subcontinent. Species‐wide genetic diversity was estimated over approximately 2800 km from Sri Lanka to Himalayas, by sequencing 16 protein‐coding nuclear loci. We found 13 SNPs in 6592 bp (π = 0.00045) indicating low species‐wide nucleotide diversity. Three genetic lineages were strongly differentiated; however, only one fixed difference among them suggests recent divergence. This is consistent with a scenario of metapopulation dynamics that homogenizes genetic diversity across the species' range. Ultimately, low standing genetic variation may hamper a species' ability to track environmental change and render social inbreeding spiders ‘evolutionary dead‐ends’.     

Molecular diversity after a range expansion in heterogeneous environments

Recent range expansions have probably occurred in many species, as they often happen after speciation events, after ice ages, or after the introduction of invasive species. While it has been shown that range expansions lead to patterns of molecular diversity distinct from those of a pure demographic expansion, the fact that many species do live in heterogeneous environments has not been taken into account. We develop here a model of range expansion with a spatial heterogeneity of the environment, which is modeled as a gamma distribution of the carrying capacities of the demes. By allowing temporal variation of these carrying capacities, our model becomes a new metapopulation model linking ecological parameters to molecular diversity. We show by extensive simulations that environmental heterogeneity induces a loss of genetic diversity within demes and increases the degree of population differentiation. We find that metapopulations with low average densities are much more affected by environmental heterogeneity than metapopulations with high average densities, which are relatively insensitive to spatial and temporal variations of the environment. Spatial heterogeneity is shown to have a larger impact on genetic diversity than temporal heterogeneity. Overall, temporal heterogeneity and local extinctions are not found to leave any specific signature on molecular diversity that cannot be produced by spatial heterogeneity.     

Evolution of Pathogen Specialisation in a Host Metapopulation: Joint Effects of Host and Pathogen Dispersal

Metapopulation processes are important determinants of epidemiological and evolutionary dynamics in host-pathogen systems, and are therefore central to explaining observed patterns of disease or genetic diversity. In particular, the spatial scale of interactions between pathogens and their hosts is of primary importance because migration rates of one species can affect both spatial and temporal heterogeneity of selection on the other. In this study we developed a stochastic and discrete time simulation model to specifically examine the joint effects of host and pathogen dispersal on the evolution of pathogen specialisation in a spatially explicit metapopulation. We consider a plant-pathogen system in which the host metapopulation is composed of two plant genotypes. The pathogen is dispersed by air-borne spores on the host metapopulation. The pathogen population is characterised by a single life-history trait under selection, the infection efficacy. We found that restricted host dispersal can lead to high amount of pathogen diversity and that the extent of pathogen specialisation varied according to the spatial scale of host-pathogen dispersal. We also discuss the role of population asynchrony in determining pathogen evolutionary outcomes.     

Effects of cognitive abilities on metapopulation connectivity

Connectivity among demes in a metapopulation depends on both the landscape's and the focal organism's properties (including its mobility and cognitive abilities). Using individual‐based simulations, we contrast the consequences of three different cognitive strategies on several measures of metapopulation connectivity. Model animals search suitable habitat patches while dispersing through a model landscape made of cells varying in size, shape, attractiveness and friction. In the blind strategy, the next cell is chosen randomly among the adjacent ones. In the near‐sighted strategy, the choice depends on the relative attractiveness of these adjacent cells. In the far‐ sighted strategy, animals may additionally target suitable patches that appear within their perceptual range. Simulations show that the blind strategy provides the best overall connectivity, and results in balanced dispersal. The near‐sighted strategy traps animals into corridors that reduce the number of potential targets, thereby fragmenting metapopulations in several local clusters of demes, and inducing sink–source dynamics. This sort of local trapping is somewhat prevented in the far‐sighted strategy. The colonization success of strategies depends highly on initial energy reserves: blind does best when energy is high, near‐sighted wins at intermediate levels, and far‐sighted outcompetes its rivals at low energy reserves. We also expect strong effects in terms of metapopulation genetics: the blind strategy generates a migrant‐pool mode of dispersal that should erase local structures. By contrast, near‐ and far‐sighted strategies generate a propagule‐pool mode of dispersal and source–sink behavior that should boost structures (high genetic variance among‐ and low variance within local clusters of demes), particularly if metapopulation dynamics is also affected by extinction–colonization processes. Our results thus point to important effects of the cognitive ability of dispersers on the connectivity, dynamics and genetics of metapopulations.     

Trans-oceanic host dispersal explains high seabird tick diversity on Cape Verde islands

Parasites represent ideal models for unravelling biogeographic patterns and mechanisms of diversification on islands. Both host-mediated dispersal and within-island adaptation can shape parasite island assemblages. In this study, we examined patterns of genetic diversity and structure of Ornithodoros seabird ticks within the Cape Verde Archipelago in relation to their global phylogeography. Contrary to expectations, ticks from multiple, geographically distant clades mixed within the archipelago. Trans-oceanic colonization via host movements probably explains high local tick diversity, contrasting with previous research that suggests little large-scale dispersal in these birds. Although host specificity was not obvious at a global scale, host-associated genetic structure was found within Cape Verde colonies, indicating that post-colonization adaptation to specific hosts probably occurs. These results highlight the role of host metapopulation dynamics in the evolutionary ecology and epidemiology of avian parasites and pathogens.     

Population diversity in salmon: linkages among response,genetic and life history diversity

Response diversity and asynchrony are important for stability and resilience of meta‐populations, however little is known about the mechanisms that might drive such processes. In salmon populations, response diversity and asynchrony have been linked to the stability of their meta‐populations and the fisheries that integrate across them. We examined how population diversity influenced response diversity and asynchrony in 42 populations of Chinook salmon from the Fraser River, British Columbia. We examined diversity in the survival responses to large‐scale ocean climate variables for populations that differed in life history. Different life‐histories responded differently to ocean environmental conditions. For instance, an increase of offshore temperature was associated with decreased survival for a population with ocean rearing juveniles but increased survival for a population with stream rearing juveniles. In a second analysis, we examined asynchrony in abundance between populations, which we then correlated with life history, spatial, and genetic diversity. Populations that were more genetically distant had the most different population dynamics. Collectively, these results suggest that fine‐scale population diversity can contribute to the asynchrony and response diversity that underpins the stability of fisheries or metapopulation dynamics, and emphasize the need to manage and conserve this scale of population diversity.     

Contemporary gene flow and the spatio-temporal genetic structure of subdivided newt populations (Triturus cristatus, T. marmoratus)

Gene flow and drift shape the distribution of neutral genetic diversity in metapopulations, but their local rates are difficult to quantify. To identify gene flow between demes as distinct from individual migration, we present a modified Bayesian method to genetically test for descendants between an immigrant and a resident in a nonmigratory life stage. Applied to a metapopulation of pond-breeding European newts (Triturus cristatus, T. marmoratus) in western France, the evidence for gene flow was usually asymmetric and, for demes of known census size (N), translated into maximally seven reproducing immigrants. Temporal sampling also enabled the joint estimation of the effective demic population size (Ne) and the immigration rate m (including nonreproductive individuals). Ne ranged between 4.1 and 19.3 individuals, Ne/N ranged between 0.05 and 0.65 and always decreased with N; m was estimated as 0.19-0.63, and was possibly biased upwards. We discuss how genotypic data can reveal fine-scale demographic processes with important microevolutionary implications.     

The dynamics of two diffusively coupled predator-prey populations

I analyze the dynamics of predator and prey populations living in two patches. Within a patch the prey grow logistically and the predators have a Holling type II functional response. The two patches are coupled through predator migration. The system can be interpreted as a simple predator-prey metapopulation or as a spatially explicit predator-prey system. Asynchronous local dynamics are presumed by metapopulation theory. The main question I address is when synchronous and when asynchronous dynamics arise. Contrary to biological intuition, for very small migration rates the oscillations always synchronize. For intermediate migration rates the synchronous oscillations are unstable and I found periodic, quasi-periodic, and intermittently chaotic attractors with asynchronous dynamics. For large predator migration rates, attractors in the form of equilibria or limit cycles exist in which one of the patches contains no prey. The dynamical behavior of the system is described using bifurcation diagrams. The model shows that spatial predator-prey populations can be regulated through the interplay of local dynamics and migration.     

Demographic and genetic factors shaping contemporary metapopulation effective size and its empirical estimation in salmonid fish

The preservation of biodiversity requires an understanding of the maintenance of its components, including genetic diversity. Effective population size determines the amount of genetic variance maintained in populations, but its estimation can be complex, especially when populations are interconnected in a metapopulation. Theory predicts that the effective size of a metapopulation (meta-N(e)) can be decreased or increased by population subdivision, but little empirical work has evaluated these predictions. Here, we use neutral genetic markers and simulations to estimate the effective size of a putative metapopulation in Atlantic salmon (Salmo salar). For a weakly structured set of rivers, we find that meta-N(e) is similar to the sum of local deme sizes, whereas higher genetic differentiation among demes dramatically reduces meta-N(e) estimates. Interdemic demographic processes, such as asymmetrical gene flow, may explain this pattern. However, simulations also suggest that unrecognized population subdivision can also introduce downward bias into empirical estimation, emphasizing the importance of identifying the proper scale of distinct demographic and genetic processes. Under natural patterns of connectivity, evolutionary potential may generally be maintained at higher levels than the local population, with implications for conservation given ongoing species declines and habitat fragmentation.     

PHASE III OF WRIGHT'S SHIFTING BALANCE PROCESS AND THE VARIANCE AMONG DEMES IN MIGRATION RATE

Interdemic selection by the differential migration of individuals out from demes of high fitness and into demes of low fitness (Phase III) is one of the most controversial aspects of Wright's Shifting Balance Theory. I derive a relationship between Phase III migration and the interdemic selection differential, S, and show its potential effect on F_ST. The relationship reveals a diversifying effect of interdemic selection by Phase III migration on the genetic structure of a metapopulation. Using experimental metapopulations, I explored the effect of Phase III migration on F_ST by comparing the genetic variance among demes for two different patterns of migration: (1) island model migration and (2) Wright's Phase III migration. Although mean migration rates were the same, I found that the variance among demes in migration rate was significantly higher with Phase III than with island model migration. As a result, F_ST for the frequency of a neutral marker locus was higher with Phase III than it was with island model migration. By increasing F_ST, Phase III enhanced the genetic differentiation among demes for traits not subject to interdemic selection. This feature makes Wright's process different from individual selection which, by reducing effective population size, decreases the genetic variance within demes for all other traits. I discussed this finding in relation to the efficacy of Phase III and random migration for effecting peak shifts, and the contribution of genes with indirect effects to among‐deme variation.     

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.