Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Genome-Wide Distribution and Organization of Microsatellites in Plants: An Insight into Marker Development in Brachypodium

  • Humira Sonah ,

    Contributed equally to this work with: Humira Sonah, Rupesh K. Deshmukh

    Affiliation National Research Centre on Plant Biotechnology, Indian Agricultural Research Institute (IARI), New Delhi, India

  • Rupesh K. Deshmukh ,

    Contributed equally to this work with: Humira Sonah, Rupesh K. Deshmukh

    Affiliations National Research Centre on Plant Biotechnology, Indian Agricultural Research Institute (IARI), New Delhi, India, Swami Ramanand Teerth Marathwada University, Nanded, Maharashtra, India

  • Anshul Sharma,

    Affiliation National Research Centre on Plant Biotechnology, Indian Agricultural Research Institute (IARI), New Delhi, India

  • Vinay P. Singh,

    Affiliation National Research Centre on Plant Biotechnology, Indian Agricultural Research Institute (IARI), New Delhi, India

  • Deepak K. Gupta,

    Affiliation National Research Centre on Plant Biotechnology, Indian Agricultural Research Institute (IARI), New Delhi, India

  • Raju N. Gacche,

    Affiliation Swami Ramanand Teerth Marathwada University, Nanded, Maharashtra, India

  • Jai C. Rana,

    Affiliation National Bureau of Plant Genetic Resources, Regional Station, Shimla (HP), India

  • Nagendra K. Singh,

    Affiliation National Research Centre on Plant Biotechnology, Indian Agricultural Research Institute (IARI), New Delhi, India

  • Tilak R. Sharma

    trsharma@nrcpb.org

    Affiliation National Research Centre on Plant Biotechnology, Indian Agricultural Research Institute (IARI), New Delhi, India

Abstract

Plant genomes are complex and contain large amounts of repetitive DNA including microsatellites that are distributed across entire genomes. Whole genome sequences of several monocot and dicot plants that are available in the public domain provide an opportunity to study the origin, distribution and evolution of microsatellites, and also facilitate the development of new molecular markers. In the present investigation, a genome-wide analysis of microsatellite distribution in monocots (Brachypodium, sorghum and rice) and dicots (Arabidopsis, Medicago and Populus) was performed. A total of 797,863 simple sequence repeats (SSRs) were identified in the whole genome sequences of six plant species. Characterization of these SSRs revealed that mono-nucleotide repeats were the most abundant repeats, and that the frequency of repeats decreased with increase in motif length both in monocots and dicots. However, the frequency of SSRs was higher in dicots than in monocots both for nuclear and chloroplast genomes. Interestingly, GC-rich repeats were the dominant repeats only in monocots, with the majority of them being present in the coding region. These coding GC-rich repeats were found to be involved in different biological processes, predominantly binding activities. In addition, a set of 22,879 SSR markers that were validated by e-PCR were developed and mapped on different chromosomes in Brachypodium for the first time, with a frequency of 101 SSR markers per Mb. Experimental validation of 55 markers showed successful amplification of 80% SSR markers in 16 Brachypodium accessions. An online database ‘BraMi’ (Brachypodium microsatellite markers) of these genome-wide SSR markers was developed and made available in the public domain. The observed differential patterns of SSR marker distribution would be useful for studying microsatellite evolution in a monocot–dicot system. SSR markers developed in this study would be helpful for genomic studies in Brachypodium and related grass species, especially for the map based cloning of the candidate gene(s).

Introduction

Microsatellites or simple sequence repeats (SSRs) are co-dominant, abundant, multi-allelic, and uniformly distributed over the genome, and can be detected by simple reproducible assays [1]. These important features have made microsatellites the markers of choice for marker-assisted plant breeding, DNA fingerprinting of genetic resources, molecular mapping and map based cloning of specific genes. Microsatellite markers have also been used in several studies to define conserved regions among related species [2][4]. Initially, SSR markers were developed from expressed sequence tags (ESTs) and bacterial artificial chromosome (BAC) end sequences in most plant species. However, whole genome sequencing has led to the identification of numerous SSR markers that are distributed over the entire genomes of rice and Arabidopsis. The mechanism of microsatellite evolution and their genome-wide distribution, however, are still not well studied in plants mostly due to the lack of genomic information. The recently sequenced genomes of Brachypodium [5], Populus [6], Medicago (www.jgi.doe.gov) and sorghum [7] along with the already well-characterized sequenced genomic information available for rice [8] and Arabidopsis [9], will facilitate the comparative genomics studies in plants. Besides understanding genome organization, sequenced genomes can be effectively used for the generation of molecular markers and their cross species utilization, specifically for those species where very little or no genomic information is available.

Thousands of SSR markers have been developed in rice and their application in other grass species has been proven and well documented [10]-[12]. However, it is postulated that the Brachypodium genome would exhibit a higher level of collinearity to the genomes of temperate cereal crops as compared to the rice genome. Therefore, SSR markers developed in Brachypodium may be utilized more effectively in wheat compared to rice SSR markers. Moreover, several features such as a typical plant structure, growth habit, rapid generation time, self pollination and a compact genome size (∼300 Mb) make Brachypodium an excellent model system for functional and structural genomic studies in cereals and grasses [13][14]. Several efforts have been made to make Brachypodium a more useful and convenient model system for genomic studies [14]. However, to develop any model system for genomic studies, a dense molecular genetic linkage map and genome-wide distributed molecular markers are required. Recently, a genetic linkage map of 139 marker loci for Brachypodium was developed using F2 population derived from a cross between diploid lines Bd3-1 and Bd21 [15]. The map was constructed using SSR markers derived from the EST, BAC end sequences and information available on conserved orthologous sequence (COS) in other grass species. Currently, efforts are underway to develop a dense molecular genetic linkage map for Brachypodium using different types of markers. However, a genome-wide resource of a large number of informative molecular markers is required to supplement this effort.

The present study deals with a genome-wide comparative analysis of microsatellite distribution in the nuclear and chloroplast genomes of monocot and dicot species, the development of genome-wide SSR markers, and the validation of a subset of these new markers in Brachypodium.

Results and Discussion

Microsatellite distribution in monocot and dicot species

A total of 797,863 SSRs were identified among six plant genomes — three monocot (Brachypodium, sorghum and rice) and three dicot (Arabidopsis, Medicago and Populus) plant species (Table 1). Among the six genomes analyzed, the maximum number and frequency of SSRs were obtained from Populus followed by Medicago whereas the sorghum genome had the lowest frequency. The frequency of SSRs was considerably higher among dicots compared to monocots. Among monocots, the frequency of SSRs in the rice genome was nearly twice that in sorghum and Brachypodium genomes (Table 1). Since the six selected plant species belonged to very diverse groups of monocots and dicots, the distribution pattern of SSR motifs, with specific sequences in these genomes, was not uniform. However, the overall pattern of SSR motifs of particular lengths was similar. Mono-nucleotide repeats dominated over other type of repeats in all the six plant species. However, the frequency of SSRs decreased stepwise with increase in motif length (mono- to hexa- nucleotide repeats) except in Brachypodium where the frequency of tri- nucleotide repeats was higher than that of the di-nucleotide repeats (Figure S1). Mono-nucleotide repeats were found to be minimum (43%) in sorghum and maximum (79%) in Medicago genomes. While the mono-, di- and tri-nucleotide repeats mostly contributed to the major proportion of SSRs, a very small share was contributed by tetra-, penta- and hexa-nucleotide repeats. A maximum of 5.4% contribution of tetra-, penta- and hexa- nucleotide repeats, was observed in the sorghum genome. A similar trend was observed for other genomes studied in the present investigation (Table 1).

thumbnail
Table 1. Distribution of microsatellite with respect to motif length and genome size in monocot and dicot plant species.

https://doi.org/10.1371/journal.pone.0021298.t001

Among the two types of mono-nucleotide repeats, (A/T)n was the most abundant in all the plant species while (G/C)n was comparatively scarce (Table 2). In the mono-nucleotide repeats category, the maximum (99%) A/T repeats were present in the Arabidopsis genome and the minimum (78%) in the Brachypodium genome. In the di-nucleotide repeat category, the distribution of SSRs in different motif types was not uniform and the most frequent motif type was different for each plant species. For example, AG/CT repeats were more frequent in Brachypodium and rice with 50.7% and 41.9% frequency, respectively; whereas AT/AT repeats were more frequent in Populus (60.5%) and Medicago (59.9%). In rice, both AG/CT and AT/AT repeats dominated other di-nucleotide repeats. Interestingly, the CG/CG motif contributed less than 0.5% in dicots, whereas it was 3.1%–7.0% in all di-nucleotide repeats identified in the monocots. The analysis of mono- and di-nucleotide repeats concluded that CG-rich motifs were least preferred in both monocot and dicot genomes. However, for tri-nucleotide repeats the AGC/CGT, AGG/CCT and CCG/CGG were observed more frequently in all the monocot species, whereas A/T-rich repeats, such as AAC/GTT, AAG/CTT and AAT/ATT, were preferred in dicots (Table 2). The frequency of tetra-, penta- and hexa-nucleotide repeats was very low in all the plant genomes investigated in the present study and their motif-wise distribution was not significant across the genomes.

thumbnail
Table 2. Frequency of different types of motifs in a class of microsatellites with mono-, di-, and tri-nucleotide repeats analysed in monocot and dicot plant genomes.

https://doi.org/10.1371/journal.pone.0021298.t002

The dominant occurrence of repeat motifs, of a particular sequence and length, in plant genomes is the outcome of selection pressures applied on that specific motif during evolution. The molecular mechanisms for the origin of microsatellites are not completely understood. The most common mutational mechanism affecting microsatellites is replication slippage, a process involving addition or removal of one or more motif repeats; however other mechanisms, such as unequal crossing over, nucleotide substitutions, or duplication events, have also been considered to be responsible for microsatellite variations [16][18]. However, these theories cannot explain the species-specific accumulation of particular motif repeats observed in the present study. Other factors, such as codon preference, DNA replication and the mismatch repair system, as well as structural and functional attributes of genomes that are unique to the species or for the particular taxon, may be responsible for the unique microsatellite distribution patterns in plant genomes. Moreover, the SSR length, motif structure and G/C content of a genome are considered to be factors influencing microsatellite evolution [19][21]. Polymorphism among SSRs is a repeat length polymorphism due to repeat elongation/shortening events, which indicates that such processes are important factors for molecular evolution. The repeat elongation/shortening processes also lead to increase in biological complexity, which is a characteristic of biological evolution. It is known that SSRs within genes are substantially involved in the regulation of evolutionary processes as SSRs in the protein-coding regions can lead to a gain or loss of gene function. Earlier, sequence variations in genomes, particularly in microsatellite distribution, were supported by the theory of stabilization patterns and potential secondary structures, as well as factors such as the mismatch repair enzymes [22][24]. All these theories, which deal with a particular factor being responsible for sequence preference, were suggested based on very limited knowledge and lack of genome-wide information on a large variety of genomes. Till date, very little work has been done to propose a genome-wide mechanism for the selection of microsatellite motifs with a particular sequence. These mechanisms may be further illustrated with available genomic resources, and the data presented in this paper would definitely help in the understanding of microsatellite evolution in the genomes of plant species.

Microsatellite distribution in coding regions

Microsatellites were identified in the coding DNA sequences (CDS) of six plant species to study the pattern of distribution in the coding regions of monocots and dicots. A total of 36,585 SSRs were identified in the 238,798 CDS of about 269.3 Mb size data for all the six plant species included in this study. Interestingly, the frequency of SSRs observed in the CDS region (CDS-SSRs) of monocots was twice that observed in the CDS of dicots (Table S1). The highest frequency (203.7 SSR/Mb) of CDS-SSRs was identified in rice followed by sorghum (181.1 SSR/Mb) whereas the lowest frequency (68.1 SSR/Mb) was observed in Populus. Tri-nucleotide repeats were found to be most abundant among the microsatellites in the coding region of plant genomes (Figure 1) and contributed to about 93% of SSRs in monocots and about 76% of SSRs in dicots. Such an accumulation of tri-nucleotide repeats in the coding regions was mostly due to the triplet-repeat nature of the codon. Mono-nucleotide repeats contributed about 2% of such SSRs in monocots while in dicots it was 14.2%; this variation may be due to the more frequent occurrence of A/T repeats in dicots. Moreover, G/C-rich repeats in the CDS region of monocots were identified with much more frequency than in dicots (Table 3). In the category of mono-nucleotide repeats, A/T repeats dominated over G/C repeats in both monocots and dicots except in rice where G/C repeats contributed to 53.1% of mono-repeat SSRs. Although in monocots G/C repeats were slightly less than A/T repeats, these were very infrequent in dicots, and as low as 3.4% in Arabidopsis (Table 3). Interestingly, in the category of di-nucleotide repeats, GC/CG repeats were predominant in monocots with an average of 49.9%, while they were completely absent in dicots. The AG/CG repeats accounted for an average of 41.1% SSRs in monocots and were also major contributors in dicots with an average of 73.5% di-nucleotide SSRs. Tri-nucleotides with repeat motifs CCG/GGC were dominant only in monocots and contributed to 51.5% of the total tri-nucleotide repeats identified in monocots, whereas only 1.9% of these repeats were present in dicots. Tri-nucleotides with repeat motifs AAG/CTT accounted for 29.5% of the total tri-nucleotide repeats in dicots. The microsatellite distribution pattern in the CDS region was found to be very unique for the monocot and dicot species. Functional annotation of Brachypodium CDS sequences containing CG-rich repeats revealed that about 50% of these genes were involved in binding activity (Table S2).

thumbnail
Figure 1. Distribution of microsatellites in coding DNA sequences (CDS) of six plant species with respect to motif length.

Microsatellites were identified with criteria of mono- to hexa-nucleotides motifs using MISA software tool, and the minimum repeat unit was defined as 10 for mono-, 6 for di-, and 5 for tri-, tetra-, penta-, and hexa-nucleotides.

https://doi.org/10.1371/journal.pone.0021298.g001

thumbnail
Table 3. Frequency of microsatellite motifs in coding DNA sequence of six plant species.

https://doi.org/10.1371/journal.pone.0021298.t003

Monocots and dicots are thought to have diverged from a common ancestor approximately 200 million years ago [25]. In several comparative genomic studies, Arabidopsis and rice have been considered as models for dicots and monocots, respectively. However, numerous interesting findings have emerged while comparing these two genomes; for example, rice genes are longer and GC-rich than Arabidopsis genes [26]. Though Arabidopsis has the smallest genome among the dicot species, it is thought to have evolved by chromosomal duplication; while the rice genome, which is comparatively larger than the Arabidopsis genome, showed more duplication [26][29]. The GC-rich monocot genomes may have microsatellites with GC-rich motifs whereas dicots lack GC-rich motifs. The relationship between microsatellite evolution and chromosomal duplications has not been well studied. The duplicated regions are thought to have different selection pressures than other regions, which may be a reason for motif preference and frequency in monocots and dicots. Such a biased selection of SSRs was observed in Populus where most of the SSRs in the coding regions are missing at the duplicated chromosomal segment mostly due to loss of corresponding genes [30]. This emphasizes the role of microsatellites in gene and genome evolution. Although, a systematic study on a number of genomes is required to make any definite conclusions, recent developments in sequencing technology and the availability of an increasing number of genome sequences for analysis would definitely provide a basis for the study of microsatellite evolution in plants [31].

SSRs frequency in plant chloroplast genomes

A total of 337 SSRs were identified for the chloroplast genome of the six plant species analyzed in this study. The highest frequency of SSRs was identified in the chloroplast genome of Populus followed by Medicago and Arabidopsis (Table S3). Compared with dicots, monocots had very infrequent SSRs in the chloroplast genome. Most of the SSRs identified in the plant chloroplast genome were mono-nucleotide repeats that contributed about 92.5% of the total SSRs (Figure 2). In the mono-nucleotide category of repeats, A/T contributed 97.4% of the total repeats. For nuclear genomes, G/C repeats were predominantly found in the sorghum and Brachypodium chloroplast genomes. Although, the number of di- and tri-nucleotide repeats identified in the chloroplast genome is not enough to compare patterns in the chloroplasts of monocots and dicots, in a broader context the chloroplasts of dicots were richer in SSRs with di- and tri-nucleotide repeats, which were otherwise lacking in monocots.

thumbnail
Figure 2. Distribution of microsatellites in the chloroplast genome of six plant species with respect to motif type.

Microsatellites were identified with criteria of mono- to hexa-nucleotides motifs using MISA software tool, and the minimum repeat unit was defined as 10 for mono-, 6 for di-, and 5 for tri-, tetra-, penta-, and hexa-nucleotides.

https://doi.org/10.1371/journal.pone.0021298.g002

SSRs frequency in Brachypodium genomes

A total of 51,875 SSRs were identified in the 271 Mb sequence of the Brachypodium genome with a maximum of 30,573 mono-, followed by 10,625 di- and 9,407 tri-nucleotide repeats (Table 4). The penta- and hexa-nucleotide repeats were only 1,270, which represented 2.5% of the total SSRs identified in the Brachypodium genome. Chromosome 4 of Brachypodium contained the maximum frequency (197 SSRs/Mb) of SSRs, whereas chromosome 5 contained the minimum frequency (175 SSRs/Mb) of SSRs. The frequency of SSRs in the different chromosomes was almost uniform and the overall frequency of SSRs in Brachypodium was 191 SSRs/Mb.

thumbnail
Table 4. Chromosome-wide distribution of microsatellites in the B. distachyon genome.

https://doi.org/10.1371/journal.pone.0021298.t004

Interestingly, the frequency of SSRs observed on the short arm of chromosome 5 was much lower than that on the long arm. This low frequency of SSRs in the short arm was common for all types of motifs (Figure S2). The short arm of chromosome 5 (Bd5s) has several features that are different from the rest of the chromosomes [5]. These include a low gene density (roughly half of the rest of the chromosomes); a high LTR retrotransposon density with the youngest intact Gypsy elements; and the lowest solo LTR density. These attributes may be responsible for the low frequency of microsatellites in Bd5s.

Development of genome-wide SSR markers and their validation in Brachypodium genotypes

A total of 27,329 SSR markers (including 7,225 class I and 20,104 class II) were successfully designed. Of these 22,879 (83.7%) SSR markers were validated by e-PCR (Table 5, Table S4). A subset of 44 (80%) markers was amplified in 16 Brachypodium genotypes with prominent PCR products of expected size (Figure 3, Table S5).

thumbnail
Figure 3. PCR-amplified product of marker (a) B2M569 and (b) B2M1 for the 16 Brachypodium genotypes resolved on the 3% metaphor agarose gel.

https://doi.org/10.1371/journal.pone.0021298.g003

thumbnail
Table 5. Summary of the primers designed for microsatellites of different motif types in B. distachyon genome.

https://doi.org/10.1371/journal.pone.0021298.t005

The frequency of SSR markers was 101 per Mb, covering the entire genome with very fine gaps of less than 10 Kb (Figure 4). Most of these gaps were found at the centromeric regions. The primer set also included 2,810 compound SSRs, which would provide better polymorphism in Brachypodium (Table 5). The rationale for categorizing SSRs as Class I and Class II was that the SSRs with the larger number of repeats (Class I) were found to be more polymorphic than those with lesser number of repeats (Class II) in animals [32][33] as well as in plants [34]. These results from animals and plants were concluded on the basis of polymorphism observed on gel electrophoresis. However, recent comparisons made between available Indica and Japonica rice genome sequences revealed that Class II SSRs were more polymorphic than Class I SSRs, but due to their small size difference polymorphism could not be achieved on agarose gel [12]. However, polymorphism can be achieved by running PCR product on polyacrylamide gel electrophoresis (PAGE) or MetaPhor® Agarose. Besides this, advancements in techniques and frequent use of capillary electrophoresis can help exploit the power of Class II SSR markers. The list of Class I and Class II markers designed in Brachypodium, are provided separately in the public domain (http://125.18.242.19/plantgenomedb/brachypodium_markersearch.jsp).

thumbnail
Figure 4. Distribution of SSR markers developed for the five chromosomes of Brachypodium distachyon genome.

Markers were presented on the left hand side of the chromosome bar as per the physical location and named (e.g. B01d12) using prefix B for Brachypodium followed by chromosome number then d for distachyon, followed by the number of SSRs as per their order on the chromosome. The left vertical bar is the physical length ruler scaling in Mb.

https://doi.org/10.1371/journal.pone.0021298.g004

Brachypodium is a member of the Poaceae family, which comprises over 600 genera and more than 10,000 species that dominate many ecological and agricultural systems. Brachypodium is related to major cereal grain crops, such as wheat, barley, oat, maize, rice, rye, sorghum and millet. It has a small genome (∼272 Mbp), small physical stature, is self-fertile, has a short lifecycle and requires normal growth conditions; it also has an efficient transformation system [35][36]. All these qualities make it an excellent model organism for genomic research in temperate grasses and cereals. However, successful implementation of Brachypodium as a model plant for genomic studies in grasses requires genome-wide well-characterized molecular markers. But, only 398 SSR makers have been developed from BAC end and EST databases in Brachypodium [37]. The EST-derived SSR markers represent the coding region that would be important for gene tagging and cloning, however due to their low potential for polymorphism, it had marginal utility in mapping and diversity related studies [38]-[40]. In the present investigation, the large number of SSR markers developed from the genome sequence of Brachypodium would serve as an important genomic resource for use in many cereal crops. We have made a maiden attempt to provide a large list of SSR markers along with their primer information in the public domain.

Database of SSR markers and genome-wide localization

To facilitate access to and utilization of SSR markers in Brachypodium, we developed a Brachypodium microsatellite marker (BraMi) database, which includes details of 27,392 SSR markers. The database is hosted on the NRCPB server and is accessible at http://125.18.242.19/plantgenomedb/brachypodium_markersearch.jsp (Figure S3). The database has searchable tools to get information for all SSR markers, and useful web links of other databases, websites and research institutes related to Brachypodium research.

In the present investigation, we made genome-wide comparisons of SSR frequencies in monocot and dicot plants with respect to repeat motif length and sequence. Among mono-nucleotide repeats, (A/T)n was most abundant among all the plant species compared to (G/C)n. Surprisingly in the Arabidopsis genome, 99% of the mono-repeats were (A/T)n. Interestingly, since the CG/CG motif in the di- nucleotide repeats category contributed to less than 0.5% SSRs in dicots, it may be concluded that the CG-rich motif was least preferred in both monocot and dicot genomes. Several interesting features revealed in this study would definitely help enhance the understanding of microsatellite evolution in plants and its relationship to the divergence of monocots and dicots. In addition, the genome-wide marker resource developed in this study would be helpful for genomic studies in Brachypodium and other related grass species.

Materials and Methods

Identification of microsatellites

Genomic (pseudo molecules) and CDS sequences of Brachypodium, sorghum, Populus and Medicago were downloaded from the Phytozome database (www.phytozome.net/). Similarly the rice genome sequence from TIGR database (http://rice.plantbiology.msu.edu/) and Arabidopsis from TAIR database (www.arabidopsis.org/) were retrieved in batches. The Perl script MIcroSAtelitte (MISA) was used to identify microsatellites in all these genomes (http://pgrc.ipk-gatersleben.de/misa/). To identify the presence of SSRs, only 1 to 6 nucleotides motifs were considered, and the minimum repeat unit was defined as 10 for mono-, 6 for di-, and 5 for tri-, tetra-, penta- and hexa-nucleotides. Compound SSRs were defined as ≥2 SSRs interrupted by ≤100 bases.

Designing SSR based primers for Brachypodium

The SSR information generated by MISA was used for designing primers flanking the repeats. To design primers flanking the microsatellite loci, two perl scripts were used as interface modules for the program-to-program data interchange between MISA and the primer designing software Primer3. Primer pairs were designed from the flanking sequences of SSRs using primer3_core (www.broadinstitute.org/genome_software/other/primer3.html) in batch mode via the p3_in.pl and p3_out.pl perl scripts available within the MISA package. The primer designing parameters were: 100–280 bpd amplicon size, 60°C optimal annealing temperature, 20 bp optimal primer length and 50% optimal GC content. Three sets of primer pairs were designed for each SSR to provide alternatives if amplification was unsuccessful. Redundancy of markers was checked by e-PCR by keeping a margin of 50 bp for product size; no mismatches were allowed in primer binding and the word size was kept at 11 bp [41]. The nomenclature of the primer set was B01d12 using prefix B for Brachypodium followed by chromosome number then d for distachyon, followed by SSR number as per their order from north to south on each chromosome.

All SSR markers were grouped into Class I (≥20) and Class II (12–19 bp) types. Graphical presentation of SSR distribution in different chromosomes of Brachypodium was made using MapChart 2.2 software (http://www.biometris.wur.nl/uk/Software/MapChart).

Statistical analysis and functional annotation

All data were tabulated and statistical analysis was performed using SPSS 10 software package (www.spss.com). All possible SSR types were analyzed for their abundance and density per Mb for both nuclear and chloroplast genomes as well as for coding sequences. The protein sequences of genes containing GC-rich repeats were downloaded from the Phytozome database (www.phytozome.net/) and functional annotation of these genes in the Brachypodium genome was performed by using Gene Ontology Tools.

Validation of SSR markers for amplification

A set of 16 Brachypodium genotypes was used for the validation of 55 SSR markers selected from 5 Mb intervals on all the five chromosomes (Table S5). Genomic DNA from the 16 genotypes was isolated from young leaves. The PCR reactions of 10 µl volume containing 20 ng of genomic DNA, 5 pmole each of forward and reverse primers, 0.1 mM dNTPs, 1x PCR buffer (10 mM Tris, pH 8.0, 50 mM KCl and 50 mM ammonium sulphate), 1.8 mM MgCl2, and 0.2 unit of Taq DNA polymerase was performed in a thermal cycler. The cycling conditions involved initial denaturation at 94°C for 4 min, followed by 35 cycles of denaturation at 94°C for 1 min, primer annealing at 55–60°C for 1 min, and primer extension at 72°C for 1 min. A final extension at 72°C for 7 min was done and products stored at 4°C until electrophoresis. The PCR products were resolved by electrophoresis in 3% Agarose gels in 1x TBE buffer and visualized by ethidium bromide staining.

Supporting Information

Figure S1.

Genome-wide distribution of microsatellites in monocot and dicot plant species with respect to motif length. Microsatellites were identified with criteria of mono- to hexa-nucleotides motifs using MISA software tool, and the minimum repeat unit was defined as 10 for mono-, 6 for di-, 5 for tri-, tetra-, penta-, and hexa-nucleotides.

https://doi.org/10.1371/journal.pone.0021298.s001

(DOC)

Figure S2.

Distribution of SSRs motif with mono-, di-, tri-, tetra-, penta-, and hexa-repeats on the chromosome 5 of Brachypodium genome.

https://doi.org/10.1371/journal.pone.0021298.s002

(DOC)

Figure S3.

Brachypodium microsatellite marker (BraMi) database includes details of 27,392 SSR markers which is (a) searchable by marker Id, Location and Key words and (b) provides all the details of SSR marker including primer sequences, annealing temperature and genomic location on pseudomolecule

https://doi.org/10.1371/journal.pone.0021298.s003

(DOC)

Table S1.

Distribution of microsatellite in coding DNA sequence (CDS) of monocot and dicot plant genomes.

https://doi.org/10.1371/journal.pone.0021298.s004

(DOC)

Table S2.

Functional annotation of coding GC-rich repeat containing genes in Brachypodium genome.

https://doi.org/10.1371/journal.pone.0021298.s005

(XLS)

Table S3.

Distribution of microsatellite with different motifs in chloroplast genome of six plant species.

https://doi.org/10.1371/journal.pone.0021298.s006

(DOC)

Table S4.

List of genome wide SSR markers designed for Brachypodium.

https://doi.org/10.1371/journal.pone.0021298.s007

(XLS)

Table S5.

Details of SSR markers chosen at 5 Mb intervals from the Brachypodium genome used to validate for the PCR amplification.

https://doi.org/10.1371/journal.pone.0021298.s008

(DOC)

Author Contributions

Wrote the paper: TRS. Compiled first draft of the MS: HS RKD. Designed the project: TRS. Provided input in designing the project: NKS RNG. Designed the SSR marker development pipeline using Perl Scripts and core primer3 software: HS RKD. Developed local database to categorized markers in genic and intergenic: HS RKD. Conducted comparative analysis, mined SSR, designed primers: HS RKD. Validated markers in wet lab experiments: HS RKD AS. Involved in database development: VPS DKG. Maintained and provided the seeds of all Brachypodium genotypes for validation of markers: JCR.

References

  1. 1. Powell W, Machray G, Provan J (1996) Polymorphism revealed by simple sequence repeats. Trends Plant Sci 1: 215–222.
  2. 2. Gupta PK, Rustogi S, Sharma S, Singh R, Kumar N, et al. (2003) Transferable EST-SSR markers for the study of polymorphism and genetic diversity in bread wheat. Mol Genet Genomics 270: 315–23.
  3. 3. Saha MC, Mian MA, Eujayl I, Zwonitzer JC, Wang L, et al. (2004) Tall fescue EST-SSR markers with transferability across several grass species. Theor Appl Genet 109: 783–91.
  4. 4. Asp T, Frei UK, Didion T, Nielsen KK, Lübberstedt T (2007) Frequency, type, and distribution of EST-SSRs from three genotypes of Lolium perenne, and their conservation across orthologous sequences of Festuca arundinacea, Brachypodium distachyon, and Oryza sativa. BMC Plant Biol 127: 36.
  5. 5. The International Brachypodium Initiative (2010) Genome sequencing and analysis of the model grass Brachypodium distachyon. Nature 463: 763–768.
  6. 6. Tuskan GA, Difazio S, Jansson S, Bohlmann J, Grigoriev I, et al. (2006) The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science 313(5793): 1596–604.
  7. 7. Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, et al. (2009) The Sorghum bicolor genome and the diversification of grasses. Nature 457: 551–556.
  8. 8. International Rice Genome Sequencing Project (2005) The map-based sequence of the rice genome. Nature 436: 793–800.
  9. 9. Arabidopsis Genome Initiative (2000) Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature 408: 796–815.
  10. 10. McCouch SR, Teytelman L, Xu Y, et al. (2002) Development and mapping of 2240 new SSR markers for rice (Oryza sativa L). DNA Res 9: 199–207.
  11. 11. Parida SK, Rajkumar KA, Dalal V, Singh NK, Mohapatra T (2006) Unigene derived microsatellite markers for the cereal genomes. Theor Appl Genet 112: 808–817.
  12. 12. Singh H, Deshmukh RK, Singh A, Singh AK, Gaikwad K, et al. (2010) Highly variable SSR markers suitable for rice genotyping using agarose gels. Mol Breeding 25(2): 359–364.
  13. 13. Draper J, Mur LAJ, Jenkins G, Ghosh-Biswas GC, Bablak P, et al. (2001) Brachypodium distachyon. A new model system for functional genomics in grasses. Plant Physiology 127(4): 1539–1555.
  14. 14. Garvin DF, Gu YQ, Hasterok R, Hazen SP, Jenkins G, et al. (2008) Development of Genetic and Genomic Research Resources for Brachypodium distachyon, a New Model System for Grass Crop Research. Crop Science 48: S-69–S-84.
  15. 15. Garvin DF, Mckenzie N, Vogel JP, Mockler TC, Blankenheim ZJ, et al. (2010) : An SSR-based genetic linkage map of the model grass Brachypodium distachyon. Genome 53(1): 1–13.
  16. 16. Levinson G, Gutman , GA (1987) Slipped-strand mispairing: a major mechanism for DNA sequence evolution. Mol Biol Evol 4: 203–221.
  17. 17. Schlötterer C, Tautz D (1992) Slippage synthesis of simple sequence DNA. Nucleic Acids Res 20: 211–215.
  18. 18. Hancock JM (1999) Microsatellites and other simple sequences: genomic context and mutational mechanisms. 1-9. New York: Oxford University Press. In Microsatellites: Evolution and applications Ed Goldstein D, Schlötterer C,.
  19. 19. Chakraborty R, Kimmel M, Stivers DN, Davison LJ, Deka R (1997) Relative mutation rates at di-, tri-, and tetranucleotide microsatellite loci. Proc Natl Acad Sci U S A 94: 1041–1046.
  20. 20. Anderson TJ, Su XZ, Roddam A, Day KP (2000) Complex mutations in a high proportion of microsatellite loci from the protozoan parasite Plasmodium falciparum. Mol Ecol 9: 1599–1608.
  21. 21. Whittaker JC, Harbord RM, Boxall N, Mackay I, Dawson G, et al. (2003) Likelihood-based estimation of microsatellite mutation rates. Genetics 164: 781–787.
  22. 22. Modrich P, Lahue R (1996) Mismatch repair in replication fidelity, genetic recombination, and cancer biology. Ann Rev Biochem 65: 101–133.
  23. 23. Harr B, Todorova J, Schlötterer C (2002) Mismatch repair-driven mutational bias in D melanogaster. Mol Cell 10: 199–205.
  24. 24. Riley DE, Jeon JS, Krieger JN (2007) Simple repeat evolution includes dramatic primary sequence changes that conserve folding potential. Biochem Biophys Res Commun 355: 619–625.
  25. 25. Wolfe KH, Gouy M, Yang YW, Sharp PM, Li WH (1989) Date of the monocot-dicot divergence estimated from chloroplast DNA sequence data. Proc Natl Acad Sci 86(16): U S A: 6201–6205.
  26. 26. Yu J, Hu S, Wang J, Wong GKS, Li S, et al. (2002) A draft sequence of the rice genome (Oryza sativa L. ssp. indica). Science 296: 79–92.
  27. 27. Vision TJ, Brown DG, Tanksley SD (2000) The origins of genomic duplications in Arabidopsis. Science 290: 2114–2117.
  28. 28. Blanc G, Hokamp K, Wolfe KH (2003) A recent polyploidy superimposed on older large-scale duplications in the Arabidopsis genome. Genome Res 13: 137–144.
  29. 29. Goff SA, Ricke D, Lan TH, Presting G, Wang RL, et al. (2002) A draft sequence of the rice genome (Oryza sativa L. ssp japonica). Science 296: 92–100.
  30. 30. Li S, Yin T, Wang M, Tuskan GA (2011) Characterization of microsatellites in the coding regions of the Populus genome Mol Breeding 27: 59–66.
  31. 31. Sonah H, Deshmukh RK, Singh VP, Gupta DK, Singh NK, et al. (2010) Genomic resources in horticultural crops: Status, Utility and Challenges. Biotechnological Advances 29: 199–209.
  32. 32. Weber JL, Wong C (1993) Mutation of human short tandem repeats. Hum Mol Genet 2: 1123–1128.
  33. 33. Xu X, Peng M, Fang Z, Xu X (2000) The direction of microsatellite mutations is dependent upon allele length. Nat Genet 24: 396–399.
  34. 34. Temnykh S, DeClerck G, Lukashova A, Lipovich L, Cartinhour S, et al. (2001) Computational and experimental analysis of microsatellites in rice (Oryza sativa L.): frequency, length variation, transposon associations, and genetic marker potential. Genome Res 11(8): 1441–52.
  35. 35. Vogel J, Hill T (2008) High-efficiency Agrobacterium-mediated transformation of Brachypodium distachyon inbred line Bd21-3. Plant Cell Rep 27: 471–478.
  36. 36. Alves SC, Worland B, Thole V, Snape JW, Bevan MW, et al. (2009) A protocol for Agrobacterium-mediated transformation of Brachypodium distachyon community standard line Bd21. Nature Protocols 4: 638–649.
  37. 37. Vogel J, Tuna M, Budak H, Huo N, Gu Y, et al. (2009) Development of SSR markers and analysis of diversity in Turkish populations of Brachypodium distachyon. BMC Plant Biol 9: 88.
  38. 38. Parida SK, Yadava DK, Mohapatra T (2010) Microsatellites in Brassica unigenes: relative abundance, marker design, and use in comparative physical mapping and genome analysis. Genome 53: 55–67.
  39. 39. Cho YG, Ishii T, Temnykh S, Chen X, Lipovich L, et al. (2000) Diversity of microsatellites derived from genomic libraries and GenBank sequences in rice (Oryza sativa L). Theor Appl Genet 100: 713–722.
  40. 40. Chabane K, Ablett GA, Cordeiro GM, Valkoun J, Henry RJ (2005) EST versus genomic derived microsatellite markers for genotyping wild and cultivated barley. Genet Resour Crop Evol 52: 903–909.
  41. 41. Schuler GD (1997) Sequence Mapping by Electronic PCR. Genome Res 7: 541–550.