Widespread gene flow and high genetic variability in populations of water voles Arvicola terrestris in patchy habitats

Theory predicts that the impact of gene flow on the genetic structure of populations in patchy habitats depends on its scale and the demographic attributes of demes (e.g. local colony sizes and timing of reproduction), but empirical evidence is scarce. We inferred the impact of gene flow on genetic structure among populations of water voles Arvicola terrestris that differed in average colony sizes, population turnover and degree of patchiness. Colonies typically consisted of few reproducing adults and several juveniles. Twelve polymorphic microsatellite DNA loci were examined. Levels of individual genetic variability in all areas were high ( H _O= 0.69–0.78). Assignments of juveniles to parents revealed frequent dispersal over long distances. The populations showed negative F _IS values among juveniles, F _IS values around zero among adults, high F _ST values among colonies for juveniles, and moderate, often insignificant, F _ST values for parents. We inferred that excess heterozygosity within colonies reflected the few individuals dispersing from a large area to form discrete breeding colonies. Thus pre-breeding dispersal followed by rapid reproduction results in a seasonal increase in differentiation due to local family groups. Genetic variation was as high in low-density populations in patchy habitats as in populations in continuous habitats used for comparison. In contrast to most theoretical predictions, we found that populations living in patchy habitats can maintain high levels of genetic variability when only a few adults contribute to breeding in each colony, when the variance of reproductive success among colonies is likely to be low, and when dispersal between colonies exceeds nearest-neighbour distances.     

Population genetic structure of the butterfly Melitaea didyma (Nymphalidae) along a northern distribution range border

The population genetic structure of the butterfly Melitaea didyma was studied along the northern distribution range border in Central Germany by means of allozyme electrophoresis. Individuals were sampled from a total of 21 habitat patches from four regions, and two provinces. Sampling was designed to estimate local vs. regional differentiation. High levels of variability were found, H _e= 0.14–0.21. The mean expected sample heterozygosity from one region, Mosel, was significantly lower than from the Hammelburg region, H _e= 0.17 and 0.19, respectively. Two hierarchical levels of genetic differentiation were found. Within regions individuals sampled from different patches behaved as belonging to one population with high levels of gene flow (Hammelburg F _ST= 0.015, Mosel F _ST= 0.044), though local isolation barriers did create a substructuring of these populations. The inbreeding coefficients, F _IS, were constant over all sample levels, suggesting a similar distribution of habitat patches within regions. Between regions gene flow was limited. An isolation by distance analysis indicated that the hierarchical structure, at the provincial level, may be breaking down due to isolation of regional populations. A more general observation was that the sampling design may greatly have influenced the estimation of genetic differentiation. Depending on which samples were included, overall F _ST estimates ranged from 0.059–0.090.     

Microsatellite analysis of genetic population structure in an endangered salmonid: the coastal cutthroat trout (Oncorhynchus clarki clarki)

The genetic population structure of coastal cutthroat trout ( Oncorhynchus clarki clarki ) in Washington state was investigated by analysis of variation in allele frequencies at six highly polymorphic microsatellite loci for 13 anadromous populations, along with one outgroup population from the Yellowstone subspecies ( O. clarki bouvieri) (mean heterozygosity = 67%; average number of alleles per locus = 24). Tests for genetic differentiation revealed highly significant differences in genotypic frequencies for pairwise comparisons between all populations within geographical regions and overall population subdivision was substantial ( F _ST = 0.121, R _ST = 0.093), with 44.6% and 55.4% of the among-population diversity being attributable to differences between streams ( F _SR = 0.054) and between regions ( F _RT = 0.067), respectively. Analysis of genetic distances and geographical distances did not support a simple model of isolation by distance for these populations. With one exception, neighbour-joining dendrograms from the Cavalli-Sforza and Edwards' chord distances and maximum likelihood algorithms clustered populations by physiogeographic region, although overall bootstrap support was relatively low (53%). Our results suggest that coastal cutthroat trout populations are ultimately structured genetically at the level of individual streams. It appears that the dynamic balance between gene flow and genetic drift in the subspecies favours a high degree of genetic differentiation and population subdivision with the simultaneous maintenance of high heterozygosity levels within local populations. Results are discussed in terms of coastal cutthroat trout ecology along with implications for the designation of evolutionarily significant units pursuant to the US Endangered Species Act of 1973 and analogous conservation units.     

Stability of genetic structure and effective population size inferred from temporal changes of microsatellite DNA polymorphisms in the land snail Helix aspersa (Gastropoda: Helicidae)

Temporal evolution of genetic variability may have far-reaching consequences for a diverse array of evolutionary processes. Within the polders of the Bay of Mont-Saint-Michel (France), populations of the land snail Helix aspersa are characterized by a metapopulation structure with occasional extinction processes resulting from farming practices. A temporal survey of genetic structure in H . aspersa was carried out using variability at four microsatellite loci, in ten populations sampled two years apart. Levels of within-population genetic variation, as measured by allelic richness, H _e or F _is , did not change over time and similar levels of population differentiation were demonstrated for both sampling years. The extent of genetic differentiation between temporal samples of the same population established (i) a stable structure for six populations, and (ii) substantial genetic changes for four populations. Using classical F -statistics and a maximum likelihood method, estimates of the effective population size ( N _e) illustrated a mixture of stable populations with high N _e, and unstable populations characterized by very small N _e estimates (of 5–11 individuals). Owing to human disturbances, intermittent gene flow and genetic drift are likely to be the predominant evolutionary processes shaping the observed genetic structure. However, the practice of multiple matings and sperm storage is likely to provide a reservoir of variability, minimizing the eroding genetic effects of population size reduction and increasing the effective population size.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society , 2004, 82 , 89–102.     

Comparison of quantitative and molecular genetic variation of native vs. invasive populations of purple loosestrife (Lythrum salicaria L., Lythraceae)

Study of adaptive evolutionary changes in populations of invasive species can be advanced through the joint application of quantitative and population genetic methods. Using purple loosestrife as a model system, we investigated the relative roles of natural selection, genetic drift and gene flow in the invasive process by contrasting phenotypical and neutral genetic differentiation among native European and invasive North American populations ( Q _ST −  F _ST analysis). Our results indicate that invasive and native populations harbour comparable levels of amplified fragment length polymorphism variation, a pattern consistent with multiple independent introductions from a diverse European gene pool. However, it was observed that the genetic variation reduced during subsequent invasion, perhaps by founder effects and genetic drift. Comparison of genetically based quantitative trait differentiation ( Q _ST) with its expectation under neutrality ( F _ST) revealed no evidence of disruptive selection ( Q _ST >  F _ST) or stabilizing selection ( Q _ST <  F _ST). One exception was found for only one trait (the number of stems) showing significant sign of stabilizing selection across all populations. This suggests that there are difficulties in distinguishing the effects of nonadaptive population processes and natural selection. Multiple introductions of purple loosestrife may have created a genetic mixture from diverse source populations and increased population genetic diversity, but its link to the adaptive differentiation of invasive North American populations needs further research.     

Allozyme Diversity in Populations of Cymbidium goeringii (Orchidaceae)

Abstract: Using 14 allozyme loci, we investigated levels of genetic diversity within populations, and degree of genetic divergence among 24 populations of Cymbidium goeringii (Orchidaceae) in Korea and Japan. Cymbidium goeringii maintains high levels of genetic diversity both at population (mean expected heterozygosity, H _e = 0.238) and species levels (0.260). Means of H _e found in 24 populations were not significantly different from each other. About 90 % of the total variation in the species is common to all populations (mean G _ST = 0.108). No unique allele was found in any population. The indirect estimate of gene flow based on the mean G _ST was high ( Nm = 2.06). Nei's genetic identities for pairs of populations had high values (mean = 0.974 [SD = 0.013]). The Mantel-Z test showed a significant correlation between genetic distance and geographic distance. However, the mean G _ST value between 17 populations in Korea and seven Japanese populations was relatively low (0.029), even though the land connection between the southern Korean peninsula and southern Japanese archipelagos has not existed since the middle Pleistocene. Large numbers of small seeds of C. goeringii might travel long distances by wind from populations to populations both in Korea and Japan, increasing genetic diversity within populations and maintaining low genetic differentiation among populations.     

Genetic consequences of an invasion through a patchy environment — the cynipid gallwasp Andricus quercuscalicis (Hymenoptera: Cynipidae)

Over the last 300–400 years, the cynipid gallwasp Andricus quercuscalicis has invaded northern and western Europe following human introduction of an obligate host plant, the Turkey oak ( Quercus cerris ) from south-eastern Europe. In the introduced range, distances between Turkey oak patches are greater and the numbers of oaks in each patch are far lower than in its native range. These changes in spatial distribution of Turkey oak are predicted to result in high genetic subdivision of A. quercuscalicis in invaded areas relative to its native range. Allozyme electrophoresis was used to examine genetic variation in 823 gall wasps from 39 populations of A. quercuscalicis. No new electrophoretic variants were found in the invaded range and both allelic diversity and mean heterozygosity decreased significantly with distance from the native range. Spatial autocorrelation analysis and values of Wright's F_st indicate that differences in allele frequences between populations were substantially greater in the invaded range than in the native range. Spatial autocorrelation analysis also suggests that changes in allele frequencies across Europe are unlikely to be the results of selection, but rather of strong directional migration followed by limited gene flow between populations. Patterns of genetic differentiation across Europe suggest that populations of A. quercuscalicis have been founded sequentially from the east through a process of random genetic subsampling and not by source-sink colonization directly from the native range.     

Genetic diversity and genetic structure of populations of Ranunculus japonicus Thunb. (Ranunculaceae)

Abstract In order to clarify the genetic diversity and population structure of Ranunculus japonicus , allozymic analysis was conducted on 60 populations in southwestern Japan. Considerable genetic variati ons were detected among the populations of R. japonicus . The genetic diversities within species ( H _es = 0.215) and within populations ( H _ep = 0.172) were slightly higher than those of other perennial herbs with widespread distribution and outcrossing plants. Significantly higher values of fixation index were detected in some populations, which might have arisen from restricted mating partners. The majority of genetic variation (approx. 80%) resided within a population and a moderate level of genetic differentiation ( G _ST = 0.203) was observed among populations. The F _ST value (0.203) suggests the existence of a substantial population structure in this species. The highly significant correlation between geographic distance and F _ST values indicates that isolation by distance has played an important role in the construction of the genetic structure of this species.     

Spatial and temporal population structure of sea trout at the Island of Gotland, Sweden, delineated from mitochondrial DNA

In a study of the genetic relationships among 879 anadromous brown trout Salmo trutta from 13 streams at the Island of Gotland, Sweden, using RFLP analysis of a mitochondrial DNA segment (NADH dehydrogenase-1 gene), six haplotypes were detected. Significant genetic divergence was observed among streams as well as between cohorts within streams. Approximately 8–9% of the total variation was due to differences between populations, and 4–5% was explained by differences between cohorts within populations. The female effective population size ( N _ef) was assessed from temporal haplotype frequency differences between consecutive cohorts; the estimated average N _ef over all populations was just below 30, suggesting that these populations were effectively quite small. With such small effective sizes the populations are expected to loose genetic variability quickly, but the observed levels do not appear particularly low. This indicates that female migration between streams occurs. The observed level of differentiation does not support the presumption that a particular pre-smolt migratory behaviour observed in Gotland streams, with large portions of fry leaving for the sea soon after hatching, results in a reduced homing ability. From a conservation management perspective the Gotland brown trout streams should be regarded as a population system where the vitality and survival of brown trout in one stream is dependent on the opportunity of contact and exchange of individuals from other streams.     

Local genetic structure in red grouse (Lagopus lagopus scoticus): evidence from microsatellite DNA markers

Allelic variation at seven hypervariable tri- and tetranucleotide microsatellite loci was used to determine levels of population differentiation between 14 populations of red grouse ( Lagopus lagopus scoticus ) in northeast Scotland, UK. Despite the potential for long-distance dispersal in grouse, and a semicontinuous habitat, significant population divergence was observed (mean R _ST = 0.153; P < 0.01) and an isolation-by-distance effect detected (Mantel test: P < 0.001). Examination of the spatial trend in principal component scores derived from allele frequencies among populations highlighted a barrier to gene flow that was confounding a simple isolation-by-distance effect. This barrier corresponded to an area of unsuitable habitat for grouse associated with a river system that bisected the study area. Mean genetic relatedness was higher for males than for females in all but one of the study populations, suggesting that the territorial behaviour and natal philopatry displayed by cocks have a manifold effect in generating the observed spatial genetic structure. Lower female relatedness values suggest a higher level of female-mediated gene flow, which is sufficient to prevent the loss of genetic variation from within populations and the onset of inbreeding effects. The potential consequences of local subdivision for red grouse populations are discussed.     

Genetic differentiation associated with commercial traffic in the Polynesian mosquito, Aedes polynesiensis Marks 1951

The population structure of the Polynesian mosquito Aedes polynesiensis was investigated using electrophoretic data from two polymorphic protein loci. Considerable differentiation was observed both within and between islands in different archipelagos (Society, Tuamotu, Austral). Gene flow evaluated by F_st estimates was independent of geographic distance between islands but related to commercial traffic intensity. The results are discussed in view of recent findings on the variability of susceptibility to insecticides and of suitability as a vector for the nematode Wuchereria bancrofti.     

Genetic variation and stock structure of school mackerel and spotted mackerel in northern Australian waters

The total mean sample heterozygosity calculated from eight polymorphic loci was 0·172 (0·047 S.E.) (F_st 0·025) for school mackerel Scomberomorus queenslandicus , and for spotted mackerel S. munroi was 0·110 (0·074 S.E.) (F_st 0·038). There was no evidence of temporal variation for either species as significant genetic differences were not detected between months or year classes within areas. Spatially, school mackerel have a complex stock structure, with stocks being associated with large embayments. In contrast, spotted mackerel appear to comprise a single stock in Australian east coast waters. Both species showed a significant pattern of stock structure between Australian east coast and northern (Arafura Sea) samples.     

Low population genetic structuring of two cryptic bat species suggests their migratory behaviour in continental Europe

Although two cryptic pipistrelle bat species, Pipistrellus pipistrellus and Pipistrellus pygmaeus , belong among the most common bat species in Europe, it is still unclear whether they can migrate over long distances between summer and winter roosts. Long-distance migratory species may be expected to show low levels of genetic structuring in large areas due to regular mixing of the gene pool by mating that occurs during migration and/or hibernation. Conversely, the dispersal of gametes in sedentary species is spatially restricted, populations are more genetically structured, and isolation by relatively short distance is visible. By analysing diversity of highly variable microsatellites within and among summer colonies of both studied species in central Europe, we found that differentiation between populations is very weak. Both classical F _ST and Bayesian clustering approach failed to detect genetic structure among colonies and there was no significant isolation-by-distance pattern. The analyses of relatedness, however, revealed that individuals within colonies are more related than random suggesting philopatry of at least one sex. The results were very similar for the two species. The high level of gene flow among central European populations, even on large geographic distances, is discussed in relation with migrations, dispersal, and mating behaviour.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 103–114.     

A nuclear DNA basis for shelf- and bank-scale population structure in northwest Atlantic cod (Gadus morhua): Labrador to Georges Bank

Variation at five microsatellite DNA loci scored in ≈ 1300 individuals provided evidence of genetic structure among 14 cod ( Gadus morhua ) populations spanning the range of the species in the northwest (NW) Atlantic. Using D _A and D _SW measures of genetic distance, as well as F _ST and R _ST measures of population structure, differences are revealed among populations at continental shelf scales (NE Newfoundland Shelf, Grand Banks, Flemish Cap, Scotian Shelf, Georges Bank) where regions are separated by submarine saddles, channels and trenches. However, we also provide evidence of genetic structure at spawning-bank scales consistent with variation in oceanographic features and in the spatiotemporal distribution of spawning, each of which may represent barriers to gene flow among geographically contiguous populations inhabiting a highly advective environment. The differences described are consistent with postdispersal spawning fidelity to natal areas, a behaviour that may be facilitated by topographically induced gyre-like circulations that can act as retention mechanisms. Significant degrees of substructure among neighbouring and contiguous cod populations may be most easily explained by the associated oceanographic features and processes that conceivably form the template for the evolution of the structure. We suggest that bathymetric and hydrodynamic structure represents a rational starting point for developing hypotheses to examine the processes that lead to the genetic structuring of marine fish species.     

Is the last glaciation the only relevant event for the present genetic population structure of the meadow brown butterfly Maniola jurtina (Lepidoptera: Nymphalidae)?

Phylogeographical studies are available for a considerable number of European species, but few analyses exist for temperate species with very large and fairly continuous populations that are also absent from Northern Europe. Therefore, we studied the butterfly Maniola jurtina as a model for this group. The species has two major genetic lineages (mean genetic distance between lineages: 0.033; F _CT: 0.052), most probably evolving in glacial differentiation centres in the western and eastern Mediterranean. The onset of this differentiation might have been the beginning of the last glacial stage maximum some 40 kyr bp . A hybrid zone between these two lineages exists in western Central Europe. No genetic substructures have been found within the two lineages ( F _SC: 0.017) and average genetic distances are very small. Therefore, it is highly probable that postglacial expansion was of the phalanx type. There is, at most, very limited differentiation at regional and local scales. However, the genetic diversity within populations is high (means: A : 2.68; H _E: 17.2%; P : 78%), as would be predicted for such a common species. Comparison of these results with a published allozyme analysis revealed a similar phylogeographical pattern, but lower genetic diversity in the latter. Morphological patterns of wings and genitalia show similar geographical patterns as allozyme data.   © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 85 , 419–431.     

Did Polyommatus icarus (Lepidoptera: Lycaenidae) have distinct glacial refugia in southern Europe? Evidence from population genetics

Pleistocene climatic oscillations strongly influenced the genetic composition of many species which are often divided into several genetic lineages. In this context, we studied the allozymes of a common and widely distributed butterfly, the common blue Polyommatus icarus, over a large part of Europe. The species had a rather high genetic diversity within populations with a strikingly high mean number of alleles per locus (2.98). In contrast, differentiation between populations was very low ( F _ST: 0.0187). Only a marginal trend of decline in genetic diversity from the south to the north was observed. Isolation-by-distance existed on a European scale ( r =  0.826), but not at a regional level. Regional differentiation between populations in western Germany was extremely low ( F _ST: 0.0041). It is probable that P. icarus was widely distributed in the Mediterranean region during the last ice age and expanded into central Europe in the postglacial period without major genetic erosion. Moderate present and past gene flow in an intact metapopulation structure may have occurred on local, regional and perhaps even continental scales.  © 2003 The Linnean Society of London, Biological Journal of the Linnean Society , 2003, 80 , 529–538.     

Geographic patterns of genetic differentiation and plumage colour variation are different in the pied flycatcher (Ficedula hypoleuca)

The pied flycatcher is one of the most phenotypically variable bird species in Europe. The geographic variation in phenotypes has often been attributed to spatial variation in selection regimes that is associated with the presence or absence of the congeneric collared flycatcher. Spatial variation in phenotypes could however also be generated by spatially restricted gene flow and genetic drift. We examined the genetic population structure of pied flycatchers across the breeding range and applied the phenotypic Q _ST ( P _ST)– F _ST approach to detect indirect signals of divergent selection on dorsal plumage colouration in pied flycatcher males. Allelic frequencies at neutral markers were found to significantly differ among populations breeding in central and southern Europe whereas northerly breeding pied flycatchers were found to be one apparently panmictic group of individuals. Pairwise differences between phenotypic ( P _ST) and neutral genetic distances ( F _ST) were positively correlated after removing the most differentiated Spanish and Swiss populations from the analysis, suggesting that genetic drift may have contributed to the observed phenotypic differentiation in some parts of the pied flycatcher breeding range. Differentiation in dorsal plumage colouration however greatly exceeded that observed at neutral genetic markers, which indicates that the observed pattern of phenotypic differentiation is unlikely to be solely maintained by restricted gene flow and genetic drift.     

Reproductive biology of a mycoheterotrophic species, Burmannia wallichii (Burmanniaceae)

The breeding system of a mycoheterotrophic species, Bumannia wallichii (Miers) Hook. f. (Burmanniaceae), is assessed using field observations, floral anatomy, pollen histochemistry and isozyme electrophoresis. The structure of the flower effectively prevents pollinators from accessing the stamens, and no pollinators were observed visiting natural populations in Hong Kong. The pollen is starch-rich but lipid-poor. Analysis of isozyme variation revealed a small proportion of polymorphic loci, a low level of total heterozygosity ( H _T= 0.1972), and a very high overall fixation index of populations ( F _IT= 0.6602). These results strongly suggest that B. wallichii is predominantly self-pollinated. The coefficient of gene differentiation ( G _ST) is low (0.05), although it is suggested that this is due to genetic homogenization of populations consonant with the formation of numerous tiny 'diaspora' seeds that are wind dispersed.     

Genetic structure of Lake Magadi tilapia populations

The control region of the mitochondrial DNA haplotype frequencies were significantly different among the two separate lagoon populations of Oreochromis alcalicus grahami in Lake Magadi and of O. a. alcalicus from lake Natron, and DNA fingerprint similarity indices were significantly higher for intra-population comparisons of the two Magadi lagoon populations and the Lake Natron population than the inter-population similarity indices among these populations. A modified F _st measure indicated population sub-division and the phylogeographic partitioning of the VNTR fragments observed were unique to specific populations further indicating substantial genetic differentiation. The lagoon populations within Lake Magadi demonstrated the same degree of genetic differentiation as either of these populations did to the outgroup (the Lake Natron population). There appears to be limited gene flow between Lake Magadi tilapia populations and this population structure has important implications for protecting locally adapted populations within this unique ecosystem.     

Genetic Diversity and Population Structure of Dioscorea tokoro Makino, a Dioecious Climber

Abstract Dioscorea tokoro Makino is a herbaceous climber species widespread in East Asia. Genetic structure of a natural population of D. tokoro was examined employing starch gel electrophoresis of allozymes. Genotypes of seven loci were studied for 1,128 individuals. Twenty-six populations located mainly in the Kinki district of Japan were subgrouped into four large clusters by the geographical distribution of alleles. The D. tokoro population was revealed to contain greater total genetic diversity ( H_T =0.282) and higher intrapopulational genetic diversity ( H_S =0.258) than other outcrossing species for which data are available. On the other hand, interpopulational differentiation ( G_ST =0.096) was smaller than in other outcrossers. For the heterozygosity deficiency observed ( F_IT =0.125), population subdivision ( F_ST =0.096) and inbreeding within the population ( F_IS =0.067) were revealed to contribute to the same extent. From these F -statistics, the migration rate among subpopulations and the rate of between-relative matings were estimated. Overall results on the genetic structure of the D. tokoro population indicated a high gene flow among its subpopulations, and this may be the consequence of its life form as a climber and its habitat in a disturbed environment. During the study, the geographical cline of Pgi allele frequencies was observed. This finding was supposed to be the result of the selection imposed on Pgi by the temperature differences between localities.