Microsatellite isolation and marker development in carrot - genomic distribution, linkage mapping, genetic diversity analysis and marker transferability across Apiaceae
© Cavagnaro et al; licensee BioMed Central Ltd. 2011
Received: 1 April 2011
Accepted: 1 August 2011
Published: 1 August 2011
The Apiaceae family includes several vegetable and spice crop species among which carrot is the most economically important member, with ~21 million tons produced yearly worldwide. Despite its importance, molecular resources in this species are relatively underdeveloped. The availability of informative, polymorphic, and robust PCR-based markers, such as microsatellites (or SSRs), will facilitate genetics and breeding of carrot and other Apiaceae, including integration of linkage maps, tagging of phenotypic traits and assisting positional gene cloning. Thus, with the purpose of isolating carrot microsatellites, two different strategies were used; a hybridization-based library enrichment for SSRs, and bioinformatic mining of SSRs in BAC-end sequence and EST sequence databases. This work reports on the development of 300 carrot SSR markers and their characterization at various levels.
Evaluation of microsatellites isolated from both DNA sources in subsets of 7 carrot F2 mapping populations revealed that SSRs from the hybridization-based method were longer, had more repeat units and were more polymorphic than SSRs isolated by sequence search. Overall, 196 SSRs (65.1%) were polymorphic in at least one mapping population, and the percentage of polymophic SSRs across F2 populations ranged from 17.8 to 24.7. Polymorphic markers in one family were evaluated in the entire F2, allowing the genetic mapping of 55 SSRs (38 codominant) onto the carrot reference map. The SSR loci were distributed throughout all 9 carrot linkage groups (LGs), with 2 to 9 SSRs/LG. In addition, SSR evaluations in carrot-related taxa indicated that a significant fraction of the carrot SSRs transfer successfully across Apiaceae, with heterologous amplification success rate decreasing with the target-species evolutionary distance from carrot. SSR diversity evaluated in a collection of 65 D. carota accessions revealed a high level of polymorphism for these selected loci, with an average of 19 alleles/locus and 0.84 expected heterozygosity.
The addition of 55 SSRs to the carrot map, together with marker characterizations in six other mapping populations, will facilitate future comparative mapping studies and integration of carrot maps. The markers developed herein will be a valuable resource for assisting breeding, genetic, diversity, and genomic studies of carrot and other Apiaceae.
The Apiaceae family (order Apiales) contains a number of important vegetable and spice crop species including carrot, celery, fennel, cilantro, parsley and parsnip. Carrot (Daucus carota L.) is the most economically important species in the Apiaceae, with more than 21 million tons produced yearly worldwide (data for year 2003; http://faostat.fao.org/faostat). Carrot consumption has steadily increased in recent decades for several reasons, including a heightened awareness of its health-promoting attributes (carrots are the richest source of provitamin A carotenoids in the U.S. diet ), development of fresh-cut carrot products convenient for consumers, and adaptation of improved cultivars for warmer production areas .
Despite their economic importance, carrot and other Apiaceae have relatively underdeveloped molecular resources . Three unsaturated linkage maps, mainly based on anonymous dominant AFLP markers, have been constructed in carrot [4, 5] to assist breeding of this species. Inheritance studies on natural carotenoid mutants have identified factors conditioning root pigment accumulation , and both simply inherited pigment traits and QTL have been included in these maps [4, 5, 7]. The more recent addition of 22 genes from the carotenoid biosynthetic pathway , as well as Transposon-Display markers , onto one of the maps has increased both maker informativeness and coverage. The latter map, which was constructed using an F2 family derived from a cultivated × wild carrot cross, is considered the carrot reference map because it harbors important phenotypic traits, has fairly good coverage and includes the largest number of informative markers. However, direct comparisons between the reference and other carrot maps are currently difficult due to the lack of common markers across maps. This fact seriously limits the usefulness of these maps for assisting breeding, especially considering that a number of important carrot traits, including nematode resistance , anthocyanin pigmentation and reduced-sugar accumulation , have been mapped in unrelated genetic backgrounds. The lack of common markers across carrot maps is mainly due to an insufficient availability of informative and robust PCR-based markers. The carotenoid genes mapped by Just et al.  are not easily transferred to other maps with different genetic backgrounds, since they were mainly mapped by sequence detection of single polymorphism nucleotides (SNPs), rather than by fragment length polymorphisms, and because sequence conservation was very high (i.e., lack of polymorphism) in some carotenoid genes . Thus, for these genes, little SNP polymorphism may be expected. The development of polymorphic and robust PCR-based markers in carrot, such as microsatellites, would facilitate their inclusion in different maps, thus serving as anchoring points for map integration. This would immediately increase map marker density, SSR-tagging relevant phenotypic traits and, perhaps, facilitate applications such as positional gene-cloning. In addition, other carrot genetic research studies, such as analysis of genetic diversity and phylogenetic reconstructions, previously approached using anonymous dominant AFLP markers  or laborious time-consuming codominant RFLPs , would also benefit from the development of microsatellite markers.
Microsatellites, or simple sequence repeats (SSRs), are the marker of choice in many molecular genetic applications including mapping, fingerprinting, genetic diversity, population structure analysis, gene flow and germplasm conservation studies. Their widespread adoption is due to several desirable characteristics: they are codominant, frequently and evenly distributed throughout genomes, selectively neutral, highly reproducible and rely on simple polymerase chain reaction (PCR) technology. In addition they are ubiquitous (SSRs can be found in nuclear and mitochondrial genomes of all organisms, as well as in plastid genomes) and hypervariable. The latter property is attributed to a high mutation rate of these repeats resulting from DNA polymerase slippage during DNA synthesis . This mutational mechanism generates gains or losses of one or a few repeat units in the microsatellite, which accumulate more rapidly than point mutations and InDels , leading to a high number of alleles per locus. In plants, the high polymorphism found in microsatellites has allowed the detection of variability in species otherwise characterized by low levels of genetic diversity .
In Apiaceae, very few publicly available SSRs have been reported previously, and these were developed from carrot (9 SSRs ) and celery (11 SSRs ), the two most economically important species in the family. The availability of a large set of SSRs in carrot is likely to benefit research in other Apiaceae, since significant marker transferability has been observed across related taxa . This is of particular interest to research groups working in minor crops or species with limited research funds; many laboratories have sufficient resources and expertise for running SSR-based PCR analyses, although perhaps not for the isolation and characterization of new loci.
In this study we report on the development of 300 new carrot SSR markers. Further characterization of these loci includes analysis of SSR distributions in genomic and EST sequence, linkage mapping onto the carrot reference map, evaluation of their mapping potential in subsets of seven carrot F2 mapping populations, evaluation of their potential for assessing genetic diversity among Daucus carota accessions, and evaluation of SSR marker transferability across 24 Apiaceae taxa.
Distribution of microsatellites in genomic and EST sequence
Distribution of microsatellites GSSRs (BSSRs), and ESSRs) of carrot*
Mean # of repeat units/mean SSR length (bp)
Most frequent motifs‡
Mean # of repeat units/mean SSR length (bp)
Most frequent motifs‡
Mean # of repeat units/mean SSR length (bp)
Most frequent motifs‡
AT(46), AC, AG
AAC, AGT, AAG
AAG(27), AAC(17), AGT(13), ACT(12), AAG(27)
AAG(25), ACT(15), AGG(10), AGC(10)
ACAT (59), AAAC
Total perfect SSRs
Total seq. (Mbp)ζ
GC content (%)
Comparisons among microsatellites isolated from the different sequence datasets (GSSRs, BSSRs, and ESSRs) revealed that di-, tri-, and tetranucleotide repeats accounted for 82.6% to 97.6% of the SSRs among the three types evaluated, with di- and tetranucleotide motifs accounting for most of the GSSRs, while trinucleotide motifs were most common among BSSRs and ESSRs (Table 1).
Within genomic DNA, variation was also found between the two sequence datasets examined. Overall, GSSRs had more repeat units (7.9 versus 4.4) and consequently were longer (23.1 bp versus 13.9 bp) than BSSRs (P < 0.0001). The larger number of repeat units in GSSRs compared to BSSRs was evident and significant (P < 0.01) for di, tri, and tetranucleotides, whereas penta to octanucleotides had the same mean number of repeats in both datasets. SSR length had a similar relationship for these repeat types in both datasets. With regard to repeat types, GSSRs yielded a higher proportion of di and tetranucleotides, as compared to BSSRs, whereas trinucleotides and penta to octanucleotides were more frequent in the latter group. In addition, GSSRs included a significant fraction (16.5%) of often long, compound microsatellites, such repeats being nearly 7 times more frequent in this group compared to BSSRs.
Departure of dinucleotide and trinucleotide genomic (BAC end, BSSR) and EST (ESSR) motifs from expected distributions*
Placement and chi-square significance of microsatellite distribution inside and outside of ORFs from GSSRs, BSSRs, and ESSRs of carrot*
p (df = 1)
p (df = 1)
p (df = 1)
The distribution of SSR types in EST sequences was especially variable, with trinucleotides predominating in ORFs, representing more than 65% of the SSRs found inside ORFs (data not presented) and 26% of all ESSRs, and dinucleotides and tetranucleotides predominating in non-protein coding sequences of the ESTs (i.e., mainly UTRs) (Figure 1C). In EST ORFs, the most frequent trinucleotide motifs were AAG, ACT, AGG, AGC, and ACC. In the dinucleotides-rich UTR region of ESTs, AG and AC motifs were three times more frequent than in protein-coding regions of ESTs (Figure 1F).
SSR marker development
Primer pairs for 156 GSSR and 144 BSSR loci were designed. Of these 300 primers pairs, 243 flanked single SSRs (202 perfect repeats and 41 compound SSRs) and 57 flanked multiple SSRs (i.e., the template sequence flanked by the primers annealing sites included more than one SSR). Markers with single perfect repeats included 59 dinucleotides, 56 trinucleotides, 64 tetranucleotides, 10 pentanucleotides, 7 hexanucleotides, 3 heptanucleotides, and 1 each of mono, octa and nonanucleotides. Further information on these 300 SSR markers is presented in Additional File 1 - Table S1, including primer sequence, annealing temperature, repeat motif and its position in template sequence, expected amplicon length, and the template DNA sequence carrying the SSR (for developing alternative primers if desired). For a number of microsatellites detected computationally, we were not able to design primers because they either lacked suitable flanking sequences or the total sequence length was too short (< 120 bp).
Marker polymorphism analyses in carrot F2 families
Carrot F2 mapping populations used to evaluate SSR markers
Phenotypic segregating traits of interest
B493 × QAL
QTL for carotenoid accumulation
(22 structural carotenoid genes mapped)
HCM × Brasilia
QTL for carotenoid accumulation
Santos and Simon (2002) 
B9304 × YC7262
Y 2 - Differential xylem/phloem carotene levels
Rs - Sugar type (reducing/non-reducing) in roots
P 1 - Purple/yellow pigment accumulation in roots
Vivek and Simon (1999) 
Biennial5 × Criolla INTA
Vern1 - Vernalization requirement
Rf1 - Fertility restoration
[(SNts × Camberly) ×
(Turkish × 2566B)]
Y 2 - Differential xylem/phloem carotene levels
Rs - Sugar type (reducing/non-reducing) in roots
P 1 - Purple/yellow pigment accumulation in roots
P 2 - Purple/green petiole
(PI173687 × B10138) ×
(PI173687 × B493)
Y 2 - Differential xylem/phloem carotene levels
P 1 - Purple/yellow pigment accumulation in roots
(Ping Ding × YC7262) × 8542
Mj 2 - Resistance to Meloidogyne javanica
Substantial variation in the degree of polymorphism was found between the two sets of markers (GSSRs and BSSRs). In general, GSSRs were more polymorphic than BSSRs. Depending on the F2 family, 21.7 - 35% (mean 29.6%) of the GSSRs, and 6.3 - 17.4% (mean 12.8%) of the BSSRs, were polymorphic. Overall, nearly 77% of GSSRs and 52% of BSSRs were polymorphic in at least one F2 family. The mean polymorphism index (PI), which takes in account the type of polymorphism (e.g., dominant, codominant), was significantly higher (P < 0.001) for GSSRs (23.6%) compared to BSSRs (9.8%), regardless of the mapping population.
Markers harboring multiple perfect microsatellites in their amplicon sequence (i.e., "multiple SSRs/amplicon") with 11-15 repeat units were the most polymorphic markers, followed by long dinucleotides and tetranucleotides with more than 11 repeat units. With the exception of trinucleotides, for which there was no clear increase in polymorphism with increased repeat number, the SSR markers considered altogether revealed a clear positive relationship between the two variables, with a nearly 3-fold increase in polymorphism when comparing the lowest (< 6; PI = 9%) and highest (> 15; PI = 26%) repeat number classes.
Marker transferability across Apiaceae
A total of 300 SSR loci were assessed across 23 Apiaceae accessions for a total of 6,900 primer/accession combinations. Of these, 4,346 (63%) produced fragments within the expected size range. Combinations that produced fragments outside the expected size range (i.e., > 100 bp larger or smaller than the original carrot sequence) were considered non-specific amplifications and regarded as negative results. This range was arbitrarily selected to simplify the analysis, especially in the cases where more than two bands were amplified. All successful amplifications were obtained at annealing temperatures between 0 and -2° of the recommended value. Alternative PCR protocols, such as touchdown, did not significantly improve amplification success.
Consistent with the previous analysis, the number of SSR primer pairs that produced expected-size amplicons across all the Daucus carota accessions (8), non-carota Daucus (16) and Apiaceae (23) was 200, 23, and 8, respectively. The performance of each SSR marker across the 23 Apiaceae accessions is presented in Additional File 1 - Table S4.
Summary of the parental (QAL and B493) maps of carrot
Total # markers mapped
# markers common to QAL and B493 map
# markers mapped on the QAL map
# mapped markers segregating only from QAL
# SSR markers mapped
Other mapped markers *
LG size (cM)
Mean marker distance (cM) ǂ
# markers mapped on the B493 map
# mapped markers segregating only from B493
# SSR markers mapped
Other mapped markers *
LG size (cM)
Mean marker distance (cM) ǂ
8A, 1T, 1G, 1S
20A, 2T, 1G, 1S
10A, 3T, 2G
11A, 4T, 2G
13A, 2T, 1G
6A, 2T, 1G
13A, 4T, 5G, 2S
11A, 3T, 5G, 2S
12A, 8T, 5G
11A, 1T, 6G
7A, 2T, 6G
12A, 2T, 2G, 1S
12A, 6T, 2G, 1S
10A, 4T, 1G
12A, 2T, 1G
99A, 27T, 23G, 4S
97A, 23T, 23G, 4S
Overall, after mapping the SSR loci, the linkage map of the wild carrot QAL contains 202 molecular markers (69 codominant and 133 dominant) covering 1,120.8 cM, with an average distance between adjacent markers of 5.8 cM, whereas the cultivated B493 map harbors 193 markers (69 codominant and 124 dominant) covering 1273.2 cM, with a 6.9 cM average marker distance. Thus, although the parental B493 map includes fewer markers, it has a larger total map length than the QAL map. A paired t-test revealed a significantly higher (P = 0.047) mean recombination fraction between adjacent markers (6.9 for B493 and 5.8 for QAL) when comparing the two parental maps. Although marginally significant, the higher mean recombination found in B493 may help explain its larger observed total map length.
Because in a very recent study  the linkage groups (LGs) from this map were integrated with actual chromosomes by means of flourescent in situ hybridization (FISH) mapping of BAC clones anchored by LG-specific markers (including some SSRs), the LGs in Figure 4 were named, ordered and oriented north-south according to the corresponding chromosomes. By convention, chromosomes are numbered consecutively from longest to shortest, and they are oriented with their short and long arms following north-south directions. Thus, correspondences between our LG designations and those from previous maps [5, 8, 9] are as follows: LGs 1, 3, 5, 6, and 7 correspond to former LGs 1, 8, 2, 3, and 5, respectively, with conserved orientation, whereas LGs 2, 4, 8, and 9, in this study (Figure 4) correspond to former LGs 4, 6, 9, and 7, respectively, with inverted north-south orientations.
The 38 fully informative (i.e., codominant) SSR markers consistently mapped to homologous LG pairs of the cultivated (B493) and wild (QAL) carrot genomes (i.e., all SSRs that mapped on the same LG in one parent also did so in the other parent). In addition, for these common markers, complete conservation of locus order was observed along the individual LGs between the two maps (Figure 4). This was also true for other -previously mapped- codominant markers (e.g., carotenoid genes, SCARs). Thus, no evidence of rearrangement of chromosomal blocks between the wild and the domesticated carrot subspecies was found.
Within each LG, SSRs were frequently positionally associated to genes. More than 40% of the SSR markers of both maps (19 and 18 SSRs in the QAL and B493 maps, respectively) mapped within 10 cM from previously mapped genes, whereas 16 of the 22 (~73%) carotenoid genes - in both, QAL and B493 maps- had 1 or more SSRs within a 10 cM distance, suggesting that these repeats are relatively frequent in genic regions of the carrot genome.
SSR diversity in Daucus carota accessions
General diversity statistics for 10 SSR loci evaluated in 65 Daucus carota accessions
Number of alleles
Allele size range (bp)
Frequency and distribution of SSRs in carrot genomic and EST sequence
Microsatellite density in genomic DNA of carrot was estimated by analysis of 1.74 Mbp of BAC end sequence (the GSSRs dataset was excluded from this analysis because it derived from an SSR-enriched library and, therefore, its analysis would result in an overestimation of the SSR density in genomic sequence). Carrot had a rather low SSR density (134.5 SSRs/Mbp) compared to other species. SSRs analyses -using the same search parameters and program as with carrot- in the complete genome sequences of four model species revealed SSR densities of 370, 507, 529, and 508 SSRs/Mbp in Arabidopsis thaliana, grapevine, rice, and poplar, respectively. The lower SSR density in carrot compared to these species cannot be attributed to differences in the source of genomic sequence (BAC ends versus whole genomes) since analyses of BAC end sequence (BES) datasets from these and other species were also much more dense in microsatellites than carrot BES (data not presented). Similarly, transcript sequences of carrot (214.8 SSRs/Mbp), although more dense in SSRs than their genomic counterparts, were also less frequent in these repeats compared to ESTs of Arabidopsis (358 SSRs/Mbp), grapevine (247 SSRs/Mbp), poplar (425 SSRs/Mbp), soybean (403 SSRs/Mbp), rice (739 SSRs/Mbp), and sorghum (646 SSRs/Mbp).
Carrot trinucleotides were more frequent in transcripts than in genomic DNA. In addition, within BSSRs trinucleotide repeats occurred preferentially inside ORFs, and accounted for ~ 50% of the total SSRs found in these protein coding regions. The abundance of these repeats in ESTs and in ORFs is consistent with the notion that protein-coding sequences tolerate better frame-shift mutations (InDels) of 3 bp -or multiples of 3 bp- than other InDel lengths. Thus, trinucleotide repeats within coding sequences may translate fully functional proteins with a few extra (or fewer) aminoacids, whereas InDels of other lengths would translate abnormal, often deleterious, proteins. Consistent with our results, an overrepresentation of trinucleotides in protein-coding sequences has been reported previously in numerous plant species [20–24], as well as in other eukaryotes including humans, primates, rodents and insects [25, 26]. The relative abundance of trinucleotides over other SSR types has been attributed not only to negative selection against frame-shift mutations in the coding regions but also to positive selection for specific single amino-acid stretches .
DNA polymerase slippage is the main mutational mechanism leading to changes in microsatellite length . These changes in SSR size are most often gradual and step-wise since polymerase slippage only generates gains or losses of one or a few repeat unit(s) . Thus, the fact that SSRs in carrot transcripts generally had fewer repeat units than SSRs in genomic sequence, even for trinucleotide repeats (trinucleotides were twice as frequent in ESTs compared to genomic data), suggests a negative selection pressure against microsatellite size increase in protein-coding sequences.
The non-random distributions of motif sequences among dinucleotide and trinucleotide SSRs of carrot included a higher than expected incidence of (AT)n repeats in genomic DNA (BAC ends), like that of several plant species including soybean, Arabidopsis and rice , but unlike the (AC)n predominant motif among dinucleotides in humans . In contrast, the (AT)n motif was less often observed in ESTs than expected, while (AG)n and (CT)n were more common than expected. This may suggest different constraints for repeat motifs across diverse organisms.
Marker development and analyses in F2 families
In this study, two different strategies were used for isolating and developing carrot SSR markers. The hybridization-based approach, as described by Glenn and Schable , yielded microsatellites (GSSRs) that were, in average, significantly longer (23.1 bp versus 13.9 bp) and had more repeat units (7.9 versus 4.4) than SSRs from BAC end sequences (BSSRs). These differences are, most likely, due to differences in the two strategies used. DNA library enrichment methods based on hybridization capture [28–31] are generally designed to yield a higher proportion of SSRs with large number of repeat units, targeting mainly long perfect repeats. Under this system, long DNA stretches of perfect repeats are hybridized more efficiently to the microsatellite probes and they are retained at a higher rate, compared to short repeats, during the washing steps, thus, increasing the relative proportion of long microsatellite sequences in cloned colonies . Conversely, the BSSRs set represents a random sample -without enrichment for length, repeat type or sequence motif- from genomic DNA. Because of this, they provide a more reliable picture of the microsatellite distribution in the carrot genome. Longer and more repetitive SSRs have been obtained through hybridization-based methods compared to sequence-searches in other plant species, regardless of the type of DNA examined (i.e., genomic or ESTs), including Brassica[32, 33], cotton , wheat and rice .
The differences in repeat number and length between GSSRs and BSSRs have important implications for marker potentiality, particularly with regard to polymorphism. In general, GSSRs were significantly more polymorphic than BSSRs, considering both the polymorphism index (PI) (23.6% versus 9.8%) and the percentage of polymorphic markers (77% versus 52%), and these differences were associated to a higher repeat number and length in the GSSRs group (as suggested by the significant positive correlations obtained between both variables and PI). Developments of SSR markers from other plant species, including cotton , barley  and pine , have also noted positive relationships between SSR polymorphism and number of repeat units. Together, these results are consistent with studies reporting that both SSR polymorphism and SSR mutation rate have a positive relationship with repeat number [38–40]. Concordantly, positive and significant relationships have also been found between repeat length and mutation rate in human , fruit fly  and yeast  microsatellites. These studies indicate that polymerase slippage, the main mutational mechanism in microsatellites , increases with higher repeat number and length, leading to a higher diversity in longer, more repetitive SSRs, as observed in the present study. However, contrary to these and our results, studies using markers developed from other plants, such as Brassica and pearl millet , have reported lack of correlation between size of the SSR, both measured by length (bp) and repeat number, and detection of polymorphic loci. As pointed out in the latter two studies, SSR evolutionary age is a key factor for SSR diversity (i.e., recently evolved microsatellites would have fewer polymorphisms because of fewer occasions for mutation, even if they are relatively long) and this may help explain the lack of association found by them. In addition, most of the above studies (including ours) cannot rule out the possibility that InDels at regions other than the SSR motifs may account for some of the polymorphisms, thus influencing the expected relationship between length and polymorphism.
A major interest for evaluating the SSR markers in the carrot F2 populations was to assess their potential for mapping. Linkage maps using some of these F2s have already been constructed (see Table 4) and others are underway (Simon, personal communication). These maps include different phenotypic traits of interest (Table 4) and -before this study- they were mainly constructed using anonymous dominant markers, such as AFLPs and RAPDs, with only very few markers, or none, in common, thus, making their comparative analyses and/or integration difficult. The present work identified 123 SSRs (87 GSSRs and 36 BSSRs) that were polymorphic in two or more mapping populations, suggesting that these common markers may serve as anchoring points for merging carrot maps. Besides the inclusion of 56 SSR markers onto the carrot reference map (see below), work is underway in our lab to include these polymorphic SSRs in other maps with different genetic backgrounds (see Table 4). The integration of carrot linkage maps would enhance their usefulness for assisting breeding of this species, by increasing marker saturation nearby genes of interest and thereby facilitating applications like positional gene cloning, among others.
From our evaluation in seven carrot F2 families, 196 SSR markers (65%) were polymorphic in at least one mapping population. Because the PCR amplicons were size-separated using high-resolution agarose gel electrophoresis, which can only resolve fragments with size differences of at least 3 bp, a fraction of the markers evaluated in some populations, generated ambiguous band patterns. Although they may have been polymorphic, the bands were too close together in the gel to unambiguously score, and were classified as monomorphic (i.e., only unambiguously polymorphic and scorable markers were classified as "dominant" or "codominant" in Additional File 1 - Table S2). Thus, if other fragment separation systems, with better resolution, are used, such as separation of fluorescently-labeled fragments through capillary electrophoresis, the number of polymorphic markers may be expanded in some populations, particularly in cases of dinucleotide SSR markers varying in a single repeat unit.
High PCR amplification efficiencies were found in the F2 families for both sets of markers, GSSRs (83%) and BSSRs (87%). Comparable amplification efficiencies have been found in other plant species with SSR markers developed using hybridizations-based methods (~ 90% ) and sequence-based searches (85% ).
Transfer success of carrot SSRs across Apiaceae
The availability of SSR loci for economically important species has increased interest in primer transferability to related taxa, especially for species in which molecular resources are limited. In Apiaceae, only a few publicly available SSRs have been reported previously, and these were developed from carrot (9 SSRs ) and celery (11 SSRs ), the two most economically important crop species in the family. Results from this study indicate that a significant fraction of carrot SSRs transfer successfully across Apiaceae. Locus amplification success was detected in 91 to 224 markers across 15 non-carrot Apiaceae species, including economically important crops like parsley (131 SSRs), celery (133 SSRs) and cilantro (91 SSRs). Prospects of a broader utilization of these markers beyond carrot include their application in taxonomic, population, and conservation studies as well as for mapping and assisting breeding in crop species.
It is, however, important to bear in mind that when using SSR markers across distantly-related species the amplification of a PCR product does not necessarily imply locus conservation, since size homoplasy, i.e. convergence in size of non-homologous fragments, may occur. Considering the possibility of this source of confusion, verification of the PCR product identity by sequencing has been suggested previously, particularly when working across genera and if there is uncertainty regarding the size range of the amplicons obtained . However, verification through sequencing may not be necessary if working within the same genus as the species from which the SSRs markers were developed . Thus, the use of carrot SSR markers for studies in non-Daucus Apiaceae should include verification, by sequencing, of the homology to the carrot SSR product sequence (see Additional File 1 - Table S1).
Transfer of carrot SSRs across Daucus species (carrot accessions excluded) was, in general, less successful than SSR transfer rate at the subgenus level reported for other species, whereas transfer of carrot SSRs across-genera was relatively higher than found in other plants. According to a previous review of SSR cross-transferability in plants , the average transferability across species in the same genus was 76.4%, and across related genera was 35.2%. We found these values to be 58.3% across Daucus species and 41% across the Apiaceae. However, it should be noted that SSR transfer success varied greatly across the different reports for both within-same genus (4.7 - 100%) and across different genera (0 - 71.4%) . The huge variation found across these studies likely reflects differences in phylogenetic distance (and thus, also in conservation of sequences at priming sites) between the source and target taxa within each family, as well as differences in the number of taxa and SSR loci analyzed, and in the type of sequences used for marker development. For example, EST-derived SSRs are more conserved and thus they transfer across genera more readily than genomic SSRs ), among other factors.
Our data (Figure 3) suggest generally a higher rate of success in amplifying carrot SSRs in plants more closely related to carrot. This should not be surprising since closer-related taxa have higher overall sequence homology which translates to more conserved SSR flanking regions and, therefore, easier transferability of primer pairs. Negative relationships between SSR transfer success and phylogenetic distance between source and target taxa have been widely observed in many plant families [18, 46].
The potential usefulness of SSR markers for diversity and phylogenetic studies in Apiaceae will depend, to a great extent, on the possibility that markers successfully amplify across different species and on the ability of the marker to detect polymorphism among the taxa. To have a preliminary picture of how suitable the SSR markers developed in this work may be for these applications, we investigated interspecific SSR variation among non-carrot species by analysis of amplicons sizes in the agarose gel images. Thus, for each SSR, the total number of different alleles in the non-carrot species dataset was recorded (Additional File 1: Table S4). Only SSRs that successfully amplified products in at least 80% of the non-carrot species (i.e., SSRs that generated amplicons in at least 12 of the 15 non-carrot species used to assess marker transferability) were considered. Overall, our results revealed 88 SSRs that generated amplicons in most (> 80%) outside-carrot species. Of these, 40 SSRs (29%) produced 3-9 different alleles (with an average of 4.9 alleles/SSR) in the non-carrot group. It should be noted that our calculation of 4.9 alleles/SSR in these selected markers is conservative, due to the low resolution of agarose gels which do not allow discrimination of different alleles varying in one or a few repeats. These results suggest that a significant proportion of the SSR markers developed herein may be suitable for addressing taxonomic or phylogenetic questions within Apiaceae.
Further analysis of the 88 SSRs that produced amplicons in the majority of the non-carrot taxa revealed interesting differences between the two SSR datasets. Although more BSSRs than GSSRs (52 and 36 markers, respectively) amplified successfully in most non-carrot taxa, GSSRs were much more polymorphic than BSSRs at the interspecific level. For example, among GSSRs 28 markers produced 3 or more different alleles (with a range of 3-9 and mean of 5.5), whereas only 12 BSSRs generated 3 or more alleles/SSR (with a range of 3-6 and mean of 3.6). It is likely that the generally higher polymorphism of GSSRs compared to BSSRs at the inter-specific level, which is in agreement with our results for both sets of markers in the carrot F2s, may be also due to the higher number of repeat units present in GSSRs.
SSR linkage mapping
Prior to this work, important advances were made in the construction of carrot genetic maps with a range of molecular marker systems. Although some RFLPs  and a few SCAR and gene-specific markers were mapped , the most extensive genetic mapping data in carrot has been generated mainly with dominant AFLP, RAPD and Transposon-display (TD) markers [4, 5, 8–10]. While RFLPs are useful for comparative mapping purposes, high throughput genotyping and probe handling are difficult. Similarly, the carotenoid genes mapped by Just et al.  are not as easily transferred to other mapping backgrounds since their analysis relied in most cases on SNPs, due to the lack of larger polymorphisms (e.g., InDels) in these genes that can be scored as easily as SSRs. On the other hand, AFLP, RAPD and TD markers, while providing a relatively large number of markers per assay and good genome coverage, have limited information content and are not of much use for comparative mapping purposes and for validating QTL across pedigrees [8–10]. The addition of 55 SSR markers to the carrot reference linkage map together with detailed characterization of this novel set of 300 SSRs in subsets of six other mapping populations should allow significant advances in carrot comparative mapping and map-integration. The fact that most of the mapped SSRs were codominant (38 SSRs) in the B493 × QAL-derived population, with 2-8 informative markers per linkage group, together with the identification of putative codominant SSRs in other mapping populations adds extra value to the data published here for pursuing these goals. The inclusion of SSRs in linkage maps with additional pedigrees is currently underway.
The parental B493 map has a slightly larger total map length than the QAL map. Although the higher mean recombination found in B493 may help explain its larger map length, other factors -e.g., related to the type of markers used- may also cause this effect. Different recombination frequencies can be obtained just by sampling of the different markers, as well as errors derived from calculations of genetic distances from dominant markers data.
In the current map, we have modified linkage group designations and orientations, in accordance to recent cytogenetic data concerning the integration of carrot LGs with actual chromosomes . Following standard conventions, consecutive numbers were assigned to the LGs in decreasing order of chromosome length (i.e., LG1 corresponds to the longest chromosome), and four LGs were inverted in their north-south orientations to agree with the standard short arm/long arm presentation of their corresponding chromosomes. It must, however, be noted that although all the LGs could be unequivocally associated to chromosomes, and thus their number designations are correct and complete, unambiguous LG orientations could only be defined for six of the nine LGs. Thus, chromosomes 4, 6, and 9 in the current map, which correspond to former LGs 6, 3, and 7, respectively, could not be unequivocally oriented, because a single anchored BAC probe was used for LG-chromosome integrations. Thus, their orientations were not modified from previous map versions [5, 8, 9]. However, the possibility remains that future cytogenetic data (for example by FISH analysis with several BAC probes SSR-anchored to these LGs) may reveal different orientations for these LGs. These modifications based on recent cytogenetic data, and the addition of 55 new SSR markers, add value to the updated reference carrot linkage map presented herein. Overall, the current maps involve 193-202 mapped loci, including 69 highly informative markers which consist of SSR, carotenoid gene and SCAR markers, spanning 1,121-1,273 cM, making it the most comprehensive genetic linkage map in the Apiaceae to date.
The SSR loci mapped across all 9 LGs in both parental maps, and they were distributed fairly evenly within most individual LGs, thus recommending their usefulness as anchor points for merging carrot maps. In addition, such dispersed map distribution of the SSR loci, has allowed us to develop BAC FISH probes carrying SSR sequences mapped to specific LGs. These SSR-anchored probes were used for integrating some LGs of carrot with chromosomes by FISH mapping .
The positional association observed between SSRs and previously mapped genes suggests that these tandem repeats are frequent in genic regions of the genome. This is in agreement with results of Morgante et al. , demonstrating higher microsatellite frequencies in the transcribed and non-repetitive fractions of plant genomes. One SSR (gssr112) and two SSRs (gssr12 and gssr119) in LG7 and LG5, respectively, were located in the vicinity of two highly-significant quantitative trait loci (QTL) for total root carotene accumulation. These correspond to the Y and Y 2 loci, respectively, described by Buishand and Gabelman . Microsatellites gssr12 and gssr119, although not tightly linked to Y 2 , may be useful for marker assisted selection, either as a complement of, or as an alternative to the lack of amplification of other robust more-closely linked markers, such as Y2mark.
This work reports on the development and characterization at various levels of a novel set of 300 carrot SSR markers. Analysis of the distribution of SSR motifs, repeat lengths and polymorphism across genomic and EST sequences, as well as concerning different SSR isolation methods (hybridization-based versus sequence SSR-mining) may help decide over future strategies for developing valuable SSR markers in this and other species. The genetic mapping of 55 SSR loci onto the reference carrot linkage map, distributed throughout all 9 linkage groups, together with the characterization of the entire set of markers in 6 other mapping populations, should facilitate comparative mapping studies and integration of carrot maps. Particularly important for these purposes are the 38 codominant SSRs added to the map, which resulted in more than doubling the original number of informative markers in this -or any other- carrot map reported to date. In addition, SSR evaluations in carrot-related taxa indicates that a significant fraction of the carrot SSRs transfer successfully across Apiaceae, with heterologous amplification success rate decreasing with the target-species evolutionary distance from carrot. Nonetheless, a fairly large number of potentially useful SSR markers were identified for non-carrot Daucus (128-224 SSRs) and non-Daucus Apiaceae (91-134 SSRs) species, increasing the prospects of their successful utilization in other Apiaceae. In addition, allelic diversity at selected SSR loci was evaluated using 65 D. carota accessions. In this germplasm, the microsatellites proved to be highly polymorphic, with an average of 19 alleles/locus and 0.84 expected heterozygosity. The marker resources developed in this work should be a valuable tool for carrot breeding and genetics.
SSR identification and marker development
Two different approaches were used to isolate carrot genomic SSRs: 1) Construction and sequence analysis of a carrot (inbred line B493) genomic DNA library enriched for SSR loci (GSSRs) and 2) Bioinformatic mining for SSR motifs in a 1.7 Mbp BAC-end sequence (BES) database (BSSRs). GSSRs were developed at the Savannah River Ecology Laboratory, University of Georgia, using a hybridization capture approach for genomic library enrichment, as described by Glenn and Schable . The DNA clones were sequenced from both directions using standard Sanger cycle-sequencing, and SSRs were detected using the program MISA . The same software was used for the identification of BSSRs in 2,696 carrot BES (NCBI acc. # FJ147695-FJ150390) derived from inbred line B8503 . Only SSRs with repeat length ≥12 nt and 3 or more repeat units were considered. Primer pairs flanking 156 GSSRs and 144 BSSRs were designed with Primer 3 (v.0.4.0).
For comparison purposes only (no markers were developed) a 3.82 Mbp EST dataset generated from 10-week old carrot (inbred B493) roots, was mined for microsatellites using the same programs and parameters described above. The EST dataset comprised 7,285 unique transcripts, 4,044 contigs and 3,241 singlets, which derived from initial analyses and processing (i.e., cleaning, trimming and assembly) of 18,044 Sanger sequence reads. SSRs in the EST dataset with a repeat length ≥12 nt and 3 or more repeat units were included. The resulting data were compared with microsatellites found in genomic DNA sequence (i.e., GSSRs, BSSRs).
For comparisons with carrot, the complete genome sequence of Arabidopsis thaliana L. (119.2 Mbp; accs. NC003070-71, NC003074-76), rice (Oryza sativa L.; 370.8 Mbp; accs. NC008394-405), grapevine (Vitis vinifera L.; 303.1 Mbp; accs. NC012007-25), and poplar (Populus trichocarpa L.; 307 Mbp; accs. NC008467-85) were downloaded from the National Center for Biotechnology (NCBI) database (Genomes Section), and mined for SSRs using the same search parameters and software. As a source of transcript sequences we used plant gene indices of The Institute of Genome Research (TIGR), which are non-redundant (unigenes) EST collections . Thus, gene indices of Arabidopsis (74.8 Mbp; AGI.release_13), poplar (67.6 Mbp; PPLGI.release_4), grapevine (81.4 Mbp; VVGI.release_6), Medicago truncatula (51.9 Mbp; MTGI.release_9), soybean (51.3 Mbp; GMGI.release_13), rice (158.2 Mbp; OGI.release_17) and sorghum (32.4 Mbp; SBGI.release_9) were downloaded from the Gene Index databases (http://compbio.dfci.harvard.edu/tgi/) and searched for SSRs.
A custom Perl program performed a computational ORF detection, which was a search in each possible reading frame for an ATG start codon followed by a stop codon (TAA, TAG, or TGA) at a distance of 100 nt or greater, with no intervening start or stop codons in that reading frame. If two ORFs in different reading frames overlapped, the longer ORF was selected and the shorter ORF was disregarded. Details on the program used for finding ORFs are included in Additional File 2. SSRs, detected by MISA, were categorized as being either inside or outside ORFs. SSRs bridging an in-frame to out-of-frame boundary were discarded from further analysis.
The t-test statistic was used to compare SSR frequencies among the datasets using the program STATGRAPHICS Centurion XV. Evaluation of differences in repeat number across the GSSR, BSSR, and ESSR datasets included chi-square goodness of fit tests to compare observed SSR distributions within each dataset with regard to 1) distribution across sequence motif, 2) distribution of repeat motif and 3) distribution of SSR as inside or outside identified ORFs, using only those clearly categorized. Posterior probability distributions for GSSRs, BSSRs, and ESSRs were calculated separately. Sequence motif distribution posterior probability was calculated from the overall base composition of each dataset. Posterior probability distributions for SSR placement inside or outside ORF regions was based on the ORF sequence distribution within each dataset.
PCR conditions and electrophoresis
PCR reactions were performed in 15 μl volume containing 7.15 μl water, 1.5 μl 10 × DNA polymerase buffer, 1.2 μl dNTPs (2.5 mM each), 1 μl of each primer at 5 μM, 0.15 μl Taq Polymerase at 10 u/μl (Promega, Madison, Wisconsin, USA) and 3 μl of genomic DNA. Thermocyclers were programmed as follows: initial denaturation at 94°C for 3 min., followed by 40 cycles of 94°C for 20 sec., appropriate annealing temperature for 1.0 min., and 72°C for 1.0 min., and a final step at 72°C for 5.0 min. Electrophoresis was carried out for 4-5 hours at 200 V on 4.5% high-resolution agarose (GenePure Hi-Res Agarose, ISC Bioexpress, Kaysville, UT) TAE gels supplemented with 4 ul (5.0 mg/ml) of ethidium bromide for each 100 ml of TAE.
Different methods for marker generation and analyses, including primer labeling, PCR conditions and separation of amplicons, were used for the genetic diversity analyses (described in detail in section "diversity analysis").
Marker analyses in carrot F2 families
Since all carrot linkage maps reported to date were constructed using predominantly anonymous dominant markers such as AFLPs and this has severely limited map merging  microsatellite markers were developed to serve as anchor points across carrot maps. All SSR primer pairs (156 GSSRs, 144 BSSRs) were evaluated in samples (16 DNAs) from 7 carrot F2 mapping populations, as well as in the parental DNAs (when available). Information regarding the populations is presented in Table 4. Markers were evaluated based on their PCR amplification efficiency (i.e., whether an amplicon of expected size was generated) and polymorphism. For the latter, a polymorphism index (PI) was developed according to the formula: PI = [(2C + D)/(7 - nd) × 2] × 100, where 'C' is the number of populations for which the markers was codominant, 'D' is the number of populations for which the marker was dominant, 'nd' populations without information on the performance of the marker.
Simple regression analyses were performed between PI and other characteristics of the SSR markers (e.g., number of repeat units, repeat length, etc.) to investigate possible microsatellite features associated with polymorphism. For this purpose, the program STATGRAPHICS Centurion XV was used.
Marker transferability across Apiaceae
To evaluate the potential utilization of SSRs within Daucus carota as well as in other carrot-related taxa, the GSSRs and BSSRs were tested in a sample of 23 Apiaceae accessions including 8 accessions of carrot (D. carota) (1 inbred line, 1 wild carrot, and 6 commercial cultivars), 8 accessions of non-carrot Daucus species, and 7 accessions of non-Daucus Apiaceae species. Accessions that produced an amplicon of the approximate length expected for that particular SSR in carrot were considered as successful PCR amplifications. PCR products were generated and resolved as described previously.
Plant materials and DNA extraction
A genetic linkage map was constructed using a subset of 103 individuals from a previous F2 mapping population derived from the cross between the white-root wild carrot Queen Annes Lace (QAL) and the cultivated orange-root carrot inbred line B493. Details concerning the development of the mapping population, plant cultivation, DNA extraction, and detection/scoring of previously mapped AFLP, SCAR, DcM TD (DcMaster Transposon Display), and gene specific markers were described before [5, 8, 9, 52].
Generation and analysis of marker data
PCR reactions were performed in the same way as the SSR reactions described above.
For purpose of marker genotyping, the F2 DNAs were analyzed in parallel with controls DNAs of QAL and B493. The genotypes of polymorphic SSRs were recorded as follows: homozygous maternal (B493) "A", homozygous paternal (QAL) "B", heterozygous "H", not A "C", not B "D", and missing data "-". The degree of segregation distortion associated with newly identified SSRs was determined by marker data comparison against the expected ratio of 1:2:1 (A:H:B) for codominant, and 3:1 (C:A or D:B) for dominant markers for an F2 using chi-square tests, where significant distortion was declared at P < 0.01 . Separate maps were constructed for each parent to avoid problems related to the use of repulsion phase dominant markers, as described previously . Dominant markers from a single parent linked in coupling were used in conjunction with all codominant markers. Linkage maps were constructed with MapMaker/EXP 3.0 , where markers were associated with the 'group' command at LOD = 4.0 and a maximum recombination frequency of 0.30. Markers within a group were ordered using 'three point' analysis followed by the 'order' command. Remaining markers were located using the 'try' command, and the map order was re-tested using the 'ripple' command. Recombination frequencies were converted to centimorgans (cM) using the Kosambi function.
SSR genetic diversity in Daucus carota
Plant materials and DNA extraction
A total of 65 Daucus carota accessions were used in this study, including 50 cultivated carrots (D.c. spp. sativa) and 15 wild relatives of carrot (12 accs. of D.c. spp. carota, and one accession each of D.c. spp. azoricus, D.c. spp. hispanicus, D.c. spp. drepanensis) (Additional File 1 - Table S5). The cultivated carrot accessions included 8 inbred lines and 42 open-pollinated cultivars with representatives from 14 of the 16 European primary cultivars , and the predominant cultivars from North America ('Imperator') and Asia ('Kuroda'). For evaluation of morphological traits the carrot accessions were grown in commercial carrot fields in Wisconsin and California. For DNA extractions, plants were grown in greenhouses in Madison, Wisconsin, and genomic DNA was extracted from single plants as described by Murray and Thompson .
PCR amplification and fragment analysis
Ten unlinked microsatellite loci (gssr3, gssr4, gssr6, gssr9, gssr16, gssr35, gssr65, gssr107, gssr85, and gssr111) were used to investigate genetic diversity among cultivated and wild Daucus accessions. Each amplification reaction was performed in a 10 μl final volume, and included ~ 50 ng of DNA template, 0.65 U of Taq DNA polymerase (EconoTaq™), 1 × PCR buffer with 1.5 mM MgCl2 (EconoTaq™), 0.15 μM of reverse primer, 0.15 μM of M13-tailed forward primer (a universal M13 tail, labelled either with 6-FAM, HEX or NED fluorochromes, was added to the 5'-end of the forward primer to enable fluorescent labelling of the amplicons as described by Shuelke ), 100 μM each of the deoxynucleotidetriphosphates (dNTPs), 0.1 mg/ml bovine serum albumin (BSA), and 1% (w/v) of polyvinylpyrrolidone (PVP). The amplification conditions were 94°C for 4 min; 40 cycles of 94°C for 20 sec, 55-60°C for 1 min, 72°C for 1 min; and a final extension step of 72°C for 5 min. Estimation of amplicon lengths and microsatellite genotyping was performed at the University of Wisconsin Biotechnology Center using an ABI 3730xl capillary sequencer and GeneMarker software version 1.5 (SoftGenetics, State College, Pennsylvania).
General statistics, such as number of alleles (NA), allele length, and expected heterozygosity (H e ), were calculated for these markers using the software Arlequin ver 3.1 .
amplified fragment length polymorphism
restriction fragment length polymorphism
expressed sequence tag
million base pairs
polymerase chain reaction
single nucleotide polymorphism
random amplified polymorphic DNA
simple sequence repeats
quantitative trait loci
fluorescent in situ hybridization
bacterial artificial chromosomes
sequence characterized amplified regions
SSRs from BAC end sequence
SSRs from EST libraries
SSRs from genomic library enrichement.
The authors acknowledge Anne E. Atkins for valuable assistance on experiments concerning marker evaluations. This work was supported by the Agricultural Research Service, United States Department of Agriculture.
- Simon PW, Pollak LM, Clevidence BA, Holden JM, Haytowitz DB: Plant breeding for human nutrition. Plant Breed Rev. 2008, 31: 325-392.Google Scholar
- Rubatzky VE, Quiros CF, Simon PW: Carrots and related Umbelliferae. 1999, New York: CABI PublishingGoogle Scholar
- Bhasi A, Senalik D, Simon PW, Kumar B, Manikandan V, Philip P, Senapathy P: RoBuST: An integrated resource of genomics information for plants in the root and bulb crop families Apiaceae and Alliaceae. BMC Plant Biol. 2010, 10: 161-10.1186/1471-2229-10-161.View ArticlePubMedPubMed CentralGoogle Scholar
- Vivek BS, Simon PW: Linkage relationships among molecular markers and storage root traits of carrot (Daucus carota L. ssp. sativus). Theor Appl Genet. 1999, 99: 58-64. 10.1007/s001220051208.View ArticleGoogle Scholar
- Santos CAF, Simon PW: QTL analyses reveal clustered loci for accumulation of major provitamin A carotenes and lycopene in carrot roots. Mol Genet Genomics. 2002, 268: 122-129. 10.1007/s00438-002-0735-9.View ArticlePubMedGoogle Scholar
- Simon PW: Domestication, historical development, and modern breeding of carrot. Plant Breed Rev. 2000, 19: 157-190.Google Scholar
- Just BJ, Santos CA, Yandell BS, Simon PW: Major QTL for carrot color are positionally associated with carotenoid biosynthetic genes and interact epistatically in a domesticated × wild carrot cross. Theor Appl Genet. 2009, 119: 1155-1169. 10.1007/s00122-009-1117-z.View ArticlePubMedGoogle Scholar
- Just BJ, Santos CAF, Fonseca MEN, Boiteux LS, Oloizia BB, Simon PW: Carotenoid biosynthesis structural genes in carrot (Daucus carota): isolation, sequence-characterization, single nucleotide polymorphism (SNP) markers and genome mapping. Theor Appl Genet. 2007, 114: 693-704. 10.1007/s00122-006-0469-x.View ArticlePubMedGoogle Scholar
- Grzebelus D, Jagosz B, Simon PW: The DcMaster transposon display maps polymorphic insertion sites in the carrot (Daucus carota L.) genome. Gene. 2007, 390: 67-64. 10.1016/j.gene.2006.07.041.View ArticlePubMedGoogle Scholar
- Boiteux LS, Belter JG, Roberts PA, Simon PW: RAPD linkage map of the genomic region encompassing the rootknot nematode (Meloidogyne javanica) resistance locus in carrot. Theor Appl Genet. 2000, 100: 439-446. 10.1007/s001220050057.View ArticleGoogle Scholar
- Bradeen JM, Bach IC, Briard M, Le Clerc V, Grzebelus D, Senalik DA, Simon PW: Molecular diversity analysis of cultivated carrot (Daucus carota L.) and wild Daucus populations reveals a genetically nonstructured composition. J Amer Soc Hort Sci. 2002, 127: 383-391.Google Scholar
- Vivek BS, Simon PW: Phylogeny and relationships in Daucus based on restriction fragment length polymorphisms (RFLPs) of the chloroplast and mitochondrial genomes. Euphytica. 1999, 105: 183-189. 10.1023/A:1003446301145.View ArticleGoogle Scholar
- Schlötterer C, Tautz D: Slippage synthesis of simple sequence DNA. Nucleic Acids Res. 1992, 20: 211-215. 10.1093/nar/20.2.211.View ArticlePubMedPubMed CentralGoogle Scholar
- Weber JL, Wong C: Mutation of human short tandem repeats. Human Mol Genetics. 1993, 2: 1123-1128. 10.1093/hmg/2.8.1123.View ArticleGoogle Scholar
- Peakall R, Gilmore S, Keys W, Morgante M, Rafalski A: Cross species amplification of soybean (Glycine max) simple sequence repeats (SSRs) within the genus and other legume genera: implications for the transferability of SSRs in plants. Mol Biol Evolution. 1998, 15: 1275-1287.View ArticleGoogle Scholar
- Niemann M, Westphal L, Wricke G: Analysis of microsatellite markers in carrot (Daucus carota L. sativus). J Appl Genet. 1997, 38A: 20-27.Google Scholar
- Acquadro A, Magurno F, Portis E, Lanteri S: dbEST-derived microsatellite markers in celery (Apium graveolens L. var. dulce). Mol Ecology Notes. 2006, 6: 1080-1082. 10.1111/j.1471-8286.2006.01440.x.View ArticleGoogle Scholar
- Rossetto M: Sourcing SSR markers from related plant species. Plant genotyping: the DNA fingerprinting of plants. Edited by: Henry RJ. 2001, Wallingford, UK: CAB International, 211-224.View ArticleGoogle Scholar
- Iovene M, Cavagnaro PF, Senalik D, Buell CR, Jiang J, Simon PW: Comparative FISH mapping of Daucus species (Apiaceae family). Chromosome Res. 2011, 19: 493-506. 10.1007/s10577-011-9202-y.View ArticlePubMedGoogle Scholar
- Gao L, Tang J, Li H, Jia J: Analysis of microsatellites in major crops assessed by computational and experimental approaches. Mol Breeding. 2003, 12: 245-261. 10.1023/A:1026346121217.View ArticleGoogle Scholar
- Morgante M, Hanafey M, Powell W: Microsatellites are preferentially associated with nonrepetitive DNA in plant genomes. Nat Genet. 2002, 30: 194-200. 10.1038/ng822.View ArticlePubMedGoogle Scholar
- Mun JH, Kim DJ, Choi HK, Gish J, Debellé F, Mudge J, Denny R, Endré G, Saurat O, Dudez AM, Kiss GB, Roe B, Young ND, Cook DR: Distribution of microsatellites in the genome of Medicago truncatula: a resource of genetic markers that integrate genetic and physical maps. Genetics. 2006, 172: 2541-2555.View ArticlePubMedPubMed CentralGoogle Scholar
- Scott KD, Eggler P, Seaton G, Rossetto M, Ablett EM, Lee LS, Henry RJ: Analysis of SSRs derived from grape ESTs. Theor Appl Genet. 2000, 100: 723-726. 10.1007/s001220051344.View ArticleGoogle Scholar
- Yu J-K, La Rota M, Kantety RV, Sorrells ME: EST Derived SSR Markers for comparative mapping in wheat and rice. Mol Genet Genomics. 2004, 271: 742-751.View ArticlePubMedGoogle Scholar
- Subramanian S, Mishra RK, Singh L: Genome-wide analysis of microsatellite repeats in humans: their abundance and density in specific genomic regions. Genome Biol. 2003, 4: R13-10.1186/gb-2003-4-2-r13.View ArticlePubMedPubMed CentralGoogle Scholar
- Toth G, Gaspari Z, Jurka J: Microsatellites in different eukaryotic genomes: survey and analysis. Genome Res. 2000, 10: 967-981. 10.1101/gr.10.7.967.View ArticlePubMedPubMed CentralGoogle Scholar
- Kruglyak SR, Durrett T, Schug MD, Aquadro CF: Equilibrium distributions of microsatellite repeat length resulting from a balance between slippage events and point mutations. Proc Natl Acad Sci USA. 1998, 95: 10774-10778. 10.1073/pnas.95.18.10774.View ArticlePubMedPubMed CentralGoogle Scholar
- Glenn TC, Schable NA: Isolating microsatellite DNA loci. Meth Enzymol. 2005, 395: 202-222.View ArticlePubMedGoogle Scholar
- Armour JAL, Neumann R, Gobert S, Jefferys AJ: Isolation of human simple repeat loci by hybridization selection. Human Mol Gen. 1994, 3: 599-605. 10.1093/hmg/3.4.599.View ArticleGoogle Scholar
- Kandpal RP, Kandpal G, Weissman SM: Construction of libraries enriched for sequence repeats and jumping clones, and hybridization selection for region-specific markers. Proc Natl Acad Sci USA. 1994, 91: 88-92. 10.1073/pnas.91.1.88.View ArticlePubMedPubMed CentralGoogle Scholar
- Kijas JMH, Fowler JCS, Garbett CA, Thomas MR: Enrichment of microsatellites from the citrus genome using biotinylated oligonucleotide sequences bound to streptavidin-coated magnetic particles. BioTechniques. 1994, 16: 656-662.PubMedGoogle Scholar
- Iniguez-Luy FL, Lukens L, Farnham MW, Amasino RM, Osborn TC: Development of public immortal mapping populations, molecular markers and linkage maps for rapid cycling Brassica rapa and B. oleracea. Theor Appl Genet. 2009, 120: 31-43. 10.1007/s00122-009-1157-4.View ArticlePubMedGoogle Scholar
- Lowe AJ, Moule C, Trick M, Edwards KJ: Efficient large-scale development of microsatellites for marker and mapping applications in Brassica crop species. Theor Appl Genet. 2004, 108: 1103-1112. 10.1007/s00122-003-1522-7.View ArticlePubMedGoogle Scholar
- Park Y-H, Alabady MS, Ulloa M, Sickler B, Wilkins TA, Yu J, Stelly DM, Kohel RJ, El-Shihy OM, Cantrell RG: Genetic mapping of new cotton fiber loci using EST-derived microsatellites in an interspecific recombinant inbred (RIL) cotton population. Mol Genet Genomics. 2005, 274: 428-441. 10.1007/s00438-005-0037-0.View ArticlePubMedGoogle Scholar
- Qureshi SN, Saha S, Kantety RV, Jenkins JN: EST-SSR: A new class of genetic markers in cotton. J Cotton Sci. 2004, 8: 112-123.Google Scholar
- Saghai Maroof MA, Biyashev RM, Yang GP, Zhang Q, Allard RW: Extraordinarily polymorphic microsatellite DNA in barley: Species diversity, chromosomal locations, and population dynamics. Proc Natl Acad Sci USA. 1994, 91: 5466-5470. 10.1073/pnas.91.12.5466.View ArticlePubMedPubMed CentralGoogle Scholar
- Fisher PJ, Richardson TE, Gardner RC: Characteristics of single- and multi-copy microsatellites from Pinus radiata. Theor Appl Genet. 1998, 96: 1069-1076. 10.1007/s001220050841.View ArticleGoogle Scholar
- Weber JL: Informativeness of human (dC-dA)n.(dG-dT)n polymorphisms. Genomics. 1990, 7: 524-530. 10.1016/0888-7543(90)90195-Z.View ArticlePubMedGoogle Scholar
- Ellegren H: Heterogeneous mutation processes in human microsatellite DNA sequences. Nat Genet. 2000, 24: 400-402. 10.1038/74249.View ArticlePubMedGoogle Scholar
- Brandström M, Ellegren H: Genome-wide analysis of microsatellite polymorphism in chicken circumventing the ascertainment bias. Genome Res. 2008, 18: 881-887. 10.1101/gr.075242.107.View ArticlePubMedPubMed CentralGoogle Scholar
- Schlötterer C, Ritter R, Harr B, Brem G: High mutation rate of a long microsatellite allele in Drosophila melanogaster provides evidence for allele-specific mutation rates. Mol Biol Evol. 1998, 15: 1269-1274.View ArticlePubMedGoogle Scholar
- Wierdl M, Dominska M, Petes TD: Microsatellite instability in yeast: dependence on the length of the microsatellite. Genetics. 1997, 146: 769-779.PubMedPubMed CentralGoogle Scholar
- Budak H, Pedraza F, Cregan PB, Baenziger PS, Dweikat I: Development and utilization of SSRs to estimate the degree of genetic relationships in a collection of pearl millet germplasm. Crop Sci. 2003, 43: 2284-2290. 10.2135/cropsci2003.2284.View ArticleGoogle Scholar
- Suwabe K, Iketani H, Nunome T, Hirai KM: Isolation and characterization of microsatellites in Brassica rapa L. Theor Appl Genet. 2002, 104: 1092-1098. 10.1007/s00122-002-0875-7.View ArticlePubMedGoogle Scholar
- Rossetto MF, Harriss CL, Mclauchlan A, Henry RJ, Baverstock PR, Lee LS: Interspecific amplification of tea tree (Melaleuca alternifolia--Myrtaceae) microsatellite loci: potential implications for conservation studies. Australian J Botany. 2000, 48: 367-373. 10.1071/BT98084.View ArticleGoogle Scholar
- Arnold C, Rossetto M, McNally J, Henry RJ: The application of SSRs characterized for grape (Vitis vinifera) to conservation studies in Vitaceae. Am J Botany. 2002, 89: 22-28. 10.3732/ajb.89.1.22.View ArticleGoogle Scholar
- Buishand JG, Gabelman WH: Investigations on the inheritance of color and carotenoid content in phloem and xylem of carrot roots (Daucus carota L.). Euphytica. 1979, 28: 611-632. 10.1007/BF00038928.View ArticleGoogle Scholar
- Bradeen JM, Simon PW: Conversion of an AFLP fragment linked to the carrot Y2 locus to a simple, codominant, PCR-based marker form. Theor Appl Genet. 1998, 97: 960-967. 10.1007/s001220050977.View ArticleGoogle Scholar
- Thiel T, Michalek W, Varshney RK, Graner A: Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.). Theor Appl Genet. 2003, 106: 411-422.PubMedGoogle Scholar
- Cavagnaro PF, Chung SM, Szklarczyk M, Grzebelus D, Senalik D, Atkins AE, Simon PW: Characterization of a deep-coverage carrot (Daucus carota L.) BAC library and initial analysis of BAC-end sequences. Mol Genet Genomics. 2009, 281: 273-288. 10.1007/s00438-008-0411-9.View ArticlePubMedGoogle Scholar
- Quackenbush J, Cho J, Lee D, Liang F, Holt I, Karamycheva S, Parvizi B, Pertea G, Sultana R, White J: The TIGR Gene Indices: analysis of gene transcript sequences in highly sampled eukaryotic species. Nucleic Acids Res. 2001, 29: 159-164. 10.1093/nar/29.1.159.View ArticlePubMedPubMed CentralGoogle Scholar
- Santos CAF, Simon PW: Merging carrot linkage groups based on conserved dominant AFLP markers in F2 populations. J Am Soc Hortic Sci. 2004, 129: 211-217.Google Scholar
- Vuylsteke M, Mank R, Antosine R, Bastiaans E, Senior ML, Stuber CW, Melchinger AE, Lubbersted T, Xia XC, Stam P, Zabeau M, Kuiper M: Two high-density AFLP linkage maps of Zea mays L.: analysis of distribution of AFLP markers. Theor Appl Genet. 1999, 99: 921-935. 10.1007/s001220051399.View ArticleGoogle Scholar
- Lander E, Green P, Abrahamson J, Barlow A, Daly M, Lincoln S, Newburg L: MAPMAKER: an interactive computer package for constructing primary genetic linkage maps of experimental and natural populations. Genomics. 1987, 1: 174-181. 10.1016/0888-7543(87)90010-3.View ArticlePubMedGoogle Scholar
- Murray M, Thompson W: Rapid isolation of high-molecular weight plant DNA. Nucleic Acids Res. 1980, 8: 4321-4325. 10.1093/nar/8.19.4321.View ArticlePubMedPubMed CentralGoogle Scholar
- Shuelke M: An economic method for the fluorescent labeling of PCR fragments. Nat Biotechnol. 2000, 18: 233-234. 10.1038/72708.View ArticleGoogle Scholar
- Excoffier L, Laval G, Schneider S: Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evol Bioinformatics Online. 2005, 1: 47-50.Google Scholar
- Just BJ: Genetic mapping of carotenoid pathway structural genes and major gene QTLs for carotenoid accumulation in wild and domesticated carrot (Daucus carota L.). PhD thesis. 2004, University of Wisconsin- Madison, Department of HorticultureGoogle Scholar
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.