Assessing genetic diversity of Citrus by DArT_seq™ genotyping

In citrus despite the diversity among cultivated genotypes related to morphological, physiological, and agronomic traits, low polymorphism is detected at the molecular marker level, and the number of DNA-markers technologies has been insufficient to identify genetic diversity among all citrus varieties. However, DArT_seq? markers can overcome this limitation and identify genetic diversity in several species of citrus. This work developed and applied DArT_seq? platform to study the relationships among species, cultivars, and hybrids of citrus. DArT_seq? yielded a total of 37,260 polymorphic presence/absence markers that generated a dendrogram of similarity among the studied genotypes. The results confirmed the relationships of sweet orange, mandarin, and citron species. DArT_seq? markers showed extensive genetic variation among citrus species and can be applied to Citrus genetic diversity studies.     

Genetic Diversity and Population Structure Analysis of European Hexaploid Bread Wheat (Triticum aestivum L.) Varieties

Progress in plant breeding is facilitated by accurate information about genetic structure and diversity. Here, Diversity Array Technology (DArT) was used to characterize a population of 94 bread wheat (Triticum aestivum L.) varieties of mainly European origin. In total, 1,849 of 7,000 tested markers were polymorphic and could be used for population structure analysis. Two major subgroups of wheat varieties, GrI and GrII, were identified using the program STRUCTURE, and confirmed by principal component analysis (PCA). These subgroups were largely separated according to origin; GrI comprised varieties from Southern and Eastern Europe, whereas GrII contained mostly modern varieties from Western and Northern Europe. A large proportion of the markers contributing most to the genetic separation of the subgroups were located on chromosome 2D near the Reduced height 8 (Rht8) locus, and PCR-based genotyping suggested that breeding for the Rht8 allele had a major impact on subgroup separation. Consistently, analysis of linkage disequilibrium (LD) suggested that different selective pressures had acted on chromosome 2D in the two subgroups. Our data provides an overview of the allele composition of bread wheat varieties anchored to DArT markers, which will facilitate targeted combination of alleles following DArT-based QTL studies. In addition, the genetic diversity and distance data combined with specific Rht8 genotypes can now be used by breeders to guide selection of crossing parents.     

Use of SRAP markers to assess genetic diversity and population structure of wild, cultivated, and ornamental pomegranates (Punica granatum L.) in different regions of Iran

Sequence-related amplified polymorphism (SRAP) was used to assess the genetic diversity of 63 cultivated, wild, and ornamental pomegranate genotypes from five different geographical regions of Iran. A total of 250 fragments were amplified using 13 primer combinations; among these, 133 bands (53?%) were polymorphic. The average PIC value was 0.28 over all PCs. The genetic distance among genotypes ranged from 0.10 to 0.37 with an average of 0.24. Cluster analysis using the neighbor-joining (NJ) method suggested there are close relationships between ornamental and some wild genotypes. Although AMOVA results revealed significant differences in the genetic diversity among the regions (P?=?0.0048), the genetic variation was mainly caused by variation of intra regions. The results indicated low genetic differentiation (Fst?=?0.025) and high gene flow (Nm?=?2.28) among regions. These results confirmed that SRAP markers could be powerful tools and an effective marker system for determining the genetic diversity and population genetic structure of the pomegranate.     

Application of EST-derived microsatellite markers for analysis of genetic variation in tall fescue and its comparison with morphological markers

Analysis of genetic diversity in germplasm collections is an important component of crop improvement programs. This study was conducted to analyze genetic variation and to classify tall fescue genotypes based on phenotypic evaluation and EST-SSR molecular markers. Twenty-five genotypes were assessed based on phenotypic and 42 EST-SSR molecular markers according to a completely randomized block design with three replications during eight years (2007–2014). Results indicated that the effect of year, genotype and their interaction were significant for all of the measured traits. Both morphological and molecular assessments showed considerable genetic variation among genotypes. The estimates of broad-sense heritability (h²_b) were moderate to high (h²_b = 42.1–78.4) for the traits studied. Based on EST-SRR analysis, a total number of 229 alleles were detected with an average of 4.58 alleles per marker. Average PIC value was 0.49 with a range of 0.014 for NFA140 to 0.95 for NFA047. Phenotypic evaluations and EST-SSR molecular marker classified genotypes into 3 and 7 clusters, respectively which mainly supported geographical origins. The general correspondence was observed between morphological and molecular classification. Therefore, combining the molecular markers with morphological responses could be more beneficial to describe genetic variation and distinguish superior genotypes for future breeding programs.     

Genetic variation, population structure, and linkage disequilibrium in European elite germplasm of perennial ryegrass

Perennial ryegrass (Lolium perenne L.) is a highly valued temperate climate grass species grown as forage crop and for amenity uses. Due to its outbreeding nature and recent domestication, a high degree of genetic diversity is expected among cultivars. The aim of this study was to assess the extent of linkage disequilibrium (LD) within European elite germplasm and to evaluate the appropriate methodology for genetic association mapping in perennial ryegrass. A high level of genetic diversity was observed in a set of 380 perennial ryegrass elite genotypes when genotyped with 40 SSRs and 2 STS markers. A Bayesian structure analysis identified two subpopulations, which were confirmed by principal coordinate analysis (PCoA). One subpopulation consisted mainly of genotypes originating from the UK, while germplasm mostly from Continental Europe was grouped into the second subpopulation. LD (r²) decay was rapid and occurred within 0.4 cM across European varieties, when population structure was taken into consideration. However, an extended LD of up to 6.6 cM was detected within the variety Aberdart. High genetic diversity and rapid LD decay provide means for high resolution association mapping in elite materials of perennial ryegrass. However, different strategies need to be applied depending on the material used. Genome-wide association study (GWAS) with several hundred markers can be applied within synthetic varieties to identify large (up to 10 cM) genomic regions affecting trait variation. A combination of available and novel DNA markers is needed to achieve resolution required for GWAS in elite breeding materials. An even higher marker density of several million SNPs might be needed for GWAS in diverse ecotype collections, potentially resulting in quantitative trait polymorphism (QTP) identification.     

SSR analysis of 38 genotypes of soybean (Glycine Max (L.) Merr.) genetic diversity in India

Sixteen polymorphic Simple sequence repeat (SSR) markers were used to determine the genetic diversity and varietal identification among 38 soybean (Glycine max (L.) Merr.) genotypes which are at present under seed multiplication chain in India. A total of 51 alleles with an average of 2.22 alleles per locus were detected. The polymorphic information content (PIC) among genotypes varied from 0.049 (Sat_243 and Satt337) to 0.526 (Satt431) with an average of 0.199. The pair wise genetic similarity between soybean varieties varied from 0.56 to 0.97 with an average of 0.761. These 16 SSR markers successfully distinguished 12 of the 38 soybean genotypes. These results suggest that used SSR markers are efficient for measuring genetic diversity and relatedness as well as identifying varieties of soybeans. Diverse genetic materials may be used for genetic improvements of soybean genotypes.     

Promoting Utilization of Saccharum spp. Genetic Resources through Genetic Diversity Analysis and Core Collection Construction

Sugarcane (Saccharum spp.) and other members of Saccharum spp. are attractive biofuel feedstocks. One of the two World Collections of Sugarcane and Related Grasses (WCSRG) is in Miami, FL. This WCSRG has 1002 accessions, presumably with valuable alleles for biomass, other important agronomic traits, and stress resistance. However, the WCSRG has not been fully exploited by breeders due to its lack of characterization and unmanageable population. In order to optimize the use of this genetic resource, we aim to 1) genotypically evaluate all the 1002 accessions to understand its genetic diversity and population structure and 2) form a core collection, which captures most of the genetic diversity in the WCSRG. We screened 36 microsatellite markers on 1002 genotypes and recorded 209 alleles. Genetic diversity of the WCSRG ranged from 0 to 0.5 with an average of 0.304. The population structure analysis and principal coordinate analysis revealed three clusters with all S. spontaneum in one cluster, S. officinarum and S. hybrids in the second cluster and mostly non-Saccharum spp. in the third cluster. A core collection of 300 accessions was identified which captured the maximum genetic diversity of the entire WCSRG which can be further exploited for sugarcane and energy cane breeding. Sugarcane and energy cane breeders can effectively utilize this core collection for cultivar improvement. Further, the core collection can provide resources for forming an association panel to evaluate the traits of agronomic and commercial importance.     

Identification of Genetic Loci Associated with Quality Traits in Almond via Association Mapping

To design an appropriate association study, we need to understand population structure and the structure of linkage disequilibrium within and among populations as well as in different regions of the genome in an organism. In this study, we have used a total of 98 almond accessions, from five continents located and maintained at the Centro de Investigación y Tecnología Agroalimentaria de Aragón (CITA; Spain), and 40 microsatellite markers. Population structure analysis performed in ‘Structure’ grouped the accessions into two principal groups; the Mediterranean (Western-Europe) and the non-Mediterranean, with K = 3, being the best fit for our data. There was a strong subpopulation structure with linkage disequilibrium decaying with increasing genetic distance resulting in lower levels of linkage disequilibrium between more distant markers. A significant impact of population structure on linkage disequilibrium in the almond cultivar groups was observed. The mean r² value for all intra-chromosomal loci pairs was 0.040, whereas, the r² for the inter-chromosomal loci pairs was 0.036. For analysis of association between the markers and phenotypic traits, five models comprising both general linear models and mixed linear models were selected to test the marker trait associations. The mixed linear model (MLM) approach using co-ancestry values from population structure and kinship estimates (K model) as covariates identified a maximum of 16 significant associations for chemical traits and 12 for physical traits. This study reports for the first time the use of association mapping for determining marker-locus trait associations in a world-wide almond germplasm collection. It is likely that association mapping will have the most immediate and largest impact on the tier of crops such as almond with the greatest economic value.     

European Phaseolus coccineus L. landraces: Population Structure and Adaptation,as Revealed by cpSSRs and Phenotypic Analyses

Relatively few studies have extensively analysed the genetic diversity of the runner bean through molecular markers. Here, we used six chloroplast microsatellites (cpSSRs) to investigate the cytoplasmic diversity of 331 European domesticated accessions of the scarlet runner bean (Phaseolus coccineus L.), including the botanical varieties albiflorus, bicolor and coccineus, and a sample of 49 domesticated and wild accessions from Mesoamerica. We further explored the pattern of diversity of the European landraces using 12 phenotypic traits on 262 individuals. For 158 European accessions, we studied the relationships between cpSSR polymorphisms and phenotypic traits. Additionally, to gain insights into the role of gene flow and migration, for a subset of 115 accessions, we compared and contrasted the results obtained by cpSSRs and phenotypic traits with those obtained in a previous study with 12 nuclear microsatellites (nuSSRs). Our results suggest that both demographic and selective factors have roles in the shaping of the population genetic structure of the European runner bean. In particular, we infer the existence of a moderate-to-strong cytoplasmic bottleneck that followed the expansion of the crop into Europe, and we deduce multiple domestication events for this species. We also observe an adaptive population differentiation in the phenology across a latitudinal gradient, which suggests that selection led to the diversification of the runner bean in Europe. The botanical varieties albiflorus, bicolor and coccineus, which are based solely on flower colour, cannot be distinguished based on these cpSSRs and nuSSRs, nor according to the 12 quantitative traits.     

Lack of structure in the gene pool of the highly polyploid ornamental chrysanthemum

Selection of clonally propagated chrysanthemums is mostly performed on F1 hybrids using phenotypic characteristics without the use of molecular information. We applied 448 amplified fragment length polymorphism markers to a set of 81 accessions, mainly from the European gene pool, covering the different horticultural types (cut, pot and garden varieties) and originating from the most important European chrysanthemum breeders. The average pairwise genetic similarity of 0.69 was moderate to rather high. Neighbour-joining clustering resulted in no grouping of the accessions, either by their common origin or their horticultural type, or by similarities in important phenotypic characteristics. The structure of the dendrogram could not be supported by bootstrap analysis. Furthermore, network analysis using SplitsTree, principal coordinate analysis via DARwin or analysis of the population with structure did not differentiate reliable and invariable clusters. Therefore, we tested the marker saturation by plotting the mean coefficient of variation for every pairwise similarity of the bootstrap analysis against the different numbers of markers. We showed that the number of markers is sufficient for a precise estimate of genetic similarity and that the lack of bootstrap support is not due to a low genetic diversity or a lack of marker information, but most likely resulted from the breeding history of the cultivars, involving repeated backcrosses and the exchange of genotypes between breeders.     

Patterns of subspecies genetic diversity among oilseed <Emphasis Type="Italic">Brassica rapa</Emphasis> as revealed by agro-morphological traits and SSR markers

Brassica rapa (2n = 20, AA genome) is an important oil yielding species of the family Brassicaceae and characterized by wide range of genetic and morphological subtypes suitable for cultivation under diverse agro-climatic regions of India. In this study, genetic diversity among three subspecies of B. rapa including yellow sarson, toria and outlier brown sarson was estimated using various agro-morphological traits and simple sequence repeat (SSR) markers. Maximum variability was recorded for siliqua angle (Coefficient of variation = 30.9%), followed by seeds/siliqua (CV = 18.8%), leaf length (CV = 10%) and plant height (CV = 16.8%). Principal component analysis explained more than 50% of the total observed morphological variability for first two components. Of the 107 SSR markers tested, 80 generated reproducible, clear and distinct amplicons of which, 65 (81.25%) were found polymorphic. The number of alleles at each locus ranged from 2 to 7, with an average of 3.03 alleles per marker. A total of 197 alleles were detected at 65 SSR loci with average PIC value of 0.457 and a mean resolving power of 3.04. Neighbor-Joining cluster analysis based on morphological traits and SSR markers separately classified all the 28 genotypes into five major groups. The population structure analysis resulted into three sub-populations with certain extent of admixture among the earlier established taxonomic sub-groups. Twenty-three unique alleles were detected in thirteen B. rapa varieties. The clustering analysis and principal coordinate analysis outlined the genetic relationships among different varieties belonging to the three subspecies of B. rapa. Genetically diverse genotypes as illustrated by score plots and from the clustering patterns brought out the wide range of diversity present among B. rapa genotypes and the underlying options available for selecting parental genotypes for hybridization and developing high yielding cultivars suitable for Indian conditions.     

Diversification of indigenous gene- pool by using exotic germplasm in lentil (Lens culinaris Medikus subsp. culinaris)

Genetic diversity was studied among 21 accessions of lentil using SSR markers and morphological traits in order to assess the diversification of Indian gene-pool of lentil through introgression of exotic genes and introduction of germplasm. Among these , 16 genotypes either had ‘Precoz’ gene, an Argentine line in their pedigree or genes from introduced lines from ICARDA. Sixty five SSR markers and eight phenotypic traits were used to analyse the level of genetic diversity in these genotypes. Forty three SSR markers (66 %) were polymorphic and generated a total of 177 alleles with an average of 4.1 alleles per SSR marker. Alleles per marker ranged from 2 to 6. The polymorphic information content ranged 0.33 to 0.80 with an average of 0.57, suggesting that SSR markers are highly polymorphic among the studied genotypes. Genetic dissimilarity based a dendrogram grouped these accessions into two main clusters (cluster I and cluster II) and it ranged 33 % to 71 %, suggesting high level of genetic diversity among the genotypes. First three components of PCA based morphological traits explained higher variance (95.6 %) compared to PCA components based on SSR markers (42.7 %) of total genetic variance. Thus, more diversity was observed for morphological traits and genotypes in each cluster and sub-cluster showed a range of variability for seed size, earliness, pods/plant and plant height. Molecular and phenotypic diversity analysis thus suggested that use of germplasm of exotic lines have diversified the genetic base of lentil germplasm in India. This diversified gene-pool will be very useful in the development of improved varieties of lentil in order to address the effect of climate change, to adapt in new cropping systems niches such as mixed cropping, relay cropping, etc. and to meet consumers’ preference.     

Population structure,morphological and genetic diversity within and among melon (Cucumis melo L.) landraces in Iran

ISSR markers were applied to evaluate the genetic diversity and differentiation of 270 individuals of 27 Iranian C. melo landraces of various varietal groups include vars. inodorous, cantalupensis, reticulatus, ameri, dudaim. Genetic diversity among the studied genotypes obtained by GeneAlex analysis (H?=?0.08, I?=?0.12, Na?=?0.77, PPL?=?22.6%). Cluster analysis divided Iranian melon landraces into two main cluster. Non-sweet genotype (dudaim group) was well separated from sweet genotypes (inodorous, ameri, reticulatus, cantalupensis). The most similar genotypes were BANI and TONI (0.95) and the most dissimilar ones were GER and TS (0.58). AMOVA result showed that the percentage of genetic variation among and within Iranian melon is 69% and 31%, respectively. All landraces evaluated based on 10 morphological traits which revealed the diversity of melon varietal groups. Bayesian analysis assigned ten landraces to Pop 1, eight landraces to Pop 2 and nine melon landraces to Pop 3. Bayesian and UPGMA cluster analyses demonstrated the almost related results. Our results indicated that ISSR markers technique alongside polyacrylamide gel analysis could be helpful to discriminate varieties of melon.     

Introgression from modern hybrid varieties into landrace populations of maize (Zea mays ssp. mays L.) in central Italy

Landraces are domesticated local plant varieties that did not experience a deliberate and intensive selection during a formal breeding programme. In Europe, maize landraces are still cultivated, particularly in marginal areas where traditional farming is often practiced. Here, we have studied the evolution of flint maize landraces from central Italy over 50 years of on-farm cultivation, when dent hybrid varieties were introduced and their use was widespread. We have compared an 'old' collection, obtained during the 1950s, before the introduction of hybrids, and a recent collection of maize landraces. For comparison, a sample of maize landraces from north Italy, and of improved germplasm, including hybrids and inbred lines were also used. A total of 296 genotypes were analysed using 21 microsatellites. Our results show that the maize landraces collected in the last 5–10 years have evolved directly from the flint landrace gene pool cultivated in central Italy before the introduction of modern hybrids. The population structure, diversity and linkage disequilibrium analyses indicate a significant amount of introgression from hybrid varieties into the recent landrace populations. No evidence of genetic erosion of the maize landraces was seen, suggesting that in situ conservation of landraces is an efficient strategy for preserving genetic diversity. Finally, the level of introgression detected was very variable among recent landraces, with most of them showing a low level of introgression; this suggests that coexistence between different types of agriculture is possible, with the adoption of correct practices that are aimed at avoiding introgression from undesired genetic sources.     

鹰嘴豆种质资源农艺性状遗传多样性分析

以100份鹰嘴豆种质资源为材料,应用聚类分析和主成分分析方法,对15个主要农艺性状的遗传多样性进行分析。结果表明,参试材料存在广泛的遗传多样性。其中,多样性指数最高的是株高,其次是百粒重;性状变异系数最大的是单株荚数,其次是单株粒重;基于各种质间形态标记的遗传差异,将100份鹰嘴豆种质聚类并划分为4大类群。第Ⅰ类群可作为选育丰产中粒型和株高适中的品种,第Ⅱ类群可作为选育矮秆耐密及特异粒色(型)品种,第Ⅲ类群丰产性较差可作为选育子粒球型、光滑的品种,第Ⅳ类群可作为选育大粒型、适宜机械化收获的品种。9个数量性状的主成分分析结果表明,前4个主成分累计贡献率达73.91%,各主成分性状载荷值反映了主要数量性状的育种选择潜力。综合分析种质资源农艺性状,为鹰嘴豆的有效利用提供一定的科学依据。     

Genetic diversity and recombination within populations of Fusarium pseudograminearum from western Canada.

Genetic diversity within populations of Fusarium pseudograminearum isolated from wheat grains from the Canadian provinces of Alberta and Saskatchewan was investigated. Three restriction enzymes (EcoRI, HaeIII, and PstI) were used to carry out restriction analysis of the nuclear ribosomal DNA (nrDNA) intergenic spacer region (IGS region) and eight primers were used to generate inter-simple sequence-repeat (ISSR) molecular markers. Our study indicated substantially high genetic diversity within these two populations, but low genetic differentiation and frequent gene flow among populations. The IGS data showed no genetic distinction between the two Alberta populations and only minor genetic differentiation between the Saskatchewan and Alberta populations. Analysis of molecular variance indicated that most genetic variability resulted from differences among isolates within populations. Multilocus linkage disequilibrium analysis suggested a panmictic population genetic structure and the occurrence of significant recombination in F. pseudograminearum. Regular gene flow and random mating between isolates from different populations could result in novel genotypes with both improved pathological and biological traits.     

基于主成分分析法的墨兰品种观赏价值评价

以福州地区栽培多年的30个墨兰品种为材料,对其植株高、花色、瓣型等20个观赏性状进行测定和观察,分析墨兰品种观赏性状的变异范围、性状间的相关性以及主成分因子,并利用主成分分析法对30个墨兰品种的观赏价值进行评价。结果表明,不同品种间的16个数量性状均存在不同程度的变异,花部形态存在较大差异,叶片形态和数量变异小;4个质量性状中,花色和叶艺的遗传多样性较高。叶部形态和花部形态的观赏性状间存在显著相关性。20个观赏性状指标可以分为7个主成分因子,分别为株高因子、瓣型因子、花色因子、叶片数因子、花径因子、叶宽因子和叶艺因子,7个主成分因子的方差贡献率累计达全部性状信息的87.35%。通过计算30个品种的重要主成分值并对其进行排序,筛选出4个观赏性状优良的品种,分别为‘白墨’、‘宝山爪’、‘龙梅’和‘金华山’。     

Relationships between four measures of genetic distance and breeding behavior in spring wheat

We estimated the genetic distances among 10 spring wheat genotypes based on pedigree data, morphological traits and AFLP markers, used individually and combined with morphological traits, to find the best predictors of general- and specific-combining abilities among parental genotypes. Ten wheat parents were crossed in a diallel form, disregarding reciprocal hybrids, totaling 45 combinations. The F(1) hybrids, F(2) populations and parents were evaluated in the field in 2007. The experimental plots consisted of 20 plants for F(1) hybrids and 40 plants for parental and F(2) populations. All methods (pedigree data, AFLP markers and morphological traits, used individually and combined) were found to be useful for the assessment of genetic diversity. The significant coefficient correlations ranged from low (0.45) to moderate (0.67) between the distance measures and hybrid performance. There was significant agreement between the distance measures based on AFLP markers vs morphological traits + AFLP markers (r = 0.47) and between pedigree data vs morphological traits + AFLP markers (r = 0.43). The pedigree distance was positively associated with traits 100-kernel weight and grain yield per plant in F(1) (correlations of 0.67 and 0.62, respectively) and F(2) (correlations of 0.62 and 0.59, respectively) generations. These correlation values indicate that the genetic distance, based on pedigree data, could replace diallel crosses for the selection of parents with higher combining ability and with moderate reliability.     

Genome-wide association mapping of yield components and drought tolerance-related traits in cotton

Drought causes serious yield losses in cotton production throughout the world. Association mapping allows identification and localization of the genes controlling drought-related traits which will be helpful in cotton breeding. In the present study, genetic diversity analysis and association mapping of yield and drought traits were performed on a panel of 99 upland cotton genotypes using 177 SSR (simple sequence repeat) markers. Yield parameters and drought tolerance-related traits were evaluated for two seasons under two watering regimes: water-stressed and well-watered. The traits included seed cotton yield (SCY), lint yield (LY), lint percentage (LP), water-use efficiency (WUE), yield potential (YP), yield reduction (YR), yield index (YI), drought sensitivity index (DSI), stress tolerance index (STI), harmonic mean (HM), and geometric mean productivity (GMP). The genotypes with the least change in seed cotton yield under drought stress were Zeta 2, Delcerro, Nazilli 87, and DAK 66/3 which were also the most water-use efficient cultivars. The average genetic diversity of the panel was 0.38. The linkage disequilibrium decayed relatively rapidly at 20–30 cM (r²?≥?0.5). We identified 30 different SSR markers associated with the traits. Fifteen and 23 SSR markers were linked to the traits under well-watered and water-stress conditions, respectively. To our knowledge, most of these quantitative yield and drought tolerance-associated loci were newly identified. The genetic diversity and association mapping results should facilitate the development of drought-tolerant cotton lines with high yield in molecular breeding programs.     

High genetic diversity and differentiation of an extremely narrowly distributed and critically endangered decaploid rose (<Emphasis Type="Italic">Rosa praelucens</Emphasis>): implications for its conservation

Rosa praelucens is a critically endangered decaploid alpine rose with an extremely narrow geographic distribution in Northwestern Yunnan, China. We sampled almost all the extant individuals (527 individuals in 31 natural locations and 56 individuals preserved in three local living collections) to assess the genetic variation and to probe the genetic connectivity among the individuals and populations based on three cpDNA intergenic spacers and six fluorescent amplified fragment length polymorphism (AFLP) markers. The morphological traits from seven populations were also measured. R. praelucens exhibited high levels of morphological variation, genetic diversity, and differentiation. The extant individuals were clustered into eight groups in neighbor-net networks, and subsequent Bayesian analysis assigned them into three larger gene pools, both in accordance with their morphological traits, especially flower color. The living collections embraced two private cpDNA haplotypes and included three out of the species’ total eight AFLP genotypes. Rhizome clonal growth, decaploid, and mixed breeding system may largely contribute to high genetic diversity and differentiation in R. praelucens. We concluded that the endangered status of R. praelucens may mainly be due to habitat fragmentation and loss and inherent reproductive difficulties, rather than low genetic diversity. The populations contributing higher cpDNA genetic diversity, representing more AFLP genotypes, and encompassing private cpDNA haplotypes should be given conservation priority by creating plant-micro reserves. The living collections should also be targeted for further ex situ conservation, population recovery, and reintroduction of R. praelucens plants.