Glacial survival does not matter – II: RAPD phylogeography of Nordic Saxifraga cespitosa

It has been suggested that many arctic-alpine plant species have limited dispersal ability and cannot have arrived in Scandinavia and the arctic archipelago of Svalbard by long-distance dispersal after a total glaciation. It has therefore been proposed that such species must have survived the entire glaciation(s) in ice-free refugia in southern Norway, northern Norway and Svalbard. We investigated random amplified polymorphic DNA (RAPD) variation among 28 populations from Norway and Svalbard of one of these arctic-alpine 'short-distance dispersers', the selfing polyploid Saxifraga cespitosa . In an analysis of molecular variance ( AMOVA ), more variation was found among populations within the three postulated refugia regions (45%) than among these regions (25%). Spatial autocorrelation (Mantel) analyses showed that the genetic distance monotonously increased with increasing geographical distance. In UPGMA and PCO analyses, the populations from Norway and Svalbard formed a south–north cline that continued across the Barents Sea barrier. The results suggest that there has been recent dispersal among the three postulated refugia regions and thus that postglacial dispersal into these refugia regions from other distant areas also must represent a possibility. The observed geographical pattern of the genetic variation may have been established after expansion from different source areas outside the North European ice sheet and/or from different refugia areas 'within' the ice sheet, but it is probably not possible to distinguish among these alternatives. The results for S. cespitosa are consistent with a dynamic late- and postglacial scenario with extensive plant dispersal, and support the conclusion from our previous study of the outbreeding Saxifraga oppositifolia ; the hypothesis of glacial survival in Norway and/or Svalbard is superfluous.     

Molecular diversity and derivations of populations of Silene acaulis and Saxifraga oppositifolia from the high Arctic and more southerly latitudes

A survey of allozyme diversity within and between populations of Silene acaulis from Spitsbergen, Norway, Iceland and Scotland, showed that populations from the high Arctic (Spitsbergen, > 76°N) contained high levels of diversity and were genetically similar to populations from more southern locations. Indirect measures of gene flow (Nm), calculated from Wrigh's F indicated that there had been extensive gene flow between Spitsbergen and some Norwegian populations. A restriction site analysis of chloroplast DNA (cpDNA) in S. acaulis revealed that all populations contained a single identical cpDNA haplotype, except one population from Norway which also contained a second haplotype. In contrast, five different cpDNA haplotypes were distinguished in a more limited survey of cpDNA variation in Saxifraga oppositifolia, with all five haplotypes present in one of two Spitsbergen populations surveyed. The contrasting cpDNA results for the two species suggest that whereas high-Arctic populations of Silene acaulis have most likely been derived from immigrants which arrived from the south after the last glacial period, high-Arctic populations of Saxifraga oppositifolia may be derived, in part, from ancient northern stocks which survived the last glaciation in high-Arctic refugia.     

Molecular analysis of the Pleistocene history of Saxifraga oppositifolia in the Alps

A recent circumpolar survey of chloroplast DNA (cpDNA) haplotypes identified Pleistocene glacial refugia for the Arctic-Alpine Saxifraga oppositifolia in the Arctic and, potentially, at more southern latitudes. However, evidence for glacial refugia within the ice sheet covering northern Europe during the last glacial period was not detected either with cpDNA or in another study of S. oppositifolia that surveyed random amplified polymorphic DNA (RAPD) variation. If any genotypes survived in such refugia, they must have been swamped by massive postglacial immigration of periglacial genotypes. The present study tested whether it is possible to reconstruct the Pleistocene history of S. oppositifolia in the European Alps using molecular methods. Restriction fragment length polymorphism (RFLP) analysis of cpDNA of S. oppositifolia, partly sampled from potential nunatak areas, detected two common European haplotypes throughout the Alps, while three populations harboured two additional, rare haplotypes. RAPD analysis confirmed the results of former studies on S. oppositifolia; high within, but low among population genetic variation and no particular geographical patterning. Some Alpine populations were not perfectly nested in this common gene pool and contained private RAPD markers, high molecular variance or rare cpDNA haplotypes, indicating that the species could possibly have survived on ice-free mountain tops (nunataks) in some parts of the Alps during the last glaciation. However, the overall lack of a geographical genetic pattern suggests that there was massive immigration of cpDNA and RAPD genotypes by seed and pollen flow during postglacial times. Thus, the glacial history of S. oppositifolia in the Alps appears to resemble closely that suggested previously for the species in northern Europe.     

Genetic Variation of Sesleria albicans Kit. ex Schultes (Poaceae): Lack of Evidence for Glacial Relict Endemism in Central Europe

Abstract: The populations of Sesleria albicans Kit. ex Schultes in central Europe are supposed to be glacial relicts. The genetic structure of these populations can provide evidence for the reliability of this hypothesis, since a distinct geographical clustering of populations and high levels of interpopulational genetic variation can be expected in the case of relict endemism. In the study presented here, we used Random Amplified Polymorphic DNA (RAPD) analysis to describe the genetic structure of 25 Sesleria albicans populations from west and central Germany, southwest Germany and the Alps. In the RAPD analysis 344 fragments could be amplified, of which 95.9 % were polymorphic. The percentage of polymorphic bands per population ranged from 29.7 % to 56.7 % and correlated significantly with population size, indicating a higher level of inbreeding in small and isolated populations. An analysis of molecular variance (AMOVA) revealed that 61.22 % of the variation resided within populations, whereas 34.25 % was found among populations within regions and only 4.53 % among regions. Genetic distance among populations (φ_ST) varied between 0.63 and 0.13, and in a Neighbour-Joining dendrogram based on the genetic distances between the investigated individuals, all populations could be discriminated from each other. The populations from the northern and southern Alps were clearly separated from the central European populations. These, however, did not form geographic clusters. Considering these results, Sesleria albicans showed a high level of intra- but only a low level of interpopulational genetic variation and very weak regional differentiation. The genetic pattern detected in this study, therefore, gives no evidence for glacial relict endemism of Sesleria albicans in central Europe.     

Genetic variation in the Afromontane tree Prunus africana,an endangered medicinal species

An understanding of the patterns of variation within and among populations of tropical trees is essential for devising optimum genetic management strategies for their conservation and sustainable utilization. Here, random amplified polymorphic DNA (RAPD) analysis was used to partition variation within and among 10 populations of the endangered Afromontane medicinal tree, Prunus africana, sampled from five countries across the geographical range of the species (Cameroon, Ethiopia, Kenya, Madagascar and Uganda). Analysis of molecular variance ( AMOVA ) employed 48 RAPD markers and revealed most variation among countries (66%, P < 0.001). However, variation among individuals within populations and among populations within Cameroon and Madagascar was also highly significant. Analysis of population product frequency data indicated Ugandan material to be more similar to populations from Cameroon than populations from Kenya and Ethiopia, while Malagash populations were most distinct. The implications of these findings for determining appropriate approaches for conservation of the species, particularly in Cameroon and Madagascar, are discussed.     

Isozyme analysis of theCerastium alpinum-C. arcticum complex (Caryophyllaceae) supports a splitting ofC. arcticum Lange

As part of a larger investigation of theC. alpinum-C. arcticum complex of arctic and North Atlantic areas, isozyme variation ofC. alpinum, C. arcticum, and related taxa was analysed. A total of 124 multilocus phenotypes was divided into more or less distinct groups by numerical analyses. Most groups correspond well to previously recognized taxa. However, what has traditionally been considered asC. arcticum was divided into two distinct groups, consisting of northern (Svalbard, Greenland) and more southern (Norway, Iceland) populations, respectively. The division ofC. arcticum into two taxa is also supported by other kinds of data and the two taxa probably deserve species rank. Serpentine plants from Shetland had multilocus phenotypes similar to those ofC. arcticum from Iceland and should be included in the southern taxon.     

Genetic variation in Fitzroya cupressoides (alerce), a threatened South American conifer

Fitzroya cupressoides (alerce, Cupressaceae) is a large and exceptionally long-lived conifer, endemic to a restricted area of southern Chile and neighbouring areas of Argentina. As a result of its high economic value, the species has been severely exploited for timber, and remnant populations are fragmented and often highly disturbed. The species is thought to have undergone a major range contraction during the last glaciation. In order to assess the extent of genetic variation using DNA markers within and between populations of this species, samples were obtained from throughout the natural range and analysed for random amplified polymorphic DNA (RAPD) variation. Eight 10-mer and three 15-mer primers were used to produce a total of 54 polymorphic bands. Shannon's diversity estimates were calculated to provide an estimate of the degree of variation within each population. Values varied from 0.343 to 0.636 with only the lowest value differing significantly from the others (Spop = 0.547). This indicated that there is a significant degree of variation within each population, and did not provide evidence for genetic 'bottle-neck' effects within the species. A pairwise distance measure calculated from the RAPD data was used as an input for principal coordinate (PCO) and AMOVA analyses. The first three principal coordinates of RAPD distances described 8.3, 5.9 and 5.4% of the total variance, respectively, and a degree of clustering of samples according to their geographical origin was detectable. AMOVA analysis indicated that although most of the variation (85.6%) was found within populations, a significant proportion (P < 0.002) was attributable to differences between populations. An UPGMA dendrogram constructed using phi ST values derived from AMOVA produced a pattern broadly similar to that produced by the PCO, highlighting differences between three main groups of populations within Chile: those from the northern Coastal Range, the southern Coastal Range and Central Depression, and the Andes. Populations from Argentina also emerged as significantly different from those in Chile. These results are interpreted in the context of the postglacial history of the species, and their implications for the development of conservation strategies for Fitzroya are discussed.     

Genetic structure and diversity in Ramonda myconi (Gesneriaceae): effects of historical climate change on a preglacial relict species

The importance of the Mediterranean Basin as a long-term reservoir of biological diversity has been widely recognized, although much less effort has been devoted to understanding processes that allow species to persist in this area. Ramonda myconi (Gesneriaceae) is a Tertiary relict plant species restricted to the NE Iberian Peninsula. We used RAPD and chloroplast markers to assess the patterns of genetic structure in eight mountain regions covering almost the full species range, to identify the main historical processes that have shaped its current distribution and to infer the number and location of putative glacial refugia. While no cpDNA variation was detected, the species had relatively high levels of RAPD variation. Maximum levels of diversity were found within populations (71%), but there was also a significant differentiation between geographical regions (20%) and among populations within regions (9%). A spatial AMOVA identified three main groups of populations, corresponding to previously recognized centers of endemism and species richness. In addition, we found a marked geographical pattern of decreasing genetic diversity and increasing population differentiation from west to east. Our results support a complex phylogeographic scenario in the Iberian Peninsula of "refugia-within-refugia" and suggest that the higher diversity observed in western regions might be associated with prolonged and more stable climatic conditions in this area during the Quaternary.     

Genetic variation of ticks (Ixodes ricinus L.) in the Lithuanian and Norwegian populations

RAPD markers were used to measure the genetic diversity of 119 individuals of Ixodes ricinus collected from Lithuania and Norway. The samples were analysed within and also between the populations. We analysed 74 loci in each of 6 populations. Our results show high levels of diversity within the populations. The percentage of polymorphic loci of the six analysed populations: Birzai, Vilnius, Kretinga, Tjore, Kjosvik and Odderoya were 68.9%, 58.1%, 78.38%, 62.2%, 44.6% and 68.9%, respectively. The percentage of polymorphic loci in the Lithuanian populations was 93.2%, and in the Norwegian populations 81.08%. The genetic distance ranged from 0.019 to 0.079 within Norwegian populations and from 0.005 to 0.0967 within Lithuanian populations and between the countries from 0.022 to 0.146. The genetic variation of I. ricinus among Norwegian populations was lower than among Lithuanian populations. The highest part of genetic variation in I. ricinus ticks depends on variation within Kretinga (Lithuania) and Odderoya (Norway) populations situated in coastal areas where many migratory and sea birds are aggregated.     

Genetic variation of Saxifraga paniculata Mill. (Saxifragaceae): molecular evidence for glacial relict endemism in central Europe

Saxifraga paniculata Mill. is an amphi-atlantic species that is widely distributed in the Alps but also occurs with disjunct populations in central Europe. These isolated populations are considered to be glacial relicts. In the study presented here, we analysed the genetic variation within and between populations of Saxifraga paniculata in central Europe to test whether the molecular data give evidence for the supposed relict endemism. We used RAPD analysis and the analysis of molecular variance (AMOVA) technique to investigate 30 populations of Saxifraga paniculata from different locations in central Europe and adjacent geographical areas. In the RAPD analysis, 319 fragments could be amplified of which 91.2% were polymorphic. The percentage of polymorphic bands within populations correlated significantly with population size, indicating a higher level of inbreeding in small populations. AMOVA revealed 43.60% of the variation within populations, only 15.45% between populations within locations, and 40.95% between locations. Genetic distance between locations (Φ_ST) correlated significantly with the geographical distance between locations, giving evidence for geographical isolation. In an UPGMA cluster analysis, based on RAPD data, the populations were grouped together according to their location. The results exhibit a strong genetic differentiation which is obviously due to genetic drift in the isolated populations. Our study therefore gives evidence for glacial relict endemism of Saxifraga paniculata in central Europe.  © 2003 The Linnean Society of London, Biological Journal of the Linnean Society , 2003, 80, 11–21.     

Trans-Atlantic dispersal and phylogeography of Cerastium arcticum (Caryophyllaceae) inferred from RAPD and SCAR markers

Cerastium arcticum is an autogamous pioneer species with a distribution limited to the North Atlantic region. It has been suggested that such species must have survived in ice-free refugia on both sides of the Atlantic throughout the last, or even several, of the Pleistocene glaciations, because they lack special adaptations for long-distance dispersal. To address the possibility for recent trans-Atlantic dispersal of C. arcticum, we analyzed random amplified polymorphic DNA (RAPD) and sequence characterized amplified region (SCAR) differentiation among 26 populations of this high-polyploid species. Three SCAR markers were obtained that verified the main patterns identified in the RAPD analysis. Eighty-four multilocus RAPD phenotypes were observed in the 126 plants analyzed, based on 35 polymorphic markers. Multivariate analyses and analyses of molecular variance (AMOVAs) identified two highly divergent groups of populations: one arctic group (western and eastern Greenland, and the archipelagos of Svalbard and Franz Josef Land) and one nonarctic group (southern and northern Norway, and Iceland), indicating that C. arcticum is composed of two lineages with different evolutionary histories. However, there was little geographic structuring within each lineage, in spite of the fact that both lineages are disjunctly distributed across the Atlantic. Occurrence of very similar, in some cases even identical RAPD multilocus phenotypes on both sides of the Atlantic in this autogamous allopolyploid is most probably caused by postglacial dispersal. The present geographic distribution of C. arcticum may thus have been established after trans-Atlantic expansion from two Weichselian refugia, one for each evolutionary lineage. Unexpectedly, the level of intrapopulational variation increased towards the north. This may reflect that interpopulational migration is most extensive in the treeless arctic environment, where the species has a more continuous distribution than in the more southerly areas.     

Random amplified polymorphic DNA variation among remnant big bluestem (Andropogon gerardii vitman) populations from Arkansas' grand prairie

Random amplified polymorphic DNA (RAPD) analysis was used to characterize genetic diversity and genetic distinctiveness of Andropogon gerardii from remnant Arkansas prairies. Six oligonucleotide primers, which generated 37 RAPD bands, were used to analyse 30-32 plants from six Grand Prairie populations, Baker Prairie (Arkansas Ozarks), two Illinois prairies and two cultivars. Genetic diversity of the Arkansas remnants ranged from 82.7 to 99.3%, with 89% of the total genetic variation within and 11% among populations. The partitioning of genetic variation was consistent with that reported for other outcrossing perennial grasses, using the more conservative allozyme markers. Principal component analysis indicated a northern and southern association within Arkansas' Grand Prairie. Although there was no genetic structuring at the landscape level, the Illinois prairies and cultivars were different from all Arkansas prairies tested. There was significant within-population structuring in four of the seven Arkansas remnants, with a negative relationship between genetic similarity and geographical distance. The three nonstructured populations were from a linear railroad remnant, suggesting different population-level dynamics from nonlinear prairies. The results of this study indicated that small isolated remnant big bluestem populations were not genetically depauperate and that genetic relationships among populations could not be predicted solely on geographical proximity.     

Diversity and genetic differentiation among populations of Indian and Kenyan tea (Camellia sinensis (L.) O. Kuntze) revealed by AFLP markers

AFLP markers were successfully employed to detect diversity and genetic differentiation among Indian and Kenyan populations of tea (Camellia sinensis (L.) O. Kuntze). Shannon's index of diversity was used to partition the total phenotypic variation into between and within population components. On average, most of the diversity was detected within populations, with 79% of the variation being within and 21% being between populations of Indian and Kenyan tea. A dendrogram constructed on the basis of band sharing distinctly separated the three populations of tea into China type (sinensis), Assam type (assamica) and Cambod type (assamica ssp. lasiocalyx) in a manner consistent with the present taxonomy of tea, the known pedigree of some of the genotypes and their geographical origin. Principal coordinate (PCO) analysis grouped Assam genotypes both from India and Kenya supporting the suggestion that the Kenyan clones have been derived from collections made in this region. The China types were more dispersed on the PCO plot which is a reflection of wider genetic variation. As would be expected, clones collected from the same region exhibited less overall genetic variation. AFLP analysis discriminated all of the tested genotypes from India and Kenya, even those which cannot be distinguished on the basis of morphological and phenotypic traits. Received: 2 May 1996 / Accepted: 14 June 1996     

High RAPD but no cpDNA sequence variation in the endemic and endangered plant, Heptacodium miconioides Rehd. (Caprifoliaceae)

Heptacodium miconioides Rehd. is an endangered species endemic to China and has suffered rapid decrease of distribution range and population size. This species has been disappeared in central China where the modal specimen was collected. We analyzed the genetic variation of the remaining populations to reveal whether the genetic diversity also suffered decrease and to provide some suggestions for conservation. All the nine known remaining populations were sampled. Genetic variation was analyzed based on RAPD markers and two fragments of cpDNA sequence, intergenic spacers of petG-trnP and trnS-trnG. No variation was observed in the two fragments of cpDNA sequence. However, the species exhibited high level of RAPD variation compared to other threatened or rare plants. Measures of genetic diversity within populations were strongly related to the log of estimated population size, indicating that large populations usually have more genetic diversity than that of small ones. About 25% of the variation was partitioned among populations. Significant relationship was observed between differentiation and geographical distance, indicating a pattern of isolation-by-distance. Given for few populations remaining, all the populations should be protected and urgent efforts be paid on the small populations to avoid their local extinction.     

Analysis of the variability of Drosophila azteca and D. athabasca populations revealed by randomly amplified polymorphic DNA

The disribution ranges of Drosophila azteca and D. athabasca overlap in northen California and southern Oregon. Seven populations, four of which are located in this area, were studied. Large random amplified polymorphic DNA (RAPD) variation was found within species; nevertheless, more than half the primers used in the study yielded greater diofference between than within species. A nested analysis of molecular variance (AMOVA) showed that the variance between populations within species was significantly greater than zero for 55% of the oligonucleotides used, which provided evidence for an underlying geographical structure of these populations. Specimens of D. azteca and D. athabasca from Salem (OR), where both species were collected together, presented the highest similarity encountered between species.     

Genetic variation and biogeography of Mertensia maritima (Boraginaceae)

To elucidate the population dynamics and biogeographic connections of Mertensia maritima in the North Atlantic we studied the genetic variation within and among 16 local populations in the Oslofjord region (south-east Norway), Finnmark (northern Norway), Iceland and Svalbard, using starch gel isozyme electrophoresis. We also tested the self-compatibility experimentally. Mertensia maritima was found to be self-compatible and able to self-pollinate by autodeposition. The total level of genetic variation was extremely low. All the screened individuals were homozygous at all loci. All individuals from the populations of mainland Norway and Iceland were isozymatically identical, whereas individuals from Svalbard differed from these at three loci (AAT-1, PGM-3 and SKD-2). The genetic identity between any population at Svalbard and any population on the mainland was 0.86. These genetic patterns may have resulted from a combination of inbreeding, frequent local population extinctions, and long-distance dispersal. The isozymatic distinctness of the Svalbard populations is consistent with the arctic distribution of the morphologically based var. tenella , here raised to subspecific rank as Mertensia maritima ssp. tenella.     

Molecules and morphology in concert. II. The Festuca brachyphylla complex (Poaceae) in Svalbard

We used a combined molecular and morphological approach to unravel variation in the autogamous Festuca brachyphylla polyploid complex in the arctic archipelago of Svalbard. Forty populations were analyzed for random amplified polymorphic DNAs (RAPDs) and 46 morphological characters. Eighteen RAPD multilocus phenotypes were observed in the 86 plants analyzed, based on 30 polymorphic markers. Multivariate analyses of the RAPD data revealed four distinct groups of multilocus phenotypes; in contrast, the variation was more or less continuous in multivariate analyses of the morphological data. However, we identified several individual morphological characters that unambiguously discriminated among the four groups of RAPD multilocus phenotypes. Analysis of type material suggests that the four groups in Svalbard can be referred to Festuca baffinensis, F. brachyphylla, F. hyperborea, and F. edlundiae. This study shows that concerted analysis of molecules and morphology is a powerful tool in low-level taxonomy.     

Genetic variation in Antarctic populations of the moss Sarconeurum glaciale

Sixty-six isolates of the moss Sarconeurum glaciale were collected from sites in continental Antarctica at Ross Island, southern Victoria Land and the Vestfold Hills. Genetic variation within and among the populations was estimated using isozymes and random amplified polymorphic DNA (RAPD) technology. Isozyme results only reproducibly showed variation between the populations with one enzyme; RAPDs indicated significantly higher levels of genetic variability within and among the Vestfold Hills samples than in the Ross Sea region samples. A dendrogram produced from the RAPD bands suggested that the Ross Island and southern Victoria Land samples form one population, and those from the Vestfold Hills form a separate and more variable population, possibly resulting from separate colonisation events on the continent. Received: 15 March 1996 / Accepted: 1 May 1997     

Multilocus phylogeography and population structure of common eiders breeding in North America and Scandinavia

Aim Glacial refugia during the Pleistocene had major impacts on the levels and spatial apportionment of genetic diversity of species in northern latitude ecosystems. We characterized patterns of population subdivision, and tested hypotheses associated with locations of potential Pleistocene refugia and the relative contribution of these refugia to the post‐glacial colonization of North America and Scandinavia by common eiders (Somateria mollissima). Specifically, we evaluated localities hypothesized as ice‐free areas or glacial refugia for other Arctic vertebrates, including Beringia, the High Arctic Canadian Archipelago, Newfoundland Bank, Spitsbergen Bank and north‐west Norway. Location Alaska, Canada, Norway and Sweden. Methods Molecular data from 12 microsatellite loci, the mitochondrial DNA (mtDNA) control region, and two nuclear introns were collected and analysed for 15 populations of common eiders (n = 716) breeding throughout North America and Scandinavia. Population genetic structure, historical population fluctuations and gene flow were inferred using F‐statistics, analyses of molecular variance, and multilocus coalescent analyses. Results Significant inter‐population variation in allelic and haplotypic frequencies were observed (nuclear DNA F_ST = 0.004–0.290; mtDNA Φ_ST = 0.051–0.927). Whereas spatial differentiation in nuclear genes was concordant with subspecific designations, geographic proximity was more predictive of inter‐population variance in mitochondrial DNA haplotype frequency. Inferences of historical population demography were consistent with restriction of common eiders to four geographic areas during the Last Glacial Maximum: Belcher Islands, Newfoundland Bank, northern Alaska and Svalbard. Three of these areas coincide with previously identified glacial refugia: Newfoundland Bank, Beringia and Spitsbergen Bank. Gene‐flow and clustering analyses indicated that the Beringian refugium contributed little to common eider post‐glacial colonization of North America, whereas Canadian, Scandinavian and southern Alaskan post‐glacial colonization is likely to have occurred in a stepwise fashion from the same glacial refugium. Main conclusions Concordance of proposed glacial refugia used by common eiders and other Arctic species indicates that Arctic and subarctic refugia were important reservoirs of genetic diversity during the Pleistocene. Furthermore, suture zones identified at MacKenzie River, western Alaska/Aleutians and Scandinavia coincide with those identified for other Arctic vertebrates, suggesting that these regions were strong geographic barriers limiting dispersal from Pleistocene refugia.     

Taxonomic relationships in East AsianVicia species with unijugate leaves based on random amplified polymorphic DNA markers

Random amplified polymorphic DNA(RAPD) markers were investigated to clarify the taxonomic positions ofVicia linearifolia andV. bifolia, and to assess the genomic diversity among the 9 populations ofV. unijuga, each of which represents a geographical variation or infraspecific taxa in southern Korea. These species are characterized by unijugate leaves in East Asia and have been controversial as to infra-or interspecific classification. The polymorphic markers among the populations examined were observed for fifteen decamer primers. The degree of band sharing was used to calculate genetic similarity between populations, and a phenogram using UPGMA cluster analysis was generated based on the Dice similarity coefficient. The taxa studied were divided into two main groups and the populations ofV. unijuga were all grouped together in the phenogram. The genetic similarities ofV. unijuga were very high among the populations and did not show distinctions between the infraspecific taxa, although the populations of Mt. Odae and adjacent areas in eastern Korea were different from others of the species.V. linearifolia fell within the range of the genomic variation among the populations ofV. unijuga, whileV. bifolia was grouped withV. venosa var.cuspidata having multijugate leaves rather thanV. unijuga. The result from studying RAPD markers suggested thatV. linearifolia should be integrated intoV. unijuga and that species with unijugate leaves ofV. bifolia andV. unijuga are polyphyletic.