Modern State of Populations of Endemic <Emphasis Type="Italic">Oxytropis</Emphasis> Species from Baikal Siberia and Their Phylogenetic Relationships Based on Chloroplast DNA Markers

The genetic diversity and population structure of the endemic species of Baikal Siberia Oxytropis triphylla, O. bargusinensis, and O. interposita were studied for the first time on the basis of the nucleotide polymorphism of intergenic spacers psbA–trnH, trnL–trnF, and trnS–trnG of chloroplast DNA. All populations of these species were characterized by a high haplotype (0.762–0.924) and relatively low nucleotide (0.0011–0.0022) diversity. Analysis of the distribution of variability in O. triphylla and O. bargusinensis showed that there was no significant genetic differentiation between populations of each species; the gene flow was 4.43 and 8.91, respectively. The high level of genetic diversity in the studied populations indicates a relatively stable state of these populations. A study of the phylogenetic relationships of closely related species confirms the concept of the origin of O. bargusinensis and O. tompudae as a result of intersectional hybridization of the species of the sections Orobia and Verticillares.     

Pollen-mediated gene flow promotes low nuclear genetic differentiation among populations of <Emphasis Type="Italic">Cycas debaoensis</Emphasis> (Cycadaceae)

Cycas debaoensis is a critically endangered cycad species endemic to China. This species is found on two kinds of habitats according to the edaphic differences, sand and karst. A previous chloroplast DNA (cpDNA) study indicated that C. debaoensis had low genetic variation within populations and high genetic differentiation among populations. Because maternally inherited cpDNA does not fully characterize genetic structure of the species, we screened seven low-copy nuclear genes and 17 nuclear microsatellite loci to detect the nuclear genetic diversity, differentiation, and the population structure of C. debaoensis. The nuclear genes revealed higher level of genetic diversity. There were both the same and region-specific haplotypes or alleles between the karst and sand regions. Nuclear gene flow among all the populations was much greater than that of cpDNA, which indicated that pollen-mediated gene flow was much greater than seed-mediated gene flow. This promoted low nuclear genetic differentiation among populations of C. debaoensis. The study suggests that both genetic and anthropogenic disturbances have resulted in the critically endangered status of C. debaoensis.     

Molecular genetic characteristic of dinucleotide microsatellite loci in parthenogenetic lizards <Emphasis Type="Italic">Darevskia unisexualis</Emphasis>

In the present study, the first molecular genetic investigation of dinucleotide (GT) _n microsatellite loci in parthenogenetic lizards Darevskia unisexualis was performed. New polymorphic locus, Du214, (GenBank Ac. No. EU252542) was identified and characterized in detail. It was demonstrated that allele of this locus differed in the size and structure of microsatellite locus, as well as in point mutations, the combinations of which enabled the isolation of stabile fixed double nucleotide substitutions A-A (alleles 2 and 4) and G-T (alleles 1, 3, 5, and 6). Double nucleotide substitutions described were also identified in the orthlogous loci of the parental species genomes, D. raddei (G-T) and D. valentine (A-A). Based on the analysis of allele distribution pattern at this locus in all populations of parthenospecies D. unisexualis, mathematic model was elaborated and realized. Using this model, frequencies of allelic variants for all populations of the species of interest were calculated and population genetic structure of D. unisexualis was characterized. Genetic contribution of each population to the species gene pool was determined. The data obtained demonstrated that microsatellite variation was one of the factors of clonal and genetic diversity of a parthenospecies.     

Isolation and Characterization of Microsatellite Markers for <Emphasis Type="Italic">Viola websteri</Emphasis> (Violaceae) and Cross-Species Amplification within the Genus <Emphasis Type="Italic">Viola</Emphasis>

We isolated and characterized microsatellite loci in Viola websteri (Violaceae), an endangered species from Korea and endemic to Northeast Asia. A total of 27 microsatellite loci were developed and tested in Korean and Chinese populations. The number of alleles per locus varied from two to eight. The observed and expected heterozygosities within two populations were 0.000 to 1.000 and 0.080 to 0.816, respectively. Korean and Chinese populations were clearly distinguished by the private alleles from 16 loci. A total of 21 loci out of the 27 developed loci were successfully cross-amplified in 39 other Viola species. We believe that these microsatellite loci will be useful for future studies on genetic diversity and population structure of V. websteri, as well as other Viola species.     

Limited diversity associated with duplicated class II MHC-<Emphasis Type="Italic">DRB</Emphasis> genes in the red squirrel population in the United Kingdom compared with continental Europe

The red squirrel (Sciurus vulgaris) population in the United Kingdom has declined over the last century and is now on the UK endangered species list. This is the result of competition from the eastern grey squirrel (S. carolinensis) which was introduced in the 19th century. However, recent evidence suggests that the rate of population decline is enhanced by squirrelpox disease, caused by a viral infection carried asymptomatically by grey squirrels but to which red squirrels are highly susceptible. Population genetic diversity provides some resilience to rapidly evolving or exotic pathogens. There is currently no data on genetic diversity of extant UK squirrel populations with respect to genes involved in disease resistance. Diversity is highest at loci involved in the immune response including genes clustered within the major histocompatibility complex (MHC). Using the class II DRB locus as a marker for diversity across the MHC region we genotyped 110 red squirrels from locations in the UK and continental Europe. Twenty-four Scvu-DRB alleles at two functional loci; Scvu-DRB1 and Scvu-DRB2, were identified. High levels of diversity were identified at both loci in the continental populations. In contrast, no diversity was observed at the Scvu-DRB2 locus in the mainland UK population while a high level of homozygosity was observed at the Scvu-DRB1 locus. The red squirrel population in the UK appears to lack the extensive MHC diversity associated with continental populations, a feature which may have contributed to their rapid decline.     

Population genetic diversity and geographical differentiation of MHC class II DAB genes in the vulnerable Chinese egret (<Emphasis Type="Italic">Egretta eulophotes</Emphasis>)

Major histocompatibility complex (MHC) genes are excellent markers for the study of adaptive genetic variation occurring over different geographical scales. The Chinese egret (Egretta eulophotes) is a vulnerable ardeid species with an estimated global population of 2600–3400 individuals. In this study, we sampled 172 individuals of this egret (approximately 6 % of the global population) from five natural populations that span the entire distribution range of this species in China. We examined their population genetic diversity and geographical differentiation at three MHC class II DAB genes by identifying eight exon 2 alleles at Egeu-DAB1, eight at Egeu-DAB2 and four at Egeu-DAB3. Allelic distributions at each of these three Egeu-DAB loci varied substantially within the five populations, while levels of genetic diversity varied slightly among the populations. Analysis of molecular variance showed low but significant genetic differentiation among five populations at all three Egeu-DAB loci (haplotype-based ?_ST: 0.029, 0.020 and 0.042; and distance-based ?_ST: 0.036, 0.027 and 0.043, respectively; all P < 0.01). The Mantel test suggested that this significant population genetic differentiation was likely due to an isolation-by-distance pattern of MHC evolution. However, the phylogenetic analyses and the Bayesian clustering analysis based on the three Egeu-DAB loci indicated that there was little geographical structuring of the genetic differentiation among five populations. These results provide fundamental population information for the conservation genetics of the vulnerable Chinese egret.     

Molecular evolution of glutathione peroxidase (<Emphasis Type="Italic">Gpx</Emphasis>) family genes in desert poplar (<Emphasis Type="Italic">Populus euphratica</Emphasis> Oliv.)

Populus euphratica Oliv. is a poplar species that is distributed mainly in deserts, making it an interesting model in which to investigate molecular mechanisms underlying different stress responses. Here, we used molecular population genetic methods to detect potential selection in candidate genes belonging to the P. euphratica glutathione (GSH) peroxidase (Gpx) family, which are associated with an enzymatic mechanism that combats oxidative damage caused by reactive oxygen species (ROS) in plant cells; earlier studies have shown that Gpx proteins play important roles in coping with increased ROS levels during biotic and abiotic stresses in plants. We analyzed the nucleotide diversity and divergence patterns of five loci encoding Gpx genes, and 16 reference loci used as controls, in order to detect departures from the neutral expectation. Gpx1 has an excess of mid-frequency alleles; high intraspecific nucleotide diversity, distributed in the upper tail of the simulated neutral model; and extensive LD, showing strong evidence of balancing selection/local adaptation. The Gpx3.2 gene exhibits very low nucleotide diversity and divergence, suggesting that it has evolved under strong purifying selection. We failed to detect any evidence for natural selection at the other loci (Gpx2, Gpx4, and Gpx5) compared with the reference loci. The results show that nucleotide diversity and/or divergence differ greatly between members of the Gpx gene family, resulting from differential selective pressure acting on genes, and that adaptive evolution influenced the distribution of P. euphratica in desert regions.     

Genetic diversity and differentiation of Siberian spruce populations at nuclear microsatellite loci

The results of the study of 21 populations of Siberian spruce (Picea obovata Ledeb.) from different parts of the species natural range by microsatellite (SSR) analysis of nuclear DNA are presented. Using nine loci developed for Picea abies (L.) Karst. and Picea glauca (Moench) Voss and detecting variation in Picea obovata, the parameters of intra- and interpopulation genetic diversity, as well as the degree of population differentiation, were determined. It was demonstrated that the population of Siberian spruce in the study was characterized by a relatively high average level of intrapopulation variability (H_o = 0.408; H_e = 0.423) and low interpopulation differentiation (F_st = 0.048, P = 0.001) at this class of DNA markers. The genetic distance between populations ranged from 0.009 to 0.167, averaging 0.039. The isolated Magadan population, located in the extreme Northeast of Russia at a considerable distance from the main species range and characterized by the lowest genetic diversity among the studied populations, was maximally differentiated from the rest of the spruce populations. In addition, the steppe Ubukun population from Buryatia and the population from the Bogd Khan Uul Biosphere Reserve, Mongolia, were considerably different in the genetic structure from most populations of Siberian spruce, although to a lesser extent than the Magadan population. These findings are consistent with the results of previous studies of this species carried out using allozyme and microsatellite loci of chloroplast DNA and point to the prospects of using nuclear microsatellites as DNA markers to analyze the population genetic structure of Siberian spruce.     

Genetic Diversity and Phylogenetic Relationships of Two Closely Related Northeast China <Emphasis Type="Italic">Vicia</Emphasis> Species Revealed with RAPD and ISSR Markers

RAPD and ISSR analyses revealed genetic diversity and relationships among 11 populations of two closely related northeast China Vicia species, Vicia ramuliflora and V. unijuga. Both methods yielded similar and complementary results, showing high genetic diversity. Vicia ramuliflora had 100% polymorphic loci in both RAPD and ISSR, and V. unijuga had 100% polymorphic loci for RAPD and 98.96% for ISSR. Genetic differentiation was moderate among populations of each species. Genetic variation was distributed mainly within populations for the two species. The high level of gene flow was important for the allocation of genetic variation. The UPGMA dendrogram and principal coordinates analysis at the level of individuals and populations showed that V. ramuliflora and V. unijuga were more closely related than either of them was to the outgroup species, V. cracca. The small molecular variance of V. ramuliflora and V. unijuga supports the conclusion that these two species had a common ancestor.     

Nucleotide polymorphisms of candidate genes of adaptive significance in the ural populations of <Emphasis Type="Italic">Larix sibirica</Emphasis> Ledeb.

The objectives of conservation and sustainable forest management require in depth study of genomes of woody plants and definition of their intraspecific genetic diversity. In recent years, an approach was developed based on the study of “candidate genes” that can potentially be involved in the formation of adaptive traits. In this study, we investigated nucleotide polymorphism of several adaptive candidate genes in the populations of Siberian larch (Larix sibirica Ledeb.) in the Urals. Representatives of this genus are among the most valuable and widely distributed forest tree species in Russia. From ten selected gene loci in the genome of L. sibirica, we isolated and investigated three loci, one of which (ABA-WDS) was sequenced in L. sibirica for the first time. The total length of the analyzed sequence in each individual amounted to 2865 bp. The length of locus alignment was from 360 bp to 1395 bp. In total, we identified 200 polymorphic positions. The most conservative is locus 4CL1-363, and the most polymorphic is locus sSPcDFD040B03103-274. The studied populations of L. sibirica are characterized by a high level of nucleotide polymorphism in comparison with other species and genuses (Picea, Pinus, Pseudotsuga, Abies) conifers plants (Hd = 0.896; π = 0.007; θ_W = 0.015). The most selectively neutral polymorphism (D _T =–0.997) was attributed to locus 4CL1-363, and polymorphism with high probability of adaptability (D _T =–1.807) was determined for the ABA-WDS locus. We identified 54 SNP markers, only five of which were nonsynonymous (9.26%) replacements. The average frequency of SNPs in the three studied loci of L. sibirica was one SNP in 53 bp. We detected unique SNP markers for eight populations, which could potentially be used to identify populations. Populations that are characterized by the highest number of unique SNP markers can be recommended for selection in order to preserve the gene pool of the species.     

Genetic diversity and population structure in the narrow endemic Chinese walnut <Emphasis Type="Italic">Juglans hopeiensis</Emphasis> Hu: implications for conservation

The conservation of narrow endemic species relies on accurate information regarding their population structure. Juglans hopeiensis Hu (Ma walnut), found only in Hebei province, Beijing, and Tianjin, China, is a threatened tree species valued commercially for its nut and wood. Sequences of two maternally inherited mitochondrial markers and two maternally inherited chloroplast intergenic spacers, three nuclear DNA sequences, and allele sizes from 11 microsatellites were obtained from 108 individuals of J. hopeiensis, Juglans regia, and Juglans mandshurica. Haplotype networks were constructed using NETWORK. Genetic diversity, population differentiation, and analysis of molecular variance (AMOVA) were used to determine genetic structure. MEGA was used to construct phylogenetic trees. Genetic diversity of J. hopeiensis was moderate based on nuclear DNA, but low based on uniparentally inherited mitochondrial DNA and chloroplast DNA. Haplotype networks showed that J. hopeiensis haplotypes were different than haplotypes found in J. regia and J. mandshurica. Allelic variants in nuclear genes that were shared among J. hopeiensis populations were not found in J. regia or J. mandshurica. Sampled populations of J. hopeiensis showed clear genetic structure. The maximum parsimony (MP) tree showed J. hopeiensis to be distinct from J. mandshurica but threatened by hybridization with J. regia and J. mandshurica. J. hopeiensis populations are strongly differentiated from sympatric Juglans species, but they are threatened by small population sizes and hybridization.     

Genetic structure and differentiation of populations of the tetraploid species <Emphasis Type="Italic">Oxytropis chankaensis</Emphasis> (Fabaceae)

The population genetic variation of the tetraploid species Oxytropis chankaensis Jurtz. (Fabaceae), a local endemic of the western coast of Khanka Lake (Primorye), was examined. Five populations were analyzed using 28 isozyme loci encoding 16 enzyme systems. Significant allelic heterogeneity among the populations was found for six out of twelve polymorphic loci. The heterozygosity of the samples (total sample size 294 plants) H _e = 0.301 was considerable higher than the mean values in populations of endemic species (0.076). Based on the results of this study, we identified two groups of O. chankaensis populations (southern and northern), in spite of the absence of marked hiatus between them. Of special interest is the population from Przhewalski Spit, which is a natural reserve of genetic diversity of the species and the putative center of formation of the autotetraploid O. chankaensis.     

Intraspecific Variation in Mitogenomes of Five <Emphasis Type="Italic">Crassostrea</Emphasis> Species Provides Insight into Oyster Diversification and Speciation

A large number of Crassostrea oysters are found in Asia-Pacific. While analyses of interspecific variation have helped to establish historical relationships among these species, studies on intraspecific variation are necessary to understand their recent evolutionary history and current forces driving population biology. We resequenced 18 and analyzed 31 mitogenomes of five Crassostrea species from China: Crassostrea gigas, Crassostrea angulata, Crassostrea sikamea, Crassostrea ariakensis, and Crassostrea hongkongensis. Our analysis finds abundant insertions, deletions, and single-nucleotide polymorphisms in all species. Intraspecific variation varies greatly among species with polymorphic sites ranging from 54 to 293 and nucleotide diversity ranging from 0.00106 to 0.00683. In all measurements, C. hongkongensis that has the narrowest geographic distribution exhibits the least sequence diversity; C. ariakensis that has the widest distribution shows the highest diversity, and species with intermediate distribution show intermediate levels of diversity. Low sequence diversity in C. hongkongensis may reflect recent bottlenecks that are probably exacerbated by human transplantation. High diversity in C. ariakensis is likely due to divergence of northern and southern China populations that have been separated without gene flow. The significant differences in mitogenome diversity suggest that the five sister species of Crassostrea have experienced different evolutionary forces since their divergence. The recent divergence of two C. ariakensis populations and the C. gigas/angulata species complex provides evidence for continued diversification and speciation of Crassostrea species along China’s coast, which are shaped by unknown mechanisms in a north–south divide.     

Genetic diversity and population structure of <Emphasis Type="Italic">Taxus cuspidata</Emphasis> Sieb. et Zucc. ex Endl. (Taxaceae) in Russia according to data of the nucleotide polymorphism of intergenic spacers of the chloroplast genome

The genetic diversity and population structure of Taxus cuspidata Sieb. et Zucc. ex Endl. on the Russian part of its range was studied for the first time on the basis of the nucleotide polymorphism of intergenic spacers psbA–trnH, trnL–trnF, and trnS–trnfM of chloroplast DNA. A high level of gene (h = 0.807) and nucleotide (π = 0.0227) diversity was revealed. The data of AMOVA showed that the interpopulation component accounted for 12% of genetic variability (F _ST = 0.12044, P = 0.0000). We revealed 15 haplotypes, four of which were shown to be unique. The presence of common haplotypes in the majority of populations, the absence of phylogenetic structure, and low values or even the absence of nucleotide divergence show that the T. cuspidata populations studied are fragments of the once common ancestor population. Geographical isolation, which resulted from climatic changes in the Pleistocene–Holocene, as well as from human activities, did not produce a significant effect on the genetic structure of the species.     

Meta-analysis reveals asymmetric reduction in the genetic diversity of introduced populations of exotic insects

Factors promoting the invasion success of introduced populations have been receiving increased attention in studies of biological invasions. Previous reports have indicated that successful invasions may be attributable to reduced genetic diversity in the invasive species. However, there is large variation in the magnitude and direction of the impact of exotic species that have remained unexplained. Here, we present a structured meta-analysis of papers investigating the genetic diversity of native and introduced populations of exotic insects using nuclear microsatellites and mitochondrial DNA sequences. The results indicate that invasion by exotic insects had an overall reducing effect on the genetic diversity of the invading population, with nonzero effect sizes for the number of alleles (NA), observed heterozygosity (Ho), expected heterozygosity (He) and nucleotide diversity (Nd). However, when analyzing different orders (e.g., Lepidoptera, Hemiptera), the effect sizes of NA, Ho and Nd in Lepidoptera were found to bracket zero, as did the effect size of He in Hemiptera. These results suggest an asymmetric reduction in the genetic diversity of introduced populations of exotic insects, indicating diverse mechanisms underlying their successful invasion.     

Novel polymorphic microsatellite loci for distinguishing rock bass (<Emphasis Type="Italic">Ambloplites rupestris</Emphasis>), Roanoke bass (<Emphasis Type="Italic">Ambloplites cavifrons</Emphasis>), and their hybrids

The rock bass (Ambloplites rupestris) is a popular sport-fish native to the Mississippi and Great Lakes basins of North America. The species has been widely introduced outside its native range, including into Atlantic-slope streams of Virginia where it may hybridize with an imperiled, similar-looking congener, the Roanoke bass (Ambloplites cavifrons). In this study, we identified and evaluated novel molecular markers to facilitate identification of these species and study the extent of hybridization. Using molecular libraries developed from A. rupestris, we identified a suite of candidate nuclear microsatellite loci, synthesized primer sets, and tested these markers for amplification and polymorphism in populations of both species. We then calculated standard diversity statistics within and differentiation statistics between species, the latter providing an indication of marker power for distinguishing the species and their hybrids. Additionally, we evaluated our efficiency for identifying hybrids by classifying simulated genotypes of known ancestry. Eleven loci were polymorphic (2–22 alleles per locus) and reliably amplified in both species. Multilocus genetic differentiation between A. cavifrons and A. rupestris was quite high (F _ST  = 0.66; D _LR  = 19.3), indicating the high statistical power of this marker set for species and hybrid identification. Analyses of simulated data suggested these markers reliably distinguish between hybrids and non-hybrids, as well as between F1 hybrids and backcrossed individuals. This panel of 11 loci should prove useful for understanding patterns of hybridization between A. rupestris and A. cavifrons. As the first microsatellite markers developed for Ambloplites, these markers also should prove broadly useful for population genetic studies of this genus.     

Analysis of Allelic Diversity and Genetic Relationships Among Cultivated Mangosteen (<Emphasis Type="Italic">Garcinia mangostana</Emphasis> L.) in Java,Indonesia Using Microsatellite Markers and Morphological Characters

The genetic variation and relationships of the mangosteen (Garcinia mangostana L.) were observed across mangosteen populations in Java, Indonesia using newly identified microsatellite loci and morphological characters. In this study, we developed an improved protocol to isolate microsatellite loci, named Selective Repeats from AFLP Sequence, by using a hybridized membrane. Twenty microsatellite loci were evaluated using 78 individuals from five mangosteen populations, and we successfully amplified four closely related Garcinia species, including G. malaccensis, G. hombroniana, G. celebica, and G. porrecta. Eight loci were monomorphic and the others were polymorphic. Sixty-nine alleles were found, with 3.491 per locus on an average. Genetic diversity (H?) was calculated with an average across loci within population (H?_S) as 0.39, an average loci across many populations (H?_T) as 0.444, and genetic differentiation (F?_ST) as 0.147. Furthermore, based on morphological characters, mangosteen individuals from four populations including Leuwiliang, Wanayasa, Puspahiang, and Kaligesing also had morphologically distinct fruit weight, rind weight, and rind thickness among populations. The study also elucidated the dispersal pattern of mangosteen in Java; the source of the genotype mangosteen in Java population was the Wanayasa population. In addition, we found evidence of tetraploidy in mangosteen. These results have potential applications in future breeding, conservation studies, and genetic assessment of mangosteen and their closely related species.     

Genetic structure of the Russian populations of <Emphasis Type="Italic">Pyrenophora tritici-repentis</Emphasis>, determined by using microsatellite markers

The population genetic structure of plant pathogenic fungus Pyrenophora tritici-repentis was examined using microsatellite (SSR) markers. According to the geographical origin of the pathogen populations, they were designated as North Caucasian (S, 33 isolates), northwest (Nw, 39), and Omsk (Om, 43). The populations were analyzed at the nine most polymorphic SSR loci, at which 75 alleles were identified. To characterize the genetic variation within and between populations, the AMOVA algorithm as implemented in the Arlequin v. 3.5 software program was used. The number of alleles per locus ranged from 5 to 12 and their sizes varied within the range from 180 to 400 bp. The mean gene diversity at SSR loci was high for all populations (H = 0.58–0.75). The populations were considerably different in the frequencies of individual alleles of the SSR loci. Most isolates in the populations were represented by unique haplotypes. The within-population variation of the isolates at molecular markers was 86.4%; among the populations, 13.6%. Substantial interpopulation differences were found between the Om and S (F_st = 0.16) and between the Om and Nw (F_st = 0.20) populations, whereas between the S and Nw populations, these differences were small (F_st = 0.05). Thus, it was demonstrated that the population of P. tritici-repentis from Omsk oblast had the independent status of the geographical population; northwest and North Caucasian populations differed in the allelic diversity of SSR loci, and despite the low F_st value (0.05), they also belonged to independent geographical populations.     

<Emphasis Type="Italic">Bruguiera hainesii</Emphasis>, a critically endangered mangrove species,is a hybrid between <Emphasis Type="Italic">B. cylindrica</Emphasis> and <Emphasis Type="Italic">B. gymnorhiza</Emphasis> (Rhizophoraceae)

Bruguiera hainesii (Rhizophoraceae) is one of the two Critically Endangered mangrove species listed in the IUCN Red List of Threatened Species. Although the species is vulnerable to extinction, its genetic diversity and the evolutionary relationships with other Bruguiera species are not well understood. Also, intermediate morphological characters imply that the species might be of hybrid origin. To clarify the genetic relationship between B. hainesii and other Bruguiera species, we conducted molecular analyses including all six Bruguiera species using DNA sequences of two nuclear genes (CesA and UNK) and three chloroplast regions (intergenic spacer regions of trnL-trnF, trnS-trnG and atpB-rbcL). For nuclear DNA markers, all nine B. hainesii samples from five populations were heterozygous at both loci, with one allele was shared with B. cylindrica, and the other with B. gymnorhiza. For chloroplast DNA markers, the two haplotypes found in B. hainesii were shared only by B. cylindrica. These results suggested that B. hainesii is a hybrid between B. cylindrica as the maternal parent and B. gymnorhiza as the paternal one. Furthermore, chloroplast DNA haplotypes found in B. hainesii suggest that hybridization has occurred independently in regions where the distribution ranges of the parental species meet. As the IUCN Red List of Threatened Species currently excludes hybrids (except for apomictic plant hybrids), the conservation status of B. hainesii should be reconsidered.