Skip to main content

Distinct patterns of SSR distribution in the Arabidopsis thalianaand rice genomes

Abstract

Background

Simple sequence repeats (SSRs) in DNA have been traditionally thought of as functionally unimportant and have been studied mainly as genetic markers. A recent handful of studies have shown, however, that SSRs in different positions of a gene can play important roles in determining protein function, genetic development, and regulation of gene expression. We have performed a detailed comparative study of the distribution of SSRs in the sequenced genomes of Arabidopsis thaliana and rice.

Results

SSRs in different genic regions - 5'untranslated region (UTR), 3'UTR, exon, and intron - show distinct patterns of distribution both within and between the two genomes. Especially notable is the much higher density of SSRs in 5'UTRs compared to the other regions and a strong affinity towards trinucleotide repeats in these regions for both rice and Arabidopsis. On a genomic level, mononucleotide repeats are the most prevalent type of SSRs in Arabidopsis and trinucleotide repeats are the most prevalent type in rice. Both plants have the same most common mononucleotide (A/T) and dinucleotide (AT and AG) repeats, but have little in common for the other types of repeats.

Conclusion

Our work provides insight into the evolution and distribution of SSRs in the two sequenced model plant genomes of monocots and dicots. Our analyses reveal that the distributions of SSRs appear highly non-random and vary a great deal in different regions of the genes in the genomes.

Background

Simple sequence repeats (SSRs) are tandem repeat nucleotides (oftentimes defined as being between 1 and 6 base-pairs (bp)) in DNA sequences. They can be found in any genome (both eukaryote and prokaryote) and in any region (protein coding regions and non-coding regions). Historically, SSRs were used often as genetic markers, helping to classify and identify various species. However, recent research has shown that SSRs have many important functions in terms of development, gene regulation, and evolution.

The locations of SSRs appear to determine the types of functional role SSRs might play, and changes in SSRs in different genetic locations can lead to changes in the phenotypes of an organism [1]. SSRs in coding regions can determine whether or not a gene gets activated or whether the protein product is truncated [1]. For instance, expansion of CAG repeats in the coding region of HD genes in humans can lead to Huntington's disease, most likely through activation of so-called 'toxic' proteins. The development of the nervous system in Drosophila appears to be associated with length variation of trinucleotide repeats in genes involved in developmental control [2]. Most recently, Fondon and Garner [3] have shown that the fast morphological evolution in domesticated dogs is due to the contraction/expansion of SSRs in the coding regions of the Alx-4 and Runx-2 genes.

SSRs in other genic regions can have large effects on organisms as well. For example, SSRs in 5'-untranslated regions (UTRs) have an effect on gene transcription and/or regulation [1]. The human calmodulin-1 (hCALM1) gene has a CAG repeat in a 5'UTR that when deleted causes a decrease in expression by 45% [4]. Intron SSRs can affect gene transcription, regulation, mRNA splicing, and gene silencing [1]. For example, the first intron of the gene encoding tyrosine hydroxylase contains a TCAT repeat that acts as a transcription regulatory element [5]. SSRs found in 3'-UTRs are involved in gene silencing and transcription slippage [1], as in the case of a CTG expansion in a kinase gene that causes myotonic dystrophy type 1 through transcription slippage [6].

The functional study of SSRs has been largely restricted to animals. In plants, the majority of research used SSRs as genetic markers to study populations and genetic diversity [7–9] and to determine sex in dioecious plants [10]. Several studies have been done to characterize the distribution of SSRs in Arabidopsis. For example, Casacuberta et al. [11] examined the abundant types of mono- and dinucleotide repeats in coding sequences of the unfinished Arabidopsis genome. Zhang et al. [12] did a more comprehensive survey of SSRs in Arabidopsis and showed that SSRs in general were more favored in upstream regions of genes and that trinucleotide repeats were the most common repeats found in the coding regions. The purpose of this paper is to compare SSRs between the two plant species: Arabidopsis thaliana and rice (Oryza sativa). These two plants have their entire genomes largely sequenced, so in-depth comparisons of SSRs can be made for not only their entire genome, but specific regions as well, such as exons, introns, and UTRs.

Results

The coding regions (exons) of 26,416 genes in Arabidopsis thaliana and 57,915 genes in rice were analyzed. The 57,915 rice genes include 14,273 transposable element (TE)-related genes. Excluding these genes from analyses does not change our results qualitatively, and, therefore, we report here only the results for the 57,915 genes. The 5'UTR regions of 16,355 genes and the 3'UTR regions of 17,617 genes were used for Arabidopsis. For rice these values were 12,907 and 14,839, respectively. The lower number of UTR sequences than coding regions is due to the fact that UTRs are curated only when there is full-length cDNA or expressed sequence tag (EST) evidence supporting the annotation [13]. Therefore, although the number of UTRs is reduced, we have high data quality. Finally, the intron data of 21,157 genes in Arabidopsis and 45,633 genes in rice were analyzed as well. The total lengths of these regions are shown in Table 1. Detailed in the following sections are the SSR amounts for the various regions. Tables 2 and 3 list the most common repeat types, including all circular permutations and their complements, similar to previous analyses done on this topic [12, 14].

Table 1 Total lengths of the studied regions and the amounts of SSRs therein
Table 2 Densities of the most abundant SSR (mono- and dinucleotide) types in different regions
Table 3 Densities of the most abundant SSR (tri- and tetranucleotide) types in different regions

Whole genome SSRs

The Arabidopsis genome contains a total of 104,102 SSRs (Table 1). The genome average SSR density is thus approximately 875 per mega-base (MB). SSRs with periods of 1 to 10 (mono-, di-, tri-, and so on) account for 33.9% (35,256 of 104,102), 12.8% (13,295), 17.5% (18,244), 5.5% (5,731), 7.8% (8,097), 8.9% (9,261), 5.2% (5,392), 3.5% (3,612), 3.6% (3,720), and 1.4% (1,494), respectively.

In comparison, the rice genome contains a total of 298,819 SSRs (Table 1). The genome average SSR density is approximately 807/MB. SSRs with periods of 1 to 10 account for 18.3% (54,809), 14.7% (43,949), 23.9% (71,373), 7.9% (23,756), 8.3% (24,718), 11.9% (35,602), 4.4% (13,216), 3.9% (11,628), 4.6% (13,854), and 2% (5,914), respectively. Figure 1 shows the corresponding SSR densities.

Figure 1
figure 1

Comparison of whole genome SSR densities. A comparison of the SSR densities (for SSRs of period 1 to 10) in the whole genome of Arabidopsis and rice.

Exon SSRs

The exon regions in Arabidopsis contain a total of 12,168 SSRs (Table 1). The average SSR density is thus approximately 334/MB. SSRs with periods 1 to 10 account for 1% (121), 2.2% (264), 65.4% (7,961), 1% (121), 2.4% (296), 17.3% (2,102), 1.6% (193), 1.6% (192), 6.9% (844), and 0.6% (74), respectively.

In comparison, the rice exon regions contain a total of 55,338 SSRs (Table 1). The average SSR density is approximately 658/MB. SSRs with periods of 1 to 10 account for 0.6% (354), 2.4% (1,353), 64% (35,437), 1.2% (637), 2.5% (1,409), 18.6% (10,268), 1.5% (856), 1.7% (929), 6.9% (3,791), and 0.5% (304), respectively. Figure 2 shows the corresponding SSR densities.

Figure 2
figure 2

Comparison of exon SSR densities. A comparison of the SSR densities (for SSRs of period 1 to 10) in the coding (exon) regions of Arabidopsis and rice.

Intron SSRs

The intron regions in Arabidopsis contain a total of 17,756 SSRs (Table 1). The average SSR density is approximately 815/MB. SSRs with periods of 1 to 10 account for 39.5% (7,011), 15.1% (2,674), 11.2% (1,981), 7.2% (1,280), 8.1% (1,432), 7.9% (1,410), 4.5% (805), 3% (538), 2.4% (421), and 1.1% (204), respectively.

In comparison, the rice intron regions contain a total of 87,529 SSRs (Table 1). The average SSR density is approximately 634/MB. SSRs with periods of 1 to 10 account for 21% (18,360), 11.2% (9,794), 28.8% (25,169), 7.3% (6,379), 7% (6,160), 11.9% (10,410), 3.8% (3,297), 3.2% (2,805), 4.5% (3,932), and 1.4% (1,223), respectively. Figure 3 shows the corresponding SSR densities.

Figure 3
figure 3

Comparison of intron SSR densities. A comparison of the SSR densities (for SSRs of period 1 to 10) in the intron regions of Arabidopsis and rice.

5'UTR SSRs

The 5'UTR regions in Arabidopsis contain a total of 6,146 SSRs (Table 1). The average SSR density is approximately 2,364/MB. SSRs with periods of 1 to 10 account for 18.4% (1,130), 28.1% (1,727), 29.7% (1,827), 4.5% (275), 5% (310), 6.8% (417), 2.6% (160), 2.2% (137), 1.8% (110), and 0.9% (53), respectively.

In comparison, the rice 5'UTR regions contain a total of 12,310 SSRs (Table 1). The average SSR density is approximately 3971/MB. SSRs with periods of 1 to 10 account for 6.2% (769), 12.6% (1,552), 42.1% (5,179), 8.5% (1,046), 12.8% (1,575), 11.1% (1,367), 2.4% (297), 1.8% (222), 1.7% (209), 0.8% (94), respectively. Figure 4 shows the corresponding SSR densities.

Figure 4
figure 4

Comparison of 5'-UTR SSR densities. A comparison of the SSR densities (for SSRs of period 1 to 10) in the 5'-UTR regions of Arabidopsis and rice.

3'UTR SSRs

The 3'UTR regions in Arabidopsis contain a total of 4,910 SSRs (Table 1). The average SSR density is approximately 982/MB. SSRs with periods of 1 to 10 account for 34.6% (1,700), 12.5% (615), 17.5% (858), 8.8% (434), 8.4% (411), 7.8% (383), 4.2% (208), 3.3% (160), 2.2% (107), 0.7% (34), respectively.

In comparison, the rice 3'UTR regions contain a total of 5,658 SSRs (Table 1). The average SSR density is approximately 832/MB. SSRs with periods of 1 to 10 account for 26.8% (1,515), 10.6% (602), 18.4% (1,042), 14.6% (827), 8.9% (502), 8.3% (472), 4.8% (270), 3.6% (206), 2.9% (162), and 1.1% (60), respectively. Figure 5 shows the corresponding SSR densities.

Figure 5
figure 5

Comparison of 3'UTR SSR densities. A comparison of the SSR densities (for SSRs of period 1 to 10) in the 3'-UTR regions of Arabidopsis and rice.

Observed versus expected densities of SSRs in different regions

The expected numbers of mono-, di-, tri-, and tetranucleotide SSRs were calculated using the de Wachter formula, as detailed in the methods section. Table 4 shows the observed and expected densities for exon, 5'UTR, 3'-UTR, and intron regions. For both Arabidopsis and rice, the 5'UTR, 3'UTR, and intron regions show similar patterns: the observed densities of mono-, di-, tri-, and tetranucleotide SSRs are much higher than those expected. In contrast, for the exon regions, in Arabidopsis, only trinucleotide SSRs are more abundant than expected, all other types of SSRs (mono-, di-, and tetranucleotides) being present at a much lower frequency than expected; in rice, the observed densities of all types of SSRs except tetranucleotide SSRs are higher than those expected.

Table 4 Observed and expected densities of SSRs in different genic regions

GO categories of genes with most repeats

In Arabidopsis, the average amount of repeats per gene is approximately two. The ten most common Gene Ontology (GO) categories for the genes with high SSR densities are: chloroplast (GO ID: 0009507), nucleus (GO ID: 0005634), mitochondria (GO ID: 0005739), extracellular region (GO ID: 0005576), transcription factor activity (GO ID: 0003700), nucleotide binding (GO ID: 0005524), DNA binding (GO ID: 0003677), structural constituent of cell wall (GO ID: 0005199), cell wall (GO ID: 0005618), and hydrolase activity (GO ID: 0004553). The hypergeometric tests show that there is no statistically significant enrichment of SSRs in genes belonging to these GO categories (P > 0.05).

In rice, the average amount of repeats per gene is approximately three. The 10 most common GO categories for the genes with high SSR densities are: transferase activity (GO ID: 0016740), hydrolase activity (GO ID: 0016787), catalytic activity (GO ID: 0003824), response to stress (GO ID: 0006950), membrane (GO ID: 0016020), protein binding (GO ID: 0005515), binding (GO ID: 0005488), DNA binding (GO ID: 0003677), response to biotic stimulus (GO ID: 0009607), and kinase activity (GO ID: 0016301). Among these GO categories, the hypergeometric tests show that SSRs are significantly enriched in genes with GO categories of DNA binding (P = 1.93e-71), response to stress (P = 2.2e-48), and binding (P = 4.47e-46).

Amino acid runs in coding regions

In coding regions, trinucleotide repeats are in fact amino acid runs. Because each amino acid is encoded by one or more synonymous codon, we are interested in how trinucleotide SSRs have contributed to the single amino acid runs. Specifically, we denote the amino acid runs as 'homogeneous runs' if they are trinucleotide repeats, in which case only one codon is used for the amino acid runs. We created a perl script that we used to analyze the proteomes of Arabidopsis and rice and calculated the amounts of amino acid runs of length ≥5 [15, 16].

A total of 7,258 amino acid runs (we require at least 5 of the same amino acid in a row) were found in the 26,416 protein sequences in Arabidopsis. The five most frequent types of amino acid run are (Table 5): serine with 1,997 runs (27.5%); proline with 865 runs (11.9%); glycine with 853 runs (11.8%); glutamic acid with 831 runs (11.4%); and glutamine with 451 runs (6.2%).

Table 5 Numbers of amino acid runs and homogeneous runs

A total of 28,367 amino acid runs were found in the 57,915 protein sequences in rice. The five most frequent types of amino acid run are (Table 5): alanine with 7,477 runs (26.4%); glycine with 6,349 runs (22.4%); proline with 3,727 runs (13.1%); serine with 2,862 runs (10.1%); and arginine with 1,636 (5.8%).

As expected, for both Arabidopsis and rice, the proportion of homogeneous runs decreases as the number of synonymous codons increases. Interestingly, we found that the proportions of homogeneous runs in rice are always much higher than that in Arabidopsis for all amino acids except aspartic acid and asparagine (Table 5).

The difference in the distribution of amino acid runs could be due to the fact that different amounts of genes from rice and Arabidopsis were analyzed. To examine this issue, we analyzed the amino acid runs for only the orthologous genes between rice and Arabidopsis. The orthologous genes were downloaded from the Gramene website [17]. The complete list of genes is included in Additional data file 1. Altogether we analyzed 10,519 pairs of orthologous genes and found that the results yielded similar values, with the most frequent types of amino acid runs remaining the same and the proportions staying consistent with the results for all genes (Table 6).

Table 6 Distribution of amino acid runs of only the orthologous genes between Arabidopsis and rice, in comparison with the total genes in the two species

Discussion

Monocots and dicots are thought to have diverged from a common species approximately 200 million years ago [18]. Arabidopsis and rice are the representative species in their respective groups whose genomes, because of their small sizes, have been largely sequenced. Arabidopsis has been traditionally used as a model plant species, and rice has gathered much attention due to its significance in being one of the major food resources in the world.

Comparative analyses of the Arabidopsis and rice genomes have yielded a number of insights about the two plants. First, since Arabidopsis and rice have 5 and 12 chromosomes, respectively, it has been commonly thought that monocots have undergone genome duplication after the split of the monocot and dicot species [18]. However, studies of both the Arabidopsis and rice genomes suggest a different story: Arabidopsis, despite having the smallest genome among the dicots, might have gone through several rounds of genome duplication [19–21]. In contrast, the rice genome shows no distinct pattern of genome duplication, instead appearing to be more the product of gradual small scale duplications and loss of duplicated genes [22–24].

Second, the Arabidopsis genome is compact, with an average gene size of approximately 2.4 kb [22]. In contrast, the rice genes are on average four times larger (approximately 9.9 kb) [25]. The much larger average gene size in rice seems to be due to the larger introns [22, 25]. Third, the gene sets in the Arabidopsis and rice genomes appear highly asymmetric to each other: approximately 80% to 90% of the Arabidopsis genes have rice homologs, yet only 49.4% to 71% of the rice genes have Arabidopsis homologs [22, 23, 25]; therefore, many genes in the rice genome might be monocot specific. In fact, these genes also do not have homologs in other sequenced genomes, including Drosophila melanogaster, Caenorhabditis elegans, Saccharomyces cerevisiae, and Schizosaccharomyces pombe [22, 23]. Unfortunately, most of these annotated rice genes have no known functions. The question remains as to how so many rice or monocot specific genes have come into existence and what their functional and evolutionary significance is.

Fourth, the G+C content of the Arabidopsis genes is rather homogenous; in contrast, the G+C content of the rice genes decreases from the 5'UTR to the 3'UTR by several percent to approximately 25% [22, 26]. This gradient of G+C content in rice genes is still present when comparing rice genes with the corresponding orthologs in Arabidopsis [22, 26]. Here, we further examined the nucleotide components in different regions of the Arabidopsis and rice genes. The genome average G+C content is 36% in Arabidopsis and 43.6% in rice, and the higher average G+C content in rice than in Arabidopsis is consistently observed for all other genic regions (5'UTR, 3'UTR, introns, and exons). In rice, the average G+C content ranking is 5'UTR (55.7%) > exons (53.2%) > introns (43.8%) > 3'UTR (40.2%). In Arabidopsis, the G+C content ranking is exons (44.2%) > 5'UTR (38.3%) > 3'UTR (33.8%) > introns (32.5%).

The gradient of G+C content along genes has also been observed in other monocots [26]. This suggests two likely evolutionary scenarios. One is that the ancestral species of monocots and dicots had no G+C gradient along genes and the genome wide G+C incline in monocots formed at early stages after the divergence of monocots and dicots, since otherwise we would have to assume that many monocot species evolved this trait independently. The second scenario is that the ancestral species of monocots and dicots possessed this trait and the dicot species has lost this trait subsequently. One can examine this issue using a proper outgroup species. It remains an open question as to what evolutionary transitions correspond to this genomic trait, if it was acquired in monocots, and the significance of having distinct G+C content in different genic regions.

Fifth, the evolution of SSRs, which is examined in this study in detail, shows several similarities and differences between Arabidopsis and rice. In the following, we discuss the similarities and differences both within and between the two genomes.

The results on the SSR distribution show that for both species, the majority of the SSRs are mono-, di-, tri-, tetra-, and pentanucleotides, accounting for up to approximately 80% of all the SSRs found in various regions and the genomes (Figure 6 and 7). For both species, the distribution of SSRs in the 5'-UTRs and exons show patterns distinct from the other genic regions and the entire genomes. Introns and 3'-UTRs have a similar SSR distribution to the whole genome for SSRs with periods of 1 to 10. The discrepancies between Arabidopsis and rice SSR distribution are most pronounced for SSRs with periods of 1 to 4 (Figures 1 to 5).

Figure 6
figure 6

Arabidopsis cumulative SSR distribution. Graph showing the cumulative distribution of SSR percentages for periods 1 to 10 in Arabidopsis.

Figure 7
figure 7

Rice cumulative SSR distribution. Graph showing the cumulative distribution of SSR percentages for periods 1 to 10 in rice.

Remarkably, the average SSR density (measured by the count of SSRs/MB) among all regions for both Arabidopsis and rice is highest for the 5'UTRs (Figure 8). A comparison between Arabidopsis and rice shows that the average repeat densities in the two genomes is similar for introns, the 3'UTR, and the whole genome, but not for exons and the 5'UTR. The SSR densities in the 5'UTR and exon regions show an almost twofold difference between the two plants, which is the combined result of much higher densities for SSRs of period 3 to 6 in the 5'UTR and much higher densities for SSRs of period 3 and period 6 in the exon regions in rice than in Arabidopsis (Figures 2 and 4). Taken together, the 5'UTR and exons stand out as the regions that differ from the remaining regions, a pattern unnoticed before, and deserve further studies on the role of SSRs in them.

Figure 8
figure 8

Comparison of SSR densities in different regions. A comparison of the SSR densities across various genic regions.

In most regions and when doing a comparison of the whole genomes, Arabidopsis shows a great affinity for mononucleotide repeats compared to rice (Figures 1 to 5). The comparison between the whole genomes of rice and Arabidopsis clearly shows that Arabidopsis has a higher percentage of mononucleotide repeats (33.9%) than rice (18.3%). The only regions where mononucleotide repeats are not as high are the 5'UTR and coding regions. Rice seems to have a greater affinity for trinucleotide repeats, with an exceptionally high density of approximately 1,670/MB in the 5'UTR, even higher than in exon regions (approximately 421.4/MB). In fact, trinucleotide SSRs are the major type of SSR in all regions except 3'UTRs (Figures 1 to 5).

Yet despite these differences in which type of SSR is most common for each organism, rice and Arabidopsis show similar distribution of the SSR types (period) in their coding regions. Both have a majority of trinucleotide repeats (Arabidopsis 65.4%, rice 64%) and a lack of any other types other than those divisible by three (the latter account for only 10.4% in Arabidopsis, and 10.4% in rice). The high percentage of SSRs with period divisible by three is expected because of the nature of translation and how it relies on triplet codons. It also corresponds with previous research that has shown that tri- and hexanucleotide repeats are the most common in the coding regions of eukaryotes [27, 28]. The rather similar distribution of SSRs of period 1 to 10 in the coding regions of the two genomes could be partially explained by the observation that approximately 80% to 90% of the predicted Arabidopsis proteins show homology with the predicted rice proteins [22, 23, 25].

Apart from the similar distribution of the SSRs in coding regions, the two plants show at least two major differences. First, the densities for the tri- and hexanucleotide SSRs are much higher in rice (trinucleotide, approximately 421.4/MB; hexanucleotide, approximately 122.1/MB) than in Arabidopsis (trinucleotide, approximately 218.7/MB; hexanucleotide, 57.7/MB). Second, when comparing the amino acid runs, while both plants contain many runs of glycine and proline (accounting for either the second or third highest amounts), they differ in what amino acid occurs in the highest amount of runs. Serine is in the highest amount for Arabidopsis and alanine is in the highest amount for rice. However, rice still has a large amount of serine repeats (the fourth largest amount of runs when comparing all of the amino acid runs), while Arabidopsis has very few alanine runs. Our initial hypothesis was that this seems to be consistent with the observation that Arabidopsis shows homology to rice but rice does not show as much homology to Arabidopsis [22, 26]. However, we observed the same patterns when we limited our analysis to only orthologous genes between rice and Arabidopsis, suggesting that the difference in coding regions between rice and Arabidopsis in amino acid runs is not the result of differences in gene content.

Over its long history, Arabidopsis has undergone at least three polyploidy events [19], leaving it with many duplicates throughout its genome. In fact, over 37% of genes are part of gene families that contain more than five members. Genes that have duplicates more often than chance are involved in signal transduction and transcription, especially in the nucleus and plasma membrane [29]. According to the data we have analyzed, these types of genes also contain an abundance of SSRs.

In terms of SSR types, a few trends can be observed. Both the rice and Arabidopsis genomes have A/T as the most frequent mononucleotide repeats, and the most common dinucleotide repeats are similar between the two genomes as well. However, the most common tri- and tetranucleotide SSR types are mostly different between the two species. What is observable in Arabidopsis is that most of the larger SSRs (≥3) consist of long strings of either A or T broken by an additional nucleotide (such as AAAAG or TTTTTC). This shows again a tendency towards mononucleotide repeats. In rice, there is a trend in the trinucleotide repeats which, with little exception, consist of various combinations of C and G. Common SSR types consist of two Cs and one G or two Cs and one T in various combinations (Table 3).

Conclusion

The SSRs show distinct patterns of distribution among different regions of the genes in both Arabidopsis and rice genomes. The amounts of differences in the SSRs between the two genomes are the combined results of ancient species divergences and the individual evolution of these plants. Considering the big discrepancy in gene content between the two plants [22, 23], we found the differences between SSRs in Arabidopsis and rice less surprising, because of their much higher mutations rates than regular genes [1]. However, the potential functional significance of the SSR changes is an important issue that is yet to be determined.

Materials and methods

We downloaded the data for Arabidopsis from the TAIR website [30] and the data for rice from the TIGR website [31]. For both species, the sequences have already been curated based on their genic locations, including the Arabidopsis and rice complete genomic sequences, coding regions (exons of the genes), introns, 5'UTR, 3'UTR, and protein sequences.

We applied mreps to search for SSRs [32]; mreps is a tool that was specifically developed to identify repeats in DNA sequences. The algorithm consists of a combinatorial and a heuristic treatment to determine SSRs. In the combinatorial step, the maximum runs of tandem repeats are found within the given error threshold. Then in the heuristic treatment, the best candidate SSR is determined for each run and overlapping repeats are merged. The results are also filtered to account for statistically expected results. Afterwards, these results are gathered and iterated for all resolution values until the final resolution value is reached.

We considered only the perfect SSRs with a length of longer than 10 bp. Throughout the paper, we have used the convention of mreps and refer to the size of a repeat unit as 'period'. For example, mononucleotide SSRs are SSRs of period 1. Because our analyses revealed that simple repeats with periods greater than 10 are rare, we focused on SSRs with periods 1 to 10. Note that a common and arbitrary definition of SSRs is simple repeats with periods of 1 to 6.

Using perl scripts, we sorted the mreps data into files where each SSR was organized by the locus in which it was contained. To examine how the observed numbers of SSRs compared to the expected numbers of SSRs in different genic regions, we calculated the expected number of SSRs using the following formula [33]:

N(M t ) = p(M)t[1 - p(M)][N'(1 - p(M) + 2L]

N' = N - tL - 2L + 1

In this formula, M is the repeat unit (repeat type), N(M t ) is the expected number of times that, in a DNA segment of length N, we find t consecutive Ms, L is the length of M, and p(M) is the probability of M (obtained by multiplying the probability of each nucleotide contained within the repeat unit M).

To examine whether SSRs are associated with gene function, we grouped all genes based on their GO [34] categories. The GO data show in what category of protein the gene product of each gene falls. Each gene can belong to multiple categories and these categories are organized into three main categories: molecular function, biological process, and cellular component. Every gene has at least one subcategory from these three main categories assigned to it, with some having additional subcategories as well. For Arabidopsis, we downloaded a complete listing of the GO data from the TAIR website. For rice, we downloaded a complete listing of the GO data from the TIGR website. We applied hypergeometric tests to formally examine whether any particular gene functions (GO categories) have statistically significant SSR enrichment [35]. Suppose that we have a total of n genes, among which there are m genes that have high SSR densities. Suppose further that there are r genes that are in a given GO category, of which k genes are in the high-SSR-density class. The following hypergeometric test gives the significance of enrichment of SSRs for this specific GO category:

Additional data files

The following additional data are available with the online version of this paper. Additional data file 1 contains a list of the orthologous genes that were analyzed.

References

  1. Li YC, Korol AB, Fahima T, Nevo E: Microsatellites within genes: Structure, function, and evolution. Mol Biol Evol. 2004, 21: 991-1007. 10.1093/molbev/msh073.

    Article  PubMed  CAS  Google Scholar 

  2. Karlin S, Burge C: Trinucleotide repeats and long homopeptides in genes and proteins associated with nervous system disease and development. Proc Natl Acad Sci USA. 1996, 93: 1560-1565. 10.1073/pnas.93.4.1560.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  3. Fondon JW, Garner HR: Molecular origins of rapid and continuous morphological evolution. Proc Natl Acad Sci USA. 2004, 101: 18058-18063. 10.1073/pnas.0408118101.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  4. Toutenhoofd SL, Garcia F, Zacharias DA, Wilson RA, Strehler EE: Minimum CAG repeat in the human calmodulin-1 gene 5' untranslated region is required for full expression. Biochim Biophys Acta. 1998, 1398: 315-320. 10.1016/S0167-4781(98)00056-6.

    Article  PubMed  CAS  Google Scholar 

  5. Meloni R, Albanese V, Ravassard P, Treilhou F, Mallet J: A tetranucleotide polymorphic microsatellite, located in the first intron of the tyrosine hydroxylase gene, acts as a transcription regulatory element in vitro. Hum Mol Genet. 1998, 7: 423-428. 10.1093/hmg/7.3.423.

    Article  PubMed  CAS  Google Scholar 

  6. Ranum LPW, Day JW: Dominantly inherited, non-coding microsatellite expansion disorders. Curr Opin Genet Dev. 2002, 12: 266-271. 10.1016/S0959-437X(02)00297-6.

    Article  PubMed  CAS  Google Scholar 

  7. Portis E, Acquadro A, Comino C, Mauromicale G, Saba E, Lanteri S: Genetic structure of island populations of wild cardoon [Cynara cardunculus L. var. sylvestris (Lamk) Fiori] detected by AFLPs and SSRs. Plant Sci. 2005, 169: 199-210. 10.1016/j.plantsci.2005.03.014.

    Article  CAS  Google Scholar 

  8. Lu H, Redus MA, Coburn JR, Rutger JN, McCouch SR, Tai TH: Population structure and breeding patterns of 145 US rice cultivars based on SSR marker analysis. Crop Sci. 2005, 45: 66-76.

    Article  CAS  Google Scholar 

  9. Saini N, Jain N, Jain S, Jain RK: Assessment of genetic diversity within and among Basmati and non-Basmati rice varieties using AFLP, ISSR and SSR markers. Euphytica. 2004, 140: 133-146. 10.1007/s10681-004-2510-y.

    Article  CAS  Google Scholar 

  10. Rode J, In-Chol K, Saal B, Flachowsky H, Kriese U, Weber WE: Sex-linked SSR markers in hemp. Plant Breeding. 2005, 124: 167-170. 10.1111/j.1439-0523.2005.01079.x.

    Article  CAS  Google Scholar 

  11. Casacuberta E, Puigdomenech P, Monfort A: Distribution of microsatellites in relation to coding sequences within the Arabidopsis thaliana genome. Plant Sci. 2000, 157: 97-104. 10.1016/S0168-9452(00)00271-5.

    Article  PubMed  CAS  Google Scholar 

  12. Zhang LD, Yuan DJ, Yu SW, Li ZG, Cao YF, Miao ZQ, Qian HM, Tang KX: Preference of simple sequence repeats in coding and non-coding regions of Arabidopsis thaliana. Bioinformatics. 2004, 20: 1081-1086. 10.1093/bioinformatics/bth043.

    Article  PubMed  CAS  Google Scholar 

  13. Yuan QP, Ouyang S, Liu J, Suh B, Cheung F, Sultana R, Lee D, Quackenbush J, Buell CR: The TIGR rice genome annotation resource: annotating the rice genome and creating resources for plant biologists. Nucleic Acids Res. 2003, 31: 229-233. 10.1093/nar/gkg059.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  14. Jurka J, Pethiyagoda C: Simple repetitive DNA-sequences from primates - compilation and analysis. J Mol Evol. 1995, 40: 120-126. 10.1007/BF00167107.

    Article  PubMed  CAS  Google Scholar 

  15. Faux NG, Bottomley SP, Lesk AM, Irving JA, Morrison JR, de la Banda MC, Whisstock JC: Functional insights from the distribution and role of homopeptide repeat-containing proteins. Genome Res. 2005, 15: 537-551. 10.1101/gr.3096505.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  16. Fiebig A, Kimport R, Preuss D: Comparisons of pollen coat genes across Brassicaceae species reveal rapid evolution by repeat expansion and diversification. Proc Natl Acad Sci USA. 2004, 101: 3286-3291. 10.1073/pnas.0305448101.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  17. Gramene. [http://www.gramene.org]

  18. Wolfe KH, Gouy ML, Yang YW, Sharp PM, Li WH: Date of the monocot dicot divergence estimated from chloroplast DNA-sequence data. Proc Natl Acad Sci USA. 1989, 86: 6201-6205.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  19. Vision TJ, Brown DG, Tanksley SD: The origins of genomic duplications in Arabidopsis. Science. 2000, 290: 2114-2117. 10.1126/science.290.5499.2114.

    Article  PubMed  CAS  Google Scholar 

  20. Blanc G, Hokamp K, Wolfe KH: A recent polyploidy superimposed on older large-scale duplications in the Arabidopsis genome. Genome Res. 2003, 13: 137-144. 10.1101/gr.751803.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  21. Arabidopsis Genome Initiative: Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature. 2000, 408: 796-815. 10.1038/35048692.

    Article  Google Scholar 

  22. Yu J, Hu SN, Wang J, Wong GKS, Li SG, Liu B, Deng YJ, Dai L, Zhou Y, Zhang XQ, et al: A draft sequence of the rice genome (Oryza sativa L. ssp indica). Science. 2002, 296: 79-92. 10.1126/science.1068037.

    Article  PubMed  CAS  Google Scholar 

  23. Goff SA, Ricke D, Lan TH, Presting G, Wang RL, Dunn M, Glazebrook J, Sessions A, Oeller P, Varma H, et al: A draft sequence of the rice genome (Oryza sativa L. ssp japonica). Science. 2002, 296: 92-100. 10.1126/science.1068275.

    Article  PubMed  CAS  Google Scholar 

  24. Blanc G, Wolfe KH: Widespread paleopolyploidy in model plant species inferred from age distributions of duplicate genes. Plant Cell. 2004, 16: 1667-1678. 10.1105/tpc.021345.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  25. Matsumoto T, Wu JZ, Kanamori H, Katayose Y, Fujisawa M, Namiki N, Mizuno H, Yamamoto K, Antonio BA, Baba T, et al: The map-based sequence of the rice genome. Nature. 2005, 436: 793-800. 10.1038/nature03895.

    Article  Google Scholar 

  26. Wong GKS, Wang J, Tao L, Tan J, Zhang JG, Passey DA, Yu J: Compositional gradients in Gramineae genes. Genome Res. 2002, 12: 851-856. 10.1101/gr.189102.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  27. Metzgar D, Bytof J, Wills C: Selection against frameshift mutations limits microsatellite expansion in coding DNA. Genome Res. 2000, 10: 72-80.

    PubMed  CAS  PubMed Central  Google Scholar 

  28. 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.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  29. Lockton S, Gaut BS: Plant conserved non-coding sequences and paralogue evolution. Trends Genet. 2005, 21: 60-65. 10.1016/j.tig.2004.11.013.

    Article  PubMed  CAS  Google Scholar 

  30. The Arabidopsis Information Resource (TAIR). [http://www.arabidopsis.org]

  31. TIGR Rice Genome Annotation. [http://rice.tigr.org]

  32. Kolpakov R, Bana G, Kucherov G: mreps: efficient and flexible detection of tandem repeats in DNA. Nucleic Acids Res. 2003, 31: 3672-3678. 10.1093/nar/gkg617.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  33. de Wachter R: The number of repeats expected in random nucleic-acid sequences and found in genes. J Theor Biol. 1981, 91: 71-98. 10.1016/0022-5193(81)90375-1.

    Article  PubMed  CAS  Google Scholar 

  34. Gene Ontology. [http://www.geneontology.org]

  35. Feller W: An Introduction to Probability Theory and its Applications. 1968, New York: John Wiley and Sons Inc

    Google Scholar 

Download references

Acknowledgements

We thank the two anonymous reviewers for helpful comments. The work was supported by a startup fund to L.Z. at Virginia Tech.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Liqing Zhang.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Reprints and permissions

About this article

Cite this article

Lawson, M.J., Zhang, L. Distinct patterns of SSR distribution in the Arabidopsis thalianaand rice genomes. Genome Biol 7, R14 (2006). https://doi.org/10.1186/gb-2006-7-2-r14

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/gb-2006-7-2-r14

Keywords