Comparison of genomic SSR and EST-SSR markers for estimating genetic diversity in cucumber

Thirteen genomic microsatellite (gSSR) and sixteen expressed sequence tag (EST)-SSR (eSSR) markers were compared to estimate genetic diversity among 29 cucumber (Cucumis sativus L.) accessions. gSSR markers detected mean 4.46 alleles with a mean polymorphic information content (PIC) of 0.664, against eSSR markers with mean 3.38 alleles and a mean PIC of 0.397. gSSRs amplified more null alleles than eSSRs. Genetic diversity within the accession set was estimated by construction of dendrograms using gSSR or eSSR data. There was a clear consistency between gSSR and eSSR trees in terms of positioning of most cucumber germplasms. gSSR markers could separate various types of cucumber germplasms on the whole, although clustering of some accessions was not based on their geographical origins in eSSR tree. eSSR markers identified an independent sub-cluster containing five accessions resistant to downy mildew, suggesting a probable relationship between eSSRs and disease-resistance trait in cucumber. The Mantel test between gSSR and eSSR matrices revealed a good fit correlation (r = 0.836). The general dendrogram constructed using the combined data of gSSRs and eSSRs was similar to those obtained separately with each marker.     

Development and characterization of microsatellite markers from tropical forage Stylosanthes species and analysis of genetic variability and cross-species transferability

A limited number of functional molecular markers has slowed the desired genetic improvement of Stylosanthes species. Hence, in an attempt to develop simple sequence repeat (SSR) markers, genomic libraries from Stylosanthes seabrana B.L. Maass & 't Mannetje (2n=2x=20) using 5' anchored degenerate microsatellite primers were constructed. Of the 76 new microsatellites, 21 functional primer pairs were designed. Because of the small number of primer pairs designed, 428 expressed sequence tag (EST) sequences from seven Stylosanthes species were also examined for SSR detection. Approximately 10% of sequences delivered functional primer pairs, and after redundancy elimination, 57 microsatellite repeats were selected. Tetranucleotides followed by trinucleotides were the major repeated sequences in Stylosanthes ESTs. In total, a robust set of 21 genomic-SSR (gSSR) and 20 EST-SSR (eSSR) markers were developed. These markers were analyzed for intraspecific diversity within 20 S. seabrana accessions and for their cross-species transferability. Mean expected (He) and observed (Ho) heterozygosity values with gSSR markers were 0.64 and 0.372, respectively, whereas with eSSR markers these were 0.297 and 0.214, respectively. Dendrograms having moderate bootstrap value (23%-94%) were able to distinguish all accessions of S. seabrana with gSSR markers, whereas eSSR markers showed 100% similarities between few accessions. The set of 21 gSSRs, from S. seabrana, and 20 eSSRs, from selected Stylosanthes species, with their high cross-species transferability (45% with gSSRs, 86% with eSSRs) will facilitate genetic improvement of Stylosanthes species globally.     

Outlier Loci and Selection Signatures of Simple Sequence Repeats (SSRs) in Flax (Linum usitatissimum L.)

Genomic microsatellites (gSSRs) and expressed sequence tag-derived SSRs (EST-SSRs) have gained wide application for elucidating genetic diversity and population structure in plants. Both marker systems are assumed to be selectively neutral when making demographic inferences, but this assumption is rarely tested. In this study, three neutrality tests were assessed for identifying outlier loci among 150 SSRs (85 gSSRs and 65 EST-SSRs) that likely influence estimates of population structure in three differentiated flax sub-populations (F _ST?=?0.19). Moreover, the utility of gSSRs, EST-SSRs, and the combined sets of SSRs was also evaluated in assessing genetic diversity and population structure in flax. Six outlier loci were identified by at least two neutrality tests showing footprints of balancing selection. After removing the outlier loci, the STRUCTURE analysis and the dendrogram topology of EST-SSRs improved. Conversely, gSSRs and combined SSRs results did not change significantly, possibly as a consequence of the higher number of neutral loci assessed. Taken together, the genetic structure analyses established the superiority of gSSRs to determine the genetic relationships among flax accessions, although the combined SSRs produced the best results. Genetic diversity parameters did not differ statistically (P?>?0.05) between gSSRs and EST-SSRs, an observation partially explained by the similar number of repeat motifs. Our study provides new insights into the ability of gSSRs and EST-SSRs to measure genetic diversity and structure in flax and confirms the importance of testing for the occurrence of outlier loci to properly assess natural and breeding populations, particularly in studies considering only few loci.     

SSR and morphological trait based population structure analysis of 130 diverse flax (Linum usitatissimum L.) accessions

A total of 130 flax accessions of diverse morphotypes and worldwide origin were assessed for genetic diversity and population structure using 11 morphological traits and microsatellite markers (15 gSSRs and 7 EST–SSRs). Analysis performed after classifying these accessions on the basis of plant height, branching pattern, seed size, Indian/foreign origin into six categories called sub-populations viz. fibre type exotic, fibre type indigenous, intermediate type exotic, intermediate type indigenous, linseed type exotic and linseed type indigenous. The study assessed different diversity indices, AMOVA, population structure and included a principal coordinate analysis based on different marker systems. The highest diversity was exhibited by gSSR markers (SI = 0.46; He = 0.31; P = 85.11). AMOVA based on all markers explained significant difference among fibre type, intermediate type and linseed type populations of flax. In terms of variation explained by different markers, EST-SSR markers (12%) better differentiated flax populations compared to morphological (9%) and gSSR (6%) markers at P = 0.01. The maximum Nei's unbiased genetic distance (D = 0.11) was observed between fibre type and linseed type exotic sub-populations based on EST-SSR markers. The combined structure analysis by using all markers grouped Indian fibre type accessions (63.4%) in a separate cluster along with the Indian intermediate type (48.7%), whereas Indian accessions (82.16%) of linseed type constituted an independent cluster. These findings were supported by the results of the principal coordinate analysis. Morphological markers employed in the study found complementary with microsatellite based markers in deciphering genetic diversity and population structure of the flax germplasm.     

Simple-sequence repeat markers used in merging linkage maps of melon (Cucumis melo L.)

A set of 118 simple sequence repeat (SSR) markers has been developed in melon from two different sources: genomic libraries (gSSR) and expressed sequence-tag (EST) databases (EST-SSR). Forty-nine percent of the markers showed polymorphism between the Piel de Sapo (PS) and PI161375 melon genotypes used as parents for the mapping populations. Similar polymorphism levels were found in gSSR (51.2%) and EST-SSR (45.5%). Two populations, F₂ and a set of double haploid lines (DHLs), developed from the same parent genotypes were used for map construction. Twenty-three SSRs and 79 restriction fragment length polymorphisms (RFLPs), evenly distributed through the melon genome, were used to anchor the maps of both populations. Ten cucumber SSRs, 41 gSSRs, 16 EST-SSR, three single nucleotide polymorphism (SNP) markers, and the Nsv locus were added in the DHL population. The maps developed in the F₂ and DHL populations were co-linear, with similar lengths, except in linkage groups G1, G9, and G10. There was segregation distortion in a higher proportion of markers in the DHL population compared with the F₂, probably caused by selection during the construction of DHLs through in vitro culture. After map merging, a composite genetic map was obtained including 327 transferable markers: 226 RFLPs, 97 SSRs, three SNPs, and the Nsv locus. The map length is 1,021 cM, distributed in 12 linkage groups, and map density is 3.11 cM/marker. SSR markers alone cover nearly 80% of the map length. This map is proposed as a basis for a framework melon map to be merged with other maps and as an anchor point for map comparison between species of the Cucurbitaceae family.Electronic Supplementary Material Supplementary material is available for this article at     

Patterns of simple sequence repeats in cultivated blueberries (Vaccinium section Cyanococcus spp.) and their use in revealing genetic diversity and population structure

Blueberry (Vaccinium section Cyanococcus) is an important small fruit crop native to North America. Relationships among the primary species known as ‘blueberry’ has been a source of speculation due to the out-breeding nature of the crop, the use of intra-specific hybridization and open-pollinated selection in breeding programs, and the lack of genomic resources to adequately address the issue. The objectives of this study were to characterize simple sequence repeats (SSRs) from an emerging genomic draft sequence, develop useful molecular markers and provide an in-depth analysis of genetic diversity and population structure in blueberry using a broad range of cultivated blueberry accessions representing multiple species, ploidy levels and sources of origin. Genomic scaffolding was assembled from the whole-genome sequencing of a diploid V. corymbosum accession ‘W8520.’ From the assembled 358 Mb sequence, a total number of 43,594 SSRs were identified (122 per Mb). Among genomic regions, SSRs were longest and occurred most frequently in predicted 5′ untranslated regions (5′ UTR), while SSRs were shortest and least common in the predicted coding sequences. AG/CT and AAG/CTT were the most frequent motifs while CG/CG and CCG/CGG motifs were the rarest di- and trinucleotide motifs, respectively. For analysis of genetic diversity and population structure, 42 genomic SSR and EST–SSR markers with an average of 14.2 alleles and 56.0 allele phenotypes per locus were used to genotype a diverse blueberry population of 150 accessions. Cluster analysis grouped the 150 accessions in a manner consistent with known information regarding species, ploidy levels and pedigree. The analysis of population structure among blueberry accessions revealed inter- and intra-specific levels of stratification. Rabbiteye blueberry (V. virgatum) represents a genetically distinct subgroup within Cyanococcus. Three additional subpopulations were detected among highbush varieties that are largely attributable to distinctions between northern and southern highbush and founder effects of a single cultivar (‘Weymouth’). The identification of substructure that correlates with known pedigree information, and the availability of new genomic molecular markers will facilitate future evolutionary and genetic studies in blueberry.     

Comparative analysis of population genetic structure in Athyrium distentifolium (Pteridophyta) using AFLPs and SSRs from anonymous and transcribed gene regions

To examine the performance and information content of different marker systems, comparative assessment of population genetic diversity was undertaken in nine populations of Athyrium distentifolium using nine genomic and 10 expressed sequence tag (EST) microsatellite (SSR) loci, and 265 amplified fragment length polymorphism (AFLP) loci from two primer combinations. In range-wide comparisons (European vs. North American populations), the EST-SSR loci showed more reliable amplification and produced more easily scorable bands than genomic simple sequence repeats (SSRs). Genomic SSRs showed significantly higher levels of allelic diversity than EST-SSRs, but there was a significant correlation in the rank order of population diversities revealed by both marker types. When AFLPs, genomic SSRs, and EST-SSRs are considered, comparisons of different population diversity metrics/markers revealed a mixture of significant and nonsignificant rank-order correlations. However, no hard incongruence was detected (in no pairwise comparison of populations did different marker systems or metrics detect opposingly significant different amounts of variation). Comparable population pairwise estimates of F(ST) were obtained for all marker types, but whilst absolute values for genomic and EST-SSRs were very similar (F(ST) = 0.355 and 0.342, respectively), differentiation was consistently higher for AFLPs in pairwise and global comparisons (global AFLP F(ST) = 0.496). The two AFLP primer combinations outperformed 18 SSR loci in assignment tests and discriminatory power in phenetic cluster analyses. The results from marker comparisons on A. distentifolium are discussed in the context of the few other studies on natural plant populations comparing microsatellite and AFLP variability.     

Newly developed SSR markers reveal genetic diversity and geographical clustering in spinach (<Emphasis Type="Italic">Spinacia oleracea</Emphasis>)

Spinach is a popular leafy green vegetable due to its nutritional composition. It contains high concentrations of vitamins A, E, C, and K, and folic acid. Development of genetic markers for spinach is important for diversity and breeding studies. In this work, Next Generation Sequencing (NGS) technology was used to develop genomic simple sequence repeat (SSR) markers. After cleaning and contig assembly, the sequence encompassed 2.5% of the 980 Mb spinach genome. The contigs were mined for SSRs. A total of 3852 SSRs were detected. Of these, 100 primer pairs were tested and 85% were found to yield clear, reproducible amplicons. These 85 markers were then applied to 48 spinach accessions from worldwide origins, resulting in 389 alleles with 89% polymorphism. The average gene diversity (GD) value of the markers (based on a GD calculation that ranges from 0 to 0.5) was 0.25. Our results demonstrated that the newly developed SSR markers are suitable for assessing genetic diversity and population structure of spinach germplasm. The markers also revealed clustering of the accessions based on geographical origin with clear separation of Far Eastern accessions which had the overall highest genetic diversity when compared with accessions from Persia, Turkey, Europe, and the USA. Thus, the SSR markers have good potential to provide valuable information for spinach breeding and germplasm management. Also they will be helpful for genome mapping and core collection establishment.     

Two highly validated SSR multiplexes (8‐plex) for Euphrates' poplar,Populus euphratica (Salicaceae)

Multiplex PCR amplification of microsatellites has significantly increased the throughput and decreased the costs of genotyping. We have developed two highly polymorphic microsatellite multiplexes for Populus euphratica, the only tree species found in desert regions of Western China and adjacent Central Asian countries. The first of these multiplex kits comprises an eight‐Plex of genomic SSRs (gSSRs) obtained from published databases. The second comprises an eight‐plex of newly designed EST‐SSRs (eSSRs) based on expressed sequence tags for P. euphratica. Both kits were tested on a sample of 170 individuals from four populations. The gSSRs exhibited slightly more polymorphism than the eSSRs. The new multiplex protocols yielded consistent results in the hands of multiple researchers, demonstrating their robustness. The 16 loci used in the kits exhibited a high transferability rate (82.0%) in eight other poplar species belonging to five different sections of the genus. Both kits should therefore be useful for further investigations of population genetics in P. euphratica and related species. Our results indicate that it is essential to follow recently established recommendations when developing microsatellite markers, including verifying the amplification efficiency, detecting null alleles and carefully measuring error rates.     

Genetic assessment of safflower (Carthamus tinctorius L.) collection with microsatellite markers acquired via pyrosequencing method

A genetic evaluation of safflower germplasm collections derived from different geographical regions and countries will provide useful information for sustainable conservation and the utilization of genetic diversity. However, the molecular marker information is limited for evaluation of genetic diversity of safflower germplasm. In this study, we acquired 509 putative genomic SSR markers for sufficient genome coverage using next‐generation sequencing methods and characterized thirty polymorphic SSRs in safflower collection composed of 100 diverse accessions. The average allele number and expected heterozygosity were 2.8 and 0.386, respectively. Analysis of population structure and phylogeny based on thirty SSR profiles revealed genetic admixture between geographical regions contrary to genetic clustering. However, the accessions from Korea were genetically conserved in distinctive groups in contrast to other safflower gene pool. In conclusion, these new genomic SSRs will facilitate valuable studies to clarify genetic relationships as well as conduct population structure analyses, genetic map construction and association analysis for safflower.     

Development of eSSR-Markers in Setaria italica and Their Applicability in Studying Genetic Diversity,Cross-Transferability and Comparative Mapping in Millet and Non-Millet Species

Foxtail millet ( Setaria italica L.) is a tractable experimental model crop for studying functional genomics of millets and bioenergy grasses. But the limited availability of genomic resources, particularly expressed sequence-based genic markers is significantly impeding its genetic improvement. Considering this, we attempted to develop EST-derived-SSR (eSSR) markers and utilize them in germplasm characterization, cross-genera transferability and in silico comparative mapping. From 66,027 foxtail millet EST sequences 24,828 non-redundant ESTs were deduced, representing ~16 Mb, which revealed 534 (~2%) eSSRs in 495 SSR containing ESTs at a frequency of 1/30 kb. A total of 447 pp were successfully designed, of which 327 were mapped physically onto nine chromosomes. About 106 selected primer pairs representing the foxtail millet genome showed high-level of cross-genera amplification at an average of ~88% in eight millets and four non-millet species. Broad range of genetic diversity (0.02–0.65) obtained in constructed phylogenetic tree using 40 eSSR markers demonstrated its utility in germplasm characterizations and phylogenetics. Comparative mapping of physically mapped eSSR markers showed considerable proportion of sequence-based orthology and syntenic relationship between foxtail millet chromosomes and sorghum (~68%), maize (~61%) and rice (~42%) chromosomes. Synteny analysis of eSSRs of foxtail millet, rice, maize and sorghum suggested the nested chromosome fusion frequently observed in grass genomes. Thus, for the first time we had generated large-scale eSSR markers in foxtail millet and demonstrated their utility in germplasm characterization, transferability, phylogenetics and comparative mapping studies in millets and bioenergy grass species.     

Development of genome-wide SSR markers in melon with their cross-species transferability analysis and utilization in genetic diversity study

Simple sequence repeats (SSRs) are one of the most informative and widely used molecular markers in plant research. The melon draft genome has provided a powerful tool for SSR marker development in this species in which there are still not enough SSR markers. We therefore developed genome-wide SSR markers from melon, which were used for genetic diversity analysis in melon accessions and comparative mapping with cucumber and watermelon. A total of 44,265 microsatellites from the melon genome were characterized, of which 28,570 SSR markers were developed. In silico PCR analysis with these SSR markers identified 4002 and 1085 with one amplicon in cucumber and watermelon genome, respectively. With these cross-species transferable melon SSR markers, the chromosome synteny between melon and cucumber as well as watermelon was established, which revealed complicated mosaic patterns of syntenic blocks among them. We experimentally validated 384 SSR markers, from which 42 highly informative SSR ones were selected for genetic diversity and population structure analysis among 118 melon accessions. The large number of melon SSR markers developed in this study provides a valuable resource for genetic linkage map construction, molecular mapping, and marker-assisted selection (MAS) in melon. Furthermore, the cross-species transferable SSR markers could also be useful in various molecular marker-related studies in other closely related species in Cucurbitaceae family in which draft genomes are not yet available.     

Utility of EST-derived SSRs as population genetics markers in a beetle

Microsatellite, or simple sequence repeat (SSR), loci can be identified by mining expressed sequence tag (EST) databases, and where these are available, marker development time and expense can be decreased considerably over conventional strategies of probing the entire genome. However, it is unclear whether they provide information on population structure similar to that generated by anonymous genomic SSRs. We performed comparative population genetic analyses between EST-derived SSRs (EST-SSRs) and anonymous SSRs developed from genomic DNA for the same set of populations of the insect Diabrotica virgifera, a beetle in the family Chrysomelidae. Compared with noncoding, nontranscribed regions, EST-SSRs were generally less polymorphic but had reduced occurrence of null alleles and greater cross-species amplification. Neutrality tests suggested the loci were not under positive selection. Across all populations and all loci, the genomic and EST-SSRs performed similarly in estimating genetic diversity, F(IS), F(ST), population assignment and exclusion tests, and detection of distinct populations. These findings, therefore, indicate that the EST-SSRs examined can be used with confidence in future genetic studies of Diabrotica populations and suggest that EST libraries can be added as a valuable source of markers for population genetics studies in insects and other animals.     

Mapping quantitative trait loci for important agronomic traits in finger millet (<Emphasis Type="Italic">Eleusine coracana</Emphasis>) mini core collection with genomic and genic SSR markers

Allele identification for agro-morphological traits and stress resistance is a major concern across the globe for improving productivity of finger millet. Here, we used 46 genomic and 58 genic simple sequence repeats (SSRs) markers in a set of 66 accessions used to constitute a global mini-core collection for analysing their genetic structure as a population and establishing association among markers and twenty morphological traits including resistance to finger blast. Phenotypic data revealed a wide range of variation for all traits except flag leaf width and flag leaf sheath width. We got amplification of 81 alleles by the 31 genomic SSRs at an average of 2.61 alleles per locus. Polymorphism information content (PIC) values varied from 0.21 to 0.75 and average gene diversity was 0.49. Structure analysis of the population using the genomic SSR data divided the accessions into two clusters where Indian and exotic accessions were grouped in separate clusters. Genic SSRs which were associated with blast resistance genes, amplified 36 alleles at an average of 2 alleles per locus. PIC values ranged from 0.32 to 0.37 and average gene diversity was 0.45. Population structure analysis using data from these SSRs grouped the accessions into three clusters, which broadly correspond to their reaction to blast disease. Twenty-two significant associations were found using the GLM approach for 20 agro-morphological traits both in 2012 and 2014, while, 7 and 5 significant marker-trait associations were identified using MLM in 2012 and 2014 respectively. The SSR markers FMBLEST35 and FMBLEST36 designed from the Pi21 gene sequence of rice were found to be associated with blast disease resistance in finger millet indicating that the gene homologues play a significant role in an important role for neck blast resistance.     

Comparative analysis of genetic diversity in the mangrove species Avicennia marina (Forsk.) Vierh. (Avicenniaceae) detected by AFLPs and SSRs

Avicennia marina is an important mangrove species with a wide geographical and climatic distribution which suggests that large amounts of genetic diversity are available for conservation and breeding programs. In this study we compare the informativeness of AFLPs and SSRs for assessing genetic diversity within and among individuals, populations and subspecies of A. marina in Australia. Our comparison utilized three SSR loci and three AFLP primer sets that were known to be polymorphic, and could be run in a single analysis on a capillary electrophoresis system, using different- colored fluorescent dyes. A total of 120 individuals representing six populations and three subspecies were sampled. At the locus level, SSRs were considerably more variable than AFLPs, with a total of 52 alleles and an average heterozygosity of 0.78. Average heterozygosity for AFLPs was 0.193, but all of the 918 bands scored were polymorphic. Thus, AFLPs were considerably more efficient at revealing polymorphic loci than SSRs despite lower average heterozygosities. SSRs detected more genetic differentiation between populations (19 vs 9%) and subspecies (35 vs 11%) than AFLPs. Principal co-ordinate analysis revealed congruent patterns of genetic relationships at the individual, population and subspecific levels for both data sets. Mantel testing confirmed congruence between AFLP and SSR genetic distances among, but not within, population comparisons, indicating that the markers were segregating independently but that evolutionary groups (populations and subspecies) were similar. Three genetic criteria of importance for defining priorities for ex situ collections or in situ conservation programs (number of alleles, number of locally common alleles and number of private alleles) were correlated between the AFLP and SSR data sets. The congruence between AFLP and SSR data sets suggest that either method, or a combination, is applicable to expanded genetic studies of mangroves. The codominant nature of SSRs makes them ideal for further population-based investigations, such as mating-system analyses, for which the dominant AFLP markers are less well suited. AFLPs may be particularly useful for monitoring propagation programs and identifying duplicates within collections, since a single PCR assay can reveal many loci at once. Received: 3 October 2000 / Accepted: 19 February 2001     

Development of genomic simple sequence repeat markers in faba bean by next-generation sequencing

Faba bean (Vicia faba L.) is an important food legume crop with a huge genome. Development of genetic markers for faba bean is important to study diversity and for molecular breeding. In this study, we used Next Generation Sequencing (NGS) technology for the development of genomic simple sequence repeat (SSR) markers. A total of 14,027,500 sequence reads were obtained comprising 4,208 Mb. From these reads, 56,063 contigs were assembled (16,367 Mb) and 2138 SSRs were identified. Mono and dinucleotides were the most abundant, accounting for 57.5 % and 20.9 % of all SSR repeats, respectively. A total of 430 primer pairs were designed from contigs larger than 350 nucleotides and 50 primers pairs were tested for validation of SSR locus amplification. Nearly all (96 %) of the markers were found to produce clear amplicons and to be reproducible. Thirty-nine SSR markers were then applied to 46 faba bean accessions from worldwide origins, resulting in 161 alleles with 87.5 % polymorphism, and an average of 4.1 alleles per marker. Gene diversity (GD) of the markers ranged from 0 to 0.48 with an average of 0.27. Testing of the markers showed that they were useful in determining genetic relationships and population structure in faba bean accessions.     

Development of SSR markers for studies of diversity in the genus <Emphasis Type="Italic">Fagopyrum</Emphasis>

The numbers of SSR markers and their utilization have not been determined and investigated as extensively in Fagopyrum species as compared to other crop species. The current report presents 136 new SSR markers in Fagopyrum esculentum ssp. esculentum and their application to related species in the genus Fagopyrum. Of the 136 SSRs, 10 polymorphic SSR markers were utilized in a genetic diversity analysis of a common buckwheat population consisting of 41 accessions of diverse origin. The study showed observed (H _O) and expected (H _E) heterozygosities ranging from 0.071 to 0.924 (mean = 0.53) and from 0.073 to 0.902 (mean = 0.412), respectively. Forty-one of the 136 SSRs amplified sequences in other Fagopyrum species, including the cymosum and urophyllum groups. The phylogenetic relationships revealed using the SSRs was consistent with results obtained using other marker systems, with one exception. The sequence and diversity information obtained using these new SSRs and their cross-transferability to related Fagopyrum species will increase our understanding of genetic structures and species relationships within the Fagopyrum genus.     

High SSR diversity but little differentiation between accessions of Nordic timothy (Phleum pratense L.)

A large collection of genebank accessions of the hexaploid outcrossing forage grass species timothy (Phleum pratense L.) was for the first time analysed for SSR diversity on individual, population and regional level. Timothy is the most important forage grass species in the Nordic countries. Eighty-eight timothy accessions from Nordic countries and eight accessions around Europe were analysed with recently developed simple sequence repeat (SSR) markers. Timothy proved to be very polymorphic: the 13 selected SSRs amplified a total of 499 polymorphic alleles, the number of alleles per SSR locus varying from 15 to 74. Taking all SSR alleles together, the observed number in each accession ranged from 95 to 203. Levels of diversity were found to be significantly different between countries, vegetation zones and different cultivar types. However, the differentiation between accessions was low: most of the variation (94%) in the studied timothy material was due to variation within accessions and only 5% was between accessions and 1% between countries. Lack of geographical differentiation may reflect the outcrossing and hexaploid nature of timothy. Our results showed that neutral SSR markers are suitable for demonstrating levels of diversity but not alone adequate to resolve population structure in timothy. Nordic timothy material seems to be diverse enough for breeding purposes and no decline in the level of diversity was observed in varieties compared to wild timothy populations. Challenges in analysing SSR marker data in a hexaploid outcrosser were discussed.