Genomic diversity,population structure,and migration following rapid range expansion in the Balsam Poplar,Populus balsamifera

Rapid range expansions can cause pervasive changes in the genetic diversity and structure of populations. The postglacial history of the Balsam Poplar, Populus balsamifera, involved the colonization of most of northern North America, an area largely covered by continental ice sheets during the last glacial maximum. To characterize how this expansion shaped genomic diversity within and among populations, we developed 412 SNP markers that we assayed for a range‐wide sample of 474 individuals sampled from 34 populations. We complemented the SNP data set with DNA sequence data from 11 nuclear loci from 94 individuals, and used coalescent analyses to estimate historical population size, demographic growth, and patterns of migration. Bayesian clustering identified three geographically separated demes found in the Northern, Central, and Eastern portions of the species’ range. These demes varied significantly in nucleotide diversity, the abundance of private polymorphisms, and population substructure. Most measures supported the Central deme as descended from the primary refuge of diversity. Both SNPs and sequence data suggested recent population growth, and coalescent analyses of historical migration suggested a massive expansion from the Centre to the North and East. Collectively, these data demonstrate the strong influence that range expansions exert on genomic diversity, both within local populations and across the range. Our results suggest that an in‐depth knowledge of nucleotide diversity following expansion requires sampling within multiple populations, and highlight the utility of combining insights from different data types in population genomic studies.     

The effects of subdivision on the genetic divergence of populations and species

Abstract. An island model of migration is used to study the effects of subdivision within populations and species on sample genealogies and on between-population or between-species measures of genetic variation. The model assumes that the number of demes within each population or species is large. When populations (or species), connected either by gene flow or historical association, are themselves subdivided into demes, changes in the migration rate among demes alter both the structure of genealogies and the time scale of the coalescent process. The time scale of the coalescent is related to the effective size of the population, which depends on the migration rate among demes. When the migration rate among demes within populations is low, isolation (or speciation) events seem more recent and migration rates among populations seem higher because the effective size of each population is increased. This affects the probability of reciprocal monophyly of two samples, the chance that a gene tree of a sample matches the species tree, and relative likelihoods of different types of polymorphic sites. It can also have a profound effect on the estimation of divergence times.     

Gene genealogies in a metapopulation

A simple genealogical process is found for samples from a metapopulation, which is a population that is subdivided into a large number of demes, each of which is subject to extinction and recolonization and receives migrants from other demes. As in the migration-only models studied previously, the genealogy of any sample includes two phases: a brief sample-size adjustment followed by a coalescent process that dominates the history. This result will hold for metapopulations that are composed of a large number of demes. It is robust to the details of population structure, as long as the number of possible source demes of migrants and colonists for each deme is large. Analytic predictions about levels of genetic variation are possible, and results for average numbers of pairwise differences within and between demes are given. Further analysis of the expected number of segregating sites in a sample from a single deme illustrates some previously known differences between migration and extinction/recolonization. The ancestral process is also amenable to computer simulation. Simulation results show that migration and extinction/recolonization have very different effects on the site-frequency distribution in a sample from a single deme. Migration can cause a U-shaped site-frequency distribution, which is qualitatively similar to the pattern reported recently for positive selection. Extinction and recolonization, in contrast, can produce a mode in the site-frequency distribution at intermediate frequencies, even in a sample from a single deme.     

A first step towards inferring levels of long-distance dispersal during past expansions

Improving the realism of spatially explicit demographic models is important for better inferring the history of past populations and for understanding the genetic bases of adaptation and speciation. One particular type of demographic event to take into account is long-distance dispersal (LDD). The goals of this study are to explore the impact of various levels of LDD on genetic diversity and to show to what extent LDD levels can be correctly inferred from multilocus data sets using an approximate Bayesian computation approach. We therefore incorporated LDD into a 2D stepping stone forward simulation framework coupled to a coalescent backward simulation step to generate genetic diversity at 100 microsatellite markers under various demographic conditions relevant to recent human evolution. Our results confirm that LDD considerably increases genetic diversity within demes and decreases levels of diversity between demes. By controlling the spatial occurrence of LDD, it appears that LDD events occurring during a phase of range expansion into new territories are more important in maintaining genetic diversity than those occurring in the wake of the expansion or when colonization is over. We also show that it is possible to infer whether LDD has occurred during a range expansion, but our results suggest that one can only approximately estimate the extent of LDD based on genetic summary statistics.     

F(st) in a Hierarchical Island Model

It is shown that in a hierarchical island model, in which demes within a neighborhood exchange migrants at a much higher rate than do demes in different neighborhoods, hierarchical F statistics introduced by S. Wright can indicate the extent of gene flow within and between neighborhoods. At equilibrium, the within-neighborhood inbreeding coefficient, FSN, is approximately 1/(1 + 4Nm1) where N is the deme size and m1 is the migration rate among demes in the same neighborhood. The between-neighborhood inbreeding coefficient, FNT, is approximately 1/(1 + 4Ndm2) where d is the number of demes in a neighborhood and m2 is the migration rate among demes in different neighborhoods.     

AQUASPLATCHE: a program to simulate genetic diversity in populations living in linear habitats

Classical models of structured populations do not apply well to species leaving in semilinear habitats such as freshwater fishes, since these habitats are intermediate between one‐dimensional and two‐dimensional stepping‐stone models. In order to investigate the genetic diversity of such populations, we have developed a new simulation program called aquasplatche . It starts by dividing a user‐defined vectorized network into segments of arbitrary length, each segment hosting a single deme. The program then proceeds by simulating the colonization of the environment from an arbitrary source, recording the evolution of the deme densities and the migration events between adjacent demes over time. This demographic history is then used to generate genetic data of population samples located in various segments of the network, using a backward coalescent framework. Different versions of aquasplatche are freely available on http://cmpg.unibe.ch/software/aquasplatche .     

The coalescent in an island model of population subdivision with variation among demes

A simple genealogical structure is found for a general finite island model of population subdivision. The model allows for variation in the sizes of demes, in contributions to the migrant pool, and in the fraction of each deme that is replaced by migrants every generation. The ancestry of a sample of non-recombining DNA sequences has a simple structure when the sample size is much smaller than the total number of demes in the population. This allows an expression for the probability distribution of the number of segregating sites in the sample to be derived under the infinite-sites mutation model. It also yields easily computed estimators of the migration parameter for each deme in a multi-deme sample. The genealogical process is such that the lineages ancestral to the sample tend to accumulate in demes with low migration rates and/or which contribute disproportionately to the migrant pool. In addition, common ancestor or coalescent events tend to occur in demes of small size. This provides a framework for understanding the determinants of the effective size of the population, and leads to an expression for the probability that the root of a genealogy occurs in a particular geographic region, or among a particular set of demes.     

Historical range expansion determines the phylogenetic diversity introduced during contemporary species invasion

For a species rapidly expanding its geographic range, such as during biological invasion, most alleles in the introduced range will have their evolutionary origins in the native range. Yet, the way in which historical processes occurring over evolutionary time in the native range contribute to the diversity sampled during contemporary invasion is largely unknown. We used chloroplast DNA (cpDNA) gene genealogies and coalescent methods to study two congeneric plants, Silene latifolia and S. vulgaris. We examined how phylogenetic diversity was shaped by demographic growth and historical range expansions in the native European range, and how this history affected the diversity sampled during their recent invasion of North America. Genealogies from both species depart from neutrality, likely as a result of demographic expansion in the ancestral range, the timing of which corresponds to shortly after each species originated. However, the species differ in the spatial distribution of cpDNA lineages across the native range. Silene latifolia shows a highly significant phylogeographic structure that most likely reflects different avenues of the post-glacial expansion into northern Europe from Mediterranean refugia. By contrast, cpDNA lineages in S. vulgaris have been widely scattered across Europe during, or since, the most recent post-glacial expansion. These different evolutionary histories resulted in dramatic differences in how phylogenetic diversity was sampled during invasion of North America. In S. latifolia, relatively few, discrete invasion events from a structured native range resulted in a rather severe genetic bottleneck, but also opportunities for admixture among previously isolated lineages. In S. vulgaris, lack of genetic structure was accompanied by more representative sampling of phylogenetic diversity during invasion, and reduced potential for admixture. Our results provide clear insights into how historical processes may feed forward to influence the phylogenetic diversity of species invading new geographic ranges.     

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.     

Spatial environmental variation can select for evolvability

Previous studies have shown that temporally fluctuating environments can create indirect selection for modifiers of evolvability. Here, we use a simple computational model to investigate whether spatially varying environments (multiple demes with limited migration among them, and a different, static selective optimum in each) can also create indirect selection for increased evolvability. The answer is surprisingly complicated. Spatial variation in the environment can sharply reduce the survival rate of migrants, because migrants may be maladapted to their new deme, relative to incumbents. The incumbent advantage can be removed by occasional extinctions in single demes. After all incumbents in a particular deme die, incoming migrants from other demes will, on average, be similarly maladapted to the new environment. This sets off a race to adapt rapidly. Over many extinction events, and the subsequent invasions by maladapted immigrants into a new environment, indirect selection for the ability to adapt rapidly, also known as high evolvability, may result.     

Effect of population structure on the amount of polymorphism and the fixation probability under overdominant selection

Under overdominant selection, mutants substantially contribute to increase the amount of polymorphism. It is also known that under neutrality as the migration rates among demes decrease in a subdivided population, the amount of polymorphism increases along with the increase of the effective population size, N(e). In this study, under overdominant selection the effect of population subdivision on the amount of polymorphism was investigated using the diffusion approximation and the low migration approximation. It was shown that if selection is medium or strong (e.g., N(T)s > 1, where N(T) is the population size and s is the selective advantage of heterozygotes), the nucleotide diversity, pi, decreases along with the decrease of Nm against the increase of N(e), where N is the size of demes and m is the migration rate per deme. In addition, the ratio of the nucleotide diversity to the evolutionary rate also decreases along with the decrease of Nm. In some cases the ratio becomes smaller than that expected under neutrality as Nm decreases.     

Fixation probability for a beneficial allele and a mutant strategy in a linear game under weak selection in a finite island model

The effect of population structure on the probability of fixation of a newly introduced mutant under weak selection is studied using a coalescent approach. Wright's island model in a framework of a finite number of demes is assumed and two selection regimes are considered: a beneficial allele model and a linear game among offspring. A first-order approximation of the fixation probability for a single mutant with respect to the intensity of selection is deduced. The approximation requires the calculation of expected coalescence times, under neutrality, for lineages starting from two or three sampled individuals. The results are obtained in a general setting without assumptions on the number of demes, the deme size or the migration rate, which allows for simultaneous coalescence or migration events in the genealogy of the sampled individuals. Comparisons are made with limit cases as the deme size or the number of demes goes to infinity or the migration rate goes to zero for which a diffusion approximation approach is possible. Conditions for selection to favor a mutant strategy replacing a resident strategy in the context of a linear game in a finite island population are addressed.     

The redivision of geographic population and genetic structure of Eriocheir in the West-Pacific Ocean

Mitten crabs (Eriocheir) were widely distributed in the West-Pacific Ocean. In response to change of temperature during the Pleistocene, they were forced into repeated cycle of glacial retreat and interglacial expansion. In the study, genetic diversity and phylogeographic structure of Eriocheir were analysed by mtDNA markers. All results supported the redivision of four geographic populations of mitten crabs corresponding to geographic distribution. And introgression of the sympatric populations might have happened for some shared haplotypes. Results of "star-like" network, unimodal curve of mismatch distributions and significantly negative values of neutral test suggested a recent population expansion.     

Coalescent simulations of Yakut mtDNA variation suggest small founding population

The Yakuts are a Turkic‐speaking population from northeastern Siberia who are believed to have originated from ancient Turkic populations in South Siberia, based on archaeological and ethnohistorical evidence. In order to better understand Yakut origins, we modeled 25 demographic scenarios and tested by coalescent simulation whether any are consistent with the patterns of mtDNA diversity observed in present‐day Yakuts. The models consist of either two simulated demes that represent Yakuts and a South Siberian ancestral population, or three demes that also include a regional Northeast Siberian population that served as a source of local gene flow into the Yakut deme. The model that produced the best fit to the observed data defined a founder group with an effective female population size of only 150 individuals that migrated northwards approximately 1,000 years BP and who experienced significant admixture with neighboring populations in Northeastern Siberia. These simulation results indicate a pronounced founder effect that was primarily kin‐structured and reconcile reported discrepancies between Yakut mtDNA and Y chromosome diversity levels. Am J Phys Anthropol, 2009. © 2009 Wiley‐Liss, Inc.     

Allopatry increases the balance of phylogenetic trees during radiation under neutral speciation

The shape of a phylogenetic tree is defined by the sequence of speciation events, represented by its branching points, and extinctions, represented by branch interruptions. In a neutral scenario of parapatry and isolation by distance, species tend to branch off the original population one after the other, leading to highly unbalanced trees. In this case the degree of imbalance, measured by the normalized Sackin index, grows linearly with species richness. Here we claim that moderate values of imbalance for trees with large number of species can occur if the geographic distribution involves more than one deme (allopatry) and speciation is parapatric within demes. The combined values of balance (normalized Sackin index) and species richness provide an estimate of how many demes were involved in the process if it happened in such neutral scenario. We also show that the spatial division in demes moderately slows down the diversification process, portraying a neutral mechanism for structuring the branch length distribution of phylogenetic trees.     

Postglacial population expansion of Dacrydium pectinatum(Podocarpaceae)in Hainan,southern China,based on cpDNA tmL-F nocoding sequence data

This study determined the sequences of chloro plast DNA(cpDNA)trnL-F non-coding regions of indiViduals of a tropical coniferous species,Dacrydium pectinatum,collected from 12 natural populations located in Hainan Province,southern China.Sequence length varied from 868 bp to 876 bp,indicating length polymorphism.Base composition in the sequences was high in A+T content between 64.17% and 64.95%.and no recombination event occurred (Rm=0).Thirty haplotypes were identified based on statistical parsimony algorithm by running the TCS program.Populations of D.pectinatum in Hainan were lacking genetic differentiation.Such a deduction was supported by the observed FST values(0.00),AMOVA(24.17% of molecular variance attributed to difference among populations,P＞0.05),high values of Nm(ranging from 1.92 to 2.50)and the branching structure in neighbor-joining(NJ)tree conin the TCS network of trnL-F haplotypes,and majority of the haplotypes coalesced near the tips in NJ tree.Gene genealogies of cpDNA haplotypes proposed a recent population expansion of D.pectinatum in Hainan,which was further supported by the results from Tajima's D test and mismatch distribution analysis.Our data.in COnjunction with geological and palynological evidences,showed that in the Holocene,due to global warming,refugee populations of D.pectinatum in Hainan might experience a range expansion.     

Populational heritability: empirical studies of evolution in metapopulations

Using demes from experimental metapopulations of the flour beetle, Tribolium castaneum, we investigated phase 3 of Wright's shifting balance process. Using parent demes of high, intermediate, and low mean fitness, we experimentally modeled migration of varying amounts from demes of high mean fitness into demes of lower mean fitness (like phase 3) as well as the reciprocal (the opposite of phase 3). In natural populations, some migration among demes occurs independently of deme fitness. In this case, demes of high mean fitness are likely to receive migrants from demes of lower mean fitness; these effects might limit the effectiveness of phase 3 but have not been studied experimentally. We estimated the populational heritability of mean fitness by the regression of offspring deme means on the weighted parental means and found moderate levels of demic heritability one (0.641-0.690) and two (0.518-0.552) generations after migration. We discuss our findings in relation to the role of interdemic migration in "adaptive peak shifts" in metapopulations and the controversies over group selection and the units of inheritance.     

A populationwide coalescent analysis of Icelandic matrilineal and patrilineal genealogies: evidence for a faster evolutionary rate of mtDNA lineages than Y chromosomes

Historical inferences from genetic data increasingly depend on assumptions about the genealogical process that shapes the frequencies of alleles over time. Yet little is known about the structure of human genealogies over long periods of time and how they depart from expectations of standard demographic models, such as that attributed to Wright and Fisher. To obtain such information and to examine the recent evolutionary history of mtDNA and Y-chromosome haplotypes in the Icelandic gene pool, we traced the matrilineal and patrilineal ancestry of all 131,060 Icelanders born after 1972 back to two cohorts of ancestors, one born between 1848 and 1892 and the other between 1798 and 1742. This populationwide coalescent analysis of Icelandic genealogies revealed highly positively skewed distributions of descendants to ancestors, with the vast majority of potential ancestors contributing one or no descendants and a minority of ancestors contributing large numbers of descendants. The expansion and loss of matrilines and patrilines has caused considerable fluctuation in the frequencies of mtDNA and Y-chromosome haplotypes, despite a rapid population expansion in Iceland during the past 300 years. Contrary to a widespread assumption, the rate of evolution caused by this lineage-sorting process was markedly faster in matrilines (mtDNA) than in patrilines (Y chromosomes). The primary cause is a 10% shorter matrilineal generation interval. Variance in the number of offspring produced within each generation was not an important differentiating factor. We observed an intergenerational correlation in offspring number and in the length of generation intervals in the matrilineal and patrilineal genealogies, which was stronger in matrilines and thus contributes to their faster evolutionary rate. These findings may have implications for coalescent date estimates based on mtDNA and Y chromosomes.