Gene flow and melanism in garter snakes revisited: a comparison of molecular markers and island vs. coalescent models

Within populations, the stochastic effect of genetic drift and deterministic effect of natural selection are potentially weakened or altered by gene flow among populations. The influence of gene flow on Lake Erie populations of the common garter snake has been of particular interest because of a discontinuous colour pattern polymorphism (striped vs. melanistic) that is a target of natural selection. We reassessed the relative contributions of gene flow and genetic drift using genetic data and population size estimates. We compared all combinations of two marker systems and two analytical approaches to the estimation of gene flow rates: allozymes (data previously published), microsatellite DNA (new data), the island model ( F _ST-based approach), and a coalescence-based approach. For the coalescence approach, mutation rates and sampling effects were also investigated. While the two markers produced similar results, gene flow based on F _ST was considerably higher (Nm > 4) than that from the coalescence-based method (Nm < 1). Estimates of gene flow are likely to be inflated by lack of migration-drift equilibrium and changing population size. Potentially low rates of gene flow (Nm < 1), small population size at some sites, and positive correlations of number of microsatellite DNA alleles and island size and between M , mean ratio of number of alleles to range in allele size, and island size suggest that in addition to selection, random genetic drift may influence colour pattern frequencies. © 2003 The Linnean Society of London, Biological Journal of the Linnean Society , 2003, 79, 389–399.     

Genetic structuring in fluctuating populations

The effect of genetic drift in spatially distributed dispersal-linked and density-regulated populations is studied in a classical one-locus two-allele system. We analyse emergence of genetic differentiation assuming random drift only, where the noise-like variability is due to demographic stochasticity. We find emergence of clusters of sub-units with local allele fixation and persistence of both alleles in lengthy simulations. We demonstrate that local allele fixation (extending over a number of adjoining spatial sub-units) – without global loss of alleles – may occur when the carrying capacities of local patches are small, under a full range population dynamic regimes, when dispersal rate is small, and when redistribution (through dispersal) does not act as global mixer. These results are novel. The key to the observations is that drift is simultaneously influenced by distance-dependent dispersal, demographic stochasticity and autocorrelated population fluctuations due to delayed-density dependence. These are standard elements of contemporary population models in spatially structured context. With stable large populations, no stochasticity and dispersal limited to neighbours only, our model collapses to the stepping-stone model, while with dispersal being random and global, the model collapses to Wright's island model.     

Microsatellites and the Genetics of Highly Selfing Populations in the Freshwater Snail Bulinus Truncatus

Hermaphrodite tropical freshwater snails provide a good opportunity to study the effects of mating system and genetic drift on population genetic structure because they are self-fertile and they occupy transient patchily distributed habitats (ponds). Up to now the lack of detectable allozyme polymorphism prevented any intrapopulation studies. In this paper, we examine the consequences of selfing and bottlenecks on genetic polymorphism using microsatellite markers in 14 natural populations (under a hierarchical sampling design) of the hermaphrodite freshwater snail Bulinus truncatus. These population genetics data allowed us to discuss the currently available mutation models for microsatellite sequences. Microsatellite markers revealed an unexpectedly high levels of genetic variation with <=41 alleles for one locus and gene diversity of 0.20-0.75 among populations. The values of any estimator of F(is) indicate high selfing rates in all populations. Linkage disequilibria observed at all loci for some populations may also indicate high levels of inbreeding. The large extent of genetic differentiation measured by F(st), R(st) or by a test for homogeneity between genic distributions is explained by both selfing and bottlenecks. Despite a limited gene flow, migration events could be detected when comparing different populations within ponds.     

Imprints from genetic drift and mutation imply relative divergence times across marine transition zones in a pan-European small pelagic fish (Sprattus sprattus)

Geographic distributions of most temperate marine fishes are affected by postglacial recolonisation events, which have left complex genetic imprints on populations of marine species. This study investigated population structure and demographic history of European sprat (Sprattus sprattus L.) by combining inference from both mtDNA and microsatellite genetic markers throughout the species' distribution. We compared effects from genetic drift and mutation for both genetic markers in shaping genetic differentiation across four transition zones. Microsatellite markers revealed significant isolation by distance and a complex population structure across the species' distribution (overall θ(ST)=0.038, P<0.01). Across transition zones markers indicated larger effects of genetic drift over mutations in the northern distribution of sprat contrasting a stronger relative impact of mutation in the species' southern distribution in the Mediterranean region. These results were interpreted to reflect more recent divergence times between northern populations in accordance with previous findings. This study demonstrates the usefulness of comparing inference from different markers and estimators of divergence for phylogeographic and population genetic studies in species with weak genetic structure, as is the case in many marine species.     

Experimental estimation of mutation rates in a wheat population with a gene genealogy approach

Microsatellite markers are extensively used to evaluate genetic diversity in natural or experimental evolving populations. Their high degree of polymorphism reflects their high mutation rates. Estimates of the mutation rates are therefore necessary when characterizing diversity in populations. As a complement to the classical experimental designs, we propose to use experimental populations, where the initial state is entirely known and some intermediate states have been thoroughly surveyed, thus providing a short timescale estimation together with a large number of cumulated meioses. In this article, we derived four original gene genealogy-based methods to assess mutation rates with limited bias due to relevant model assumptions incorporating the initial state, the number of new alleles, and the genetic effective population size. We studied the evolution of genetic diversity at 21 microsatellite markers, after 15 generations in an experimental wheat population. Compared to the parents, 23 new alleles were found in generation 15 at 9 of the 21 loci studied. We provide evidence that they arose by mutation. Corresponding estimates of the mutation rates ranged from 0 to 4.97 x 10(-3) per generation (i.e., year). Sequences of several alleles revealed that length polymorphism was only due to variation in the core of the microsatellite. Among different microsatellite characteristics, both the motif repeat number and an independent estimation of the Nei diversity were correlated with the novel diversity. Despite a reduced genetic effective size, global diversity at microsatellite markers increased in this population, suggesting that microsatellite diversity should be used with caution as an indicator in biodiversity conservation issues.     

Microsatellite genetic variation in small and isolated populations of Magnolia sieboldii ssp. japonica

Magnolia sieboldii ssp. japonica, distributed mainly in western Japan, is restricted to high elevation areas (1000-2000 m above sea level) and usually forms small isolated populations. Four microsatellite loci were assayed for 19 populations from six regions spanning the range of distribution, and the levels and distribution of genetic variation were estimated. All four loci were variable, with a total of 39 alleles, but the overall level of microsatellite genetic variation was low, especially compared with a related species, M. obovata. Genetic structure in M. sieboldii was characterised by low intrapopulational genetic variation (A = 3.74 and H(o) = 0.366 on average) and high genetic differentiation even among regional populations. Highly significant isolation-by-distance (IBD) models at the short distance were detected. Genetic drift and limited gene flow was considered to be important in determining the genetic structure within regions. Total genetic differentiation was remarkably high (F(ST) = 0.488 and R(ST) = 0.538), suggesting genetic barriers among regions. Neighbour-joining dendrograms relating the 19 populations, and further analysis on the IBD models, revealed that a stepwise mutation model was more suited than an infinite allele model to explain the genetic differentiation among regions. It is suggested that mutation at microsatellite loci might be influential in generating the genetic differentiation among regions. These results showed the potential of hypervariable microsatellite loci to evaluate the effects of genetic drift and population isolation within regions, and to detect genetic distinctiveness, in spite of the loss of overall genetic variation in M. sieboldii.     

Spatial and temporal patterns of neutral and adaptive genetic variation in the endangered African wild dog (Lycaon pictus)

Deciphering patterns of genetic variation within a species is essential for understanding population structure, local adaptation and differences in diversity between populations. Whilst neutrally evolving genetic markers can be used to elucidate demographic processes and genetic structure, they are not subject to selection and therefore are not informative about patterns of adaptive variation. As such, assessments of pertinent adaptive loci, such as the immunity genes of the major histocompatibility complex (MHC), are increasingly being incorporated into genetic studies. In this study, we combined neutral (microsatellite, mtDNA) and adaptive (MHC class II DLA‐DRB1 locus) markers to elucidate the factors influencing patterns of genetic variation in the African wild dog (Lycaon pictus); an endangered canid that has suffered extensive declines in distribution and abundance. Our genetic analyses found all extant wild dog populations to be relatively small (N_e < 30). Furthermore, through coalescent modelling, we detected a genetic signature of a recent and substantial demographic decline, which correlates with human expansion, but contrasts with findings in some other African mammals. We found strong structuring of wild dog populations, indicating the negative influence of extensive habitat fragmentation and loss of gene flow between habitat patches. Across populations, we found that the spatial and temporal structure of microsatellite diversity and MHC diversity were correlated and strongly influenced by demographic stability and population size, indicating the effects of genetic drift in these small populations. Despite this correlation, we detected signatures of selection at the MHC, implying that selection has not been completely overwhelmed by genetic drift.     

Genetic Polymorphism and Bottleneck Analysis of Balkhi,Hashtnagri, and Michni Sheep Populations Using Microsatellite Markers

Pakistan is rich in sheep genetic resources. Balkhi, Hashtnagri, and Michni are neighboring sheep populations found in Khyber Pakhtunkhwa province of Pakistan. In this study, we analyzed the genetic structures and bottleneck incidents within these sheep populations using 31 microsatellite DNA markers. Total numbers of 116, 100, and 95 alleles, with average numbers of 3.20, 3.26, and 3.74 alleles per locus were observed, respectively, in Balkhi, Hashtnagri, and Michni population. Mean observed heterozygosity was 0.402 in Balkhi, 0.416 in Hashtnagri, and 0.522 in Michni population. All the three sheep populations showed significantly high inbreeding. Michni population was found to be in mutation drift equilibrium, showing the absence of genetic bottleneck. The data of Balkhi and Hashtnagri indicated the presence of genetic bottleneck in these populations. These results suggest a moderate level of genetic diversity within Michni population that may be useful for breed improvement programs. Hashtnagri and Balkhi populations having low within breed genetic variability may contain some valuable characteristics that need to be conserved.     

Molecular and genetic characterization of the allelic variants of Du215, Du281, Du323, and Du47G microsatellite loci in parthenogenetic lizard <Emphasis Type="Italic">Darevskia armeniaca</Emphasis> (Lacertidae)

A key issue in the study of unisexual (parthenogenetic) vertebrate species is the determination of their genetic and clonal diversity. In pursuing this aim, various markers of nuclear and mitochondrial genomes can be used. The most effective genetic markers include microsatellite DNA, characterized by high variability. The development and characterization of such markers is a necessary step in the genetic studies of parthenogenetic species. In the present study, using locus-specific PCR, for the first time, an analysis of allelic polymorphism of four microsatellite loci is performed in the populations of parthenogenetic species Darevskia armeniaca. In the studied populations, allelic variants of each locus are identified, and the nucleotide sequences of each allele are determined. It is demonstrated that allele differences are associated with the variation in the structure of microsatellite clusters and single nucleotide substitutions at fixed distances in flanking DNA regions. Structural allele variations form haplotype markers that are specific to each allele and are inherited from their parental bisexual species. It is established which of the parental alleles of each locus were inherited by the parthenogenetic species. The characteristics of the distribution and frequency of the alleles of microsatellite loci in the populations of D. armeniaca determining specific features of each population are obtained. The observed heterozygosity of the populations at the studied loci and the mutation rates in genome regions, as well as Nei’s genetic distances between the studied populations, are determined, and the phylogenetic relationships between them are established.     

Long-term population genetic structure of an invasive urochordate: the ascidian Botryllus schlosseri

The accelerated pace of marine biological invasions raises questions pertaining to genetic traits and dynamics underlying the successful establishment of invasive species. Current research stresses the importance of multiple introductions and prolonged gene flow as the main sources for genetic diversity, which, along with genetic drift, affect invasive species success. We here attempt to determine the relative contribution of gene flow and mutation rates as sources of genetic variability using the invasive tunicate Botryllus schlosseri as a model. The study was performed over a 13-year period in the Santa Cruz Harbor, California. With a characteristic life history of five generations/year, the Santa Cruz Botryllus population has already experienced approximately 155 generations since the onset of its invasion. The results (278 specimens, 127 scored alleles, five microsatellite loci) support limited gene flow rate (2.89?×?10^?3) and relative genetic isolation. Furthermore, the study population was found to be influenced by both, genetic drift and a high mutation rate (2.47?×?10^?2). These findings were supported by high fluctuations in the frequencies of microsatellite alleles, the appearance of new alleles and the loss of others. The balance between genetic drift and a high mutation rate is further elucidated by the high, stable level of genetic variation. We suggest that rapid mutation rates at the microsatellite loci reflect genome-wide phenomena, helping to maintain high genetic variability in relatively isolated populations. The potential adaptability to new environments is discussed.     

Demographic history and genetic diversity in West Indian Coereba flaveola populations

The bananaquit (Coereba flaveola) has been well studied throughout the Caribbean region from a phylogenetic perspective. However, data concerning the population genetics and long-term demography of this bird species are lacking. In this study, we focused on three populations within the Lesser Antilles and one on Puerto Rico and assessed genetic and demographic processes, using five nuclear and two mitochondrial markers. We found that genetic diversity of bananaquits on Puerto Rico exceeds that on the smaller islands (Dominica, Guadeloupe and Grenada); this might reflect either successive founder events from Puerto Rico to Grenada, or more rapid drift in smaller populations subsequent to colonization. Population growth rate estimates showed no evidence of rapid expansion and migration was indicated only between populations from the closest islands of Dominica and Guadeloupe. Overall, our results suggest that a "demographic fission" model, considering only mutation and drift, but without migration, can be applied to these bananaquit populations in the West Indies.     

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.     

Comparative Genetic Diversity of Wild and Captive Populations of the Bare-Faced Curassow (Crax fasciolata) Based on Cross-Species Microsatellite Markers: Implications for Conservation and Management

The bare-faced curassow (Crax fasciolata) is a large Neotropical bird that suffers anthropogenic pressure across much of its range. A captive population is maintained for conservation management, although there has been no genetic screening of stocks. Based on the six microsatellite markers developed for Crax globulosa, the genetic variability of C. fasciolata and possible differences between a wild and a captive population were investigated. Only three loci were polymorphic, with a total of 27 alleles. More than half of these alleles were private to the wild (n = 8) or captive (n = 7) populations. Significant deviations from Hardy–Weinberg equilibrium were restricted to the captive population. Despite the number of private alleles, genetic drift has probably promoted differentiation between populations. Our results indicate that wild C. fasciolata populations are genetically impoverished and structured, but species-specific microsatellite markers will be necessary for a more reliable assessment of the species’ genetic diversity.     

Genetic divergence in the small Indian mongoose (Herpestes auropunctatus), a widely distributed invasive species

The combination of founder events, random drift and new selective forces experienced by introduced species typically lowers genetic variation and induces differentiation from the ancestral population. Here, we investigate microsatellite differentiation between introduced and native populations of the small Indian mongoose (Herpestes auropunctatus). Many expectations based on introduction history, such as loss of alleles and relationships among populations, are confirmed. Nevertheless, when applying population assignment methods to our data, we observe a few specimens that are incorrectly assigned and/or appear to have a mixed ancestry, despite estimates of substantial population differentiation. Thus, we suggest that population assignments of individuals should be viewed as tentative and that there should be agreement among different algorithms before assignments are applied in conservation or management. Further, we find no congruence between previously reported morphological differentiation and the sorting of microsatellite variation. Some introduced populations have retained much genetic variation while others have not, irrespective of morphology. Finally, we find alleles from the sympatric grey mongoose (Herpestes edwardsii) in one small Indian mongoose within the native range, suggesting an alternative explanation for morphological differentiation involving a shift in female preferences in allopatry.     

Contrasting patterns of mitochondrial and microsatellite population structure in fragmented populations of greater prairie-chickens

Greater prairie-chickens (Tympanuchus cupido pinnatus) were once found throughout the tallgrass prairie of midwestern North America but over the last century these prairies have been lost or fragmented by human land use. As a consequence, many current populations of prairie-chickens have become isolated and small. This fragmentation of populations is expected to lead to reductions in genetic variation as a result of random genetic drift and a decrease in gene flow. As expected, we found that genetic variation at both microsatellite DNA and mitochondrial DNA (mtDNA) markers was reduced in smaller populations, particularly in Wisconsin. There was relatively little range-wide geographical structure (FST) when we examined mtDNA haplotypes but there was a significant positive relationship between genetic (FST) and geographical distance (isolation by distance). In contrast, microsatellite DNA loci revealed significant geographical structure (FST) and a weak effect of isolation by distance throughout the range. These patterns were much stronger when populations with reduced levels of genetic variability (Wisconsin) were removed from the analyses. This suggests that the effects of genetic drift were stronger than gene flow at microsatellite loci, whereas these forces were in range-wide equilibrium at mtDNA markers. These differences between the two molecular markers may be explained by a larger effective population size (Ne) for mtDNA, which is expected in species such as prairie-chickens that have female-biased dispersal and high levels of polygyny. Our results suggest that historic populations of prairie-chickens were once interconnected by gene flow but current populations are now isolated. Thus, maintaining gene flow may be important for the long-term persistence of prairie-chicken populations.     

The impacts of inbreeding,drift and selection on genetic diversity in captive breeding populations

The goal of captive breeding programmes is often to maintain genetic diversity until re‐introductions can occur. However, due in part to changes that occur in captive populations, approximately one‐third of re‐introductions fail. We evaluated genetic changes in captive populations using microsatellites and mtDNA. We analysed six populations of white‐footed mice that were propagated for 20 generations using two replicates of three protocols: random mating (RAN), minimizing mean kinship (MK) and selection for docility (DOC). We found that MK resulted in the slowest loss of microsatellite genetic diversity compared to RAN and DOC. However, the loss of mtDNA haplotypes was not consistent among replicate lines. We compared our empirical data to simulated data and found no evidence of selection. Our results suggest that although the effects of drift may not be fully mitigated, MK reduces the loss of alleles due to inbreeding more effectively than random mating or docility selection. Therefore, MK should be preferred for captive breeding. Furthermore, our simulations show that incorporating microsatellite data into the MK framework reduced the magnitude of drift, which may have applications in long‐term or extremely genetically depauperate captive populations.     

Microsatellite DNA variability in the populations of muskoxen <Emphasis Type="Italic">Ovibos moschatus</Emphasis> transplanted into the Russian North

The muskoxen populations introduced to the Taimyr Peninsula and Wrangel Island in 1974 to 1975 were examined for sequence variation at seven microsatellite loci. Donor material originated from the populations of Banks Island (Canada) and Eastern Greenland. Relative to the allele frequencies, both introduced populations demonstrated rather strong deviation from the populations of the native range. At the same time, population allelic structures evidenced that they were closer to the Greenland populations. Estimates of genetic diversity at microsatellite loci (expected heterozygosity and the allele number) in the introduced muskoxen were found to be high for populations originating from a small number of founder individuals. In the immigrants, linkage disequilibrium and deviation of the genotype frequencies from the Hardy-Weinberg proportions were observed, which was mainly caused by the deficit of heterozygotes. The same pattern was also typical of native populations and was explained in terms of specific population structure and demographic processes. The latter were manifested as a periodic decline of the effective population size, resulting in the prevailing influence of genetic drift and inbreeding. The consequences of genetic drift were not as dramatic, as could be expected, which may be explained by a high mutation rate of neutral microsatellite loci and fast growth of the new populations.     

Recent demographic bottlenecks are not accompanied by a genetic signature in banner-tailed kangaroo rats (Dipodomys spectabilis)

Single-sample methods of bottleneck detection are now routine analyses in studies of wild populations and conservation genetics. Three common approaches to bottleneck detection are the heterozygosity excess, mode-shift, and M-ratio tests. Empirical groundtruthing of these methods is difficult, but their performances are critical for the accurate reconstruction of population demography. We use two banner-tailed kangaroo rat (Dipodomys spectabilis) populations from southeastern Arizona (USA) that are known to have experienced recent demographic reductions to search for genetic bottleneck signals with eight microsatellite loci. Over eight total sample-years, neither population showed a genetic bottleneck signature. M-ratios in both populations were large, stable, and never fell below a critical significance value (Mc). The mode shift test did not detect any distortion of allele frequencies, and tests of heterozygosity excess were not significant in postbottleneck samples when we used standard microsatellite mutation models. The genetic effects of bottlenecks like those experienced by our study populations should be strongly influenced by rates of mutation and migration. We used genetic parentage data to estimate a relatively high mutation rate in D. spectabilis (0.0081 mutants/generation/locus), but mutation alone is unlikely to explain the temporal distribution of rare alleles that we observed. Migration (gene flow) is a more likely explanation, despite prior mark-recapture analysis that estimated very low rates of interpopulation dispersal. We interpret our kangaroo rat data in light of the broader literature and conclude that in natural populations connected by dispersal, demographic bottlenecks may prove difficult to detect using molecular genetic data.     

Low genetic variation support bottlenecks in Scandinavian red deer

Loss of genetic variation from genetic drift during population bottlenecks has been shown for many species. Red deer (Cervus elaphus) may have been exposed to bottlenecks due to founder events during postglacial colonisation in the early Holocene and during known population reductions in the eighteenth and nineteenth centuries. In this study, we assess loss of genetic variation in Scandinavian red deer due to potential bottlenecks by comparing microsatellite (n = 14) and mitochondrial DNA variation in the Norwegian and Swedish populations with the Scottish, Lithuanian and Hungarian populations. Bottlenecks are also assessed from the M ratio of populations, heterozygosity excess and from hierarchical Bayesian analyses of their demographic history. Strong genetic drift and differentiation was identified in both Scandinavian populations. Microsatellite variation was lower in both Scandinavian populations compared with the other European populations and mitochondrial DNA variation was especially low in the Swedish population where only one unique haplotype was observed. Loss of microsatellite alleles was demonstrated by low M ratios in all populations except the Hungarian. M ratios’ were especially low in the Scandinavian populations, indicating additional or more severe bottlenecks. Heterozygosity excess compared with the expectation from the number of observed microsatellite alleles suggested a recent bottleneck of low severity in the Norwegian population. Hierarchical Bayesian coalescent analyses consistently yielded estimates of a large ancestral and a small current population size in all investigated European populations and suggested the onset of population decline to be between 5,000 and 10,000 years ago, which coincide well with postglacial colonisation.     

The population genetic effects of ancestry and admixture in a subdivided cattle breed

The genetic structure of the Dexter, a minority cattle breed with complex demographic history, was investigated using microsatellite markers and a range of statistical approaches designed to detect both admixture and genetic drift. Modern representatives of two putative ancestral populations, the Devon and Kerry, together with the different populations of the Dexter, which have experienced different demographic histories, were analysed. Breed units showed comparatively high levels of genetic variability ( H _E = 0.63–0.68); however, distinct genetic subgroups were detected within the Dexter, which could be attributed to known demographic events. Much lower diversity was identified in three small, isolated Dexter populations ( H _E = 0.52–0.55) and higher differentiation ( F _ST > 0.13) was found. For one of these populations, where strong selection has taken place, we also found evidence of a demographic bottleneck. Three methods for quantifying breed admixture were applied and substantial method-based variation in estimates for the genetic contribution of the two proposed ancestral populations for each subdivision of the Dexter was found. Results were consistent only in the case of a group consisting of selected Traditional Dexter animals, where the ancestor of the modern Kerry breed was also determined as the greater parental contributor to the Dexter. The inconsistency of estimation of admixture proportions between the methods highlights the potentially confounding role of genetic drift in shaping small population structure, and the consequences of accurately describing population histories from contemporary genetic data.