Genetic differentiation within and between two habitats.

We investigate the usefulness of analyses of population differentiation between different ecological types, such as host races of parasites or sources and sink habitats. To that aim, we formulate a model of population structure involving two classes of subpopulations found in sympatry. Extensions of previous results for Wright''s F-statistics in island and isolation-by-distance models of dispersal are given. It is then shown that source and sinks cannot in general be distinguished by F-statistics nor by their gene diversities. The excess differentiation between two partially isolated classes with respect to differentiation within classes is shown to decrease with distance, and for a wide range of parameter values it should be difficult to detect. In the same circumstances little differentiation will be observed in "hierarchical" analyses between pools of samples from each habitat, and differences between levels of differentiation within each habitat will only reflect differences between levels of gene diversity within each habitat. Exceptions will indicate strong isolation between the different classes or habitat-related divergent selection.     

THE RELEVANCE OF GENE FLOW IN METAPOPULATION DYNAMICS OF AN OCEANIC ISLAND ENDEMIC,OLEA EUROPAEA SUBSP. GUANCHICA

Theoretical and empirical studies suggest that geographical isolation and extinction–recolonization dynamics are two factors causing strong genetic structure in metapopulations, but their consequences in species with high dispersal abilities have not been tested at large scales. Here, we investigated the effect of population age structure and isolation by distance in the patterns of genetic diversity in a wind‐pollinated, zoochorous tree (Olea europaea subsp. guanchica) sporadically affected by volcanic events across the Canarian archipelago. Genetic variation was assessed at six nuclear microsatellites (nDNA) and six chloroplast fragments (cpDNA) in nine subpopulations sampled on four oceanic islands. Subpopulations occurring on more recent substrates were more differentiated than those on older substrates, but within‐subpopulation genetic diversity was not significantly different between age groups for any type of marker. Isolation‐by‐distance differentiation was observed for nDNA but not for cpDNA, in agreement with other metapopulation studies. Contrary to the general trend for island systems, between‐island differentiation was extremely low, and lower than differentiation between subpopulations on the same island. The pollen‐to‐seed ratio was close to one, two orders of magnitude lower than the average estimated for other wind‐pollinated, animal‐dispersed plants. Our results showed that population turnover and geographical isolation increased genetic differentiation relative to an island model at equilibrium, but overall genetic structure was unexpectedly weak for a species distributed among islands. This empirical study shows that extensive gene flow, particularly mediated by seeds, can ameliorate population subdivision resulting from extinction–recolonization dynamics and isolation by distance.     

Gene Flow and Genetic Differentiation

A brief analysis is presented for the effects of gene flow upon genetic differentiation within and between populations generated by mutation and drift. Previous results obtained with the "island" model are developed into a form that lends itself to biological interpretation. Attention is focused upon the effective local population size and the ratio of the genetic identity of two genes in different populations to that of two genes in the same population. The biological significance of this ratio, which is independent of population size, is discussed. Similarities between the results of this model and those of the "stepping-stone" model are noted.     

Patterning of geographic variation in Middle Pleistocene Homo frontal bone morphology

A quantitative assessment of the frontal bone morphology of a sample of Middle Pleistocene hominins was undertaken in order to address questions regarding their population structure and evolutionary history. Outline tracings of the frontal bones of forty-seven fossil crania were obtained, and size-standardized measurements were then computed using an Elliptical Fourier analysis of these tracings. Principal component scores of the Fourier harmonic amplitudes were derived and served as a quantitative representation of the morphology of the frontal bone. Morphological, geographical, and temporal distance matrices were then constructed between each pair of fossils. A partial Mantel matrix correlation test was performed between morphological and geographical distance matrices, controlling for temporal distance, in order to determine if the pattern of geographical differentiation in features of the frontal bone of mid-Pleistocene Homo followed that of an isolation-by-distance model of population structure. The results of the partial Mantel tests indicate that the overall patterning of differentiation in the features of the frontal bone cannot best be explained by a population structure shaped by isolation-by-distance. Additionally, various aspects of the frontal bone quantified here follow different patterns of geographical differentiation, suggesting that a mosaic pattern of evolution holds true for characters within one cranial region and not just for those between regions.     

Ancestral processes in population genetics-the coalescent

A special stochastic process, called the coalescent, is of fundamental interest in population genetics. For a large class of population models this process is the appropriate tool to analyse the ancestral structure of a sample of n individuals or genes, if the total number of individuals in the population is sufficiently large. A corresponding convergence theorem was first proved by Kingman in 1982 for the Wright-Fisher model and the Moran model. Generalizations to a large class of exchangeable population models and to models with overlying mutation processes followed shortly later. One speaks of the "robustness of the coalescent, as this process appears in many models as the total population size tends to infinity. This publication can be considered as an introduction to the theory of the coalescent as well as a review of the most important "convergence-to-the-coalescent-theorems. Convergence theorems are not only presented for the classical exchangeable haploid case but also for larger classes of population models, for example for diploid, two-sex or non-exchangeable models. A review-like summary of further examples and applications of convergence to the coalescent is given including the most important biological forces like mutation, recombination and selection. The general coalescent process allows for simultaneous multiple mergers of ancestral lines.     

Rapid genetic divergence in postglacial populations of threespine stickleback (Gasterosteus aculeatus): the role of habitat type, drainage and geographical proximity

The goal of this study was to evaluate the role of common ancestry, and of geographical or reproductive isolation, in genetic divergence in populations of threespine sticklebacks (Gasterosteus aculeatus). Using seven DNA microsatellite loci we compared the effects of habitat type, drainage system and geographical proximity on genetic distance among 16 populations situated in an area in Schleswig-Holstein (Germany) that became deglaciated approximately 12 000 years ago. Stickleback population structure correlated only weakly with drainage system, whereas the primary divergence was among habitat types. Phylogenetic analysis revealed that lake (n = 7) and river (n = 5) populations formed two distinct clades (Cavalli-Sforza's and Edwards' chord distance, 82-100% bootstrap support) at approximately equal genetic distances to a third clade, comprising putative estuarine (n = 4) ancestors. Allele frequencies in lake and river populations represented different subsets of the genetically more diverse estuarine populations. In nested amovas approximately twice the genetic variance was distributed among lake vs. river vs. estuarine populations as compared with the combined effects of drainage system and geographical distance. Limited gene flow between habitat types must have been established after postglacial colonization, suggesting ecological hybrid inferiority or behavioural mating barriers between ecotypes. Within estuarine and lake populations, population differentiation followed an isolation-by-distance model. Given the high observed heterozygosities within the 16 study populations (HO = 0.65-0.87), the mean divergence between lake and river population pairs (FST = 0.18 +/- 0.007) would be reached after 300-6000 generations in a stepwise mutation model, depending on the size of N(e). This demonstrates both the utility of hypervariable microsatellites for detecting recent population divergences and the danger of operating at temporal or spatial scales which are beyond their resolution.     

Ecological correlates of population genetic structure: a comparative approach using a vertebrate metacommunity

Identifying ecological factors associated with population genetic differentiation is important for understanding microevolutionary processes and guiding the management of threatened populations. We identified ecological correlates of several population genetic parameters for three interacting species (two garter snakes and an anuran) that occupy a common landscape. Using multiple regression analysis, we found that species interactions were more important in explaining variation in population genetic parameters than habitat and nearest-neighbour characteristics. Effective population size was best explained by census size, while migration was associated with differences in species abundance. In contrast, genetic distance was poorly explained by the ecological correlates that we tested, but geographical distance was prominent in models for all species. We found substantially different population dynamics for the prey species relative to the two predators, characterized by larger effective sizes, lower gene flow and a state of migration-drift equilibrium. We also identified an escarpment formed by a series of block faults that serves as a barrier to dispersal for the predators. Our results suggest that successful landscape-level management should incorporate genetic and ecological data for all relevant species, because even closely associated species can exhibit very different population genetic dynamics on the same landscape.     

Quantifying gene flow from spatial genetic structure data in a metapopulation of Chamaecrista fasciculata (Leguminosae)

Abstract An extensive allozyme survey was conducted within a natural "meta" population of the native North American annual legume, Chamaecrista fasciculata (Leguminosae) to quantify genetic structure at different spatial scales. Gene flow was then estimated by a recently developed indirect method based on a continuous population model, using pairwise kinship coefficients between individuals. The indirect estimates of gene flow, quantified in terms of neighborhood size, with an average value on the order of 150 individuals, were concordant among different spatial scales (subpopulation, population, metapopulation). This gene-flow value lies within the range of direct estimates previously documented from observations of pollen and seed dispersal for the same metapopulation. Monte Carlo simulations using the direct measures of gene flow as parameters further demonstrated that the observed spatial pattern of allozyme variation was congruent with a model of isolation by distance. Combining previously published estimates of pollen dispersal distances with kinship coefficients from this study, we quantified biparental inbreeding relative to either a single subpopulation or the whole metapopulation. At the level of a neighborhood, little biparental inbreeding was observed and most departure from Hardy-Weinberg genotypic proportions was explained by self-fertilization, whereas both selfing and biparental inbreeding contributed to nonrandom mating at the metapopulation level. Gene flow was also estimated from indirect methods based on a discontinuous population structure model. We discuss these results with respect to the effect of a patchy population structure on estimation of gene flow.     

Dissimilarity of individual microsatellite profiles under different mutation models: Empirical approach

Microsatellites (simple sequence repeats, SSRs) still remain popular molecular markers for studying neutral genetic variation. Two alternative models outline how new microsatellite alleles evolve. Infinite alleles model (IAM) assumes that all possible alleles are equally likely to result from a mutation, while stepwise mutation model (SMM) describes microsatellite evolution as stepwise adding or subtracting single repeat units. Genetic relationships between individuals can be analyzed in higher precision when assuming the SMM scenario with allele size differences as a proxy of genetic distance. If population structure is not predetermined in advance, an empirical data analysis usually includes (a) estimating proximity between individual SSR profiles with a selected dissimilarity measure and (b) determining putative genetic structure of a given set of individuals using methods of clustering and/or ordination for the obtained dissimilarity matrix. We developed new dissimilarity indices between SSR profiles of haploid, diploid, or polyploid organisms assuming different mutation models and compared the performance of these indices for determining genetic structure with population data and with simulations. More specifically, we compared SMM with a constant or variable mutation rate at different SSR loci to IAM using data from natural populations of a freshwater bryozoan Cristatella mucedo (diploid), wheat leaf rust Puccinia triticina (dikaryon), and wheat powdery mildew Blumeria graminis (monokaryon). We show that inferences about population genetic structure are sensitive to the assumed mutation model. With simulations, we found that Bruvo's distance performs generally poorly, while the new metrics are capturing the differences in the genetic structure of the populations.     

Microsatellite analysis of genetic structure in the mangrove species Avicennia marina (Forsk.) Vierh. (Avicenniaceae)

The level of genetic variation throughout the entire worldwide range of the mangrove species Avicennia marina (Forsk.) Vierh. was examined using microsatellite markers. Three microsatellite loci detected high levels of allelic diversity (70 alleles in total), essential for an accurate estimation of population genetic parameters. The informativeness of the microsatellite loci tended to increase with increasing average number of repeats. The levels of heterozygosity detected for each population, over all loci, ranged from 0.0 to 0.8, with an average of 0.407, indicating that some populations had little or no genetic variation, whereas others had a large amount. Populations at the extremes of the distribution range showed reduced levels of heterozygosity, and significant levels of inbreeding. This is not unexpected as these populations may be subject to founder effects and environmental constraints. The presence of genetic structure was tested in A. marina populations using three models: (i) a single panmictic model; (ii) the discrete subpopulation model; and (iii) the isolation by distance model. The discrete subpopulations model was supported by the overall measures of population differentiation based on the infinite alleles model (F-statistics), and the stepwise mutation model (R statistics). In addition, an analysis of molecular variance (AMOVA), using both theoretical models, found that most of the variation was between populations (41-71%), and within individuals in the total population (31-49%). There was little variation among individuals within populations (0-10%). There was no significant isolation by distance. The high levels of genetic differentiation observed among populations of A. marina may be due to environmental and ecological factors, particularly past sea level and climatic changes.     

Joint estimation of migration rate and effective population size using the island model

Using the island model of population demography, I report that the demographic parameters migration rate and effective population size can be jointly estimated with equilibrium probabilities of identity in state calculated using a sample of genotypes collected at a single point in time from a single generation. The method, which uses moment-type estimators, applies to dioecious populations in which females and males have identical demography and monoecious populations with no selfing and requires that offspring genotypes are sampled following reproduction and prior to migration. I illustrate the estimation procedure using the infinite-island model with no mutation and the finite-island model with three kinds of mutation models. In the infinite-island model with no mutation, the estimators can be expressed as simple functions of estimates of the F-statistic parameters F(IT) and F(ST). In the finite-island model with mutation among k alleles, mutation rate, migration rate, and effective population size can be simultaneously estimated. The estimates of migration rate and effective population size are somewhat robust to violations in assumptions that may arise in empirical applications such as different kinds of mutation models and deviations from temporal equilibrium.     

Contrasting microsatellite variation between subalpine and western larch, two closely related species with different distribution patterns

Subalpine larch (Larix lyallii Parl.) and western larch (Larix occidentalis Nutt.) represent two closely related species with contrasting abundance and distribution patterns in Western North America. Genetic diversity at seven informative microsatellite loci was determined for 19 populations of subalpine larch and nine populations of western larch. Contrasting genetic diversity and patterns of population differentiation were observed between the two species. The overall within-population genetic diversity parameters were lower in subalpine larch (A = 3.2; A(P) = 3.6; H(E) = 0.418) than in western larch (A(P) = 5.51; H(E) = 0.580), a pattern that is likely related to historical or demographic factors. No evidence of interspecific hybridization was observed. Significantly more population differentiation (theta = 0.15; R(ST) = 0.07), consistent with more restricted gene flow, was observed for subalpine larch as compared to western larch (theta = 0.05; R(ST) = 0.04). Under the assumption of an infinite allele mutation model, 12 of the 19 subalpine larch populations showed signs of deviation from the mutation-drift equilibrium, which suggests Holocene population bottlenecks and fluctuations in effective population size for this species. None of the western larch populations deviated significantly from the mutation-drift equilibrium. For both species, Mantel's test revealed a significant positive relationship between geographical and genetic distances indicative of isolation by distance. A similar geographical structure was detected in both species, suggesting at least two genetically distinct glacial populations in each species. The various implications for gene conservation are discussed.     

Genetic population structure and contemporary dispersal patterns of a recent European invader, the Chinese mitten crab, Eriocheir sinensis

Genetic studies of recently established populations are challenging because the assumption of equilibrium underlying many analyses is likely to be violated. Using microsatellites, we investigated determinants of genetic structure and migration among invasive European-Chinese mitten crab populations, applying a combination of traditional population genetic analyses and nonequilibrium Bayesian methods. Consistent with their recent history, invasive populations showed much lower levels of genetic diversity than a native Chinese population, indicative of recent bottlenecks. Population differentiation was generally low but significant and especially pronounced among recently established populations. Significant differentiation among cohorts from the same geographical location (River Thames) suggests the low effective population size and associated strong genetic drift that would be anticipated from a very recent colonization. An isolation-by-distance pattern appears to be driven by an underlying correlation between geographical distance and population age, suggesting that cumulative homogenizing gene flow reduces founder bottleneck-associated genetic differentiation between longer-established populations. This hypothesis was supported by a coalescent analysis, which supported a drift + gene flow model as more likely than a model excluding gene flow. Furthermore, admixture analysis identified several recent migrants between the UK and Continental European population clusters. Admixture proportions were significantly predicted by the volume of shipping between sites, indicating that human-mediated transport remains a significant factor for dispersal of mitten crabs after the initial establishment of populations. Our study highlights the value of nonequilibrium methods for the study of invasive species, and also the importance of evaluating nonequilibrium explanations for isolation by distance patterns.     

Fine-scale spatial genetic structure and within population male-biased gene-flow in the grasshopper <Emphasis Type="Italic">Mioscirtus wagneri</Emphasis>

Dispersal is a life history trait that plays a key role in population dynamics, determining gene flow and influencing the size, structure and persistence of populations. For these reasons, the study of the genetic consequences of dispersal can be considered a central topic in both conservation and population genetics. In this study we examine the patterns of fine-scale genetic structure within two populations of the grasshopper Mioscirtus wagneri (Orhoptera: Acrididae). For this purpose, we have used seven species-specific microsatellite markers to type 266 individuals from two populations (Peña Hueca and El Salobral) located in Central Spain. We have found subtle genetic differentiation between some sampling patches and significant kinship structures up to 25 m distance which were particularly patent for females. In Peña Hueca locality, patterns of isolation-by-distance at both the patch scale and the individual level have also revealed an association between genetic differentiation/similarity and geographical distance in females but not in males. Overall, these data suggest a fine-scale spatial genetic substructure in the studied populations which seems to be mainly driven by female philopatry. Such pattern of within population genetic structure together with the inferred restricted dispersal distances is likely to contribute to reduce effective population sizes and inter-population gene flow. This can erode genetic variability and limit the colonization ability of this orthoptera, factors which can ultimately compromise the long-term persistence of their small size and isolated populations.     

Influence of mutational and sampling factors on the estimation of demographic parameters in a "continuous" population under isolation by distance

In numerous species, individual dispersal is restricted in space so that "continuous" populations evolve under isolation by distance. A method based on individual genotypes assuming a lattice population model was recently developed to estimate the product Dsigma2, where D is the population density and sigma2 is the average squared parent-offspring distance. We evaluated the influence on this method of both mutation rate and mutation model, with a particular reference to microsatellite markers, as well as that of the spatial scale of sampling. Moreover, we developed and tested a nonparametric bootstrap procedure allowing the construction of confidence intervals for the estimation of Dsigma2. These two objectives prompted us to develop a computer simulation algorithm based on the coalescent theory giving individual genotypes for a continuous population under isolation by distance. Our results show that the characteristics of mutational processes at microsatellite loci, namely the allele size homoplasy generated by stepwise mutations, constraints on allele size, and change of slippage rate with repeat number, have little influence on the estimation of Dsigma2. In contrast, a high genetic diversity (approximately 0.7-0.8), as is commonly observed for microsatellite markers, substantially increases the precision of the estimation. However, very high levels of genetic diversity (>0.85) were found to bias the estimation. We also show that statistics taking into account allele size differences give unreliable estimations (i.e., high variance of Dsigma2 estimation) even under a strict stepwise mutation model. Finally, although we show that this method is reasonably robust with respect to the sampling scale, sampling individuals at a local geographical scale gives more precise estimations of Dsigma2.     

THE INFLUENCE OF SELF-FERTILIZATION AND POPULATION DYNAMICS ON THE GENETIC STRUCTURE OF SUBDIVIDED POPULATIONS: A CASE STUDY USING MICROSATELLITE MARKERS IN THE FRESHWATER SNAIL BULINUS TRUNCATUS

The distribution of neutral genetic variability within and among sets of populations results from the combined actions of genetic drift, migration, extinction and recolonization processes, mutation, and the mating system. We here analyzed these factors in 38 populations of the hermaphroditic snail Bulinus truncatus. The sampling area covered a large part of the species range. The variability was analyzed using four polymorphic microsatellite loci. A very large number of alleles (up to 55) was found at the level of the whole study. Observed heterozygote deficiencies within populations are consistent with very high selfing rates, generally above 0.80, in all populations. These should depress the variability within populations, because of low effective size, genetic hitchhiking, and background selection, whatever the model of mutation assumed. However, that some populations exhibit much more variability than others suggests that historical demographic processes (e.g., population size variation, bottlenecks, or founding events) may play a significant role. A hierarchical analysis of the distribution of the variability across populations indicates a strong pattern of isolation by distance, whatever the geographical scale considered. Our analysis also illustrates how the mutation rate may affect population differentiation, as different mutation rates result in different levels of homoplasy at microsatellite loci. The effects of both genetic drift and gene flow vary with the temporal and spatial scales considered in B. truncatus populations.     

Low but significant genetic differentiation underlies biologically meaningful phenotypic divergence in a large Atlantic salmon population

Despite decades of research assessing the genetic structure of natural populations, the biological meaning of low yet significant genetic divergence often remains unclear due to a lack of associated phenotypic and ecological information. At the same time, structured populations with low genetic divergence and overlapping boundaries can potentially provide excellent models to study adaptation and reproductive isolation in cases where high‐resolution genetic markers and relevant phenotypic and life history information are available. Here, we combined single nucleotide polymorphism (SNP)‐based population inference with extensive phenotypic and life history data to identify potential biological mechanisms driving fine‐scale subpopulation differentiation in Atlantic salmon (Salmo salar) from the Teno River, a major salmon river in Europe. Two sympatrically occurring subpopulations had low but significant genetic differentiation (F_ST = 0.018) and displayed marked differences in the distribution of life history strategies, including variation in juvenile growth rate, age at maturity and size within age classes. Large, late‐maturing individuals were virtually absent from one of the two subpopulations, and there were significant differences in juvenile growth rates and size at age after oceanic migration between individuals in the respective subpopulations. Our findings suggest that different evolutionary processes affect each subpopulation and that hybridization and subsequent selection may maintain low genetic differentiation without hindering adaptive divergence.     

Geographical and genetic distances among zooplankton populations in a set of interconnected ponds: a plea for using GIS modelling of the effective geographical distance

In systems of interconnected ponds or lakes, the dispersal of zooplankton may be mediated by the active population component, with rivulets and overflows functioning as dispersal pathways. Using a landscape-based approach, we modelled the effective geographical distance among a set of interconnected ponds (De Maten, Genk, Belgium) in a Geographic Information System (GIS) environment. The first model (the Landscape Model; LM) corrects for the presence of direct connections among ponds and was based on the existing landscape structure (i.e. network of connecting elements among ponds, travelling distance and direction of the current). A second model (the Flow Rate Model; FRM) also incorporated field data on flow rates in the connecting elements as the driving force for the passive dispersal of the active zooplankton population component. Finally, the third model (the Dispersal Rate Model; DRM) incorporated field data on zooplankton dispersal rates. An analysis of the pattern of genetic differentiation among Daphnia ambigua populations inhabiting 10 ponds in the pond complex reveals that the effective geographical distance as modelled by the flow rate and the dispersal rate model provide a better approximation of the true rates of genetic exchange among populations than mere Euclidean geographical distances or the landscape model that takes into account solely the presence of physical connections.     

Arctic fox Vulpes lagopus population structure: circumpolar patterns and processes

Movement is a prominent process shaping genetic population structure. In many northern mammal species, population structure is formed by geographic distance, geographical barriers and various ecological factors that influence movement over the landscape. The Arctic fox Vulpes lagopus is a highly mobile, opportunistic carnivore of the Arctic that occurs in two main ecotypes with different ecological adaptations. We assembled microsatellite data in 7 loci for 1834 Arctic foxes sampled across their entire distribution to describe the circumpolar population structure and test the impact of (1) geographic distance, (2) geographical barriers and (3) ecotype designation on the population structure. Both Structure and Geneland demonstrated distinctiveness of Iceland and Scandinavia whereas low differentiation was observed between North America–northern Greenland, Svalbard and Siberia. Genetic differentiation was significantly correlated to presence of sea ice on a global scale, but not to geographical distance or ecotype designation. However, among areas connected by sea ice, we recorded a pattern of isolation by distance. The maximum likelihood approach in Migrate suggested that connectivity across North America–northern Greenland and Svalbard was particularly high. Our results demonstrate the importance of sea ice for maintaining connectivity between Arctic fox populations and we therefore predict that climate change will increase genetic divergence among populations in the future.     

Genetic structure of Araucaria angustifolia (Araucariaceae) populations in Brazil: implications for the in situ conservation of genetic resources

The distribution of the genetic variation within and among natural populations of A. ANGUSTIFOLIA growing in different regions in Brazil was assessed at microsatellite and AFLP markers. Both markers revealed high gene diversity ( H = 0.65; AR = 9.1 for microsatellites and H = 0.27; P = 77.8 % for AFLPs), moderate overall differentiation ( RST = 0.13 for microsatellites and FST = 0.10 for AFLPs), but high divergence of the northernmost, geographically isolated population. In a Bayesian analysis, microsatellite data suggested population structure at two levels: at K = 2 and at K = 3 in agreement to the geographical distribution of populations. This result was confirmed by the UPGMA dendrogram based on microsatellite data (bootstrap support > 95 %). Non-hierarchical AMOVA revealed high variation among populations from different A POSTERIORI defined geographical groups. The genetic distance between sample locations increased with geographical distance for microsatellites ( R = 0.62; P = 0.003) and AFLPs ( R = 0.32; P = 0.09). This pattern of population differentiation may be correlated with population history such as geographical isolation and postglacial colonization of highlands. Implications of the population genetic structure for the conservation of genetic resources are discussed.