The hierarchical spatial distribution of chloroplast DNA polymorphism across the introduced range ofSilene vulgaris

Phylogeography of a widely distributed plant species reveals cryptic genetic lineages with parallel phenotypic responses to warming and drought conditions

To predict how widely distributed species will perform under future climate change, it is crucial to understand and reveal their underlying phylogenetics. However, detailed information about plant adaptation and its genetic basis and history remains scarce and especially widely distributed species receive little attention despite their putatively high adaptability. To examine the adaptation potential of a widely distributed species, we sampled the model plant Silene vulgaris across Europe. In a greenhouse experiment, we exposed the offspring of these populations to a climate change scenario for central Europe and revealed the population structure through whole-genome sequencing. Plants were grown under two temperatures (18°C and 21°C) and three precipitation regimes (65, 75, and 90 mm) to measure their response in biomass and fecundity-related traits. To reveal the population genetic structure, ddRAD sequencing was employed for a whole-genome approach. We found three major genetic clusters in S. vulgaris from Europe: one cluster comprising Southern European populations, one cluster of Western European populations, and another cluster containing central European populations. Population genetic diversity decreased with increasing latitude, and a Mantel test revealed significant correlations between F ST and geographic distances as well as between genetic and environmental distances. Our trait analysis showed that the genetic clusters significantly differed in biomass-related traits and in the days to flowering. However, half of the traits showed parallel response patterns to the experimental climate change scenario. Due to the differentiated but parallel response patterns, we assume that phenotypic plasticity plays an important role for the adaptation of the widely distributed species S. vulgaris and its intraspecific genetic lineages.

Genetic Diversity of Invasive Spartina alterniflora Loisel. (Poaceae) Introduced Unintentionally Into Japan and Its Invasion Pathway

Among invasive species, aquatic plants pose serious threats to local biodiversity and ecosystem functions. Spartina alterniflora Loisel. (Poaceae), native to the eastern United States, was introduced unintentionally into Japan (Aichi and Kumamoto Prefectures) at around 2010. This invasive species could easily and rapidly spread to estuarine areas of Japan via vigorous trade and transport, making the prediction of its future invasion necessary. Here, the distribution and structure of the genetic variation of S. alterniflora in Japan were examined using chloroplast DNA (cpDNA) and microsatellite genotyping analyses for clarifying its invasion route and process. According to the cpDNA analysis, S. alterniflora populations in Japan had a single haplotype (haplotype C4) that is the most dominant genotype around the Florida Peninsula, the region of its origin, and is also widely found in the introduced populations in the East Asia. Microsatellite analysis also showed a loss of genetic diversity in Japanese S. alterniflora populations (allelic richness (A _R) = 1.20-1.39) compared with that in its native region (A _R = 4.58-4.59), suggesting a founder effect on S. alterniflora that might have occurred after invasion of the species into Japan. The principal coordinate analysis and The STRUCTURE analysis indicated that no gene mixing among Japanese local populations (Aichi, northern and southern Kumamoto) was observed, indicating that Spartina invasion occurred independently into these regions. Among the three regions, trading between the ports of northern Kumamoto and the U.S. was obviously lower than trading with China. We concluded that invasive S. alterniflora might have independently invaded Japan at different times through an East Asia route, particularly via China (i.e., secondary introduction). Therefore, it is important to strengthen the quarantine control on the importation of commodities, especially of transport vehicles at potential donor spots (i.e., border control/border biosecurity system), and to share information networks on invasive species between each region/port for minimizing further risks of biological species such as Spartina.

Recurrent gene duplication in the angiosperm tribe Delphinieae (Ranunculaceae) inferred from intracellular gene transfer events and heteroplasmic mutations in the plastid matK gene

The study of intracellular gene transfer may allow for the detection of interesting evolutionary processes such as ancient polyploidization. We compared 24 plastid genomes (plastomes) from tribe Delphinieae, one from tribe Nigelleae and one from tribe Ranunculeae, including five newly sequenced genomes. The functional transfers of the plastids rpl32 and rps16 to the nucleus in tribe Delphinieae were identified. Unexpectedly, we discovered multiple divergent copies of the nuclear-encoded plastid rpl32 in the genus Aconitum. Phylogenetic and synonymous substitution rate analyses revealed that the nuclear-encoded plastid rpl32 underwent two major duplication events. These ancient gene duplication events probably occurred via multiple polyploidization events in Aconitum between 11.9 and 24.7 Mya. Furthermore, our sequence rate analysis indicated that the eight plastid-encoded rpl subunits in Aconitum had a significantly accelerated evolutionary rate compared to those in other genera, suggesting that highly divergent paralogs targeted to the plastid may contribute to an elevated rate of evolution in plastid rpl genes. In addition, heteroplasmy of the plastid matK from two Aconitum species suggested the existence of potentially functional plastid maturases in its plastome. Our results provide insight into the evolutionary history of the tribe Delphinieae.

Geographic distribution of genetic variation among native and introduced populations of Chinese tallow tree, Triadica sebifera (Euphorbiaceae)

PREMISE OF THE STUDY

Invasive plants often display genetically determined variation in patterns of growth and resource allocation between native and introduced genotypes, as well as among genotypes within different regions of the introduced range. We examined patterns of genetic variation within and among native and introduced populations of the tetraploid Chinese tallow tree (Triadica sebifera, Euphorbiaceae) to determine whether nonselective evolutionary processes or the introduction history could contribute to previously observed phenotypic differences between native and introduced populations as well as among introduced populations.

METHODS

We used six microsatellite markers to study 12 native populations in China, 51 introduced populations in the southeastern USA, and one introduced population in Australia.

KEY RESULTS

Genetic diversity was greater within and among native populations than introduced populations. Within the southeastern USA, populations in Georgia and South Carolina differed substantially in their genetic composition and had greater genetic diversity than the rest of the southeastern USA. Greater genetic similarity between some populations in the native range and introduced range indicate a common provenance for Georgia and South Carolina populations that could have come from any of several western or southern Chinese populations and a different provenance for other southeastern USA populations and the Australian population, which were most similar to more northeastern Chinese populations.

CONCLUSIONS

Differences among introduced populations in potentially adaptive traits (e.g., herbivore tolerance, herbivore resistance, growth rates) may result in part from the introduction history, in particular from differences present among source populations in the native range.

Multiple introductions boosted genetic diversity in the invasive range of black cherry (Prunus serotina; Rosaceae)

BACKGROUND AND AIMS

Black cherry (Prunus serotina) is a North American tree that is rapidly invading European forests. This species was introduced first as an ornamental plant, then it was massively planted by foresters in many countries, but its origins and the process of invasion remain poorly documented. Based on a genetic survey of both native and invasive ranges, the invasion history of black cherry was investigated by identifying putative source populations and then assessing the importance of multiple introductions on the maintenance of gene diversity.

METHODS

Genetic variability and structure of 23 populations from the invasive range and 22 populations from the native range were analysed using eight nuclear microsatellite loci and five chloroplast DNA regions.

KEY RESULTS

Chloroplast DNA diversity suggests there were multiple introductions from a single geographic region (the north-eastern United States). A low reduction of genetic diversity was observed in the invasive range for both nuclear and plastid genomes. High propagule pressure including both the size and number of introductions shaped the genetic structure in Europe and boosted genetic diversity. Populations from Denmark, The Netherlands, Belgium and Germany showed high genetic diversity and low differentiation among populations, supporting the hypothesis that numerous introduction events, including multiple individuals and exchanges between sites, have taken place during two centuries of plantation.

CONCLUSIONS

This study postulates that the invasive black cherry has originated from east of the Appalachian Mountains (mainly the Allegheny plateau) and its invasiveness in north-western Europe is mainly due to multiple introductions containing high numbers of individuals.

Recombination and linkage disequilibrium among mitochondrial genes in structured populations of the gynodioecious plant Silene vulgaris

The impact of intergenic recombination on the population genetics of plant mitochondrial genomes is unknown. In an effort to study this in the gynodioecious plant Silene vulgaris three-locus PCR/RFLP genotypes (based on the mitochondrial genes atpA, cox1, and cob) were determined for 239 individuals collected from 20 North American populations. Seventeen three-locus PCR/RFLP genotypes were found. Recombination was indicated by observation of each of the four two-locus genotypes possible when the two most common alleles are considered for each of two loci. Based on these common alleles the absolute values of standardized linkage disequilibrium |D'| between pairs of loci range from 0.17 to 0.78. This indicates modest disequilibrium, rather than the maximum value expected in the absence of recombination |D'=1|, or the linkage equilibrium expected if recombination is pervasive (D'=0). Values of D' did not depend on which pair of loci contributed alleles to the analysis. The direction of D' obtained for the common atpA and cox1 alleles was comparable in sign and magnitude to that obtained by examining similar information obtained in a prior study of European samples. All three loci indicated a high degree of population structure (average FST=0.63), which would limit the within-population genetic diversity required for intergenic recombination to create novel genotypes, if most mating is local. Thus, population structure acts as a constraint on the approach to linkage equilibrium.

Spatial analysis of nuclear and cytoplasmic DNA diversity in wild sea beet (Beta vulgaris ssp. maritima) populations: do marine currents shape the genetic structure?

Patterns of seed dispersal in the wild sea beet (Beta vulgaris ssp. maritima) are predicted to be influenced by marine currents because populations are widely distributed along the European Atlantic coast. We investigated the potential influence of marine currents on the pattern of spatial genetic structuring in natural populations of sea beet. Populations were located along the French coasts of the Anglo-Norman gulf that features peculiar marine currents in the Channel. Thirty-three populations were sampled, among which 23 were continental and 10 were insular populations located in Jersey, Guernsey and Chausey, for a total of 1224 plants genotyped. To validate the coastal topography influence and the possibility of marine current orientated gene flow on the genetic features of sea beet populations, we assessed patterns of genetic structuring of cytoplasmic and nuclear diversity by: (i) searching for an isolation-by-distance (IBD) pattern using spatial autocorrelation tools; (ii) using the Monmonier algorithm to identify genetic boundaries in the area studied; and (iii) performing assignment tests that are based on multilocus genotype information to ascertain population membership of individuals. Our results showed a highly contrasted cytoplasmic and nuclear genetic differentiation and highlighted the peculiar situation of island populations. Beyond a classical isolation-by-distance due to short-range dispersal, genetic barriers fitting the orientation of marine currents were clearly identified. This suggests the occurrence of long-distance seed dispersal events and an asymmetrical gene flow separating the eastern and western part of the Anglo-Norman gulf.

Testing the role of genetic factors across multiple independent invasions of the shrub Scotch broom (Cytisus scoparius)

Knowledge of the introduction history of invasive plants informs on theories of invasiveness and assists in the invasives management. For the highly successful invasive shrub Scotch broom, Cytisus scoparius, we analysed a combination of nuclear and chloroplast microsatellites for eight native source regions and eight independent invasion events in four countries across three continents. We found that two exotic Australian populations came from different sources, one of which was derived from multiple native populations, as was an invasive sample from California. An invasive population from New Zealand appeared to be predominantly sourced from a single population, either from the native or exotic ranges. Four invasive populations from Chile were genetically differentiated from the native range samples analysed here and so their source of introduction could not be confirmed, but high levels of differentiation between the Chilean populations suggested a combination of different sources. This extensive global data set of replicated introductions also enabled tests of key theories of invasiveness in relation to genetic diversity. We conclude that invasive populations have similar levels of high genetic diversity to native ranges; levels of admixture may vary across invasive populations so admixture does not appear to have been an essential requirement for invasion; invasive and native populations exhibit similar level of genetic structure indicating similar gene flow dynamics for both types of populations. High levels of diversity and multiple source populations for invasive populations observed here discount founder effects or drift as likely explanations for previously observed seed size differences between ranges. The high levels of genetic diversity, differential and source admixture identified for most exotic populations are likely to limit the ability to source biocontrol agents from the native region of origin of invasive populations.

Inheritance of chloroplast DNA is not strictly maternal in Silene vulgaris (Caryophyllaceae): evidence from experimental crosses and natural populations

Chloroplast DNA (cpDNA) is maternally inherited in the majority, but not all, of angiosperm species. The mode of inheritance of cpDNA is a critical determinant of its molecular evolution and of its population genetic structure. Here, we present the results of investigations of the inheritance of cpDNA in Silene vulgaris, a plant used in a variety of studies in which cpDNA is an important component. PCR/RFLP markers were used to compare mother and offspring cpDNA genotypes sampled from two natural populations, and mother, father, and offspring genotypes obtained from controlled greenhouse crosses. Ten of 215 offspring cpDNA genotypes studied in the controlled crosses and three of 156 offspring from natural populations did not match that of the mother, demonstrating rare nonmaternal inheritance. That the chloroplast genome is occasionally transmitted through pollen is discussed in the context of using S. vulgaris cpDNA as a marker in studies of seed dispersal and when considering the joint evolution of the chloroplast and mitochondrial genomes.

Historical range expansion determines the phylogenetic diversity introduced during contemporary species invasion

For a species rapidly expanding its geographic range, such as during biological invasion, most alleles in the introduced range will have their evolutionary origins in the native range. Yet, the way in which historical processes occurring over evolutionary time in the native range contribute to the diversity sampled during contemporary invasion is largely unknown. We used chloroplast DNA (cpDNA) gene genealogies and coalescent methods to study two congeneric plants, Silene latifolia and S. vulgaris. We examined how phylogenetic diversity was shaped by demographic growth and historical range expansions in the native European range, and how this history affected the diversity sampled during their recent invasion of North America. Genealogies from both species depart from neutrality, likely as a result of demographic expansion in the ancestral range, the timing of which corresponds to shortly after each species originated. However, the species differ in the spatial distribution of cpDNA lineages across the native range. Silene latifolia shows a highly significant phylogeographic structure that most likely reflects different avenues of the post-glacial expansion into northern Europe from Mediterranean refugia. By contrast, cpDNA lineages in S. vulgaris have been widely scattered across Europe during, or since, the most recent post-glacial expansion. These different evolutionary histories resulted in dramatic differences in how phylogenetic diversity was sampled during invasion of North America. In S. latifolia, relatively few, discrete invasion events from a structured native range resulted in a rather severe genetic bottleneck, but also opportunities for admixture among previously isolated lineages. In S. vulgaris, lack of genetic structure was accompanied by more representative sampling of phylogenetic diversity during invasion, and reduced potential for admixture. Our results provide clear insights into how historical processes may feed forward to influence the phylogenetic diversity of species invading new geographic ranges.

Horizontal gene transfer of a plastid gene in the non-photosynthetic flowering plants Orobanche and Phelipanche (Orobanchaceae)

Plastid sequences are among the most widely used in phylogenetic and phylogeographic studies in flowering plants, where they are usually assumed to evolve like non-recombining, uniparentally transmitted, single-copy genes. Among others, this assumption can be violated by intracellular gene transfer (IGT) within cells or by the exchange of genes across mating barriers (horizontal gene transfer, HGT). We report on HGT of a plastid region including rps2, trnL-F, and rbcL in a group of non-photosynthetic flowering plants. Species of the parasitic broomrape genus Phelipanche harbor two copies of rps2, a plastid ribosomal gene, one corresponding to the phylogenetic position of the respective species, the other being horizontally acquired from the related broomrape genus Orobanche. While the vertically transmitted copies probably reside within the plastid genome, the localization of the horizontally acquired copies is not known. With both donor and recipient being parasitic plants, a possible pathway for the exchange of genetic material is via a commonly attacked host.

Nonneutral evolution of organelle genes in Silene vulgaris

Knowledge of mitochondrial gene evolution in angiosperms has taken a dramatic shift within the past decade, from universal slow rates of nucleotide change to a growing realization of high variation in rates among lineages. Additionally, evidence of paternal inheritance of plant mitochondria and recombination among mitochondrial genomes within heteroplasmic individuals has led to speculation about the potential for independent evolution of organellar genes. We report intraspecific mitochondrial and chloroplast sequence variation in a cosmopolitan sample of 42 Silene vulgaris individuals. There was remarkably high variation in two mitochondrial genes (atp1, atp9) and additional variation within a third gene (cob). Tests for patterns of nonneutral evolution were significant for atp1 and atp9, indicative of the maintenance of balanced polymorphisms. Two chloroplast genes (matK, ndhF) possessed less, but still high, variation and no divergence from neutral expectations. Phylogenetic patterns of organelle genes in both the chloroplast and mitochondria were incongruent, indicating the potential for independent evolutionary trajectories. Evidence indicated reassociation among cytoplasmic genomes and recombination between mitochondrial genes and within atp1, implying transient heteroplasmy in ancestral lineages. Although the mechanisms for long-term maintenance of mitochondrial polymorphism are currently unknown, frequency-dependent selection on linked cytoplasmic male sterility genes is a potential candidate.

Correlates of genetic differentiation and isolation by distance in 17 congeneric Silene species

The contemporary pattern of intraspecific genetic variation can indicate the relative role of gene flow and local differentiation in shaping the evolutionary history and future trajectory of a species. To assess the recent influence of contrasting life history and demographic characteristics on genetic structure within a group of closely related species, patterns of genetic differentiation (F(ST) and related statistics) and isolation by distance (IBD) were compared among 17 congeneric herbaceous plant species. Data came from 35 published studies of 16 species, and a previously unpublished analysis of chloroplast genetic variation in the rare endemic Silene rotundifolia. Among-population genetic variance was most strongly influenced by the type of genetic marker used; cytoplasmic markers showed larger values than allozyme and anonymous nuclear markers. Other independently significant factors were geographical range size and, for allozyme studies, reproductive system; in particular, endemism and hermaphroditism were associated with higher among-population genetic variance, whereas large native geographical range and dioecy were associated with lower among-population variance. Over equivalent spatial scales, dioecious populations also showed weaker IBD than hermaphrodites, perhaps because increased population transience and/or variance in the spatial pattern of gene flow are more closely associated with dioecy in this genus. Invasive populations had both highly variable among-population genetic variance, and no evidence for IBD, consistent with nonequilibrium conditions. Other analysed factors including predominant pollinator had no discernable influence on genetic structure or patterns of IBD. In general, this comparative approach appears to be valuable for synthesizing the complementary information provided by F-statistics and IBD, and for indicating the relative importance of particular biological factors in shaping genetic variation within different species of a closely related plant group.

Inbreeding effects on progeny sex ratio and gender variation in the gynodioecious Silene vulgaris (Caryophyllaceae)

In gynodioecious species, sex expression is generally determined through cytoplasmic male sterility genes interacting with nuclear restorers of the male function. With dominant restorers, there may be an excess of females in the progeny of self-fertilized compared with cross-fertilized hermaphrodites. Moreover, the effect of inbreeding on late stages of the life cycle remains poorly explored. Here, we used hermaphrodites of the gynodioecious Silene vulgaris originating from three populations located in different valleys in the Alps to investigate the effects of two generations of self- and cross-fertilization on sex ratio and gender variation. We detected an increase in females in the progeny of selfed compared with outcrossed hermaphrodites and inbreeding depression for female and male fertility. Male fertility correlated positively with sex ratio differences between outbred and inbred progeny, suggesting that dominant restorers are likely to influence male fertility qualitatively and quantitatively in S. vulgaris. We argue that the excess of females in the progeny of selfed compared with outcrossed hermaphrodites and inbreeding depression for gamete production may contribute to the maintenance of females in gynodioecious populations of S. vulgaris because purging of the genetic load is less likely to occur.

Multiple origin of metallicolous populations of the pseudometallophyte Arabidopsis halleri (Brassicaceae) in central Europe: the cpDNA testimony

The population structure of the pseudo-metallophyte herb, Arabidopsis halleri, was studied using Polymerase Chain Reaction-Restriction Fragment Length Polymorphism (PCR-RFLP) on chloroplast DNA (cpDNA). The history of metallicolous (M) populations showing increased zinc tolerance was investigated. Eight primer-enzyme combinations out of 72 tested were applied to a total of 625 individuals from 28 widespread populations, 14 of them being M. Eleven distinct chlorotypes were found: five were common to nonmetallicolous (NM) and M populations, whereas six were only observed in one edaphic type (five in NM and one in M). No difference in chlorotype diversity between edaphic types was detected. Computed on the basis of chlorotype frequencies, the level of population differentiation was high but remained the same when taking into account levels of molecular divergence between chlorotypes. Isolation by distance was largely responsible for population differentiation. Geographically isolated groups of M populations were more genetically related to their closest NM populations than to each other. Our results suggest that M populations have been founded separately from distinct NM populations without suffering founding events and that the evolution towards increased tolerance observed in the distinct M population groups occurred independently.

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.

Mitochondrial haplotypes and the New Zealand origin of clonal European Potamopyrgus, an invasive aquatic snail

The small aquatic snail Potamopyrgus antipodarum is an important invading species in Europe, Australia and North America. European populations are generally believed to derive from accidental introductions from New Zealand, probably dating back to the mid-19th century. We have employed mitochondrial DNA sequences to test the proposed New Zealand origin of European Potamopyrgus, and to learn more about its genealogical history. Using a 481-bp region of the 16S ribosomal RNA gene, we identified 17 distinct haplotypes among 65 snails from New Zealand. In marked contrast, only two haplotypes were found across all European samples, which cover a large geographical area. Importantly, these two haplotypes are shared with snails from the North Island of New Zealand. Due to sampling limitations we cannot rule out a South Island origin for one of the haplotypes, but our results clearly demonstrate the New Zealand origin of European populations. The marked divergence among the two European haplotypes implies the successful colonization by two distinct mitochondrial lineages, which is consistent with previous data based on nuclear markers.

Molecular evidence for multiple introductions of garlic mustard (Alliaria petiolata, Brassicaceae) to North America

Invasive species offer excellent model systems for studying rapid evolutionary change. In this context, molecular markers play an important role because they provide information about pathways of introduction, the amount of genetic variation introduced, and the extent to which founder effects and inbreeding after population bottlenecks may have contributed to evolutionary change. Here, we studied microsatellite variation in eight polymorphic loci among and within 27 native and 26 introduced populations of garlic mustard (Alliaria petiolata), a European herb which is a current serious invader in North American deciduous forests. Overall, introduced populations were genetically less diverse. However, considerable variability was present and when compared to the probable source regions, no bottleneck was evident. Observed heterozygosity was very low and resulted in high inbreeding coefficients, which did not differ significantly between native and introduced populations. Thus, selfing seems to be equally dominant in both ranges. Consequently, there was strong population differentiation in the native (F(ST) = 0.704) and the introduced (F(ST) = 0.789) ranges. The high allelic diversity in the introduced range strongly suggests multiple introductions of Alliaria petiolata to North America. Out of six European regions, the British Isles, northern Europe, and central Europe had significantly higher proportions of alleles, which are common to the introduced range, and are therefore the most probable source regions. The genetic diversity established by multiple introductions, and the lack of inbreeding depression in this highly selfing species, may have contributed to the invasion success of Alliaria petiolata.

Comparison between mitochondrial and chloroplast DNA variation in the native range of Silene vulgaris

A detailed survey of mitochondrial and chloroplast diversity in eight populations of Silene vulgaris from Central Europe was conducted for comparison with previously published data on diversity from S. vulgaris populations in the introduced range. Mitochondrial DNA (mtDNA) variation around the coxI gene was assessed with Southern blotting/restriction fragment length polymorphism methods. Chloroplast variation was assessed by sequencing the intergenic spacer separating the trnH and psbA genes. Thirty mtDNA haplotypes and 24 chloroplast DNA (cpDNA) haplotypes were found within 86 individuals. The overall genetic diversity h (0.941 for mitochondrial, and 0.893 for chloroplast markers) and within-population diversity were higher than reported in previous population studies of S. vulgaris in the USA and Europe. The frequency of private alleles was surprisingly high - more than 90% for both kinds of markers. Most of our populations were large and located in relatively undisturbed meadows, whereas surveys in Virginia consisted of smaller roadside populations. The slow rate of population turnover in European populations is discussed as a factor responsible for the relatively high diversity of S. vulgaris in undisturbed areas of its native range. Association between mtDNA and cpDNA haplotypes was also demonstrated. Finally, gender and mtDNA haplotype were associated in the Alps populations, where females were very rare.