Specialization clines in the pollination systems of agaves (Agavaceae) and columnar cacti (Cactaceae): A phylogenetically controlled meta-analysis

The biogeography of plant-animal interactions is a novel topic on which many disciplines converge (e.g., reproductive biology, biogeography, and evolutionary biology). Narrative reviews have indicated that tropical columnar cacti and agaves have highly specialized pollination systems, while extratropical species have generalized systems. However, this dichotomy has never been quantitatively tested. We tested this hypothesis using traditional and phylogenetically informed meta-analysis. Three effect sizes were estimated from the literature: diurnal, nocturnal, and hand cross-pollination (an indicator of pollen limitation). Columnar cactus pollination systems ranged from purely bat-pollinated in the tropics to generalized pollination, with diurnal visitors as effective as nocturnal visitors in extratropical regions; even when phylogenetic relatedness among species is taken into account. Metaregressions identified a latitudinal increase in pollen limitation in columnar cacti, but this increase was not significant after correcting for phylogeny. The currently available data for agaves do not support any latitudinal trend. Nectar production of columnar cacti varied with latitude. Although this variation is positively correlated with pollination by diurnal visitors, it is influenced by phylogeny. The degree of specificity in the pollination systems of columnar cacti is heavily influenced by ecological factors and has a predictable geographic pattern.     

Pollination ecology of Stenocereus queretaroensis (Cactaceae), a chiropterophilous columnar cactus, in a tropical dry forest of Mexico

Flowers of columnar cacti are animal-pollinated, often displaying a chiropterophylic syndrome. This study examined if the columnar cactus Stenocereus queretaroensis, a tropical species endemic to western Mexico, is bat-pollinated, by studying its pollination biology and the foraging behavior of potential pollinators. Flowers were produced in winter through spring, peaking in April. Anthesis was nocturnal, and stigma and anther turgidity began around 2200 hours. Production of nectar secretion and highest sugar concentration and energy supply were nocturnal, peaking between 2200 and 2400 hours. Manual auto-pollination and exclusion experiments showed that self-pollination yielded no fruits, while nocturnal pollinators resulted in high fruit set and seed set compared to diurnal pollination treatments. The nectar-feeding bat Leptonycteris curasoae (Phyllostomidae) was the main nocturnal pollinator with the highest effective pollination. Peak bat visitation coincided with peaks in nectar production. The high abundance of L. curasoae throughout the 4-yr study, suggests that it is a seasonally reliable pollinator for this columnar cactus. While pollination syndromes have been increasingly called into question in recent years, this study suggests that at least for this system, there is a fairly close fit between pollinator and pollination syndrome.     

Comparative pollination biology of Venezuelan columnar cacti and the role of nectar-feeding bats in their sexual reproduction

The floral biology, reproductive system, and visitation behavior of pollinators of four species of columnar cacti, Stenocereus griseus, Pilosocereus moritzianus, Subpilocereus repandus, and Subpilocereus horrispinus, were studied in two arid zones in the north of Venezuela. Our results support the hypothesis that Venezuelan species of columnar cacti have evolved toward specialization on bat pollination. Additional information on the floral biology of a fifth species, Pilosocereus lanuginosus, was also included. All species showed the typical traits that characterize the pollination syndrome of chiropterophily. All species but Pilosocereus moritzianus were obligate outcrossers. Nectar and pollen were restricted to nocturnal floral visitors. Two species of nectar-feeding bats, Leptonycteris curasoae Miller and Glossophaga longirostris Miller, were responsible for practically all the fruit set in these cacti. Frequency of bat visitation per flower per night was highly variable within and between species of cactus, with average frequencies varying between 27 and 78 visits/flower/night. In general terms, the pattern of floral visitation through the night was significantly correlated with the pattern of nectar production and nectar sugar concentration for all species of cactus. Under natural pollination, fruit:flower ratios varied from 0.46 in Subpilocereus repandus to 0.76 in Stenocereus griseus. The efficiency of bat pollination in terms of seed:ovule ratio was high in all species, varying between 0.70 and 0.94.     

Analysis of form and function in North American columnar cacti (tribe Pachycereeae)

Simple and multivariate linear models were used to demonstrate the influence of mechanical design and climate on stem morphology and branching architecture in 25 species of North American columnar cacti. The effect of phyletic inertia was tested by the method of independent contrasts. Stem girth was found to increase significantly slower with increased height within taxa (cross-sectional stem area ;ps [plant height] 0.603), than across taxon (cross-sectional stem area ;ps [plant height] 1.451). Juveniles are shown to be mechanically overbuilt and subsequently grow into more slender adult forms determined in part by structural limitations and the optimization of other stem functions. We make a structural analogy of relatively rigid columnar cacti to concrete columns and compare plants and models with similar growth forms lacking woody skeletons (barrel cacti). Taxa with woody support achieved a surface-to-volume ratio six times greater than taxa without woody support. Across taxon, cooler winter temperatures were associated with larger stem girths, and greater annual precipitation was associated with less frequent branching. The relationship between total plant surface and volume approaches isometry within taxa, but across taxon average individuals are scaled replicates. We hypothesize that architecture and average plant height are adjusted, in an evolutionary sense, to maintain geometric similitude between surface and volume along a climatic gradient.     

Genetic relationships among Hylocereus and Selenicereus vine cacti (Cactaceae): evidence from hybridization and cytological studies

BACKGROUND AND AIMS: Hylocereus and Selenicereus are native to tropical and sub-tropical America. Based on its taxonomic status and crossability relations it was postulated that H. megalanthus (syn. S. megalanthus) is an allotetraploid (2n = 4x = 44) derived from natural hybridization between two closely related diploid taxa. The present work aimed at elucidating the genetic relationships between species of the two genera. METHODS: Crosses were performed and the putative hybrids were analysed by chromosome counts and morphological traits. The ploidy level of hybrids was confirmed by fluorescent in situ hybridization (FISH) of rDNA sites. Genomic in situ hybridization (GISH) was used in an attempt to identify the putative diploid genome donors of H. megalanthus and an artificial interploid hybrid. KEY RESULTS: Reciprocal crosses among four diploid Hylocereus species (H. costaricensis, H. monacanthus (syn. H. polyrhizus), H. undatus and Hylocereus sp.) yielded viable diploid hybrids, with regular chromosome pairing. Reciprocal crosses between these Hylocereus spp. and H. megalanthus yielded viable triploid, pentaploid, hexaploid and aneuploid hybrids. Morphological and phenological traits confirm the hybrid origin. In situ detection of rDNA sites was in accord with the ploidy status of the species and hybrid studied. GISH results indicated that overall sequence composition of H. megalanthus is similar to that of H. ocamponis and S. grandiflorus. High sequence similarity was also found between the parental genomes of H. monacanthus and H. megalanthus in one triploid hybrid. CONCLUSIONS: The ease of obtaining partially fertile F1 hybrids and the relative sequence similarity (in GISH study) suggest close genetic relationships among the taxa analysed.     

Population variability in the genusTrichocereus (Cactaceae) in central Chile

Trichocereus has 4 species in central Chile. Population analyses indicate that introgressive hybridization has occurred betweenT. chilensis andT. litoralis in the Valparaíso and Santiago coast ranges. A taxonomic key is presented.     

Patterns of abundance and population structure of Pachycereus pringlei (Cactaceae), a columnar cactus of the Sonoran Desert

Understanding the mechanisms that determine the distribution and abundance of plants is a major problem in ecology. However, very few studies have explored the factors controlling the abundance of columnar cacti throughout their range of distribution. In this paper, we describe the density and size structure of 26 populations of Pachycereus pringlei throughout its distribution range in the Sonoran Desert. Major differences in abundance were detected between island and mainland and peninsular areas, with islands sustaining significantly larger densities than mainland and peninsular populations. Within peninsular populations, the abundance was negatively associated with latitude and positively associated with annual and seasonal rainfall. In contrast, the abundance in mainland populations showed neither latitudinal trend nor an association with rainfall. In peninsular populations, mean height and basal diameter of branched plants showed a negative association with population density whereas mainland populations showed no significant association. None of the populations exhibited a population structure that fitted the log-normal distribution expected for young, growing populations with constant recruitment. Insular, peninsular and␣mainland populations showed a population structure with an uneven size distribution typical of populations experiencing regeneration pulses.     

Stem biomechanics of three columnar cacti from the Sonoran Desert

The allometric relationship of stem length L with respect to mean stem diameter D was determined for 80 shoots of each of three columnar cactus species (Stenocereus thurberi, Lophocereus schottii, and S. gummosus) to determine whether this relationship accords with that predicted by each of three contending models purporting to describe the mechanical architecture of vertical shoots (i.e., geometric, stress, and elastic similitude, which predict L proportional to D(alpha), with alpha = 1/1, 1/2, and 2/3, respectively). In addition, anatomical, physical, and biomechanical stem properties were measured to determine how the stems of these three species maintain their elastic stability as they increase in size. Reduced major axis regression of L with respect to D showed that alpha = 2.82 ± 0.14 for S. thurberi, 2.32 ± 0.19 for L. schottii, and 4.21 ± 0.31 for S. gummosus. Thus, the scaling exponents for the allometry of L differed significantly from that predicted by each of the three biomechanical models. In contrast, these exponents were similar to that for the allometry previously reported for saguaro. Analyses of biomechanical data derived from bending tests performed on 30 stems selected from each of the three species indicated that the bulk stem tissue stiffness was roughly proportional to L2, while stem flexural rigidity (i.e., the ability to resist a bending force) scaled roughly as L3. Stem length was significantly and positively correlated with the volume fraction of wood, while regression analysis of the pooled data from the three species (i.e., 90 stems) indicated that bulk tissue stiffness scaled roughly as the 5/3-power of the volume fraction of wood in stems. These data were interpreted to indicate that wood served as the major stiffening agent in stems and that this tissue accumulates at a sufficient rate to afford unusually high scaling exponents tot stem length with respect to stem diameter (i.e., disproportionately large increments of stem length with respect to increments in stem diameter). Nevertheless, the safety factor against the elastic failure of stems (computed on the basis of the critical buckling height divided by actual stem length) decreased with increasing stem size tot each species, even though each species maintained an average safety factor equal to two. We speculate that the apparent upper limit to plant height calculated for each species may serve as a biomechanical mechanism for vegetative propagation and the establishment of dense plant colonies by means of extreme stem flexure and ultimate breakage, especially for S. gummosus.     

Distribution and habitat in Mexico of Dactylopius Costa (Hemiptera: Dactylopiidae) and their cacti hosts (Cactaceae: Opuntioideae)

The distribution pattern of species of the genus Dactylopius Costa in Mexico was analyzed in relation to the distribution of their host plants (subfamily Opuntioideae) to evaluate the specificity of the insect-host association. The distribution of Dactylopius currently recognized is narrower than that of its hosts and probably is not representative. Therefore, a broader distribution of the Dactylopius species in correspondence with those of their hosts was hypothesized. Insects and their hosts were collected and georeferenced in 14 states of Mexico from 2005 to 2007. The distribution areas, maps, and habitat characteristics of Dactylopius, Opuntia sensu stricto, Nopalea and Cylindropuntia were determined on the basis of field collections and examination of museum collections. This information was complemented with information from the exhaustive examination of microscope slides from a local insect collection, plants from local herbaria, and literature reviews. The current distribution of the genus Dactylopius and its hosts included 22 and 25 states of Mexico, respectively, and Dactylopius had a continuous distribution according to its hosts, broader than recognized hitherto. The new georeferenced records of the five Mexican Dactylopius species are reported. Insects with morphological characteristics of D. confusus combined with those of D. salmianus were identified, as well as insects with characteristics of D. opuntiae combined with those of D. salmianus. These records suggest that the number of local Dactylopius species could be higher than previously thought or that possible new processes of hybridization between native and introduced species may be occurring.     

华木莲居群遗传结构与保护单元

华木莲(Sinomanglietia glauca)仅分布于江西宜春和湖南永顺, 是我国一级重点保护植物。前人采用RAPD、ISSR以及叶绿体SSCP(single-stranded conformation polymorphism)标记对华木莲进行了居群遗传学研究, 但未包括后发现的湖南居群或未检出居群内遗传变异。为了全面检测华木莲遗传多样性及其空间分布格局, 并据此确定保护单元, 本研究采用细胞核微卫星标记对华木莲所有4个居群共77个个体进行了居群遗传学分析。结果表明, 华木莲具有较低的遗传多样性(平均等位基因数A = 2.604, 平均期望杂合度H_E = 0.423)和较高的遗传分化(F_ST = 0.425)。STRUCTURE和主成分分析(Principal Coordinated Analysis, PCA)将4个居群首先分为江西、湖南两组, 江西的2个居群实际上是同一个自然繁育居群, 而湖南的2个居群则为2个分化明显的自然繁育居群。研究还发现湖南居群存在明显的杂合子过剩现象, 可能是小居群内随机因素造成的。研究结果表明华木莲可能在近期历史上遭受过强烈的瓶颈效应, 导致种群缩小、遗传多样性丧失和居群分化加剧, 需要加强对其进化潜力的保护。在制定保护措施时, 需要考虑其较高的遗传分化水平, 根据遗传结构可以将其划分为3个保护单元, 即湖南居群和江西居群分别为2个进化显著单元, 湖南居群进一步划分为2个管理单元(分别为朗溪乡云盘村和小溪乡鲁家村居群)。     

Population genetic structure of North American thinhorn sheep (Ovis dalli)

The thinhorn sheep (Ovis dalli ssp.) provides a rare example of a North American large mammal that occupies most of its native range and maintains close to ancestral population size. There are currently two recognized subspecies, Dall's sheep (O. d. dalli) and Stone's sheep (O. d. stonei), the validity of which remains uncertain. We investigated the spatial genetic structure of thinhorn sheep populations representing both subspecies by genotyping individuals (n = 919) from across the species range at 12 variable microsatellite loci. We found high levels of genetic diversity within (HE = 0.722) and significant genetic structure among the 24 sampled areas (FST = 0.160). Genetic distance measures and Bayesian clustering analyses revealed the presence of at least eight subpopulations that are delineated by mountain range topology. A strong overall pattern of isolation-by-distance is evident across the sampling range (r = 0.75, P < 0.001) suggesting limited dispersal and extensive philopatry. Partial Mantel tests of this relationship showed mountain range distinctions represent significant barriers to gene flow (P = 0.0001), supporting the Bayesian analyses. Genetic structure was more strongly pronounced in southern Yukon and Alaska than elsewhere. We also show evidence for genetic differences between the two currently recognized thinhorn subspecies.     

Allozyme diversity and genetic structure of the leafy cactus (Pereskia guamacho [Cactaceae])

We examined levels of genetic variation and genetic structure in the leafy cactus (Pereskia guamacho) in arid and semiarid zones in Venezuela. We surveyed genetic diversity within 17 populations using 19 allozyme loci. Genetic diversity was relatively high at both the species (P(s) = 89%, A(s) = 3.26, AP(s) = 3.53, H(es) = 0.24) and population (P(p) = 63%, A(p) = 1.90, AP(p) = 2.42, H(ep) = 0.20) levels. A significant deficit of heterozygote individuals was detected within populations in the Paraguana Peninsula region (F(IS) = 0.301). Relatively low levels of population differentiation were detected at macrogeographic (G(ST) = 0.112) and regional levels (G(ST) = 0.044 for peninsula region and G(ST) = 0.074 for mainland region), suggesting substantial genetic exchange among populations; however, gene flow in this species seems to be regulated by the distance among populations. Overall, estimates of genetic diversity found in P. guamacho are concordant with the pattern observed for other cacti surveyed, namely high levels of polymorphism and genetic diversity with one common allele and several rare alleles per locus. Differences in gene dispersal systems between this species and other cacti studied were not reflected in the patterns of genetic diversity observed. The concentration of the highest estimates of genetic variation in northwestern Venezuela suggests a potential reservoir of plant genetic diversity within xerophilous ecosystems in northern South America.     

Population genetic structure of the tropical tree species Aegiphila sellowiana (Lamiaceae)

The Tibagi River, located in southern Brazil, is associated with a significant degree of environmental heterogeneity, along its 550 km extension. There is concern about the integrity of this river's ecosystem, as human interference has been increasing. Aegiphila sellowiana (Lamiaceae) is an important pioneer tree species, commonly found near rivers; the fruit is consumed by avifauna. We studied this species along three ecological gradients, comprising the upper, middle, and lower regions of the Tibagi River basin. The genetic structure of nine subpopulations of A. sellowiana distributed along these gradients was investigated using RAPDs. Moderate levels of gene diversity (ranging from 0.091 to 0.132) were identified, inferred by a traditional approach and a Bayesian model-based method. The F-statistic, G(ST) parameters and molecular variance analysis showed high genetic differentiation among the three regions (39.5 to 50.26%). Analysis of molecular variance revealed high levels of genetic variation between populations (50.26%), while lower values of genetic variation (ranging from 9.56 to 16.35%) were seen between subpopulations within the upper, middle, and lower regions of the Tibagi River basin. The validity of these results was confirmed by principal coordinate analysis. Linear regression analysis showed significant correlations (r = 0.621, P = 0.0001) between the genetic and geographical distances. The differences observed in genetic variation between regions are probably due to habitat fragmentation; for conservation purposes, we recommend that at least one subpopulation from each region of the Tibagi River should be maintained.     

Population genetic structure of orchid bees (Euglossini) in anthropogenically altered landscapes

Habitat degradation and fragmentation are widespread phenomena in tropical regions. Negative effects on the biota are numerous, ranging from interruption of gene flow among populations, to the loss of genetic diversity within populations, to a decline in species richness over time. Orchid bees (Hymenoptera: Apidae: Euglossini) are of major conservation interest due to their function as pollinators of numerous orchid species and other tropical plants. Here, we used microsatellite markers to investigate the effects of geographic distance and habitat fragmentation on gene flow among populations. Populations of Euglossa dilemma in three geographic regions??the Yucat??n peninsula (Mexico), Veracruz (Mexico), and Florida (USA)??were genetically structured predominantly across the regions, with the strength of differentiation among populations being positively correlated with geographic distance. Within geographic regions only little substructure was found, suggesting that dispersal is substantial in the absence of geographic or ecological barriers. In a second study, patterns of genetic differentiation among eight species of Euglossa were not related to habitat fragmentation following deforestation in southern Mexico (Veracruz). Specifically, most bee populations in the 9,800?ha forest remnant of Los Tuxtlas (Volcano San Martin) were neither differentiated from, nor had less genetic diversity than, populations in near-continuous forest separated from Los Tuxtlas by 130?km of agricultural land. Either occasional long distance dispersal across open areas has buffered the expected genetic effects of fragmentation, or the history of fragmentation in southern Mexico is too recent to have caused measurable shifts in allelic composition.     

Population genetic structure of Yushania niitakayamensis (Bambusoideae, Poaceae) in Taiwan

Yushania niitakayamensis is distributed in Taiwan, south-west China and northern Philippines. In Taiwan, the species occurs in the central mountain ranges from 1500 to 3500 m in altitude. Morphological variation, especially in terms of plant height, is large, with plants ranging from 10 cm to 5 m in height. The species appears to spread mainly by rhizomes and flowers rarely, leading to the prediction that most populations are comprised of a single or a few clonal genotypes and that the observed morphological variation is primarily due to phenotypic plasticity. The purpose of the present study was to investigate the genetic structure of this species on Mt Hohuan in central Taiwan. Ten plants from a single clone and ten plants of unknown genetic background were surveyed at one site in order to select RAPD primers useful for clone identification. Plants at a second site were collected at 1-m intervals across a 50-m transect through the population. Plants at one extreme (exposed portion) of the transect were approximately 15–30 cm in height, whereas plants up to 410 cm in height were found at the other shaded end of the transect. Comparison of amplification profiles for 12 primers revealed that in contrast to our predictions of genetic uniformity, many samples had reproducibly different RAPD amplification profiles, with the 51 samples representing 31 clones. These data imply that the clone size is relatively small, and the population is actually highly diverse genetically. The genetic variation in this population may be due to a higher frequency of sexual reproduction during the evolutionary history of the species and/or a high somatic mutation rate for RAPD loci in clones of Yushania.     

Population genetic structure of the South American species Hypochaeris lutea (Asteraceae)

The genus Hypochaeris has a recent evolutionary history caused by long‐distance dispersal in conjunction with adaptive radiation in the South American continent. Hypochaeris lutea is a perennial herb that grows mostly at altitudes of around 1000 m in cold swamps of the southern regions of Brazil. We investigated the amplified fragment length polymorphism (AFLP) in 270 individuals representing 11 Brazilian populations of H. lutea to elucidate the population genetic structure of this species. The frequencies of polymorphic loci and gene diversity ranged from 83.42% to 91.66% and from 0.26 to 0.34, respectively. Analysis of molecular variance revealed that most of the genetic variability was found within (76.67%) rather than among (23.3%) populations, agreeing with the pattern of genetic distribution within and among populations observed in other allogamous species of Hypochaeris. A Mantel test showed no correlation between genetic and geographic distances when all populations were considered. Simulations performed using a Bayesian approach consistently identified two clusters with different admixture proportions of individuals, as also revealed by a UPGMA dendrogram of populations. The pattern of genetic structure observed in H. lutea is consistent with a process of successive colonization events by long‐distance dispersal resembling the rapid and recent radiation that has been proposed to explain the origin of the South American species of Hypochaeris.     

Population genetic structure of Ceriops tagal (Rhizophoraceae) in Thailand and China

Knowledge of the amount and patterns of genetic variation within and among populations of mangrove trees is essential for devising optimum genetic management strategies for their conservation and sustainable utilization. Ceriops tagal is a widespread viviparous mangrove. Genetic diversity in the species was examined with inter-simple sequence repeat (ISSR). Nine natural populations were collected from Thailand and China. The estimates of genetic variation were extremely low (H_T = 0.0179 ± 0.005, H_S = 0.0084 ± 0.001), and only 47% of the total gene diversity was maintained within populations (G_ST = 0.529). The eastern coastal populations of Thailand were more similar to populations from China than to populations from the western coastline of Thailand. A high level of Nei's genetic identity exists between populations of C. tagal (I = 0.989), suggesting their common ancestry. The low levels of genetic diversity in the species may result from a series of genetic bottlenecks during several glacial epochs.     

Pollination biology of two columnar cacti (Neobuxbaumia mezcalaensis and Neobuxbaumia macrocephala) in the Tehuacan Valley,central Mexico

We document the pollination biology and mating systems of Neobuxbaumia mezcalaensis and Neobuxbaumia macrocephala, two Mexican giant columnar cacti. These two species form mixed forests in the western Tehuacan Valley, Mexico. The flowers of both N. mezcalaensis and N. macrocephala exhibit nocturnal anthesis, are self-incompatible, and are pollinated primarily by three species of nectar-feeding bats (Choeronycteris mexicana, Leptonycteris curasoae, and Leptonycteris nivalis). Neobuxbaumia mezcalaensis is androdioecious, a breeding system that appears to be uncommon among Cactaceae. Neobuxbaumia macrocephala is hermaphroditic. We hypothesize that columnar cacti show a geographical dichotomy in floral biology specialization that probably can be related to predictability in pollinator abundance.     

Digestibility of columnar cacti pollen grains in the glosophagine bats Glossophaga longirostris and Leptonycteris curasoae (Chiroptera: Phyllostomidae)

We examined the protoplasmic assimilation of columnar cacti pollen grains in two species of Venezuelan desert glosophagine bats, Glossophaga longirostris and Leptonycteris curasoae, by determining the amount of empty (digested) pollen grains found in their fecal samples. To determine the amount of empty pollen grains, the fecal samples were stained to differ between empty and full (non-digested) pollen grains. The number of empty and full pollen grains observed in the fecal samples were corrected using the amount of aborted pollens present (before anthesis) in flowers of the columnar cacti species (Subpilocereus repandus, Stenocereus griseus and Pilosocereus tillianus) used by bats as food in the study site; G. longirostris and L. curasoae digested 64.2% and 71.3% of all the pollens fed, respectively. These high values confirm the importance of pollen in the diet of these bats, given its high nitrogen level.     

Population genetic structure and colonization history of short ninespine sticklebacks (Pungitius kaibarae)

The contemporary distribution and genetic structure of a freshwater fish provide insight into its historical geodispersal and geographical isolation following Quaternary climate changes. The short ninespine stickleback, Pungitius kaibarae, is a small gasterosteid fish occurring in freshwater systems on the Korean Peninsula and in southeast Russia. On the Korean Peninsula, P. kaibarae populations are distributed in three geographically separated regions: the NE (northeast coast), SE (southeast coast), and a limited area in the ND (Nakdong River). In this study, we used mitochondrial loci and microsatellites to investigate the evolutionary history of P. kaibarae populations by assessing their pattern of genetic structure. Our analyses revealed a marked level of divergence among three regional populations, suggesting a long history of isolation following colonization, although ND individuals showed relatively higher genetic affinity to populations from SE than those from NE. The populations from NE showed a great degree of interpopulation differentiation, whereas populations from SE exhibited only weak genetic structuring. Upon robust phylogenetic analysis, P. kaibarae formed a monophyletic group with Russian P. sinensis and P. tymensis with strong node confidence values, indicating that P. kaibarae populations on the Korean Peninsula originated from the southward migration of its ancestral lineage around the middle Pleistocene.