Genetic variation of a global germplasm collection of chickpea (<Emphasis Type="Italic">Cicer arietinum</Emphasis> L.) including Italian accessions at risk of genetic erosion

Chickpea (Cicer arietinum L.) is one of the most important legumes worldwide. We addressed this study to the genetic characterization of a germplasm collection from main chickpea growing countries. Several Italian traditional landraces at risk of genetic erosion were included in the analysis. Twenty-two simple sequence repeat (SSR) markers, widely used to explore genetic variation in plants, were selected and yielded 218 different alleles. Structure analysis and hierarchical clustering indicated that a model with three distinct subpopulations best fits the data. The composition of two subpopulations, named K1 and K2, broadly reflects the commercial classification of chickpea in the two types desi and kabuli, respectively. The third subpopulation (K3) is composed by both desi and kabuli genotypes. Italian accessions group both in K2 and K3. Interestingly, this study highlights genetic distance between desi genotypes cultivated in Asia and Ethiopia, which respectively represent the chickpea primary and the secondary centres of diversity. Moreover, European desi are closer to the Ethiopian gene pool. Overall, this study will be of importance for chickpea conservation genetics and breeding, which is limited by the poor characterization of germplasm collection.     

Investigation of Chinese Wolfberry (<Emphasis Type="Italic">Lycium</Emphasis> spp.) Germplasm by Restriction Site-Associated DNA Sequencing (RAD-seq)

Chinese wolfberry (Lycium spp.) is an important edible and medicinal plant, with a long cultivation history. The genetic relationships among wild Lycium species and landraces have been unclear for a number of reasons, which has hindered the breeding of modern Chinese wolfberry cultivars. In this study, we collected 19 accessions of Chinese wolfberry germplasm, and constructed the genetic relationship based on RAD-seq markers. We obtained 30.32 Gb of clean data, with the average value of each sample being 1.596 Gb. The average mapping rate was 85.7%, and the average coverage depth was 6.76 X. The phylogeny results distinguished all accessions clearly. All the studied landraces shared their most recent common ancestor with L. barbarum, which indicated that L. barbarum may be involved in cultivation of these landraces. The relationship of some landraces, namely the ‘Ningqi’ series, ‘Qingqi-1’ and ‘Mengqi-1,’ has been supported by the phylogeny results, while the triploid wolfberry was shown to be based on a hybrid between ‘Ningqi-1’ and a tetraploid wolfberry. This study uncovered the genetic background of Chinese wolfberry, and developed the foundation for species classification, accession identification and protection, and the production of hybrid cultivars of wolfberry.     

<Emphasis Type="Italic">Miscanthus</Emphasis> sp.: Genetic Diversity and Phylogeny in China

Miscanthus genetic resources are widely distributed throughout China. However, genetic studies on Miscanthus lagged far behind other crops (e.g., sorghum, maize). To establish the comprehensive genetics knowledge of Miscnathus in China, here we report the genetic and phylogenetic diversity of 174 domestic Miscanthus accessions, along with an external Miscanthus × giganteus control. Cytological observations and flow cytometry analyses indicated that there were two major Miscanthus cytotypes in China: diploid (86.86%) and tetraploid (12.57%) without triploid. A total of 108 polymorphic loci generated from 25 SSR primers were used to evaluate the genetic variation. Large variations in genetic similarity coefficients (GSCs), ranging from 0.08 to 0.97 with a mean value of 0.39, were observed between these Miscanthus accessions. Our phylogenetic data revealed that these accessions were clustered into four main clades: M. section Miscanthus, M. section Diandranthus, M. section Triarrhena, and hybrids. The average percentage of polymorphic loci (P), gene diversity (H), and Shannon’s diversity index (I) among Miscanthus species are 70.93%, 0.22, and 0.34, respectively. These were consistent with the analysis of molecular variance (AMOVA) results, showing that 85% of genetic variation was found within clades. This study investigated the clear phylogenetic relationship of Miscanthus species in China, which will be valuable for further utilization of the germplasm in genetic improvement and hybrid breeding of Miscanthus.     

Complete chloroplast genome sequence of <Emphasis Type="Italic">Capsicum</Emphasis><Emphasis Type="Italic">baccatum</Emphasis> var. <Emphasis Type="Italic">baccatum</Emphasis>

Molecular markers derived from the complete chloroplast genome can provide effective tools for species identification and phylogenetic resolution. Complete chloroplast (cp) genome sequences of Capsicum species have been reported. We herein report the complete chloroplast genome sequence of Capsicum baccatum var. baccatum, a wild Capsicum species. The total length of the chloroplast genome is 157,145 bp with 37.7 % overall GC content. One pair of inverted repeats, 25,910 bp in length, was separated by a small single-copy region (17,974 bp) and large single-copy region (87,351 bp). This region contains 86 protein-coding genes, 30 tRNA genes, 4 rRNA genes, and 11 genes contain one or two introns. Pair-wise alignments of chloroplast genome were performed for genome-wide comparison. Analysis revealed a total of 134 simple sequence repeat (SSR) motifs and 282 insertions or deletions variants in the C. baccatum var. baccatum cp genome. The types and abundances of repeat units in Capsicum species were relatively conserved, and these loci could be used in future studies to investigate and conserve the genetic diversity of the Capsicum species.     

Assessment of Genetic Diversity in <Emphasis Type="Italic">Secale cereale</Emphasis> Based on SSR Markers

The primary aim of this study was to estimate genetic diversity among Secale cereale L. accessions using 22 previously published simple sequence repeat (SSR) markers. The plant material included 367 rye accessions comprising historical and contemporary cultivars, cultivated materials, landraces, and breeding strains from the Polish breeding company Danko. The studied accessions represented a wide geographical diversity. Several methods were employed to analyze genetic diversity among the Secale cereale L. accessions and to determine population structure: principal coordinate analysis (PCoA), neighbor-joining (NJ), and Bayesian clustering. We also defined a core collection of 25 rye accessions representing over 93 % of SSR alleles. The results of these analyses showed that accessions from the rye gene bank are clearly divergent in comparison with materials received directly from European breeding companies. Our findings suggest also that the genetic pool of current rye cultivars is becoming narrower during breeding processes. The selected panel of SSR markers performed well in detection of genetic diversity patterns and can be recommended for future germplasm characterization studies in rye.     

Genetic diversity and population structure of pummelo (<Emphasis Type="Italic">Citrus maxima</Emphasis>) germplasm in China

Pummelo (Citrus maxima) is one of the basic species of Citrus. It has been cultivated for about 4000 years in China, and therefore, there are abundant germplasm during the long time of culture. However, there is still a lack of a detailed study of the genetic characteristics of pummelo population. In this study, genetic diversity and population structure among 274 pummelo accessions collected in China were analyzed using 31 nuclear simple sequence repeat (nSSR) markers. The observed heterozygosity was calculated as 0.325 and genetic differentiation Fst as 0.077. Genetic structure analysis divided the whole germplasm into three subpopulations, Pop-a, Pop-b, and Pop-c. Pop-a was composed of accessions mostly from Southeast China, Pop-b was composed of accessions from the central region of South China, and Pop-c was composed of accessions from Southwest China. Meanwhile, the analysis of principal coordinate analysis and neighbor-join tree supported the viewpoint of three subpopulations, and then the possible dispersal routes of pummelos in China were proposed. This study provides an insight into the genetic diversity, facilitates future genome-wide association studies, and promotes the breeding program of pummelo as well.     

The southwestern origin and eastward dispersal of pear (<Emphasis Type="Italic">Pyrus pyrifolia</Emphasis>) in East Asia revealed by comprehensive genetic structure analysis with SSR markers

Pyrus pyrifolia is considered one of the most important cultivated Pyrus species. Hundreds of landraces and bred cultivars have been developed through the natural and artificial hybridizations necessary due to self-incompatibility. In this study, the genetic diversity of 478 Pyrus accessions, including Chinese landraces, bred cultivars, and wild samples, as well as introduced pear cultivars from Japan and Korea, was investigated with a set of 17 simple sequence repeat (SSR) markers distributed across all 17 linkage groups of the pear genome. A total of 121 alleles were detected, including 4 rare alleles with a frequency lower than 5%. Diversity statistics indicated a high level of genetic variation as quantified by the average values of the observed heterozygosity, the expected heterozygosity, and Wright’s fixation index, at 0.76, 0.78, and 0.02, respectively. Population structure and discriminant analysis of principal component analysis implied extensive genetic communication between sand pears in China and revealed four contiguous geographical clusters with overlapping geographical regions. The diversity of the four clusters and approximate Bayesian computation (ABC) indicated that sand pear spread from west to east along the Pearl River and Yangtze River valleys. High diversity and polyphyletic genetic components of cultivars in southwestern China further support southwestern China as the probable center of divergence for Pyrus species. A core collection of 80 out of 470 cultivars was selected, accounting for about 17% of accessions, and capturing 91% of all alleles, including all rare alleles. Our research provides a comprehensive understanding of sand pear germplasm in East Asia and constructs a preliminary core collection, which will be useful for association genetics studies, germplasm conservation, and breeding programs.     

Genetic structure of cultivated <Emphasis Type="Italic">Zanthoxylum</Emphasis> species investigated with SSR markers

Zanthoxylum is an economically and ecologically important genus of the Rutaceae family, of which Z. bungeanum and Z. armatum have a long history of cultivation in China. However, how the natural processes such as selection and drift and agriculture practices have influenced the genetic variation of cultivated Zanthoxylum species during long-term domestication remains elusive. Herein, we determined the population genetic structure of current widely cultivated Zanthoxylum species, Z. bungeanum and Z. armatum. Microsatellite markers revealed a high level of genetic variation and significant genetic differentiation for both species despite Z. bungeanum showed higher genetic diversity than Z. armatum. AMOVA indicated that most of the genetic variation exists within individuals rather than among provenances for both species. Population structure analyses generated three distinct groups within the entire accessions. All Z. bungeanum accessions were distinguished into two major geographic groups, north and south groups, with Qinling Mountains as the main geographic barrier to gene flow while a significant genetic differentiation was observed between cultivated and wild Z. armatum accessions. Mantel test of Z. bungeanum displayed a significant correlation between genetic and geographic distances within each inferred group but no correlation between genetic and geographic distance was observed when comparing genetic and geographic distances focusing only on pairwise of north vs. south provenances, ruling out the hypothesis that gene flow between north and south provenances followed an isolation-by-distance model. Our research provided a fundamental genetic profile that will improve the conservation and responsible exploitation of the extant germplasm of Zanthoxylum.     

Analysis of microsatellite loci variability in rare and endemic species <Emphasis Type="Italic">Allium regelianum</Emphasis> A.K. Becker ex Iljin

SSR analysis of rare and endemic species Allium regelianum, which grows in the south of Russia, was performed for the first time. Variability analysis of 88 accessions of A. regelianum was carried out using four highly polymorphic microsatellite loci (PIC value ranged from 0.55 to 0.72). SSR-analysis made it possible to revealed polymorphism within and among the populations of A. regelianum from Volgograd region. Analysis of Wright’s F-statistics and the analysis of molecular variance showed that more than 90% of total genetic variation of the species was attributed to the differences within the populations and less than 10% of the differences were found among the populations. Cluster analysis of 46 accessions of A. regelianum from three populations of Volgograd region and principal coordinate analysis of all accessions did not reveal differentiation among the populations.     

Development of cost-effective single nucleotide polymorphism marker assays for genetic diversity analysis in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Competitive allele-specific PCR (KASPar) assay is a user-friendly system that provides flexibility in the numbers of single nucleotide polymorphisms (SNPs) and genotypes. Based on Illumina-GA-IIx genomic data from 10 genotypes with a broad genetic background, 3183 SNPs were selected for KASPar assays development, and 568 were finally converted and selected for Brassica rapa germplasm characterization (17.8%) on the basis of reproducibility, missing data rate, and uniform genetic distribution. High levels of polymorphism of these markers across 231 B. rapa genotypes were verified, illustrating by high polymorphic information content (averaged 0.34), minor allele frequency (0.37), genetic diversity (0.45), and the low observed heterozygosity (0.10). Based on the SNP dataset, structure and principal coordinates analysis, and neighbor-joining phylogenetic methods were used to examine the population structure and gave highly consistent results. The 231 accessions were divided into the four primary subspecies, representing 99 accessions from B. rapa ssp. pekinensis, 85 from B. rapa ssp. chinensis, 30 from B. rapa ssp. rapifera, and 17 from B. rapa ssp. oleifera and were further subdivided into 12 lower-order clusters according to different morphotypes. The genetic variability and pairwise fixation index analysis revealed that the ssp. pekinensis accessions possess the most extensive genetic variation among the four subspecies. The KASPar system is highly useful for validating SNPs and will be valuable for genetics research and breeding applications in B. rapa.     

Microsatellite based linkage disequilibrium analyses reveal <Emphasis Type="Italic">Saltol</Emphasis> haplotype fragmentation and identify novel QTLs for seedling stage salinity tolerance in rice (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.)

A set of 84 diverse rice genotypes were assessed for seedling stage salt tolerance and their genetic diversity using 41 polymorphic SSR markers comprising of 19 Saltol QTL linked and 22 random markers. Phenotypic screening under hydroponics identified three indica landraces (Badami, Shah Pasand and Pechi Badam), two Oryza rufipogon accessions (NKSWR2 and NKSWR17) and one each of Basmati rice (Seond Basmati) and japonica cultivars (Tompha Khau) as salt tolerant, having similar tolerance as of Pokkali and FL478. Among the salt tolerant genotypes, biomass showed positive correlation with shoot fresh weight and negative association with root and shoot Na⁺ content. The results indicated repression of Na⁺ loading within the tolerant plants. Linkage disequilibrium (LD) of the Saltol linked markers was weak, suggestive of high fragmentation of Pokkali haplotype, a result of evolutionary active recombination events. Poor haplotype structure of the Saltol region, may reduce its usefulness in marker assisted breeding programmes, if the target foreground markers chosen are wide apart. LD mapping identified eight robust marker-trait associations (QTLs), of which RM10927 was found linked to root and shoot Na⁺ content and RM10871 with shoot Na⁺/K⁺ ratio. RM271 on chromosome 10, an extra Saltol marker, was found associated to root Na⁺/K⁺ ratio. This marker showed a distinct allele among O. rufipogon accessions. There were also other novel loci detected on chromosomes 2, 5 and 10 influencing salt tolerance in the tested germplasm. Although Saltol remained as the key locus, the role of other genomic regions cannot be neglected in tailoring seedling stage salt tolerance in rice.     

Assessment of genetic diversity in <Emphasis Type="Italic">Mucuna</Emphasis> species of India using randomly amplified polymorphic DNA and inter simple sequence repeat markers

Genus Mucuna which is native to China and Eastern India comprises of perennial climbing legume with long slender branches, trifoliate leaves and bear green or brown pod covered with soft or rigid hairs that cause intense irritation. The plants of this genus are agronomically and economically important and commercially cultivated in India, China and other regions of the world. The high degrees of taxonomical confusions exist in Mucuna species that make authentic identification and classification difficult. In the present study, the genetic diversity among the 59 accessions of six species and three varieties of M. pruriens has been assessed using DNA fingerprinting based molecular markers techniques namely randomly amplified polymorphic DNA (RAPD), inter simple sequence repeats (ISSR) and combined dataset of RAPD and ISSR. Also, genetic relationship among two endemic species of Mucuna namely M. imbricata and M. macrocarpa and two varieties namely IIHR hybrid (MHR) and Dhanwantari (MD) with other species under study was investigated by using cluster analysis and principal coordinate analysis. The cluster analysis of RAPD, ISSR and combined dataset of RAPD and ISSR clearly demonstrated the existence of high interspecific variation than intra-specific variation in genus Mucuna. The utility and efficacy of RAPD and ISSR for the study of intra species and interspecies genetic diversity was evident from AMOVA and PCoA analysis. This study demonstrates the genetic diversity in Mucuna species and indicates that these markers could be successfully used to assess genetic variation among the accessions of Mucuna species.     

Assessment of genetic diversity among okra genotypes using SSR markers

Okra (Abelmoschus esculentus) is an important nutritious vegetable. Despite its high economic and industrial value, very little attention has been paid to assess genetic diversity of okra at molecular level. For effective conservation and proper deployment of germplasm, a study on diversity analysis of okra germplasm was conducted with DNA markers. Microsatellite/Simple sequence repeat (SSR) markers were utilized to evaluate the genetic diversity among 96 accessions of Abelmoschus, of which 92 accessions were of A. esculentus and one accession each of A. tuberculatus, A. moschatus, A. moschatus subspecies tuberosus and A. manihot. A set of 40 SSR primers were tested, of which 30 primers gave reproducible amplification which were used further for diversity analysis. With a mean of 7.1 bands per SSR, DNA amplification with 30 SSRs generated a total 213 bands, of which 60.66 % were recorded polymorphic. Polymorphic information content ranged between 0.11 and 0.80 with an average of 0.52, indicating that the majority of primers were informative. The Jaccard’s coefficient ranged from 0.107 to 0.969. The UPGMA analysis grouped Abelmoschus genotypes into three main clusters at a cut-off of 0.20. Results of present study revealed that sufficient variation exists among the studied accessions and GAO-5 which was found highly diverse can be exploited for okra improvement. The outcome of present research would assist to make use of Ablemoschus germplasm for okra breeding.     

Screening of rice landraces (<Emphasis Type="Italic">Oryza sativa</Emphasis> L.) for seedling stage salinity tolerance using morpho-physiological and molecular markers

Traditional rice landraces of coastal area in Bangladesh are distinct regarding their phenotype, response to salt stress and yield attributes. With characterization of these landraces, suitable candidate genes for salinity tolerance could be identified to introgress into modern rice varieties. Therefore, the aim of this experiment was to uncover prospective rice landraces tolerant to salinity. Relying on morphological, biochemical and molecular parameters 25 rice genotypes were tested for salt tolerance at germination and seedling stage. At germination stage 0 and 12 dSm^?1 salinity were imposed on rice genotypes. Ward’s cluster analysis divided rice genotypes into three clusters (susceptible, moderately tolerant and tolerant) based on the physiological indices. The tolerant rice landraces to salinity were Sona Toly, Nakraji and Komol Bhog. At seedling stage screening was performed following IRRI standard protocol at 12 dSm^?1 salinity level. Based on all morphological and biochemical parameters Komol Bhog was identified as the highly salinity tolerant landrace while Bolonga, Sona Toly, Dud Sail, Tal Mugur and Nakraji were found as tolerant to salinity. Molecular characterization using two simple sequence repeats (SSR) markers, viz. RM121 and RM337 displayed Bolonga, Til Kapor, Panbra, Sona Toly, Bina Sail, Komol Bhog, Nakraji, Tilkapur, Gajor Goria and Gota were tolerant landraces through genetic similarity in dendrogram. These identified salt-resistant landraces can be used as promising germplasm resources for breeding salt-tolerant high-yielding rice varieties in future.     

Molecular genetic diversity and association mapping of morphine content and agronomic traits in Turkish opium poppy (<Emphasis Type="Italic">Papaver somniferum</Emphasis>) germplasm

As the sole plant source of many potent alkaloids, opium poppy (Papaver somniferum L.) is an important medicinal crop. Nevertheless, few studies have characterized opium poppy germplasm with crop-specific molecular markers. Because Turkey is a diversity center for opium poppy, Turkish germplasm is a valuable genetic resource for association mapping studies aimed at identifying QTLs controlling morphine content and agronomic traits. In this study, the morphological diversity and molecular diversity of 103 Turkish opium poppy landraces and 15 cultivars were analyzed. Potentially useful morphological variation was observed for morphine content, plant height, and capsule index. However, the landraces exhibited limited breeding potential for stigma number, and seed and straw yields. Both morphological and molecular analyses showed distinct clustering of cultivars and landraces. In addition, a total of 164 SSR and 367 AFLP polymorphic loci were applied to an opium poppy association mapping panel composed of 95 opium poppy landraces which were grown for two seasons. One SSR and three AFLP loci were found to be significantly associated with morphine content (P < 0.01 and LD value (r ²) = 0.10–0.32), and six SSR and 14 AFLP loci were significantly associated with five agronomic traits (plant height, stigma number, capsule index, and seed and straw yields) (P < 0.01 and LD value (r ²) = 0.08–0.35). This is the first report of association mapping in this crop. The identified markers provide initial information for marker-assisted selection of important traits in opium poppy breeding.     

Natural variation and genetic make-up of leaf blade area in spring barley

Key message

GWAS analysis for leaf blade area (LA) revealed intriguing genomic regions associated with putatively novel QTL and known plant stature-related phytohormone and sugar-related genes.

Abstract

Despite long-standing studies in the morpho-physiological characters of leaf blade area (LA) in cereal crops, advanced genetic studies to explore its natural variation are lacking. The importance of modifying LA in improving cereal grain yield and the genes controlling leaf traits have been well studied in rice but not in temperate cereals. To better understand the natural genetic variation of LA at four developmental stages, main culm LA was measured from 215 worldwide spring barleys including 92 photoperiod-sensitive accessions [PHOTOPERIOD RESPONSE LOCUS 1 (Ppd-H1)] and 123 accessions with reduced photoperiod sensitivity (ppd-H1) locus under controlled greenhouse conditions (long-day; 16/8 h; ~ 20/~ 16 °C day/night). The LA of Ppd-H1-carrying accessions was always smaller than in ppd-H1-carrying accessions. We found that nine SNPs from the Ppd-H1 gene were present in the collection of which marker 9 (M9; G/T in the CCT-domain) showed the most significant and consistent effect on LA at all studied developmental stages. Genome-wide association scans (GWAS) showed that the accessions carrying the ppd-H1 allele T/M9 (late heading) possessed more genetic variation in LA than the Ppd-H1 group carrying G/M9 (early heading). Several QTL with major effects on LA variation were found close to plant stature-related heading time, phytohormone- and sugar-related genes. The results provide evidence that natural variation of LA is an important source for improving grain yield, adaptation and canopy architecture of temperate cereals.

     

Simple Sequence Repeat and <Emphasis Type="Italic">S</Emphasis>-locus Genotyping to Explore Genetic Variability in Polyploid <Emphasis Type="Italic">Prunus spinosa</Emphasis> and <Emphasis Type="Italic">P. insititia</Emphasis>

Polyploid Prunus spinosa (2n = 4×) and P. insititia (2n = 6×) represent enormous genetic potential in Central Europe, which can be exploited in breeding programmes. In Hungary, 17 cultivar candidates were selected from wild-growing populations including 10 P. spinosa, 4 P. insititia and three P. spinosa × P. domestica hybrids (2n = 5×). Their taxonomic classification was based on their phenotypic characteristics. Six simple sequence repeats (SSRs) and the multiallelic S-locus genotyping were used to characterize genetic variability and reliable identification of the tested accessions. A total of 98 SSR alleles were identified, which presents 19.5 average allele number per locus, and each of the 17 genotypes could be discriminated based on unique SSR fingerprints. A total of 23 S-RNase alleles were identified. The complete and partial S-genotype was determined for 8 and 9 accessions, respectively. The identification of a cross-incompatible pair of cultivar candidates and several semi-compatible combinations help maximize fruit set in commercial orchards. Our results indicate that the S-allele pools of wild-growing P. spinosa and P. insititia are overlapping in Hungary. A phylogenetic and principal component analysis confirmed the high level of diversity and genetic differentiation present within the analysed genotypes and helped clarify doubtful taxonomic identities. Our data confirm that S-locus genotyping is suitable for diversity studies in polyploid Prunus species. The analysed accessions represent huge genetic potential that can be exploited in commercial cultivation.     

Genetic diversity of the African poplar (<Emphasis Type="Italic">Populus ilicifolia</Emphasis>) populations in Kenya

We evaluated the genetic diversity of the African poplar (Populus ilicifolia) populations found in Kenya compared with reference samples of five poplar species from North America and one species introduced in Kenya from India (KEFRI-Kenya). Amplified fragment length polymorphism (AFLP) was used with the objective of providing important information for breeding and in situ/ex situ conservation of this species. Samples collected from three locations along the species’ natural range (Athi, Ewaso Nyiro, and Tana rivers) were compared with four samples of locally planted Populus deltoides stand introduced from India and ten reference samples from North America. Six AFLP primer combinations produced 521 clear bands for analysis. The percentage polymorphic loci were lowest in Tana (20.4 %) and highest in Athi (40.6 %). The average heterozygosity across the studied populations was between 0.07 and 0.3. AMOVA revealed more genetic variation partitioning within population (87 %; P?<?0.01) than among populations (13 %; P?<?0.01) suggesting significant genetic variation between populations. Further, UPGMA delineation showed two clusters of the Tana, Athi, and Ewaso Nyiro populations clustered together compared to the North America and India/KEFRI reference samples. Moreover, the study showed that the Athi population is more diverse than those of Tana and Ewaso Nyiro and may be important for conservation, domestication, and improvement studies. The genetic differentiation (F _ST ?=?0.134) among Kenyan P. ilicifolia populations suggests limited possibility of gene flow between these populations.