Skip to main content
  • Research article
  • Open access
  • Published:

Coevolution between simple sequence repeats (SSRs) and virus genome size

Abstract

Background

Relationship between the level of repetitiveness in genomic sequence and genome size has been investigated by making use of complete prokaryotic and eukaryotic genomes, but relevant studies have been rarely made in virus genomes.

Results

In this study, a total of 257 viruses were examined, which cover 90% of genera. The results showed that simple sequence repeats (SSRs) is strongly, positively and significantly correlated with genome size. Certain repeat class is distributed in a certain range of genome sequence length. Mono-, di- and tri- repeats are widely distributed in all virus genomes, tetra- SSRs as a common component consist in genomes which more than 100 kb in size; in the range of genome < 100 kb, genomes containing penta- and hexa- SSRs are not more than 50%. Principal components analysis (PCA) indicated that dinucleotide repeat affects the differences of SSRs most strongly among virus genomes. Results showed that SSRs tend to accumulate in larger virus genomes; and the longer genome sequence, the longer repeat units.

Conclusions

We conducted this research standing on the height of the whole virus. We concluded that genome size is an important factor in affecting the occurrence of SSRs; hosts are also responsible for the variances of SSRs content to a certain degree.

Background

Viruses are small infectious agents, which are found wherever there is a life and have probably existed since living cells first evolved[1, 2]. There are millions of virus types[3]. Wherein, those virus species which have been reported were sorted into dsDNA, ssDNA, dsDNA-RT, ssRNA-RT, dsRNA, (−)ssRNA and (+)ssRNA viruses based on their genome types; they can also be sorted into algae, archaea, bacteria, fungi, invertebrates, plants, protozoa and vertebrates viruses based on the general host categories according to the ICTV (International Committee on the Taxonomy of Viruses)[4]. These viruses can infect all types of organisms including archaea, bacteria, plants and animals[5]. Many common human diseases are caused by viruses, such as common cold, influenza, chickenpox, cold scores, etc. In addition, many serious diseases such as ebola, AIDS, avian influenza and SARS are also caused by viruses. What's more, many genotypes of viruses are responsible for cancers, for example, human papillomavirus, hepatitis B virus, hepatitis C virus, Epstein-Barr virus, Kaposi's sarcoma-associated herpesvirus and human T-lymphotropic virus, and so on (http://en.wikipedia.org/wiki/Virus). Though there are three main theories on the origin of virus: regressive, cellular and coevolution origin theory, it is still unclear how viruses originated because they do not like other organisms forming fossils[6, 7]. So studying viruses via molecular information has been the most useful means in investigating how they arose and evolved[6, 8–10]. Success of viral genome researches will promote our understandings and solutions of numerous problems, including their origin, evolution, infection mechanism, disease treatment, etc. The genome sizes (defined as haploid DNA content) of viruses vary greatly between species. The smallest viral genomes — the ssDNA circoviruses, family Circoviridae — code for only two proteins and have a genome size of only 2 kb; the largest — miniviruses have genome sizes of over 1.2 Mb and code for over one thousand proteins[11, 12]. Two main mechanisms have been implicated in changes of genome size: one is the accumulation of transposable elements[13, 14]; the other is the accumulation of tandemly repetitive sequences[15].

Simple sequence repeats (SSRs), also known as microsatellites, generally defined as simple sequences of 1–6 nucleotides that are repeated multiple times and are present in both coding and non-coding regions of the genome[16, 17]. SSRs are ubiquitous and highly abundant in eukaryotic[18–21] and prokaryotic genomes[22, 23]. DNA repeats are primarily expanded by three models: replication, repair and recombination[24]. Meiotic recombination plays a key role in the maintenance of sequence diversity in the human genome, and SSRs have been reported to be hot spots for recombination as well as sites for random integration[25, 26]. Thus, alterations in SSRs lie at the center of DNA evolution and sequence diversity that drives adaptation; on the other hand, changes in repetitive sequences can result in deleterious effects on gene expression and function, leading to diseases[17]. The instability of SSRs was identified to be a pathway to lead to colorectal cancer[27]. It is now accepted that unstable maintenance of microsatellites occurs in about 15% of sporadic colorectal cancers[28, 29]. Microsatellite instability is also frequently associated with other diseases such as ovarian cancers, malignant tumors of endometrium[30], small intestine[29], stomach[31], skin[32] and brain, etc. The features of microsatellite instability observed in bacteria, yeast, mice and man can provide general clues as to how genomes evolve and how certain instability could contribute to human disease[17]. Some pathogens use SSRs in a strategy that counteracts the host immune response by increasing the antigenic variance of the pathogen population[33].

Genome sequences with diverse lengths make it possible to investigate the relationship between genome size and accumulation of SSRs in all virus genera whose complete genome sequences have been reported. Therefore, scatter plots and regression analysis were performed to survey the correlation between repetitiveness (SSRs occurrence as well as SSRs length) and genome size. Distributions of different repeat classes were also surveyed among virus genomes of various sizes. While, relative abundance and relative density were examined to make the SSRs comparison parallel among differently sized species genomes; principal component analysis (PCA) was designed to investigate which repeat class(es) made a greater contribution to the variance among virus species as well as the relationships between repeat classes.

Methods

Genome sequences

The Eighth Report of ICTV (International Committee on Taxonomy of Viruses) provided information on 3 orders, 73 families, 9 subfamilies, 287 genera and 1938 virus species[4]; wherein 257 genera have been reported on complete genome sequences on NCBI and one typical species was identified as the representative for each genus according to the Listing in Taxonomic Order (http://ictvdb.bio-mirror.cn/Ictv/index.htm). Therefore, the 257 genome sequences were selected as samples for the analysis of relationship between SSRs distribution and genome size in the level of the whole virus. All the genome sequences were downloaded in both Genbank and FASTA formats from the NCBI (ftp://ncbi.nlm.nih.gov/genbank/). Sequences obtained include DNA and RNA, so both T and U bases were represented with T. Some genomes were segmented, multipartite and consist of two or more segments with various sizes (Additional file 1).

SSRs extraction

SSRs were identified and localized using the software SSR Identification Tool (SSRIT), which identifies perfect di-, tri-, tetra-, penta- and hexanucleotide repeats. We have considered only those repeats, wherein the motif was repeated more than 3 times for further analysis. Mononucleotide repeats (with a repeat length of 6 nt) were identified using the tool IMEX (Imperfect Microsatellite Extractor), which can extract perfect microsatellites as well as imperfect microsatellites. Here we presented the data for all perfect repeat types. No distinctions between the occurrence of repeats in coding and noncoding regions were made, the rationale for this decision was that the coding regions often account for the large proportion (mean value approximately 90%); while the sequences of noncoding regions are usually very short; moreover, the overlap phenomenon is very common in virus genomes, and many of the details were presented inAdditional file 1.

Relative abundance and relative density

These total numbers have been normalized by using relative abundance and relative density of SSRs to allow the comparisons to be parallel among genome sequences with different sizes. Relative abundance was calculated by dividing the number of SSRs by kilo base pair (kb) of sequences; and relative density (bp/kb) was calculated by dividing the total sequences analyzed (kb) by the number of base pairs of sequence contributed by each SSR.

PCA

Principal Components Analysis (PCA) is a well known statistical technique which has wide ranging applications. The main goal of PCA is to reduce the dimensionality by decomposing the total variances observed in an original data set. That is to say, we use PCA method to transform a set of original variables into a set of new and uncorrelated variables. The mathematic principle of PCA method lies in coordinate conversion. Consequently, PC (principal component) is a linear combination of the original variables.

Mathematical model. If the sample size isn, and each sample has P observed index (X 1,  X 2, …, Xp), we can get the following matrix of the original dataset:

X = x 11 x 12 … x 1 p x 21 x 22 … x 2 p ⋮ ⋮ ⋮ x n 1 x n 2 … x np = X 1 , X 2 , … , X P

Wherein, X i = x 1 i x 2 i ⋮ x ni , i = 1 , 2 , … , p .

Making linear combinations using the p variables (X 1,  X 2, …, Xp) of the original data matrix X:

Y = { Y 1 = e 11 X 1 + e 21 X 2 + … + e p 1 X p Y 2 = e 12 X 1 + e 22 X 2 + … + e p 2 X p … Y p = e 1 p X 1 + e 2 p X 2 + … + e pp X p

Hence, Y i = e 1 i X 1 + e 2 i X 2 + … + e pi X p , i = 1 , 2 , … , p

Here, Yi is the principal component, but it must meet the following conditions: (1) e 1 i 2 + e 2 i 2 + … + e pi 2 = 1 , i = 1 , 2 , … , p ; (2) there is no correlation between Yi and Yj ( i ≠ j , i , j = 1 , 2 , … , p ); (3) the variance of Yi is the maximum during  Y i , Y i + 1 , … , Y p ; (4) V a r Y 1 + V a r Y 2 + … + V a r Y p = V a r X 1 + V a r X 2 + … + V a r X p

Geometric meaning. Supposing that the sample contains n individuals, each individual has two variables X 1,  X 2, and in addition, variables subject to the normal distribution. That is, we discuss the geometric meaning of PCA by using bivariate normally distributed variables. Therefore, scatters of sample are roughly distributed in the shape of ellipse (Figure1). Then orthogonally rotate the original plane rectangular coordinates composed of X 1 and X 2 with an angle θ, thus, two original correlated variables (X 1,  X 2) were transformed into two integrated and uncorrelated variables (Y 1,  X 2), and the correlation between the original and new axes is as following:

{ F 1 = X 1 cos θ + X 2 sin θ F 2 = − X 1 sin θ + X 2 cos θ
Figure 1
figure 1

Geometric meaning of PCA explained by using bivariate normally distributed variables. Scatters of sample are distributed in the shape of ellipse roughly, then orthogonally rotate the original plane rectangular coordinates composed of X 1 and X 2 with an angle θ. By now, two original correlated variables(X 1 , X 2 )were transformed into two integrated and uncorrelated variables (Y1,Y 2 ). Because the variance of the original variables is greater in Y 1 axis than in Y 2 axis, so the minimum of information will be lost if integrated variable Y 1 is used for replacing all original variables. Hence,Y 1 is defined as the first principal component; in contrast, variance of variables is smaller in Y 2 axis, and it can explain minor information relative to Y 1 , soY 2 is called the second principal component.

Because the variance of the original variables is greater in Y 1 axis than in Y 2 axis, so a minimum of information will be lost if integrated variable Y 1 is used for replacing all original variables. Hence, Y 1 is defined as the first principal component; in contrast, the variance of variables is smaller in Y 2 axis, and it can explain minor information relative to Y 1, so Y 2 is called the second principal component.

Results

To obtain an expansive and unbiased data set, all virus genera with complete genome sequences reported on NCBI were scanned for SSRs analysis; wherein, one typical species was selected as the representative for each genus according to the ICTVdb (http://ictvdb.bio-mirror.cn/Ictv/index.htm). Therefore, we analyzed perfect SSRs over 6 bp long, from the 257 completely sequenced virus genomes. While, the genome size varies widely, ranging from 1682 bp (S170-(−)ssRNA-31, Hepatitis delta virus, NC_001653) to 407339 bp (S42-dsDNA-42, Emiliania huxleyi virus 86, NC_007346) (Additional file 1).

Relationship between SSRs and genome size

We constructed two sets of scatter plots and then performed regression analysis of SSRs (occurrence and length) versus complete genome size for all analyzed viruses to examine the relationship between SSRs and genome size. Above all, scatter plots were made, in which, genome size was taken as an independent variable, and all analyzed data were split into two groups (genome > 30000 bp and ≤ 30000 bp) to make the scatters and curves natural and visible (Figures2,3); and then 10 curves (linear, logarithmic, inverse, quadratic, cubic, compound, power, S, growth and exponential) were fitted according to their respective mathematical models by using the software SPSS 17.0. Parameter estimates and visual inspection showed that goodness fit of data varies greatly to different models; nevertheless, curves with the best goodness of fit were picked out for correlation analysis between SSRs (occurrence and length) and genome size (Figures2,3). The number of repeat arrays varies from 4 in Nodamura virus genome (S206-(+)ssRNA-36) to 3823 in Amsacta moorei entomopoxvirus 'L' genome (S33-dsDNA-33) (Additional file 2). The power function model provides the best fitted values towards all studied SSRs occurrence and genome size by regression analysis, and results display a very strong and significant positive relationship between the occurrence of SSRs and genome size clearly (R2 = 0.919, P < 0.001) (Figure2A). Power function and cubic model best fit for the data of genome > 30000 bp and ≤ 30000 bp group, respectively (Figure2B,C). Clearly, the SSRs occurrence is strongly, significantly and positively related to the genome size in both genome > 30000 bp (R2 = 0.815, P < 0.001) and ≤ 30000 bp (R2 = 0.718, P < 0.001) group. Especially in the group of genome ≤ 30000 bp, the values of SSR occurrences fluctuate with a relatively narrow range. An exceptional case is worth noting. One point of the scatter plot locating far above the fitted curve represents the value of SSRs in Amsacta moorei entomopoxvirus 'L' genome (S33-dsDNA-33, NC_002520) with the size of 232392 bp, in which the SSRs occurrence is a total of 3823, far more than SSRs in any other analyzed virus genome.

Figure 2
figure 2

Regression analysis of relationship between SSRs occurrence and genome size. (A) Scatter plot of SSRs occurrences in all analyzed virus genomes. (B) Scatter plot of SSRs occurrences in analyzed virus genomes > 30000 bp. (C) Scatter plot of SSRs occurrences in analyzed virus genomes < 30000 bp.

Figure 3
figure 3

Regression analysis of relationship between SSRs length and genome size.

The length of SSRs varies from 27 bp in Nodamura virus genome (S206-(+)ssRNA-36) to 26829 bp in Amsacta moorei entomopoxvirus 'L' genome (S33-dsDNA-33); and the percentage of SSRs varies from 0.59% in Nodamura virus genome (S206-(+)ssRNA-36) to 11.54% in Amsacta moorei entomopoxvirus 'L' genome (S33-dsDNA-33) (Additional file 3). Similarly, we investigated the correlation between SSRs length and genome size. Figure3 showed that the distribution of SSRs length is similar to the SSRs occurrence in differently-sized genomes, and it indicated that SSRs length is also significantly and positively correlated with the genome size to all analyzed data (R2 = 0.915, P < 0.001), to genome >30000 bp group (R2 = 0.818, P < 0.001) and to genome ≤ 30000 bp (R2 = 0.705, P < 0.001) group. Likewise, Amsacta moorei entomopoxvirus 'L' genome (S33-dsDNA-33, NC_002520) shows features out of the ordinary, with the total SSRs length of 26829 bp and SSRs percentage of 11.54%, occupying the number-one spot in length and percentage of SSRs among all analyzed virus genomes. Except that, other points float up and down the curve with a small range (Figure3). The above results indicated that genome size is an important factor in affecting repetitiveness of microsatellites in viruses.

Relationship between repeat class and genome size

We surveyed the distribution of different SSR classes in virus genomes to investigate the relationship between repeat classes (mono-, di-, tri-, tetra-, penta- and hexa-) and genome sequence length. The data of genome size < 2 kb group are not in our consideration here, because too small sample sizes lead to statistical insignificance. Data presents such a trend that, for the same repeat class, the ratio of genomes with corresponding SSRs to all genomes increases with the genome sequence growing, although the genome distribution is uneven among different genome ranges (Table1). For example, the ratio of genomes with hexanucleotide SSRs is 0 in group of 2 ~ 5 kb, and it is 1.1% in 5 ~ 10 kb, 2.6% in 10 ~ 20 kb, 6.7% in 30 ~ 100 kb and 63.9% in > 100 kb group, respectively. For the same range of genome sizes, tendency seems to be that the ratio decreases with the increase of the length of repeat unit. For example, in the genome range of 10 ~ 30 kb, the ratio is 100% (mono-), 100% (di-), 98.7% (tri-), 19.2% (tetra-), 2.6% (penta-) and 2.6% (hexa-), respectively. Observed value per virus genome showed a rising trend with the increase of the genome sequence. Additionally, long repeat units such as penta- and hexa- SSRs were rarely, or even not, observed in small genomes, and certain repeat unit class distributed in genomes with a certain range of sequence length. All mono- and di- repeats were observed in analyzed genomes except Duck hepatitis B virus (S103-dsDNA-RT-2), Cryphonectria parasitica mitovirus 1 (S174-(+)ssRNA-4) and Nodamura virus (S206-(+)ssRNA-36) in which mono- repeats were not found; tri- repeats seemed to widely distribute in all virus genomes; and tetra- SSRs, as a common component, consist in genomes with size more than 100 kb (94.4% of the virus genomes contain tetra- in group of genome >100 kb); In contrast, it is rarely observed in genomes with size < 100 kb; and genomes containing penta- and hexa- SSRs are not more than 50% in < 100 kb group. Moreover, the number of tetra-, penta- and hexa- SSRs is very small in genome range of < 100 kb (Table1). Results indicated that the correlation is strong between length of repeat unit and genome size. The longer the genome sequence, the longer repeat units. For the same repeat unit class such as mononucleotide SSRs, the number of SSRs increases with the genome length increasing. It confirmed a preference that SSRs tend to accumulate in larger virus genomes.

Table 1 Distribution of repeat classes in different ranges of genome size

Relative abundance and relative density of SSRs

Because of the irregular sizes of analyzed virus genomes, we calculated the relative abundance and relative density of SSRs to make the comparison of SSRs abundance parallel among differently-sized genomes. Frequency of virus genomes with the SSRs relative abundance of 2.0 ~ 6.0 is quite high with the value of 212 (82.8% of all analyzed viruses). Wherein, 108 genomes (42.2% of all analyzed viruses) were found to have the SSRs relative abundance of 3.0 ~ 4.5. However, genomes with the SSRs relative abundance of < 2.0 and > 6.0 are relatively fewer (with the total number of 44, accounting for 17.2% of all analyzed viruses) (Figure4,Additional file 4). Paralleling, frequency of genomes is relatively high in the SSRs relative density range of 12 ~ 44 bp/kb with the genome number of 226 (88.3% of all analyzed viruses), and 147 genomes (57.4%) have the SSRs relative density among 16 ~ 32 bp/kb; moreover, 85 genomes (33.2%) have the SSRs density of 20 ~ 28 bp/kb (Figure5,Additional file 5). The relationship between SSRs relative abundance, relative density and genome size were investigated respectively. Scatter plots showed that the correlations between the SSRs relative abundance and genome size and between the relative density and genome size are quite weak (Additional file 6,Additional file 7). The results indicated that the genome size has slightly affected the relative abundance and relative density of SSRs in virus genomes. Chen et al.[34] also found that the relative abundance and relative density of SSRs were not significantly related to genome size. On the contrary, SSRs are distributed in the virus genomes with a certain proportion.

Figure 4
figure 4

Histogram of SSRs relative abundance. The horizontal axis represents the relative abundance of SSRs in all analyzed virus genomes. The vertical axis represents the genome frequency with the corresponding SSRs relative abundance. The definition of relative abundance of SSRs can be seen in MATERIAL AND METHODS.

Figure 5
figure 5

Histogram of SSRs relative density.

PCA applying to SSRs study

PCA was used to examine which factor(s) primarily lead(s) to differences in SSRs abundance among the virus species. The sample with the size of 257 (n = 257 virus genomes) contains 6 variables (p = 6, including the percentages of mono-, di-, tri-, tetra-, penta-, hexa-, respectively). Di- SSRs is the most and hexa- SSRs is the least on average, but the standard deviation is very large for each repeat unit class among the virus genomes (Additional file 8). Even so, correlation is still strong and extremely significant between the original variables (Additional file 9). The results showed that the two principal components with eigenvalues of 4.041 and 0.811 together can account for 80.869% of all differences of SSRs abundance among viruses. Wherein, the first component can account for 67.351% and the second 13.518% of all variances, respectively. Other components played a less important role in explaining the differences of SSRs abundance among virus genomes. The comparison of the parameters' coefficients for the first and second components showed that the first component has a major loading on the difference of SSRs during analyzing genomes (Table2). The results indicated that the SSRs differences among virus genomes are mainly due to the following parameters: mono-, di-, tri- and tetra-. Wherein, the variable of di- affects the differences of SSRs among virus genomes most strongly with the loading of 0.939, followed by tri-, mono- and tetra-. In this component, penta- and hexa- played relative minor role in explaining the differences of SSRs among virus genomes. In the second component, hexa- with high positive coefficient and tetra-, penta- with negative coefficients showed different effects on SSRs abundance. In this component, hexa- played the most important role in explaining differences of SSRs abundance. Overall, the results of PCA indicated that di- affected the SSRs variances among virus genomes most strongly, followed by tri-, mono- and tetra-; and then by hexa-; penta- played the weakest role in influencing the variances of SSRs abundance among viruses.

Table 2 Loadings of variables on the first two extracted principal components

All results of Kaiser-Meyer-Olkin (KMO), Bartlett's and scree test indicated that it is significantly meaningful to analyze our data using PCA (Table2). The KMO measure with the value of 0.866 is close to 1, and Bartlett's test (< 0.001) approximates to 0, and scree plot displays the "cliff" and the "screes" vividly (Additional file 10). Moreover, the correlation is strong between the original variables (Additional file 9).

Preference of SSRs

SSRs vary greatly in repeat classes and motifs among analyzed virus genomes (Table3,Additional file 11,Additional file 12,Additional file 13 andAdditional file 14). Dinucleotide SSRs accounts for the largest proportion of 48.68% in all repeat classes, followed by mono- (37.36%) and trinucleotide SSRs (13.11%). Both A and T mono- SSRs are much more than C and G SSRs, and they make up about 16.38%, 15.54%, 2.74% and 2.69% of all SSRs in analyzed viruses respectively. AT/TA SSRs predominate in dinucleotide repeats with the proportion of 17.27%, and it is slightly more than A and T mono- SSRs (16.38%, 15.54%); other di- repeat motifs are neck and neck in occurrence, but they are all higher than C and G mono- SSRs (Table3). Repeat motif group of AAT/ATA/ATT/TAA/TAT/TTA showed the highest percentage and AGT/ACT/CTA/GTA/TAC/TAG showed the lowest percentage in tri- SSRs. Tetra-, penta- and hexanucleotide SSRs are rare, accounting for 0.5% more or less. It's abnormal that penta- SSRs are less than hexa- SSRs with 0.09%, which is approximately only one third of hexa- SSRs. However, it is usually assumed that the longer repeat unit, the lower frequency it occurred. Repeat motifs differ greatly among different virus genomes (details inAdditional file 11,Additional file 12,Additional 13,Additional file 14).

Table 3 Frequency of repeat motifs (group) in all analyzed virus genomes

Discussion

These analyses extend those in Chen et al.[34] in three ways: firstly, by using larger sample such that these analyses cover almost all taxonomic virus genera; secondly, by making the data more comprehensive because the genome size varies greatly, ranging from 1682 bp (S170-(−)ssRNA-31, Hepatitis delta virus, NC_001653) to 407339 bp (S42-dsDNA-42, Emiliania huxleyi virus 86, NC_007346), (Additional file 1); and thirdly, by applying statistically significant methods. The above extension made it possible to investigate the relationship between repetitiveness of microsatellites and genome size more fully and deeply.

The previous analysis[34] simply considered the correlation between microsatellites and genome size based on relatively small sample with 54 complete Hepatitis C virus (HCV) genomes, and they found that the number of SSRs is weakly correlated with genome size. We believe that Chen's result is lacking of statistical significance due to the relatively small sample size and uniform genome length. Here, the sample made up of 257 representative virus genome sequences was designed to investigate the relationship between SSRs and genome size on the level of the whole virus. The result of our data showed a very strong and significant positive relationship between the occurrence, or length of SSRs and genome size with the value of R2 = 0.919, P < 0.001 (Figure2A) and R2 = 0.915, P < 0.001 (Figure3A), respectively. That is, the longer the virus genome sequence, the more SSRs extracted. Hancock[15, 35, 36] confirmed that the simple sequence repeats were positively and significantly correlated with the genome size in both archaea and eubacteria, and SSRs accumulate preferentially in organisms with larger genomes. Moreover, there is evidence proved that short SSRs (1–4 bp length) exist in reduced genomes, but long SSRs (5–11 bp length) consist in larger genomes in prokaryotes[23]. The overall level of repetition in genomes is related to genome size and to the degree of repetition, and the entire genome accepts simple sequences in a concerted manner when its size increases[36, 37]. A relative scarcity of repeating DNA is a major factor in causing the relatively compact size of the avian genome[38, 39]. What's more, differences in genome size account for approximately 10% of the variance in genomic repetition in archaea and eubacteria[15], suggesting that other factors can also play important roles. DNA structure and base-stacking determined the number and length distributions of microsatellites in vertebrate genomes over evolutionary time[18]. Hosts are responsible for the variances of SSRs content to a certain degree. For example, with the similar genome size, viruses infecting vertebrates and invertebrates tend to be higher than viruses attacking bacteria in SSRs content, relative abundance and relative density of SSRs overall (Additional file 15). This can be explained by the following statements. Genomes of reptiles are estimated to consist of about 30-50% repeats, birds have been estimated to consist of 15-20% of repeats[40, 41], Mus musculus of 26.1%[42, 43], and 44.9% of human genome were occupied by repeats[44, 45]. While SSR tracts make up 2.4% of the E. coli genome[46], significantly less than vertebrates'. SSRs have been reported to be hot spots for recombination as well as sites for random integration[25, 26]. Thus, the increase of viral SSRs content is maybe due to combining partial genome sequences of hosts in the process of infecting vertebrates and invertebrates. As we know, hosts evolved a number of defense systems in response to the challenge from parasites. Meanwhile, the parasites evolved multiple counter-defense mechanism as well under the selection pressure from hosts. Bacteria have developed CRISPR/Cas (CRISPR, Clustered regularly interspaced short palindromic repeats; Cas, CRISPR-associated) immune system to defend against bacteriophages by cleaving their DNA[47]. Antagonistic coevolution between bacteria and their ubiquitous parasites, bacteriophage (phage), is well known[48, 49]. The genomic regions of CRISPR/Cas are hot spot of recombination, and CRISPR/Cas modules underwent rapid evolution in natural environments because of recurrent selection pressure exerted by coevolving viruses[50]. Meanwhile, viruses may combine partial CRISPR/Cas sequence in response to the counter-defense of bacteria. Therefore, it is no coincidence that SSRs content is high in both viruses that infect vertebrates and invertebrates and these hosts themselves. The recombination enhanced the virus's ability of infection and anti-immunity to a certain extent. Evolutionarily speaking, it is the result of selection in the process of interaction between viruses and hosts. It has proposed that reduced genome size represents an adaptation to the high rate of oxidative metabolism in birds, which results primarily from the demands of flight, and the relatively small genome size of birds in general may reflect the selective pressure to minimize the amount of repetitive DNA[51, 52].

Overall, the longer genome sequence, the stronger capability the genome holding long SSRs. Each type of repeat unit is distributed in a certain length range of genomes. Mono- and di- SSRs were observed in almost all analyzed virus genomes; tri- repeats appeared to widely distribute in all virus genomes but it's number is obviously less than mono- and di- SSRs; tetra- SSRs as a common component consist in genomes with size more than 100 kb (94.4% of the genomes contain tetra- SSRs in group of genome > 100 kb). In contrast, it is relatively rare in genomes with the size < 100 kb; genomes containing penta- and hexa- SSRs are not more than 50% in < 100 kb group. Moreover, the number of tetra-, penta- and hexa- SSRs is very small (Table1). Dinucleotide and trinucleotide SSRs were observed in all analyzed HIV genomes (genome size approximately 9 kb), but almost no tetra-, penta- and hexanucleotide SSRs were found[53]. Tetranucleotide SSRs are contained in 26.7% of the analyzed Potyvirus genomes (genome size approximately 10 kb), but the number of tetranucleotide SSRs is small[54]. The data of tetra-, penta- and hexanucleotide SSRs are also rare in Mycoplasma, but they are relatively sufficient in bacterial[46, 55], fungal[56], plant[57], vertebrates[39, 41] and human[58, 59]. Those results confirmed that SSRs distribution is closely related to the genome size, indeed. The accumulation of simple sequence repeats would be attributed to the results of selection in the process of evolution. It has been well known that viruses such as influenza virus, hepatitis virus and human immunodeficiency virus (HIV) have a higher mutation rate to resist drugs, vaccines and so on during the process of replication and (or) recombination, which is one of the reasons for curing flu, hepatitis and acquired immunodeficiency syndrome (AIDS) with difficulty. Moreover, viruses lack complete repair mechanisms. Therefore, long SSRs can be poorly found in viruses. In the opinion of Mrázek et al.[23], small genomes have a strong negative selection against long SSRs due to their strong constraints against expansion.

Conclusions

Genome size is an important factor in affecting the occurrence and the total length of SSRs, moreover, there is a positive correlation between them. Additionally, hosts are also responsible for the variances of SSRs content to a certain degree. For example, with similar genome sizes, viruses infecting vertebrates and invertebrates tend to be higher than viruses attacking bacteria in SSRs content, relative abundance and relative density of SSRs, overall. We inferred that maybe viruses combined partial genome sequences of hosts in infecting, resulting in relative large genome and high content of SSRs. Evolutionarily speaking, it is the result of selection in the process of interaction between viruses and hosts. Virus is a group of parasite, so studying of SSRs in viruses is helpful to the research of many etiopathogenesis of its hosts.

Misc

Xiangyan Zhao, Yonglei Tian, both authors contributed equally to this work.

References

  1. Gao L, Qi J: Whole genome molecular phylogeny of large dsDNA viruses using composition vector method. BMC Evol Biol. 2007, 7: 41-10.1186/1471-2148-7-4.

    Article  PubMed Central  PubMed  Google Scholar 

  2. Iyer LM, Balaji S, Koonin EV, Aravind L: Evolutionary genomics of nucleo-cytoplasmic large DNA viruses. Virus Res. 2006, 117: 156-184. 10.1016/j.virusres.2006.01.009.

    Article  CAS  PubMed  Google Scholar 

  3. Breitbart M, Rohwer F: Here a virus, there a virus, everywhere the same virus?. Trends Microbiol. 2005, 13: 278-284. 10.1016/j.tim.2005.04.003.

    Article  CAS  PubMed  Google Scholar 

  4. Fauquet CM, Mayo MA, Maniloff J, Desselberger U, Ball LA: Virus Taxonomy, VIIIth Report of the ICTV. 2005, Elsevier/Academic Press, London

    Google Scholar 

  5. Koonin EV, Senkevich TG, Dolja VV: The ancient Virus World and evolution of cells. 2006, 1: 29-doi:10.1186/1745-6150-1-29.

    Google Scholar 

  6. Holland J, Spindler K, Horodyski F, Grabau E, Nichol S, VandePol S: Rapid evolution of RNA genomes. Science. 1982, 215: 1577-1585. 10.1126/science.7041255.

    Article  CAS  PubMed  Google Scholar 

  7. Holmes EC: The evolution of viral emergence. Proc Natl Acad Sci USA. 2006, 103: 4803-4804. 10.1073/pnas.0601166103.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Domingo E: Viruses at the edge of adaptation. Virology. 2000, 270: 251-253. 10.1006/viro.2000.0320.

    Article  CAS  PubMed  Google Scholar 

  9. Elena SF, Lenski RE: Evolution experiments with microorganisms: the dynamics and genetic bases of adaptation. Nat Rev Genet. 2003, 4: 457-469.

    Article  CAS  PubMed  Google Scholar 

  10. Sanjuán R, Nebot MR, Chirico N, Mansky LM, Belshaw R: Viral mutation rates. J Virol. 2010, 84: 9733-9748. 10.1128/JVI.00694-10.

    Article  PubMed Central  PubMed  Google Scholar 

  11. Mrázek J, Karlin S: Distinctive features of large complex virus genomes and proteomes. Proc Natl Acad Sci USA. 2007, 104: 5127-5132. 10.1073/pnas.0700429104.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Van Etten JL, Lane LC, Dunigan DD: DNA Viruses: The Really Big Ones (Giruses). Annu Rev Microbiol. 2010, 64: 83-99. 10.1146/annurev.micro.112408.134338.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Bennetzen JL: Transposable element contributions to plant gene and genome evolution. Plant Mol Biol. 2000, 42: 251-269. 10.1023/A:1006344508454.

    Article  CAS  PubMed  Google Scholar 

  14. International Human Genome Sequencing Consortium: Initial sequencing and analysis of the human genome. Nature. 2001, 409: 860-921. 10.1038/35057062.

    Article  Google Scholar 

  15. Hancock JM: Genome size and the accumulation of simple sequence repeats: implications of new data from genome sequencing projects. Genetica. 2002, 115: 93-103. 10.1023/A:1016028332006.

    Article  CAS  PubMed  Google Scholar 

  16. 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  CAS  PubMed  Google Scholar 

  17. Kovtun IV, McMurray CT: Features of trinucleotide repeat instability in vivo. Cell Res. 2008, 18: 198-213. 10.1038/cr.2008.5.

    Article  CAS  PubMed  Google Scholar 

  18. Bacolla A, Larson JE, Collins JR, Li J, Milosavljevic A, Stenson PD, Cooper DN, Wells RD: Abundance and length of simple repeats in vertebrate genomes are determined by their structural properties. Genome Res. 2008, 18: 1545-1553. 10.1101/gr.078303.108.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Rajendrakumar P, Biswal AK, Balachandran SM, Srinivasarao K, Sundaram RM: Simple sequence repeats in organellar genomes of rice: frequency and distribution in genic and intergenic regions. Bioinformatics. 2007, 23: 1-4. 10.1093/bioinformatics/btl547.

    Article  CAS  PubMed  Google Scholar 

  20. Sharma PC, Grover A, Kahl G: Mining microsatellites in eukaryotic genomes. Trends Biotechnol. 2007, 25: 490-498. 10.1016/j.tibtech.2007.07.013.

    Article  CAS  PubMed  Google Scholar 

  21. Temnykh S, DeClerck G, Lukashova A, Lipovich L, Cartinhour S, McCouch S: Computational and experimental analysis of microsatellites in rice (Oryza sativa L.): frequency, length variation, transposon associations, and genetic marker potential. Genome Res. 2001, 11: 1441-1452. 10.1101/gr.184001.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Mrázek J: Analysis of distribution indicates diverse functions of simple sequence repeats in mycoplasma genomes. Mol Biol Evol. 2006, 23: 1370-1385. 10.1093/molbev/msk023.

    Article  PubMed  Google Scholar 

  23. Mrázek J, Guo X, Shah A: Simple sequence repeats in prokaryotic genomes. Proc Natl Acad Sci USA. 2007, 104: 8472-8477. 10.1073/pnas.0702412104.

    Article  PubMed Central  PubMed  Google Scholar 

  24. Mirkin SM: Expandable DNA repeats and human disease. Nature. 2007, 447: 932-940. 10.1038/nature05977.

    Article  CAS  PubMed  Google Scholar 

  25. Jeffreys AJ, Holloway JK, Kauppi L, May CA, Neumann R, Slingsby MT, Webb AJ: Meiotic recombination hot spots and human DNA diversity. Philos Trans R Soc Lond B Biol Sci. 2004, 359: 141-152. 10.1098/rstb.2003.1372.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  26. Yant SR, Wu X, Huang Y, Garrison B, Burgess SM, Kay MA: High-resolution genome-wide mapping of transposon integration in mammals. Mol Cell Biol. 2005, 6: 2085-2094.

    Article  Google Scholar 

  27. Söreide K, Janssen EAM, Söiland H, Körner H, Baak JPA: Microsatellite instability in colorectal cancer. Brit J Surg. 2006, 93: 395-406. 10.1002/bjs.5328.

    Article  PubMed  Google Scholar 

  28. Chan TL, Yuen ST, Kong CK, Chan YW, Chan AS, Ng WF, Tsui WY, Lo MW, Tam WY, Li VS, Leung SY: Heritable germline epimutation of MSH2 in a family with hereditary nonpolyposis colorectal cancer. Nat Genet. 2006, 38: 1178-1183. 10.1038/ng1866.

    Article  CAS  PubMed  Google Scholar 

  29. Suter CM, Martin DIK, Ward RL: Germline epimutation of MLH1 in individuals with multiple cancers. Nat Genet. 2004, 36: 497-501. 10.1038/ng1342.

    Article  CAS  PubMed  Google Scholar 

  30. Wijnen J, de Leeuw W, Vasen H, van der Klift H, Møller P, Stormorken A, Meijers-Heijboer H, Lindhout D, Menko F, Vossen S, et al: Familial endometrial cancer in female carriers of MSH6 germline mutations. Nat Genet. 1999, 23: 142-144. 10.1038/13773.

    Article  CAS  PubMed  Google Scholar 

  31. Bacani J, Zwingerman R, Nicola DN, Spencer S, Wegrynowski T, Mitchell K, Hay K, Redston M, Holowaty E, Huntsmanet D: Tumor microsatellite instability in early onset gastric cancer. J Mol Diagn. 2005, 7: 465-477. 10.1016/S1525-1578(10)60577-6.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  32. Castiglia D, Pagani E, Alvino E, Vernole P, Marra G, Cannavò E, Jiricny J, Zambruno G, D'Atri S: Biallelic somatic inactivation of the mismatch repair gene MLH1 in a primary skin melanoma. Gene Chromosome Canc. 2003, 37: 165-175. 10.1002/gcc.10193.

    Article  CAS  Google Scholar 

  33. Rocha EPC: An appraisal of the potential for illegitimate recombination in bacterial genomes and its consequences: from duplications to genome reduction. Genome Res. 2003, 13: 1123-1132. 10.1101/gr.966203.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Chen M, Tan Z, Zeng G: Microsatellite is an important component of complete Hepatitis C virus genomes. Infect Genet Evol. 2011, 11: 1646-1654. 10.1016/j.meegid.2011.06.012.

    Article  CAS  PubMed  Google Scholar 

  35. Hancock JM: The contribution of slippage-like processes to genome evolution. J Mol Evol. 1995, 41: 1038-1047.

    Article  CAS  PubMed  Google Scholar 

  36. Hancock JM: Simple sequences in a 'minimal' genome. Nat Genet. 1996, 14: 14-15. 10.1038/ng0996-14.

    Article  CAS  PubMed  Google Scholar 

  37. Hancock JM: Simple sequences and the expanding genome. BioEssays. 1996, 18: 421-425. 10.1002/bies.950180512.

    Article  CAS  PubMed  Google Scholar 

  38. Hughes AL, Piontkivska H: DNA repeat arrays in chicken and human genomes and the adaptive evolution of avian genome size. BMC Evol Biol. 2005, 5: 12-10.1186/1471-2148-5-12.

    Article  PubMed Central  PubMed  Google Scholar 

  39. Wicker T, Robertson JS, Schulze SR, Feltus FA, Magrini V, Morrison JA, Mardis ER, Wilson RK, Peterson DG, Paterson AH, et al: The repetitive landscape of the chicken genome. Genome Res. 2005, 15: 126-136. 10.1101/gr.2438005.

    Article  PubMed Central  PubMed  Google Scholar 

  40. Epplen JT, Leipoldt M, Engel W, Schmidtke J: DNA sequence organization in avian genomes. Chromosoma. 1978, 69: 307-321. 10.1007/BF00332134.

    Article  CAS  PubMed  Google Scholar 

  41. Brandström M, Ellegren H: Genome-wide analysis of microsatellite polymorphism in chicken circumventing the ascertainment bias. Genome Res. 2008, 18: 881-887. 10.1101/gr.075242.107.

    Article  PubMed Central  PubMed  Google Scholar 

  42. Waterston RH, Lindblad-Toh K, Birney E, et al: Initial sequencing and comparative analysis of the mouse genome. Nature. 2002, 420: 520-562. 10.1038/nature01262.

    Article  CAS  PubMed  Google Scholar 

  43. MGSC genome assembly release 3. ftp://ftp.ncbi.nih.gov/genomes/M_musculus/ARCHIVE/MGSCv3_Release3/Assembled_Chromosomes

  44. Warburton PE, Hasson D, Guillem F, Lescale C, Jin X, Abrusan G: Analysis of the largest tandemly repeated DNA families in the human genome. BMC genomics. 2008, 9: 533-10.1186/1471-2164-9-533.

    Article  PubMed Central  PubMed  Google Scholar 

  45. Ames D, Murphy N, Helentjaris T, Sun N, Chandler V: Comparative analyses of human single- and multilocus tandem repeats. Genetics. 2008, 179: 1693-1704. 10.1534/genetics.108.087882.

    Article  PubMed Central  PubMed  Google Scholar 

  46. Gur-Arie R, Cohen CJ, Eitan Y, Shelef L, Hallerman EM, Kashi Y: Simple sequence repeats in Escherichia coli: abundance, distribution, composition, and polymorphism. Genome Res. 2000, 10: 62-71.

    PubMed Central  CAS  PubMed  Google Scholar 

  47. Garneau JE, Dupuis M-È, Villion M, Romero DA, Barrangou R, Boyaval P, Fremaux C, Horvath P, Magadán AH, Moineau S: The CRISPR/Cas bacterial immune system cleaves bacteriophage and plasmid DNA. Nature. 2010, 468: 67-72. 10.1038/nature09523.

    Article  CAS  PubMed  Google Scholar 

  48. Gómez P, Buckling A: Bacteria-phage antagonistic coevolution in soil. Science. 2011, 332: 106-109. 10.1126/science.1198767.

    Article  PubMed  Google Scholar 

  49. Pal C, Maciá MD, Oliver A, Schachar I, Buckling A: Coevolution with viruses drives the evolution of bacterial mutation rates. Nature. 2007, 450: 1079-1081. 10.1038/nature06350.

    Article  CAS  PubMed  Google Scholar 

  50. Takeuchi N, Wolf YI, Makarova KS, Koonin EV: Nature and Intensity of Selection Pressure on CRISPR-Associated Genes. J Bacteriol. 2012, 194: 1216-1225. 10.1128/JB.06521-11.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  51. Gregory TR: A bird's-eye view of the C-value enigma: Genome size, cell size, and metabolic rate in the class Aves. Evolution. 2002, 56: 121-130.

    Article  CAS  PubMed  Google Scholar 

  52. Szarski H: Cell size and nuclear DNA content in vertebrates. Int Rev Cytol. 1976, 44: 93-111.

    Article  CAS  PubMed  Google Scholar 

  53. Chen M, Tan Z, Jiang J, Li M, Chen H, Shen G, Yu R: Similar distribution of simple sequence repeats in diverse completed Human Immunodeficiency Virus Type 1 genomes. FEBS Lett. 2009, 583: 2959-2963. 10.1016/j.febslet.2009.08.004.

    Article  CAS  PubMed  Google Scholar 

  54. Zhao X, Tan Z, Feng H, Yang R, Li M, Jiang J, Shen G, Yu R: Microsatellites in different Potyvirus genomes: Survey and analysis. Gene. 2011, 488: 52-56. 10.1016/j.gene.2011.08.016.

    Article  CAS  PubMed  Google Scholar 

  55. Chen M, Zeng G, Tan Z, Jiang M, Zhang J, Zhang C, Lu L, Lin Y, Peng J: Compound microsatellites in complete Escherichia coli genomes. FEBS Lett. 2011, 585: 1072-1076. 10.1016/j.febslet.2011.03.005.

    Article  CAS  PubMed  Google Scholar 

  56. Karaoglu H, Lee CM, Meyer W: Survey of simple sequence repeats in completed fungal genomes. Mol Biol Evol. 2005, 22: 639-649.

    Article  CAS  PubMed  Google Scholar 

  57. Hong CP, Piao ZY, Kang TW, Batley J, Yang TJ, Hur YK, Bhak J, Park BS, Edwards D, Lim YP: Genomic Distribution of Simple Sequence Repeats in Brassica rapa. Mol Cells. 2007, 23: 349-356.

    CAS  PubMed  Google Scholar 

  58. Kelkar YD, Tyekucheva S, Chiaromonte F, Makova KD: The genome-wide determinants of human and chimpanzee microsatellite evolution. Genome Res. 2008, 18: 30-38.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  59. Mayer C, Leese F, Tollrian R: Genome-wide analysis of tandem repeats in Daphnia pulex–a comparative approach. BMC genomics. 2010, 11: 277-10.1186/1471-2164-11-277.

    Article  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements

We would like to thank Chuansheng He for the language editing and 3 anonymous reviewers for constructive comments on the earlier version of the manuscript. This work was supported by the AQSIQ Scientific Program of China [2007IK255]; National Scientific and Technique Program of China [2006BAD08A13]; Hunan Scientific Program of China [2008CKI3070] and Changsha Scientific Program of China [2011 K1113021/11].

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Zhongyang Tan or Mingfu Li.

Additional information

Competing interest

The authors declare that they have no competing interests.

Authors' contributions

ZT and ML conceived and designed this study. XZ and YT performed and drafted manuscript. HF, QO, YT and YN participated in the data processing. RY, JJ, GS and RY involved in revising the manuscript critically for important intellectual content. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: List of the basic information of all analyzed viruses.(PDF 2 MB)

Additional file 2: Occurrence of SSRs in analyzed virus genomes.(PDF 736 KB)

Additional file 3: Length (bp) of SSRs in analyzed virus genomes.(PDF 601 KB)

Additional file 4: Relative abundance of SSRs in analyzed virus genomes.(PDF 794 KB)

Additional file 5: Relative density of SSRs in analyzed virus genomes.(PDF 615 KB)

12864_2012_4772_MOESM6_ESM.pdf

Additional file 6: Scatter plots of SSRs relative abundance versus genome size. (A)Scatter plot of SSRs relative abundances in all analyzed virus genomes. (B) Scatter plot of SSRs relative abundances in analyzed virus genomes with size of < 30000 bp. (C) Scatter plot of SSRs relative abundances in analyzed virus genomes with size of > 30000 bp. (PDF 56 KB)

12864_2012_4772_MOESM7_ESM.pdf

Additional file 7: Scatter plots of SSRs relative density versus genome size. (A) Scatter plot of SSRs relative densities in all analyzed virus genomes. (B) Scatter plot of SSRs relative densities in analyzed virus genomes with size of < 30000 bp. (C) Scatter plot of SSRs relative densities in analyzed virus genomes with size of > 30000 bp. (PDF 97 KB)

Additional file 8: Descriptive statistics of SSRs variables.(PDF 12 KB)

Additional file 9: Matrix of correlation coefficients and 1-tailed tests between SSRs.(PDF 21 KB)

12864_2012_4772_MOESM10_ESM.pdf

Additional file 10: Scree plot. It displays the "cliff" and the "screes" vividly, which can be visually proved that the applicability of PCA is very good to the current data set. (PDF 33 KB)

Additional file 11: Occurrence of mono- SSRs in analyzed virus genomes.(PDF 693 KB)

Additional file 12: Occurrence of di- SSRs in analyzed virus genomes.(PDF 524 KB)

Additional file 13: Occurrence of tri- SSRs in analyzed virus genomes.(PDF 617 KB)

Additional file 14: Occurrence of tetra-, penta- and hexa- SSRs in analyzed virus genomes.(PDF 435 KB)

Additional file 15: Hosts of analyzed virus genomes.(PDF 849 KB)

Authors’ original submitted files for images

Rights and permissions

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

Reprints and permissions

About this article

Cite this article

Zhao, X., Tian, Y., Yang, R. et al. Coevolution between simple sequence repeats (SSRs) and virus genome size. BMC Genomics 13, 435 (2012). https://doi.org/10.1186/1471-2164-13-435

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-13-435

Keywords