- Methodology article
- Open Access
Use of a multi-way method to analyze the amino acid composition of a conserved group of orthologous proteins in prokaryotes
BMC Bioinformatics volume 7, Article number: 257 (2006)
Amino acids in proteins are not used equally. Some of the differences in the amino acid composition of proteins are between species (mainly due to nucleotide composition and lifestyle) and some are between proteins from the same species (related to protein function, expression or subcellular localization, for example). As several factors contribute to the different amino acid usage in proteins, it is difficult both to analyze these differences and to separate the contributions made by each factor.
Using a multi-way method called Tucker3, we have analyzed the amino composition of a set of 64 orthologous groups of proteins present in 62 archaea and bacteria. This dataset corresponds to essential proteins such as ribosomal proteins, tRNA synthetases and translational initiation or elongation factors, which are common to all the species analyzed. The Tucker3 model can be used to study the amino acid variability within and between species by taking into consideration the tridimensionality of the data set. We found that the main factor behind the amino acid composition of proteins is independent of the organism or protein function analyzed. This factor must be related to the biochemical characteristics of each amino acid. The difference between the non-ribosomal proteins and the ribosomal proteins (which are rich in arginine and lysine) is the main factor behind the differences in amino acid composition within species, while G+C content and optimal growth temperature are the main factors behind the differences in amino acid usage between species.
We show that a multi-way method is useful for comparing the amino acid composition of several groups of orthologous proteins from the same group of species. This kind of dataset is extremely useful for detecting differences between and within species.
Amino acids are not used equally in proteins. In addition to the physical, chemical and biochemical differences between amino acids, which may explain why some amino acids are more used than others, there are also differences in the amino acid composition that are associated with protein function or some characteristic of the species to which they belong. Between proteins from the same species, differences in amino acid composition have been associated with expressivity , the chromosomal position of their genes [2, 3], hydrophobicity and the number of transmembrane regions , subcellular localization  and protein function such as differences between ribosomal and non-ribosomal proteins . The availability of complete proteomes from a large number of prokaryotic species has enabled a global comparison of their amino acid compositions to be made and some of the causes of amino acid variability between species to be identified [7–11]. Between species, nucleotide bias and organism lifestyle, especially whether the organisms are mesophilic or thermophilic, are the two main factors behind the variability in amino acid composition [7–11] and even in synonymous codon usage . The methodology used in these analyses involves comparing the means and evaluating which differences are statistically significant or using principal component analysis (PCA) or correspondence analysis to determine which factors influence the amino acid composition of proteins.
Most analyses that compare the amino acid composition between species use whole proteomic averages. However, it could be interesting to compare the amino acid composition of the same group of proteins from different species to check whether the predicted bias is widespread and affects virtually all genes within a genome. Any differences found will then definitely be due to the use of different amino acids between species, not to differences in protein functions. Several authors have used a set of homologous genes from different genomes to show that the observed amino acid bias is a general trend for the evolution of proteins [7, 13]. Here we present a similar approach in which we compare the amino acid composition of a set of 64 groups of orthologous proteins common to 62 archaeal and bacterial species. The additional values of our approach, however, are the large number of species used, the fact that the sets of orthologous sequences used correspond to different functional classes, and the procedure used. This procedure consisted of a multi-way method called Tucker3 that analyzes the amino acid variability within and between species using a three-dimensional matrix as an input. In this matrix, the rows represented the frequency with which each amino acid is used, the columns represented groups of orthologous proteins, and the third dimension represented the 62 species analyzed. Briefly, the Tucker3 algorithm decomposes the three-dimensional matrix into a matrix of residuals, three component matrices A, B, C called loadings matrices, and a 3-way core array. This core array defines how the individual loadings vectors interact, providing which factor combinations best represent the data set in terms of explained variability. The order of the core array, i.e. the number of components that are calculated from each loadings matrix, needs to be determined from a priori knowledge of the data or by evaluating models with different combinations and choosing the order that provides the most accurate model. See the Methods section for a more detailed explanation of the Tucker3 algorithm.
To find the optimal order of the Tucker3 model, we studied several combinations of different orders. We found that the 5 × 3 × 3 model was the most suitable because the other models did not improve in fit when higher orders were used. Table 1 shows the elements of the core array obtained in the model of order 5 × 3 × 3 that describe the relationships between the three modes (one for amino acid composition, one for protein function and one for the organisms). The elements are sorted by the percentage of fraction variance that they represent. Only those with the five highest values are shown. With only these five elements, more than 98% of the summed fraction of the variability of the dataset is represented. For all of these five elements, we have also found a biological interpretation and a correlation with some biological character.
There is a general amino acid usage that is independent of the function or organism analyzed
The [1,1,1] factor combination (which means the first component of the amino acid composition, protein function and organism loadings matrices, respectively) explains almost 80% of the variability of the dataset. This value is much higher than the other elements of the core array. Figure 1 shows the representation of the [1,1,1] against the [2,2,2] combination (the [2,2,2] combination is discussed in the next section). Points representing protein function and organism (grey squares and blue circles, respectively) have a similar position in the [1,1,1] combination and appear as a vertical line. The projections of the amino acid variables onto the horizontal axis show that the amino acids are not used equally. The amino acids on the right (A, V, L, K, G, R, E and I) are the ones used with the highest frequency in the data set. The amino acids on the left (C and W) are the ones used with the lowest frequency. This amino acid composition is very similar to the composition of databases of protein sequences. So, for example, the projections of the twenty amino acids onto the [1,1,1] element correlates very well (r = 0.917) with the amino acid abundance extracted from Swiss-Prot release 46.5. These results show that there is a general amino acid usage in the proteins analyzed and that this is independent of the protein function or organism to which they belong. This important observation is only possible if a multi-way algorithm is applied to a three-dimensional dataset.
In the dataset analyzed, ribosomal vs non-ribosomal is the main factor behind variability in the amino acid usage between different protein functions
The [3,2,1] factor combination has a variance fraction of 6.69%. This represents the second element of the core array with the highest variance fraction. This element is represented in figure 2. The loadings values of the organisms, represented by green triangles in figure 2, are localized at the same coordinates. This means that the differences in amino acid composition related to protein function are a general trend and are independent of the organisms analyzed.
Ribosomal and non-ribosomal proteins are represented by red squares and blue circles, respectively. Each group forms a different cluster, which shows that these proteins have different amino acid usages. The good correlation (r = 0.902) between the amino acid differences between ribosomal proteins and non-ribosomal proteins and the projections of the amino acid variables onto the horizontal axis of figure 2 shows that, except for S2, L10 and L29, the ribosomal proteins are distinguishable from the non-ribosomal proteins if we compare their amino acid compositions. Figure 2 shows that the ribosomal proteins use the basic amino acids lysine (K) and arginine (R) and the small hydrophobic amino acids glycine (G) and valine (V) with the greatest frequency and use leucine (L) and the negatively charged amino acids glutamate (E) and aspartate (D) with the lowest frequency. The ribosomal proteins L10, L29, S2 and S14, however, do not cluster with the other ribosomal proteins in figure 2. S2 and L10 appear in the non-ribosomal cluster, and S14 and L29 have a slightly different amino acid composition from those of the other proteins (see figure 2). These differences may be due to the characteristics of these ribosomal proteins in the position or role when the ribosome is formed.
G+C composition and optimal growth temperature are the main factors behind variability in amino acid usage between organisms
Figure 3 represents the variability in amino acid usage between organisms. The position of the loadings related to protein function shows that the differences in amino acid composition related to organisms are independent of the proteins analyzed.
The [2,1,2] combination is the element of the core array with the third highest percentage of variance fraction. The position of the organisms in this element (i.e. the projections of this variable onto the horizontal axis) correlates very well (r = -0.901) with the G+C content of the organisms. This means that the organisms represented on the left of figure 3 (e.g. Halobacterium sp, Hbs) have the highest G+C content and those on the right (e.g. Methanococcus jannaschii, Mja) have the lowest. Because of the genetic code, the organisms with highest G+C contents use the amino acids glycine (G), alanine (A), arginine (R) and proline (P) [7, 8] with the greatest frequency. These amino acids are encoded by codons with a G or C in the first and second codon positions. This is confirmed in figure 3, where the lowest values for the amino acid loadings in the [2,1,2] factor combination correspond to A, R, G and P. The position of the amino acid valine (V) in this axis, close to the position of proline (P), is also interesting. Valine is encoded by GTN codons, but there are more valines in G+C rich species than expected . This may be due to the many conservative replacements between valine and isoleucine (encoded by AT [A,T,C]) forced by positive GC pressure . On the other hand, figure 3 also shows that lysine (K), isoleucine (I) and asparagine (N) – three amino acids encoded by codons with T or A in the first and second codon positions – are the amino acids that are most used in species with the lowest G+C contents.
The [5,1,3] combination is also represented in figure 3. The positions of the organisms on the y axis correlate well (r = -0.840) with their optimal growth temperatures. Thermophiles (species with an optimal growth temperature above 60°C) and mesophiles form two clusters in this figure, so they can be distinguished by using their amino acid composition. From the position of the loadings values, we can deduce that the amino acids preferred by thermophiles are mainly glutamate (E) and valine (V). The positions of Halobacterium sp (Hbs) and Methanosarcina acetivorans (Mac) in figure 3 are interesting. The optimal growth temperatures of these Euryarchaea are below 40°C, so they cannot be classified as thermophiles. However, they have a similar amino acid composition to that of thermophiles and cluster with them in figure 3.
Finally, the [2,2,2] factor combination shown in figure 1 explains only 0.66% of the variability of the dataset. This factor is a combination of the effects of the G+C content of the organisms and the variability in the G+C contents of the various orthologous genes. Although the G+C content of a gene depends mainly on the G+C content of the organism to which it belongs, there are also variations in the G+C content of genes from the same organism. One of the reasons for these G+C variations within organisms may be compositional amino acid constraints. This is the case of ribosomal proteins, which, because of their compositional amino acid constraints, have a lower G+C content at the second codon position than non-ribosomal proteins. This effect is reflected in the correlation (r = -0.680) between the projections onto the vertical axis of figure 1 of the variables associated with the protein functions and the mean G+C content at the second codon position of each group of orthologous genes.
The high percentage of the variance explained by the [1,1,1] factor combination shows that the main determinant of the amino acid composition of proteins is independent of the protein function or organism to which they belong. The different uses of amino acids may therefore be due to differences in several biochemical characteristics. However, which amino acid properties influence their usage in proteins is still unknown. Moreover, Jordan et al.  have shown that the amino acid composition of proteins is not in equilibrium. By comparing sets of orthologous proteins of closely related genomes from 15 species representing the three domains of life and comparing the fluxes of reciprocal substitutions caused by single-nucleotide replacements, these authors found that cysteine, methionine, histidine, serine and phenylalanine are strong 'gainers' (i.e. their frequency is increasing), and proline, alanine, glutamate and glycine are strong 'losers' (i.e. their frequency is decreasing) . Except for methionine, gainers tend to be under-represented and losers are over-represented . This loser-rich and gainer-poor amino acid composition may be due to the order in which amino acids were recruited into the genetic code . The correlation between the general amino acid frequencies that we observe (the projection of the amino acids onto the x axis in figure 1) and the rate of gain or loss defined by Jordan et al.  is only -0.39. The correlation, however, is -0.68 when we compare the general amino acid frequencies and a consensus chronology of incorporation of amino acids into the genetic code defined by Trifonov . This relatively high correlation value means that the order of recruitment of the amino acids into the genetic code can be an additional factor that influences the different use of the amino acids. However, because Trifonov used the amino acid composition of extant proteins as one of the 60 criteria to obtain his consensus chronology of amino acids , the above correlation is not unexpected and must be interpreted with caution.
In addition to this general amino acid composition, there are obviously differences in the amino acid composition of proteins due to the function or organism to which they belong. The difference between ribosomal and non-ribosomal proteins is the main factor behind the amino acid usage within species in the data set we analyzed. Shape and charge complementarity rather than sequence-specific interactions are responsible for the specific interactions of most ribosomal proteins with RNA . Because of these interactions, ribosomal proteins prefer positively charged amino acids and avoid negatively charged ones . The mapping of conserved arginines and lysines onto the ribosome structure has revealed that these charged residues frequently form surface patches that reflect RNA-binding sites . The ribosomal proteins L10, L29, S2 and S14, however, do not cluster with the other ribosomal proteins in figure 2. S2 and L10 appear in the non-ribosomal cluster, and S14 and L29 have a slightly different amino acid composition from the other proteins (see figure 2). These differences may be due to the characteristics of these ribosomal proteins in the position or role when the ribosome is formed. Although S2 is one of the largest ribosomal proteins in the 30S subunit, it is very loosely attached to this subunit (only seven out of 236 residues contact with the rRNA) and has the lowest percentage of arginine and lysine in it . It is not unusual, therefore, for S2 to clusters with non-ribosomal proteins. With approximately 60 amino acids, S14 and L29 are the smallest ribosomal proteins in the data set. The short sequence of these proteins influences their amino acid composition and both appear as outliers in figure 2. However, projection of the S14 protein onto the horizontal axis of figure 2 shows that, despite its short length, this protein has some characteristics of the majority of ribosomal proteins. S14 is completely devoid of any globular domain, and most of the protein has an extended coil structure . Although S14 is involved in intimate protein-protein interactions, almost its entire length is involved in RNA contacts and its arginine and lysine contents are similar to those of most ribosomal proteins . S14 is therefore indistinguishable from most ribosomal proteins in the x-axis projection of figure 2. On the other hand, L29 interacts with the L23 protein and with only one of the six domains of 23S rRNA . This characteristic, and its short length, may therefore explain the position of L29 in figure 2.
G+C content and optimal growth temperature are the two factors that most influence differences in amino acid composition between organisms. Analysis of the optimal temperatures of the enzymes extracted from hyperthermophilic organisms showed that thermal resistance was an intrinsic property of these enzymes . Comparative analysis of the amino acid composition of orthologous proteins from several mesophilic and thermophilic organisms indicated some amino acid substitutions that are preferred in thermophiles . However, the small number of sequences analyzed and the fact that factors other than temperature can affect the amino acid composition of proteins revealed the inconsistency of theses results . Comparison of the first completely sequenced genomes of several thermophiles and mesophiles showed that proteins from thermophiles contain higher levels of both charged and hydrophobic residues and lower levels of polar and uncharged ones . Once more complete genomes were sequenced, new analyses were performed using different methods and different datasets [8–10, 21–25]. Although these studies show several discrepancies in the role of each amino acid, there is a consensus that glutamate (E) and, to a lesser extent, valine (V) are the amino acids that are more represented in thermophiles than in mesophiles. These were also the amino acids that were most represented in thermophiles when our method was used.
There are greater discrepancies, however, over which amino acids are used with the lowest frequency in thermophiles or with the highest frequency in mesophiles. For example, Singer and Hickey  found that these amino acids were A, H, Q and T; Kreil and Ouzounis  found that they were Q and T; and Tekaia and coworkers  found only Q. These discrepancies indicate that hyperthermophilic and mesophilic enzymes may be very similar – their difference being that hyperthermophilic enzymes are more rigid than mesophilic enzymes . To increase their rigidity, hyperthermophilic enzymes may adopt several strategies but a common rule could be that more charged residues are found in hyperthermophilic proteins, mostly at the expense of uncharged polar residues . Computational, biochemical, and structural evidence now supports the hypothesis that ion pair formation, hydrogen bonds, and hydration, rather than hydrophobic interactions, play important roles in the stabilization of enzymes from extremophiles . Also, we cannot talk of a common amino acid usage in mesophiles because an adaptation to live at intermediate temperatures is unnecessary. When comparing the amino acid compositions of thermophilic and mesophilic proteins, therefore, different datasets and methods obtain different results.
The use of certain amino acids with higher or lower frequencies in thermophiles is important for the thermal stability of their enzymes. However, other factors may contribute to survival at high temperatures. Thermophilic archaea, for example, may be protected by their unique membrane lipids, the use of a reverse gyrase that introduces positive supercoils , a DNA repair system [28, 29] and the presence of special DNA-binding proteins . One of these thermophilic-specific proteins may be highly basic histone-like proteins that wind and compact DNA into a nucleosome-like structure and thus protect them from heat denaturation . Loss of some of these factors may lead to a lesser ability to grow at high temperatures. This could be the case of the Euryarcheota Halobacterium sp (Hbs) and M. acetivorans (Mac), two archaea whose optimal growth temperature is below 40°C but that cluster with other thermophilic species in figure 3. The amino acid compositions of these two Euryarchaeota, which are very similar to those of other thermophiles, may be a trace of their past ability to grow at high temperatures. A thermophile-specific NTPase found in 13 thermophilic genomes and absent in 52 mesophilic genomes is present in M. acetivorans . This suggests that M. acetivorans facultatively could be thermophilic . Although the phylogenetic position of these two archaea and our analysis of the amino acid composition suggest a recent transversion to mesophily in Halobacterium sp and M. acetivorans, this hypothesis is speculative and needs to be supported by stronger evidence. In this sense, it would be useful to identify proteins present in all thermophilic Euryarchaeota but not in mesophilic Euryarchaeota. One of these proteins is a dsDNA-binding protein called Alba (short for "acetylation lowers binding affinity"), which is present in several thermophilic archaea but not in Halobacterium sp or M. acetivorans . The correlation of Alba with growth at high temperatures hints at a role for Alba in DNA protection and stability under these conditions . Interestingly, it has been suggested that this protein constrains negative DNA supercoils in a temperature-dependent fashion, which suggests that it may function in chromosomal organization and accessibility .
The relationship between genomic G+C content and optimal growth temperature in prokaryotes has been debated recently in the literature [34–37]. Because G:C pairs in DNA are more thermally stable than A:T pairs, it has been suggested that a high G+C content may be a selective response to high temperature. In this sense, a significant correlation has been observed between optimal growth temperature and the G+C content of structural RNAs [35, 36]. When open reading frames are analyzed, some studies have concluded that there is no correlation between G+C content and optimal growth temperature [34–36] and others have found a positive correlation among some families of prokaryotes . If this correlation exists, it could be argued that the G+C-content dependence observed in the amino acid composition of prokaryotes is a consequence of their thermophily-dependence. In our dataset, the G+C content at the third codon position and the optimal growth temperature do not correlate significantly (r = 0.081). In addition, the results obtained with the tucker3 algorithm indicate that these two variables are independent. The amino acid variation associated with G+C content and optimal growth temperature corresponds to the second and fifth factor of the amino acid loadings matrix, respectively. Because the principal components obtained with the tucker3 model were constrained to be orthogonal, it can be concluded that the two factors are independent. The correlation observed therefore for the second and fifth factors of the amino acid loadings matrix is only 6.03E-5. Similar arguments can be applied to the second and third factors (those associated with differences in G+C content and optimal growth temperature, respectively) of the organism's loadings matrix. Moreover, the amino acid preferred by thermophiles is glutamate, which is an amino acid encoded by intermediate-GC content codons (GA [A,G]). All this evidence suggests that the amino acid variations related to variations of G+C content and optimal growth temperature are independent and that the observed G+C-dependence is not a consequence of a thermophily dependence.
We have shown that a multi-way method can be used to analyze differences in the amino acid composition within and between species. This method determines the relative influence of the various factors behind the heterogeneity of amino acid composition in proteins. Also, using a dataset consisting of a group of orthologous proteins present in all the species analyzed ensures that the differences in the amino acid composition between species are related to an intrinsic property of their proteins.
Definition of the data set
Sixty-four orthologous groups of proteins present in 62 archaea and bacteria were imported from the COG database . This data set was chosen because it represents a group of proteins that are present in all the genomes (except for the eukaryotic species and the archaea Sulfolobus solfataricus) in the August 2003 version of the COG database . The organisms analyzed [see Additional file 1] comprise various taxonomic groups of bacteria and archaea with different growth temperatures (from mesophiles to hyperthermophiles) and a wide range of G+C contents. The 64 groups of orthologous proteins analyzed [see Additional file 2] have essential functions such as ribosome formation, tRNA synthesis or translation initiation and are therefore conserved in all the organisms we analyzed. Protein sequences of each COG family for the 62 species analyzed were extracted using our own PERL programs. When paralogous sequences were detected, only the largest proteins were retained. In three cases a fusion protein had two of the activities analyzed. These proteins were the HP1198 and jhp1121 from H. pylori 26695 and H. pylori J99, respectively, which contain the two fused beta subunits of a DNA-directed RNA polymerase (COG0085 and COG0086), and the MTH39 protein from M. thermoautotrophicum, which contains fused the ribosomal proteins L13 and S9 (COG0102 and COG0103). These fusion proteins were computationally cut into two proteins using information from the comparison of orthologous sequences. So as not to be predicted as horizontally transferred, all the protein sequence in the final dataset was checked using the HGT-Database [39, 40]. Using our own PERL programs, we calculated the amino acid composition of each protein of the 64 orthologous groups from the 62 previously defined species. Our data set was therefore a three-dimensional matrix of order (20 × 64 × 62), in which the rows represented the frequency of use of each amino acid, the columns represented the group of orthologous proteins, and the third dimension represented the 62 species analyzed. To analyze the amino acid variability within and between species, the three-dimensional matrix was column mean-centered and analyzed using a multi-way method called Tucker3  developed for the MATLAB (The Mathworks, Natick, MA, USA) environment.
Analysis of the 3-way data array using the Tucker3 method
In standard multivariate data analysis, data are arranged in a two-way structure (i.e. a table or a matrix) that contains objects and variables. These tables are analyzed with a method such as PCA or correspondence analysis, which enables large amounts of data to be condensed to a few representative variables (called principal components or factors). The projections of objects and variables in the representative principal factors are called scores and loadings matrices, respectively. These matrices are used to identify the similarities and dissimilarities between objects and variables in the data under investigation. In this sense, the scores and loadings plots of different principal components are indispensable for finding patterns and clusters between variables and identifying the variables responsible for the formation of these clusters. They are also indispensable for identifying the variables, or combinations of variables, responsible for the maximum variability between objects. Sometimes, however, the structure of the data set is such that a standard two-way table is not enough to describe it. In these cases, a third mode needs to be added to represent the data set, which can be imagined as a parallelepiped of size r1 × r2 × r3. To apply standard PCA, these three-way data arrays must be unfolded to obtain a two-way data table. This can be done in several ways, depending on what one would like to focus on. Multi-way analysis is the natural extension of multivariate analysis, in which data are arranged in three or more dimensions. There are several multi-way models for analyzing three-way data sets. Here we focus on the Tucker3 model. The origin of this method lies in psychometrics and the pioneering work of Tucker . The algorithm solution for estimating the model was later considerably improved by Kroonenberg and de Leeuw .
Figure 4 shows the basis of the Tucker3 model, which takes as a starting point a 3-way data array X of order (20 × 64 × 62). The tucker3 algorithm decomposes the data array into a matrix E of residuals, a 3-way core array G of order (w1 × w2 × w3), and three component matrices A, B, C called loadings matrices. In order to use less computation time and make the results easier to interpret, the component matrices of each loadings matrix are usually constrained to be orthogonal. Mathematically, it can be said that
where a il , b jm and c kn denote elements of the component matrices A, B and C of orders 20 × w1, 64 × w2 and 62 × w3 respectively, g lmn denotes the elements (l, m, n) of the w1 × w2 × w3 core array G, and e ijk denotes the error term for element x ijk and is an element of the 20 × 64 × 62 array E. See  for more details on Tucker's models.
The number of factors in each of the three modes, i.e. w1, w2 and w3, is determined by the analyst from a priori knowledge of the data or by evaluating models with different combinations and choosing the order that provides the most accurate model of X . The elements of the core array define how individual loadings vectors in the different modes interact. For example, the [1,1,1] factor combination means the first component of the amino acid composition, the protein function and the organism modes, respectively. The core array therefore provides a way to interpret the solutions since its squared entries represent the relative importance of each individual factor combination in terms of explained variability. A highly effective way to analyze these combinations of factors is to represent them in a series of bi- or triplots. In our case, triplots consist of the superposition of the representation of two of the principal components for each loadings matrix. These plots can be interpreted as the plot of scores and loadings in standard PCA. By analysing points of the same type (amino acids, functions or organisms), we can analyse the variation observed and the clusters formed. Also, by analysing the mutual positions of points of different types, we can identify which variables are responsible for the variation. Though several successful applications of multi-way models have been demonstrated in quite different areas [43, 44], these models have not received much attention in the analysis of proteomic or genomic data. In this paper we show that this kind of algorithm is useful for comparing the amino acid composition of several groups of orthologous proteins.
Robustness of the model and reliability of the results
Robustness measures the ability of a method to remain unaffected by small variations in variables or procedure. To validate the robustness of our model and check that the results were not a product of chance, we have verified that minor changes in the modelling procedure do not affect our conclusions. Specifically we have verified that the model: (i) was unaffected by small changes in the structural model, and (ii) was unaffected by slight changes such as the elimination at random of some organisms or group of orthologous sequences.
Another important matter is the reliability of the results. Our results have been shown to be reliable by the percentage of summed variance fraction (98.60%) correlated with some biological character, the high value of these correlations, and the consistency of our results with previous studies and our knowledge of the system analyzed.
Rispe C, Delmotte F, van Ham RCHJ, Moya A: Mutational and selective pressures on codon and amino acid usage in Buchnera , endosymbiotic bacteria of aphids. Genome Res 2004, 14: 44–53. 10.1101/gr.1358104
Mackiewicz P, Gierlik A, Kowalczuk M, Dudek MR, Cebrat S: How does replication-associated mutational pressure influence amino acid composition of proteins? Genome Res 1999, 9: 409–416.
Rocha EPC, Danchin A, Viari A: Universal replication biases in bacteria. Mol Microbiol 1999, 32: 11–16. 10.1046/j.1365-2958.1999.01334.x
Krogh A, Larsson B, von Heijne G, Sonnhammer ELL: Predicting transmembrane protein topology with a hidden markov model: application to complete genomes. J Mol Biol 2001, 305: 567–580. 10.1006/jmbi.2000.4315
Fujiwara Y, Asogawa M: Prediction of subcellular localizations using amino acid composition and order. Genome Informatics 2001, 12: 103–112.
Lin K, Kuang Y, Joseph JS, Kolatkar PR: Conserved codon composition of ribosomal protein coding genes in Escherichia coli, Mycobacterium tuberculosis and Saccharomyces cerevisiae: lessons from supervised machine learning in functional genomics. Nucleic Acids Res 2002, 30: 2599–2607. 10.1093/nar/30.11.2599
Singer GAC, Hickey DA: Nucleotide bias causes a genomewide bias in the amino acid composition of proteins. Mol Biol Evol 2000, 17: 1581–1588.
Kreil DP, Ouzounis CA: Identification of thermophilic species by the amino acid compositions deduced from their genomes. Nucleic Acids Res 2001, 29: 1608–1615. 10.1093/nar/29.7.1608
Tekaia F, Yeramian E, Dujon B: Amino acid composition of genomes, lifestyles of organisms, and evolutionary trends: a global picture with correspondence analysis. Gene 2002, 297: 51–60. 10.1016/S0378-1119(02)00871-5
Lobry JR, Chessel D: Internal correspondence analysis of codon and amino-acid usage in thermophilic bacteria. J Appl Genet 2003, 44: 235–261.
Pe'er I, Felder CE, Man O, Silman I, Sussman JL, Beckmann JS: Proteomic signatures: amino acid and oligopeptide compositions differentiate among phyla. Proteins 2004, 54: 20–40. 10.1002/prot.10559
Lynn D, Singer GAC, Hickey DA: Synonymous codon usage is subject to selection in thermophilic bacteria. Nucleic Acids Res 2002, 30: 4272–4277. 10.1093/nar/gkf546
Gu X, Hewett-Emmett D, Li W-H: Directional mutational pressure affects the amino acid composition and hydrophobicity of proteins in bacteria. Genetica 1998, 102/103: 383–391. 10.1023/A:1017028102013
Jordan IK, Kondrashov FA, Adzhubei IA, Wolf YI, Koonin EV, Kondrashov AS, Sunyaev S: A universal trend of amino acid gain and loss in protein evolution. Nature 2005, 433: 633–638. 10.1038/nature03306
Trifonov EN: The triplet code from first principles. J Biomol Struct Dyn 2004, 22: 1–11.
Brodersen DE, Clemons WM, Carter AP, Wimberly BT, Ramakrishnan V: Crystal structure of the 30 S Ribosomal subunit from Thermus thermophilus : Structure of the proteins and their interactions with 16 S RNA. J Mol Biol 2002, 316: 725–768. 10.1006/jmbi.2001.5359
Ban N, Nissen P, Hansen J, Moore PB, Steitz TA: The complete atomic structure of the large ribosomal subunit at 2.4 A Resolution. Science 2000, 289: 905–920. 10.1126/science.289.5481.905
Vieille C, Zeikus GY: Hyperthermophilic enzymes: sources, uses and molecular mechanisms for thermostability. Microbiol. Mol Biol Rev 2001, 65: 1–43. 10.1128/MMBR.65.1.1-43.2001
Bohm G, Jaenicke R: Relevance of sequence statistics for the properties of extremophilic proteins. Int J Pept Protein Res 1994, 43: 97–106.
Deckert G, Warren PV, Gaasterland T, Young WG, Lenox AL, Graham DE, Overbeek R, Snead MA, Keller M, Aujay M, Huber R, Feldman RA, Short JM, Olsen GJ, Swanson RV: The complete genome of the hyperthermophilic bacterium Aquifex aeolicus . Nature 1998, 392: 353–358. 10.1038/32831
Cambillau C, Claverie J-M: Structural and genomic correlates of hyperthermostability. J Biol Chem 2000, 275: 32383–32386. 10.1074/jbc.C000497200
Farias ST, Bonato MCM: Preferred amino acids and thermostability. Genetics Mol Res 2003, 2: 383–393.
Nakashima H, Fukuchi S, Nishikawa K: Compositional changes in RNA, DNA and proteins for bacterial adaptation to higher and lower temperatures. J Biochem 2003, 133: 507–513. 10.1093/jb/mvg067
Saunders NFW, Thomas T, Curmi PM, Mattick JS, Kuczek E, Slade R, Davis J, Franzmann PD, Boone D, Rusterholtz K, Feldman R, Gates C, Bench S, Sowers K, Kadner K, Aerts A, Dehal P, Detter C, Glavina T, Lucas S, Richardson P, Larimer F, Hauser L, Land M, Cavicchioli R: Mechanisms of thermal adaptation revealed from the genomes of the Antarctic archaea Methanogenium frigidum and Methanococcoides burtonii . Genome Res 2003, 13: 1580–1588. 10.1101/gr.1180903
Singer GAC, Hickey DA: Thermophilic prokaryotes have characteristic patterns of codon usage, amino acid composition and nucleotide content. Gene 2003, 317: 39–47. 10.1016/S0378-1119(03)00660-7
Rossi M, Ciaramella M, Cannio R, Pisani FM, Moracci M, Bartolucci S: Extremophiles 2002. J Bacteriol 2003, 185: 3683–3689. 10.1128/JB.185.13.3683-3689.2003
Forterre P: A hot story from comparative genomics: reverse gyrase is the only hyperthermophile-specific protein. Trends Genet 2003, 18: 236–237. 10.1016/S0168-9525(02)02650-1
Makarova KS, Aravind L, Grishin NV, Rogozin IB, Koonin EV: A DNA repair system specific for thermophilic Archaea and bacteria predicted by genomic context analysis. Nucleic Acids Res 2002, 30: 482–496. 10.1093/nar/30.2.482
Guy CP, Majernik AI, Chong JPJ, Bolt EL: A novel nuclease-ATPase (Nar71) from archaea is part of a proposed thermophilic DNA repair system. Nucleic Acids Res 2004, 32: 6176–6186. 10.1093/nar/gkh960
Klinger C, Robbach M, Howe R, Kaufmann M: Thermophile-specific proteins: the gene product of aq_1292 from Aquifex aeolicus is an NTPase. BMC Biochemistry 2003, 4: 12. 10.1186/1471-2091-4-12
White MF: Archaeal DNA repair: paradigms and puzzles. Biochemical Society Transactions 2003, 31: 690–693. 10.1042/BST0310690.
Jelinska C, Conroy MJ, Craven CJ, Hounslow AM, bullough PA, Waltho JP, Taylor GL, White MF: Obligate heterodimerization of the archaeal Alba2 protein with Alba1 provides a mechanism for control of DNA packaging. Structure 2005, 13: 963–971. 10.1016/j.str.2005.04.016
Xue H, Guo R, Wen Y, Liu D, Huang L: An abundant DNA binding protein from the hyperthermophilic archaeon Sulfolobus shibatae affects DNA supercoiling in a temperature-dependent fashion. J Bacteriol 2000, 182: 3929–3933. 10.1128/JB.182.14.3929-3933.2000
Wang H-C, Susko E, Roger AJ: On the correlation between genomic G+C content and optimal growth temperature in prokaryotes: Data quality and confounding factors. Biochem Biophys Res Commun 2006, 342: 681–684. 10.1016/j.bbrc.2006.02.037
Galtier N, Lobry JR: Relationships between genomic G+C content, RNA secondary structures, and optimal growth temperature in prokaryotes. J Mol Evol 1997, 44: 632–636. 10.1007/PL00006186
Hurst LD, Merchant AR: High guanine-cytosine content is not an adaptation to high temperature: a comparative analysis amongst prokaryotes. Proc R Soc Lond B 2001, 268: 493–497. 10.1098/rspb.2000.1397
Musto H, Naya H, Zavala A, Romero H, Alvarez-Valin F, Bernardi G: Correlations between genomic GC levels and optimal growth temperatures in prokaryotes. FEBS Lett 2004, 573: 73–77. 10.1016/j.febslet.2004.07.056
Tatusov RL, Koonin EV, Lipman DJ: A genomic perspective on protein families. Science 1997, 278: 631–637. 10.1126/science.278.5338.631
Garcia-Vallve S, Romeu A, Palau J: Horizontal gene transfer in bacterial and archaeal complete genomes. Genome Res 2000, 10: 1719–1725. 10.1101/gr.130000
Garcia-Vallve S, Guzman E, Montero MA, Romeu A: HGT-DB: a database of putative horizontally transferred genes in prokaryotic complete genomes. Nucleic Acids Res 2003, 31: 187–189. 10.1093/nar/gkg004
Kroonenberg PM, de Leeuw J: Principal component analysis of three-mode data by means of alternating leats squares algorithms. Psychometrika 1980, 45: 69–97. 10.1007/BF02293599
Tucker L: Some mathematical notes on three-mode factor analysis. Psychometrika 1966, 31: 279–311. 10.1007/BF02289464
Henrion R: N-way principal component analysis. Theory, algorithms and applications. Chemom Intell Lab Syst 1994, 25: 1–23. 10.1016/0169-7439(93)E0086-J
Andersson CA, Munck L, Henrion R, Henrion G: Analysis of N-dimensional data arrays from fluorescence spectroscopy of an intermediary sugar product. Fresenius J Anal Chem 1997, 359: 138–142. 10.1007/s002160050550
This work has been financed by the following projects: "Ministerio de Ciencia y Tecnología" (Ref. BIO2003-07672), "Fundación BBVA, I Convocatoria de Ayudas a la Investigación en Bioinformática (Ref.:BIO 04)" and "FIS (Ref. C03/O8)". We thank Kevin Costello of the Language Service of the Rovira i Virgili University for his help with writing the manuscript and Ricard Boque from the Chemometrics, Qualimetrics and Nanosensors Group for his comments on the manuscript.
AP helped to design the study, carried out all the analyses and helped to draft the manuscript. SG-V conceived the study, collected the sequences of the orthologous groups of proteins, supervised the work and wrote the paper. Both authors read and approved the final manuscript.
Electronic supplementary material
Additional File 1: Optimal growth temperatures, G+C content and position in the x and y axes of figure 3 of the organisms analyzed. Thermophilic organisms (defined as species with an optimal growth temperature above 60°C) are shown in red. (PDF 135 KB)
Additional File 2: COG family, function and position on the x and y axis of figure 2 of the genes common to all species analyzed. (PDF 117 KB)
About this article
Cite this article
Pasamontes, A., Garcia-Vallve, S. Use of a multi-way method to analyze the amino acid composition of a conserved group of orthologous proteins in prokaryotes. BMC Bioinformatics 7, 257 (2006). https://doi.org/10.1186/1471-2105-7-257
- Ribosomal Protein
- Amino Acid Composition
- Codon Position
- Factor Combination
- Optimal Growth Temperature