Molecular evidence for a natural diploid hybrid between <Emphasis Type="Italic">Miscanthus</Emphasis><Emphasis Type="Italic">sinensis</Emphasis> (Poaceae) and <Emphasis Type="Italic">M. sacchariflorus</Emphasis>

The hybrid origin of Miscanthus purpurascens has previously been proposed, primarily because of its intermediate morphology. In this study, phylogenies based on the DNA sequences from the internal transcribed spacer region of nuclear ribosomal DNA (nrDNA ITS), on the DNA sequences of the trnL intron and trnL-F intergenic spacer of chloroplast DNA, and on amplified fragment length polymorphism (AFLP) fingerprinting confirm that M. purpurascens originated through homoploid hybridization between M. sinensis and M. sacchariflorus. Two different types of ITS sequences were identified from almost all plants of M. purpurascens. One type was found to be closely related to M. sinensis and the other to M. sacchariflorus. Miscanthus purpurascens was found to possess many M. sinensis- and M. sacchariflorus-specific AFLP bands but no band specific to itself. Clustering with the Unweighted Pair Group Method with Arithmetic Mean and principal coordinate analysis based on the AFLP data also demonstrated that M. purpurascens is an approximate intermediate of the two species. In addition, M. purpurascens has the plastid genome of M. sinensis or M. sacchariflorus, suggesting that either species could be its maternal parent. All specimens of M. purpurascens and its coexisting parental species are identified as diploids (2n = 2x = 38). Possible mechanisms of natural hybridization, hybrid status, chloroplast DNA recombination, and evolutionary implications of this hybridization are also discussed.     

Genetic diversity and reproductive traits in triploid and tetraploid populations of <Emphasis Type="Italic">Gladiolus tenuis</Emphasis> (Iridaceae)

Most perennial herbaceous plants are able to reproduce vegetatively as well as sexually. Sometimes, such plants may lose the capacity for sexual reproduction. We have studied the case of sterility in triploid populations (2n = 3x = 45) of Gladiolus tenuis M.Bieb. in a considerable part of its area of distribution. Initially, we recorded the presence of a large clone of G. tenuis to the east of the Volga River, as a result of isozyme analysis. We also used AFLP fingerprinting to genotype 55 samples from 10 populations of G. tenuis and one population of the closely related G. imbricutus L. This analysis revealed an extremely low genetic diversity in sterile triploid populations of G. tenuis and a rather high genetic diversity in fertile tetraploid populations (2n = 4x = 60) over most of the area of this species. Genetic distances between fertile and sterile populations of G. tenuis were similar to those between different species of gladioli. It appears that a single sterile genotype has spread vegetatively over 800 km, propagating by daughter corms. The study of the reproductive features of G. tenuis suggests that the cause of sterility may be self-incompatibility between individuals of the clone.     

High genetic diversity in an endangered medicinal plant, <Emphasis Type="Italic">Saussurea involucrata</Emphasis> (<Emphasis Type="Italic">Saussurea</Emphasis>, Asteraceae), in western Tianshan Mountains,China

Saussurea involucrata (Asteraceae) is a medicinal and second-degree national priority endangered plant that is mainly distributed in the high latitude region of the western Tianshan Mountains. The population is fragmented and isolated, and extensive human impact merits a suitable and specific conservation strategy, which can be compiled based on the genetic diversity, population structure, and demographic history. Phylogeographic studies were conducted on a total of five natural populations and 150 individuals were sampled. Data from three cpDNA intergenic spacer regions (trnL-F, matK, and ndhF-rpl32) and nrDNA ITS sequences showed that twelve haplotypes in cpDNA and five haplotypes in nrDNA indicated high genetic diversity among populations sampled (H _T?=?0.820 and 0.756) and within populations sampled (H _S?=?0.792 and 0.721). Additionally, the high genetic diversity did not mirror genetic structure in either cpDNA (F _ST?=?0.03153, G _ST?>?N _ST, p?<?0.05) or nrDNA (F _ST?=?0.03666, meaningless G _ST and N _ST). Two groups (north and south) were determined for a SAMOVA analysis. Based on this analysis, the demographic history was conducted with a Bayesian Skyline Plot and Isolation with Migration analysis, which showed sustainable and stable extension without a marked bottleneck. Divergence time was indicated at c. 6.25 Mya (90%HPD: 15.30–0.22 Mya) in the Miocene, which is consistent with the formation of the Kelasu section of Tianshan. The southern populations in the Bayanbulak and Gonglu regions require additional attention and transplanting would be an effective way to restore rare cpDNA haplotypes, increase effective population size, and migration rate. Our results suggested that in situ conservation of S. involucrata in western Tianshan should be the main strategy for protection and recovery of the species.     

Molecular phylogenetic analysis of key <Emphasis Type="Italic">Jatropha</Emphasis> species inferred from nrDNA ITS and chloroplast (<Emphasis Type="Italic">trn</Emphasis>L-F and <Emphasis Type="Italic">rbc</Emphasis>L) sequences

The genus Jatropha (Euphorbiaceae) contains species that are of significant economic and ornamental value. However, Jatropha breeding material is rather limited due to incomplete information regarding phylogenetic relationships among germplasm resources. Phylogenetic analyses were performed based on the internal transcribed spacer of nuclear ribosomal DNA (nrDNA ITS), two chloroplast regions (trnL-F and rbcL), and the combined (ITS+trnL-F+rbcL) dataset among twenty-five specimens representing six key Jatropha species. Phylogenetic relationships of Jatropha were well resolved between subgenus Curcas and subgenus Jatropha, and demonstrated the intermediate position of section Polymorphae among sections of both subgenera. Jatropha curcas and J. integerrima demonstrated a close phylogenetic relationship. The molecular data agreed with the morphological classification that recognized J. multifida and J. podagrica in sec. Peltatae. The distinct intraspecific divergence that occurred in J. curcas could be attributed to restricted gene flow caused by geographical isolation and different ecological conditions. Phylograms produced with trnL-F and rbcL sequence data suggested slow rates of sequence divergence among Jatropha spp., while the ITS gene tree had good resolution suggesting high genetic variation of ITS among Jatropha species.     

<Emphasis Type="Italic">Allophoma hayatii</Emphasis> sp. nov., an undescribed pathogenic fungus causing dieback of <Emphasis Type="Italic">Lantana camara</Emphasis> in Iran

Dieback and canker diseases are a major problem in ornamental shrubs and trees of Ahvaz, southwestern Iran. Symptomatic stems and branches were collected from two urban parks in the downtown regions of Ahvaz to identify the dieback-causing agents of Lantana camara. Accordingly, two isolates of a new species, Allophoma hayatii sp. nov., were obtained, which are described and illustrated. This species is identified based on morphological characters and analyses of nucleotide sequences of four regions, including internal transcribed spacer 1 and 2 and 5.8S nrDNA (ITS), partial large subunit 28S nrDNA (LSU-D1/D2), a partial sequence of the β-tubulin (tub2) and part of the RNA polymerase II (rpb2). The isolates of A. hayatii generated a well-supported clade in the trees constructed from multi-locus phylogenetic analysis, distinct from other previously known species of Allophoma. Pathogenicity of both isolates was verified by the inoculation of stem fragments of L. camara. These findings confirm A. hayatii as the causal agents of dieback and canker of L. camara in Ahvaz.     

Differentiation of the endemic Greek genus <Emphasis Type="Italic">Hymenonema</Emphasis> and its relatives of subtribe Scolyminae (Compositae,Cichorieae) based on a multilocus species tree reconstruction

Hymenonema (Compositae, tribe Cichorieae) together with the genera Catananche, Gundelia, and Scolymus forms the subtribe Scolyminae. It is endemic to Greece and consists of two species, Hymenonema laconicum and Hymenonema graecum, which occur in the south Peloponnisos and central Aegean area, respectively. The present contribution aims at a phylogenetic reconstruction of evolutionary relationships among the 12 species of the subtribe, focusing on the temporal and spatial framework for its evolution. The phylogenetic relationships among the members of Scolyminae were inferred from molecular data based on the multi-copy region of the nrDNA internal transcribed spacers ITS1 and ITS2, two intergenic spacers of the cpDNA (trnL-trnF, rpl32-trnL), and one single-copy nuclear region (D10). The gene trees were reconstructed using Bayesian phylogenetic methods. All gene trees support the monophyly of Hymenonema and the sister-group relationship with the genus Scolymus. The further sister-group relationship of this group (Hymenonema–Scolymus) with Catananche is also supported by nrDNA and cpDNA analyses. Finally, a species tree (inferred in a Bayesian coalescent framework) was reconstructed and dates the divergence time between the two Hymenonema species to the Pleistocene (around 1.3 Ma ago). Maximum likelihood-based biogeographical reconstructions suggest a Miocene (pre-Messinian) differentiation of the subtribe on the northern Tethyan platform, followed by Miocene/Pliocene dispersal events to the western Mediterranean and North-African platforms and final, small-scale vicariance events within the genera in the Pleistocene.     

Integrated gene mapping and synteny studies give insights into the evolution of a sex proto-chromosome in <Emphasis Type="Italic">Solea senegalensis</Emphasis>

The evolution of genes related to sex and reproduction in fish shows high plasticity and, to date, the sex determination system has only been identified in a few species. Solea senegalensis has 42 chromosomes and an XX/XY chromosome system for sex determination, while related species show the ZZ/ZW system. Next-generation sequencing (NGS), multi-color fluorescence in situ hybridization (mFISH) techniques, and bioinformatics analysis have been carried out, with the objective of revealing new information about sex determination and reproduction in S. senegalensis. To that end, several bacterial artificial chromosome (BAC) clones that contain candidate genes involved in such processes (dmrt1, dmrt2, dmrt3, dmrt4, sox3, sox6, sox8, sox9, lh, cyp19a1a, amh, vasa, aqp3, and nanos3) were analyzed and compared with the same region in other related species. Synteny studies showed that the co-localization of dmrt1-dmrt2-drmt3 in the largest metacentric chromosome of S. senegalensis is coincident with that found in the Z chromosome of Cynoglossus semilaevis, which would potentially make this a sex proto-chromosome. Phylogenetic studies show the close proximity of S. senegalensis to Oryzias latipes, a species with an XX/XY system and a sex master gene. Comparative mapping provides evidence of the preferential association of these candidate genes in particular chromosome pairs. By using the NGS and mFISH techniques, it has been possible to obtain an integrated genetic map, which shows that 15 out of 21 chromosome pairs of S. senegalensis have at least one BAC clone. This result is important for distinguishing those chromosome pairs of S. senegalensis that are similar in shape and size. The mFISH analysis shows the following co-localizations in the same chromosomes: dmrt1-dmrt2-dmrt3, dmrt4-sox9-thrb, aqp3-sox8, cyp19a1a-fshb, igsf9b-sox3, and lysg-sox6.     

Morphological and molecular investigations of <Emphasis Type="Italic">Gagea</Emphasis> (Liliaceae) in southeastern Kazakhstan with special reference to putative altitudinal hybrid zones

The most important center of speciation in the genus Gagea is thought to be in Central Asia. Here, we focus on species diversity in southeastern Kazakhstan (around Almaty, Ili-Alatau range of the Western Tian-Shan mountains). We studied an elevational transect, reaching from lowland steppes to the alpine zone (500–2750 m a. s. l.), and carried out detailed morphological and molecular investigations for populations of Gagea spp. Nine species were identified in different altitudinal zones; one of these (Gagea almaatensis) is described as new to science. We could detect two altitudinal contact zones between closely related species: G. filiformis and G. granulosa (sect. Minimae), and G. almaatensis and G. kuraminica (sect. Gagea). Morphological and molecular investigations (ITS data and cpDNA networks) indicate ongoing hybridization of co-occurring G. filiformis into G. granulosa and putative bidirectional hybridization events between G. almaatensis and G. kuraminica.     

<Emphasis Type="Italic">Viola woosanensis</Emphasis>, a recurrent spontaneous hybrid between <Emphasis Type="Italic">V. ulleungdoensis</Emphasis> and <Emphasis Type="Italic">V. chaerophylloides</Emphasis> (Violaceae) endemic to Ulleung Island,Korea

Ulleung Island is an oceanic volcanic island in Korea, which has never been connected to the adjacent continent. Previous studies highlighted Ulleung Island as an excellent system to study the pattern and process of early stages of flowering plant evolutions on oceanic island. The predominant mode of speciation in flowering plants on Ulleung Island appears to be anagenesis. However, the potentially important role of hybrid speciation among incompletely reproductively isolated lineages cannot be ruled out. Viola woosanensis (Violaceae) is of purportedly hybrid origin between V. ulleungdoensis (i.e., formerly recognized as V. selkirkii in Ulleung Island) and V. chaerophylloides, based on morphology. To examine the origin of V. woosanensis, we sampled a total of 80 accessions, including V. woosanensis and its putative parental species and sequenced nrDNA ITS, and four highly variable chloroplast noncoding regions (trnL-trnF, rpl16 intron, atpF-atpH, and psbA-trnH). Representative species of Viola from Korea were also included in the phylogenetic analyses (maximum parsimony, maximum likelihood, and Bayesian inference). Additive polymorphic sites in the nrDNA ITS regions were confirmed by cloning amplicons from representative species. The molecular data strongly supported the hybrid origin of V. woosanensis, and the maternal and paternal parent were determined to be V. ulleungdoensis and V. chaerophylloides, respectively. The presence of two parental ribotypes in V. woosanensis (with the exception in one population) was confirmed by cloning, suggesting V. woosanensis is primarily the F₁ generation. No trace of backcrossing and introgression with its parents was detected due to low fertility of hybrid species. We found a multiple and unidirectional hybrid origin of V. woosanensis. Additional studies are required to determine which factors contribute to asymmetric gene flow of Viola species in Ulleung Island.     

Conservation genetics and geographic patterns of genetic variation of endangered shrub <Emphasis Type="Italic">Ammopiptanthus</Emphasis> (Fabaceae) in northwestern China

Phylogeographic patterns of Ammopiptanthus in northwestern China were examined with internal transcribed spacer (ITS) and three chloroplast intergenic spacers (trnH–psbA, trnL–trnF, and trnS–trnG). Two ITS genotypes (a–b) and 8 chloroplast haplotypes (A–H) were detected. Both ITS genotypes and chloroplast lineages were split in two geographic regions: western Xinjiang and the Alxa Desert. This lineage split was also supported by AMOVA analysis and the Mantel test. AMOVA showed that 89.81 % of variance in Ammopiptanthus occurred between the two geographic regions, and correlation between genetic distances and geographical distances was significant (r = 0.757, p < 0.0001). All populations in western Xinjiang shared haplotype A with high frequency, and range expansion was strongly supported by negative Fu’s F_S value, and mismatch distribution analysis, whereas populations in the Alxa Desert had higher genetic diversity and structure. We speculate that the cold and dry climate during the early Quaternary fragmented habitats of the species, limiting gene flow between regions, and interglacial periods most likely led to the range expansion in western Xinjiang. The low genetic diversity of Ammopiptanthus indicate a significant extinction risk, and protective measures should be taken immediately.     

Fingerprinting by amplified fragment length polymorphism (AFLP) and barcoding by three plastidic markers in the genus <Emphasis Type="Italic">Wolffiella</Emphasis> Hegelm

Amplified fragment length polymorphism (AFLP) fingerprinting and three different plastidic DNA regions (rpl16, rps16, atpF-atpH) were used to investigate species identity in the genus Wolffiella. For this purpose, clones (67 in total) belonging to all ten species were selected. Almost all the species were represented by more than one clone. The fingerprinting technique, AFLP, clearly distinguished the species, W. caudata, W. gladiata, W. neotropica, W. rotunda, and W. welwitschii. Apart from confirming the molecular identity of these five species, the plastidic markers could delineate two additional species, W. hyalina and W. denticulata, although the conclusion concerning the latter is restricted by the availability of only one clone. The efficiency of the plastid-derived markers in identifying the number of species-specific clusters followed the sequence rps16 > rpl16 > atpF-atpH. The species W. lingulata, W. oblonga, and W. repanda could not be identified by any of the molecular methods presented here, but could be strictly defined on a morphological basis. In several clones, high amounts of genetic admixtures between different species were detected. Further, simulation studies demonstrated that these clones are genetic hybrids. This might be one of the obstacles in molecular identification of species in the genus Wolffiella.     

Species diversity of <Emphasis Type="Italic">Pseudocercospora</Emphasis> from Far East Asia

This study reflects on the monophyly of, and species diversity within, the genus Pseudocercospora in Far East Asia. Morphological characteristics and phylogenetic analyses of Pseudocercospora species were based on type specimens and ex-type cultures, which were collected from Japan and Taiwan. A phylogenetic tree was generated from multi-locus DNA sequence data of the internal transcribed spacer regions of the nrDNA cistron (ITS), partial actin (actA), and partial translation elongation factor 1-alpha (tef1), as well as the partial DNA-directed RNA polymerase II second largest subunit (rpb2). Based on these results, Pseudocercospora amelanchieris on Amelanchier and Ps. iwakiensis on Ilex were newly described from Japan, and a further 22 types (incl. two neo-, five lecto-, and 15 epitypes), were designated. The genus Pseudocercospora as presently circumscribed was found to be monophyletic, while the secondary barcodes, actA, tef1, and rpb2 were shown to be well suited to delimitate species within the genus.     

Drivers of natural spread of invasive <Emphasis Type="Italic">Impatiens parviflora</Emphasis> differ between life-cycle stages

Impatiens parviflora is one of the most widely spread invasive species in central Europe, yet the factors affecting its spread are still subject to discussion. The aim of this study was to determine which factors affect the spread of I. parviflora. This was achieved by observing the natural spread of the species on 15 permanent transects, in six different habitats, found within the Czech Republic from 2012 to 2016. The transects were divided into 321 plots, whilst data on environmental conditions, as well as of the spread of I. parviflora and its performance, were collected in each plot. The results showed that individual stages of the I. parviflora life-cycle were affected by individual environmental conditions to different extents. The most important factor preventing seedling emergence and establishment was a high cover of herb layer. It did not, however, affect survival of older plants. Thus, I. parviflora can grow in sites with dense cover of herb layer in case the cover formed after I. parviflora seedlings established. Juvenile mortality was the highest in sites with low nutrient levels and low soil moisture. Canopy openness had a negative effect on I. parviflora performance. Impatiens parviflora performed better in neutral soils, in comparison to acidic soils. Oak-hornbeam forests were the most suitable habitat for I. parviflora, followed by acidophilous oak and mixed coniferous forests. However, I. parviflora was able to penetrate into even species-rich habitats, such as thermophilous oak forests, as well as steppe grasslands on rocks. This makes it a potential threat to biodiversity. Only heathlands found on former pastures proved to be unsuitable for I. parviflora, as these remained uninvaded until the end of the study.     

Characterisation of seven <Emphasis Type="Italic">Inocybe</Emphasis> ectomycorrhizal morphotypes from a semiarid woody steppe

Ectomycorrhizas (ECM) of Inocybe species (Inocybaceae, Basidiomycota) formed by three host plant species (Populus alba, Salix rosmarinifolia and Pinus nigra) in a semiarid woody steppe of Hungary were studied. To identify the fungal partners, we performed phylogenetic analyses of nucleotide sequences for the internal transcribed spacer region of nuclear DNA (nrDNA ITS) together with sequences gained from public databases. Seven Inocybe ectomycorrhiza morphotypes were morpho-anatomically characterised. Five morphotypes were identified (I. phaeoleuca, I. psammophila, I. semifulva, I. splendens and I. subporospora), whereas two morphotypes represented unidentified Inocybe species. Differences were discernible among the morphotypes, and they showed general anatomical characteristics of Inocybe ECM, such as the slightly organised plectenchymatic mantle (types A, B and E and the gelatinous C). The ECM of I. subporospora and I. phaeoleuca were detected from the introduced Pinus nigra. These two fungi are probably native to the area but capable of forming a novel ectomycorrhizal association with the invasive host.     

Phylogenetic relationships in Iranian <Emphasis Type="Italic">Scutellaria</Emphasis> (Lamiaceae) based on nuclear ribosomal ITS and chloroplast <Emphasis Type="Italic">trn</Emphasis>L-F DNA data

The genus Scutellaria (Lamiaceae, Scutellarioidae) includes taxa with a wide range of pharmacological traits, which make this genus economically important. Previous taxonomical studies on Scutellaria were mainly based on morphological and chromosomal characters. The infrageneric classification of the genus has been the subject of a long controversy among authors, and delimitation of taxa remains problematic. In this study, the phylogeny of the genus was evaluated using nrDNA ITS and trnL-F sequences, with a focus on the 42 Iranian taxa belonging to subgenera Apeltanthus and Scutellaria. In both ITS and trnL-F trees, there were two main clades within the genus corresponding to the two subgenera Scutellaria and Apeltanthus. In the subg. Scutellaria, the sect. Anaspis was monophyletic, but the monophyly of sect. Scutellaria was not supported. In this subgenus, S. galericulata was located away from other species in both trees, supported the isolation of this species from the sect. Scutellaria and inclusion in the sect. Galericularia as was previously considered by some authors. In subg. Apeltanthus, sect. Apeltanthus and sect. Lupulinaria were not monophyletic; furthermore, the monophyly of subsect. Lupulinaria was not supported. Based on the results of this study, it can be suggested that the sectional classification of both subgenera Apeltanthus and Scutellaria should be revised.     

Polyphasic characterisation of three novel species of <Emphasis Type="Italic">Paraboeremia</Emphasis>

Paraboeremia was recently introduced for a distinct lineage in the family Didymellaceae. Currently, three species are included, i.e. P. adianticola, P. putaminum and P. selaginellae, all of which are plant pathogens. Paraboeremia is morphologically similar to Phoma but phylogenetically distinct. In this paper, three new species, i.e. Paraboeremia camelliae isolated from Camellia sp., P. litseae from Litsea sp., and P. oligotrophica from cave limestone, are described, illustrated and compared with closely related taxa. Phylogenetic analysis based on the multi-locus sequences of the internal transcribed spacer regions 1 and 2 and 5.8S nuclear ribosomal RNA gene (ITS), partial large subunit 28S nrDNA region (LSU), partial β-tubulin (TUB2) gene and RNA polymerase II (RPB2) gene regions confirmed the distinction of these species in Paraboeremia. These three new species were discovered from habitats and hosts that are previously unknown from this genus.     

Spatial patterns and intraspecific diversity of the glacial relict legume species <Emphasis Type="Italic">Vavilovia formosa</Emphasis> (Stev.) Fed. in Eurasia

Vavilovia formosa is one of five genera in tribe Fabeae, (Fabaceae, Leguminosae) with close phylogenetic relationships to Pisum. It grows in subalpine and alpine levels in Armenia, Azerbaijan, Georgia, Iran, Iraq, Lebanon, Russia and Turkey and is recognized as an endangered and protected plant. This study was conducted to reveal its intraspecific variability, as well as to predict the past, extant and future species distribution range. We analysed 51 accessions with common phylogenetic markers (trnF-trnL, trnS-trnG, matK, rbcL, psbA-trnH and ITS). These represent in total up to 2551 bp of chloroplast and 664 bp of nuclear sequences per sample. Two populations from Turkey and Armenia were analysed for genetic diversity by AFLP. Leaf morphometry was conducted on 1457 leaflets from 43 specimens. Extracted bioclimatic parameters were used for niche-modelling approach. Analysis of cpDNA revealed two haplotypes, 12 samples from Armenia, Daghestan, Nakhichevan and Iran belonged to H1 group, while 39 samples of all Turkish and part of Armenian were in H2 group. The mean intrapopulation diversity based on AFLP was low (H _E = 0.088) indicating limited outcrossing rate. A significantly positive correlation between geographical latitude and leaf area (\(\rho\) = 0.527, p < 0.05) was found. Niche modelling has shown temporal variation of predicted occurrence across the projected time periods. Vavilovia formosa has suffered a range reduction following climate warming after last glacial maximum, which classify this species as cold-adapted among the Fabeae species as well as a glacial relict.     

Investigating the genetic diversity and differentiation patterns in the <Emphasis Type="Italic">Penstemon scariosus</Emphasis> species complex under different sample sizes using AFLPs and SSRs

Habitat fragmentation due to anthropogenic activities is the major cause of biodiversity loss. Endemic and narrowly distributed species are the most susceptible to habitat degradation. Penstemon scariosus is one of many species whose natural habitat is vulnerable to industrialization. All varieties of P. scariosus (P. scariosus var. albifluvis, P. scariosus var. cyanomontanus, P. scariosus var. garrettii, P. scariosus var. scariosus) have small distribution ranges, but only P. scariosus var. albifluvis is being considered for listing under the Endangered Species Act. We used eight microsatellites or simple sequence repeats (SSRs) loci and two amplified fragment length polymorphism (AFLP) primer combinations to investigate the population genetic structure and diversity of P. scariosus varieties. Moreover, we compared the utility of the two marker systems in conservation genetics and estimated an appropriate sample size in population genetic studies. Genetic differentiation among populations based on F_st ranged from low to moderate (F_st?=?0.056–0.157) and from moderate to high when estimated with D_es (D_es?=?0.15–0.32). Also, AMOVA analysis shows that most of the genetic variation is within populations. Inbreeding coefficients (F_is) were high in all varieties (0.20–0.56). The Bayesian analysis, STRUCTURE, identified three clusters from SSR data and four clusters from AFLPs. Clusters were not consistent between marker systems and did not represent the current taxonomy. MEMGENE revealed that a high proportion of the genetic variation is due to geographic distance (R²?=?0.38, P?=?0.001). Comparing the genetic measurements from AFLPs and SSRs, we found that AFLP results were more accurate than SSR results across sample size when populations were larger than 25 individuals. As sample size decreases, the estimates become less stable in both AFLP and SSR datasets. Finally, this study provides insight into the population genetic structure of these varieties, which could be used in conservation efforts.