Contrasting patterns of historical colonization in white oaks (Quercus spp.) in California and Europe

Phylogeography allows the inference of evolutionary processes that have shaped the current distribution of genealogical lineages across a landscape. In this perspective, comparative phylogeographical analyses are useful in detecting common historical patterns by either comparing different species within the same area within a continent or by comparing similar species in different areas. Here, we analyse one taxon (the white oak, genus Quercus, subgenus Quercus, section Quercus) that is widespread worldwide, and we evaluate its phylogeographical pattern on two different continents: western North America and Western Europe. The goals of the present study are: (i) to compare the chloroplast genetic diversity found in one California oak species vs. that found in the extensively studied European oak species (in France and the Iberian Peninsula); (ii) to contrast the geographical structure of haplotypes between these two taxa and test for a phylogeographical structure for the California species. For this purpose, we used the same six maternally inherited chloroplast microsatellite markers and a similar sampling strategy. The haplotype diversity within site as well as the differentiation among sites was alike in both taxa, but the Californian species has higher allelic richness with a greater number of haplotypes (39 vs. 11 in the European white oak complex). Furthermore, in California these 39 haplotypes are distributed locally in patches while in the European oaks haplotypes are distributed into lineages partitioned longitudinally. These contrasted patterns could indicate that gene movement in California oak populations have been more stable in response to past climatic and geological events, in contrast to their European counterparts.     

Is Cupressus sempervirens native in Italy? An answer from genetic and palaeobotanical data

This study represents the first large-scale analysis using nuclear molecular markers to assess genetic diversity and structure of Cupressus sempervirens L.. Genetic and fossil data were combined to infer the possible role of human activity and evolutionary history in shaping the diversity of cypress populations. We analysed 30 populations with six polymorphic nuclear microsatellite markers. Dramatic reductions in heterozygosity and allelic richness were observed from east to west across the species range. Structure analysis assigned individuals to two main groups separating central Mediterranean and eastern populations. The two main groups could be further divided into five subgroups which showed the following geographical distributions: Turkey with the Greek islands Rhodes and Samos, Greece (Crete), Southern Italy, Northern Italy, Tunisia with Central Italy. This pattern of genetic structure is also supported by samova and Barrier analyses. Palaeobotanical data indicated that Cupressus was present in Italy in the Pliocene, Pleistocene and Holocene. Furthermore, our molecular survey showed that Italian cypress populations experienced bottlenecks that resulted in reduced genetic diversity and allelic richness and greater genetic differentiation. Recent colonization or introduction may also have influenced levels of diversity detected in the Italian populations, as most individuals found in this range today have multilocus genotypes that are also present in the eastern range of the species. The data reveal a new interpretation of the history of cypress distribution characterized by ancient eastern populations (Turkey and Greek islands) and a mosaic of recently introduced trees and remnants of ancient, depauperate populations in the central Mediterranean range.     

Patterns of seed dispersal and pollen flow in Quercus garryana (Fagaceae) following post-glacial climatic changes

Aim We examined the genetic structure of Quercus garryana to infer post‐glacial patterns of seed dispersal and pollen flow to test the hypotheses that (1) peripheral populations are genetically distinct from core populations and from one another; (2) genetic diversity declines towards the poleward edge of the species’ range; and (3) genetic diversity in the chloroplast genome, a direct measure of seed dispersal patterns, declines more sharply with increasing latitude than diversity in the nuclear genome. We address our findings in the context of known historical oak distribution from pollen core data derived from previously published research. Location Oak–savanna ecosystems from southern Oregon, USA (core populations/non‐glaciated range) northward to Vancouver Island, British Columbia, Canada (peripheral populations/glaciated range). Methods We genotyped 378 trees from 22 sites with five chloroplast and seven nuclear microsatellite loci. For both sets of markers, we estimated genetic diversity and differentiation using an analysis of molecular variance and generated Mantel correlograms to detect genetic and geographical distance correlations. For the nuclear markers, we also used a Bayesian approach to infer population substructure. Results There was a large degree of population differentiation revealed by six chloroplast haplotypes, with little (≤ 3) or no haplotype diversity within sites. Peripheral island locations shared the same, maternally inherited chloroplast haplotype, whereas locations in mainland Washington had greater haplotype diversity. In contrast, genetic diversity of the nuclear markers was high at all locations sampled. Populations clustered into two groups and were significantly positively correlated over large spatial scales (≤ 200 km), although allele richness decreased significantly with latitude. Population substructure was observed between core and peripheral populations because rare alleles were absent in peripheral localities and common allele frequencies differed. Main conclusions The observed pattern of chloroplast haplotype loss at the northern periphery suggests restricted seed dispersal events from mainland sites to peripheral islands. This pattern was unexpected, however, as refugial oak populations remained near the current post‐glacial range even during the Last Glacial Maximum. Using nuclear markers, we found high within‐population diversity and population differentiation only over large spatial scales, suggesting that pollen flow is relatively high among populations.     

Evidence for multiple sources of invasion and intraspecific hybridization in Brachypodium sylvaticum (Hudson) Beauv. in North America

We compared the levels and distribution of genetic diversity in Eurasian and North American populations of Brachypodium sylvaticum (Huds.) Beauv. (false brome), a newly invasive perennial bunchgrass in western North America. Our goals were to identify source regions for invasive populations, determine the number of independent invasion events, and assess the possibility that postinvasion bottlenecks and hybridization have affected patterns of genetic diversity in the invaded range. We tested the hypothesis that this Eurasian grass was accidentally introduced into two areas in Oregon and one site in California by examining nuclear microsatellites and chloroplast haplotype variation in 23 introduced and 25 native populations. In the invaded range, there was significantly lower allelic richness (R(S)), observed heterozygosity (H(O)) and within-population gene diversity (H(S)), although a formal test failed to detect a significant genetic bottleneck. Most of the genetic variation existed among populations in the native range but within populations in the invaded range. All of the allelic variation in the invaded range could be explained based on alleles found in western European populations. The distribution of identified genetic clusters in the North American populations and the unique alleles associated with them is consistent with two historical introductions in Oregon and a separate introduction to California. Further analyses of population structure indicate that intraspecific hybridization among genotypes from geographically distinct regions of western Europe occurred following colonization in Oregon. The California populations, however, are more likely to be derived from one or perhaps several genetically similar regions in the native range. The emergence and spread of novel recombinant genotypes may be facilitating the rapid spread of this invasive species in Oregon.     

Genetic diversity and differentiation of Michelia maudiae (Magnoliaceae) revealed by nuclear and chloroplast microsatellite markers

Michelia maudiae Dunn. is a Magnoliaceae species threatened by habitat destruction and over-exploitation. Genetic diversity and differentiation, population contribution to total diversity and allelic/haplotypic richness, and the relative importance of pollen- and seed-mediated gene flow were investigated in nine populations (192 individuals) of M. maudiae using nuclear and chloroplast microsatellites to further our understanding of the genetic structure and evolutionary history of this tree species and to provide a genetic perspective for its conservation. The species had strong pollen mediated gene flow in the past. The ratio of pollen to seed gene flow was 25.4. Three clusters from the western, central, and eastern China were identified by both chloroplast and nuclear microsatellites. Western populations at Xiaodanjiang and Daoxian were phylogenetically divergent from the remaining populations and might be particularly important for the conservation of this species. The populations of Xiaodanjiang, Daoxian, and Minjiangyuan made positive contribution to the total diversity and allelic/haplotypic richness, and were worthy of being conserved with priority. In the central cluster, population at Laopengding should be protected since it harbored the greatest genetic diversity.     

Tracking a genetic signal of extinction-recolonization events in a neotropical tree species: Vouacapoua americana Aublet in French Guiana

Drier periods from the late Pleistocene and early Holocene have been hypothesized to have caused the disappearance of various rainforest species over large geographical areas in South America and restricted the extant populations to mesic sites. Subsequent improvement in climatic conditions has been associated with recolonization. Changes in population size associated with these extinction-recolonization events should have affected genetic diversity within species. However, these historical hypotheses and their genetic consequences have rarely been tested in South America. Here, we examine the diversity of the chloroplast and nuclear genomes in a Neotropical rainforest tree species, Vouacapoua americana (Leguminosae, Caesalpinioideae) in French Guiana. The chloroplast diversity was analyzed using a polymerase chain reaction-restriction fragment length polymorphism method (six pairs of primers) in 29 populations distributed over most of French Guiana, and a subset of 17 populations was also analyzed at nine polymorphic microsatellite loci. To determine whether this species has experienced extinction-recolonization, we sampled populations in areas supposedly not or only slightly affected by climatic changes, where the populations would not have experienced frequent extinction, and in areas that appear to have been recently recolonized. In the putatively recolonized areas, we found patches of several thousands of hectares homogeneous for chloroplast variation that can be interpreted as the effect of recolonization processes from several geographical origins. In addition, we observed that, for both chloroplast and nuclear genomes, the populations in newly recolonized areas exhibited a significantly smaller allelic richness than others. Controlling for geographic distance, we also detected a significant correlation between chloroplast and nuclear population differentiation. This result indicates a cytonuclear disequilibrium that can be interpreted as a historical signal of a genetic divergence between fragmented populations. In conclusion, the spatial genetic structure of contemporary V. americana populations shows evidence that this species has experienced large extinction-recolonization events, which were possibly caused by past climatic change.     

Population differentiation of sessile oak at the altitudinal front of migration in the French Pyrenees

To assess the effects of altitude on the level and structure of genetic diversity, a genetic survey was conducted in 12 populations of sessile oak (Quercus petraea) located between 130 and 1660 m in two parallel valleys on the northern side of the Pyrenees Mountains. Genetic diversity was monitored at 16 nuclear microsatellite loci and 5 chloroplast DNA (cpDNA) markers. The cpDNA survey suggested that extant populations in both valleys shared the same source populations from the plain. There was no visible trend of nuclear genetic diversity along altitude, even if indirect estimates of effective population sizes revealed a consistent reduction at higher altitudes. Population differentiation, although low, was mostly present among populations of the same valleys and reached similar levels than differentiation across the range of distribution of sessile oak. Contribution to the overall differentiation in the valleys was mostly due to the genetic divergence of the highest populations and the altitudinal variation of allelic frequencies at a few loci. Bayesian inference of migration between groups of populations showed that gene flow is preferentially unidirectional from lower altitudes in one valley to other groups of populations. Finally, we found evidence of clonal reproduction in high altitude populations. The introgression of Quercus robur and Quercus pubescens was also more frequent at the altitudinal margin suggesting that this mechanism may have contributed to the present migration and adaptation of Q. petraea and may also facilitate its future upslope shift in the context of climate change.     

The importance of sexual and asexual reproduction in the recent evolution of Allium vineale

In the weedy plant species Allium vineale (wild garlic), individuals may simultaneously produce sexually and asexually derived offspring, by seed and bulbils, respectively. In this study, genetic and genotypic diversity was determined in samples from 14 European A. vineale populations using nuclear (RAPD) and cytoplasmic (PCR-RFLP of cpDNA) markers to investigate the importance of the different reproductive modes. In the whole sample, 77 nuclear multilocus genotypes and four chloroplast haplotypes (chlorotypes) were found. Populations exhibited a high degree of subdivision for nuclear and cytoplasmic markers as estimated from hierarchical F-statistics; at the same time, identical chlorotypes could be found in populations separated by large distances. Genotypic diversity was significantly lower than expected under free recombination in almost all populations, indicating that recruitment into populations is mostly by asexually produced offspring. Nevertheless, within each chlorotype, the distribution of markers from pairs of nuclear loci was incompatible with a purely clonal structure, suggesting that many multilocus genotypes have originated by sexual recombination rather than by mutation within asexual lineages. It is argued that the weedy habit of A. vineale is likely to have favored bulbil reproduction, whereas sexually generated genotypes may have facilitated local adaptation during the species' expansion across Europe.     

Geographical patterns of turnover and nestedness-resultant components of allelic diversity among populations

The analysis of geographical patterns in population divergence has always been a powerful way to infer microevolutionary processes involved in population differentiation, and several approaches have been used to investigate such patterns. Most frequently, multivariate spatial patterns of population differentiation are analyzed by computing pairwise genetic distances or F_ST (or related statistics, such as ?_ST from AMOVA), which are then correlated with geographical distances or landscape features. However, when calculating distances, especially based on presence-absence of alleles in local populations, there would be a confounding effect of allelic richness differences in the population differentiation. Moreover, the relative magnitude of these components and their spatial patterns can help identifying microevolutionary processes driving population differentiation. Here we show how recent methodological advances in ecological community analyses that allows partitioning dissimilarity into turnover (turnover) and richness differences, or nestedness-resultant dissimilarity, can be applied to allelic variation data, using an endemic Cerrado tree (Dipteryx alata) as a case study. Individuals from 15 local populations were genotyped for eight microsatellite loci, and pairwise dissimilarities were computed based on presence-absence of alleles. The turnover of alleles among populations represented 69?% of variation in dissimilarity, but only the richness difference component shows a clear spatial structure, appearing as a westward decrease of allelic richness. We show that decoupling richness difference and turnover components of allelic variation reveals more clearly how similarity among populations reflects geographical patterns in allelic diversity that can be interpreted in respect to historical range expansion in the species.     

Inferring the complex origins of horticultural invasives: French broom in California

Investigating the origins of invasive populations provides insight into the evolutionary and anthropogenic factors underlying invasions, and can inform management decisions. Invasive species introduced for horticultural purposes often have complex origins typified by multiple introductions of species, cultivars, and genotypes, and interspecific and intraspecific hybridizations in introduced ranges. Such complex introduction histories may result in complex genetic signatures in the invaded range, making inferences about origins difficult, particularly when all putative sources cannot be sampled. In this study, we inferred the origins of the invasive French broom complex in California using 12 nuclear microsatellite markers. We characterized the genetic diversity and population structure of invasive and horticultural brooms in their invaded range in California and of Genista monspessulana in its native Mediterranean range. Overall, no significant differences in allelic richness, observed heterozygosity, inbreeding, or genetic structure were observed between the invaded and native ranges, but differences existed among populations within ranges. Bayesian STRUCTURE analysis revealed three genetic clusters in the French broom complex. Nearly all native G. monspessulana assigned highly to a single cluster. Many invasives assigned to a second cluster that contained Genista canariensis, Genista stenopetala, and ornamental sweet broom, and the remaining invasives assigned to a third cluster that also contained some G. monspessulana individuals from Sardinia and Corsica. Admixture between the second and third clusters was detected. Approximate Bayesian Computation analysis of six alternative scenarios supported the hypothesis that some invasive French broom is derived from an unsampled population branching from ornamental sweet broom. A combination of factors, including multiple introductions, escapes from cultivation, and inter-taxon hybridization, likely contribute to the invasive success of French broom in California and may have important implications for management, in particular biological control.     

Mapping genetic diversity of cherimoya (Annona cherimola Mill.): application of spatial analysis for conservation and use of plant genetic resources

There is a growing call for inventories that evaluate geographic patterns in diversity of plant genetic resources maintained on farm and in species' natural populations in order to enhance their use and conservation. Such evaluations are relevant for useful tropical and subtropical tree species, as many of these species are still undomesticated, or in incipient stages of domestication and local populations can offer yet-unknown traits of high value to further domestication. For many outcrossing species, such as most trees, inbreeding depression can be an issue, and genetic diversity is important to sustain local production. Diversity is also crucial for species to adapt to environmental changes. This paper explores the possibilities of incorporating molecular marker data into Geographic Information Systems (GIS) to allow visualization and better understanding of spatial patterns of genetic diversity as a key input to optimize conservation and use of plant genetic resources, based on a case study of cherimoya (Annona cherimola Mill.), a Neotropical fruit tree species. We present spatial analyses to (1) improve the understanding of spatial distribution of genetic diversity of cherimoya natural stands and cultivated trees in Ecuador, Bolivia and Peru based on microsatellite molecular markers (SSRs); and (2) formulate optimal conservation strategies by revealing priority areas for in situ conservation, and identifying existing diversity gaps in ex situ collections. We found high levels of allelic richness, locally common alleles and expected heterozygosity in cherimoya's putative centre of origin, southern Ecuador and northern Peru, whereas levels of diversity in southern Peru and especially in Bolivia were significantly lower. The application of GIS on a large microsatellite dataset allows a more detailed prioritization of areas for in situ conservation and targeted collection across the Andean distribution range of cherimoya than previous studies could do, i.e. at province and department level in Ecuador and Peru, respectively.     

Across-species patterns of genetic variation in forest trees of Central Europe

The study focuses on geographical patterns of genetic variation at allozyme loci common for four main tree species of Central Europe (Norway spruce, silver fir, common beech and sessile oak). Moving-window averaging of four indicators of allelic richness and diversity (proportion of polymorphic loci, mean number of alleles per locus, effective number of alleles and expected heterozygosity) with window size of 50 × 50 km was used to identify the patterns. Moreover, local genetic divergence was assessed using the G _ST (Nei, Molecular population genetic and evolution, Amsterdam and Oxford, North-Holland, 1975) and D _j (Gregorius and Roberds, Theor Appl Genet 71:826–834, 1986) statistics for common beech and silver fir, where raw genotype data were available. Spatial patterns of diversity and allelic richness were quite similar. Romanian Carpathians were identified as the most important hotspot of genetic diversity and evolutionary divergence in Central Europe. Implications for genetic conservation are briefly discussed.     

The rise and fall of isolation by distance in the anadromous brook charr (Salvelinus fontinalis Mitchill)

Geographic patterns of genetic diversity depend on a species' demographic properties in a given habitat, which may change over time. The rates at which patterns of diversity respond to changes in demographic properties and approach equilibrium are therefore pivotal in our understanding of spatial patterns of diversity. The brook charr Salvelinus fontinalis is a coastal fish exhibiting limited marine movements, such that a stable one-dimensional isolation-by-distance (IBD) pattern should be observed over the whole range. Its range, however, recently shifted northward such that northern populations may still be in the process of reaching equilibrium. We investigated variation in IBD patterns, genetic divergence, and allelic richness at six microsatellite markers in 2087 anadromous brook charr from 59 rivers along the most likely postglacial colonization route. We observed a decrease in allelic richness, together with an increase in differentiation and a decrease in IBD in the most recently colonized northern populations, as expected following recent colonization. Contrary to expectation, however, similar patterns were also observed at the southernmost part of the range, despite the fact that these populations are not considered to be newly colonized. We propose that the loss of dispersal capabilities associated with anadromy may have caused the southernmost populations to evolve relatively independently of one another. This study thus demonstrated that changes in a species' geographic range and dispersal capabilities may contribute to shaping geographic patterns of genetic diversity.     

Genetically depauperate but widespread: the case of an emblematic Mediterranean pine

Genetic variation is generally considered a prerequisite for adaptation to new environmental conditions. Thus the discovery of genetically depauperate but geographically widespread species is unexpected. We used 12 paternally inherited chloroplast microsatellites to estimate population genetic variation across the full range of an emblematic circum-Mediterranean conifer, stone pine (Pinus pinea L.). The same chloroplast DNA haplotype is fixed in nearly all of the 34 investigated populations. Such a low level of variation is consistent with a previous report of very low levels of diversity at nuclear loci in this species. Stone pine appears to have passed through a severe and prolonged demographic bottleneck, followed by subsequent natural- and human-mediated dispersal across the Mediterranean Basin. No other abundant and widespread plant species has as little genetic diversity as P. pinea at both chloroplast and nuclear markers. However, the species harbors a nonnegligible amount of variation at adaptive traits. Thus a causal relationship between genetic diversity, as measured by marker loci, and the evolutionary precariousness of a species, cannot be taken for granted.     

High genetic differentiation of red gorgonian populations from the Atlantic Ocean and the Mediterranean Sea

Patterns of genetic variation within a species may be used to infer past events in the evolutionary history of marine species. In the present study we aimed to compare the genetic diversity of the red gorgonian Paramuricea clavata in the Atlantic Ocean and the Mediterranean Sea. For genetic markers we used microsatellites and a mitochondrial gene fragment. Our results revealed a distinct genetic composition and diversity between the Mediterranean and the Atlantic. The Mediterranean samples had higher microsatellite heterozygosity, allelic richness and private allelic richness. The hypotheses that can explain these patterns are the isolation of Atlantic populations and/or a founder effect. Additionally, a clear difference was obtained from the mitochondrial locus, since sequences from Atlantic and Mediterranean samples diverged by 1%, which is high for soft corals.     

Influence of late Quaternary climate change on present patterns of genetic variation in valley oak,Quercus lobata Née

Phylogeography and ecological niche models (ENMs) suggest that late Quaternary glacial cycles have played a prominent role in shaping present population genetic structure and diversity, but have not applied quantitative methods to dissect the relative contribution of past and present climate vs. other forces. We integrate multilocus phylogeography, climate‐based ENMs and multivariate statistical approaches to infer the effects of late Quaternary climate change on contemporary genetic variation of valley oak (Quercus lobata Née). ENMs indicated that valley oak maintained a stable distribution with local migration from the last interglacial period (~120 ka) to the Last Glacial Maximum (~21 ka, LGM) to the present compared with large‐scale range shifts for an eastern North American white oak (Quercus alba L.). Coast Range and Sierra Nevada foothill populations diverged in the late Pleistocene before the LGM [104 ka (28–1622)] and have occupied somewhat distinct climate niches, according to ENMs and coalescent analyses of divergence time. In accordance with neutral expectations for stable populations, nuclear microsatellite diversity positively correlated with niche stability from the LGM to present. Most strikingly, nuclear and chloroplast microsatellite variation significantly correlated with LGM climate, even after controlling for associations with geographic location and present climate using partial redundancy analyses. Variance partitioning showed that LGM climate uniquely explains a similar proportion of genetic variance as present climate (16% vs. 11–18%), and together, past and present climate explains more than geography (19%). Climate can influence local expansion–contraction dynamics, flowering phenology and thus gene flow, and/or impose selective pressures. These results highlight the lingering effect of past climate on genetic variation in species with stable distributions.     

Influence of environmental heterogeneity on genetic diversity and structure in an endemic southern Californian oak

Understanding how specific environmental factors shape gene flow while disentangling their importance relative to the effects of geographical isolation is a major question in evolutionary biology and a specific goal of landscape genetics. Here, we combine information from nuclear microsatellite markers and ecological niche modelling to study the association between climate and spatial genetic structure and variability in Engelmann oak (Quercus engelmannii), a wind-pollinated species with high potential for gene flow. We first test whether genetic diversity is associated with climatic niche suitability and stability since the Last Glacial Maximum (LGM). Second, we use causal modelling to analyse the potential influence of climatic factors (current and LGM niche suitability) and altitude in the observed patterns of genetic structure. We found that genetic diversity is negatively associated with local climatic stability since the LGM, which may be due to higher immigration rates in unstable patches during favourable climatic periods and/or temporally varying selection. Analyses of spatial genetic structure revealed the presence of three main genetic clusters, a pattern that is mainly driven by two highly differentiated populations located in the northern edge of the species distribution range. After controlling for geographic distance, causal modelling analyses showed that genetic relatedness decreases with the environmental divergence among sampling sites estimated as altitude and current and LGM niche suitability. Natural selection against nonlocal genotypes and/or asynchrony in reproductive phenology may explain this pattern. Overall, this study suggests that local environmental conditions can shape patterns of genetic structure and variability even in species with high potential for gene flow and relatively small distribution ranges.     

Species richness and evolutionary niche dynamics: a spatial pattern-oriented simulation experiment

Evolutionary processes underlying spatial patterns in species richness remain largely unexplored, and correlative studies lack the theoretical basis to explain these patterns in evolutionary terms. In this study, we develop a spatially explicit simulation model to evaluate, under a pattern-oriented modeling approach, whether evolutionary niche dynamics (the balance between niche conservatism and niche evolution processes) can provide a parsimonious explanation for patterns in species richness. We model the size, shape, and location of species' geographical ranges in a multivariate heterogeneous environmental landscape by simulating an evolutionary process in which environmental fluctuations create geographic range fragmentation, which, in turn, regulates speciation and extinction. We applied the model to the South American domain, adjusting parameters to maximize the correspondence between observed and predicted patterns in richness of about 3,000 bird species. Predicted spatial patterns, which closely resemble observed ones (r2=0.795), proved sensitive to niche dynamics processes. Our simulations allow evaluation of the roles of both evolutionary and ecological processes in explaining spatial patterns in species richness, revealing the enormous potential of the link between ecology and historical biogeography under integrated theoretical and methodological frameworks.     

Invasion of Lepidium draba (Brassicaceae) in the western United States: distributions and origins of chloroplast DNA haplotypes

Advances in phylogeography are of great value for understanding the population structure and origins of invasive genotypes. Such insights provide constructive information for current or future biological control research efforts. In this study, we investigated a highly variable chloroplast DNA (cpDNA) marker for populations of the weed Lepidium draba (Brassicaceae) in its native Eurasian and invasive US ranges. We sequenced DNA from 684 individuals from Eurasia and the US and found 41 different haplotypes. Our comparative study between the native and invasive ranges showed a 33% reduction in allelic richness (A) and a 7% reduction in haplotype diversity (h) since introduction into the US. Most genetic variation in the native range was observed within geographical regions and populations, not between regions, and this result was similar for the invasive range. Assignment tests indicated the most likely origins of many invasive haplotypes. Some of these occurred in western Europe, supporting an expanded native range that had been proposed for the species. Exact locations were identified for a diverse set of invasive haplotypes which can be used in ongoing host-specificity tests of potential biological control agents.