The coalescence time of sampled genes in the structured coalescent model

The aim of this article is to investigate the distribution of the coalescence time (T) for sampled genes in the structured coalescent. We obtain some exact solutions for small samples and approximate distributions for n sampled genes in strong and weak migration. We also conduct computer simulation to evaluate efficiencies of these approximations and show the dependency of the distribution of the coalescence time on the geographical structure and the intensity of migration. In a panmictic population, we prove that the conditional distribution of the coalescence time given the number of segregating sites (S) among sampled genes is given by the weighted mean of the convolution of gamma distributions. We also study the joint distribution of T and S in the structured coalescent model and show some exact solutions.     

The strong-migration limit for the genealogical process in geographically structured populations

In this article we assume that the entire population is subdivided into a finite number of panmictic colonies, each of which consists of a respective number of haploid individuals. We also assume that random genetic drift occurs in each colony and migration among colonies, which is independent of time and ergodic. We study the genealogical process of sampled genes from geographically structured populations. We prove that if the actual total population number is replaced by the effective population number, the mean coalescence time converges to that in a panmictic population in the strong migration limit. We also obtain the geographical distribution of the common ancestor.     

Genealogy of Neutral Genes and Spreading of Selected Mutations in a Geographically Structured Population

In a geographically structured population, the interplay among gene migration, genetic drift and natural selection raises intriguing evolutionary problems, but the rigorous mathematical treatment is often very difficult. Therefore several approximate formulas were developed concerning the coalescence process of neutral genes and the fixation process of selected mutations in an island model, and their accuracy was examined by computer simulation. When migration is limited, the coalescence (or divergence) time for sampled neutral genes can be described by the convolution of exponential functions, as in a panmictic population, but it is determined mainly by migration rate and the number of demes from which the sample is taken. This time can be much longer than that in a panmictic population with the same number of breeding individuals. For a selected mutation, the spreading over the entire population was formulated as a birth and death process, in which the fixation probability within a deme plays a key role. With limited amounts of migration, even advantageous mutations take a large number of generations to spread. Furthermore, it is likely that these mutations which are temporarily fixed in some demes may be swamped out again by non-mutant immigrants from other demes unless selection is strong enough. These results are potentially useful for testing quantitatively various hypotheses that have been proposed for the origin of modern human populations.     

Geographical invariance and the strong-migration limit in subdivided populations

Invariance under population subdivision and the strong-migration limit are investigated for digenic samples in neutral models. The monoecious, diploid population is subdivided into a finite number of panmictic colonies that exchange gametes. The backward migration matrix is arbitrary, but time independent and ergodic (i.e., irreducible and aperiodic). Results are derived for the distribution of the place and time of coalescence, for the probability of identity in the model of infinitely many alleles, and for the distribution of the number of nucleotide differences in the model of infinitely many sites without recombination. Received: 5 October 1999 / Revised version: 1 February 2000 / Published online: 4 July 2000     

The number of segregating sites in a sample of DNA sequences from a geographically structured population

 The distribution of the number of segregating sites among randomly sampled DNA sequences from a geographically structured population is studied. We assume the infinitely-many-sites model of neutral genes and no recombination. Employing the genealogical process, we derive an equation for the generating function of the distribution of the number of segregating sites. First we study the strong-migration limit and prove that the distribution converges to that for a panmictic population. We also study the case of two sampled DNA sequences in the d-dimensional torus model with homogeneous migration. Received 13 July 1995; received in revised form 21 April 1997     

THE EFFECT OF COLLECTIVE DISPERSAL ON THE GENETIC STRUCTURE OF A SUBDIVIDED POPULATION

Correlated dispersal paths between two or more individuals are widespread across many taxa. The population genetic implications of this collective dispersal have received relatively little attention. Here we develop two‐sample coalescent theory that incorporates collective dispersal in a finite island model to predict expected coalescence times, genetic diversities, and F‐statistics. We show that collective dispersal reduces mixing in the system, which decreases expected coalescence times and increases F_ST. The effects are strongest in systems with high migration rates. Collective dispersal breaks the invariance of within‐deme coalescence times to migration rate, whatever the deme size. It can also cause F_ST to increase with migration rate because the ratio of within‐ to between‐deme coalescence times can decrease as migration rate approaches unity. This effect is most biologically relevant when deme size is small. We find qualitatively similar results for diploid and gametic dispersal. We also demonstrate with simulations and analytical theory the strong similarity between the effects of collective dispersal and anisotropic dispersal. These findings have implications for our understanding of the balance between drift–migration–mutation in models of neutral evolution. This has applied consequences for the interpretation of genetic structure (e.g., chaotic genetic patchiness) and estimation of migration rates from genetic data.     

Expected coalescence times and segregating sites in a model of glacial cycles

The climatic fluctuations of the Quaternary have influenced the distribution of numerous plant and animal species. Several species suffer population reduction and fragmentation, becoming restricted to refugia during glacial periods and expanding again during interglacials. The reduction in population size may reduce the effective population size, mean coalescence time and genetic variation, whereas an increased subdivision may have the opposite effect. To investigate these two opposing forces, we proposed a model in which a panmictic and a structured phase alternate, corresponding to interglacial and glacial periods. From this model, we derived an expression for the expected coalescence time and number of segregating sites for a pair of genes. We observed that increasing the number of demes or the duration of the structured phases causes an increase in coalescence time and expected levels of genetic variation. We compared numerical results with the ones expected for a panmictic population of constant size, and showed that the mean number of segregating sites can be greater in our model even when population size is much smaller in the structured phases. This points to the importance of population structure in the history of species subject to climatic fluctuations, and helps explain the long gene genealogies observed in several organisms.     

The distribution of the coalescence time and the number of pairwise nucleotide differences in a model of population divergence or speciation with an initial period of gene flow

This paper is concerned with a model of “isolation with an initial period of migration”, where a panmictic ancestral population split into n descendant populations which exchanged migrants symmetrically at a constant rate for a period of time and subsequently became completely isolated. In the limit as the population split occurred an infinitely long time ago, the model becomes an “isolation after migration” model, describing completely isolated descendant populations which arose from a subdivided ancestral population. The probability density function of the coalescence time of a pair of genes and the probability distribution of the number of pairwise nucleotide differences are derived for both models. Whilst these are theoretical results of interest in their own right, they also give an exact analytical expression for the likelihood, for data consisting of the numbers of nucleotide differences between pairs of DNA sequences where each pair is at a different, independent locus. The behaviour of the distribution of the number of pairwise nucleotide differences under these models is illustrated and compared to the corresponding distributions under the “isolation with migration” and “complete isolation” models. It is shown that the distribution of the number of nucleotide differences between a pair of DNA sequences from different descendant populations in the model of “isolation with an initial period of migration” can be quite different from that under the “isolation with migration model”, even if the average migration rate over time (and hence the total number of migrants) is the same in both scenarios. It is also illustrated how the results can be extended to other demographic scenarios that can be described by a combination of isolated panmictic populations and “symmetric island” models.     

Multilocus dynamics under haploid selection

 A general haploid selection model with arbitrary number of multiallelic loci and arbitrary linkage distribution is considered. The population is supposed to be panmictic. A dynamically equivalent diploid selection model is introduced. There is a position effect in this model if the original haploid selection is not multiplicative. If haploid selection is additive then the fundamental theorem is established even with an estimate for the change in the mean fitness. On this basis exponential convergence to an equilibrium is proved. As rule, the limit states are single-gamete ones. If, moreover, linkage is tight, then the single-gamete state with maximal fitness attracts the population for almost all initial states. Received 27 November 1995; received in revised form 17 January 1996     

Across a migratory divide: divergent migration directions and non‐breeding grounds of Eurasian reed warblers revealed by geolocators and stable isotopes

Migratory divides represent narrow zones of overlap between parapatric populations with distinct migration directions and, consequently, expected divergent non‐breeding distributions. The composition of the mixed population at a migratory divide and the corresponding non‐breeding ranges remain, however, unknown for many Palaearctic‐African migrants. Here, we used light‐level geolocation to track migration direction and non‐breeding grounds of Eurasian reed warblers Acrocephalus scirpaceus from three breeding populations across the species’ migratory divide. Moreover, by using feathers grown at non‐breeding grounds, we quantified stable isotope composition for individuals with known southwestern (SW) and southeastern (SE) migration directions. On a larger sample per population, we then assessed the proportions of SW‐ and SE‐migrating phenotypes in each of the three populations. All tracked reed warblers from Germany and two thirds of the birds tagged from the Czech population headed initially SW. Nevertheless, about one third of the birds from the Czech site migrated towards SE. No tracking data have been obtained for the Bulgarian population. The initial migration direction determined by geolocators was a strong predictor of the non‐breeding region, with SW migrants staying in west Africa and SE migrants in central Africa. Feather δ³⁴S and δ¹⁵N values confirmed the predominance of SW migrants in the German population, the co‐occurrence of SW and SE migrants in the Czech population, and indicated a high (72%) proportion of SE migrants in the Bulgarian population. Thus, the combined approach of geolocator tracking and stable isotopic assignments provided clear evidence for the existence of a migratory divide in the southeast of central Europe and predicted non‐breeding range in central and central‐eastern Africa for the eastern population.     

Estimating the minimum rate of genetic transformation in bacteria

Extensive data from multilocus electrophoresis are available for many bacterial populations. In some cases, for example Neisseria gonorrhoeae, these data are consistent with the population being in linkage equilibrium. This raises the following question. What frequency of transformation, or other means of genetic recombination, is needed, relative to mutation, to produce apparent panmixis? Simulation of a finite-population model suggests that, if transformation is at least twenty times as frequent as mutation, the population structure will be indistinguishable from a panmictic one, using the best available data sets. That is, relatively infrequent transformation is sufficient to produce approximate linkage equilibrium.     

Coalescence time for two genes from a subdivided population

In this paper a new form of the solution for the Laplace transform and moments of the distribution of the waiting time for two genes to coalescence is presented. The two genes are sampled from a subdivided population where migration rates between populations are constant in time. Equal subpopulation size is not assumed. For the special case of an island model with equal migration rates between islands, the Laplace transform of the coalescence time and the first and second moments are found explicitly. The new form of the solutions allows numerical calculation. The connection of how the results relate to a panmictic population when migration rates are large is illustrated using strong-migration-limit theory. Received: 19 April 1999 / Revised version: 22 March 2001 / Published online: 19 September 2001     

The Impact of Sampling Schemes on the Site Frequency Spectrum in Nonequilibrium Subdivided Populations

Using coalescent simulations, we study the impact of three different sampling schemes on patterns of neutral diversity in structured populations. Specifically, we are interested in two summary statistics based on the site frequency spectrum as a function of migration rate, demographic history of the entire substructured population (including timing and magnitude of specieswide expansions), and the sampling scheme. Using simulations implementing both finite-island and two-dimensional stepping-stone spatial structure, we demonstrate strong effects of the sampling scheme on Tajima's D (D_T) and Fu and Li's D (D_FL) statistics, particularly under specieswide (range) expansions. Pooled samples yield average D_T and D_FL values that are generally intermediate between those of local and scattered samples. Local samples (and to a lesser extent, pooled samples) are influenced by local, rapid coalescence events in the underlying coalescent process. These processes result in lower proportions of external branch lengths and hence lower proportions of singletons, explaining our finding that the sampling scheme affects D_FL more than it does D_T. Under specieswide expansion scenarios, these effects of spatial sampling may persist up to very high levels of gene flow (Nm > 25), implying that local samples cannot be regarded as being drawn from a panmictic population. Importantly, many data sets on humans, Drosophila, and plants contain signatures of specieswide expansions and effects of sampling scheme that are predicted by our simulation results. This suggests that validating the assumption of panmixia is crucial if robust demographic inferences are to be made from local or pooled samples. However, future studies should consider adopting a framework that explicitly accounts for the genealogical effects of population subdivision and empirical sampling schemes.     

The two-locus ancestral graph in a subdivided population: convergence as the number of demes grows in the island model

We study the ancestral recombination graph for a pair of sites in a geographically structured population. In particular, we consider the limiting behavior of the graph, under Wrights island model, as the number of subpopulations, or demes, goes to infinity. After an instantaneous sample-size adjustment, the graph becomes identical to the two-locus graph in an unstructured population, but with a time scale that depends on the migration rate and the deme size. Interestingly, when migration is gametic, this rescaling of time increases the population mutation rate but does not affect the population recombination rate. We compare this to the case of a partially-selfing population, in which both mutation and recombination depend on the selfing rate. Our result for gametic migration holds both for finite-sized demes, and in the limit as the deme size goes to infinity. However, when migration occurs during the diploid phase of the life cycle and demes are finite in size, the population recombination rate does depend on the migration rate, in a way that is reminiscent of partial selfing. Simulations imply that convergence to a rescaled panmictic ancestral recombination graph occurs for any number of sites as the number of demes approaches infinity.Send offprint request to: Sabin LessardS. Lessard was supported by grants from the Natural Sciences and Research Council of Canada, the Fonds Québécois de la Recherche sur la Nature et les Technologies, and the Université de Montréal.J. Wakeley was supported by a Career Award (DEB-0133760) and by a grant (DEB-9815367) from the National Science Foundation.     

The diffusion model for migration and selection in a dioecious population

The diffusion approximation is derived for migration and selection at a multiallelic locus in a dioecious population subdivided into a lattice of panmictic colonies. Generations are discrete and nonoverlapping; autosomal and X-linked loci are analyzed. The relation between juvenile and adult subpopulation numbers is very general and includes both soft and hard selection; the zygotic sex ratio is the same in every colony. All the results hold for both adult and juvenile migration. If ploidy-weighted average selection, drift, and diffusion coefficients are used, then the ploidy-weighted average allelic frequencies satisfy the corresponding partial differential equation for a monoecious population. The boundary conditions and the unidimensional transition conditions for coincident discontinuities in the carrying capacity and migration rate extend identically. The previous unidimensional formulation and analysis of symmetric, nearest-neighbor migration of a monoecious population across a geographical barrier is generalized to symmetric migration of arbitrary finite range, and the transition conditions are shown to hold for a dioecious population. Thus, the entire theory of clines and of the wave of advance of favorable alleles is applicable to dioecious populations.This work was supported by National Science Foundation grant BSR-9006285     

NEUTRAL THEORY OF QUANTITATIVE GENETIC VARIANCE IN AN ISLAND MODEL WITH LOCAL EXTINCTION AND COLONIZATION

Additive genetic variance maintained by mutation in a selectively neutral quantitative character is analyzed for an ideal population distributed on n islands, each with local effective size N, that exchange migrants at a small rate, m. In a stable population structure, the expected genetic variance maintained within islands is identical to that in a panmictic population of the same total size, regardless of the migration rate (m > 0). This result contrasts with Wright's classical conclusion, based on inbreeding coefficients, that at least one immigrant per island every other generation (Nm > ½) is necessary for the genetic variance within local populations to approach that under panmixia. The expected genetic variance maintained among islands is inversely proportional to m and increases with the number of islands, but is independent of N. Local extinction and colonization diminish the genetic variance maintained within islands by reducing the effective size of island populations through the founder effect, although the expected genetic variance within islands is nearly as large as that in a panmictic population of the same total effective size. If the founders of new colonies originate from more than one island, rates of local extinction and colonization larger than about twice the migration rate will substantially reduce the genetic variance maintained among islands. These results indicate the importance of mutation and migration in maintaining quantitative genetic variance within small local populations.     

The strong-migration limit in geographically structured populations

Summary Some strong-migration limits are established for geographically structured populations. A diploid monoecious population is subdivided into a finite number of colonies, which exchange migrants. The migration pattern is fixed and ergodic, but otherwise arbitrary. Generations are discrete and nonoverlapping; the analysis is restricted to a single locus. In all the limiting results, an effective population number N _e ( N _T ) appears instead of the actual total population number N _T . 1. If there is no selection, every allele mutates at rate u to types not preexisting in the population, and the (finite) subpopulation numbers N _i are very large, then the ultimate rate and pattern of convergence of the probabilities of allelic identity are approximately the same as for panmixia. If, in addition, the N _i are proportional to 1/u, as N _T 8, the equilibrium probabilities of identity converge to the panmictic value. 2. With a finite number of alleles, any mutation pattern, an arbitrary selection scheme for each colony, and the mutation rates and selection coefficients proportional to 1/N _T , let P _j be the frequency of the allele A _j in the entire population, averaged with respect to the stationary distribution of the backward migration matrix M. As N _T 8, the deviations of the allelic frequencies in each of the subpopulations from P _j converge to zero; the usual panmictic mutation-selection diffusion is obtained for P _j , with the selection intensities averaged with respect to the stationary distribution of M. In both models, N _e = N _T and all effects of population subdivision disappear in the limit if, and only if, migration does not alter the subpopulation numbers.Supported by the National Science Foundation (Grant No. DEB77-21494)     

The distribution of the coalescence time and the number of pairwise nucleotide differences in the "isolation with migration" model

This paper is concerned with the “isolation with migration” model, where a panmictic ancestral population gave rise to a symmetric n-island model, time τ ago. Explicit analytical expressions are derived for the probability density function of the coalescence time of a pair of genes sampled at random from the same subpopulation or from different subpopulations, and for the probability distribution of the number of pairwise nucleotide differences.     

Postglacial history of a widespread conifer produces inverse clines in selective neutrality tests

Deviations of the site frequency spectrum of mutations (SFS) from neutral expectations may be caused by natural selection or by demographic processes such as population subdivision or temporal changes in population size. As most widespread temperate and boreal tree species have expanded from glacial refugia in the past 13 000 years, colonization bottlenecks associated with this migration may have left variable demographic signatures among geographic populations corresponding to distance from the refugia. To determine whether the signature of postglacial re‐colonization has skewed the SFS in the widely distributed conifer Sitka spruce (Picea sitchensis (Bong.) Carr.), we re‐sequenced 153 nuclear genes in six populations from across the species range. We found that while the SFS for the pooled sample produced negative values for Tajima’s D and Fay and Wu’s H, these statistics exhibited strong clinal variation when populations were analysed separately (R² = 0.84, P = 0.007 for Tajima’s D and R² = 0.65, P = 0.033 for Fay and Wu’s H). When historical bottlenecks of varying age were simulated using approximate Bayesian computation, distance of populations from the southern range limit explained most of the variation in bottleneck timing among populations (R² = 0.89, P = 0.003). These data suggest that sequential population bottlenecks during postglacial re‐colonization have resulted in diverse among‐population signatures within the contemporary SFS in Sitka spruce, with rare variants more common in the south, and medium‐frequency variants more common in the north. Our results also emphasize the need to consider sampling strategy and to explore population‐specific null demographic models in surveys of nucleotide variation in widely distributed species.