ISOLATION BY DISTANCE

Evolutionary processes in a continental island system: molecular phylogeography of the Aegean Nigella arvensis alliance (Ranunculaceae) inferred from chloroplast DNA

Continental shelf island systems, created by rising sea levels, provide a premier setting for studying the effects of past fragmentation, dispersal, and genetic drift on taxon diversification. We used phylogeographical (nested clade) and population genetic analyses to elucidate the relative roles of these processes in the evolutionary history of the Aegean Nigella arvensis alliance (= 'coenospecies'). We surveyed chloroplast DNA (cpDNA) variation in 455 individuals from 47 populations (nine taxa) of the alliance throughout its core range in the Aegean Archipelago and surrounding mainland areas of Greece and Turkey. The study revealed the presence of three major lineages, with largely nonoverlapping distributions in the Western, Central, and Eastern Aegean. There is evidence supporting the idea that these major lineages evolved in situ from a widespread (pan-Aegean) ancestral stock as a result of multiple fragmentation events, possibly due to the influence of post-Messinian sea flooding, Pleistocene eustatic changes and corresponding climate fluctuations. Over-sea dispersal and founder events appear to have played a rather insignificant role in the group's history. Rather, all analytical approaches identified the alliance as an organism group with poor seed dispersal capabilities and a susceptibility to genetic drift. In particular, we inferred that the observed level of cpDNA differentiation between Kikladian island populations of Nigella degenii largely reflects population history, (viz. Holocene island fragmentation) and genetic drift in the near absence of seed flow since their time of common ancestry. Overall, our cpDNA data for the N. arvensis alliance in general, and N. degenii in particular, indicate that historical events were important in determining the phylogeographical patterns seen, and that genetic drift has historically been relatively more influential on population structure than has cytoplasmic gene flow.

Colonization and dispersal in a social species, the Bechstein's bat (Myotis bechsteinii)

Metapopulation genetic models consider that colonization and dispersal are distinct behaviours. However, whether colonization and dispersal indeed reflect different biological processes in nature is unclear. One possibility to test this assumption is to assess patterns of autosomal and mitochondrial genetic structure in species with strict female philopatry, such as the communally breeding Bechstein's bat. In this species, mitochondrial DNA can spread only when females establish new colonies, and autosomal DNA is transmitted among colonies only when females mate with solitary males born in foreign colonies. Investigating the genetic structure among 37 colonies, we found that autosomal genes followed an island model on a regional scale and a model of isolation by distance on a larger geographical scale. In contrast, mitochondrial genetic structure revealed no pattern of isolation by distance at a large scale but exhibited an effect of ecological barriers on a regional scale. Our results provide strong empirical evidence that colonization and dispersal do not follow the same behavioural rules in this bat, supporting the assumption of metapopulation genetic models.

Marine landscapes and population genetic structure of herring (Clupea harengus L.) in the Baltic Sea

Numerically small but statistically significant genetic differentiation has been found in many marine fish species despite very large census population sizes and absence of obvious barriers to migrating individuals. Analyses of morphological traits have previously identified local spawning groups of herring (Clupea harengus L.) in the environmentally heterogeneous Baltic Sea, whereas allozyme markers have not revealed differentiation. We analysed variation at nine microsatellite loci in 24 samples of spring-spawning herring collected at 11 spawning locations throughout the Baltic Sea. Significant temporal differentiation was observed at two locations, which we ascribe to sympatrically spawning but genetically divergent 'spawning waves'. Significant differentiation was also present on a geographical scale, though pairwise F(ST) values were generally low, not exceeding 0.027. Partial Mantel tests showed no isolation by geographical distance, but significant associations were observed between genetic differentiation and environmental parameters (salinity and surface temperature) (0.001 < P < or = 0.099), though these outcomes were driven mainly by populations in the southwestern Baltic Sea, which also exhibits the steepest environmental gradients. Application of a novel method for detecting barriers to gene flow by combining geographical coordinates and genetic differentiation allowed us to identify two zones of lowered gene flow. These zones were concordant with the separation of the Baltic Sea into major basins, with environmental gradients and with differences in migration behaviour. We suggest that similar use of landscape genetics approaches may increase the understanding of the biological significance of genetic differentiation in other marine fishes.

Lack of taxonomic differentiation in an apparently widespread freshwater isopod morphotype (Phreatoicidea: Mesamphisopidae: Mesamphisopus) from South Africa

The unambiguous identification of phreatoicidean isopods occurring in the mountainous southwestern region of South Africa is problematic, as the most recent key is based on morphological characters showing continuous variation among two species: Mesamphisopus abbreviatus and M. depressus. This study uses variation at 12 allozyme loci, phylogenetic analyses of 600 bp of a COI (cytochrome c oxidase subunit I) mtDNA fragment and morphometric comparisons to determine whether 15 populations are conspecific, and, if not, to elucidate their evolutionary relationships. Molecular evidence suggested that the most easterly population, collected from the Tsitsikamma Forest, was representative of a yet undescribed species. Patterns of differentiation and evolutionary relationships among the remaining populations were unrelated to geographic proximity or drainage system. Patterns of isolation by distance were also absent. An apparent disparity among the extent of genetic differentiation was also revealed by the two molecular marker sets. Mitochondrial sequence divergences among individuals were comparable to currently recognized intraspecific divergences. Surprisingly, nuclear markers revealed more extensive differentiation, more characteristic of interspecific divergences. This disparity and the mosaic pattern of differentiation may be driven by stochastic population crashes and genetic bottlenecks (caused by seasonal habitat fluctuations), coupled with genetic drift.

Unifying the spatial population dynamics and molecular evolution of epidemic rabies virus

Infectious disease emergence is under the simultaneous influence of both genetic and ecological factors. Yet, we lack a general framework for linking ecological dynamics of infectious disease with underlying molecular and evolutionary change. As a model, we illustrate the linkage between ecological and evolutionary dynamics in rabies virus during its epidemic expansion into eastern and southern Ontario. We characterized the phylogeographic relationships among 83 isolates of fox rabies virus variant using nucleotide sequences from the glycoprotein-encoding glycoprotein gene. The fox rabies virus variant descended as an irregular wave with two arms invading from northern Ontario into southern Ontario over the 1980s and 1990s. Correlations between genetic and geographic distance suggest an isolation by distance population structure for the virus. The divergence among viral lineages since the most recent common ancestor correlates with position along the advancing wave front with more divergent lineages near the origin of the epidemic. Based on divergence from the most recent common ancestor, the regional population can be partitioned into two subpopulations, each corresponding to an arm of the advancing wave. Subpopulation A (southern Ontario) showed reduced isolation by distance relative to subpopulation B (eastern Ontario). The temporal dynamics of subpopulation A suggests that the subregional viral population may have undergone several smaller waves that reduced isolation by distance. The use of integrated approaches, such as the geographical analysis of sequence variants, coupled with information on spatial dynamics will become indispensable aids in understanding patterns of disease emergence.

Landscape genetics of the blotched tiger salamander (Ambystoma tigrinum melanostictum)

The field of landscape genetics has great potential to identify habitat features that influence population genetic structure. To identify landscape correlates of genetic differentiation in a quantitative fashion, we developed a novel approach using geographical information systems analysis. We present data on blotched tiger salamanders (Ambystoma tigrinum melanostictum) from 10 sites across the northern range of Yellowstone National Park in Montana and Wyoming, USA. We used eight microsatellite loci to analyse population genetic structure. We tested whether landscape variables, including topographical distance, elevation, wetland likelihood, cover type and number of river and stream crossings, were correlated with genetic subdivision (F(ST)). We then compared five hypothetical dispersal routes with a straight-line distance model using two approaches: (i) partial Mantel tests using Akaike's information criterion scores to evaluate model robustness and (ii) the BIOENV procedure, which uses a Spearman rank correlation to determine the combination of environmental variables that best fits the genetic data. Overall, gene flow appears highly restricted among sites, with a global F(ST) of 0.24. While there is a significant isolation-by-distance pattern, incorporating landscape variables substantially improved the fit of the model (from an r2 of 0.3 to 0.8) explaining genetic differentiation. It appears that gene flow follows a straight-line topographic route, with river crossings and open shrub habitat correlated with lower F(ST) and thus, decreased differentiation, while distance and elevation difference appear to increase differentiation. This study demonstrates a general approach that can be used to determine the influence of landscape variables on population genetic structure.

Invasion genetics of the Eurasian spiny waterflea: evidence for bottlenecks and gene flow using microsatellites

The Eurasian spiny waterflea (Bythotrephes longimanus) is a predacious zooplankter that has increased its range in Europe and is rapidly invading inland water-bodies throughout North America's Great Lakes region. To examine the genetics of these invasions, we isolated five microsatellite DNA loci with between 5 and 19 alleles per locus. We sampled three populations where B. longimanus has been historically present (Switzerland, Italy, and Finland) as well as an introduced European population (the Netherlands) and three North American populations (Lakes Erie, Superior, Shebandowan). Consistent with a bottleneck during colonization (i.e. founder effect), average heterozygosities of the four European populations ranged from 0.310 to 0.599, and were higher than that of three North American populations (0.151-0.220). Pairwise F(ST) estimates among North American populations (0.002-0.063) were not significantly different from zero and were much lower than among European populations (0.208-0.474). This is consistent with a scenario of high gene flow among North American populations relative to that of European ones. Contrary to an invasion bottleneck, however, Erie and Superior populations contained similar numbers of rare alleles as European populations. Assignment tests identified several migrant genotypes in all introduced populations (the Netherlands, Erie, Superior, Shebandowan), but rarely in native ones (Switzerland, Italy and Finland). A large number of genotypes from North America were assigned to our Italian population suggesting a second, previously unidentified, invasion source somewhere in the region of northern Italy. Together, our results support an invasion bottleneck for North American populations that has been largely offset by gene flow from multiple native sources, as well as gene flow among introduced populations.

Streams over mountains: influence of riparian connectivity on gene flow in the Pacific jumping mouse (Zapus trinotatus)

In species affiliated with heterogeneous habitat, we expect gene flow to be restricted due to constraints placed on individual movement by habitat boundaries. This is likely to impact both individual dispersal and connectivity between populations. In this study, a GIS-based landscape genetics approach was used, in combination with fine-scale spatial autocorrelation analysis and the estimation of recent intersubpopulation migration rates, to infer patterns of dispersal and migration in the riparian-affiliated Pacific jumping mouse (Zapus trinotatus). A total of 228 individuals were sampled from nine subpopulations across a system of three rivers and genotyped at eight microsatellite loci. Significant spatial autocorrelation among individuals revealed a pattern of fine-scale spatial genetic structure indicative of limited dispersal. Geographical distances between pairwise subpopulations were defined following four criteria: (i) Euclidean distance, and three landscape-specific distances, (ii) river distance (distance travelled along the river only), (iii) overland distance (similar to Euclidean, but includes elevation), and (iv) habitat-path distance (a least-cost path distance that models movement along habitat pathways). Pairwise Mantel tests were used to test for a correlation between genetic distance and each of the geographical distances. Significant correlations were found between genetic distance and both the overland and habitat-path distances; however, the correlation with habitat-path distance was stronger. Lastly, estimates of recent migration rates revealed that migration occurs not only within drainages but also across large topographic barriers. These results suggest that patterns of dispersal and migration in Pacific jumping mice are largely determined by habitat connectivity.

Strong spatial genetic structure in peripheral but not core populations of Sitka spruce [Picea sitchensis (Bong.) Carr.]

We examined spatial genetic structure within eight populations of Sitka spruce classified as core or peripheral based on ecological niche, and continuous or disjunct based on species distribution. In each population, 200 trees were spatially mapped and genotyped for eight cDNA-based sequence tagged site (STS) codominant markers. Spatial autocorrelation was assessed by estimating p(ij), the average co-ancestry coefficient, between individuals within distance intervals. The distribution of alleles and genotypes within core populations was almost random, with nonsignificant co-ancestry values among trees as close as 50 m in core populations. In contrast, the distribution of alleles and genotypes within peripheral populations revealed an aggregation of similar multilocus genotypes, with co-ancestry values greater than 0.20 among trees up to 50 m apart and significant, positive values between trees up to 500 m. The relatively high density of reproductive adults in core populations may lead to highly overlapping seed shadows that limit development of spatial genetic structure. However, in peripheral populations with a lower density of adults, the distribution of alleles and genotypes was highly structured, likely due to offspring establishment near maternal trees and subsequent biparental inbreeding, as well as more recent population establishment at the leading edge of post-Pleistocene range expansion. Conserving genetic diversity in peripheral populations may require larger reserves for in situ conservation than required in core populations. These data on spatial genetic structure can be used to provide guidance for sampling strategies for both ex situ conservation and research collections.

Clonal and fine-scale genetic structure in populations of a restricted Korean endemic, Hosta jonesii (Liliaceae) and the implications for conservation

BACKGROUND AND AIMS

In plant populations the magnitude of spatial genetic structure of apparent individuals (including clonal ramets) can be different from that of sexual individuals (genets). Thus, distinguishing the effects of clonal versus sexual individuals in population genetic analyses could provide important insights for evolutionary biology and conservation. To investigate the effects of clonal spread on the fine-scale spatial genetic structure within plant populations, Hosta jonesii (Liliaceae), an endemic species to Korea, was chosen as a study species.

METHODS

Using allozymes as genetic markers, spatial autocorrelation analysis of ramets and of genets was conducted to quantify the spatial scale of clonal spread and genotype distribution in two populations of H. jonesii.

KEY RESULTS

Join-count statistics revealed that most clones are significantly aggregated at < 3-m interplant distance. Spatial autocorrelation analysis of all individuals resulted in significantly higher Moran's I values at 0-3-m interplant distance than analyses of population samples in which clones were excluded. However, significant fine-scale genetic structure was still observed when clones were excluded.

CONCLUSIONS

These results suggest that clones enhance the magnitude of spatial autocorrelation due to localized clonal spread. The significant fine-scale genetic structure detected in samples excluding clones is consistent with the biological and ecological traits exhibited by H. jonesii including bee pollination and limited seed dispersal. For conservation purposes, genetic diversity would be maximized in local populations of H. jonesii by collecting or preserving individuals that are spaced at least 5 m apart.

Genetic structure and environmental heterogeneity in the European hake (Merluccius merluccius)

This study aimed at assessing the genetic structure and the state of the stocks of the European hake (Merluccius merluccius). To this end, 15 samples were taken from the whole range of the species and analysed using allozymes. Since 11 samples were taken from the poorly studied Mediterranean Sea, the results obtained provided a complete picture of the hake's genetic structure and an initial insight into its relationships with environmental features. Atlantic and Mediterranean hake populations are separated by the Almeria-Oran front. This area has been proved to be the boundary between Atlantic and Mediterranean stocks of many marine organisms, but some doubt exists concerning the efficaciousness of the local gyres as barriers to the gene flow. Our data have evidenced a latitudinal cline at loci Gapdh and Gpi-2 within the Mediterranean Sea, with a further steep change across the Almeria-Oran front. The genetic pattern showed a strong correlation with the values of the salinity both at the surface and at -320 m and of the salinity + temperature at the surface, suggesting a role for these parameters in maintaining the genetic differentiation among the two population groups through selective processes. Finally, the levels of genetic variability were found to be slightly lower in the depleted Atlantic stock than in the Mediterranean one.

A separation-of-timescales approach to the coalescent in a continuous population

This article presents an analysis of a model of isolation by distance in a continuous, two-dimensional habitat. An approximate expression is derived for the distribution of coalescence times for a pair of sequences sampled from specific locations in a rectangular habitat. Results are qualitatively similar to previous analyses of isolation by distance, but account explicitly for the location of samples relative to the habitat boundaries. A separation-of-timescales approach takes advantage of the fact that the sampling locations affect only the recent coalescent behavior. When the population size is larger than the number of generations required for a lineage to cross the habitat range, the long-term genealogical process is reasonably well described by Kingman's coalescent with time rescaled by the effective population size. This long-term effective population size is affected by the local dispersal behavior as well as the geometry of the habitat. When the population size is smaller than the time required to cross the habitat, deep branches in the genealogy are longer than would be expected under the standard neutral coalescent, similar to the pattern expected for a panmictic population whose population size was larger in the past.

Genetic structure and migration from mainland to island populations inAbies proceraRehd

Noble fir (Abies procera Rehd) is a narrowly distributed conifer with a typical mainland-island structure of natural distribution. Here, we examined the genetic structure of populations native to the Pacific coast from Oregon to Washington (5 island and 16 mainland populations) with 14 polymorphic allozyme loci. A general method for estimating the number of unidirectional migrants from the mainland to island populations is presented in terms of the relation of average heterozygosity between the mainland and island populations. The results indicated that there were substantial island–mainland population differentiations (Fst = 0.107±0.029~0.154±0.039) but small differentiation within the mainland/submainland populations (0.037±0.008~0.054±0.010). Significant isolation by distance existed among the island–mainland populations and among the populations in Washington submainland. Four islands investigated received different numbers of migrants from the mainland/submainland. The southern island populations received a smaller number of migrants from the mainland but had greater genetic diversity, implying that there could be introgression with A. magnifica and (or) they represented possible glacial refuges and had expanded northwards after the last glaciations. The island populations close to the Pacific coast were more likely mainland-dependent.Key words: Abies procera Rehd, gene flow, island, mainland, genetic structure.

Conservation of genetic diversity in the Japanese Black cattle population by the construction of partially isolated lines

In the Japanese Black cattle population, five genetically divergent sub-populations have played important roles as suppliers of breeding stocks to the entire breed. We supposed a situation where five lines were constructed from each of the five subpopulations, and applied to this set of lines a management plan to conserve genetic diversity. Assuming that the male migration pattern among the lines followed the island model, we assessed the optimum male migration rate and required male numbers in each line, satisfying the following three conditions simultaneously: (i) the rate of inbreeding in each line was below 0.01 per generation; (ii) at least 97% of the initial genetic diversity was preserved after 10 generations; and (iii) more than 50% of the genes in an initial line were retained in the line after 10 generations. We found that approximately one breeding male should be selected per year and one breeding male should be exchanged among the lines per generation in each line to satisfy these three conditions. Numerical analysis with the migration rates actually observed among the five subpopulations demonstrated that the initial genetic differentiation among the lines was rapidly decayed by an asymmetrical migration pattern. For a successful plan, migration among lines should be strictly managed.

F-statistics under alternation of sexual and asexual reproduction: a model and data from schistosomes (platyhelminth parasites)

Accurate inferences on population genetics data require a sound underlying theoretical null model. Nearly nothing is known about the gene dynamics of organisms with complex life cycles precluding any biological interpretation of population genetics parameters. In this article, we used an infinite island model to derive the expectations of those parameters for the life cycle of a dioecious organism obligatorily alternating sexual and asexual reproductions as it is the case for schistosomes (plathyhelminth parasites). This model allowed us to investigate the effects of the degree of mixing among individuals coming from different subpopulations at each new generation (represented in the model by the migration rates before and after clonal reproductions) and the variance in the reproductive success of individuals during the clonal phase. We also consider the effects of different migration rates and degrees of clonal reproductive skew between male and female individuals. Results show that the variance in the reproductive success of clones is very important in shaping the distribution of the genetic variability both within and among subpopulations. Thus, higher variance in the reproductive success of clones generates heterozygous excesses within subpopulations and also increases genetic differentiation between them. Migration occurring before and after asexual reproduction has different effects on the patterns of F(IS) and F(ST). When males and females display different degrees of reproductive skew or migration rates, we observe differences in their respective population genetic structure. While results of the model apply to any organism alternating sexual and clonal reproductions (e.g. all parasitic trematodes, many plants, and all aphididae), we finally confront some of these theoretical expectations to empirical data from Schistosoma mansoni infecting Rattus rattus in Guadeloupe.

Measuring the genetic structure of the pollen pool as the probability of paternal identity

Contemporary pollen flow in forest plant species is measured by the probability of paternal identity (PPI) for two randomly sampled offspring, drawn from a single female, and contrasting that with PPI for two random offspring, drawn from different females. Two different estimation strategies have emerged: (a) an indirect approach, using the 'genetic structure' of the pollen received by different mothers and (b) a direct approach, based on parentage analysis. The indirect strategy is somewhat limited by the assumptions, but is widely useful. The direct approach is most appropriate where a large majority of the true fathers can be identified exactly, which is sometimes possible with high-resolution SSR markers. Using the parentage approach, we develop estimates of PPI, showing that the obvious estimates are severely biased, and providing an unbiased alternative. We then illustrate the methods with SSR data from a 36-tree isolated population of Pinus sylvestris from the Meseta region of Spain, for which categorical paternity assignment was available for over 95% of offspring. For all the females combined, we estimate that PPI=0.0425, indicating uneven male reproductive contributions. Different (but overlapping) arrays of males pollinate different females, and for the average female, PPI=0.317, indicating substantial 'pollen structure' for the population. We also relate the direct measures of PPI to those available from indirect approaches, and show that they are generally comparable.

Estimating dispersal from short distance spatial autocorrelation

A series of theoretical studies has formed a strong connection between spatial statistics observed in populations and summary measures of the amount of dispersal. Synthesized, these developments allow dispersal to be indirectly estimated from standing spatial patterns of genetic variation under a range of conditions broad enough to be likely met in most populations of either plants or animals. The spatial correlations at the shortest distances are particularly robust to range of conditions and have disproportionately high statistical power. This review integrates theoretical results in a way that maximizes robustness and flexibility in the use of short distance autocorrelation to estimate Wright's neighborhood size, or the total variance in dispersal distances. Empirical guidelines are developed that are meant to be as practical and broad as possible. The guidelines focus on Moran's I-statistics for diploid genotypes converted to allele frequencies, but are also extended to or compared with several other approaches.

Genetic differentiation and gene flow among populations of the alpine butterfly, Parnassius smintheus, vary with landscape connectivity

Levels of gene flow among populations vary both inter- and intraspecifically, and understanding the ecological bases of variation in levels of gene flow represents an important link between the ecological and evolutionary dynamics of populations. The effects of habitat spatial structure on gene flow have received considerable attention; however, most studies have been conducted at a single spatial scale and without background data on how individual movement is affected by landscape features. We examined the influence of habitat connectivity on inferred levels of gene flow in a high-altitude, meadow-dwelling butterfly, Parnassius smintheus. For this species, we had background data on the effects of landscape structure on both individual movement and on small-scale population genetic differentiation. We compared genetic differentiation and patterns of isolation by distance, based on variation at seven microsatellite loci, among three regions representing two levels of connectivity of high-altitude, nonforested habitats. We found that reduced connectivity of habitats, resulting from more forest cover at high altitudes, was associated with greater genetic differentiation among populations (higher estimated FST), a breakdown of isolation by distance, and overall lower levels of inferred gene flow. These observed differences were consistent with expectations based on our knowledge of the movement behaviour of this species and on previous population genetic analyses conducted at the smaller spatial scale. Our results indicate that the role of gene flow may vary among groups of populations depending on the interplay between individual movement and the structure of the surrounding landscape.

Genetic structure of the tick Ornithodoros coriaceus (Acari: Argasidae) in California, Nevada, and Oregon

The argasid tick Ornithodoros coriaceus (Koch) is the only confirmed vector of epizootic bovine abortion (EBA) in the United States. The disease and its tick vector have historically been reported in the foothills of the Sierra Nevada Mountains and coast ranges of California. In the past two decades, the range of EBA has apparently expanded into southern Oregon and northern Nevada. Possible explanations for this expansion include 1) increased recognition and reporting of EBA in these regions; 2) widespread movement of tick-infested and EBA-infected hosts with subsequent colonization of these regions by infected ticks; and 3) widespread movement of the EBA agent, independent of tick movements, into extant tick populations in these new regions. The current study was performed to evaluate these hypotheses by examining patterns of variability in a 420-bp segment of the 16S mitochondrial rDNA gene sequence among 210 O. coriaceus individuals from 14 sites in California, Oregon, and Nevada. Sixty-three unique haplotypes were identified in the ticks tested, with 84% of the sequence variation attributable to among-population variation and 16% to within-population variation. A majority of the haplotypes were unique to their particular collection site, whereas only four collection sites shared haplotypes. Overall, very little evidence of gene flow among tick populations was detected, making it unlikely that widespread tick movement had introduced O. coriaceus and the EBA agent into new regions.

Inbreeding and kinship in the ant Plagiolepis pygmaea

In ants the presence of multiple reproductive queens (polygyny) decreases the relatedness among workers and the brood they rear, and subsequently dilutes their inclusive fitness benefits from helping. However, adoption of colony daughters, low male dispersal in conjunction with intranidal (within nest) mating and colony reproduction by budding may preserve local genetic differences, and slow down the erosion of relatedness. Reduced dispersal and intranidal mating may, however, also lead to detrimental effects owing to competition and inbreeding. We studied mating and dispersal patterns, and colony kinship in three populations of the polygynous ant Plagiolepis pygmaea using microsatellite markers. We found that the populations were genetically differentiated, but also a considerable degree of genetic structuring within populations. The genetic viscosity within populations can be attributed to few genetically homogeneous colony networks, which presumably have arisen through colony reproduction by budding. Hence, selection may act at different levels, the individuals, the colonies and colony networks. All populations were also significantly inbred (F=0.265) suggesting high frequencies of intranidal mating and low male dispersal. Consequently the mean regression relatedness among workers was significantly higher (r = 0.529-0.546) than would be expected under the typically reported number (5-35) of queens in nests of the species. Furthermore, new queens were mainly recruited from their natal or a neighbouring related colony. Finally, the effective number of queens coincided with that found upon excavation, suggesting low reproductive skew.

Invasion genetics of a freshwater mussel (Dreissena rostriformis bugensis) in eastern Europe: high gene flow and multiple introductions

In recent years, the quagga mussel, Dreissena rostriformis bugensis, native to the Dnieper and Bug Limans of the northern Black Sea, has been dispersed by human activities across the basin, throughout much of the Volga River system, and to the Laurentian Great Lakes. We used six published microsatellite markers to survey populations throughout its native and introduced range to identify relationships among potential source populations and introduced ones. Mussels from 12 sites in Eurasia, including the central Caspian Sea and one in North America (Lake Erie), were sampled. Field surveys in the Volga River basin suggested that the species first colonized the middle reach of the river near Kubyshev Reservoir, and thereafter spread both upstream and downstream. Evidence of considerable gene flow among populations was observed and genetic diversity was consistent with a larger, metapopulation that has not experienced bottlenecks or founder effects. We propose that high gene flow, possibly due to multiple invasions, has facilitated establishment of quagga mussel populations in the Volga River system. The Caspian Sea population (D. rostriformis rostriformis (=distincta)) was genetically more distinct than other populations, a finding that may be related to habitat differences. The geographical pattern of genetic divergence is not characteristic of isolation-by-distance but, rather, of long-distance dispersal, most likely mediated by commercial ships' ballast water transfer.

Peromyscuspopulations and theirCuterebraparasites display congruent phylogeographical structure

The relationships between populations of the Deer Mouse (Peromyscus maniculatus) and the White-footed Mouse (P. leucopus) and their respectiveCuterebraparasites were examined. Population genetic structure of hosts and parasites was inferred using cytochrome oxidase mitochondrial sequences of specimens from 7 populations. Genetic analyses revealed that isolation-by-distance applies forP. maniculatusand its associated parasite (C. grisea). A significant correlation was also observed between the genetic distances of these host and parasite species. Furthermore, populations ofP. maniculatusandC. griseafrom the North and South shores of the St Lawrence River were found to be significantly different. This structure may be explained by the St Lawrence River being a dispersal barrier for both species. A robust analysis of the other species pair (P. leucopusandC. fontinella) could not be performed because of limited sample sizes.

Evolutionary divergence and possible incipient speciation in post-glacial populations of a cosmopolitan aquatic plant

Habitat configuration is expected to have a major influence on genetic exchange and evolutionary divergence among populations. Aquatic organisms occur in two fundamentally different habitat types, the sea and freshwater lakes, making them excellent models to study the contrasting effects of continuity vs. isolation on genetic divergence. We compared the divergence in post-glacial populations of a cosmopolitan aquatic plant, the pondweed Potamogeton pectinatus that simultaneously occurs in freshwater lakes and coastal marine sites. Relative levels of gene flow were inferred in 12 lake and 14 Baltic Sea populations in northern Germany using nine highly polymorphic microsatellite markers developed for P. pectinatus. We found highly significant isolation-by-distance in both habitat types (P < 0.001). Genetic differentiation increased approximately 2.5-times faster among freshwater populations compared with those from the Baltic Sea. As different levels of genetic drift or population history cannot explain these differences, higher population connectivity in the sea relative to freshwater populations is the most likely source of contrasting evolutionary divergence. These findings are consistent with the notion that freshwater angiosperms are more conducive to allopatric speciation than their life-history counterparts in the sea, the relative species poor seagrasses. Surprisingly, population pairs from different habitat types revealed almost maximal genetic divergence expected for complete reproductive isolation, regardless of their respective geographical distance. Hence, the barrier to gene flow between lake and sea habitat types cannot be due to dispersal limitation. We may thus have identified a case of rapid incipient speciation in post-glacial populations of a widespread aquatic plant.

Population differentiation and hybrid success in Campanula americana: geography and genome size

Populations within a species may diverge through genetic drift and natural selection. Few studies report on population differentiation in autopolyploids where multiple gene copies and the ratio of cytoplasmic to nuclear genes differ from diploids and may influence divergence. In autotetraploid Campanula americana we created hybrids between populations that differed in geographic proximity and genome size. Differences in genome size (up to 6.5%) did not influence hybrid performance. In contrast, hybrid performance was strongly influenced by population proximity. F1 hybrids between distant populations performed poorly relative to their parents while hybrids between proximate populations outperformed their parents. Outbreeding depression was strongest for juvenile traits. The expression of outbreeding depression often differed between reciprocal hybrids indicating interactions between nuclear and cytoplasmic genes contribute to population differentiation. Because plants were grown under greenhouse conditions, the outbreeding depression was likely due to genetic (underdominance or loss of additive-by-additive epistasis) rather than ecological factors.

Isolation by distance and post-glacial range expansion in the rough-skinned newt, Taricha granulosa

Allozymes and mitochondrial DNA sequences were used to examine the phylogeographical history of the rough-skinned newt, Taricha granulosa, in western North America. Nineteen populations were analysed for allozyme variation at 45 loci, and 23 populations were analysed for cytochrome b sequence variation. Both data sets agree that populations in the southern part of the range are characterized by isolation by distance, whereas northern populations fit the expectations of a recent range expansion. However, the northern limit of isolation by distance (and the southern limit of range expansion) is located in Oregon State by the mtDNA data, and in Washington State by the allozyme data. Nevertheless, both data sets are consistent with the known Pleistocene history of western North America, with phylogenetically basal populations in central and northern California, and a recent range expansion in the north following the retreat of the Cordilleran ice sheet 10,000 years ago. Additionally, a population in Idaho, previously considered introduced from central California based on morphometric analyses, possesses a distinct mtDNA haplotype, suggesting it could be native. The relevance of these results for Pacific Northwest biogeography is discussed.

Spatial genetic structure among and within populations of Primula sieboldii growing beside separate streams

We investigated the hierarchical genetic structure of SSR (simple sequence repeats) and cpDNA (chloroplast DNA) polymorphisms among and within populations of Primula sieboldii, a heterostylous clonal herb. Seven out of eight populations at the study site, located in a mountainous region of Nagano Prefecture, had each developed alongside a different stream, and the other occurred on a flat area 70 m from the nearest stream. The magnitude of genetic differentiation among streamside populations in maternally inherited cpDNA (Phi = 0.341) was much higher than that in biparentally inherited SSRs (Phi = 0.011). This result suggests that seed dispersal among streams was restricted, and pollen was the primary agent of gene flow among streamside populations. In contrast, genetic differentiation among subpopulations within streams were low at both markers (Phi = 0.053 for cpDNA, Phi = 0.025 for SSR). This low differentiation among subpopulations in cpDNA compared with that among streamside populations suggest that seed dispersal occur along the stream probably during flooding. This hypothesis was supported by the fact that in cpDNA haplotypes, no clear genetic structure was detected within the streamside population, while a significant genetic structure was found within 20 m in the nonstreamside population. Furthermore, within the streamside populations, two pairs of ramets with identical multilocus genotypes for eight SSR loci were distantly (> 50 m) distributed along the same streamside, suggesting dispersal of clonal propagule. Our study showed that the heterogeneity of the landscape can influence gene flow and hence spatial genetic structure.

Spatial and population genetic structure of microsatellites in white pine

We evaluated the population genetic structure of seven microsatellite loci for old growth and second growth populations of eastern white pine (Pinus strobus). From each population, located within Hartwick Pines State Park, Grayling, Michigan, USA, 120-122 contiguous trees were sampled for genetic analysis. Within each population, genetic diversity was high and inbreeding low. When comparing these populations, there is a significant, but small (less than 1%), genetic divergence between populations. Spatial distance between populations or timber harvest at the second growth site were reasonable explanations for the observed minor differences in allele frequencies between populations. Spatial autocorrelation analysis suggested that, for the old growth population, weak positive structuring at 15 m fits the isolation by distance model for a neighbourhood size of about 100 individuals. In comparison, genotypes were randomly distributed in the second growth population. Thus, logging may have decreased spatial structuring at the second growth site, suggesting that management practices may be used to alter natural spatial patterns. In addition, the amount of autocorrelation in the old growth population appears to be lower for some of the microsatellites, suggesting higher numbers of rare alleles and that higher mutation rates may have directly affected spatial statistics by reducing structure.

Opposing natural selection from herbivores and pathogens may maintain floral-color variation in Claytonia virginica (Portulacaceae)

The maintenance of floral-color variation within natural populations is enigmatic because directional selection through pollinator preferences combined with random genetic drift should lead to the rapid loss of such variation. Fluctuating, balancing, and negative frequency-dependent selection mediated through pollinators have been identified as factors that may contribute to the maintenance of floral-color variation, and recently it has been suggested that indirect responses to selection on correlated characters through agents of selection other than pollinators may substantially shape the evolution of floral traits. Here, I provide empirical support for this latter view in Claytonia virginica (Portulacaceae) through a multiseason field study, a pollen supplementation study, and an artificial herbivory experiment. These studies show that most individuals fall into one of four discrete color classes, and suggest pollinator-mediated selection for increased floral redness in concurrent years. Floral color is also an indirect target of opposing directional selection via herbivores and pathogens that fluctuates through time. Taken together, these data suggest a novel mechanism by which floral-color variation may be maintained, and illustrate the importance of an inclusive, pluralistic view of selection when investigating the evolution of complex phenotypes.

Molecular differentiation within and among island populations of the endemic plant Scalesia affinis (Asteraceae) from the Galápagos Islands

Molecular variance was estimated in seven populations of the endemic species Scalesia affinis within and among islands of the Galapagos. The analysis, based on 157 polymorphic AFLP markers, revealed a high differentiation among populations, of which most was partitioned among islands. In addition, the information content of AFLP markers was tested with sets of discriminant analyses based on different numbers of AFLP markers. This indicated that the markers were highly informative in discriminating the populations. Although one of four populations from the island Isabela was sampled from a volcano 100 km away from the remaining populations, this population resembled the others on Isabela. The partitioning of molecular variance (AFLP) resulted in two unities, one consisting of populations from Isabela and one of populations from Santa Cruz and Floreana. The differentiation in two chloroplast microsatellites was higher than for AFLP markers and equally partitioned among populations within islands as among islands. Thus, gene flow via fruits within islands is as limited as among islands. The lower differentiation within islands in the nuclear AFLP markers may thus indicate that gene flow within islands is mostly accounted for by pollen transfer. S. affinis is the only species in the genus that is not listed in 2000 IUCN Red List of Threatened Species. However, due to prominent grazing and land exploitation, some populations have recently been reduced markedly, which was reflected in lower diversity. As inbreeding depression is present in the species, the rapid bottlenecks are threats to the populations.

Habitat specialization and adaptive phenotypic divergence of anuran populations

We tested for adaptive population structure in the frog Rana temporaria by rearing tadpoles from 23 populations in a common garden experiment, with and without larval dragonfly predators. The goal was to compare tadpole phenotypes with the habitats of their source ponds. The choice of traits and habitat variables was guided by prior information about phenotypic function. There were large differences among populations in life history, behaviour, morphological shape, and the predator-induced plasticities in most of these. Body size and behaviour were correlated with predation risk in the source pond, in agreement with adaptive population divergence. Tadpoles from large sunny ponds were morphologically distinct from those inhabiting small woodland ponds, although here an adaptive explanation was unclear. There was no evidence that plasticity evolves in populations exposed to more variable environments. Much among-population variation in phenotype and plasticity was not associated with habitat, perhaps reflecting rapid changes in wetland habitats.

Landscape genetics of alpine-snowbed plants: comparisons along geographic and snowmelt gradients

The genetic structure of three snowbed-herb species (Peucedanum multivittatum, Veronica stelleri, and Gentiana nipponica) was analyzed using allozymes across nine populations arranged as a matrix of three snowmelt gradients x three geographic locations within 3 km in the Taisetsu Mountains, northern Japan. Phenologically asynchronous populations are packed within a local area in alpine snowbeds, because flowering season of alpine plants depends strongly on the timing of snowmelt. Moderate genetic differentiation was detected among local populations in every species (FST=0.03-0.07). There was a significant correlation between the geographic distance and genetic distance in the P. multivittatum populations, but not in the V. stelleri and G. nipponica populations. On the other hand, a significant correlation between the phenological distance caused by snowmelt timing and genetic distance was detected in the V. stelleri and G. nipponica populations, but not in the P. multivittatum populations. The snowmelt gradient or geographic separation influenced hierarchical genetic structure of these species moderately (FRT <0.04). Restriction of gene flow due to phenological separation and possible differential selection along the snowmelt gradient may produce genetic clines at microgeographic scale in these species.

Population genetic structure of annual and perennial populations of Zostera marina L. along the Pacific coast of Baja California and the Gulf of California

The Baja California peninsula represents a biogeographical boundary contributing to regional differentiation among populations of marine animals. We investigated the genetic characteristics of perennial and annual populations of the marine angiosperm, Zostera marina, along the Pacific coast of Baja California and in the Gulf of California, respectively. Populations of Z. marina from five coastal lagoons along the Pacific coast and four sites in the Gulf of California were studied using nine microsatellite loci. Analyses of variance revealed significant interregional differentiation, but no subregional differentiation. Significant spatial differentiation, assessed using theta(ST) values, was observed among all populations within the two regions. Z. marina populations along the Pacific coast are separated by more than 220 km and had the greatest theta(ST) (0.13-0.28) values, suggesting restricted gene flow. In contrast, lower but still significant genetic differentiation was observed among populations within the Gulf of California (theta(ST) = 0.04-0.18), even though populations are separated by more than 250 km. This suggests higher levels of gene flow among Gulf of California populations relative to Pacific coast populations. Direction of gene flow was predominantly southward among Pacific coast populations, whereas no dominant polarity in the Gulf of California populations was observed. The test for isolation by distance (IBD) showed a significant correlation between genetic and geographical distances in Gulf of California populations, but not in Pacific coast populations, perhaps because of shifts in currents during El Niño Southern Oscillation (ENSO) events along the Pacific coast.

Phenotypic divergence despite high levels of gene flow in Galápagos lava lizards (Microlophus albemarlensis)

The extent of evolutionary divergence of phenotypes between habitats is predominantly the result of the balance of differential natural selection and gene flow. Lava lizards (Microlophus albemarlensis) on the small island of Plaza Sur in the Galápagos archipelago inhabit contrasting habitats: dense vegetation on the western end of the island thins rapidly in a transitional area, before becoming absent on the eastern half. Associated with these habitats are phenotypic differences in traits linked to predator avoidance (increased wariness, sprint speed, and endurance in lizards from the sparsely vegetated habitat). This population provides an opportunity to test the hypothesis that reduced gene flow is necessary for phenotypic differentiation. There was no evidence of any differences among habitats in allele frequencies at six out of seven microsatellite loci examined, nor was there any indication of congruence between patterns of genetic variability and the change in vegetation regime. We infer that gene flow between the habitats on Plaza Sur must be sufficiently high to overcome genetic drift within habitats but that it does not preclude phenotypic differentiation.

Mutation at high rates reduces spatial structure within populations

The question of whether or not the high rates (mu) of mutation that occur for some hypervariable markers can affect commonly used empirical measures of spatial structure of genetic variation within populations is addressed. The results show that values of these measures are approximately halved when mu is 10(-2). Finest spatial-scale correlations, measured by either Moran's I-statistics or conditional kinship, are reduced by 30%-50%. When the mutation rate is 10 times lower, much smaller reductions result, e.g. averaging 7% for the finest scale correlations. Still smaller orders of magnitude of mu cause negligible changes in spatial structure, where any effects normally would not be detectable. The reductions are caused by forward mutations, and when the reductions are measured as percentages, they are nearly independent of the amount of structure produced sans mutation, except when dispersal is nearly minimal. The percent reductions are also nearly independent of the number of alleles and of back mutations, hence of the nature of the mutation process (e.g. stepwise or not). The results demonstrate that some hypervariable loci should have reduced spatial structuring, and that marker choice may affect the values observed in experimental surveys. Moreover, if fine-scale correlations are used to indirectly estimate dispersal distances, then mutation at high rates could inflate estimates, easily up to two- to three-fold.

Impacts of gene flow and logging history on the local genetic structure of a scattered tree species, Sorbus torminalis L. Crantz

Sorbus torminalis L. Crantz is a colonizing tree species usually found at low density in managed European forests. Using six microsatellite markers, we investigated spatial and temporal patterns of genetic structure within a 472-ha population of 185 individuals to infer processes shaping the distribution of genetic diversity. Only eight young stems were found to be the result of vegetative reproduction. Despite high levels of gene flow (standard deviation of gene dispersal = 360 m), marked patterns of isolation by distance were detected, associated with an aggregated distribution of individuals in approximately 100-m patches. This spatial structure of both genes and individuals is likely to result from patterns of seedling recruitment combined with low tree density. Our results suggest that landscape factors and logging cycles markedly shape the distribution of favourable sites for seedling establishment, which are then colonized by sibling cohorts as a result of joint seed transportation by frugivores. These combined genetic and demographic processes result in similar genetic structure both within and among logging units. However, conversion to high forest may enhance genetic structuring.

Isolation by distance and sharp discontinuities in gene frequencies: implications for the phylogeography of an alpine insect species, Carabus solieri

Analysis of genetic isolation by distance (IBD) is of prime importance for the study of processes responsible for spatial population genetic structure and is thus frequently used in case studies. However, the identification of a significant IBD pattern does not necessarily imply the absence of sharp discontinuities in gene frequencies. Therefore, identifying barriers to gene flow and/or secondary contact between differentiated entities remains a major challenge in population biology. Geographical genetic structure of 41 populations (1080 individuals) of an alpine insect species, Carabus solieri, was studied using 10 microsatellite loci. All populations were significantly differentiated and spatially structured according to IBD over the entire range. However, clustering analyses clearly identified three main clusters of populations, which correspond to geographical entities. Whereas IBD also occurs within each cluster, population structure was different according to which group of populations was considered. The southernmost cluster corresponds to the most fragmented part of the range. Consistently, it was characterized by relatively high levels of differentiation associated with low genetic diversity, and the slope of the regression of genetic differentiation against geographical distances was threefold those of the two other clusters. Comparisons of within-cluster and between-cluster IBD patterns revealed barriers to gene flow. A comparison of the two approaches, IBD and clustering analyses, provided us with valuable information with which to infer the phylogeography of the species, and in particular to propose postglacial colonization routes from two potential refugia located in Italy and in southeastern France. Our study highlights strongly the possible confounding contribution of barriers to gene flow to IBD pattern and emphasizes the utility of the model-based clustering analysis to identify such barriers.

Microsatellite Analysis of the Spectacled Bear (Tremarctos ornatus) Across its Range Distribution

DNA samples of the spectacled bear (Tremarctos ornatus) from five Andean countries, Venezuela, Colombia, Ecuador, Peru and Bolivia, were analyzed for nine microsatellite loci. Seven of them were polymorphic, which led us to investigate several population-genetic parameters. Private alleles and significant differences in gene frequencies were found among the populations studied, which demonstrated the extent of genetic differentiation among the spectacled bear populations. The levels of gene diversity measured with these microsatellites were rather modest in this species. Hardy-Weinberg disequilibrium was especially found for the overall and the Ecuadorian samples, and might be due to the Wahl-und effect or consanguinity. Significant genetic heterogeneity was mainly observed among the Colombian and the Ecuadorian populations. Markov chain Monte Carlo simulations clearly showed that two different gene pools were present, one present in the Venezuelan-Colombian bears and other in the Ecuadorian ones.