Email updates

Keep up to date with the latest news and content from Genome Biology and BioMed Central.

Open Access Research

Mutation patterns of amino acid tandem repeats in the human proteome

Loris Mularoni1, Roderic Guigó12 and M Mar Albà1*

Author Affiliations

1 Research Unit on Biomedical Informatics, Institut Municipal d'Investigació Mèdica, Universitat Pompeu Fabra, Barcelona 08003, Spain

2 Centre de Regulació Genòmica, Barcelona 08003, Spain

For all author emails, please log on.

Genome Biology 2006, 7:R33  doi:10.1186/gb-2006-7-4-r33


The electronic version of this article is the complete one and can be found online at: http://genomebiology.com/2006/7/4/R33


Received:3 February 2006
Revisions received:17 March 2006
Accepted:23 March 2006
Published:26 April 2006

© 2006 Mularoni et al.; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Amino acid tandem repeats are found in nearly one-fifth of human proteins. Abnormal expansion of these regions is associated with several human disorders. To gain further insight into the mutational mechanisms that operate in this type of sequence, we have analyzed a large number of mutation variants derived from human expressed sequence tags (ESTs).

Results

We identified 137 polymorphic variants in 115 different amino acid tandem repeats. Of these, 77 contained amino acid substitutions and 60 contained gaps (expansions or contractions of the repeat unit). The analysis showed that at least about 21% of the repeats might be polymorphic in humans. We compared the mutations found in different types of amino acid repeats and in adjacent regions. Overall, repeats showed a five-fold increase in the number of gap mutations compared to adjacent regions, reflecting the action of slippage within the repetitive structures. Gap and substitution mutations were very differently distributed between different amino acid repeat types. Among repeats containing gap variants we identified several disease and candidate disease genes.

Conclusion

This is the first report at a genome-wide scale of the types of mutations occurring in the amino acid repeat component of the human proteome. We show that the mutational dynamics of different amino acid repeat types are very diverse. We provide a list of loci with highly variable repeat structures, some of which may be potentially involved in disease.

Background

Single amino acid tandem repeats, also called homopolymeric amino acid tracts, are very abundant in eukaryotic proteins and are present in nearly one-fifth of human gene products [1,2]. They can be encoded by runs of a single codon or by a mixture of synonymous codons. Pure runs of the same codon will be susceptible to expansions and contractions of the core repetitive unit via slippage of trinucleotide repeat units [3,4]. In accordance, repeats that are poorly conserved in orthologous genes across different species are more often encoded by homogeneous codon tracts than repeats that are well conserved across species [2,5].

It has been proposed that the high mutability associated with slippage may provide an evolutionary advantage in the adaptation to new environments and to the rapid evolution of morphological traits [6,7]. But slippage can also have pathogenic effects: the uncontrolled expansion of trinucleotide repeats within human coding sequences is associated with several neurodegenerative disorders. Examples are Huntington's disease and dentatorubro-pallidolusyan atrophy, both associated with abnormally long expansions of CAG runs encoding poly-glutamine tracts (for reviews, see [8,9]). The high mutability of disease-associated repeats is reflected in high repeat size polymorphism levels in the human population [10,11]. Detection of highly variable amino acid tandem repeats can thus help discover new loci that may be particularly prone to suffer repeat expansions and become pathogenic.

Here we report on the mutations found in regions encoding amino acid tandem repeats in human genes using the human expressed sequence tag (EST) database. Of 115 different variants, each supported by at least 2 ESTs, almost half contain expansions or contractions of the amino acid repeat. We analyze the properties of repeats formed by different types of amino acids and identify a group of human genes that could potentially suffer expansions similar to those observed in the disease genes.

Results

Survey of polymorphic amino acid repeats

We analyzed 33,860 human peptide sequences from the Ensembl database [12] for the presence of tandem amino acid repeats of size 5 or longer; 5,467 proteins contained at least one such tandem repeat (about 16%). The most common amino acid repeat types (n > 200) were glutamic acid (888), proline (883), alanine (681), serine (623), glycine (510), leucine (392), glutamine (273) and lysine (223). The average repeat size was similar for different amino acids (in the range 5.8 to 6.8) except for glutamine, with longer repeats (average 8.7).

We mapped the human ESTs [13] to the repeat regions using TBLASTN [14]. We selected those repeat regions, including the tandem repeat and 15 nucleotides of adjacent sequence at each side, that were covered by at least 4 different ESTs with >90% identity matches. These comprised 2,227 repeat regions, about 41% of the initial ones, with an average repeat coverage of 27.4 ESTs per repeat. Within these, 115 (5.2%) showed one or several polymorphic variants, each supported by at least 2 different ESTs (Table 1). The amount of polymorphism varied between different amino acids, from 2.4% in leucine repeats to 10.2% in glutamine repeats. Considering only cases for which we had 100 or more ESTs, repeats that were polymorphic went up from 5.2% to 21% (26 out of 123).

Table 1. Human amino acid repeat variants

We detected 137 different polymorphic variants in the 115 polymorphic amino acid repeats. We classified them as those containing gaps (or indels) of which there were 60 (43.8%), and those containing only amino acid substitutions, of which there were 77 (56.2%) (Additional data file 1) We also measured the repeat codon homogeneity as the fraction of the repeat encoded by a perfect codon run. A high average codon homogeneity in the sequences encoding different types of amino acid repeats was generally associated with a high percentage of gap polymorphic variants (Table 1). Glutamine repeats showed the highest frequency of gap polymorphisms (88% of the glutamine polymorphic variants) and the highest average codon homogeneity (0.66), while proline and serine repeats showed the lowest gap polymorphism frequencies (14% and 18%, respectively) and the lowest average codon homogeneity (0.41 in both cases).

We also analyzed polymorphisms supported by ≥4 ESTs, which should be enriched in more common polymorphic variants; they comprised 38% of the polymorphism data. In this dataset, the frequency of repeat gap polymorphisms was higher than that of substitutions (35 versus 17).

Analysis of repeat adjacent sequences

We compared the repeat polymorphism levels to those of the sequences immediately adjacent to the repeats, considering, at each side of the repeat, a sequence of the same length as the corresponding repeat (Additional data file 1). The overall number of polymorphic variants was similar to that found within the repeats (4.8% versus 5.2%), but the number of polymorphisms containing gaps was remarkably lower, 8 in upstream regions and 14 in downstream regions, about 5 times less than within repeats (60 cases). In contrast, substitutions were slightly more common outside the repeats than inside them: 103 and 93 in upstream and downstream regions, respectively, compared to 77 within repeats.

Among polymorphisms supported by ≥4 ESTs, the trend was maintained for a larger number of gap polymorphisms within repeats than outside repeats (35 in respect to 3 in upstream and 11 in downstream sequences) and only a small difference for substitutions (21 in upstream and 26 in downstream sequences, in comparison to 17 within repeats).

Types of polymorphism by amino acid repeat type

We compared the number of polymorphisms involving gaps or substitutions in different amino acid repeat types and adjacent regions (Figure 1). This analysis showed that the previously observed larger number of substitutions outside the repeats with respect to inside the repeats could be mainly attributed to leucine repeats (10 and 12 polymorphic variants in upstream and downstream sequences, respectively, versus only one within the repeat). On the other hand, glutamine, alanine and glutamic acid repeats were the main contributors to the increased number of gap polymorphisms inside the repeats than outside the repeats. The ratio between gap polymorphisms and substitution polymorphisms was highest in the case of glutamine (15 versus 2) and lowest for proline (3 versus 22), indicating strong differences in the susceptibility to slippage of different repeat types.

thumbnailFigure 1. Number of polymorphic variants for regions containing different kinds of amino acid repeats. For the upstream and downstream sequences adjacent to the repeat the average value was taken. Bars indicate the actual values of both repeat adjacent sides.

Position and nature of amino acid substitutions

We investigated the frequency of the different amino acid substitutions in repeats and adjacent sequences, focusing on the eight most common amino acids forming tandem repeats (Table 2). The aim was to identify possible biases in the amino acid substitution patterns inside repeats with respect to the repeats' adjacent sequences, as this could be informative of specific selective constraints operating in the repetitive structures. The dataset analyzed comprised 79 substitutions inside repeats and 135 in adjacent regions. In the first place, we determined that the vast majority of amino acid substitutions could be explained by single non-synonymous nucleotide changes. Inspection of the types of amino acid substitutions in repeats and adjacent sequences indicated that there were no major differences between them. For example, nearly all amino acid substitutions that occurred at least five times in the adjacent sequences, representing the most common amino acid replacements, could also be observed inside the repeats. The only exception was the replacement of G by V, with seven cases in adjacent sequences versus none within repeats. Given the low number of cases, however, this observation should be treated with caution.

Table 2. Amino acid substitutions in polymorphic variants

In addition, we inspected the relative position of substitutions inside the repeats. This could be informative for biases in the positions where substitutions more often occurred. For example, an excess of substitutions at the repeat extremes could indicate a selective pressure to preserve a minimum length of the repeat. However, the observed position of amino acid substitutions was overall not significantly different from the expected distribution if substitutions were located at random (see Materials and methods and Additional data file 1). In conclusion, this analysis did not detect any specific differences in the selective constraints related to different amino acid substitutions inside or outside repeats, or in the relative position of the substitutions within the repeats.

Non-synonymous versus synonymous substitutions

Another aspect we studied was the relative frequency of synonymous and non-synonymous substitutions inside repeats and in repeat adjacent regions. We identified all the synonymous and non-synonymous nucleotide changes in the EST dataset and divided it by the total number of synonymous and non-synonymous positions analyzed (Figure 2 and Additional data file 1). In the two types of regions, the frequency of non-synonymous substitutions was lower than that of synonymous substitutions, as expected if some of the substitutions resulting in amino acid changes were negatively selected. The frequency of synonymous substitutions was very similar inside and outside the repeats: 0.015 (1.5% of sites) inside repeats and 0.014 to 0.016 (1.4% to 1.6% of sites) in repeat upstream and downstream regions, respectively. In agreement with the results obtained on amino acid substitutions, the frequency of non-synonymous substitutions was similar inside repeats (0.009, 0.9% of sites) and outside repeats (0.011; 1.1% of sites, in both repeat upstream and downstream regions). By amino acid type, only proline and glutamine repeats showed a non-synonymous substitution pattern different from their corresponding adjacent sequences. In the case of proline, the frequency of non-synonymous substitutions was 0.02 while the average of the two adjacent regions was 0.012. That is, there appeared to be an almost two-fold excess of non-synonymous changes within repeats. In the case of glutamine repeats, the opposite trend was observed, with a non-synonymous substitution frequency of 0.005 inside the repeats versus 0.011 in the adjacent regions.

thumbnailFigure 2. Frequency of synonymous and non-synonymous nucleotide substitutions for regions containing different kinds of amino acid repeats. For the upstream and downstream sequences adjacent to the repeat the average value was taken. Bars indicate the actual values of both repeat adjacent sides.

Relationship between polymorphism and codon homogeneity

We next compared the codon homogeneity values of all repeats to those of the repeats associated with gap polymorphisms or with substitution polymorphisms (Figure 3). The average value was 0.49 in all the repeats analyzed, 0.44 for those with substitution variants and 0.65 for those with gap variants. Repeats that showed gap polymorphisms had higher codon homogeneity values than average (p = 0.001, Kolmogorov-Smirnov test). Those with substitution variants, instead, were similar to the general repeat population. These results are expected if we consider that slippage will mainly act on long pure codon tracts, resulting in expansions and contractions of the repeats. Interestingly, however, the presence of a long pure codon tract does not appear to be indispensable for this type of polymorphism to occur, as in about 25% of the cases the longest pure codon run had a short size, between 1 and 3 codon repeat units.

thumbnailFigure 3. Codon homogeneity distribution of the sequence regions encoding different types of repeats: polymorphic with substitutions, polymorphic with expansions or contractions (gaps), all repeats. Codon homogeneity value intervals labeled as X-Y stand for values >X and <=Y (for example, 0-0.2 are values >0 and <= 0.2).

Repeat expansion/contraction polymorphisms

Polymorphic cases related to the expansion or contraction of repeats, those that involve gaps, are of particular interest because of the potential of these elements to cause disease. Table 3 lists genes containing this type of polymorphism for the most abundant amino acid repeat types. Among them we detected two poly-glutamine containing genes known to be associated with neurodegenerative disorders: dentatorubral-pallidoluysian atrophy protein (DRPLA) and spinocerebellar ataxia protein 6 (voltage-dependent P/Q-type calcium channel alpha-1A subunit, or CACNA1A). These two disease loci contained long runs (19 and 13 glutamines, respectively) and high codon homogeneity levels (0.79 and 1, respectively). Other genes in the list with long homopeptide runs and high codon homogeneity are thus possible candidates to be associated with disease. Among genes showing expansion/contraction polymorphisms was an abundance of transcription factors and RNA-binding proteins. Most of the polymorphic variants were one repeat unit away from the reference repeat, the maximum difference being three repeat units. The detection of longer repeat size variants would be hindered by our >90% identity EST match criteria, but, given that only two variants were found that show a 3 repeat unit size difference, these cases are expected to be rare. In 12 cases, the longest pure codon run occupied the totality of the repeat (codon homogeneity of 1). Length polymorphisms were most frequently associated with CAG (glutamine), GAG (glutamic acid) and CTG/GCT (leucine/alanine).

Table 3. Repeat gap polymorphic variants

Discussion

Databases of ESTs can be used to rapidly screen for potential polymorphisms in the products of eukaryotic genomes [15,16] and in particular are of great use for identifying microsatellite size variants [17-19]. We have explored this type of resource to obtain an overview of the polymorphisms associated with amino acid tandem repeats in human proteins, including potential expansion/contraction polymorphisms that may be associated with disease. We have focused on variants supported by at least two different EST sequences, discarding those associated with a single EST, to minimize the effect of possible errors introduced by the EST sequencing procedure. Other studies have focused on the detection of polymorphisms in highly homogeneous DNA repeats in coding sequences [19-21], or in specific amino acid repeat datasets [10,11,22]. While our results are generally consistent with these studies, we have taken a more extensive genome-wide approach, using all human sequence information currently available in databases to obtain a more complete picture of the types of mutations found in different amino acid repeat types. In this regard, the study by O'Dushlaine et al. [19] has points in common with our study, since ESTs were also used to infer patterns of copy number variation in protein coding genes in the human genome. In the former, however, repeat length polymorphism was investigated at the nucleotide level, whereas we investigate it at amino acid level. For this reason, while in [19] the analysis is based on UniGene clusters, and a representative sequence from each cluster is compared to all ESTs in the same cluster, we have based our analysis on the Ensembl set of proteins, and compared each of them against the entire set of ESTs. The two studies are complementary and a fraction of about 30% of the polymorphic variants identified in our study maps to polymorphic variants found in [19].

It has been suggested that the evolutionary dynamics of microsatellite-type structures can be explained by a balance between expansion by slippage and growth interruption by point mutation [23,24]. The different frequencies of gap and substitution mutations that can be observed in different types of repeats are, therefore, likely to reflect the different strength of these two evolutionary forces at the DNA level, coupled with the action of selection at the protein level. Many of the gap variants may have originated by trinucleotide slippage, as they show significantly higher levels of codon homogeneity, and this has been linked to increased repeat expansions [25] and to higher inter-specific repeat divergence [5]. Unequal recombination has also been suggested to result in large size differences in a number of disease-associated poly-alanine tracts [26,27], but it seems unlikely that it plays a major contribution here, as the variants we describe mostly diverge by one repeat unit and are biased toward long pure codon runs. Within human amino acid repeats it becomes clear that glutamine has a much higher propensity to suffer expansions than other types of repeats, with 88% of the polymorphic variants containing gaps (15 out of 17). On the contrary, proline repeats appear to be little exposed to this type of mutation, with only 14% of the polymorphic variants containing gaps (3 out of 22). In spite of the low expansion/contraction rate observed for proline repeats, which would seem to suggest a low rate of de novo formation of this kind of repeat, it is interesting to note that these are among the most common repeats. Their abundance may be related to a role in mediating protein-protein interactions, as proline-rich regions are often found in protein-protein interaction surfaces [28] and proline tandem repeats are strongly associated with 'protein binding' functional annotations [2].

Our analysis captures the elevated levels of repeat size polymorphism previously reported in poly-glutamine disease-associated loci [10,11]. Of the 4 different disease loci for which we have obtained coverage with ≥4 ESTs - spinocerebellar ataxia 6 (CACNA1A or SCA6), dentarubro-pallidolusyan atrophy, Huntington's disease and spinocerebellar ataxia 7 (SCA7) - we detected repeat size polymorphic variants for the first two. The lack of observed variability for Huntington's disease and SCA7 may be explained by their poor EST coverage, 5 and 6 ESTs, respectively. Glutamine repeats associated with human disease share a number of characteristics: they are highly polymorphic, they are among the longest tandem amino acid repeats in the proteome, and they are encoded by highly homogeneous codon runs. A fourth characteristic, previously reported, is that they are generally much shorter in rodent species than in humans, probably denoting a recent expansion in primates [29,30]. Interestingly, we have detected several other loci with similar characteristics, which could, therefore, be good candidates for involvement in trinucleotide expansion diseases. For example, the mRNA decapping enzyme 1B contains a poly-glutamine run of 10 units, codon homogeneity of 0.9, and has no detectable repeat in the rodent homologues. Another example is zinc finger protein 384 (nuclear matrix transcription factor 4), which contains a run of 16 repeat units, codon homogeneity of 0.88, and a shorter repeat of size 7 in both mouse and rat.

We observed that about 5.2% of the repeats (115 out of 2,227) show some kind of polymorphism but as the average EST coverage is only 27.4 ESTs per loci, many polymorphisms occurring in natural populations may have been missed. A closer estimate may be obtained using cases with an EST coverage of 100 or more ESTs per loci (123 repeats with average EST coverage 179.8). In this case, 21% of the repeats (26 out of 123) have at least one polymorphic variant. In a study based on a selection of highly homogeneous DNA repeats in human coding sequences, it was found that out of 42 repeats tested by PCR amplification from 36 individuals about 40% were polymorphic [21]. For comparative purposes, let's consider those repeats in our dataset with coverage of ≥100 ESTs and encoded by sequences containing pure codon repeats of size ≥5 (17 different ones). The polymorphism level within these repeats is 17.3% considering cases supported by ≥2 ESTs, and 35.5% considering those supported by ≥1 EST; the latter is similar to that obtained in [21]. In another study, the authors screened polymorphisms associated with human sequences coding for more than 7 alanines in 42 DNA samples, and determined that 24.5% (24 out of 98) had triplet expansion/contraction polymorphic variants [22]. Using the same size cutoff, we detected repeat size variants for 40% of poly-alanine repeats (2 out of 5) using ≥100 EST coverage (or 3 out of 38 (13%,)using ≥4 EST coverage). One of them, in 60S ribosomal protein L14 (RPL14), is common to both datasets.

The intra-specific variability within repeat structures has been compared to that in adjacent regions. In general, the number of gap polymorphisms in the repeat surrounding regions is five times lower than that within repeats, indicating a much more reduced slippage activity outside the repeats. However, the number of substitutions is, in general, comparable to that within the repeats. Considering that an important fraction of the repeats is likely to comprise neutral structures, many of which might have originated by slippage, it is somewhat surprising to observe a similar relaxed level of negative selection inside and outside the repeats. An exception is leucine, which shows a very small number of substitutions inside the repeat compared to the adjacent region. That would be consistent with the existence of stronger functional constraints inside the repeat. Leucine tandem repeats are often found at the amino terminus of transmembrane receptor proteins [2], where it has been suggested that they could function as signal peptides [1]. Other proteins, such as Toll-like receptors, contain leucine-rich regions that can be involved in the recognition of pathogens [31]. These putative functions could result in a reduced number of observed substitution polymorphisms.

A deeper insight into the substitution patterns in repeats and adjacent regions can be obtained by the analysis of the types of amino acid substitutions, as well as the frequencies of synonymous and non-synonymous substitutions, in the two types of regions. We have found that the vast majority of amino acid changes can be explained by a single nucleotide change, indicating a low incidence of multiple substitutions at the same site, as expected for intra-specific sequence comparisons. This analysis has also shown that a broad range of different amino acid replacements can be observed in the polymorphic variants of both repeats and adjacent sequences. The effect of selection can be better analyzed by comparing the non-synonymous and the synonymous substitution frequencies, as only the former will be related to selective constraints at the protein sequence level. Our results show that non-synonymous substitution frequencies are lower than synonymous ones, both in repeats and adjacent sequences, indicating that selection plays a role in shaping the amino acid content of these regions. The overall observed ratio is about 1 non-synonymous substitution for every 1.5 synonymous substitutions, which is higher than that typically observed in inter-specific comparisons [32]. The increased non-synonymous to synonymous substitution ratio in intra-specific measurements versus inter-specific ones is not unexpected in light of several recent reports [33,34], and one of the reasons could be the persistence in populations of slightly deleterious non-synonymous mutations that are yet to be lost [34]. In comparing repeats and adjacent regions, few differences in the non-synonymous versus synonymous substitution frequencies are observed, which, together with the data on amino acid substitutions, indicates that overall the selective constraints related to substitutions, at least at the intra-specific level, do not appear to be too different inside repeats and in the regions adjacent to them.

An interesting question for future studies will be to determine if similar conclusions can be derived from inter-specific comparisons. Interestingly, it has been previously noted that regions adjacent to poly-glutamine tracts in human and mouse proteins tend to show high divergence rates, particularly when repeats are not conserved between the two species [35]. This may indicate that repeats tend to originate in regions that are subjected to low selective constraints, where disruption of the structure or function of the protein will be less severe. Another interesting scenario is that many of the adjacent regions may indeed be old degenerate repeats, which, in the absence of selection, are in a rapidly evolving phase. The expected increase in available sequence and variability data will undoubtedly contribute to deepen our understanding of these highly mutagenic sequences.

Conclusion

We have identified a large number of human amino acid repeat variants and classified them according to the mutational mechanism, amino acid substitution or expansion/contraction, of the repeat. This has allowed us to quantify the mutation propensity of regions located within and outside tandem repeats and of repeats formed by different amino acid repeat types. The analysis has led to the identification of new candidate disease genes.

Materials and methods

Sequence databases

Human protein and cDNA sequences were extracted from the Ensembl database (NCBI35-based release) [12]. The number of initial peptide sequences was 33,860. The source of EST sequences was the NCBI-EST database (Feb 3 2005) at the National Center for Biotechnology Information [13], containing 5,430,499 EST human sequences.

Repeat count and analysis

We used our own programs to identify all single amino acid tandem repeats of size five or longer in the human proteins and to extract the DNA sequences encoding them. We identified repeats in 5,467 different proteins. For each repeat we stored the repeated amino acid, position in the sequence, repeat length, length of the longest pure codon run and codon(s) in the longest pure codon run(s). In specific cases we also retrieved the equivalent repeat in the mouse and rat orthologous sequences using BLASTP [14] at NCBI. For each repeat we calculated codon homogeneity as the fraction of the repeat occupied by the longest pure codon run. The non-parametric Kolmogorov-Smirnov test was used to assess the difference in the codon homogeneity values of different samples.

EST mapping

We mapped all human ESTs to the repeat regions in the reference proteins using the program TBLASTN [14]. The repeat regions included the perfect tandem repeat and 15 nucleotide sequences at each side of the repeat. We considered only EST matches that covered the entire repeat region and showed a percent identity >90%. This may hinder the detection of very divergent polymorphic variants but limits the chances of matches between unrelated sequences. For analysis we selected those repeat regions that were covered by at least 4 different ESTs (2,227 cases). We also retrieved cases covered by at least 100 different ESTs (123 cases).

Detection of polymorphic variants

Polymorphic variants were identified as changes to the original sequence supported by at least 2 independent ESTs (137 within repeats, 111 in upstream regions, 107 in downstream regions). We counted the different types of polymorphisms within the tandem repeats and in sequences of the same length immediately adjacent to the repeats. We discarded those cases where adjacent regions also contained repeats. For comparison we also analyzed polymorphic variants supported by at least 4 ESTs (52 within repeats, 24 in upstream regions, 37 in downstream regions). They were classified as variants involving expansions and/or contractions (gaps or indels) and variants involving only amino acid substitutions.

Type and location of amino acid substitutions

We counted the observed frequency of all possible types of amino acid substitutions, in the repeat and adjacent sequence polymorphic variants, for those amino acid repeat types that were most frequently found in tandem repeats (A, E, G, L, P, S, K, Q). No strong differences were observed in the two datasets. We also counted the position of substitutions within the repeats, by assigning each substitution to one of the following classes: pos = 1 (first position of the repeat), pos = 2 (second position), pos = -1 (last position), pos = -2 (position before the last one) and middle (remainder of positions). We calculated the expected values under a random distribution; for example, in a repeat of size 5, each class will have an expected value of 0.2, and in a repeat of size 6 all classes will have an expected value of 0.16 except the middle, which will have an expected value of 0.33. The total expected values for each group were compared with the observed values using a chi-square test. No significant differences were found.

Synonymous and non-synonymous nucleotide substitutions

We counted the observed number of synonymous and non-synonymous nucleotide substitutions in the non-redundant EST dataset matching the repeats and their adjacent regions. In this case we included substitutions represented by a single EST as well as by several identical ESTs to have a sufficiently large dataset to be able to obtain and compare substitution frequencies. Some of the changes could be due to sequencing errors. This type of error should affect both synonymous and non-synonymous substitution rates inside and outside the repeats in the same manner. As we still detected differences between the two types of rates, and as our main goal was to compare different regions and types of homopeptides, we used all the observed mutations in the non-redundant EST dataset. To maximize the reliability of the alignments we discarded ESTs containing gaps (3%). The dataset comprised 8,196 different non-redundant ESTs. We counted the number of synonymous and non-synonymous positions analyzed to obtain the frequency of substitutions of each kind. Overall, we analyzed 430,161 nucleotide positions, 107,845 of which were synonymous and 322,316 non-synonymous. The total number of synonymous substitutions was 1,663 (1.54% of sites) and of non-synonymous substitutions 3,458 (1.07% of sites). We also extracted the results for sequences containing each different amino acid repeat type.

Additional data files

The following additional data is available with the online version of this manuscript. Additional data file 1 contains a listing of substitution polymorphic variants within tandem amino acid repeats (subs_rep), a listing of gap polymorphic variants within tandem amino acid repeats (gaps_rep), a listing of substitution polymorphic variants in repeat adjacent sequences (subs_adj), a listing of gap polymorphic variants in repeat adjacent sequences (gaps_adj), data on observed and expected substitution positions (subs_position) and, data on synonymous and non-synonymous substitutions (nucl_subs).

Additional Data File 1. In the file, subs_rep is a list of substitution polymorphic variants within tandem amino acid repeats, gaps_rep is a list of gap polymorphic variants within tandem amino acid repeats, subs_adj is a list of substitution polymorphic variants in repeat adjacent sequences, gaps_adj is a list of gap polymorphic variants in repeat adjacent sequences, subs_position contains the observed and expected substitution positions and, nucl_subs contains data on synonymous and non-synonymous substitutions.

Format: XLS Size: 67KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Acknowledgements

We acknowledge support by the program Ramón y Cajal and Fundació ICREA (MMA) and from Universitá di Bologna (LM). We are grateful to Celine Poux, Nicolás Bellora and Domènec Farré for their useful comments. This research was funded by grants BIO2002-04426-C02-01 and BIO2003-05073 from Ministerio de Ciencia y Tecnología (Spain), and STAR European Project.

References

  1. Karlin S, Brocchieri L, Bergman A, Mrazek J, Gentles AJ: Amino acid runs in eukaryotic proteomes and disease associations.

    Proc Natl Acad Sci USA 2002, 99:333-338. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  2. Albà MM, Guigó R: Comparative analysis of amino acid repeats in rodents and humans.

    Genome Res 2004, 14:549-554. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  3. Tachida H, Iizuka M: Persistence of repeated sequences that evolve by replication slippage.

    Genetics 1992, 131:471-478. PubMed Abstract | Publisher Full Text OpenURL

  4. Li Y, Korol AB, Fahima T, Nevo E: Microsatellites within genes: structure, function and evolution.

    Mol Biol Evol 2004, 21:991-1007. PubMed Abstract | Publisher Full Text OpenURL

  5. Albà MM, Santibáñez-Koref MF, Hancock JM: Conservation of polyglutamine tract size between mouse and human depends on codon interruption.

    Mol Biol Evol 1999, 16:1641-1644. PubMed Abstract | Publisher Full Text OpenURL

  6. Kashi Y, King D, Soller M: Simple sequence repeats as a source of quantitative genetic variation.

    Trends Genet 1997, 13:74-78. PubMed Abstract | Publisher Full Text OpenURL

  7. Fondon JW 3rd, Garner HR: Molecular origins of rapid and continuous morphological evolution.

    Proc Natl Acad Sci USA 2004, 101:18058-18063. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  8. Wells RD: Molecular basis of genetic instability of triplet repeats.

    J Biol Chem 1996, 271:2875-2878. PubMed Abstract | Publisher Full Text OpenURL

  9. Gatchel JR, Zoghbi HY: Diseases of unstable repeat expansion: mechanisms and common principles.

    Nat Rev Genet 2005, 6:743-755. PubMed Abstract | Publisher Full Text OpenURL

  10. Jodice C, Giovannone B, Calabresi V, Bellocchi M, Terrenato L, Novelletto A: Population variation analysis at nine loci containing expressed trinucleotide repeats.

    Ann Hum Genet 1997, 61:425-438. PubMed Abstract | Publisher Full Text OpenURL

  11. Andrés AM, Lao O, Soldevila M, Calafell F, Bertranpetit J: Dynamics of CAG repeat loci revealed by the analysis of their variability.

    Hum Mutat 2003, 21:61-70. PubMed Abstract | Publisher Full Text OpenURL

  12. Hubbard T, Andrews D, Caccamo M, Cameron G, Chen Y, Clamp M, Clarke L, Coates G, Cox T, Cunningham F, et al.: Ensembl 2005.

    Nucl Acids Res 2005, 33:D447-D453. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Wheeler DL, Barrett T, Benson DA, Bryant SH, Canese K, Church DM, DiCuccio M, Edgar R, Federhen S, Helmberg W, et al.: Database resources of the National Center for Biotechnology Information.

    Nucl Acids Res 2005, 33:D39-45. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs.

    Nucleic Acids Res 1997, 25:3389-3402. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Picoult-Newberg L, Ideker TE, Pohl MG, Taylor SL, Donaldson MA, Nickerson DA, Boyce-Jacino M: Mining SNPs from EST databases.

    Genome Res 1999, 9:167-174. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Guryev V, Berezikov E, Malik R, Plasterk RH, Cuppen E, Guryev V: Single nucleotide polymorphisms associated with rat expressed sequences.

    Genome Res 2004, 14:1438-1443. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Wilder SP, Bihoreau MT, Argoud K, Watanabe TK, Lathrop M, Gauguier D: Integration of the rat recombination and EST maps in the rat genomic sequence and comparative mapping analysis with the mouse genome.

    Genome Res 2004, 14:758-765. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  18. La Rota M, Kantety RV, Yu JK, Sorrells ME: Nonrandom distribution and frequencies of genomic and EST-derived microsatellite markers in rice, wheat, and barley.

    BMC Genomics 2005, 6:23. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  19. O'Dushlaine CT, Edwards RJ, Park SD, Shields DC: Tandem repeat copy-number variation in protein-coding regions of human genes.

    Genome Biol 2005, 6:R69. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  20. Deka R, Guangyn S, Smelser D, Zhong Y, Kimmel M, Chakraborty R: Rate and directionality of mutations and effects of allele size constraints at anonymous, gene-associated, and disease-causing trinucleotide loci.

    Mol Biol Evol 1999, 16:1166-1177. PubMed Abstract | Publisher Full Text OpenURL

  21. Wren JD, Forgacs E, Fondon JW 3rd, Pertsemlidis A, Cheng SY, Gallardo T, Williams RS, Shohet RV, Minna JD, Garner HR: Repeat polymorphisms within gene regions: phenotypic and evolutionary implications.

    Am J Hum Genet 2000, 67:345-356. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Lavoie H, Debeane F, Trinh QD, Turcotte JF, Corbeil-Girard LP, Dicaire MJ, Saint-Denis A, Page M, Rouleau GA, Brais B: Polymorphism, shared functions and convergent evolution of genes with sequences coding for polyalanine domains.

    Hum Mol Genet 2003, 12:2967-2979. PubMed Abstract | Publisher Full Text OpenURL

  23. Kruglyak S, Durrett RT, Schug MD, Aquadro CF: Equilibrium distributions of microsatellite repeat length resulting from a balance between slippage events and point mutations.

    Proc Natl Acad Sci USA 1998, 95:10774-10778. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  24. Santibáñez-Koref MF, Gangeswaran R, Hancock JM: A relationship between lengths of microsatellites and nearby substutitution rates in mammalian genomes.

    Mol Biol Evol 2001, 18:2119-2123. PubMed Abstract | Publisher Full Text OpenURL

  25. Kunst CB, Leeflang EP, Iber JC, Arnheim N, Warren ST: The effect of FMR1 CGG repeat interruptions on mutation frequency as measured by sperm typing.

    J Med Genet 1997, 34:627-631. PubMed Abstract OpenURL

  26. Warren ST: Polyalanine expansion in synpolydactyly might result from unequal crossing-over of HOXD13.

    Science 1997, 275:408-409. PubMed Abstract | Publisher Full Text OpenURL

  27. Utsch B, Becker K, Brock D, Lentze MJ, Bidlingmaier F, Ludwig M: A novel stable polyalanine [poly(A)] expansion in the HOXA13 gene associated with hand-foot-genital syndrome: proper function of poly(A)-harbouring transcription factors depends on a critical repeat length?.

    Hum Genet 2002, 110:488-494. PubMed Abstract | Publisher Full Text OpenURL

  28. Kay BK, Williamson MP, Sudol M: The importance of being proline: the interaction of proline-rich motifs in signaling proteins with their cognate domains.

    FASEB J 2000, 14:231-241. PubMed Abstract | Publisher Full Text OpenURL

  29. Huang H, Winter EE, Wang H, Weinstock KG, Xing H, Goodstadt L, Stenson PD, Cooper DN, Smith D, Albà MM, et al.: Conservation of human disease genes in the rat genome.

    Genome Biol 2004, 5:R47. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  30. Gibbs RA, Weinstock GM, Metzker ML, Muzny DM, Sodergren EJ, Scherer S, Scott G, Steffen D, Worley KC, Burch PE, et al.: Genome sequence of the Brown Norway rat yields insights into mammalian evolution.

    Nature 2004, 428:493-521. PubMed Abstract | Publisher Full Text OpenURL

  31. Fujita M, Into T, Yasuda M, Okusawa T, Hamahira S, Kuroki Y, Eto A, Nisizawa T, Shibata K: Involvemen of leucine residues at positions 107, 112, and 115 in a leucine-rich repeat motif of human Toll-like receptor 2 in the recognition of diacylated lipoproteins and lipopeptides and Staphylococcus aureus peptidoglycans.

    J Immunol 2003, 171:3675-3683. PubMed Abstract | Publisher Full Text OpenURL

  32. Mouse Genome Sequencing Consortium: Initial sequencing and comparative analysis of the mouse genome.

    Nature 2002, 420:520-562. PubMed Abstract | Publisher Full Text OpenURL

  33. Ho S, Phillips MJ, Cooper A, Drummond AJ: Time dependency of molecular rate estimates and systematic overestimation of recent divergence times.

    Mol Biol Evol 2005, 22:1561-1568. PubMed Abstract | Publisher Full Text OpenURL

  34. Penny D: Relativity for molecular clocks.

    Nature 2005, 436:183-184. PubMed Abstract | Publisher Full Text OpenURL

  35. Hancock JM, Worthey EA, Santibáñez-Koref MF: A role for selection in regulating the evolutionary emergence of disease-causing and other coding CAG repeats in human and mice.

    Mol Biol Evol 2001, 18:1014-1023. PubMed Abstract | Publisher Full Text OpenURL