Potential for parthenogenesis of virgin females in a bisexual population of the geographically parthenogenetic mayfly Ephoron shigae (Insecta: Ephemeroptera,Polymitarcyidae)

The burrowing polymitarcyid mayfly Ephoron shigae is a geographically parthenogenetic species. Interestingly, the distributions of the bisexual and unisexual populations overlap broadly in their respective geographic ranges. In this mayfly, obligatory diploid thelytoky appears within unisexual populations. In the present study, we examined the potential for parthenogenesis or the parthenogenetic ability of females in a bisexual population aiming to understand the emergence of unisexual populations. The results obtained revealed that females in the examined bisexual populations showed a potential for diploid thelytoky as also seen in the unisexual populations, although, in females from bisexual populations, the development success rates of their unfertilized eggs were considerably lower than those of virgin females from unisexual populations. In the three bisexual reproducing species (Ephemera japonica, Ephemera strigata, and Ephemera orientalis) in the closely‐related family Ephemeridae, diploid thelytoky (i.e. tychoparthenogenesis; < 3%) was also observed. However, in this case, the parthenogenetic development success rates of unfertilized eggs were significantly lower than those of virgin females in the bisexual (Hino‐yosui Irrigation Canal) population of E. shigae. Accordingly, we suggest that parthenogenetic ability (i.e. tychoparthenogenesis or facultative parthenogenesis) in bisexual populations of E. shigae may facilitate the evolutionary transition to unisexual populations with fully obligatory parthenogenesis. © 2009 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 326–334.     

Reproductive mode of the geographic parthenogenetic mayfly Ephoron shigae, with findings from some new localities (Insecta: Ephemeroptera, Polymitarcyidae)

The burrowing polymitarcyid mayfly Ephoron shigae is distributed widely in Japan. Some populations are bisexual, others are unisexual, and the distributions of the two types overlap broadly. Experimental evidence of parthenogenetic reproduction, long suspected in unisexual populations, is presented here, based on a comparative analysis of the developmental rate of fertilized and unfertilized eggs. The developmental rate of fertilized eggs from 20 mated females in a bisexual population was 98.4% ± 0.73% (mean ± SD), and no unfertilized eggs from 20 virgin females in that population developed. The developmental rates of unfertilized eggs in two unisexual populations were 89.0% ± 4.59% and 84.2% ± 1.96%, respectively. This article presents experimental evidence of geographic parthenogenesis in E. shigae and provides support for the previous interpretation. In addition, we discuss the relationship between the sex ratio of each population and the developmental rate of fertilized versus unfertilized eggs from the females in those populations.     

Phylogeography of the East Asian polymitarcyid mayfly genus Ephoron (Ephemeroptera: Polymitarcyidae): a comparative analysis of molecular and ecological characteristics

In East Asia, the mayfly genus Ephoron, a member of Polymitarcyidae, has been classified into three species: Ephoron eophilum, Ephoron limnobium, and Ephoron shigae. Using the nuclear histone 3 and internal transcribed spacer 1 and mitochondrial 16S rRNA and cytochrome c oxidase subunit I genes, we investigated the phylogenetic relationships of these mayflies. Ephoron limnobium formed polyphyletic relationships within the clade of Japanese E. shigae and shared identical sequences with E. shigae, suggesting that E. limnobium is a synonym of E. shigae. Monophyly was supported in E. eophilum, Korean E. shigae, and E. limnobium + Japanese E. shigae. Ephoron eophilum that co‐occurs with the other Ephoron species has several peculiar traits in its adult swarming period and karyotype. This suggests that both pre‐ and post‐zygotic isolation are likely to have promoted speciation of these Ephoron species. The divergence between E. eophilum and E. shigae would have occurred from the late Late Miocene to Early Pliocene according to the Bayesian relaxed clock model calibrated using the mitochondrial substitution rate. The area currently inhabited by E. eophilum emerged above sea level some 0.1 Mya, suggesting that E. eophilum moved into this area from other areas that it occupied prior to this time. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 109 , 181–202.     

Automictic parthenogenesis of a geographically parthenogenetic mayfly,Ephoron shigae (Insecta: Ephemeroptera,Polymitarcyidae)

In the geographically parthenogenetic mayfly, Ephoron shigae, egg maturation and counts of chromosome number of unfertilized, parthenogenetic eggs were studied, in comparison with fertilized eggs from a bisexual population. The primary oocyte becomes mature through two successive maturation divisions. The first maturation division (meiotic division) takes place in the primary oocyte to produce a secondary oocyte and a first polar body. The second maturation division soon occurs in the secondary oocyte, in which the nucleus is divided into a mature egg nucleus (female pronucleus) and second polar body nuclei. The first polar body, in some cases, was successively divided into two polar bodies; in other instances, it was not divided. After the successive maturation division, the egg nucleus and the sister second polar body nucleus drew near to fuse into the zygote nucleus. The chromosome number was doubled in the zygote, and this conjugation initiates subsequent embryonic development. This suggests that, in E. shigae, the process of parthenogenetic recovery of diploidy is the automictic type categorized as the ‘terminal fusion’ pattern. © 2009 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 335–343.     

Genetic variation and Wolbachia infection status of the lycaenid butterfly Zizina emelina (Lepidoptera,Lycaenidae) from Korean populations

The lycaenid butterfly Zizina emelina emelina (de l’Orza) (previously Zizina otis emelina; Lepidoptera, Lycaenidae) has populations distributed across Japan and Korea. Zizina emelina is listed on Japan’s Red List as an endangered species, owing to the loss of its principal plant food and habitat. Researchers previously investigated the genetic diversity and structure in Japanese populations by using mitochondrial and nuclear DNA analyses alongside their Wolbachia infection status to consider the conservation strategy. In this study, we investigated the genetic variation and Wolbachia infection status of Korean populations and compared our findings with these traits in the Japanese populations. Four novel mitochondrial haplotypes were found in the Korean populations, whereas nuclear DNA analysis revealed one new haplotype and another that matched the one reported previously in Japan. The two known strains of Wolbachia in the Japanese populations were also found in the Korean ones. Mitochondrial DNA analysis demonstrated that no gene flow had occurred between the Japanese and Korean populations. The fact that they share the same nuclear haplotype and Wolbachia strains, however, indicates a low level of population differentiation between the two.     

Geographic variations in Pacific sand eels <Emphasis Type="Italic">Ammodytes personatus</Emphasis> (Ammodytidae) from Korea and Japan using multivariate morphometric analysis

Multivariate methods were used to compare meristic and morphometric characters of Ammodytes personatus from four populations in Korea (3) and Japan (1). ANCOVA showed the slopes of head length (Japan > East > South > West) and body depth (South > Japan > East > West) relative to total length to differ significantly among populations. Vertebrae and fin rays’ numbers differed significantly among populations, the former ranging between 64 and 67 (mode = 65, 66) in the East population, 62–67 (63, 64) in the West population and 62–64 (62, 63) in the South population in Korea. PCA and CDA, based on 14 morphometric characters, showed that scores plotted on axes 1-2 and 2-3 separated the three Korean populations but overlapped between the South and Japanese populations. These results suggest that the Korean A. personatus may comprise at least two populations. The text was submitted by the authors in English.     

Genetic structure of Korean populations of bumblebees Bombus ignitus (Hymenoptera: Apidae) as revealed by microsatellite markers

Bumblebee Bombus ignitus, which is indigenous to Korea, Japan, and China, has been recognized as a valuable pollinator for both crops and wild plants. Bombus ignitus has now become commercially important as a pollinator because of its use in the agricultural industry, particularly for greenhouse pollination. For long‐term management and effective conservation of B. ignitus, an understanding of the genetic structure of its naturally occurring populations is practically important. In this study, the genetic structure among the five populations of B. ignitus in South Korea was assessed using nine microsatellite loci. These markers showed high levels of genetic variability, with the number of alleles ranging from 6 to 22 (mean = 13.4) and the frequency of the most common allele ranging from 0.11 to 0.66. Only the Sabuk (SB) population showed a genetic signature of a recent bottleneck, which was further supported by the lowest level of allelic richness (AR) (mean AR = 3.944). Genetic differentiation was highly significant among all pairs of populations (P < 0.001) across the nine microsatellite markers, suggesting a lack of gene flow among those populations. Interestingly, F_ST (and R_ST and D_est) values were always greater for the Taebaek population than for the four remaining populations. The phylogenetic analysis showed evidence supporting our hypothesis that the Taebaek population is genetically more divergent than the other populations. Overall, our results suggest that the Korean populations of B. ignitus might have undergone geographic isolation and have become highly separated spatially from one another, thereby having a limited range of migration among their local habitats.     

Genetic homogeneity between Korean and Japanese populations of the broad-leaved evergreen tree Machilus thunbergii (Lauraceae): A massive post-glacial immigration through the Korea Strait or something else?

The warm-temperate vegetation of Korea, currently limited to southern coastal areas, shifted during the Last Glacial Maximum (LGM) towards glacial refugia, putatively located in southern Japan. We hypothesized two scenarios of post-glacial re-colonization of warm-temperate species in relation to current levels of genetic diversity within their populations. If extant Korean populations originated from a single source (a single glacial refugium), we expect significantly lower levels of genetic diversity relative to those from Japan due to founder effects. Alternatively, if they were derived from multiple source populations, levels of genetic diversity within Korean populations will not be significantly reduced compared to those of Japanese ones. To test which of these scenarios is more likely, we investigated the patterns of genetic diversity in 14 populations (seven from southern Korea and seven from southern Japan) of the broad-leaved evergreen tree Machilus thunbergii, employing 11 allozyme loci. High levels of genetic variation in M. thunbergii were found both in southern Korea and southern Japan, with a considerable genetic homogeneity not only between the two regions but also between populations within the two regions. These results suggest a pattern of re-colonization after the LGM fitting the second scenario (immigration from multiple refugia), probably through multiple waves and/or with large founder populations.     

Molecular and morphological differences between two geographic populations of Salanx ariakensis (Salangidae) from Korea and Japan

Geographic variations of a salangid, Salanx ariakensis, from Korea and Japan, based on partial mtDNA cyt b gene sequences and morphometric data, are presented and compared with Salanx cuvieri from China. Analyses of molecular variance (AMOVA) suggested that the two populations of S. ariakensis differed significantly from each other in pairwise fixation index (F _ST = 0.43), although their reciprocal monophyly was not constructed in the phylogenetic analysis. The phylogenetic analysis, however, supported the reciprocal monophyly of S. ariakensis (Korean population + Japanese population) and S. cuvieri. Differences in genetic structure between the two populations of S. ariakensis were suggested by diversity indices, mismatch distribution shape, and minimum spanning network. In particular, diversity indices indicated that the Korean population was historically larger and more stable than the Japanese population. In addition, both the principal component analysis (PCA) and canonical discriminant analysis (CDA) for morphometric characters, and the Kruskal–Wallis test for some meristic characters, showed that the two populations of S. ariakensis differed from each other, although morphological differences between the two populations of S. ariakensis were smaller than those between S. ariakensis and S. cuvieri. Such morphological differences were consistent with the differences in mtDNA.     

Genetic status of Asiatic black bear (Ursus thibetanus) reintroduced into South Korea based on mitochondrial DNA and microsatellite loci analysis

The Asiatic black bear is one of the most endangered mammals in South Korea owing to population declines resulting from human exploitation and habitat fragmentation. To restore the black bear population in South Korea, 27 bear cubs from North Korea and Russian Far East (Primorsky Krai) were imported and released into Jirisan National Park, a reservoir of the largest wild population in South Korea, in 2004. To monitor the success of this reintroduction, the genetic diversity and population structure of the reintroduced black bears were measured using both mitochondrial and nuclear DNA markers. Mitochondrial D-loop region DNA sequences (615 bp) of 43 Japanese black bears from previous study and 14 Southeast Asian black bears in this study were employed to obtain phylogenetic inference of the reintroduced black bears. The mitochondrial phylogeny indicated Asiatic black bear populations from Russian Far East and North Korea form a single evolutionary unit distinct from populations from Japan and Southeast Asia. Mean expected heterozygosity (H(E)) across 16 microsatellite loci was 0.648 for Russian and 0.676 for North Korean populations. There was a moderate but significant level of microsatellite differentiation (F(ST) = 0.063) between black bears from the 2 source areas. In addition, genetic evidences revealed that 2 populations are represented as diverging groups, with lingering genetic admixture among individuals of 2 source populations. Relatedness analysis based on genetic markers indicated several discrepancies with the pedigree records. Implication of the phylogenetic and genetic evidences on long-term management of Asiatic black bears in South Korea is discussed.     

MICROSATELLITE MARKERS REVEAL POPULATION GENETIC STRUCTURE OF THE TOXIC DINOFLAGELLATE ALEXANDRIUM TAMARENSE (DINOPHYCEAE) IN JAPANESE COASTAL WATERS1

This is the first report to explore the fine‐scale diversity, population genetic structure, and biogeography of a typical planktonic microbe in Japanese and Korean coastal waters and also to try to detect the impact of natural and human‐assisted dispersals on the genetic structure and gene flow in a toxic dinoflagellate species. Here we present the genetic analysis of Alexandrium tamarense (Lebour) Balech populations from 10 sites along the Japanese and Korean coasts. We used nine microsatellite loci, which varied widely in number of alleles and gene diversity across populations. The analysis revealed that Nei's genetic distance correlated significantly with geographic distance in pair‐wise comparisons, and that there was genetic differentiation in about half of 45 pair‐wise populations. These results clearly indicate genetic isolation among populations according to geographic distance and restricted gene flow via natural dispersal through tidal currents among the populations. On the other hand, high P‐values in Fisher's combined test were detected in five pair‐wise populations, suggesting similar genetic structure and a close genetic relationship between the populations. These findings suggest that the genetic structure of Japanese A. tamarense populations has been disturbed, possibly by human‐assisted dispersal, which has resulted in gene flow between geographically separated populations.     

Genetic integrity and individual identification-based population size estimate of the endangered long-tailed goral,Naemorhedus caudatus from Seoraksan National Park in South Korea,based on a non-invasive genetic approach

ABSTRACT

The long-tailed goral (also called the Amur goral) Naemorhedus caudatus (subfamily Caprinae), a vulnerable and protected species designated by IUCN and CITES, has sharply been declining in the population size and is now becoming critically endangered in South Korea. This species has been conserved as a natural monument by the Korean Cultural Heritage Administration since 1968. In this study, using 78 fecal DNA samples with a non-invasive genetic approach, we assessed the genetic integrity and individual identification-based population size for the goral population from Seoraksan National Park representing the largest wild population in Korea. Using the successfully isolated 38 fecal DNA, phylogeographic and population genetic analyses were performed with mitochondrial DNA control region (CR) sequences and nine microsatellite loci. We found seven CR haplotypes, of which five were unique to the Seoraksan population, considering previously determined haplotypes in Korean populations. The Seoraksan population showed higher haplotype diversity (0.777?±?0.062) and mean number of alleles (4.67?±?1.563) relative to southern populations in Korea reported from previous studies, with no signal of a population bottleneck. Microsatellite-based individual identification estimate based on probability of identity (PID) indicated a population size of ≥30 in this population. Altogether, we suggest that for future management efforts of this species in the Seoraksan National Park, conserving its genetic integrity as an ‘endemic’ lineage, and curbing a decrease in its number through mitigating habitat destruction might be key to secure the population for the long term.     

PHYLOGEOGRAPHY OF CUTLERIA CYLINDRICA (CUTLERIALES,PHAEOPHYCEAE) IN NORTHEASTERN ASIA,AND THE IDENTITY OF AN INTRODUCED POPULATION IN CALIFORNIA1

Cutleria cylindrica Okamura was described from Japan in 1902 and has been reported only from northwestern Asia until its relatively recent discovery in California, USA, and Baja California, Mexico. To clarify the genetic relationships within and among the disjunct populations, we carried out a molecular phylogenetic study, as well as the examination of sex ratio and the life‐history patterns, of populations in Japan, Korea, and California. Based on the DNA sequences of mitochondrial genes cox2, cox3, the open reading frame (ORF) region, and the spacer between cox3 and ORF, a total of 23 haplotypes were detected in the 85 individuals from 20 localities in Japan, Korea, and California. All localities in Japan and Korea included multiple haplotypes, but only a single haplotype was found in California. There was a positive relationship between distance and genetic divergence in Japan and Korea. The single haplotype found in California was the same as one occurring in Japan (Aomori Pref. and Fukuoka Pref.) and Korea (Daedaepo, Pusan). Both male and female gametophytes were distributed in most northeastern Asian populations. Only female gametophytes, developing parthenogenetically from female gametes, were found in California and Aomori Pref., Japan. On the basis of these results, we conclude that the disjunct population of C. cylindrica in California originated from a relatively recent introduction from Japan and shares its origin with the parthenogenetic population in the Tsugaru Strait.     

Population structure of chum salmon and selection on the markers collected for stock identification

Genetic stock identification (GSI) is a major management tool of Pacific salmon (Oncorhynchus Spp.) that has provided rich genetic baseline data of allozymes, microsatellites, and single‐nucleotide polymorphisms (SNPs) across the Pacific Rim. Here, we analyzed published data sets for adult chum salmon (Oncorhynchus keta), namely 10 microsatellites, 53 SNPs, and a mitochondrial DNA locus (mtDNA3, control region, and NADH‐3 combined) in samples from 495 locations in the same distribution range (n = 61,813). TreeMix analysis of the microsatellite loci identified the greatest convergence toward Japanese/Korean populations and suggested two admixture events from Japan/Korea to Russia and the Alaskan Peninsula. The SNPs had been purposively collected from rapidly evolving genes to increase the power of GSI. The largest expected heterozygosity was observed in Japanese/Korean populations for microsatellites, whereas it was largest in Western Alaskan populations for SNPs, reflecting the SNP discovery process. A regression of SNP population structures on those of microsatellites indicated the selection of the SNP loci according to deviations from the predicted structures. Specifically, we matched the sampling locations of the SNPs with those of the microsatellites and performed regression analyses of SNP allele frequencies on a 2‐dimensional scaling (MDS) of matched locations obtained from microsatellite pairwise F _ST values. The MDS first axis indicated a latitudinal cline in American and Russian populations, whereas the second axis showed differentiation of Japanese/Korean populations. The top five outlier SNPs included mtDNA3, U502241 (unknown), GnRH373, ras1362, and TCP178, which were identified by principal component analysis. We summarized the functions of 53 nuclear genes surrounding SNPs and the mtDNA3 locus by referring to a gene database system and propose how they may influence the fitness of chum salmon.     

Phylogeography of Pulsatilla cernua (Ranunculaceae), a grassland species,in Japan

The genetic diversity and structure of Pulsatilla cernua, a continental‐grassland relict, were investigated using variations in chloroplast DNA (cpDNA) and microsatellites of nuclear DNA. In the analyses of three cpDNA regions, 17 haplotypes were found in 24 populations of P. cernua from Japan, Korea, and Russia. Although the route and time of migration between the continent of Asia and Japan could not be well resolved, the cpDNA haplotype network suggests the existence of several ancient lineages in Japan and a recent secondary migration from Japan to the continent. Microsatellite analyses did not indicate genetic structure among the Japanese populations, indicating the existence of gene flow across the distribution area until recently. These results indicate that the present fragmentation of P. cernua in Japan may reflect a rapid, recent reduction from a previously large, continuous distribution.     

Application of Microsatellite Markers to Population Genetics Studies of Japanese Flounder Paralichthys olivaceus

We examined population genetic structure by means of microsatellite analysis among 7 Japanese flounder (Paralichthys olivaceus) populations collected from coastal sea areas around Japan. As was expected, all of the 11 microsatellite loci examined were variable in all populations (number of alleles per locus, 15.2–18.2; average of expected heterozygosity, 0.74–0.78). Eleven population pairs in 21 possible pairwise comparisons showed significant genetic heterogeneity associated with allele frequency distributions or fixation index (F _ST ). Modified Cavalli-Sforza chord distance (D _A ) and Nei's standard genetic distance (D _ST ) ranged from 0.051 to 0.090, and from 0.000 to 0.025, respectively. There was evidence that the populations assessed in this study were not drawn from a single panmictic population; however, it appears that Japanese flounder populations around Japan are not well-structured, as an estimate of the fixation index value among the 7 localities was very low (F _ST = 0.0025). Received April 27, 2001; accepted June 27, 2001.     

Phylogeography of the Korean pine (Pinus koraiensis) in northeast Asia: inferences from organelle gene sequences

Range-wide genetic variation of Korean pine (Pinus koraiensis) was assessed using maternally inherited mtDNA and paternally inherited cpDNA for 16 natural populations throughout northeast Asia in order to study its phylogeographical history during the Quaternary. The cpDNA variation indicated that there was no difference between populations on the Asian continent and those in the Japanese archipelago. In contrast, the mtDNA variation indicated that there was significant difference between the populations from the two regions, with each region having a different lineage. The continental populations exhibited no diversity in the mtDNA examined despite the species’ current extensive range and large populations. Conversely, while the Korean pine is rare in Japan, the Japanese populations exhibited greater levels of mtDNA diversity (H _T?=?0.502). The higher mtDNA diversity and evidence from numerous Korean pine macrofossil remains dated to the Pleistocene and recovered various sites in Japan suggest that the Japanese archipelago once served as a refugium to a much larger Korean pine population with a more extensive range than is the case today. The presence of the single mtDNA haplotype across the Asian continent suggests that the present widespread populations could have expanded from a single refugium population after the last glacial periods.     

Allozyme Diversity in Populations of Cymbidium goeringii (Orchidaceae)

Abstract: Using 14 allozyme loci, we investigated levels of genetic diversity within populations, and degree of genetic divergence among 24 populations of Cymbidium goeringii (Orchidaceae) in Korea and Japan. Cymbidium goeringii maintains high levels of genetic diversity both at population (mean expected heterozygosity, H _e = 0.238) and species levels (0.260). Means of H _e found in 24 populations were not significantly different from each other. About 90 % of the total variation in the species is common to all populations (mean G _ST = 0.108). No unique allele was found in any population. The indirect estimate of gene flow based on the mean G _ST was high ( Nm = 2.06). Nei's genetic identities for pairs of populations had high values (mean = 0.974 [SD = 0.013]). The Mantel-Z test showed a significant correlation between genetic distance and geographic distance. However, the mean G _ST value between 17 populations in Korea and seven Japanese populations was relatively low (0.029), even though the land connection between the southern Korean peninsula and southern Japanese archipelagos has not existed since the middle Pleistocene. Large numbers of small seeds of C. goeringii might travel long distances by wind from populations to populations both in Korea and Japan, increasing genetic diversity within populations and maintaining low genetic differentiation among populations.     

Two invasive populations of the dry rot fungus Serpula lacrymans show divergent population genetic structures

The dry rot fungus Serpula lacrymans is a devastating basidiomycete occurring in wooden constructions in temperate regions worldwide. In this study, we compare the genetic structures of two invasive populations from Europe and Japan. Microsatellite data from 14 loci and DNA sequences from four loci demonstrated that the two populations were highly differentiated. Significant isolation by distance effect was observed in Europe and Japan. Higher genetic variation was observed within the Japanese population than within the European population, corresponding with the observed higher richness of vegetative compatibility types in Japan, indicating that there has been a higher level of gene flow from the Asian source populations to Japan than to Europe. The European population is genetically more homogenous with only six detected vegetative compatibility types. Various tests indicate that both the European and the Japanese populations have gone through population bottlenecks prior to population expansion. No identical multi‐locus genotypes were observed within Japan and very few within Europe, indicating limited clonal dispersal. Deviations from Hardy Weinberg expectations were observed both in Europe and Japan and heterozygote excess were observed at several loci, especially in Europe. Possible explanations for this pattern are discussed.     

Phylogeography of East Asian <Emphasis Type="Italic">Lespedeza buergeri</Emphasis> (Fabaceae) based on chloroplast and nuclear ribosomal DNA sequence variations

The dynamic changes in land configuration during the Quaternary that were accompanied by climatic oscillations have significantly influenced the current distribution and genetic structure of warm-temperate forests in East Asia. Although recent surveys have been conducted, the historical migration of forest species via land bridges and, especially, the origins of Korean populations remains conjectural. Here, we reveal the genetic structure of Lespedeza buergeri, a warm-temperate shrub that is disjunctively distributed around the East China Sea (ECS) at China, Korea, and Japan. Two non-coding regions (rpl32-trnL, psbA-trnH) of chloroplast DNA (cpDNA) and the internal transcribed spacer of nuclear ribosomal DNA (nrITS) were analyzed for 188 individuals from 16 populations, which covered almost all of its distribution. The nrITS data demonstrated a genetic structure that followed geographic boundaries. This examination utilized AMOVA, comparisons of genetic differentiation based on haplotype frequency/genetic mutations among haplotypes, and Mantel tests. However, the cpDNA data showed contrasting genetic pattern, implying that this difference was due to a slower mutation rate in cpDNA than in nrITS. These results indicated frequent migration by this species via an ECS land bridge during the early Pleistocene that then tapered gradually toward the late Pleistocene. A genetic isolation between western and eastern Japan coincided with broad consensus that was suggested by the presence of other warm-temperate plants in that country. For Korean populations, high genetic diversity indicated the existence of refugia during the Last Glacial Maximum on the Korean Peninsula. However, their closeness with western Japanese populations at the level of haplotype clade implied that gene flow from western Japanese refugia was possible until post-glacial processing occurred through the Korea/Tsushima Strait land bridge.