Evolution of complete proteomes: guanine-cytosine pressure, phylogeny and environmental influences blend the proteomic architecture
© Chen et al.; licensee BioMed Central Ltd. 2013
Received: 21 March 2013
Accepted: 1 October 2013
Published: 3 October 2013
Guanine-cytosine (GC) composition is an important feature of genomes. Likewise, amino acid composition is a distinct, but less valued, feature of proteomes. A major concern is that it is not clear what valuable information can be acquired from amino acid composition data. To address this concern, in-depth analyses of the amino acid composition of the complete proteomes from 63 archaea, 270 bacteria, and 128 eukaryotes were performed.
Principal component analysis of the amino acid matrices showed that the main contributors to proteomic architecture were genomic GC variation, phylogeny, and environmental influences. GC pressure drove positive selection on Ala, Arg, Gly, Pro, Trp, and Val, and adverse selection on Asn, Lys, Ile, Phe, and Tyr. The physico-chemical framework of the complete proteomes withstood GC pressure by frequency complementation of GC-dependent amino acid pairs with similar physico-chemical properties. Gln, His, Ser, and Val were responsible for phylogeny and their constituted components could differentiate archaea, bacteria, and eukaryotes. Environmental niche was also a significant factor in determining proteomic architecture, especially for archaea for which the main amino acids were Cys, Leu, and Thr. In archaea, hyperthermophiles, acidophiles, mesophiles, psychrophiles, and halophiles gathered successively along the environment-based principal component. Concordance between proteomic architecture and the genetic code was also related closely to genomic GC content, phylogeny, and lifestyles.
Large-scale analyses of the complete proteomes of a wide range of organisms suggested that amino acid composition retained the trace of GC variation, phylogeny, and environmental influences during evolution. The findings from this study will help in the development of a global understanding of proteome evolution, and even biological evolution.
KeywordsProteome Evolution GC pressure Phylogeny Environmental influence
The dramatic increase of genomic data in public databases for a wide variety of species as the result of plummeting costs and newly introduced sequencing technologies, has set off a wave of post-genomic research [1, 2]. One major goal of the post-genome analyses is to characterize features such as composition, content, and functionality of genomes and/or proteomes and to decipher their connection to already known biological characteristics [1, 3]. Guanine-cytosine (GC) content is a well-known feature of genomes, similarly, amino acid composition is the distinct feature of proteomes . However, the many in silico studies on proteomes have focused mainly on the functional annotation of individual proteins in a proteome . Comparatively, the distribution of amino acids on the proteomic scale has received far less attention. Thus, it has been suggested that information about the amino acid compositions of proteomes should be provided in databases, as is now done routinely for the GC content of genomes . Proteomic architecture carries striking signatures of biological traits, and its compositional biases are closely associated with already recognized internal factors, such as GC variation, and emerging external factors, such as species phylogeny and environmental niche [4–7]. Elucidating the associations between proteomic architecture and such factors is essential to achieve a better understanding of the evolution of proteomes and even organisms.
After the seminal studies of Sueoka in 1961 , many different studies drew the recurring conclusion that amino acid composition, whether at the protein or proteomic level, was essentially under the influence of GC pressure [3, 9–13]. It is neither new nor surprising that the use of amino acids with GC-rich codons such as Ala, Arg, Gly, and Pro will increase with increasing genomic GC content, while the use of amino acids with AT-rich codons such as Asn, Ile, Lys, Phe, and Tyr will decrease [10, 12–14]. Therefore, it is believed that mutational pressure on DNA composition is a very powerful and pervasive force in long-term protein evolution . However, what response mechanisms the organisms have to counter the potentially deleterious effects of this nucleotide bias-induced change on protein function is still unclear, mainly because the effect of nucleotide bias on amino acid composition is very large and widespread. One hypothesis was that the resulting amino acid sequence changes are nonrandom because the mutational bias is strongly directional, and are not caused by natural selection acting directly on protein function. Consequently, the evolutionary dynamics of these changes cannot be described in terms of either Darwinian selection or random genetic drift . Another subsequent explanation suggested compensatory changes were made to reduce the impact of amino acid bias on protein structure and function. Previous analyses have shown that amino acid substitutions such as Arg for Lys and Ala for Ser are important for the stability of proteins in some species .
Studies of interlinkages between species phylogeny and proteomic architecture have been less conducted, primarily due to the restriction of information capacity for 20-dimensional amino acid data, far enough to fully characterize the species diversity. Nevertheless, several previous works have implied a close connection between species phylogeny and the amino acid composition information of their proteomes. For example, an in-depth comparison of proteomic composition across more than 70 genera from all three domains revealed a clear differentiation of archaea, eubacteria, and eukaryote from each other . An extensive analysis of amino acid usage in 208 proteomes showed a clearcut segregation of eukaryotes from prokaryotes . The phylogenetic groups within 1029 bacterial and archaeal strains, generally, could be discriminated from other groups based on the amino acid composition of their proteomes . However, it should be stressed that proteomic amino acid information is limited by capacity, and is capable of discrimination of species only at a general level of phylogenetic classification, such as at the domain level, rather than at a detailed classification level.
It is an emerging view that proteomic amino acid composition is influenced by external factors, namely environmental influences . Throughout evolution, organisms may have used this to increase their fitness to an environmental niche, especially under conditions such as extreme temperature, pressure, pH, and salinity. Regardless of the environmental niche, however, adaptation and maintenance of protein integrity and function seem to be fundamental to the survival of entire organisms [3, 5, 17–19]. Accordingly, many studies have attempted to relate habitats of organisms to the amino acid composition of their proteins and proteomes. At the habitat level, many studies have investigated the relationship between extreme environments, especially at extremely high temperature and salinity ranges, and the adaptive proteomes to identify amino acid signatures associated with different lifestyles [19–25]. At the phylogeny level, environmental signatures have been studied mostly in the proteomes of unicellular prokaryotes. However, a few studies on higher multicellular eukaryotes have also found that the proteomes of these organisms may be under evolutionary adaptation to their habitats. For example, proteomic comparisons among eleven endothermic and ectothermic vertebrates suggested that the amino acid composition in endothermic vertebrates was biased in the same direction as the composition in the proteomes of thermophilic prokaryotes . Analysis of the proteomes across streptophyte lineages, including charophycean algae and embryophytic plants, showed amino acid compositional shifts during streptophyte transitions to terrestrial habitats . However, previous results from direct comparisons of proteomes associated with different lifestyles can be partly misleading. For example, comparisons among temperature-sensitive species may exaggerate the influence of temperature on amino acid composition and thus mask underlying interactions between environmental niche and proteomic architecture.
The compatibility between the genetic code and amino acid composition of proteomes is also worthy of further study, because it may provide some clues about their co- adaptation. Based only on theoretical probability, amino acid frequency should coincide with the expected random distribution of the 61 sense codons of the universal genetic code. However, the observed frequencies of some amino acids show obvious biases against the expected rates. Most of the deviations can be interpreted as the consequences of selective pressures at the level of nucleotide composition or protein structure and function [7, 11]. At the genetic level, as mentioned above, variations in the nucleotide composition could shift amino acid composition noticeably. Moreover, some nucleotide combinations, such as the dinucleotide cystosine-guanine (CG), are suppressed in many genomes, which may lead to suppression of the codons that contain the suppressed nucleotides and hence to an underrepresentation of the corresponding amino acids [11, 28–31]. At the protein level, on the one hand, the frequency bias of some amino acids can be attributed to structural constraints; for example, Met is principally the initiation amino acid and thus its presence does not follow a random distribution. On the other hand, some amino acids like Cys and Pro occur at significantly less than the expected rates because of their strong influences on protein tertiary structure; that is, Cys can form disulphide bridges, and Pro can terminate helices and introduce turns . Moreover, other factors such as biosynthetic costs, metabolic efficiency, and elements of amino acids also result in their compositional bias [7, 32–35].
Unlike previous studies that compared the proteomes of unicellular prokaryotes with habitat tags, a large number of eukaryotes covering a wide taxonomy range were considered in the present study. First, the complete proteome sequences from 63 archaea, 270 bacteria, and 128 eukaryotes in public databases were surveyed for a comprehensive comparison and analysis of their organization, similarity, uniqueness, and variability at the proteomic level to help understand long-term changes in amino acid compositions, selective constraints, and pressures across the proteomes. Then, a global picture was created to help visualize how the amino acid composition of complete proteomes could have been blended mainly by the joint actions of three forces, namely GC pressure, phylogeny, and environmental influences. This study has provided new insights into how proteomes change over evolutionary time and how they are used as the blueprints for organisms.
Frequency distribution of amino acids in archaea, bacteria, and eukaryota
The differences and biases in amino acid frequency were studied among the three domains. Translation from the genetic material (DNA or mRNA sequences) to proteins (amino acid sequences) is based on the triplet genetic code ; therefore, based merely on probability distribution, the genomic amino acid frequency should be the same as the frequency of the synonymous codons in all the codons in a genome. Therefore, we assumed that there was no difference between the genomic amino acid frequency and the corresponding synonymous codon frequency. The comparison showed differences in the amino acid frequencies among the three domains. The frequencies of Asp, Glu, Gly, Ile, Tyr and Val were higher in Archaea than in Bacteria and Eukaryota, while, the dominant amino acids were Ala, Leu, and Phe in Bacteria, and Asn, Cys, Gln, His, Pro, and Ser in Eukaryota. It is worth mentioning that there seemed to be no frequency difference for Lys among the three domains. We found that the amino acid frequencies differed from the corresponding synonymous codon frequencies in the three domains and two distinguishing groups were found. One group, the “up group”, contained Ala, Asn, Asp, Glu, Lys, Met, and Phe because the frequencies of the amino acids in this group were significantly higher than their assumed frequencies based on the synonymous codon frequencies. The other group, the “down group”, included Arg, Cys, His, Pro, Ser, Thr, and Trp because the frequencies of the amino acids in this group were significantly lower than their assumed frequencies based on the synonymous codon frequencies.
Physico-chemical profile of amino acids in archaea, bacteria, and eukaryota
Based on the chemical structure of the side chains, the amino acids can be classified into aliphatic, aromatic, and heterocyclic amino acids. Our analysis showed that the aliphatic amino acids made up the largest proportion of the amino acids in all three domains (about 85%), although, of the three, Archaea had the highest proportion (Figure 3b). Bacteria and Eukaryota had the highest frequencies of aromatic and heterocyclic amino acids, respectively.
Interlinkage of amino acid frequencies in the tested proteomes
When the correlations between the frequencies of the 20 amino acid were studied, we identified two distinct amino acid groups; one that contained Ala, Arg, Gly, Pro, Trp, and Val, and another that included Asn, Lys, Ile, Phe, and Tyr. The frequencies of the amino acids were positively correlated within each group, but negatively correlated between the two groups (Additional file 1). These associations were consistent in all three domains.
Principal components of amino acid frequency matrix from the tested proteomes
Principal components of the amino acid frequency in the 20-dimensional matrices for the 461 tested species were analyzed, and the top three components, which together accounted for 71.81% of the total inertia, were listed (Additional file 2). The data showed that principal component 1 (Prin1), which made up the largest proportion, was positively and closely related to the frequencies of Ala, Arg, Gly, Pro, Trp, and Val, but negatively related to the frequencies of Asn, Ile, Lys, Phe, and Tyr. These amino acids had already been found to be correlated as mentioned in the previous section. Therefore, Prin1 was focused mainly on changes among the 11 frequency-dependent amino acids. On the other hand, the second principal component (Prin2) mainly reflected changes among Gln, His, Ser, and Val, and the third principal component (Prin3) mainly focused on changes among Cys, Leu, and Thr.
Cluster analysis of 461 species based on amino acid frequency differences
The distribution characteristics of the amino acids were compared among the 11 groups (Additional file 3). The group A species were from the Bacteria and Eukaryota domains, and the main members were Actinobacteria, plants, and protists. The proteomes of the species in this group contained the highest frequencies of Ala and Cys, and the lowest frequencies of Phe, Ile, and Tyr, and the aromatic amino acids. Group B species were all Eukaryota, covering all the tested animals and fungi, and some of the plants and protists. The proteomes of the species in this group had the highest frequencies of His, Pro, and Ser, and the lowest frequencies of Gly and Val. The group B proteomes shared a distinctive amino acid framework; for example, the proteomes contained the highest percentages of polar and uncharged amino acids and the lowest percentages of hydrophobic amino acids, and the most heterocyclic amino acids and the least aliphatic amino acids, among the 11 groups. The group C species were mainly Thermoprotei and Methanomicrobia in the Archaea domain. The proteomes of the species in this group had the lowest frequencies of Asp, His, and Gln and no obviously dominant amino acids. The group D species were all Archaea and all were Halobacteria. Interestingly, the Halobacteria proteomes showed a distinct distribution in amino acid frequency with Asp, Glu, Gly, Thr, and Val being the amino acids with highest frequencies and with many other amino acids, such as Lys, Leu, Met, and Asn, with low frequencies compared with other groups. Furthermore, the group D proteomes contained the largest percentages of charged and aliphatic amino acids among the 11 groups. The group E species were all Bacteria, and mainly from the Alphaproteobacteria class. The proteomes of the species in this group had the lowest frequency of Cys among the 11 groups. The group F was a large branch of Bacteria species mainly from the Actinobacteria and Proteobacteria. Arg was the dominant amino acid proteomes of the species in group F, and Glu and Ser were least abundant. Moreover, the group F proteomes had the largest portion of hydrophobic amino acids, but the smallest portion of polar and uncharged amino acids among the 11 groups. The groups G, H and I were three small branches of Bacteria species, mainly from Gammaproteobacteria (group G), Cyanobacteria (group H), and Firmicutes and Proteobacteria (group I). The dominant amino acid in the group G proteomes was Met, while, in the group H, Leu, Gln, and Trp were the dominant amino acids. The group H proteomes also had the lowest composition of charged amino acids among the 11 groups. Amino acid frequencies in the proteomes of the species in group I showed a balanced distribution, implying no obvious frequency preference. The group J and K branches contained species from both Archaea and Bacteria. In group J, Euryarchaeota and Firmicutes were the dominant species, and in group K Bacteroidetes/Chlorobi, Euryarchaeota, and Firmicutes were dominant. The group J proteomes had the lowest frequency of Thr among the 11 groups. The amino acid distribution in the group K proteomes was distinctive; Asn, Ile, Lys, Phe, and Tyr were dominant, while Ala, Arg, Pro and Trp, had the lowest frequencies among the 11 groups. Furthermore, the group K proteomes contained the highest percentages of aromatic amino acids and the lowest percentages of heterocyclic amino acids.
When we examined the correlations between the 20 amino acid frequencies and their corresponding synonymous codon frequencies in the 11 groups (Additional file 3), we found that the correlation was highest in group B, which contained only Eukaryota. In groups E, F, G, H, and I, all of which contained only Bacteria, the correlations were also high, while in group D, which contained only Archaea, the correlation was relatively low. In the other groups, which contained species from different domains, the correlation coefficients showed a divergent distribution. It is worth mentioning that, in group K, the overall amino acid frequencies appeared to be the least correlated with the synonymous codon frequencies among the 11 groups.
Correlation between amino acid frequencies and their synonymous codon frequencies
Although the origin of the standard genetic code with its 20 amino acids has long been debated, it may be that the genetic code and amino acid frequencies have evolved to approach an optimal codon/amino acid relation; in other words, the genetic code may be selected to match the average optimal concentration of the amino acids [35, 37–39]. Indeed, there is a clear tendency for amino acids that are encoded by more codons to be more frequent in proteomes, although, as we have shown, the strength of this relationship, as evidenced by the Pearson coefficients between amino acid frequencies and their synonymous codon frequencies, varied considerably between species (Figure 2).
For individual amino acids, the concordance was not perfect. We found that most of the individual amino acids showed an obvious bias against the rates expected from a random distribution of the 61 sense codons (Figure 1). The up group of Ala, Asn, Asp, Glu, Lys, Met, and Phe showed significantly higher frequencies than their expected rates and most of the amino acids in this group have relatively small numbers of encoding codons. The down group of Arg, Cys, His, Pro, Ser, Thr, and Trp had significantly lower frequencies than their expected frequencies and these amino acids have relatively high numbers of encoding codons. Some of the deviation from expected frequency rates can be explained as a consequence of selection pressure at the level of protein structure and function. For example, Met is principally an initiation amino acid and, therefore, its presence does not follow a random distribution, explaining why it occurs at higher than expected rates. Both Cys and Pro, strongly influence protein tertiary structure, Cys by forming disulphide bridges, and Pro by terminating helices and introducing turns . If they occurred in higher proportions, protein misfolding may result. In addition, some of this variation from expected frequency rates can be explained as the result of selection pressure at the level of the nucleic acid sequences. For example, Arg is coded by six triplet codons, four of which begin with the dinucleotide CG, which is known to be suppressed in many classes of organisms such as vertebrates, and in many diverse species from protist, dicot plants and bacteria [11, 29]. This characteristic may lead to the suppression of codons for Arg and, hence, to its underrepresentation. In their study, Bharanidharan et al.  reported that the variations of Arg, Val, Asp, Glu, Ser, and Cys from predicted values were attributable to selection pressure at the amino acid level, while the variations of Thr, Phe, Lys, and Asn arose only from mutation and selection pressure at the level of the nucleic acid sequences.
For some amino acids, however, the variations from the expected frequency cannot be easily explained by these reasons. In our view, the deviation of some amino acids may not have arisen from themselves, but from strong preferences for other amino acids. In other words, when some amino acids are favored or rejected, this bias is bound to affect the frequencies of the other amino acids, even with no bias from them.
Effect of GC pressure on the amino acid compositions and on the framework of amino acid properties in the proteomes
Surprisingly, not only was amino acid composition affected by nucleotide bias, but the effect was large and widespread among the three domains. The genomic GC content across all the tested species ranged from 25 to 75%. For each 10% increase in genomic GC content, the use of Ala, Arg, Gly, Pro, Trp, and Val increased by 2.05, 1.04, 0.85, 0.52, 0.14, and 0.59%, respectively, while the use of Asn, Lys, Ile, Phe, and Tyr decreased by 1.02, 1.62, 1.16, 0.41 and 0.47%, respectively.
Next, we performed a correlation analysis between the distribution of amino acid properties and the genomic GC content across all the species to find whether GC variation would affect the framework of the amino acid properties (Additional file 4). We found that genomic GC content showed a strong negative correlation with the aromatic amino acids and a strong positive correlation with the hydrophobic amino acids. These correlations can be explained by the amino acids in both these groups. In the aromatic group, the main contributors, Phe and Tyr, were negatively correlated with the GC content, while in the hydrophobic group, the main contributors, Ala, Pro, Trp, and Val, were positively correlated with the GC content (Figure 6). However, we also observed that the impact of genomic GC content on the global physico-chemical trends of the amino acids seemed to be limited. For each 10% increase in genomic GC content, the proportion of aromatic amino acids decreased by 0.74%, while the proportion of hydrophobic amino acids increased by 1.88%. Thus, the change of the physico-chemical groups with genomic GC content was relatively low compared with that of the individual amino acids. In addition, the genomic GC variation showed no significant correlations with other important physico-chemical groups (Additional file 4). Therefore, we hypothesized that the framework of the amino acid properties in the proteomes was in a relatively stable state, and, at least, able to withstand GC pressure. Accordingly, two challenging questions arose. One, is the relatively stable physico-chemical framework in the proteomes the material basis of species in response to the pressure for survival? And two, can the use of amino acids be substituted by amino acids with similar physico-chemical properties? Further research on the roles that the physico-chemical properties play in the protein structure and function is needed to address the first question. As for the second question, despite our proposal that genomic GC variation largely affects the use of many of the amino acids, the physico-chemical framework was relatively independent of genomic GC content because of the frequency complementation between the GC-dependent amino acid pairs in the physico-chemical groups. For example, the frequency of amino acid Arg was positively correlated with genomic GC content while the frequency of Lys was negatively correlated (Figure 6); therefore, the GC-dependent amino acids Arg and Lys are complementary pairs in charged group and probably, the use of Arg and Lys was interchangeable to a certain degree. This finding might provide some evolutionary clues on amino acid substitutions. Indeed, it is worth exploring the biological significance of a stable physico-chemical framework of amino acids by the complementation of GC-dependent amino acids.
Principal components of amino acid matrices may represent GC variation, domain discrimination, and environmental influences
To extract information about the associations between the amino acid features and the corresponding observations from the data matrices, we used principal component analysis based on the projection of the 20 dimensional frequency vectors of the 20 standard amino acids onto their salient aspects. The most significant principal component (Prin1) in the analysis accounted for 43.91% of the total information in the matrices. In Prin1, the GC-dependent amino acids, Ala, Arg, Gly, Pro, Trp, and Val, were concentrated at high values, whereas Ala, Arg, Gly, Pro, Trp, and Val displayed low values (Figure 6 and Additional file 2). Prin1 showed a highly significant correlation to the genomic GC content (R = 0.9394, p < 0.0001), indicating that, to a large extent, genomic GC content determined amino acid usage. This factorial axis was assigned as the “GC-axis”. This finding was similar to earlier observations for numerous genomes that GC content was the most important determinant of global amino acid composition despite differences in information volume [3, 16, 17, 24]. The observed consistency of the impact of genomic GC content on patterns of amino acid composition across all three domains supports the idea that genomic GC content may be a driving force in genome evolution .
In contrast, Prin2, which represented 16.84% of the total information, separated the Archaea, Bacteria, and Eukaryota domains (Figure 4). Accordingly, this factorial axes was responsible mainly for the domain discrimination. The main amino acids were Gln, His, Ser, and Val. In previous studies, the second factorial axis was assimilated to a temperature axis because it could distinguish temperature-sensitive species [3, 4, 24]. Interestingly, in previous studies, these four amino acids were identified as amino acid signatures of hyperthermophiles, indicating that Prin2 shared a similar structure despite minor differences in the amino acid signatures arising from differences in the scope of the tested species [3, 4, 17, 24]. However, it should be stressed that Prin2 was responsible for domain discrimination rather than growth temperature, because the separation was more related to taxonomic class than to the temperature preferences of the species. For example, the frontiers between Eukaryota and the other domains were clear in this factorial axis (Figure 4). Consequently, we proposed that Prin2 was subjected to the underlying amino acid preferences in the different taxonomic domains. We showed that the distribution of Gln, His, Ser, and Val had distinguishing features between the three domains (Figure 1). For example, the occurrence of Gln and Val was significantly lower in Archaea compared with the other two domains, while the occurrence of Ser was remarkably higher in Eukaryota compared with the other two domains. The occurrence of His also showed a distinguishable distribution in the three domains. Therefore, we suggested that the obvious preferences of these amino acids in Prin2 were in response to the domain discrimination.
Prin3 accounted for 11.06% of the variance in the multivariate analysis and its amino acid signature contained Cys, Leu, and Thr (Additional file 2). Prin3 appeared to be related with environmental niches. In Archaea, the Prin3 values decreased from hyperthermophiles (Prin3 values of 3.3–1.6; representative species Desulfurococcus fermentans, Pyrobaculum aerophilum, Pyrococcus abyssi, Thermococcus barophilus, and Thermoproteus tenax) to acidophiles (Prin3 values of 1.6–1.4; representative species Thermoplasma acidophilum and Sulfolobus acidocaldarius) to mesophiles or anaerobes (Prin3 values of 1.1–-0.2; representative species Methanococcus aeolicus and Methanocorpusculum labreanum) to psychrophiles (Prin3 values of around −0.3; representative species Cenarchaeum symbiosum) to halophiles (Prin3 values of − 0.3– − 3.5; Halobacteria species). Accordingly, in Bacteria, the thermophiles (representative species Anaerolinea thermophila, Meiothermus ruber, Thermodesulfatator indicus and Thermus aquaticus) also gathered at highest end, but with higher Prin3 values (5.9–3.2) than in Archaea. For the Eukaryota species, especially the higher animals and plants, it was more difficult to determine the environmental conditions that could be correlated directly with Prin3 because of the many differences in their morphology and living environments. Interestingly, previous studies have suggested that the amino acid composition of the proteomes in endothermic vertebrates was biased in the same direction as the amino acid composition in thermophilic bacteria and archaea . Together, these results might suggest that Prin3 is subjected to a complex mixture of environmental influences such as growth temperature, pH, oxygen, and solvent, which may provide important new information to refine the global picture of lifestyles and genomes.
Association of amino acid composition dendrogram with genomic GC content, phylogeny, and lifestyles
To compare amino acid composition across the species, organisms with similar compositions for all 20 amino acids were clustered using a neighbor-joining method with distances computed using the Euclidean metric on a dataset that consisted of the overall percent amino acid composition for each of the 461 organisms. The resulting dendrogram had 11 main branches (Figure 5). Comparisons between branches indicated that amino acid composition was associated closely with genomic GC content, phylogeny, and environmental factors. Generally, organisms within a cluster had a similar genomic GC content. However, it should be noted that organisms with similar GC content were not necessarily proximate neighbours. It was worth mentioning that the group F proteomes had the highest GC content (average of 63%) and group E proteomes had the lowest GC content (average of 36%).
Most Eukaryota were in group B, which reflected a distinctive feature of the amino acid architecture of their proteomes, while Archaea and Bacteria revealed general differences in their amino acid composition. The differentiation of the three domains was reflected more obviously in Prin2 (Figure 4). Previous studies have also revealed that proteomic composition could generally discriminate among archaea, bacteria, and eukaryotes [4, 6, 16]. However, interlinkages between proteomic amino acid composition and species phylogeny were limited in characterizing the species diversity. The frequency distributions of the 20 amino acids among the different taxa were compared (Additional file 5). In general, on the one hand, it was common for remote species, such as those in the groups Crenarchaeota and Synergistetes, to share a similar amino acid distribution. On the other hand, species in the phylogenetically related group might also show considerable differences. For example, species in the group Spirochaetes displayed large differences in the amino acids frequencies, as reflected by the CV values. Therefore, information about the amino acid composition of the proteomes was insufficient to reflect the detailed species classification, probably because of the limited information capacity of 20-dimensional amino acid data.
Organisms that share similar lifestyles tended to be grouped together by amino acid composition clustering. For instance, halophiles were clustered in group D while groups C and J consisted mainly of archaea and bacteria hyperthermophiles. Interestingly, a clear distinction was observed between the species in these two groups: species with a genomic GC content around 53% were clustered in group C, and species with a relatively low genomic GC content (around 44%) were clustered in group J. This finding is not surprising because genomic GC content was found to show no correlated response to optimal growth temperatures [40, 41].
Effect of genomic GC content, phylogeny, and lifestyles on frequency concordance of the genetic code and amino acids
Species with similar lifestyles tended to cluster and to have similar concordance. For example, the hyperthermophiles in group C reflected a high adaptability between overall amino acid frequencies and the encoding codon frequencies with correlation coefficients (r) around 0.7, while halophiles in group D showed medium adaptability with correlation coefficients (r) around 0.6 (Additional file 3). This finding suggested that the lifestyles of species may also affect the concordance of the genetic code and amino acids.
In this study, the overall composition, chemical architecture, frequency relationship, main component, and distribution features of the 20 standard amino acids in the proteomes of 461 species from the three domains of life were analyzed systematically. Our results and previous observations together point to the idea that proteomic architecture has been shaped by the integrated forces of GC pressure, phylogeny, and environmental niche during evolution. These findings could contribute to our understanding of proteome evolution, and even to biological evolution.
GC content is the most important determinant of global amino acid composition. In particular, GC content has significant positive effects on the frequencies of Ala, Arg, Gly, Pro, Trp, and Val, and negative effects on the frequencies of Asn, Lys, Ile, Phe, and Tyr. Except for Val, the other GC-dependent amino acids are encoded either by GC-rich or AT-rich codons. However, physico-chemical framework of proteomes is relatively independent against GC pressure through a frequency complementation of GC-dependent amino acid pairs with similar physico-chemical property. On the one hand, our results suggested that a relatively stable physico-chemical framework might be the material basis of the response of species to GC pressure during evolution. On the other hand, the results implied that the use of complementary pairs such as Arg and Lys might be in mutual substitution, which, to a certain degree, is attributable to their similar physico-chemical properties, to buffer the influence of GC pressure on protein structure and function.
Species phylogeny is the second most important factor that influenced proteomic architecture. The cardinal amino acids are Gln, His, Ser and Val, and the constituted component could discriminate Archaea, Bacteria, and Eukaryotes from each other. This finding suggested that the different species phylogenies could be traced based on the proteomic architecture during evolution as early as the divergence of the three domains. Consequently, the proteomic architecture vector might serve as an additional useful tool for understanding taxonomic differences.
Environmental niche brewing a series of habitat factors, such as growth temperature, pH, oxygen and solvent, is also a significant factor to determine the proteomic architecture especially for archaea and its amino acid signature is Cys, Leu and Thr. The environmental component could clearly differentiate lifestyle groups such as hyperthermophiles, acidophiles, mesophiles, and halophiles among the unicellular prokaryotes. However, the relatedness between environmental influences and the proteomes of multicellular eukaryotes was rather weak. It is likely that multicellular higher organisms have a greater capacity to respond to physiological constraints. In other words, external factors have limited impacts on the intracellular environment in multicellular eukaryotes compared with their impact in unicellular prokaryotes, and this ability to respond probably arose from evolutionary resistance to environmental pressure. These results may provide important new information to refine the global picture of the relationships between lifestyles and genomes.
We also found that the frequency concordance of the genetic code and amino acids was not perfect for an individual amino acid. Most of the individual amino acids showed a statistically significant bias against the rates expected from a random distribution of the 61 sense codons. In general, Ala, Asn, Asp, Glu, Lys, Met, and Phe showed significantly higher frequencies than their expected rates, while Arg, Cys, His, Pro, Ser, Thr, and Trp showed significantly lower frequencies. The deviation of some amino acids such as Arg, Cys, Met, and Pro was explained as a consequence of selection pressure at the level of the nucleic acid sequences or protein structure and function. However, we also proposed that the deviation of some amino acids may not have arisen from themselves, but from other amino acids. When certain amino acids are favored or rejected, the frequencies of other amino acids will be affected, even with no bias from them. Overall, the concordance of the amino acids was associated closely with genomic GC content, phylogeny, and lifestyles. Extreme genomic GC contents weakened the concordance and the results suggested that concordance preferentially occurred at a genomic GC content around 56.7%. The concordance in Eukaryota was significantly better than in the other two domains, which we explained as a balance between metabolic efficiency and energy cost in protein biosynthesis. The concordance also differed between species with similar habitats; for example, hyperthermophiles reflected a higher adaptability than halophiles. The observations reported here may promote the understanding of co-adaptation between proteomic architecture and the genetic code.
A total of 461 species covering the Archaea, Bacteria, and Eukaryota domains were used in this study (Additional file 6). The protein sequences from the proteomes of each of the species were extracted from the NCBI Genome database (http://www.ncbi.nlm.nih.gov/genome/browse/).
Analysis of amino acid composition
The percentages of the 20 standard amino acids in the proteomes of each species were calculated by parsing the overall protein sequences. The 20 amino acids were classified into hydrophobic (Ala, Ile, Leu, Met, Phe, Pro, Trp, and Val), charged (Arg, Asp, Glu, His, and Lys), and polar and uncharged amino acids (Asn, Cys, Gln, Gly, Ser, Thr, and Tyr) according to the polarity of their side chains. The 20 amino acids were also divided into aliphatic (Ala, Arg, Asn, Asp, Cys, Gln, Glu, Gly, Ile, Leu, Lys, Met, Ser, Thr, and Val), aromatic (Phe, Trp, and Tyr), and heterocyclic (His and Pro) amino acids based on the chemical structure of their side chains [10, 48]. The proteomes were analyzed by summing the respective contributions of the corresponding amino acids under each of the classifications.
Construction of a clustering tree
The cluster analysis of the proteomes from the 461 species was based on their different amino acid compositions. The Euclidean distance of the amino acid composition was used to reflect the distance between any two species. A neighbor-joining tree was constructed from the distance data using Mega 5.05 . To view the phylogenetic trees and draw the figures, the tree cluster data were submitted to iTOL (http://itol.embl.de/upload.cgi) .
All the data analyses in this study were performed using SAS statistic software (SAS Institute Inc., Cary, North Carolina, USA). The normality of amino acid frequencies was performed with the Shapiro-Wilk test; Pearson correlation coefficients were used to evaluate the correlation between amino acid frequencies; comparative analysis among the amino acid frequencies was performed at the 0.05 level using Duncan's multiple-range test; and principal component analysis was performed to explore the data matrix of amino acid frequencies.
This study was supported financially by the National Natural Science Foundation of China (Program No. 31171649 and 31271834), the Fundamental Research Funds for the Central Universities (Program No. 2012YB06), the National High Technology Research and Development Program of the People's Republic of China (863 Program No. 2006AA10Z1A3), and Program for New Century Excellent Talents in University of the Ministry of Education of the People's Republic of China (NCET-05-0667). The authors thank Amanda Skarlupka (University of Wisconsin-Madison, USA) for reviewing the manuscript.
- Casto AM, Amid C: Beyond the Genome: genomics research ten years after the human genome sequence. Genome Biol. 2010, 11: 309-10.1186/gb-2010-11-11-309.PubMed CentralPubMedView ArticleGoogle Scholar
- Yao T: Bioinformatics for the genomic sciences and towards systems biology. Japanese activities in the post-genome era. Prog Biophys Mol Biol. 2002, 80: 23-42. 10.1016/S0079-6107(02)00011-1.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Tekaia F, Yeramian E: Evolution of proteomes: fundamental signatures and global trends in amino acid compositions. BMC Genomics. 2006, 7: 307-10.1186/1471-2164-7-307.PubMed CentralPubMedView ArticleGoogle Scholar
- Smole Z, Nikolic N, Supek F, Smuc T, Sbalzarini IF, Krisko A: Proteome sequence features carry signatures of the environmental niche of prokaryotes. BMC Evol Biol. 2011, 11: 26-10.1186/1471-2148-11-26.PubMed CentralPubMedView ArticleGoogle Scholar
- Schmidt A, Rzanny M, Hagen M, Schutze E, Kothe E: GC content-independent amino acid patterns in Bacteria and Archaea. J Basic Microbiol. 2012, 52: 195-205. 10.1002/jobm.201100067.PubMedView ArticleGoogle Scholar
- Vieira-Silva S, Rocha EP: An assessment of the impacts of molecular oxygen on the evolution of proteomes. Mol Biol Evol. 2008, 25: 1931-1942. 10.1093/molbev/msn142.PubMed CentralPubMedView ArticleGoogle Scholar
- Sueoka N: Correlation between base composition of deoxyribonucleic acid and amino acid composition of protein. Proc Natl Acad Sci USA. 1961, 47: 1141-1149. 10.1073/pnas.47.8.1141.PubMed CentralPubMedView ArticleGoogle Scholar
- Lobry JR: Influence of genomic G + C content on average amino-acid composition of proteins from 59 bacterial species. Gene. 1997, 205: 309-316. 10.1016/S0378-1119(97)00403-4.PubMedView ArticleGoogle Scholar
- Chen W, Xie T, Shao Y, Chen F: Genomic characteristics comparisons of 12 food-related filamentous fungi in tRNA gene set, codon usage and amino acid composition. Gene. 2012, 497: 116-124. 10.1016/j.gene.2012.01.016.PubMedView ArticleGoogle Scholar
- Bharanidharan D, Bhargavi GR, Uthanumallian K, Gautham N: Correlations between nucleotide frequencies and amino acid composition in 115 bacterial species. Biochem Bioph Res Co. 2004, 315: 1097-1103. 10.1016/j.bbrc.2004.01.129.View ArticleGoogle Scholar
- Lightfield J, Fram NR, Ely B: Across bacterial phyla, distantly-related genomes with similar genomic GC content have similar patterns of amino acid usage. PLoS One. 2011, 6: e17677-10.1371/journal.pone.0017677.PubMed CentralPubMedView ArticleGoogle Scholar
- Singer GAC, Hickey DA: Nucleotide bias causes a genomewide bias in the amino acid composition of proteins. Mol Biol Evol. 2000, 17: 1581-1588. 10.1093/oxfordjournals.molbev.a026257.PubMedView ArticleGoogle Scholar
- Knight RD, Freeland SJ, Landweber LF: A simple model based on mutation and selection explains trends in codon and amino-acid usage and GC composition within and across genomes. Genome Biol. 2001, 2: 0010.Google Scholar
- Nishio Y, Nakamura Y, Kawarabayasi Y, Usuda Y, Kimura E, Sugimoto S, Matsui K, Yamagishi A, Kikuchi H, Ikeo K, Gojobori T: Comparative complete genome sequence analysis of the amino acid replacements responsible for the thermostability of Corynebacterium efficiens. Genome Res. 2003, 13: 1572-1579. 10.1101/gr.1285603.PubMed CentralPubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Dumontier M, Michalickova K, Hogue CWV: Species-specific protein sequence and fold optimizations. BMC Bioinforma. 2002, 3: 39-10.1186/1471-2105-3-39.View ArticleGoogle Scholar
- Nikolic N, Smole Z, Krisko A: Proteomic properties reveal phyloecological clusters of Archaea. PLoS One. 2012, 7: e48231-10.1371/journal.pone.0048231.PubMed CentralPubMedView ArticleGoogle Scholar
- Hickey DA, Singer GAC: Genomic and proteomic adaptations to growth at high temperature. Genome Biol. 2004, 5: 117-10.1186/gb-2004-5-10-117.PubMed CentralPubMedView ArticleGoogle Scholar
- Fukuchi S, Yoshimune K, Wakayama M, Moriguchi M, Nishikawa K: Unique amino acid composition of proteins in halophilic bacteria. J Mol Biol. 2003, 327: 347-357. 10.1016/S0022-2836(03)00150-5.PubMedView ArticleGoogle Scholar
- Bardavid RE, Oren A: The amino acid composition of proteins from anaerobic halophilic bacteria of the order Halanaerobiales. Extremophiles. 2012, 16: 567-572. 10.1007/s00792-012-0455-y.View ArticleGoogle Scholar
- Brocchieri L: Environmental signatures in proteome properties. Proc Natl Acad Sci U S A. 2004, 101: 8257-8258. 10.1073/pnas.0402797101.PubMed CentralPubMedView ArticleGoogle Scholar
- Singer GAC, Hickey DA: Thermophilic prokaryotes have characteristic patterns of codon usage, amino acid composition and nucleotide content. Gene. 2003, 317: 39-47.PubMedView ArticleGoogle Scholar
- Puigbo P, Pasamontes A, Garcia-Vallve S: Gaining and losing the thermophilic adaptation in prokaryotes. Trends Genet. 2008, 24: 10-14. 10.1016/j.tig.2007.10.005.PubMedView ArticleGoogle Scholar
- Zeldovich KB, Berezovsky IN, Shakhnovich EI: Protein and DNA sequence determinants of thermophilic adaptation. PLoS Comput Biol. 2007, 3: e5-10.1371/journal.pcbi.0030005.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang GZ, Lercher MJ: Amino acid composition in endothermic vertebrates is biased in the same direction as in thermophilic prokaryotes. BMC Evol Biol. 2010, 10: 263-10.1186/1471-2148-10-263.PubMed CentralPubMedView ArticleGoogle Scholar
- Jobson RW, Qiu YL: Amino acid compositional shifts during streptophyte transitions to terrestrial habitats. J Mol Evol. 2011, 72: 204-214. 10.1007/s00239-010-9416-1.PubMedView ArticleGoogle Scholar
- Hanai R, Wada A: The effects of guanine and cytosine variation on dinucleotide frequency and amino acid composition in the human genome. J Mol Evol. 1988, 27: 321-325. 10.1007/BF02101194.PubMedView ArticleGoogle Scholar
- Karlin S, Burge C: Dinucleotide relative abundance extremes - a genomic signature. Trends Genet. 1995, 11: 283-290. 10.1016/S0168-9525(00)89076-9.PubMedView ArticleGoogle Scholar
- De Amicis F, Marchetti S: Intercodon dinucleotides affect codon choice in plant genes. Nucleic Acids Res. 2000, 28: 3339-3345. 10.1093/nar/28.17.3339.PubMed CentralPubMedView ArticleGoogle Scholar
- Gentles AJ, Karlin S: Genome-scale compositional comparisons in eukaryotes. Genome Res. 2001, 11: 540-546. 10.1101/gr.163101.PubMed CentralPubMedView ArticleGoogle Scholar
- Alves R, Savageau MA: Evidence of selection for low cognate amino acid bias in amino acid biosynthetic enzymes. Mol Microbiol. 2005, 56: 1017-1034. 10.1111/j.1365-2958.2005.04566.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Bragg JG, Thomas D, Baudouin-Cornu P: Variation among species in proteomic sulphur content is related to environmental conditions. Proc R Soc B. 2006, 273: 1293-1300. 10.1098/rspb.2005.3441.PubMed CentralPubMedView ArticleGoogle Scholar
- Akashi H, Gojobori T: Metabolic efficiency and amino acid composition in the proteomes of Escherichia coli and Bacillus subtilis. Proc Natl Acad Sci USA. 2002, 99: 3695-3700. 10.1073/pnas.062526999.PubMed CentralPubMedView ArticleGoogle Scholar
- Jun T, Robert N, Erika W, Erin M, Brad M, Grant H, Tygue A, Kelly H, Rachel H, Manuel R, Mike W: The frequencies of amino acids encoded by genomes that utilize standard and nonstandard genetic codes. BIOS. 2010, 81: 22-31. 10.1893/011.081.0103.View ArticleGoogle Scholar
- Clancy S, Brown W: Translation: DNA to mRNA to protein. Nat Educ. 2008, 1: 1.Google Scholar
- Wu HL, Bagby S, van den Elsen JMH: Evolution of the genetic triplet code via two types of doublet codons. J Mol Evol. 2005, 61: 54-64. 10.1007/s00239-004-0224-3.PubMedView ArticleGoogle Scholar
- Tlusty T: A colorful origin for the genetic code: Information theory, statistical mechanics and the emergence of molecular codes. Phys Life Rev. 2010, 7: 362-376. 10.1016/j.plrev.2010.06.002.PubMedView ArticleGoogle Scholar
- Gilis D, Massar S, Cerf NJ, Rooman M: Optimality of the genetic code with respect to protein stability and amino-acid frequencies. Genome Biol. 2001, 2: 0049.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Hurst LD, Merchant AR: High guanine–cytosine content is not an adaptation to high temperature: a comparative analysis amongst prokaryotes. Proc R Soc B. 2001, 268: 493-497. 10.1098/rspb.2000.1397.PubMed CentralPubMedView ArticleGoogle Scholar
- Smith DR, Chapman MR: Economical evolution: Microbes reduce the synthetic cost of extracellular proteins. mBio. 2010, 1: e00131-10.PubMed CentralPubMedGoogle Scholar
- Heizer EM, Raiford DW, Raymer ML, Doom TE, Miller RV, Krane DE: Amino acid cost and codon-usage biases in 6 prokaryotic genomes: A whole-genome analysis. Mol Biol Evol. 2006, 23: 1670-1680. 10.1093/molbev/msl029.PubMedView ArticleGoogle Scholar
- Raiford DW, Heizer EM, Miller RV, Doom TE, Raymer ML, Krane DE: Metabolic and translational efficiency in microbial organisms. J Mol Evol. 2012, 74: 206-216. 10.1007/s00239-012-9500-9.PubMedView ArticleGoogle Scholar
- Raiford DW, Heizer EM, Miller RV, Akashi H, Raymer ML, Krane DE: Do amino acid biosynthetic costs constrain protein evolution in Saccharomyces cerevisiae?. J Mol Evol. 2008, 67: 621-630. 10.1007/s00239-008-9162-9.PubMedView ArticleGoogle Scholar
- Barton MD, Delneri D, Oliver SG, Rattray M, Bergman CM: Evolutionary systems biology of amino acid biosynthetic cost in yeast. PLoS One. 2010, 5: e11935-10.1371/journal.pone.0011935.PubMed CentralPubMedView ArticleGoogle Scholar
- Seligmann H: Cost-minimization of amino acid usage. J Mol Evol. 2003, 56: 151-161. 10.1007/s00239-002-2388-z.PubMedView ArticleGoogle Scholar
- Wang J, Zhu S, Xu C: Edited by: Xu C. 2002, Beijing: Higher Education Press, 124-127. Amino acids, 3, Biochemistry, Volume 1,
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.PubMed CentralPubMedView ArticleGoogle Scholar
- Letunic I, Bork P: Interactive Tree Of Life v2: online annotation and display of phylogenetic trees made easy. Nucleic Acids Res. 2011, 39: W475-W478. 10.1093/nar/gkr201.PubMed CentralPubMedView ArticleGoogle Scholar