Serial founder effects and genetic differentiation during worldwide range expansion of monarch butterflies

Range expansions can result in founder effects, increasing genetic differentiation between expanding populations and reducing genetic diversity along the expansion front. However, few studies have addressed these effects in long-distance migratory species, for which high dispersal ability might counter the effects of genetic drift. Monarchs (Danaus plexippus) are best known for undertaking a long-distance annual migration in North America, but have also dispersed around the world to form populations that do not migrate or travel only short distances. Here, we used microsatellite markers to assess genetic differentiation among 18 monarch populations and to determine worldwide colonization routes. Our results indicate that North American monarch populations connected by land show limited differentiation, probably because of the monarch''s ability to migrate long distances. Conversely, we found high genetic differentiation between populations separated by large bodies of water. Moreover, we show evidence for serial founder effects across the Pacific, suggesting stepwise dispersal from a North American origin. These findings demonstrate that genetic drift played a major role in shaping allele frequencies and created genetic differentiation among newly formed populations. Thus, range expansion can give rise to genetic differentiation and declines in genetic diversity, even in highly mobile species.     

Evidence of rapid change in genetic structure and diversity during range expansion in a recovering large terrestrial carnivore

Recovery of natural populations occurs often with simultaneous or subsequent range expansions. According to population genetic theory, genetic structuring emerges at the expansion front together with decreasing genetic diversity, owing to multiple founder events. Thereupon, as the expansion proceeds and connectivity among populations is established, homogenization and a resurgence of genetic diversity are to be expected. Few studies have used a fine temporal scale combined with genetic sampling to track range expansions as they proceed in wild animal populations. As a natural experiment, the historical eradication of large terrestrial carnivores followed by their recovery and recolonization may facilitate empirical tests of these ideas. Here, using brown bear (Ursus arctos) as model species, we tested predictions from genetic theory of range expansion. Individuals from all over Finland were genotyped for every year between 1996 and 2010 using 12 validated autosomal microsatellite markers. A latitudinal shift of about 110 km was observed in the distribution and delineation of genetic clusters during this period. As the range expansion proceeded, we found, as theory predicts, that the degree of genetic structure decreased, and that both genetic variation and admixture increased. The genetic consequences of range expansions may first be detected after multiple generations, but we found major changes in genetic composition after just 1.5 generations, accompanied by population growth and increased migration. These rapid genetic changes suggest an ongoing concerted action of geographical and demographic expansion combined with substantial immigration of bears from Russia during the recovery of brown bears within the large ecosystem of northern Europe.     

Shifting dispersal modes at an expanding species’ range margin

While it is generally recognized that noncontiguous (long‐distance) dispersal of small numbers of individuals is important for range expansion over large geographic areas, it is often assumed that colonization on more local scales proceeds by population expansion and diffusion dispersal (larger numbers of individuals colonizing adjacent sites). There are few empirical studies of dispersal modes at the front of expanding ranges, and very little information is available on dispersal dynamics at smaller geographic scales where we expect contiguous (diffusion) dispersal to be prevalent. We used highly polymorphic genetic markers to characterize dispersal modes at a local geographic scale for populations at the edge of the range of a newly invasive grass species (Brachypodium sylvaticum) that is undergoing rapid range expansion in the Pacific Northwest of North America. Comparisons of Bayesian clustering of populations, patterns of genetic diversity, and gametic disequilibrium indicate that new populations are colonized ahead of the invasion front by noncontiguous dispersal from source populations, with admixture occurring as populations age. This pattern of noncontiguous colonization was maintained even at a local scale. Absence of evidence for dispersal among adjacent pioneer sites at the edge of the expanding range of this species suggests that pioneer populations undergo an establishment phase during which they do not contribute emigrants for colonization of neighbouring sites. Our data indicate that dispersal modes change as the invasion matures: initial colonization processes appear to be dominated by noncontiguous dispersal from only a few sources, while contiguous dispersal may play a greater role once populations become established.     

Genetic signature of the colonisation dynamics along a coastal expansion front in the damselfly Coenagrion scitulum

1. Many insects are expanding their distribution range polewards as a result of climate change, which has been shown to be associated with founder effects leading to a reduction in genetic diversity and an increase in genetic differentiation. These spatial genetic patterns may arise from colonisation from a broad expansion front or a limited neighbourhood after a stepping stone model of dispersal. The temporal persistence of such founder effects are poorly understood, mainly because studies looking at the fine‐scale initial temporal dynamics of the genetic signature of a range expansion are rare. 2. Using microsatellite markers, we performed a detailed spatiotemporal genetic analysis of the range expanding damselfly Coenagrion scitulum (Rambur) along a coastal axis during the first years after colonisation. 3. A decrease was in (private) allelic richness when going northwards along the coastline, which is consistent with a scenario of cumulative founder events. In spite of the spatiotemporal dynamics in the observation records of the species along the coastline, the spatial genetic data indicated a major contribution from the broad expansion front during the colonisation of the coastline rather than a stepping‐stone colonisation process. 4. The fine‐scale temporal dynamics of the range expansion indicated the absence of persistent founder effects and instead showed considerable temporal instability in genetic indices at the more northern edge populations. This may be explained by genetic immigration and admixture from the broad expansion front in this active disperser.     

Patterns of genetic diversity reveal multiple introductions and recurrent founder effects during range expansion in invasive populations of Geranium carolinianum (Geraniaceae)

Genetic diversity, and thus the adaptive potential of invasive populations, is largely based on three factors: patterns of genetic diversity in the species'' native range, the number and location of introductions and the number of founding individuals per introduction. Specifically, reductions in genetic diversity (‘founder effects'') should be stronger for species with low within-population diversity in their native range and few introductions of few individuals to the invasive range. We test these predictions with Geranium carolinianum, a winter annual herb native to North America and invasive in China. We measure the extent of founder effects using allozymes and microsatellites, and ask whether this is consistent with its colonization history and patterns of diversity in the native range. In the native range, genetic diversity is higher and structure is lower than expected based on life history traits. In China, our results provide evidence for multiple introductions near Nanjing, Jiangsu province, with subsequent range expansion to the west and south. Patterns of genetic diversity across China reveal weak founder effects that are driven largely by low-diversity populations at the expansion front, away from the introduction location. This suggests that reduced diversity in China has resulted from successive founder events during range expansion, and that the loss of genetic diversity in the Nanjing area was mitigated by multiple introductions from diverse source populations. This has implications for the future of G. carolinianum in China, as continued gene flow among populations should eventually increase genetic diversity within the more recently founded populations.     

Fine-scale genetic structure and marginal processes in an expanding population of Biscutella laevigata L. (Brassicaceae)

Evolutionary processes acting at the expanding margins of a species' range are still poorly understood. Genetic drift is considered prevalent in marginal populations, and the maintenance of genetic diversity during recolonization might seem puzzling. To investigate such processes, a fine-scale investigation of 219 individuals was performed within a population of Biscutella laevigata (Brassicaceae), located at the leading edge of its range. The survey used amplified fragment length polymorphisms (AFLPs). As commonly reported across the whole species distribution range, individual density and genetic diversity decreased along the local axis of recolonization of this expanding population, highlighting the enduring effect of the historical colonization on present-day diversity. The self-incompatibility system of the plant may have prevented local inbreeding in newly found patches and sustained genetic diversity by ensuring gene flow from established populations. Within the more continuously populated region, spatial analysis of genetic structure revealed restricted gene flow among individuals. The distribution of genotypes formed a mosaic of relatively homogenous patches within the continuous population. This pattern could be explained by a history of expansion by long-distance dispersal followed by fine-scale diffusion (that is, a stratified dispersal combination). The secondary contact among expanding patches apparently led to admixture among differentiated genotypes where they met (that is, a reshuffling effect). This type of dynamics could explain the maintenance of genetic diversity during recolonization.     

Weak founder effects but significant spatial genetic imprint of recent contraction and expansion of European beech populations

Understanding the ecological and evolutionary processes occurring during species range shifts is important in the current context of global change. Here, we investigate the interplay between recent expansion, gene flow and genetic drift, and their consequences for genetic diversity and structure at landscape and local scales in European beech (Fagus sylvatica L.) On Mont Ventoux, South-Eastern France, we located beech forest refugia at the time of the most recent population minimum, ~150 years ago, and sampled 71 populations (2042 trees) in both refugia and expanding populations over an area of 15,000 ha. We inferred patterns of gene flow and genetic structure using 12 microsatellite markers. We identified six plots as originating from planting, rather than natural establishment, mostly from local genetic material. Comparing genetic diversity and structure in refugia versus recent populations did not support the existence of founder effects: heterozygosity (He = 0.667) and allelic richness (Ar = 4.298) were similar, and F_ST was low (0.031 overall). Still, significant spatial evidence of colonization was detected, with He increasing along the expansion front, while genetic differentiation from the entire pool (β_WT) decreased. Isolation by distance was found in refugia but not in recently expanding populations. Our study indicates that beech capacities for colonization and gene flow were sufficient to preserve genetic diversity despite recent forest contraction and expansion. Because beech has long distance pollen and seed dispersal, these results illustrate a ‘best case scenario’ for the maintenance of high genetic diversity and adaptive potential under climate-change-related range change.Subject terms: Molecular ecology, Forest ecology, Genetic variation     

Some genetic consequences of ice ages, and their role, in divergence and speciation

The genetic effects of pleistocene ice ages are approached by deduction from paleoenvironmental information, by induction from the genetic structure of populations and species, and by their combination to infer likely consequences. (1) Recent palaeoclimatic information indicate rapid global reversals and changes in ranges of species which would involve elimination with spreading from the edge. Leading edge colonization during a rapid expansion would be leptokurtic and lead to homozygosity and spatial assortment of genomes. In Europe and North America, ice age contractions were into southern refugia, which would promote genome reorganization. (2) The present day genetic structure of species shows frequent geographic subdivision, with parapatric genomes, hybrid zones and suture zones. A survey of recent DNA phylogeographic information supports and extends earlier work. (3) The grasshopper Chorthippus parallelus is used to illustrate such data and processes. Its range in Europe is divided on DNA sequences into five parapatric races, with southern genomes showing greater haplotype diversity - probably due to southern mountain blocks acting as refugia and northern expansion reducing diversity. (4) Comparison with other recent studies shows a concordance of such phylogeographic data over pleistocene time scales. (5) The role that ice age range changes may have played in changing adaptations is explored, including the limits of range, rapid change in new invasions and refugial differentiation in a variety of organisms. (6) The effects of these events in causing divergence and speciation are explored using Chorthippus as a paradigm. Repeated contraction and expansion would accumulate genome differences and adaptations, protected from mixing by hybrid zones, and such a composite mode of speciation could apply to many organisms.     

Fine-scale genetic structure arises during range expansion of an invasive gecko

Processes of range expansion are increasingly important in light of current concerns about invasive species and range shifts due to climate change. Theoretical studies suggest that genetic structuring may occur during range expansion. Ephemeral genetic structure can have important evolutionary implications, such as propagating genetic changes along the wave front of expansion, yet few studies have shown evidence of such structure. We tested the hypothesis that genetic structure arises during range expansion in Hemidactylus mabouia, a nocturnal African gecko recently introduced to Florida, USA. Twelve highly variable microsatellite loci were used to screen 418 individuals collected from 43 locations from four sampling sites across Florida, representing a gradient from earlier (～1990s) to very recent colonization. We found earlier colonized locations had little detectable genetic structure and higher allelic richness than more recently colonized locations. Genetic structuring was pronounced among locations at spatial scales of tens to hundreds of meters near the leading edge of range expansion. Despite the rapid pace of range expansion in this introduced gecko, dispersal is limited among many suitable habitat patches. Fine-scale genetic structure is likely the result of founder effects during colonization of suitable habitat patches. It may be obscured over time and by scale-dependent modes of dispersal. Further studies are needed to determine if such genetic structure affects adaptation and trait evolution in range expansions and range shifts.     

Population genetic structure and demographic history of the lone star tick,Amblyomma americanum (Ixodida: Ixodidae): New evidence supporting old records

Range expansions are a potential outcome of changes in habitat suitability, which commonly result as a consequence of climate change. Hypotheses on such changes in the geographic distribution of a certain species can be evaluated using population genetic structure and demography. In this study we explore the population genetic structure, genetic variability, demographic history of, and habitat suitability for Amblyomma americanum, a North American tick species that is a known vector of several pathogenic microorganisms. We used a double digestion restriction site‐associated DNA sequencing technique (dd‐RAD seq) and discovered 8,181 independent single nucleotide polymorphisms (SNPs) in 189 ticks from across the geographic range of the species. Genetic diversity was low, particularly when considering the broad geographic range of this species. The edge populations were less diverse than populations belonging to the historic range, possibly indicative of a range expansion, but this hypothesis was not statistically supported by a test based on genetic data. Nonetheless, moderate levels of population structure and substructure were detected between geographic regions. For New England, demographic and species distribution models support a scenario where A. americanum was present in more northern locations in the past, underwent a bottleneck, and subsequently recovered. These results are consistent with a hypothesis that this species is re‐establishing in this area, rather than one focused on range expansion from the south. This hypothesis is consistent with old records describing the presence of A. americanum in the northeastern US in the early colonial period.     

Genetic drift in range expansions is very sensitive to density dependence in dispersal and growth

Theory predicts rapid genetic drift during invasions, yet many expanding populations maintain high genetic diversity. We find that genetic drift is dramatically suppressed when dispersal rates increase with the population density because many more migrants from the diverse, high‐density regions arrive at the expansion edge. When density dependence is weak or negative, the effective population size of the front scales only logarithmically with the carrying capacity. The dependence, however, switches to a sublinear power law and then to a linear increase as the density dependence becomes strongly positive. We develop a unified framework revealing that the transitions between different regimes of diversity loss are controlled by a single, universal quantity: the ratio of the expansion velocity to the geometric mean of dispersal and growth rates at expansion edge. Our results suggest that positive density dependence could dramatically alter evolution in expanding populations even when its contribution to the expansion velocity is small.     

Genetic drift and rapid evolution of viviparity in insular fire salamanders (Salamandra salamandra)

Continental islands offer an excellent opportunity to investigate adaptive processes and to time microevolutionary changes that precede macroevolutionary events. We performed a population genetic study of the fire salamander (Salamandra salamandra), a species that displays unique intraspecific diversity of reproductive strategies, to address the microevolutionary processes leading to phenotypic and genetic differentiation of island, coastal and interior populations. We used eight microsatellite markers to estimate genetic diversity, population structure and demographic parameters in viviparous insular populations and ovoviviparous coastal and interior populations. Our results show considerable genetic differentiation (F(ST) range: 0.06-0.27), and no clear signs of gene flow among populations, except between the large and admixed interior populations. We find no support for island colonization by rafting or intentional/accidental anthropogenic introductions, indicating that rising sea levels were responsible for isolation of the island populations approximately 9000 years ago. Our study provides evidence of rapid genetic differentiation between island and coastal populations, and rapid evolution of viviparity driven by climatic selective pressures on island populations, geographic isolation with genetic drift, or a combination of these factors. Studies of these viviparous island populations in early stages of divergence help us better understand the microevolutionary processes involved in rapid phenotypic shifts.     

Inferring the mode of colonization of the rapid range expansion of a solitary bee from multilocus DNA sequence variation

Rapid geographic range expansions can have dramatic effects on the distribution of genetic diversity, both within and among populations. Based on field records collected over the past two decades in Western Europe, we report on the rapid geographic range expansion in Colletes hederae, a solitary bee species. To characterize how this expansion shaped the distribution of genetic diversity within and among populations, we performed a genetic analysis based on the sequencing of three nuclear loci (RNAp, CAD and WgL). We then simulated the evolution of DNA sequences under a spatially explicit model of coalescence to compare different hypotheses regarding the mode of colonization associated with this rapid expansion and to identify those that are most consistent with the observed molecular data. Our genetic analyses indicate that the range expansion was not associated with an important reduction in genetic diversity, even in the most recently colonized area in the United Kingdom. Moreover, little genetic differentiation was observed among populations. Our comparative analysis of simulated data sets indicates that the observed genetic data are more consistent with a demographic scenario involving relatively high migration rates than with a scenario based on a high reproduction rate associated with few migrants. In the light of these results, we discuss the factors that might have contributed to the rapid geographic range expansion of this pollen‐specialist solitary bee species across Western Europe.     

Neutral and adaptive genomic signatures of rapid poleward range expansion

Many species are expanding their range polewards, and this has been associated with rapid phenotypic change. Yet, it is unclear to what extent this reflects rapid genetic adaptation or neutral processes associated with range expansion, or selection linked to the new thermal conditions encountered. To disentangle these alternatives, we studied the genomic signature of range expansion in the damselfly Coenagrion scitulum using 4950 newly developed genomic SNPs and linked this to the rapidly evolved phenotypic differences between core and (newly established) edge populations. Most edge populations were genetically clearly differentiated from the core populations and all were differentiated from each other indicating independent range expansion events. In addition, evidence for genetic drift in the edge populations, and strong evidence for adaptive genetic variation in association with the range expansion was detected. We identified one SNP under consistent selection in four of the five edge populations and showed that the allele increasing in frequency is associated with increased flight performance. This indicates collateral, non‐neutral evolutionary changes in independent edge populations driven by the range expansion process. We also detected a genomic signature of adaptation to the newly encountered thermal regimes, reflecting a pattern of countergradient variation. The latter signature was identified at a single SNP as well as in a set of covarying SNPs using a polygenic multilocus approach to detect selection. Overall, this study highlights how a strategic geographic sampling design and the integration of genomic, phenotypic and environmental data can identify and disentangle the neutral and adaptive processes that are simultaneously operating during range expansions.     

Recovery process of genetic diversity through seed and pollen immigration at the northernmost leading-edge population of Fagus crenata

The range expansion of a plant species begins with colonization of ecological empty patches from posterior source populations. This process involves stochastic loss of genetic diversity. However, the founder population could restore genetic diversity by gene flow from posterior populations via seeds and pollen and its recovery affects evolutionary potential for species expansion. To clarify the recovery process of genetic diversity during species range expansion, gene flow via seeds and pollen was investigated at the expansion front of Fagus crenata. Based on eight nuclear microsatellite genotypes of a total of 150 individuals and 225 seeds at the northernmost leading-edge population, genetic diversity, fine-scale spatial genetic structure (FSGS), and genetic differentiation from other five northern populations were investigated. Moreover, both seed and pollen immigration and their effects on genetic diversity at different successional stages were analyzed. The leading-edge population showed lower genetic diversity and substantial genetic differentiation, reflecting its strong genetic drift. Non-significant FSGS and a negative inbreeding coefficient for mature trees may indicate that the earliest generation consisted of founders from foreign seed sources. The significant proportion of seed and pollen immigration increased the number of different alleles for later successional stages. The effective number of pollen parents from foreign sources (20.8) was markedly higher than that from the local source (2.1). These results indicated that pollen immigration incorporated new and rare alleles and increased the genetic diversity of the population. However, the proportion of foreign gene flow decreased during succession, probably due to the increased reproductive success of local individuals as they reached maturity and grew in size.     

Population genetic differences along a latitudinal cline between original and recently colonized habitat in a butterfly

BACKGROUND: Past and current range or spatial expansions have important consequences on population genetic structure. Habitat-use expansion, i.e. changing habitat associations, may also influence genetic population parameters, but has been less studied. Here we examined the genetic population structure of a Palaeartic woodland butterfly Pararge aegeria (Nymphalidae) which has recently colonized agricultural landscapes in NW-Europe. Butterflies from woodland and agricultural landscapes differ in several phenotypic traits (including morphology, behavior and life history). We investigated whether phenotypic divergence is accompanied by genetic divergence between populations of different landscapes along a 700 km latitudinal gradient. METHODOLOGY/PRINCIPAL FINDINGS: Populations (23) along the latitudinal gradient in both landscape types were analyzed using microsatellite and allozyme markers. A general decrease in genetic diversity with latitude was detected, likely due to post-glacial colonization effects. Contrary to expectations, agricultural landscapes were not less diverse and no significant bottlenecks were detected. Nonetheless, a genetic signature of recent colonization is reflected in the absence of clinal genetic differentiation within the agricultural landscape, significantly lower gene flow between agricultural populations (3.494) than between woodland populations (4.183), and significantly higher genetic differentiation between agricultural (0.050) than woodland (0.034) pairwise comparisons, likely due to multiple founder events. Globally, the genetic data suggest multiple long distance dispersal/colonization events and subsequent high intra- and inter-landscape gene flow in this species. Phosphoglucomutase deviated from other enzymes and microsatellite markers, and hence may be under selection along the latitudinal gradient but not between landscape types. Phenotypic divergence was greater than genetic divergence, indicating directional selection on some flight morphology traits. MAIN CONCLUSIONS/SIGNIFICANCE: Clinal differentiation characterizes the population structure within the original woodland habitat. Genetic signatures of recent habitat expansion remain, notwithstanding high gene flow. After differentiation through drift was excluded, both latitude and landscape were significant factors inducing spatially variable phenotypic variation.     

Long‐distance dispersal maximizes evolutionary potential during rapid geographic range expansion

Conventional wisdom predicts that sequential founder events will cause genetic diversity to erode in species with expanding geographic ranges, limiting evolutionary potential at the range margin. Here, we show that invasive European starlings (Sturnus vulgaris) in South Africa preserve genetic diversity during range expansion, possibly as a result of frequent long‐distance dispersal events. We further show that unfavourable environmental conditions trigger enhanced dispersal, as indicated by signatures of selection detected across the expanding range. This brings genetic variation to the expansion front, counterbalancing the cumulative effects of sequential founding events and optimizing standing genetic diversity and thus evolutionary potential at range margins during spread. Therefore, dispersal strategies should be highlighted as key determinants of the ecological and evolutionary performances of species in novel environments and in response to global environmental change.     

Strong genetic differentiation due to multiple founder events during a recent range expansion of an introduced wall lizard population

Biological invasions represent ideal systems for the study of evolutionary processes associated with colonization events. It has been hypothesized that the genetic diversity is generally decreasing from the centre of the range to the margins due to multiple founder events. Invasive populations offer the opportunity to test this hypothesis at a fine spatial and temporal scale. We analysed the genetic structure of a large expanding non-native population of the Common Wall Lizard (Podarcis muralis) in Passau (Germany) using thirteen microsatellite loci. We analyzed the genetic structure and levels of admixture across a transect reflecting the expansion process and tested for a loss of genetic diversity and an increase of genetic differentiation from the centre to the invasion front. Our results demonstrate that significant genetic population structure can emerge rapidly at a small spatial scale. We found a trend for an increase in genetic differentiation and a decrease in genetic diversity from the invasion centre to the expanding range margin, suggesting that genetic drift is the major factor causing this pattern. The correlation between genetic diversity and average genetic differentiation was significant among sites. We hypothesize that the territoriality of P. muralis generates sufficient rates of noncontiguous and stratified dispersal from longer established sites to maintain significant genetic diversity at the invasion front. Simultaneously, territoriality might restrict the colonization success of migrants at established sites, so that in combination with founder events a strong differentiation arises.     

Introgression in peripheral populations and colonization shape the genetic structure of the coastal shrub Armeria pungens

The coastal shrub Armeria pungens has a disjunct Atlantic-Mediterranean distribution. The historic range expansion underlying this distribution was investigated using the nuclear internal transcribed spacer region, three plastid regions (namely trnL-F, trnS-fM and matK) and morphometric data. A highly diverse ancestral lineage was identified in southwest Portugal. More recently, two areas have been colonized: (1) Corsica and Sardinia, where disjunct Mediterranean populations have been established as a result of the long-distance dispersal of Portuguese genotypes, and (2) the southern part of the Atlantic range, Gulf of Cadiz, where a distinct lineage showing no genetic differentiation among populations occurs. Genetic consequences of colonization seem to have been more severe in the Gulf of Cadiz than in Corsica-Sardinia. Although significant genetic divergence is associated with low plastid diversity in the Gulf of Cadiz, in Corsica-Sardinia, the loss of plastid haplotypes was not accompanied by divergence from disjunct Portuguese source populations. In addition, in its northernmost and southernmost populations, A. pungens exhibited evidence for ancient or ongoing introgression from sympatric congeners. Introgression might have created novel genotypes able to expand beyond the latitudinal margins of the species or, alternatively, these genotypes may be the result of surfing of alleles from other species in demographic equilibrium into peripheral populations of A. pungens. Our results highlight the evolutionary significance of genetic drift following the colonization of new areas and the key role of introgression in range expansion.     

Phylogeography and genetic structure of disjunct Salix arbutifolia populations in Japan

Disjunct geographic distributions of boreal plant species at the southern edges of their ranges are expected to result in low genetic diversity and high genetic differentiation in the disjunct populations. This prediction was tested in a riparian willow, Salix arbutifolia, distributed in the northeastern Eurasian continent and the Sakhalin, Hokkaido, and Honshu Islands, using chloroplast DNA haplotypes and nuclear microsatellite genotypes. Hokkaido and Honshu populations shared a chloroplast haplotype identical to a closely related species, S. cardiophylla. This haplotype was divergent from haplotypes in the Eurasian continent (Primorsky) and the Sakhalin Island. In the nuclear microsatellites, most Hokkaido populations were genetically closer to Primorsky populations than to Sakhalin populations in spite of the geographical vicinity between Sakhalin and Hokkaido. The unexpected genetic divergence between Sakhalin and Hokkaido implies a complicated history of migration and colonization. The most peripheral populations in Honshu had the lowest genetic diversity and were most differentiated from the others. Thus, low genetic diversity and high genetic differentiation at the range periphery suggest substantial effects of genetic drift on genetic structure in the disjunct populations of Salix arbutifolia at the southern edge of its range.