Microsatellite markers from sugarcane (Saccharum spp.) ESTs cross transferable to erianthus and sorghum

Analysis of a sugarcane (Saccharum spp.) EST (expressed sequence tag) library of 8678 sequences revealed approximately 250 microsatellite or simple sequence repeats (SSRs) sequences. A diversity of dinucleotide and trinucleotide SSR repeat motifs were present although most were of the (CGG)_n trinucleotide motif. Primer sets were designed for 35 sequences and tested on five sugarcane genotypes. Twenty-one primer pairs produced a PCR product and 17 pairs were polymorphic. Primer pairs that produced polymorphisms were mainly located in the coding sequence with only a single pair located within the 5′ untranslated region. No primer pairs producing a polymorphic product were found in the 3′ untranslated region. The level of polymorphism (PIC value) in cultivars detected by these SSRs was low in sugarcane (0.23). However, a subset of these markers showed a significantly higher level of polymorphism when applied to progenitor and related genera (Erianthus sp. and Sorghum sp.). By contrast, SSRs isolated from sugarcane genomic libraries amplify more readily, show high levels of polymorphism within sugarcane with a higher PIC value (0.72) but do not transfer to related species or genera well.     

Development of microsatellite markers and characterization of simple sequence length polymorphism (SSLP) in rice (Oryza sativa L.)

Microsatellite markers containing simple sequence repeats (SSR) are a valuable tool for genetic analysis. Our objective is to augment the existing RFLP map of rice with simple sequence length polymorphisms (SSLP). In this study, we describe 20 new microsatellite markers that have been assigned to positions along the rice chromosomes, characterized for their allelic diversity in cultivated and wild rice, and tested for amplification in distantly related species. Our results indicate that the genomic distribution of microsatellites in rice appears to be random, with no obvious bias for, or clustering in particular regions, that mapping results are identical in intersubspecific and interspecific populations, and that amplification in wild relatives ofOryza sativa is reliable in species most closely related to cultivated rice but becomes less successful as the genetic distance increases. Sequence analysis of SSLP alleles in three relatedindica varieties demonstrated the clustering of complex arrays of SSR motifs in a single 300-bp region with independent variation in each. Two microsatellite markers amplified multiple loci that were mapped onto independent rice chromosomes, suggesting the presence of duplicated regions within the rice genome. The availability of increasing numbers of mapped SSLP markers can be expected to increase the power and resolution of genome analysis in rice.     

Sequence characterization of microsatellites in diploid and polyploid Ipomoea

 The objectives of the present study were to evaluate the inheritance and nucleotide sequence profiles of microsatellite genetic markers in hexaploid sweetpotato [Ipomoea batatas (L.) Lam.] and its putative tetraploid and diploid ancestors, and to test possible microsatellite mutation mechanisms in polyploids by direct sequencing of alleles. Sixty three microsatellite loci were isolated from genomic libraries of I. batatas and sequenced. PCR primers were designed and used to characterize microsatellite loci in two hexaploid I. batatas populations, a tetraploid Ipomoea trifida population, and a diploid I. trifida population. Nine out of the sixty three primer pairs tested yielded a clearly discernible, heritable banding pattern; five showed Mendelian segregation. All other primer pairs produced either smeared banding patterns, which could not be scored, or no bands at all in I. batatas. All of the primers which produced discernible banding patterns from I. batatas also amplified products of similar size in tetraploid and diploid I. trifida accessions. The sequence analysis of several alleles in the three species showed differences due to mutations in the repeat regions consistent with small differences in the repeat number. However, in some cases insertions/deletions and base substitutions in the microsatellite flanking regions were responsible for polymorphisms in both polyploid and diploid species. These results provide strong empirical evidence that complex genetic mechanisms are responsible for SSR allelic variation in Ipomoea. Four I. batatas microsatellite loci showed polysomic segregation fitting tetraploid segregation ratios. To our knowledge this is the first report of segregation ratios for microsatellites markers in polyploids. Received: 4 January 1999 / Accepted: 4 January 1999     

Development of a microsatellite framework map providing genome-wide coverage in rice (Oryza sativa L.)

 Ninety-four newly developed microsatellite markers were integrated into existing RFLP framework maps of four rice populations, including two doubled haploid, a recombinant inbred, and an interspecific backcross population. These simple sequence repeats (SSR) were predominantly poly(GA) motifs, targetted because of their abundance in rice. They were isolated from a previously described sheared library and a newly constructed enzyme-digested library. Differences in the average length of poly(GA) tracts were observed for clones isolated from the two libraries. The length of GA motifs averaged 21 repeat units for clones isolated from the Tsp-509-digested library, while motifs averaged 17 units for clones from the sheared library. There was no evidence of clustering of microsatellite markers near centromeres or telomeres. Mapping of the 94 newly developed markers as well as of 27 previously reported microsatellites provided genome-wide coverage of the 12 chromosomes, with an average distance of 1 SSLP (simple sequence repeat polymorphism) per 16–20 cM. Received: 13 February 1997/Accepted: 28 February 1997     

Diversity of microsatellites derived from genomic libraries and GenBank sequences in rice (Oryza sativa L.)

The growing number of rice microsatellite markers warrants a comprehensive comparison of allelic variability between the markers developed using different methods, with various sequence repeat motifs, and from coding and non-coding portions of the genome. We have performed such a comparison over a set of 323 microsatellite markers; 194 were derived from genomic library screening and 129 were derived from the analysis of rice-expressed sequence tags (ESTs) available in public DNA databases. We have evaluated the frequency of polymorphism between parental pairs of six inter- subspecific crosses and one inter-specific cross widely used for mapping in rice. Microsatellites derived from genomic libraries detected a higher level of polymorphism than those derived from ESTs contained in the GenBank database (83.8% versus 54.0%). Similarly, the other measures of genetic variability [the number of alleles per locus, polymorphism information content (PIC), and allele size ranges] were all higher in genomic library-derived microsatellites than in their EST-database counterparts. The highest overall degree of genetic diversity was seen in GA-containing microsatellites of genomic library origin, while the most conserved markers contained CCG- or CAG-trinucleotide motifs and were developed from GenBank sequences. Preferential location of specific motifs in coding versus non-coding regions of known genes was related to observed levels of microsatellite diversity. A strong positive correlation was observed between the maximum length of a microsatellite motif and the standard deviation of the molecular-weight of amplified fragments. The reliability of molecular weight standard deviation (SDmw) as an indicator of genetic variability of microsatellite loci is discussed. Received: 5 May 1999 / Accepted: 16 August 1999     

镍离子对尖叶莴苣氮素吸收和生长生理的影响

以‘全年油麦菜’尖叶莴苣为试验材料,采用水培方式,研究3个浓度(0 mg·L^-1、0.1 mg·L^-1、1 mg·L^-1)Ni²⁺在22.4 mg·L^-1 N处理下对尖叶莴苣氮素吸收的生长及生理影响。结果显示：(1)尖叶莴苣根系和地上部生物量随处理时间的增加呈上升趋势。与对照T₁(0 mg·L^-1 Ni²⁺、112 mg·L^-1 N)相比,T₂处理(0 mg·L^-1 Ni²⁺、22.4 mg·L^-1 N)对尖叶莴苣根系及叶片生长具有一定抑制作用,植株鲜重、干重、根冠比、根系长度、平均直径、表面积、体积、根尖数、分根数、叶片表面积和体积在T3处理(0.1 mg·L^-1 Ni²⁺、22.4 mg·L^-1 N)下显著高于对照,T₄处理(1 mg·L^-1 Ni²⁺、22.4 mg·L^-1 N)对尖叶莴苣根系及其叶片生长具有一定促进作用,但对其根尖数和分根数表现出一定抑制性。(2)随着Ni²⁺浓度的增加,尖叶莴苣叶片叶绿素a、叶绿素b和总叶绿素含量呈先升后降的变化规律,且均在T₃处理下显著提高。(3)随着处理时间的增加,尖叶莴苣叶片的净光合速率(P_n)、气孔导度(G_s)和蒸腾速率(T_r)逐渐上升,胞间CO₂浓度(C_i)逐渐下降,且T₃处理叶片的G_s显著高于对照,其C_i最低,P_n最大。(4)施加Ni²⁺对尖叶莴苣有机酸、可溶性蛋白和可溶性糖含量以及SOD和POD活性有显著影响,在T₃处理下有机酸含量降低,可溶性糖和可溶性蛋白含量显著增加,SOD和POD活性显著提高。(5)T₃处理尖叶莴苣根系中N及叶片中B和Ca含量较高;根系中Ni含量高于叶片,T₃处理叶片中的Ni含量较低,Mg含量较高;植株体内Cu含量随Ni²⁺浓度增加而下降。研究表明,外源Ni²⁺处理能影响低氮条件下(22.4 mg·L^-1 N)尖叶莴苣幼苗生长及生理状况,适宜浓度(0.1 mg·L^-1)Ni²⁺可有效提高尖叶莴苣根系对氮素的吸收利用效率,减少氮素施用量,促进尖叶莴苣根系和地上部叶片生长,增加光合色素含量,并提高净光合速率,进而改善植株的产量和营养品质。     

硝态和铵态氮配比对水培油麦菜苗期生长及生理特性的影响

采用水培技术,以油麦菜幼苗为材料,研究不同硝铵态氮配比（NO₃^-∶NH₄⁺）对油麦菜苗期地上部和根系生长及生理特性的影响。结果表明：（1）油麦菜地上部和根系硝酸盐含量皆与营养液中NO₃^--N比例呈正相关关系,且各处理均达到无公害蔬菜的标准。（2）随着营养液中NH₄⁺-N比例的增加,油麦菜地上部有机酸含量先降低后升高,且在硝铵态氮配比为5∶5时最低,可溶性糖含量呈上升趋势,而可溶性蛋白质含量先升高后降低,在硝铵态氮配比为5∶5时最高;油麦菜根系有机酸和可溶性糖含量先升高后降低,两者分别在硝铵态氮配比为5∶5和7.5∶2.5时最高,而可溶性蛋白质含量呈下降趋势,在全NO₃^--N时最高。（3）随着营养液中NH₄⁺-N比例的增加,油麦菜地上部和根系中SOD活性先升后降,并分别在硝铵态氮配比为5∶5和7.5∶2.5时最高,而地上部和根系中MDA、脯氨酸含量和POD、CAT活性的变化趋势则与其相反。（4）随着营养液中NH₄⁺-N比例的增加,油麦菜地上部和根系干重皆先升后降,根冠比则逐渐减小;在硝铵态氮配比为7.5∶2.5时干重最大,根冠比适宜且稳定。研究表明,水培油麦菜苗期地上部和根系生长及生理特性受到氮素形态配比的显著影响,且根系的生理响应更敏感;营养液中硝铵态氮配比为7.5∶2.5时,油麦菜受胁迫程度最低,地上部和根系生长较协调,油麦菜生长和生理状况最佳。     

Development, characterization and mapping of microsatellite markers in Eucalyptus grandis and E. urophylla

 We report on the development, genetic characterization and linkage mapping of a battery of SSR (simple sequence repeat) loci in Eucalyptus grandis and E. urophylla. This study reveals the abundance of SSRs in Eucalyptus, the very high information content of these markers for mapping and individual identification, and demonstrates the feasibility of constructing a comprehensive microsatellite-based linkage map for Eucalyptus. Primer sequence for a set of 20 highly informative EMBRA (Eucalyptus microsatellites from Brazil) loci are made available together with their map position and estimates of the expected heterozygosity and allele size range in these two species. Using genomic library enrichment and anchored-PCR screening prior to sequencing, the efficiency of SSR marker locus development was 63% from sequencing data to operationally useful SSR loci. Absolute transportability between the two species and very high levels of allelic variability and expected heterozygosity (H) were seen at all SSR loci surveyed. The number of alleles per locus ranged from 9 to 26 with an average of 16.3±4.8. The average H of 15 loci was 0.86±0.04, 0.83±0.08 and 0.89±0.04, respectively, for E. urophylla, E. grandis and the combined two-species estimate. In the mapping analysis 16 out of 20 marker loci segregated in a fully informative configuration, allowing the determination of synteny of six homologous linkage groups between the two species. The availability of transportable, multiallelic, PCR-based co-dominant SSR loci represents a dramatic improvement in our ability to carry out detailed population genetic analysis and to search, understand, and manipulate allelic variation at QTLs (quantitative trait loci) in species of Eucalyptus. Received: 16 March 1998 / Accepted: 22 March 1998     

A chlorophyll fluorescence analysis of photosynthetic efficiency, quantum yield and photon energy dissipation in PSII antennae of Lactuca sativa L. leaves exposed to cinnamic acid

This study investigated the effects of cinnamic acid (CA) on growth, biochemical and physiological responses of Lactuca sativa L. CA (0.1, 0.5, 1.0 and 1.5 mM) treatments decreased plant height, root length, leaf and root fresh weight, but it did not affect the leaf water status. CA treatment (1.5 mM) significantly reduced F_v, F_m, photochemical efficiency of PSII (F_v/F_m) and quantum yield of PSII (ΦPSII) photochemistry in L. sativa. The photochemical fluorescence quenching (qP) and non-photochemical quenching (NPQ) were reduced after treatment with 1.5 mM CA. Fraction of photon energy absorbed by PS II antennae trapped by “open” PS II reaction centers (P) was reduced by CA (1.5 mM) while, portion of absorbed photon energy thermally dissipated (D) and photon energy absorbed by PSII antennae and trapped by “closed” PSII reaction centers (E) was increased. Carbon isotope composition ratios (δ¹³C) was less negative (−27.10) in CA (1.5 mM) treated plants as compared to control (−27.61). Carbon isotope discrimination (Δ¹³C) and ratio of intercellular CO₂ concentration (ci/ca) from leaf to air were also less in CA treated plants. CA (1.5 mM) also decreased the leaf protein contents of L. sativa as compared to control.     

Genetic diversity and population structure of the endangered medicinal plant Phellodendron amurense in China revealed by SSR markers

Phellodendron amurense is an endangered tree with important medicinal and economic value in China. In this study, eight nuclear SSR primer pairs were employed to assess the genetic diversity and structure of 22 natural populations, including 516 individuals. A total of 66 alleles were detected with an average of 8.3 alleles per locus ranging from 3 to 17. The expected heterozygosity (He) of each SSR locus varied from 0.347 to 0.877 (average 0.627). Analysis of molecular variance (AMOVA) revealed that the main variation component existed within populations (95.11%) rather than among populations (4.89%). The Wilcoxon's sign-rank tests did not show any recent bottleneck effect in any population. A Mantel test displayed a significant correlation between the geographic distances and genetic distances for all populations (r = 0.566, P = 0.0001), indicating conformity to the isolation by distance model. Bayesian clustering and UPGMA supported grouping the populations into two groups. The present genetic structure of P. amurense may be explained by geographical isolation. The lack of genetic structure and genetic diversity decreased with increasing latitude within the Northeast China group may be due to postglacial northward expansion from a single refugium. Proper conservation measures are proposed for this species.     

Genetic mapping of expressed sequences in onion and in silico comparisons with rice show scant colinearity

The Poales (which include the grasses) and Asparagales [which include onion (Allium cepa L.) and other Allium species] are the two most economically important monocot orders. Enormous genomic resources have been developed for the grasses; however, their applicability to other major monocot groups, such as the Asparagales, is unclear. Expressed sequence tags (ESTs) from onion that showed significant similarities (80% similarity over at least 70% of the sequence) to single positions in the rice genome were selected. One hundred new genetic markers developed from these ESTs were added to the intraspecific map derived from the BYG15-23×AC43 segregating family, producing 14 linkage groups encompassing 1,907 cM at LOD 4. Onion linkage groups were assigned to chromosomes using alien addition lines of Allium fistulosum L. carrying single onion chromosomes. Visual comparisons of genetic linkage in onion with physical linkage in rice revealed scant colinearity; however, short regions of colinearity could be identified. Our results demonstrate that the grasses may not be appropriate genomic models for other major monocot groups such as the Asparagales; this will make it necessary to develop genomic resources for these important plants. Electronic Supplementary Material Supplementary material is available for this article at     

Genetic structure of three Oryza AA genome species (O. rufipogon, O. nivara and O. sativa) as assessed by SSR analysis on the Vientiane Plain of Laos

Microsatellite (SSR) markers can reveal a high level of polymorphic loci, and are increasingly being used in population genetic structure studies. On the Vientiane plain of Laos all components of the rice crop complex exist, wild annual (O. nivara), wild perennial (O. rufipogon) and weedy relatives of rice as well as rice itself. To understand gene flow in the rice complex, the genetic structures of O. rufipogon (10 populations), O. nivara (10 populations) and O. sativa (24 samples) from across the Vientiane Plain, Laos, were compared. Higher genetic differentiation was detected among O. nivara populations (G _ST = 0.77, R _ST = 0.71) than O. rufipogon populations (G _ST = 0.29, R _ST = 0.28), whereas genetic diversity for all populations of these two wild species showed similar values (H _T = 0.77 and 0.64 in O. rufipogon and O. nivara, respectively). Based on neighbor-joining tree constructed on the basis of genetic distance (D _A), three genetic clusters were detected, corresponding to (1) O. sativa samples, (2) O. nivara populations and (3) O. rufipogon populations. Pairwise tests confirmed the genetic differentiation of the three species. Although none of the wild rice individuals used in this study had any cultivated-specific phenotypic traits, genetic admixture analysis detected more than 10% O. sativa membership in three O. rufipogon and one O. nivara populations, indicating that O. sativa alleles may cryptically persist in natural populations of O. rufipogon and O. nivara on the Vientiane Plain.     

Characterization and analysis of microsatellite loci in Elymus caninus (Triticeae: Poaceae)

 Microsatellites are highly variable DNA sequences that can be used as markers for the genetic analysis of plants. The potential of microsatellite markers for use in a genetic diversity study in Elymus species was evaluated. Genomic libraries of Elymus caninus were constructed. The libraries were screened with two dinucleotide, (GA)n and (GT)n, and two trinucleotide repeats, (TCT)n and (CAC)n. A total of 19 positive clones were found for the two dinucleotide repeats; no positive clone was found for the trinucleotide repeats. Positive clones were sequenced to confirm the presence of microsatellites and to generate polymerase chain reaction (PCR) primers based on the sequences flanking the microsatellite. All sequenced (GA)n clones have repeats of n>10; over half of the (GT)n microsatellites have n<10 repeats. Primer pairs were designed and evaluated for 8 selected microsatellites. PCR products were amplified from 15 Elymus caninus accessions. The number of alleles found for the eight loci varied from 1 for ECGA89 and ECGT35 to 13 for ECGA22, as determined by non-denaturing polyacrylamide electrophoresis. Six microsatellite loci were found to be polymorphic in E. caninus. The eight primer pairs were tested on three other species; seven were successful in amplifying DNA from Elymus alaskanus and E. mutabilis, and four amplified DNA from E. caucasicus. Based on these results, microsatellites appear to be useful markers in detecting variation in E. caninus. Received: 8 September 1997/Accepted: 6 October 1997     

Inter-simple sequence repeat (ISSR) amplification for analysis of microsatellite motif frequency and fingerprinting in rice (Oryza sativa L.)

 Inter-simple sequence repeat (ISSR) amplification was used to analyze microsatellite motif frequency in the rice genome and to evaluate genetic diversity among rice cultivars. A total of 32 primers, containing different simple sequence repeat (SSR) motifs, were tested for amplification on a panel of 59 varieties, representative of the diversity of cultivated rice (Oryza sativa L.). The ISSR analysis provided insights into the organization, frequency and levels of polymorphism of different simple sequence repeats in rice. The more common dinucleotide motifs were more amenable to ISSR analysis than the more infrequent tri-, tetra- and penta-nucleotide motifs. The ISSR results suggested that within the dinucleotide class, the poly(GA) motif was more common than the poly(GT) motif and that the frequency and clustering of specific tri- and tetra-nucleotide simple sequence repeats was variable and motif-specific. Furthermore, trinucleotide ISSR markers were found to be less polymorphic than either dinucleotide or certain tetranucleotide ISSR markers, suggesting which motifs would be better targets for microsatellite marker development. The ISSR amplification pattern was used to group the rice genotypes by cluster analysis. These results were compared to surveys of the same varieties for amplified fragment length polymorphism (AFLP), restriction fragment length polymorphism (RFLP) and isozyme markers. The ISSR fingerprint could be used to differentiate the genotypes belonging to either Japonica or Indica sub species of cultivated rice and to dissect finer levels of diversity within each subspecies. A higher percentage of polymorphic bands was produced with the ISSR technique than the AFLP method, based on a similar PCR reaction. Therefore, ISSR amplification proved to be a valuable method for determining genetic variability among rice varieties and for rapidly identifying cultivars. This efficient genetic fingerprinting technique would be useful for characterizing the large numbers of rice accessions held in national and international germplasm centers. Received: 25 May 1998 / Accepted: 17 September 1998     

Molecular diversity and association mapping of fiber quality traits in exotic G. hirsutum L. germplasm

The narrow genetic base of cultivated cotton germplasm is hindering the cotton productivity worldwide. Although potential genetic diversity exists in Gossypium genus, it is largely ‘underutilized’ due to photoperiodism and the lack of innovative tools to overcome such challenges. The application of linkage disequilibrium (LD)-based association mapping is an alternative powerful molecular tool to dissect and exploit the natural genetic diversity conserved within cotton germplasm collections, greatly accelerating still ‘lagging’ cotton marker-assisted selection (MAS) programs. However, the extent of genome-wide linkage disequilibrium (LD) has not been determined in cotton. We report the extent of genome-wide LD and association mapping of fiber quality traits by using a 95 core set of microsatellite markers in a total of 285 exotic Gossypium hirsutum accessions, comprising of 208 landrace stocks and 77 photoperiodic variety accessions. We demonstrated the existence of useful genetic diversity within exotic cotton germplasm. In this germplasm set, 11–12% of SSR loci pairs revealed a significant LD. At the significance threshold (r² ≥ 0.1), a genome-wide average of LD declines within the genetic distance at < 10 cM in the landrace stocks germplasm and > 30 cM in variety germplasm. Genome wide LD at r² ≥ 0.2 was reduced on average to  1–2 cM in the landrace stock germplasm and 6–8 cM in variety germplasm, providing evidence of the potential for association mapping of agronomically important traits in cotton. We observed significant population structure and relatedness in assayed germplasm. Consequently, the application of the mixed liner model (MLM), considering both kinship (K) and population structure (Q) detected between 6% and 13% of SSR markers associated with the main fiber quality traits in cotton. Our results highlight for the first time the feasibility and potential of association mapping, with consideration of the population structure and stratification existing in cotton germplasm resources. The number of SSR markers associated with fiber quality traits in diverse cotton germplasm, which broadly covered many historical meiotic events, should be useful to effectively exploit potentially new genetic variation by using MAS programs.     

菰核基因组SSR引物的开发及其在稻族植物中的通用性检测

利用数据库中已有的部分菰(Zizania latifolia Turcz.)核基因组序列,采用in silico方法开发其SSR引物,并选取我国不同纬度的5个菰野生种群,对合成的64对引物进行筛选。结果显示:64对引物中有15对至少在一个种群中表现出多态性;共发现84个等位基因,每个位点平均有5.6个等位基因。在5个种群中,观察杂合度为0.000~0.941,预期杂合度为0.072~0.625。种群间的基因流(Nm=0.576)水平较低导致了种群间表现出较高的遗传分化(FST=0.432)。进一步对稻族其他物种的通用性检测发现,15个多态位点中,有8个位点在亚洲栽培稻(Oryza sativa L.)中得到扩增,有9个位点在普通野生稻(O.rufipogon Griff.)中得到扩增。     

Cloning, expression and function of phosphate transporter encoded gene in Oryza sativa L.

OsPT6:1, a phosphate transporter encoding gene from the leaf samples of Oryza sativa, was identified through PCR with specifically designed primers. The phylogenetic analysis and the conserved amino acid residue site detection suggested OsPT6:1 a possible high-affinity phosphate transporter encoding gene. In situ hybridization and RT-PCR demonstrated the expression of OsPT6:1 in both roots and leaves. The peak expression signal was observed in mesophyll cells under low phosphorus (P) induction. A homologous recombination study indicated that OsPT6:1 can enhance the Pi uptake efficiency of Pichia pastoris. At the meantime, the introduction of OsPT6:1 was able to complement the Pi uptake function of yeast cells with high-affinity phosphate transporters deficient. Those results substantiated our contention that OsPT6:1 encoded a high-affinity phosphate transporter of Oryza sativa.     

Race structure within the Mesoamerican gene pool of common bean (Phaseolus vulgaris L.) as determined by microsatellite markers

Common bean (Phaseolus vulgaris L.) cultivars are distinguished morphologically, agronomically and ecologically into specific races within each of the two gene pools found for the species (Andean and Mesoamerican). The objective of this study was to describe the race structure of the Mesoamerican gene pool using microsatellite markers. A total of 60 genotypes previously described as pertaining to specific Mesoamerican races as well as two Andean control genotypes were analyzed with 52 markers. A total of 267 bands were generated with an average of 5.1 alleles per marker and 0.297 heterozygosity across all microsatellites. Correspondence analysis identified two major groups equivalent to the Mesoamerica race and a group containing both Durango and Jalisco race genotypes. Two outlying individuals were classified as potentially of the Guatemala race although this race does not have a defined structure and previously classified members of this race were classified with other races. Population structure analysis with K = 1–4 agreed with this classification. The genetic diversity based on Nei’s index for the entire set of genotypes was 0.468 while this was highest for the Durango–Jalisco group (0.414), intermediate for race Mesoamerica (0.340) and low for race Guatemala (0.262). Genetic differentiation (G _ST) between the Mesoamerican races was 0.27 while genetic distance and identity showed race Durango and Jalisco individuals to be closely related with high gene flow (N _m) both between these two races (1.67) and between races Durango and Mesoamerica (1.58). Observed heterozygosity was low in all the races as would be expected for an inbreeding species. The analysis with microsatellite markers identified subgroups, which agreed well with commercial class divisions, and seed size was the main distinguishing factor between the two major groups identified.Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.