Large Scale Geographical Mapping of Essential Oil Volatiles in Heracleum (Apiaceae): Identification of Novel Compounds and Unraveling Cryptic Variation

Conspecific populations of plants in their native range are expected to show considerable variation due to long‐term ecological and evolutionary factors. We investigated the levels of secondary metabolites in Heracleum including H. persicum a valuable medicinal plant to depict the magnitude of cryptic variation and the potential significance of novel chemical traits. The essential oil volatiles from fruits of 34 populations from different species of Heracleum in Iranian distribution range and a native of H. sphondylium and an invasive population of H. persicum from Norway were analyzed with GC/MS. Out of 48 compounds identified, a contrasting pattern in the level of two major compounds, octyl acetate and hexyl butyrate was found among all studied species. Interestingly, a significant geographic pattern was observed; the hexyl butyrate/octyl acetate ratio was high (range 1.8 – 3.2) in the northwestern Iranian populations of H. persicum compared to that in northern and central populations (range 0.3 – 0.9). Four populations from Zagros mountains also exhibited a unique composition. Anethole was found in two populations of H. persicum from central Zagros, which has not been previously reported for essential oil of fruits of Heracleum so far. The results suggest high efficiency of large scale sampling from distribution range of species in identifying novel compounds. The unique pattern of geographic structuring also provides novel information to unravel cryptic variation in Heracleum.     

Geographic characteristics of sable (Martes zibellina) distribution over time in Northeast China

Understanding historical context can help clarify the ecological and biogeographic characteristics of species population changes. The sable (Martes zibellina ) population has decreased dramatically in Northeast China since the l950s, and understanding the changes in its distribution over time is necessary to support conservation efforts. To achieve this goal, we integrated ecological niche modeling and historical records of sables to estimate the magnitude of change in their distribution over time. Our results revealed a 51.71% reduction in their distribution in 2000–2016 compared with the potential distribution in the 1950s. This reduction was related to climate change (Pearson's correlation: Bio1, ?.962, p  < .01; Bio2, ?.962, p  < .01; Bio5, .817, p  < .05; Bio6, .847, p  < .05) and human population size (?.956, p  < .01). The sable population tended to migrate in different directions and elevations over time in different areas due to climate change: In the Greater Khingan Mountains, they moved northward and to lower elevations; in the Lesser Khingan Mountains, they moved northward; and in the Changbai Mountains, they move southward and to higher elevations. Active conservation strategies should be considered in locations where sable populations have migrated or may migrate to.     

Genetic variation in Tertiary relics: The case of eastern‐Mediterranean Abies (Pinaceae)

The eastern‐Mediterranean Abies taxa, which include both widely distributed species and taxa with minuscule ranges, represent a good model to study the impacts of range size and fragmentation on the levels of genetic diversity and differentiation. To assess the patterns of genetic diversity and phylogenetic relationships among eastern‐Mediterranean Abies taxa, genetic variation was assessed by eight nuclear microsatellite loci in 52 populations of Abies taxa with a focus on those distributed in Turkey and the Caucasus. Both at the population and the taxon level, the subspecies or regional populations of Abies nordmanniana s.l. exhibited generally higher allelic richness, private allelic richness, and expected heterozygosity compared with Abies cilicica s.l. Results of both the Structure analysis and distance‐based approaches showed a strong differentiation of the two A. cilicica subspecies from the rest as well as from each other, whereas the subspecies of A. nordmanniana were distinct but less differentiated. ABC simulations were run for a set of scenarios of phylogeny and past demographic changes. For A. ×olcayana, the simulation gave a poor support for the hypothesis of being a taxon resulting from a past hybridization, the same is true for Abies equi‐trojani: both they represent evolutionary branches of Abies bornmuelleriana.     

Widely Targeted Metabolomics Analysis of Dendrobium officinale at Different Altitudes

In order to investigate the differences in the chemical composition of Dendrobium officinale at different altitudes, we collected Dendrobium officinale from different altitudes in Guizhou Province, China, and firstly determined its polysaccharide content by using a sulfuric acid-phenol color development method with reference to the Chinese Pharmacopoeia, and further determined its metabolites by using widely targeted metabolomics, and explored the differences in the chemical composition of Dendrobium officinale at different altitudes using multivariate statistical analysis. We found that the polysaccharide content was higher in plants growing at 1122 m, a total of 902 secondary metabolites were detected by wildly targeted metabolomics, and amino acids and their derivatives were more highly expressed at 1122 m, while other metabolites were more highly expressed at 835 m. Furthermore, we found that the phenolic acid compound nerugein was only present in plants at 835 m, and two lipid compounds, namely Lyso PE 20 : 4 and its isomer, were only present at 1122 m. Taken together, these results may provide a basis for the selection and clinical application of D. officinale cultivated at different altitudes.     

Diversification in the South American Pampas: the genetic and morphological variation of the widespread Petunia axillaris complex (Solanaceae)

Understanding the spatiotemporal distribution of genetic variation and the ways in which this distribution is connected to the ecological context of natural populations is fundamental for understanding the nature and mode of intraspecific and, ultimately, interspecific differentiation. The Petunia axillaris complex is endemic to the grasslands of southern South America and includes three subspecies: P. a. axillaris, P. a. parodii and P. a. subandina. These subspecies are traditionally delimited based on both geography and floral morphology, although the latter is highly variable. Here, we determined the patterns of genetic (nuclear and cpDNA), morphological and ecological (bioclimatic) variation of a large number of P. axillaris populations and found that they are mostly coincident with subspecies delimitation. The nuclear data suggest that the subspecies are likely independent evolutionary units, and their morphological differences may be associated with local adaptations to diverse climatic and/or edaphic conditions and population isolation. The demographic dynamics over time estimated by skyline plot analyses showed different patterns for each subspecies in the last 100 000 years, which is compatible with a divergence time between 35 000 and 107 000 years ago between P. a. axillaris and P. a. parodii, as estimated with the IMa program. Coalescent simulation tests using Approximate Bayesian Computation do not support previous suggestions of extensive gene flow between P. a. axillaris and P. a. parodii in their contact zone.     

Discovering variation of secondary metabolite diversity and its relationship with disease resistance in Cornus florida L.

Understanding intraspecific relationships between genetic and functional diversity is a major goal in the field of evolutionary biology and is important for conserving biodiversity. Linking intraspecific molecular patterns of plants to ecological pressures and trait variation remains difficult due to environment‐driven plasticity. Next‐generation sequencing, untargeted liquid chromatography–mass spectrometry (LC‐MS) profiling, and interdisciplinary approaches integrating population genomics, metabolomics, and community ecology permit novel strategies to tackle this problem. We analyzed six natural populations of the disease‐threatened Cornus florida L. from distinct ecological regions using genotype‐by‐sequencing markers and LC‐MS‐based untargeted metabolite profiling. We tested the hypothesis that higher genetic diversity in C. florida yielded higher chemical diversity and less disease susceptibility (screening hypothesis), and we also determined whether genetically similar subpopulations were similar in chemical composition. Most importantly, we identified metabolites that were associated with candidate loci or were predictive biomarkers of healthy or diseased plants after controlling for environment. Subpopulation clustering patterns based on genetic or chemical distances were largely congruent. While differences in genetic diversity were small among subpopulations, we did observe notable similarities in patterns between subpopulation averages of rarefied‐allelic and chemical richness. More specifically, we found that the most abundant compound of a correlated group of putative terpenoid glycosides and derivatives was correlated with tree health when considering chemodiversity. Random forest biomarker and genomewide association tests suggested that this putative iridoid glucoside and other closely associated chemical features were correlated to SNPs under selection.     

Influence of altitude on local adaptation in upland tree species from central Argentina

• Steep climatic gradients boost morphological and physiological adjustments in plants, with consequences on performance. The three principal woody species of the Sierras Grandes Mountains of central Argentina have marked differences in sapling performance along their altitudinal distribution. We hypothesize that the steep gradient of climatic conditions across the species’ altitudinal distribution promotes trait differences between populations of different altitudes that are inherited by the following generation.
• Seeds from different altitudes were exposed to three temperature regimes to assess differential germination responses. Saplings were then transplanted to a greenhouse to assess possible variations in attributes and performance after 18 months.
• The three species showed differences in germination responses to temperature among altitudes and/or in sapling attributes and performance. In Maytenus boaria and Escallonia cordobensis, germination success was higher under high temperatures for the highest‐altitude, whereas lower temperatures boosted germination of the lowest altitudes. Polylepis australis showed no differences in germination among temperature treatments. In the greenhouse, saplings of the three species from intermediate altitudes showed high performance, whereas the upper and lower populations seemed to be adjusted to tolerating more stressful conditions (i.e., lower temperatures at the upper end and water stress at the lower end), showing lower performance toward both altitudinal limits.
• These patterns agree with those described for saplings growing under field conditions, suggesting adjustments in response to environmental changes undergone by populations along the altitudinal range. The marked adjustments of populations to the local environment suggest a potentially high impact of climatic change on species distribution.

     

Reproductive biology of Syzygiella rubricaulis (Nees) Steph. (Adelanthaceae,Marchantiophyta), a liverwort disjunctly distributed in high‐altitude Neotropical mountains

Syzygiella rubricaulis is a dioecious leafy liverwort disjunctly distributed and restricted to high‐altitude mountains in the Neotropics and the Azores. This study is part of a larger project examining the phylogeography of S. rubricaulis in the Neotropics, and our main goals were to understand its reproductive biology, where sex expression occurs, if vegetative propagules are frequently found, how the sexes are distributed in populations, how frequently sporophytes are formed and what environmental conditions influence sexual expression. S. rubricaulis patches are mostly female, but all patches also contain non sex‐expressing shoots. Out of 42 patches examined, 29 (69%) were sex‐expressing: 25 were unisexual (21 female and four male) and four of mixed sex (two male‐biased and two unbiased). At shoot level, out of 4200 shoots 18% were female and 7% male; among sex‐expressing shoots, 73% were female, representing a sex ratio of 0.8 (female‐biased). We encountered a total of 33 sporophytes in six patches (in Brazil, Venezuela and Ecuador). Leaf regenerants were found in one patch in Mexico. Low rates of sporophytes were likely related to low frequencies of male shoots and large distances between the sexes. As 25% of S. rubricaulis shoots expressed sex (occasionally producing sporophytes), we suggest that short‐distance (and rarely long‐distance) spore dispersal events occur in mountainous areas on a short‐term basis. On a long‐term basis, however, these events likely contribute to dynamic exchanges among populations in the Neotropics.     

Thermal ecology of allopatric lizards (Sphenomorphus) in southeast Australia

Summary The Sphenomorphus quoyi species complex (Lacertilia: Scincidae) of southeast Australia consists of three lizard species, distributed as continuous but largely allopatric populations, and a fourth species which has a disjunct distribution at high altitudes within the range of the other three. The four species have similar Voluntary temperatures (range of body temperatures within which normal activity occurs) despite the fact that they occupy different geographical areas.Altitudinal zonation is shown in the distribution pattern of the four species with S. quoyi and S. tympanum (warm temperate) restricted to the lower altitudes. The level of the mean daily minimum and maximum temperatures are different in each of the four species distribution areas. Similarly, the area occupied by each species is characterized by different daily sunlit and shade surface temperatures and different duration of ground surface subzero temperatures. Floristic composition of the four areas and the basking sites utilized by each species are also dissimilar.The lower temperature tolerances and the acclimation range of the Critical Minimum temperatures of S. quoyi and S. tympanum (warm temperate) exclude these species from the higher altitudes with the relatively lower environmental temperature levels. The Critical temperatures of the species which are found at high altitudes, do not explain why sympatry is not found with the lowland species.     

Geographic Distribution of Natural Products Produced by the Red Alga Laurencia dendroidea J. Agardh

In order to evaluate the chemical diversity of Laurencia dendroidea J. Agardh , a widely distributed seaweed in Brazilian coast, a phytochemical study was carried out with algae collected from six different locations along the Southeast Brazilian coast. Purified compounds were identified by MS and NMR techniques. The chemical profiles of lipophilic extracts were obtained by GC/MS for each population. In total, 15 compounds were described. The sesquiterpene composition accounted for 49 – 63% of the GC/MS chromatogram area. The discrimination of three chemotypes was done by the use of HCA on GC/MS chromatograms. They were also analyzed by the PCA and, together with peak area analysis, it was possible to discriminate all populations by the main variation of elatol, obtusol, rogiolol, and triquinane. The results revealed the high diversity of sesquiterpene composition among populations of L. dendroidea. Curiously, the within and among population variation of elatol and obtusol suggested a biochemical interplay on the content of these compounds. More studies are necessary to understand the patterns of chemical diversity and compound variation within and among populations of L. dendroidea.     

Genome‐wide SNP analysis unveils genetic structure and phylogeographic history of snow sheep (Ovis nivicola) populations inhabiting the Verkhoyansk Mountains and Momsky Ridge (northeastern Siberia)

Insights into the genetic characteristics of a species provide important information for wildlife conservation programs. Here, we used the OvineSNP50 BeadChip developed for domestic sheep to examine population structure and evaluate genetic diversity of snow sheep (Ovis nivicola) inhabiting Verkhoyansk Range and Momsky Ridge. A total of 1,121 polymorphic SNPs were used to test 80 specimens representing five populations, including four populations of the Verkhoyansk Mountain chain: Kharaulakh Ridge–Tiksi Bay (TIK, n = 22), Orulgan Ridge (ORU, n = 22), the central part of Verkhoyansk Range (VER, n = 15), Suntar‐Khayata Ridge (SKH, n = 13), and Momsky Ridge (MOM, n = 8). We showed that the studied populations were genetically structured according to a geographic pattern. Pairwise F_ST values ranged from 0.044 to 0.205. Admixture analysis identified K = 2 as the most likely number of ancestral populations. A Neighbor‐Net tree showed that TIK was an isolated group related to the main network through ORU. TreeMix analysis revealed that TIK and MOM originated from two different ancestral populations and detected gene flow from MOM to ORU. This was supported by the f3 statistic, which showed that ORU is an admixed population with TIK and MOM/SKH heritage. Genetic diversity in the studied groups was increasing southward. Minimum values of observed (Ho) and expected (He) heterozygosity and allelic richness (Ar) were observed in the most northern population—TIK, and maximum values were observed in the most southern population—SKH. Thus, our results revealed clear genetic structure in the studied populations of snow sheep and showed that TIK has a different origin from MOM, SKH, and VER even though they are conventionally considered a single subspecies known as Yakut snow sheep (Ovis nivicola lydekkeri). Most likely, TIK was an isolated group during the Late Pleistocene glaciations of Verkhoyansk Range.     

Glucosinolate polymorphism in wild cabbage (Brassica oleracea) influences the structure of herbivore communities

Natural plant populations often show substantial heritable variation in chemical structure of secondary metabolites. Despite a great deal of evidence from laboratory studies that these chemicals influence herbivore behaviour and life history, there exists little evidence for the structuring of natural herbivore communities according to plant chemical profiles. Brassica oleracea (Brassicaceae) produces aliphatic glucosinolates, which break down into toxins when leaf tissue is damaged. Structural diversity in these glucosinolates is heritable, and varies considerably at two ecological scales in the UK: both within and between populations. We surveyed herbivore attack on plants producing different glucosinolates, using 12 natural B. oleracea populations. In contrast to the results of previous studies in this system, which suffered low statistical power, we found significant differential responses of herbivore species to heritable glucosinolates, both within and between plant populations. We found significant correlations between herbivore infestation rates and the presence or absence of two heritable glucosinolates: sinigrin and progoitrin. There was variation between herbivore species in the direction of response, the ecological scale at which responses were identified, and the correlations for some herbivore species changed at different times of the year. We conclude that variation in plant secondary metabolites can structure the community of herbivores that attack them, and propose that herbivore-mediated differential selection deserves further investigation as a mechanism maintaining the observed diversity of glucosinolates in wild Brassica.     

Altitudinal distribution of chytrid (Batrachochytrium dendrobatidis) infection in subtropical Australian frogs

The disappearance of amphibian populations from seemingly pristine upland areas worldwide has become a major focus of conservation efforts in the last two decades, and a parasitic chytrid fungus, Batrachochytrium dendrobatidis, is thought to be the causative agent of the population declines. We examined the altitudinal distribution of chytrid infections in three stream‐dwelling frog species (Litoria wilcoxii, L. pearsoniana and L. chloris) in southeast Queensland, Australia, and hypothesized that if B. dendrobatidis were responsible for the disappearance of high‐altitude frog populations, infection prevalence and intensity would be greatest at higher altitudes. Overall, 37.7% of the 798 adult frogs we sampled were infected with B. dendrobatidis, and infections were found in all the populations we examined. Contrary to our initial hypothesis, we found no consistent evidence that high‐altitude frogs were more likely to be infected than were lowland frogs. Further, the intensity of fungal infections (number of zoospores) on high‐altitude frogs did not differ significantly from that of lowland frogs. Batrachochytrium dendrobatidis appears to be capable of infecting frogs at all altitudes in the subtropics, suggesting that all populations are at risk of decline when conditions favour disease outbreaks. We did find evidence, however, that chytrid infections persist longer into summer in upland as compared with lowland areas, suggesting that montane amphibian populations remain susceptible to disease outbreaks for longer periods than do lowland populations. Further, we found that at high altitudes, temperatures optimal for chytrid growth and reproduction coincide with frog metamorphosis, the life‐stage at which frogs are most susceptible to chytrid infections. While these altitudinal differences may account for the differential population‐level responses to the presence of B. dendrobatidis, the reason why many of southeast Queensland's montane frog populations declined precipitously while lowland populations remained stable has yet to be resolved.     

Influence of northern limit range on genetic diversity and structure in a widespread North American tree,sugar maple (Acer saccharum Marshall)

Due to climate change, the ranges of many North American tree species are expected to shift northward. Sugar maple (Acer saccharum Marshall) reaches its northern continuous distributional limit in northeastern North America at the transition between boreal mixed‐wood and temperate deciduous forests. We hypothesized that marginal fragmented northern populations from the boreal mixed wood would have a distinct pattern of genetic structure and diversity. We analyzed variation at 18 microsatellite loci from 23 populations distributed along three latitudinal transects (west, central, and east) that encompass the continuous–discontinuous species range. Each transect was divided into two zones, continuous (temperate deciduous) and discontinuous (boreal mixed wood), based on sugar maple stand abundance. Respective positive and negative relationships were found between the distance of each population to the northern limit (D_north), and allelic richness (A_R) and population differentiation (F_ST). These relations were tested for each transect separately; the pattern (discontinuous–continuous) remained significant only for the western transect. structure analysis revealed the presence of four clusters. The most northern populations of each transect were assigned to a distinct group. Asymmetrical gene flow occurred from the southern into the four northernmost populations. Southern populations in Québec may have originated from two different postglacial migration routes. No evidence was found to validate the hypothesis that northern populations were remnants of a larger population that had migrated further north of the species range after the retreat of the ice sheet. The northernmost sugar maple populations possibly originated from long‐distance dispersal.     

Spatiotemporal and demographic variation in the diet of New Zealand lesser short‐tailed bats (Mystacina tuberculata)

Variation in the diet of generalist insectivores can be affected by site‐specific traits including weather, habitat, and season, as well as demographic traits such as reproductive status and age. We used molecular methods to compare diets of three distinct New Zealand populations of lesser short‐tailed bats, Mystacina tuberculata. Summer diets were compared between a southern cold‐temperate (Eglinton) and a northern population (Puroera). Winter diets were compared between Pureora and a subtropical offshore island population (Hauturu). This also permitted seasonal diet comparisons within the Pureora population. Lepidoptera and Diptera accounted for >80% of MOTUs identified from fecal matter at each site/season. The proportion of orders represented within prey and the Simpson diversity index, differed between sites and seasons within the Pureora population. For the Pureora population, the value of the Simpson diversity index was higher in summer than winter and was higher in Pureora compared to Eglinton. Summer Eglinton samples revealed that juvenile diets appeared to be more diverse than other demographic groups. Lactating females had the lowest dietary diversity during summer in Pureora. In Hauturu, we found a significant negative relationship between mean ambient temperature and prey richness. Our data suggest that M. tuberculata incorporate a narrower diversity of terrestrial insects than previously reported. This provides novel insights into foraging behavior and ecological interactions within different habitats. Our study is the first from the Southern Hemisphere to use molecular techniques to examine spatiotemporal variation in the diet of a generalist insectivore that inhabits a contiguous range with several habitat types and climates.     

Demographic history and the low genetic diversity in Dipteryx alata (Fabaceae) from Brazilian Neotropical savannas

Genetic effects of habitat fragmentation may be undetectable because they are generally a recent event in evolutionary time or because of confounding effects such as historical bottlenecks and historical changes in species'' distribution. To assess the effects of demographic history on the genetic diversity and population structure in the Neotropical tree Dipteryx alata (Fabaceae), we used coalescence analyses coupled with ecological niche modeling to hindcast its distribution over the last 21 000 years. Twenty-five populations (644 individuals) were sampled and all individuals were genotyped using eight microsatellite loci. All populations presented low allelic richness and genetic diversity. The estimated effective population size was small in all populations and gene flow was negligible among most. We also found a significant signal of demographic reduction in most cases. Genetic differentiation among populations was significantly correlated with geographical distance. Allelic richness showed a spatial cline pattern in relation to the species'' paleodistribution 21 kyr BP (thousand years before present), as expected under a range expansion model. Our results show strong evidences that genetic diversity in D. alata is the outcome of the historical changes in species distribution during the late Pleistocene. Because of this historically low effective population size and the low genetic diversity, recent fragmentation of the Cerrado biome may increase population differentiation, causing population decline and compromising long-term persistence.     

The Phenolics of the Ornithogalum umbellatum L. (Hyacinthaceae): Phytochemical and Ecological Characterization

Ornithogalum umbellatum L. is a widely distributed species in Europe, that exhibits considerable variability at ecological, morphological, anatomical, and karyological level. Previous reports of the chemical investigations among Ornithogalum species indicate significant diversity of the secondary metabolites, as well. Knowing that environment affects the phenolic composition in plants to a large extent, the main objective of the research was to define relationship between phytochemical and ecological characters. To estimate an environmental influence on these results, plant material was collected at four habitats that differ in ecological factors and belong to two biogeographical regions: the Balkan Peninsula and the Pannonian Plane. Measured phytochemical characters are yield of dry extract, total phenolic and flavonoid contents, the presence of selected phenolic compounds as well as free‐radical scavenging activity (neutralization of DPPH and OH radicals). Results revealed that all analyzed phytochemical parameters differ between investigated O. umbellatum samples. The moisture level of habitat has the highest correlation, either positive or negative, with most of phytochemical characters, and is followed by temperature and soil reaction. Light intensity and nitrogen level have mostly moderate correlation coefficient with phytochemical characters. More complex correlation is revealed between ecological factors and nine phenolic compounds, with three observed patterns of relationship.     

Finding the appropriate variables to model the distribution of vector‐borne parasites with different environmental preferences: climate is not enough

Understanding how environmental variation influences the distribution of parasite diversity is critical if we are to anticipate disease emergence risks associated with global change. However, choosing the relevant variables for modelling current and future parasite distributions may be difficult: candidate predictors are many, and they seldom are statistically independent. This problem often leads to simplistic models of current and projected future parasite distributions, with climatic variables prioritized over potentially important landscape features or host population attributes. We studied avian blood parasites of the genera Plasmodium, Haemoproteus and Leucocytozoon (which are viewed as potential emergent pathogens) in 37 Iberian blackcap Sylvia atricapilla populations. We used Partial Least Squares regression to assess the relative importance of a wide array of putative determinants of variation in the diversity of these parasites, including climate, landscape features and host population migration. Both prevalence and richness of parasites were predominantly related to climate (an effect which was primarily, but not exclusively driven by variation in temperature), but landscape features and host migration also explained variation in parasite diversity. Remarkably, different models emerged for each parasite genus, although all parasites were studied in the same host species. Our results show that parasite distribution models, which are usually based on climatic variables alone, improve by including other types of predictors. Moreover, closely related parasites may show different relationships to the same environmental influences (both in magnitude and direction). Thus, a model used to develop one parasite distribution can probably not be applied identically even to the most similar host–parasite systems.     

Comparative phylogeographic study of Hosta sieboldiana and Hosta albomarginata (Asparagaceae) in Japan

We analyzed variations in chloroplast DNA (cpDNA) in the widespread herbaceous species Hosta sieboldiana and Hosta albomarginata across large portions of their geographic ranges in the Japanese archipelago. Our objective was to compare the phylogeographic histories and phylogeographic structures of the two congeneric species in the Japanese archipelago. The location of the study is Japanese archipelago. We sequenced 1380 bp of noncoding cpDNA from 45 populations of H. sieboldiana (n = 362) and 55 populations of H. albomarginata (n = 436) to assess genetic variations within and among populations across almost the entire distributions of the species in Japan. Extant patterns of geographic structure were analyzed using statistical parsimony networks and spatial analysis of molecular variance (SAMOVA). We also used Monmonier's algorithm to detect genetic barriers between regions. Relationships between the populations were examined using a neighbor‐joining (NJ) method. Four haplotypes were found for H. sieboldiana, whereas eight haplotypes were identified for H. albomarginata. Total genetic haplotype diversity (hT) and within‐population haplotype diversity (hS) for H. sieboldiana were 0.352 and 0.040, respectively, while the values for H. albomarginata were 0.529 and 0.085, respectively. The population differentiations (GST) for H. sieboldiana and H. albomarginata were 0.839 and 0.886, respectively. The SAMOVA analysis revealed two clusters in H. sieboldiana and four clusters in H. albomarginata. Differentiations between and among the clusters were supported by the BARRIER analysis and the NJ tree. We detected differences in the population genetic structure between the two species. We found that H. sieboldiana had lower haplotype diversity than H. albomarginata. These results may be partially explained by the difference in ecological habitats and geographic distributions between the species. Hosta albomarginata is more widely distributed than H. sieboldiana in East Asia including Russia, and this large distribution range would enable more chances to intraspecific gene flow.