Skip to main content

Characterisation of sugar beet (Beta vulgaris L. ssp. vulgaris) varieties using microsatellite markers



Sugar beet is an obligate outcrossing species. Varieties consist of mixtures of plants from various parental combinations. As the number of informative morphological characteristics is limited, this leads to some problems in variety registration research.


We have developed 25 new microsatellite markers for sugar beet. A selection of 12 markers with high quality patterns was used to characterise 40 diploid and triploid varieties. For each variety 30 individual plants were genotyped. The markers amplified 3-21 different alleles. Varieties had up to 7 different alleles at one marker locus. All varieties could be distinguished. For the diploid varieties, the expected heterozygosity ranged from 0.458 to 0.744. The average inbreeding coefficient Fis was 0.282 ± 0.124, but it varied widely among marker loci, from Fis = +0.876 (heterozygote deficiency) to Fis = -0.350 (excess of heterozygotes). The genetic differentiation among diploid varieties was relatively constant among markers (Fst = 0.232 ± 0.027). Among triploid varieties the genetic differentiation was much lower (Fst = 0.100 ± 0.010). The overall genetic differentiation between diploid and triploid varieties was Fst = 0.133 across all loci. Part of this differentiation may coincide with the differentiation among breeders' gene pools, which was Fst = 0.063.


Based on a combination of scores for individual plants all varieties can be distinguished using the 12 markers developed here. The markers may also be used for mapping and in molecular breeding. In addition, they may be employed in studying gene flow from crop to wild populations.


Sugar beet (Beta vulgaris L.) is a crop of major importance for sugar production in temperate zones. Varieties are produced through crosses of diploid male sterile (CMS) lines with tetraploid, or increasingly, diploid pollinator lines, resulting in triploid or diploid varieties, respectively [1]. As the parental lines are mixtures of genotypes, the varieties will consist of mixtures of plants from various parental combinations. This leads to some problems in variety registration research. Variety registration is based on Distinctiveness, Uniformity, and Stability (DUS) research. Using a visual inspection of morphological characteristics, distinctiveness from other varieties is not easy to assess, for several reasons: the crop has a narrow genetic basis [2, 3], which results in varieties that are highly similar in appearance [4], the varieties are mixtures of genotypes, and breeders change the pollinator line in modern hybrids frequently to produce locally adapted hybrid varieties. For these reasons, the other two aspects of the standard DUS research, uniformity and stability, are not determined, and there are no UPOV (International Union for the Protection of New Varieties of Plants) guidelines for this crop.

The number of informative morphological characteristics is limited. Therefore, most often production-related characteristics as beet yield, sugar content and total sugar yield are included as descriptors. A preliminary characterisation ("pre-screening") of newly submitted varieties with molecular markers during the winter before sowing could be of help in the planning of the field trials and may give a first indication for distinctiveness, provided that a sufficient number of markers is used and that overall marker profile and phenotype correlate well.

Molecular markers have been used successfully for variety identification in a large number of crops, including selfing species [5, 6] and clonally propagated plants [7, 8]. In sugar beet, RFLP, RAPD, and AFLP [915] studies have been reported. Although AFLP markers are reproducible between laboratories [16, 17], data base building can be a problem as different equipment may lead to different profiles. Six co-dominant microsatellites were used to study genetic diversity in wild, cultivated, and weedy forms of Beta vulgaris[18, 19]. Rae et al. [20] developed a set of mostly dinucleotide repeats for incorporation into the linkage map of B. vulgaris, and Richards et al. [21] characterized eight new polymorphic microsatellite markers, of which five were based on trinucleotide repeats. Cureton et al. [22] developed six microsatellite markers to measure gene flow in sea beet (Beta vulgaris ssp. maritima). Laurent et al. [23] mapped a large number of genomic and EST-derived microsatellites on a genetic map of sugar beet, the majority of which were dinucleotide repeat markers.

To be useful for identification of varieties the markers should allow determining unequivocally the genotype of each plant independently. The ease and accuracy of scoring varies among microsatellite markers, with significantly more problems when applying dinucleotide repeats, due to their tendency to generate more stutter bands, which may co-migrate with neighbouring alleles. The experience in those species in which large replication studies have been set up among laboratories, is that rigorous screening of markers is necessary [5, 6]. For that reason we have developed a set of new microsatellite markers for B. vulgaris with PIG-tailed reverse primers [24] and stringent quality demands (Quality 1 or 2 of Smulders et al. [25]). We have applied this set of 12 di-, tri- and tetranucleotide repeat microsatellite markers to determine the genetic variation within and between 40 diploid and triploid varieties. Using 30 plants per variety we have generated a dataset of genotypes of 1200 plants. We analysed the data with respect to allelic diversity, and discuss applications of the markers in sugar beet, sea beets, and ruderal beets.


Microsatellite marker development

For accurate genotyping of varieties and database building, high quality microsatellite markers are needed. Therefore the isolation of microsatellites was focussed on tri- and tetranucleotide repeats, although dinucleotide repeats were isolated as well. In total 3200 clones were screened for microsatellite-containing inserts. In total 31 clones (1%) were found positive for tetranucleotide motives, 240 (7.7%) for trinucleotide repeats and 240 (7.7%) for dinucleotide repeats. For 65 unique microsatellite sequences, primer pairs were designed on the flanking regions. For each locus the amplification pattern was evaluated with respect to pattern quality and degree of polymorphism on a set of individual plants of 10 varieties originating from different breeders. Twenty-five primer pairs (39%) produced polymorphic and simple banding patterns. These primers were selected for further analysis with fluorescent primers on an ABI 3700 using the same set of test varieties. The twelve most robust markers showing no or moderate stutter bands, a low degree of differential amplification, and easy scorability, were used for genotyping the sugar beet varieties (Table 1). These 12 markers consisted of two perfect and four compound dinucleotide repeat loci, five trinucleotide loci, and one locus with both a perfect dinucleotide repeat and a perfect tetranucleotide repeat.

Table 1 Characteristics of 25 newly developed sugar beet microsatellite markers.

Alleles detected

For the evaluation of the markers 30 individual plants per variety were genotyped. Table 2 shows the number of alleles detected for each marker, which varied widely (from 3 to 21), but the effective number of alleles was quite comparable across loci (1.95-3.74; Table 2). In total 91 different alleles were detected. From the number of dropouts in amplification and the positive Fis values we deduced that null alleles may exist. Additional population-genetic parameters of these varieties are listed in Additional file 1.

Table 2 The number of microsatellite alleles detected in 30 individual plants per variety.

Variety characterization based on dominant scoring of alleles

Using the set of 12 marker loci, we found 25-38 different alleles (on average 32.3 per variety) in the 30 plants of a diploid variety and 33-46 (average 39.0) alleles in a triploid variety (Table 2). In general, individual plants from varieties reported to be diploid had only one or two alleles per locus. There were only 15 out of 330 plants from reportedly diploid varieties with three different alleles at one or two loci (5 plants each of Rebecca and Brigitta, 3 of Nemil, 1 each of HI0032 and Fortis). On average, diploid plants had 1.3 alleles per locus. Among plants of the triploid varieties, the average number of alleles per locus was 1.6. Depending on the locus, between 0 (markers bvv17 and bvv21) and 183 (bvv15) plants contained three different alleles at a single locus. Overall, 528 of the 870 plants of these varieties had three different alleles at one or more marker loci, underlining a considerable amount of genetic variation present within these plants.

Triploid varieties are produced from tetraploid males and diploid female plants. While females are always diploid and may be shared between diploid and triploid varieties, the male plants are either diploid or tetraploid and these may form genetically distinct groups. However, tetraploid lines can also easily be made from diploids. When we calculated genetic differentiation between diploid and triploid groups (330 and 870 plants, respectively) Fst = 0.1327 across all loci, ranging from 0.037 for marker bvv30 to 0.2066 for bvv23. Each of these estimates was highly significant (p < 0.001, tested by permutation of individual plants among all varieties).

This differentiation between diploid and triploid varieties could also be the result of the fact that some breeders specialise in diploid varieties, and others in triploids. If so, it would reflect differentiation among breeders rather than between ploidy levels. We therefore also tested the differentiation among breeding companies. Among breeders, we found Fst = 0.0628 ± 0.0092, which is roughly half of the difference between diploid and triploid varieties.

Genetic diversity and differentiation among varieties

A NJ tree was made using the pairwise genetic distances between varieties to visualise the genetic distances among varieties (Figure 1). It shows that the genetic distance is, on average, larger among diploid varieties. For instance, the inner part of the dendrogram contains 17 triploid varieties at relatively small distances from each other. The same pattern is visible in a PCO plot, with the triploid varieties central in the plot and the diploid varieties further from each other (Additional file 2). Triploid varieties have a higher probability of sharing alleles due to the fact that they have more gene copies, hence on average more alleles, which may explain the pattern observed.

Figure 1
figure 1

Neighbour-joining tree based on pairwise genetic distances between sugar beet varieties. The genetic distances were calculated using dominant scoring of alleles. The names of the varieties are followed by their ploidy level: 2 = diploid (2n = 2x), 3 = triploid (2n = 3x).

There is no clear structure in the genetic relatedness of varieties from particular breeding companies in the tree, except that the top branch consists exclusively of nine varieties from KWS (Ariana, Aurelia, KWS8123, KWS9226, Rebecca, Tiara, Brigitta, Lenora, and Madonna).

Overall, Fst = 0.133, but this value was lower among triploid varieties (Fst = 0.100) and much higher among diploid varieties (Fst = 0.232) (Table 3), which is consistent with the pattern observed in the dendrogram. The correlation between the values of individual marker loci for triploid and diploid varieties is relatively poor (R2 = 0.54), suggesting that the gene pool differences between triploid and diploid varieties are not evenly spread across loci.

Table 3 F-statistics of 40 varieties genotyped with 12 microsatellite markers.

The estimate of Fis for the whole dataset was negative for each of the markers (not shown), which is most likely an artefact of the dominant scoring of the markers. In theory, this can influence the Fst estimates as well. For the diploid varieties we were able to estimate the magnitude of this effect through a comparison with an analysis using codominant scoring (assuming two alleles per locus per plant and no null alleles). Table 3 (middle panel) shows that the actual Fis value varies widely among marker loci, from Fis = +0.876 (heterozygote deficiency) to Fis = -0.350 (excess of heterozygotes), with an average of Fis = 0.282 ± 0.124, Table 3). The effect on the estimation of the variation present among varieties (Fst) is limited: Fst averaged across loci is 0.232 for dominant scores (left panel) and 0.271 for codominant scores (middle panel; 17% more). The Fst estimates for most loci are close to this systematic difference of 17%, and the pairwise correlation between the values per locus is R2 = 0.91. This indicates that differentiation among diploid varieties is being estimated comparably using dominant or codominant scores.


We have developed a set of new microsatellite loci for sugar beet, which amplified 2-21 alleles per locus. This is comparable to the 2-11 alleles found by Richards et al. [21] for their microsatellite markers in a set of sugar beet and sea beet plants. Desplanque et al. [18] and Viard et al. [19] found up to 10 alleles for a marker in a single variety. This level of gene diversity does not seem to correspond with the notion of little genetic variation in the crop sugar beet due to a bottleneck during its development from wild beets [1]. The breeding system, which employs separate gene pools for paternal and maternal parents, increases the gene diversity within individual plants, and the habit of working with pools of parental plants, which contain a large amount of genetic diversity [19], may contribute to the fact that 84-92% [26] of the genetic variation of the crop is present within hybrid varieties.

Ploidy level

We have applied 12 of our markers to analyse 30 plants of each of 40 sugar beet varieties.

The markers detected only few (15/330) triploid plants in diploid varieties. The highest frequencies of triploid plants were found for two varieties (5/30 plants each for Brigitta and Rebecca). These plants are probably the result of pollination by tetraploid pollen donors from production fields for other, triploid, varieties in the neighbourhood of the seed production fields of the diploid varieties. In Europe seed production of sugar beet varieties takes place in the South-West and South-East of France, Northern Italy, and the South of Ukraine, and in these areas the distance between production fields is at least 1000 m to severely limit cross-pollination, but this cannot be avoided completely. Accidental cross-fertilization may also take place with ruderal populations in the vicinity of the seed production fields [27, 28], but this would produce diploid offspring.

Genetic differentiation

The overall genetic differentiation between diploid and triploid varieties was Fst = 0.1327 across all loci. Part of this differentiation coincides with the differentiation among breeders' gene pools, which was Fst = 0.0628. This suggests that breeders use parental lines that are, to some extent, genetically different. The latter value can be expected to gradually decrease in the future, as there have been mergers between sugar beet breeding companies in recent years, which may result in merging of the breeding programs.

When partitioning the genetic variation using F statistics, the estimate of Fis of diploid plants turned out to be highly variable among microsatellite loci: from Fis = 0.876 (large shortage of heterozygotes) to Fis = -0.35 (excess of heterozygotes). The excess of heterozygotes is not surprising as the propagation system pairs selected male-sterile (CMS) mother lines with selected father lines, with the aim of assortative mating and hybrid seed production. The shortage of heterozygotes at some marker loci may indicate selection. It may also indicate the presence of null-alleles, i.e. alleles that have gone undetected, or skewed inheritance [12]. Laurent et al. [23] found 14% skewed segregation in an F2 population, notably for markers on linkage group V [29]. Viard et al. [19, 30] found significant heterozygote deficiencies in weed beets. Fénart et al. [1] observed also significant deviations in Fis, in both directions, in wild sea beet and weed beet populations. Viard et al. [19] thought it may be related to a low frequency of self-compatibility alleles commonly used in breeding programs. This was recently confirmed by Arnaud et al. [27].

Nonetheless, Fst values of dominantly scored and codominantly scored markers (for diploid varieties) were in good concordance, indicating that regardless of the statistical analysis of the data, genetically similar and dissimilar varieties can be distinguished reliably. This is in agreement with the conclusions of De Riek et al. [24], who compared the power of these microsatellites with that of a set of AFLP markers. The differentiation among diploid varieties was quite high: Fst ranged from 0.093 to 0.421 (Table 3). The average of 0.232 is higher than Fénart et al. [1]'s estimate of Fst = 0.082 among 13 diploid sugar beet varieties using 5 microsatellite markers, which in turn was higher than the differentiation among weed beets and among sea beets. It would be interesting to determine the level of differentiation assessed with our markers among these groups of beets.


Based on a combination of scores for individual plants all varieties can be distinguished using the 12 markers employed here. However, as the varieties are mixtures of genotypes, not all individual plants can always be identified or classified unequivocally. De Riek et al. [26] compared various ways of analysing the data for eight of these varieties. They concluded that, using the data for 30 individual plants for each variety, assignment methods accomplished a very good distinction among the genetically diverse varieties. In their assignment-based method, for each individual plant the 10 most genetically similar partner plants were identified across the whole data set. The origin of these highest-ranking plants was then used to assign the plants to a particular variety. With microsatellite data, between 24 and 30 of the 30 plants analysed for each variety, were assigned correctly to this variety. The partitioning of the origin of the highest-ranking partners over all varieties in the dataset was also used to develop an assignment-based similarity measure for such sets of mixtures of genotypes, called similarity-by-assignment (Sax, y) [26].


Microsatellite markers may be used for genetic mapping and breeding purposes [29]. The markers developed here were polymorphic within all or nearly all varieties, which indicates that they may be used for mapping in most crosses in sugar beet. In addition, they may be employed in studies of crop-to-wild gene flow [1], including those in the frame of biosafety studies [31].


Plant material

For the isolation of microsatellites, genomic DNA of Beta vulgaris L. ssp. vulgaris variety Holly was used. For the characterization of varieties, 30 individual plants of 40 varieties (listed in Table 2) were analyzed (in total 1200 plants). Young leaves of a single individual were harvested, immediately frozen in liquid nitrogen and stored at -80°C until use.

DNA extraction

For the construction of a genomic library enriched for microsatellites, nuclear DNA of high quality was extracted from leaves of variety Holly according to Vosman et al. [32]. For microsatellite amplification, DNA of single individuals was extracted from freeze-dried leaves either according to Fulton et al. [33] or by a combination of this method with the Qiagen Dneasy Plant Mini kit (Westburg, The Netherlands). In the combination extraction protocol, after chloroform extraction the cleared supernatant was mixed with Qiagen binding buffer (AP3/EtOH) and applied to a DNeasy spin column (Esselink, unpublished). Subsequently, the column was washed and DNA eluted. Typical yield of this extraction protocol was 20 μg DNA per 20 mg dried weight.

Microsatellite isolation

Microsatellites were isolated from enriched small-insert genomic libraries essentially as described by Van de Wiel et al. [34] and Esselink et al. [7]. The DNA was digested with Alu I, Mbo I or Rsa I, and the enrichment was carried out using filter-immobilized synthetic dinucleotide [(GT)12, (GA)12], trinucleotide [(TCT)10, (TGT)9, (GAG)8, (GTG)8, (TGA)9, (AGT)10, (CTG)8, (CGT)8], and tetranucleotide [(GATA)8, (TGTT)8, (GTAT)8] repeats, all in separate reactions. Primers were designed on the obtained sequences using primer3

Microsatellite amplification and detection

Microsatellites were amplified in a 20 μl reaction volume containing 20 ng of genomic DNA, 2-10 pmol of each primer, 100 μM of each dNTP, 10 mM Tris-HCL pH 9.0, 20 mM (NH4)2SO4, 0.01% Tween 20, 1.5 mM MgCl2 and 0.3 Units Goldstar Taq DNA polymerase (Eurogentec, Maastricht, The Netherlands). The optimized PCR conditions used for the database construction were 94°C for 3 min. followed by 30 cycles of 94°C for 30 s, at the calculated annealing temperature for 30 s, 72°C for 60 s and a final extension at 72°C for 3 min. Unlabeled primers were obtained from Isogen (Maarssen, The Netherlands), fluorescently labelled (HEX, NED, 6-FAM) primers from Applied Biosystems (Warrington, United Kingdom). The amplification products were separated on a 6% acrylamide gel and visualized with silver staining according to Promega Silver sequence DNA sequencing system (Promega, Leiden, The Netherlands) as described [34]. Fluorescent amplification products were combined (see Table 1) and purified using Multiscreen 96-well Sephadex G50 filtration plates (Millipore). One μl of purified sample was mixed with 10 μl of formamide loading buffer containing a ROX-labelled internal lane standard. After denaturation at 95°C for 3 min, followed by quenching on ice, 1 μl samples were loaded in a capillary sequencer (3700 POP6, ABI) and run for 1.5 h. Fragment sizes were determined automatically using Genescan 1.1 (ABI). All genotypes were analyzed using Genotyper 3.5 NT (ABI).

Data analysis

A selection of 12 microsatellite markers with high quality patterns (see Table 1) was used for the characterization of the varieties. Screening of varieties in a first round revealed all existing alleles for each marker and allowed selection of a set of varieties representing all the alleles. These varieties were included in each following run and used as a reference for allele determination. In this way for each marker the alleles were assigned a name (a, b, c, etc.) based on an exact match to the length of the corresponding allele present in the reference variety, rather than as a particular length in base pairs. Only the presence of alleles was scored and recorded as a presence/absence (1/0) matrix. As a consequence, both AAB and ABB genotypes, for example, are scored and entered in the database as AB. We call this the 'allelic phenotype' [7, 8, 24, 35] after Becher et al. [36] to distinguish it from the genotype. An allelic phenotype is not the same as a genotype, as it only includes information on the presence of alleles, not on the allele frequency [26]. We report the number of alleles per locus, the effective number of alleles, and the number of allelic phenotypes. The effective number of alleles (ne) is estimated as 1/Σpi2, where pi is the frequency of the ith allele in the variety examined. We prefer calculating the effective number of alleles to the expected heterozygosity (which is 1-Σpi2). These two measures have a non-linear relationship (nE = 1/(1-Hexp)), and the effective number of alleles scales better when there are many alleles. More importantly, it is less affected by our dominant way of scoring alleles, and has a straightforward interpretation even across ploidy levels. On the basis of individual allele scores Jaccard distances were calculated. The Jaccard distance and the related Dice distance ignore absence-absence pairs, whose number may be inflated by the dominant scoring of a codominant marker. The varieties were clustered using neighbour-joining in NTSYSpc 2.1.

SpaGeDi 1.0b [37], which can handle plants of different ploidy levels, was used to calculate genetic differentiation (Fst) among varieties on the basis of the presence of alleles. The magnitude of the error in allele frequencies caused by scoring only presence/absence and ignoring all presence of more than one copy in diploid and triploid varieties, was estimated for the diploid plants through a comparison with the results of an analysis of codominantly scored data.

For the codominantly scored diploid plants also Nei's heterozygosity, gene diversity, allelic richness, and Fis values were calculated per variety, using SpaGeDi.


  1. Fénart S, Arnaud J-F, De Cauwer I, Cuguen J: Nuclear and cytoplasmic genetic diversity in weed beet and sugar beet accessions compared to wild relatives: New insights into the genetic relationships within the Beta vulgaris complex species. Theor Appl Genet. 2008, 116: 1063-1077. 10.1007/s00122-008-0735-1.

    Article  PubMed  Google Scholar 

  2. Fischer HE: Origin of the 'Weisse Schlesische Rübe' (white Silesian beet) and resynthesis of sugar beet. Euphytica. 1989, 41: 75-80. 10.1007/BF00022414.

    Article  Google Scholar 

  3. Van Geyt JPC, Lange W, Oleo M, De Bock ThSM: Natural variation within the genus Beta and its possible use for breeding sugar beet: A review. Euphytica. 1990, 49: 57-76. 10.1007/BF00024131.

    Article  Google Scholar 

  4. McGrath JM, Saccomani M, Stevanato P, Biancardi E: Beet. Genome Mapping and Molecular Breeding in Plants. Vegetables. Edited by: Kole C. 2007, Berlin: Springer, 5: 191-207. full_text.

    Google Scholar 

  5. Bredemeijer GMM, Cooke RJ, Ganal MW, Peeters R, Isaac P, Noordijk Y, Rendell S, Jackson J, Röder MS, Wendehake K, Dijcks M, Amelaine M, Wickaert R, Bertrand L, Vosman B: Construction and testing of a microsatellite database containing more than 500 tomato varieties. Theor Appl Genet. 2002, 105: 1019-1026. 10.1007/s00122-002-1038-6.

    Article  PubMed  Google Scholar 

  6. Röder MS, Wendehake K, Korzun V, Bredemeijer G, Laborie D, Bertrand L, Isaac P, Rendell S, Jackson J, Cooke RJ, Vosman B, Ganal MW: Construction and analysis of a microsatellite-based database of European wheat varieties. Theor Appl Genet. 2002, 106: 67-73.

    PubMed  Google Scholar 

  7. Esselink D, Smulders MJM, Vosman B: Identification of cut-rose (Rosa hybrida) and rootstock varieties using robust Sequence Tagged Microsatellite markers. Theor Appl Genet. 2003, 106: 277-286.

    CAS  PubMed  Google Scholar 

  8. Smulders MJM, Esselink D, Voorrips RE, Vosman B: Analysis of a database of DNA profiles of 734 Hybrid Tea Rose varieties. Acta Horticult. 2009, 836: 169-174.

    Article  Google Scholar 

  9. Barzen E, Mechelke W, Ritter E, Schulte-Kappert E, Salamini F: An extended map of the sugar beet genome containing RFLP and RAPD loci. Theor Appl Genet. 1995, 90: 189-193. 10.1007/BF00222201.

    Article  CAS  PubMed  Google Scholar 

  10. Hansen M, Kraft T, Christiansson M, Nilsson N-O: Evaluation of AFLP in Beta. Theor Appl Genet. 1999, 98: 845-852. 10.1007/s001220051143.

    Article  CAS  Google Scholar 

  11. Pillen K, Steinrücken G, Wricke G, Herrmann RG, Jung C: A linkage map of Sugar beet (Beta vulgaris L.). Theor Appl Genet. 1992, 84: 129-135. 10.1007/BF00223992.

    Article  CAS  PubMed  Google Scholar 

  12. Pillen K, Sleinrücken G, Herrmann RG, Jung C: An extended linkage map of Sugar beet (Beta vulgaris L.) including nine putative lethal genes and the Restorer gene X. Plant Breed. 1993, 111: 265-272. 10.1111/j.1439-0523.1993.tb00641.x.

    Article  CAS  Google Scholar 

  13. Schondelmaier J, Steinrücken G, Jung G: Integration of AFLP markers into a linkage map of sugar beet (Beta vulgaris L.). Plant Breed. 1996, 115: 231-237. 10.1111/j.1439-0523.1996.tb00909.x.

    Article  CAS  Google Scholar 

  14. Schumacher K, Schondelmaier J, Barzen E, Steinrücken G, Borchardt D, Weber WE, Jung C, Salamini F: Combining different linkage maps in sugar beet (Beta vulgaris L.) to make one map. Plant Breed. 1997, 116: 23-38. 10.1111/j.1439-0523.1997.tb00971.x.

    Article  Google Scholar 

  15. Saccomani M, Stevanato P, Trebbi D, McGrath JM, Biancardi E: Molecular and morpho-physiological characterization of sea, ruderal and cultivated beets. Euphytica. 2009, 169: 19-29. 10.1007/s10681-009-9888-5.

    Article  Google Scholar 

  16. Jones CJ, Edwards KJ, Castiglione S, Winfield MO, Sala F, Wiel Van de C, Bredemeijer G, Vosman B, Matthes M, Daly A, Brettschneider R, Bettini P, Buiatti M, Maestri E, Malcevschi A, Marmiroli N, Aert R, Volckaert G, Rueda J, Linacero R, Vazquez A, Karp A: Reproducibility testing of RAPD, AFLP and SSR markers in plants by a network of European laboratories. Mol Breed. 1997, 3: 381-390. 10.1023/A:1009612517139.

    Article  CAS  Google Scholar 

  17. De Riek J, Calsyn E, Everaert I, Van Bockstaele E, De Loose M: AFLP based alternatives for the assessment of Distinctness, Uniformity and Stability of sugar beet varieties. Theor Appl Genet. 2001, 103: 1254-1265. 10.1007/s001220100710.

    Article  CAS  Google Scholar 

  18. Desplanque B, Boudry P, Broomberg K, Saumitou-Laprade P, Cuguen J, Van Dijk H: Genetic diversity and gene flow between wild, cultivated and weedy forms of Beta vulgaris L. (Chenopodiaceae), assessed by RFLP and microsatellite markers. Theor Appl Genet. 1999, 98: 1194-1201. 10.1007/s001220051184.

    Article  CAS  Google Scholar 

  19. Viard F, Bernard J, Desplanque B: Crop-weed interactions in the Beta vulgaris complex at a local scale: allelic diversity and gene flow within sugar beet fields. Theor Appl Genet. 2002, 104: 688-697. 10.1007/s001220100737.

    Article  CAS  PubMed  Google Scholar 

  20. Rae SJ, Aldam C, Dominguez I, Hoebrechts M, Barnes SR, Edwards KJ: Development and incorporation of microsatellite markers into the linkage map of sugar beet (Beta vulgaris spp.). Theor Appl Genet. 2000, 100: 1240-1248. 10.1007/s001220051430.

    Article  CAS  Google Scholar 

  21. Richards CM, Brownson M, Mitchell SE, Kresovich S, Panella L: Polymorphic microsatellite markers for inferring diversity in wild and domesticated sugar beet (Beta vulgaris). Mol Ecol Notes. 2004, 4: 243-245. 10.1111/j.1471-8286.2004.00630.x.

    Article  CAS  Google Scholar 

  22. Cureton AN, Burns MJ, Ford-Lloyd BV, Newbury HJ: Development of simple sequence repeat (SSR) markers for the assessment of gene flow between sea beet (Beta vulgaris ssp. maritima) populations. Mol Ecol Notes. 2002, 2: 402-403. 10.1046/j.1471-8286.2002.00253.x.

    Article  CAS  Google Scholar 

  23. Laurent V, Devaux P, Thiel T, Viard F, Mielordt S, Touzet P, Quillet MC: Comparative effectiveness of sugar beet microsatellite markers isolated from genomic libraries and GenBank ESTs to map the sugar beet genome. Theor Appl Genet. 2007, 115: 793-805. 10.1007/s00122-007-0609-y.

    Article  CAS  PubMed  Google Scholar 

  24. Brownstein ML, Carpten JD, Smith JR: Modulation of non-templated nucleotide addition by Taq DNA polymerase: primer modifications that facilitate genotyping. BioTechniques. 20: 1004-1010.

  25. Smulders MJM, Bredemeijer G, Rus-Kortekaas W, Arens P, Vosman B: Use of short microsatellites from database sequences to generate polymorphisms among Lycopersicon esculentum cultivars and accessions of other Lycopersicon species. Theor Appl Genet. 1997, 94: 264-272. 10.1007/s001220050409.

    Article  CAS  Google Scholar 

  26. De Riek J, Everaert I, Esselink D, Calsyn E, Smulders MJM, Vosman B: Assignment tests for variety identification compared to genetic similarity-based methods using experimental datasets from different marker systems in sugar beet. Crop Sci. 2007, 47: 1964-1974. 10.2135/cropsci2006.04.0265.

    Article  Google Scholar 

  27. Arnaud J-F, Fénart S, Godé C, Deledicque S, Touzet P, Cuguen J: Fine-scale geographical structure of genetic diversity in inland wild beet populations. Mol Ecol. 2009, 18: 3201-3215. 10.1111/j.1365-294X.2009.04279.x.

    Article  PubMed  Google Scholar 

  28. Arnaud J-F, Fénart S, Cordellier M, Cuguen J: Populations of weedy crop-wild hybrid beets show contrasting variation in mating system and population genetic structure. Evol Applic. 2010, 3: 305-318. 10.1111/j.1752-4571.2010.00121.x.

    Article  Google Scholar 

  29. McGrath JM, Trebbi D, Fenwick A, Panella L, Schulz B, Laurent V, Barnes S, Murray SC: An open-source first-generation molecular genetic map from a Sugarbeet × Table Beet cross and its extension to physical mapping. Crop Sci. 2007, 47: S27-S44.

    Article  Google Scholar 

  30. Viard F, Arnaud J-F, Delescluse M, Cuguen J: Tracing back seed and pollen flow within the crop-wild Beta vulgaris complex: genetic distinctiveness vs. hot spots of hybridization over a regional scale. Mol Ecol. 2004, 13: 1357-1364. 10.1111/j.1365-294X.2004.02150.x.

    Article  CAS  PubMed  Google Scholar 

  31. Arnaud J-F, Viard F, Delescluse M, Cuguen J: Evidence for gene flow via seed dispersal from crop to wild relatives in Beta vulgaris (Chenopodiaceae): consequences for the release of genetically modified crop species with weedy lineages. Proc R Soc Lond B. 2003, 270: 1565-1571. 10.1098/rspb.2003.2407.

    Article  Google Scholar 

  32. Vosman B, Arens P, Rus-Kortekaas W, Smulders MJM: Identification of highly polymorphic DNA regions in tomato. Theor Appl Genet. 1992, 85: 239-244. 10.1007/BF00222865.

    CAS  PubMed  Google Scholar 

  33. Fulton TM, Chunwongse J, Tanksley SD: Microprep protocol for extraction of DNA from tomato and other herbaceous plants. Plant Molecular Biology Reporter. 1995, 13: 207-209. 10.1007/BF02670897.

    Article  CAS  Google Scholar 

  34. Wiel Van de C, Arens P, Vosman B: Microsatellite retrieval in lettuce (Lactuca sativa L.). Genome. 1999, 42: 139-149. 10.1139/gen-42-1-139.

    Article  PubMed  Google Scholar 

  35. Samiei L, Naderi R, Khalighi A, Shahnejat-Bushehri AA, Mozaffarian V, Esselink GD, Kazempour Osaloo S, Smulders MJM: Genetic diversity and genetic similarities between Iranian rose species. J Hortic Sci Biotech. 2010, 85: 231-237.

    CAS  Google Scholar 

  36. Becher SA, Steinmetz K, Weising K, Boury S, Peltier D, Renou J-P, Kahl G, Wolff K: Microsatellites for cultivar identification in Pelargonium. Theor Appl Genet. 2000, 101: 643-651. 10.1007/s001220051526.

    Article  CAS  Google Scholar 

  37. Hardy OJ, Vekemans X: SPAGeDi: a versatile computer program to analyse spatial genetic structure at the individual or population levels. Mol Ecol Notes. 2002, 2: 618-620. 10.1046/j.1471-8286.2002.00305.x.

    Article  Google Scholar 

Download references


This research was funded in part by the Ministry of Agriculture, Nature, and Food quality of the Netherlands under the Method development scheme for Plant Breeders Rights research.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Marinus JM Smulders.

Additional information

Authors' contributions

JdR and BV conceived and designed the study; GDE and IE performed the experiments; MJMS, IE, JdR, and GDE analyzed the data; MJMS, GDE and BV wrote the paper. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1:A table reporting gene diversity, allelic richness and Fis values per marker and variety.(XLS 43 KB)


Additional file 2:A PCO plot of the sugar beet varieties based on pairwise genetic distances between sugar beet varieties calculated using dominant scoring of alleles. Triploid varieties: open circles; diploid varieties: filled circles. (PPT 74 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Smulders, M.J., Esselink, G.D., Everaert, I. et al. Characterisation of sugar beet (Beta vulgaris L. ssp. vulgaris) varieties using microsatellite markers. BMC Genet 11, 41 (2010).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: