- Research article
- Open access
- Published:
Hypothesis: protein and RNA attributes are continuously optimized over time
BMC Genomics volume 20, Article number: 1012 (2019)
Abstract
Background
Little is known why proteins and RNAs exhibit half-lives varying over several magnitudes. Despite many efforts, a conclusive link between half-lives and gene function could not be established suggesting that other determinants may influence these molecular attributes.
Results
Here, I find that with increasing gene age there is a gradual and significant increase of protein and RNA half-lives, protein structure, and other molecular attributes that tend to affect protein abundance. These observations are accommodated in a hypothesis which posits that new genes at ‘birth’ are not optimized and thus their products exhibit low half-lives and less structure but continuous mutagenesis eventually improves these attributes. Thus, the protein and RNA products of the oldest genes obtained their high degrees of stability and structure only after billions of years while the products of younger genes had less time to be optimized and are therefore less stable and structured. Because more stable proteins with lower turnover require less transcription to maintain the same level of abundance, reduced transcription-associated mutagenesis (TAM) would fixate the changes by increasing gene conservation.
Conclusions
Consequently, the currently observed diversity of molecular attributes is a snapshot of gene products being at different stages along their temporal path of optimization.
Background
Typical protein attributes such as abundance, stability, number of protein-protein interactions, or folded vs. unfolded regions are increasingly characterized on a global level. Large-scale analyses offer the opportunity for an unbiased characterization of such innate molecular attributes. For example, protein half-life and turnover is an attribute which has been thoroughly investigated [1,2,3,4,5,6,7,8,9]. However, analyses of protein turnover in a variety of organisms including bacteria [3], yeast [2, 7], mammalian cells [1, 4, 9], whole animals [5, 8], and even humans [6] have not yielded any significant insight as to why turnover values among proteins of the same cell can vary by orders of magnitudes. Similarly, also RNA half-lives in human B-cells ranged from minutes to days [10]. Some biological correlations have been identified, for example membrane proteins exhibit a significantly higher protein turnover than cytosolic proteins and phosphorylated proteins have a higher turnover than non-phosporylated proteins [11]. Yet, a compelling protein function vs. protein turnover correlation could not be established though.
During evolution, mutations to genes lead to functional divergence of their protein products but they also affect innate molecular attributes such as turnover or stability. However, if and how continuous mutagenesis changed innate molecular attributes over time is rarely discussed [12, 13]. There have been few reports of correlations between protein age and percentage of protein disorder (negative correlation) [13] or protein age and protein packing density [12]. Here, I present over a dozen additional correlations between gene age and various molecular attributes. For example, ‘old’ proteins that already existed in unicellular organisms on average have a lower turnover than ‘young’ proteins which appeared more recently. Together with other already published correlations, some also reaffirmed here, these observations suggest that there is a continuous and gradual change of different molecular attributes over time through nonsynonymous mutations. Obviously, there are countless and diverse molecular attributes such as the propensity for protein-protein-interaction [14] or the length of poly(A) tails [15]. Consequently, there is vast literature on molecular attributes and thus it is important to note that the claim here is neither to have uncovered all possible attribute-gene age correlations nor to be the first to present them specifically. Rather, the scope of this research was to show an overarching trend of attribute optimization over time by analyzing many different molecular attributes by the same, simple correlation with gene age. I find that older genes tend to produce more stable and structured proteins and mRNAs. I present a hypothesis that suggests that such optimized molecular attributes arise from cumulative mutational drifts of old genes. Consequently, young genes produce less optimized molecules. Just how mutations that favor attribute optimization tend to accumulate and become fixated is also discussed. Rather than invoking increased cellular fitness as a driving force for selection, I postulate that optimized, stable proteins and mRNA reduce the need for transcription. In turn, less transcription reduces transcription-associated mutagenesis (TAM) at this specific gene locus so that these optimizing mutations become fixated.
Results
Molecular stability and gene conservation correlate with gene age
Triplicate, high-throughput mass spectrometry (MS) was previously used to analyze and compare protein turnover in non-dividing arrested human cervical HeLa and differentiated mouse muscle C2C12 cells [11]. MS based on SILAC (stable isotope labeling of amino acids in cell culture) labeled amino acids can be used for quantitative protein abundance comparison between samples [16]. Similar to incorporation experiments with radioisotope-labeled amino acids decades ago [17, 18], SILAC allows the analysis of turnover of thousands of proteins. Sub-saturating metabolic incorporation of ‘heavy’ isotopes produced a turnover value (heavy / unlabeled ratio after 24 h) for 4106 human and 3574 mouse proteins [11]. Here, the half-lives of proteins and their corresponding gene age were compared to reveal if there is a general correlation between the time of existence of a gene and the stability of its protein product. Gene ages were obtained from the ProteinHistorian Database [19] and genes were taxonomically grouped as unicellular organisms (u org), unicellular eukaryotes (u euk), Ophistokonta/Bilateria/Deuterostomia (OBD), chordates (chor), or mammals (mamm). This grouping was chosen to reflect major steps in evolution. So all genes for which the protein turnover and the gene age was available were assigned to one of these five groups according to the gene age and the median protein turnover for all genes in the group was determined. It is important to note that in the figures, u org, u euk, OBD, chor, and mamm, indicate the age of the genes, not the origin. Thus, ‘chor’ genes originally appeared during the time when chordates first existed. The grouping according to gene age thus allowed comparison of human genes that already existed in prokaryotes billions of years ago with those human genes that appeared more recently and existed only since the age of mammals.
For human proteins, the correlation of decreasing protein turnover with increasing protein age was significant (Spearman’s correlation, r = − 0.20, P < 0.0001) as were most differences between taxonomically grouped proteins (Fig. 1a) (ANOVA, Bonferroni post-hoc analysis; also Fig. 1c-h). A box plot of the same data is shown in Additional file 1: Figure S1a. Notably, the variance of turnover values among all proteins in each of the five groups substantially decreased with increasing protein age (Fig. 1b) indicating that turnover was more uniform for old proteins. Similarly, the protein turnover values derived from the mouse C2C12 MS experiments correlated with gene age as well (Spearman’s correlation, r = − 0.31, P < 0.0001) (Additional file 1: Figure S1b). Protein turnover of rat proteins showed the same trend (Additional file 1: Figure S1c). Analysis in Schizosaccharomyces pombe indicated that older genes produce proteins with longer half-lives (Additional file 1: Figure S1d). Moreover, yeast Saccharomyces cerevisiae proteins with prokaryotic orthologues [20] exhibited longer half-lives [2] compared to those without (52 vs. 40 min, P < 0.0001, Mann-Whitney test). These correlations are supported by a previous publication showing that human protein stability in terms of free energy folding ΔG was higher for old genes vs. young ones [21]. In summary, these data suggested that proteins from older genes on average exhibit a lower turnover than proteins from younger genes. To test if additional molecular attributes also follow such a trend, various other attributes were further examined.
Protein aggregation, once thought to be a characteristic of diverse diseases such as Alzheimer’s or Parkinson’s disease, is now considered to be more of a generic property of polypeptide chains [22]. I found that the overall strength of aggregation nucleating regions per protein significantly decreased with gene age, albeit weakly (Spearman’s correlation, r = − 0.11, P < 0.0001) (Fig. 1c). Additionally, it was also reported that proteins with high turnover were found to have an increased propensity to aggregate [23]. Together, this suggests that older proteins have a lower tendency to aggregate than younger ones.
Since it was demonstrated that the predicted extent of intrinsically unstructured protein (IUP) levels negatively correlated with protein half-lives [24], protein structure may itself be influenced by protein age. Indeed, the levels of unstructured regions in human proteins significantly decreased with increasing protein age (Spearman’s correlation, r = − 0.18, P < 0.0001) although there is a minor decrease rather than increase from chordates to mammals (Fig. 1d). When analyzing yeast Saccharomyces cerevisiae genes the same way, there was an even more pronounced correlation between protein structure and gene age (Spearman’s correlation, r = − 0.35, P < 0.0001) (Additional file 2: Figure S2a). In addition, when comparing different species, prokaryotes have been found to exhibit significantly less disordered proteins compared to eukaryotes [25, 26]. Thus, longer existing proteins are on average more structured as was independently demonstrated in a recent study [27].
Protein abundance in a data set can be approximated by summed MS peptides intensities [28, 29] as the measured peptide signal is greater for abundant proteins. In the HeLa data set, protein abundance was another molecular attribute that changed over time, as older human proteins were significantly more abundant (Spearman’s correlation r = 0.12, P < 0.0001) (Fig. 1e). Similar trends have been observed before [30] albeit often without statistics to support it. In the present study, the MS HeLa protein abundance analyses were not skewed towards high abundant proteins since abundance varied over five orders of magnitude and exhibited a bell-shaped distribution (Additional file 2: Figure S2b). Moreover, almost one hundred human transcription factors, i.e. proteins considered to be of low abundance, were detected in the HeLa lysates [11]. Similar to protein abundance, mRNA abundance was also significantly higher for old genes compared to young ones (Spearman’s correlation r = 0.34, P < 0.0001) (Additional file 2: Figure S2c).
RNA half-lives were larger for old genes compared young genes (Spearman’s correlation for human, r = 0.23, P < 0.0001; Spearman’s correlation for mouse, r = 0.26, P < 0.0001) (human in Fig. 1f, mouse in Additional file 2: Figure S2d). Analysis of RNA secondary structure in Saccharomyces cerevisiae genes revealed a higher level of structure, i.e. a higher average Parallel Analysis of RNA Structure (PARS) score [31], for older genes with prokaryotic orthologs (0.28 vs. 0.23, P < 0.0001, Mann-Whitney test) compared to younger genes without orthologs.
The translation initiation efficiency describes how well a particular mRNA assembles with 80S ribosomes which is an indication of the fraction of mRNA molecules that can be translated. Genes with higher values are thus able to produce more proteins within a certain time window. Using published data on the translation initiation efficiency in human HeLa cells [32] older genes on average were found to have higher translation efficiencies than younger ones (Fig. 1g) (Spearman’s correlation r = 0.08, P < 0.0001).
Taken together, the data demonstrate in different organisms that newly evolved genes tend to exhibit molecular attributes that are not well optimized, such as less stability and less structure of RNAs and proteins. In addition, the human gene conservation score [33] was as expected significantly lower in young compared to old genes (Spearman’s correlation, r = 0.67, P < 0.0001) (Fig. 1h). When analyzing more defined gene age classes (Additional file 2: Figure S2e), this correlation with conservation becomes even more obvious with an almost steady decrease with decreasing gene age, an observation that was not reported in previous publications [34].
Continued mutagenesis decreases individual protein turnover variance
Figure 1b demonstrated that the spread, i.e. variance of human turnover values of all proteins in a respective taxonomic group was highest for young proteins. The question then arose if the spread of measured turnover values for each specific protein also changed with age. In SILAC-MS experiments [35], the various peptides that identify a specific protein are often quantified multiple times and so one can determine the variance of the measured turnover values for each protein separately. Indeed, the protein-specific variance was highest for high turnover proteins and thus the variance indirectly related to gene age. Plotting the variance of measured peptide values versus protein turnover shows (Fig. 2a,b) that the peptide value variance was much greater for high turnover compared to low turnover proteins. These data therefore suggest that for a few low turnover proteins, the variance can decrease to a level of near uniformity. For display purposes, all proteins whose overall turnover value happened to be based on 15 peptide measurements are shown in Fig. 2c and the spread of variances recapitulates that pattern seen in Fig. 2a. The variance of three exemplary proteins (Fig. 2d) and normalization of the respective peptide values confirmed the trend (Fig. 2e) of low turnover proteins exhibiting much less variance. Therefore, each individual molecule of an old protein exists more or less the same length of time before it gets degraded. Conversely, for a protein of a younger gene, the time it takes to be degraded following its synthesis can vary substantially.
Molecular attributes change independently of function and essentiality
It has long been assumed that core cellular proteins (‘housekeeping genes’) are ubiquitously expressed and intuitively should be more conserved [36]. However, from young to old human housekeeping genes [37], not only did I observe a continuous increase in gene conservation, but also in protein structure as well as mRNA and protein turnover (Additional file 3: Figure S3a,b,c,d). The family of human transcription factors [38] exhibited the same trend (Additional file 4: Figure S4a,b,c,d). Both sets of data thus suggest that the optimization of molecular attributes over time occurs independently of essentiality and function. Also, the mean conservation score of human genes (0.71) considered to be essential [39] was significantly lower than the mean of those genes that first appeared in prokaryotes (0.82; P < 0.0001, Mann-Whitney test) or unicellular eukaryotes (0.78; P = 0.0002, Mann-Whitney test). Thus, old genes are more conserved than essential genes and consequently protein function appears to have lesser role in determining gene conservation than generally assumed. In addition, essential mouse genes [40] also showed an increase towards higher mRNA and protein stability from young genes to old ones (Additional file 4: Figure S4e,f).
GO categories
GO (gene ontology) biological process analyses of human and mouse turnover data allowed ranking of GO categories according to their median turnover values (Additional file 6: Tables S1,S2). In the present HeLa and C2C12 data sets, core cellular and metabolic categories that must have appeared early in evolution such as ‘protein synthesis’ displayed low median turnover values while categories central to the appearance of multi-cellular organisms and vertebrates exhibited high and even higher values, respectively. Analogous trends were previously also seen in other analyses such as GO categorization of protein conservation scores [41] or RNA half-lives [10]. Comparing the order of almost 600 GO categories between mouse and human indicated a high similarity (Spearman’s correlation, r = 0.67, P < 0.0001). However, because of extensive cross-annotations and functional overlap between categories, it was not possible to establish a common temporal ranking based on both species. Instead, a manually curated selection of human categories is shown in Table 1. The data (Additional file 6: Tables S1,S2) suggest that one of the first biological processes to exist was ATP production. Other basic biological processes such as translation or biosynthesis of carbohydrates also exhibited low average turnover values while biological processes related to complex multicellular organisms such as embryonic morphogenesis or axon guidance had significantly higher values. Undoubtedly, this ranking provides an appealing ‘time stamp’ for the time of emergence of biological processes although this ranking must be considered preliminary. Such an approach can clearly not be extended to individual proteins as the turnover based GO ranking only became meaningful if several proteins were included. Nevertheless, one can approximate the coarse order of many of these processes based on known cellular physiology. For example, nuclear transport as GO category must have appeared in eukaryotes not prokaryotes, while the GO categories translation and carbohydrate metabolism must have appeared in prokaryotes. Similarly, GO categories such as cell motion or response to chemical stimulus are relevant to unicellular organisms and thus appeared early in evolution while GO categories such as organ or nervous system development are key features of multi-cellular organisms and therefore must have appeared later. In conclusion, the GO data independently support the gene age – protein turnover correlation.
A hypothesis for dynamic molecular attribute optimization
The data discussed above showed that in different organisms, attributes such as half-life of proteins and RNAs significantly change over time. This temporal correlation can be interpreted in two ways. First, the molecular attributes of molecules that appeared billions of years ago are very different from those that arose only millions of years ago. This implies creation of stable, structured proteins with long mRNA half-lives in prokaryotes, but less stable and less structured proteins with shorter mRNA half-lives in mammals, and proteins with intermediate qualities in between. This broadly relates to the ‘constant restraint’ model [34, 42]. A second, alternative explanation would be that newly evolved genes, from the time of prokaryotic life to the age of mammals, always have the tendency to produce more unstructured and unstable molecules when they first appear. Continued mutagenesis then gradually optimizes the respective molecular attributes over time so that the oldest, prokaryotic genes were subjected to the most extensive optimization and thus produce the most structured proteins as well as the most stable mRNAs and proteins today. Genes that appeared later when multicellular organisms and invertebrates first inhabited Earth had less time available and therefore were less optimized compared to prokaryotic genes, but are more optimized compared to those that arose in mammals. Genes that appeared the latest at the age of mammals had the least time available for optimization and thus produce the least structured and least stable proteins and mRNAs. This hypothesis is graphically visualized in Fig. 3. Genes A, B, and C represent genes that first appeared in prokaryotes, multicellular organisms/invertebrates, or mammals, respectively, and so gene A is older than B, and B is older than C. The vertical axis indicates the degree of molecular stability, structure, and gene conservation - unstructured, unstable molecules from less conserved genes are at the top while structured, stable molecules from conserved genes are at the bottom of the axis. The stippled lines reflect the path to greater optimization of each gene over time. ‘Newborn’, imperfect genes generally start their path near the top of the axis and then gradually, but not steadily, their molecular attributes change to reach the bottom of the graph. The hypothesis thus proposes that proteins and mRNAs which were present in the first forms of unicellular life were optimized subsequently during the billions of years until today to now have a low uniform turnover and high conservation. Molecules that appeared later in evolution also change but are more ‘work in progress’ and will, like all other proteins and mRNAs, continue to be optimized. Consequently, the wide spectra of half-lives, protein disorder, and degrees of conservation we currently observe actually represent a snapshot of molecules being at different stages along their evolutionary paths.
Inter-species differences
So far, all the presented data were based on analyses of proteins or RNAs within one species with the observation that longer existing genes and their products were different to younger ones. Because of their longer existence, older genes have been subject to more mutagenesis over time and the extent of mutagenesis a gene experienced determines the extent of optimization of molecular attributes. Consequently, within one species, the products of older genes are more stable and more structured. Overall, when attribute data from two separate species could be obtained, consistently I found a larger age-attribute correlation value for those species that have a higher cumulative mutational load. For example, when comparing yeast with human proteins, yeast proteins displayed a higher correlation between gene age and unstructured protein regions. Also, there was a higher correlation for gene age and RNA turnover in mice than men. Thus, the prediction would be that proteins from species with higher mutational rates are more optimized. Comparing protein turnover in mammalian cells vs. lower organisms would be difficult as metabolic rates, cell cycles, or body temperatures can be very different. Therefore, protein turnover was compared in human vs. mouse using the two aforementioned data sets from non-dividing HeLa and C2C12 cells [11]. Mice have higher mutational rates [43] and therefore the expectation would be that the average protein turnover is decreased compared to human. The mean turnover value in human was 2.2 and in mice 2.0 (P < 0.0001, Mann-Whitney test) (Fig. 4a). Most (73%) of the 2107 protein homologs identified in both species exhibited a larger value in humans. The same trend was seen in individual subunits of several previously characterized protein complexes such as the COP9 complex (Fig. 4b, Additional file 5: Table S3, Figure S5a,b). Also, the fraction of human proteins with lower turnover values in mice was increasing with decreasing turnover values (Additional file 5: Figure S5c). Thus, the lower the human turnover value, the higher the probability that the corresponding mouse homolog displayed an even smaller value. The data therefore support the idea that increased mutagenesis leads to more optimized proteins as mouse orthologs exhibited on average lower turnover compared to human. If this were the case, one would expect that less optimized human proteins required more ‘cellular support’ to maintain proper functioning. Relative quantitative proteome comparison of mouse and human data showed that heat shock proteins were much more abundant in human compared to mouse cells (Fig. 4c). This is in line with the heat shock capacitor hypothesis that views molecular chaperones as key mediators of adaptive evolution by buffering genetic variation [44]. Of course, although the data showed highly significant trends that were as predicted, further analyses are certainly necessary to confirm the results of such inter-species comparisons. Nevertheless, data from both, intra-species and inter-species comparisons suggest that the cumulative levels of mutagenesis affect molecular attributes.
Protein turnover and cellular fitness
To assess if high turnover proteins affect cellular fitness, the protein turnover values of all genes in the ENSEMBL database known to produce lethal phenotypes were correlated. Such an unbiased approach did not reveal any obvious association between known lethal diseases and human proteins of high turnover. However, there was a clear association between lethal phenotypes and the ratio of human/mouse turnover values from orthologs detected in both species. As expected from the high correlation of turnover values between both species [11], most orthologs displayed a human/mouse turnover ratio close to 1 (blue dots, Fig. 4d). Bins of proteins (30 each) with a ratio close to 1 (bins 31–60, 61–90, and 91–120) had a total of zero, one, two, or three lethal diseases associated (brown bars and dots). But for orthologs with very high or low ratios however, there were up to eight counts of diseases. These extreme ratios and thus extreme differences between human and mouse turnover values of the same protein may be an indication that these proteins are particularly prone to mutagenesis and thus one could speculate that they have a higher chance to acquire lethal phenotypes.
Continued mutagenesis and post-translational modifications
We previously showed that phosphorylated proteins have a significantly higher turnover than non-phosphorylated proteins [11]. I now find that acetylated [45] and methylated [46] proteins have a significantly lower turnover than phosphorylated proteins (Fig. 4e). Since the data presented here suggest that continued mutagenesis decreases protein turnover, this mutagenesis may also influence the type of post-translational modifications (PTMs). The abundance of post-translational modifications in humans is phosphorylation > > acetylation > methylation [45, 47, 48]. This is opposite to prokaryotes such as Escherichia coli where the abundance of phosphorylation < acetylation [49, 50] or the spirochete Leptospira interrogans where phosphorylation < acetylation < methylation [51]. Thus, prokaryotes tend to have proteins that are methylated and acetylated while phosphorylation is predominant in humans. Similar to other molecular attributes, PTMs could potentially correlate with gene age as well. Analysis of PTMs in human proteins showed that acetylation was significantly more associated with older proteins than phosphorylation (Mann-Whitney test, P < 0.0001). For example, for very old proteins that originated during the stage of ‘unicellular organisms’, the percentage of acetylated proteins was twice as high compared to the percentage of phosphorylated proteins (Fig. 4f). Given that young proteins tend to be phosphorylated and older proteins are more likely to be acetylated, I would like to speculate that post-translational modifications are also molecular attributes that are undergoing ‘optimization’ because of continued mutagenesis. Notably, phosphorylation was found to occur mostly in unstructured regions of proteins [52] furthering the notion that younger genes produce proteins with less structure and more phosphorylation. The analyses presented here thus suggest that the existence of PTMs on proteins is ultimately transient as proteins become older. As the hydrophobicity of the actual PTM increases from phosphorylation to acetylation and then methylation, water solubility and protein-protein interactions should be affected. According to the hypothesis, proteins are therefore more likely to be phosphorylated at ‘birth’, yet they eventually lose this modification and may acquire an acetylation or methylation subsequently as mutations optimize the proteins towards higher stability. However, it is not clear at this point whether changes in PTMs relate at all to changes protein turnover. Turnover and PTMs may simply be independent manifestations of continued mutagenesis. Also, whether the changing PTMs locate to the same or different positions within the protein will have to be determined.
Discussion
Here, I presented data which show that molecules have distinct properties according to their time of birth independent of function, essentiality, or expression. Unless one were to postulate that nature ‘by design’ progressively produces increasingly unstructured and unstable proteins and RNAs, the data are best accommodated by the proposed hypothesis assuming that new genes are imperfect at the beginning and are then subjected to optimization over time. At least in respect to protein turnover, this optimization also reduces variance. The surprisingly smooth, almost steady decrease of conservation for genes that appeared later in evolution (Additional file 2: Figure S2e) is an appealing representation of such optimization. Importantly, the hypothesis does not require that the different attributes of a specific molecule must correlate. In other words, a gene that produces a stable mRNA may in turn give rise to a protein with high turnover because a mutation that changes mRNA stability may have little effect on protein turnover or vice versa. Indeed, published data showed that protein and mRNA half-lives did not correlate [53]. Moreover, the hypothesis allows for genetic selection to maintain high turnover or less structure where functionally required. For example, proteins that functionally require short half-lives such as cyclines possibly optimize only mRNA half-lives and protein structure, independent of a selective pressure that supports high protein turnover.
Protein abundance also correlated with gene age and molecular optimization may indeed increase abundance over time, similar to structure und half-lives (Fig. 1a,d). However, a different interpretation of the data would be that the increased protein abundance is a direct consequence of increased protein structure, increased protein half-lives, and increased mRNA half-lives. Especially changes in mRNA and protein stability should lead to accumulation of protein and both are expected to act independently of each other [53]. Not surprisingly, the abundance of mRNAs was also higher for old genes compared to young ones.
The observed gradual change of attributes could be the consequence of two opposing scenarios. These scenarios are highlighted below using the degree of structure of transcription factors for illustration. In one scenario, nature produced highly structured transcription factors billions of years ago but created transcription factors with much more disorder millions of years ago, and transcription factors with intermediate structure in between. Thus, the molecular attributes of transcription factors were very different depending on the time when they were ‘born’. This appears to be prevailing view of molecular attribute genesis, but an understanding of the molecular basis and the driving force behind this scenario is vague at best [27]. It is difficult to imagine that protein function is the underlying driving force responsible for producing less and less structured transcription factors over time when the principles of transcription are the same for all transcription factors, young and old. If having less structure were beneficial, why did nature only produce such transcription factors millions, but not billions of years ago?
Conversely, a second, more compelling scenario to explain the correlation of molecular attributes with gene age is the following: Rather than assuming that newly existing molecules of recent times have different attributes than those that appeared long time ago, I propose a hypothesis that posits that new genes always, from billions of years ago to millions years ago, had imperfect attributes at the beginning of their existence which were then optimized through mutagenesis over time. This implies that a transcription factor that was present billions of years ago also exhibited less protein structure billions of years ago, but nature optimized its various attributes so that the same transcription factor today is highly structured. I also described the correlation of several other molecular attributes with gene age, including half-lives of proteins and RNAs. Together, all these observations suggest that mutations not only change the function and activity of a molecule but also dynamically and continuously change its innate attributes.
How are the different optimized attributes fixated and selected for during evolution? For example, reduced turnover requires less energy for cellular homeostasis, but these changes are incremental as the energy needs for the production of one protein compared to all proteins in a cell is negligible. With more than 10,000 different protein species detected in eukaryotic cells [54], changing the turnover of a single specific protein will correspondingly affect about 1/10000 of the energy used for protein production. The degree to which this confers a competitive growth advantage is not clear. The cellular consequences of stepwise optimizing any of the other attributes also appear negligible. I therefore wish to propose an alternative explanation. Hypothetically, if an mRNA species becomes ten times more stable, ten times less transcription should be necessary to maintain equal levels of steady-state mRNA abundance. I showed here that old genes produce more stable mRNAs and proteins, and so their individual transcription rates should be reduced. Use-dependent transcription-associated mutagenesis (TAM) of genomic DNA [55] is a possible link between optimization, transcription rates, and fixation of mutations. TAM is locally altering and destabilizing the genomic DNA template through a variety of different mechanisms [56]. Thus, instead of invoking selective pressure for such mutations that optimize molecular attributes such as protein turnover, perhaps these mutations become fixated because they directly reduce mutagenesis at their own gene locus. If a random mutation produces more stable mRNA or protein, less transcription should be necessary to reach normal levels of abundance and consequently, fixation of this gene increases as it is more protected from mutations because of reduced TAM. Indeed, I found a significant negative correlation between gene conservation scores and corresponding transcriptional rates [10] (Spearman’s correlation, r = − 0.18, P < 0.0001) which suggests that reduced transcription accounts to some extent for the increased fixation of old genes. From a molecular perspective, any random mutation that reduces the necessity for transcription, i.e. by ultimately increasing mRNA and protein abundance, should principally be fixated because of TAM. Such fixation would occur ‘passively’ and not by selection based on function. To this end, it is known that old, conserved housekeeping genes generally exhibit lower evolutionary rates [57]. Since I showed here that old genes have on average more stable mRNAs and proteins, the reduced TAM should be in line with their lower evolutionary rates.
Conclusions
Molecular attributes such as protein or RNA stability often vary over several magnitudes but a conclusive explanation for this variability has been missing. A new hypothesis was introduced that views individual attributes of every molecule as being at different stages along a path of optimization based on continued mutagenesis over time. The hypothesis does not attempt to link the various attributes of the same molecule with function but rather allows attributes to be affected independently by a mutation. Thus, the biological challenge would be to maintain molecular function in the face of ever changing molecular attributes. Based on these assumptions, this hypothesis can be tested by in vitro evolution.
Methods
All analyses were based on previously published data. Statistical data analyses were performed using Prism 5.0 (GraphPad Software, San Diego, CA USA). Each gene obtained from the different molecular attributes data sets was given its published age as described in the database for human gene/protein age (‘ProteinHistorian’/ http://lighthouse.ucsf.edu/ProteinHistorian/) [19]. According to this age, genes in each data set were grouped. Each taxonomic grouping of the ProteinHistorian output is indicated in Additional file 2: Figure S2e, together with the corresponding gene ages.
Protein and peptide turnover values as well as protein abundance were from a previous publication [11]. RNA folding energies (PARS) [31] (http://genie.weizmann.ac.il/pubs/PARS10/pars10_catalogs.html), gene conservation scores (https://dsgweb.wustl.edu/hutz/index.html), and bacterial orthologues in yeast (http://makana.berkeley.edu/phylofacts/) [20] were obtained from publicly available databases. RNA half-lives were previously published [10] from which also the relative transcription rates were derived assuming that transcription rates are proportional to degradation rates at steady state. Yeast protein half-lives [2], human transcription factors [38], human housekeeping genes [57], human [39] and mouse [40] essential genes, acetylated [45] and phosphorylated [47] proteins, human protein aggregation propensity [58], human translation initiation values [32], human mRNA abundance [59] were each obtained from the indicated references. The mean normalized peptide variance was defined as (mean variance) / (mean turnover value)2 for all proteins in the respective bin (93 proteins/bin). Peptides were required to be quantified with at least three peptides during the MS experiments.
Ranking of biological process GO categories was based on the median protein turnover of each particular category. In parallel, 1000 sets of random protein turnover values were iteratively created (“bootstrapping”) where each set contained as many values as the number in the respective category and the mean of all 1000 random sets produced the bootstrap distribution. The difference between the bootstrap distribution and the mean of the category was computed in standard deviations and yielded a probability value for a null correlation [11]. An unbiased description of diseases associated with human proteins was obtained from the biomart database ‘MIM Morbid’ (www.ensemble.org). The MaxQuant software [60] was used for label-free quantitative analysis [61] of mouse and human proteomes [11].
Availability of data and materials
For every figure, all source data have been published before and can be obtained from the indicated publications. Processed source data that support the findings of this study are available from the corresponding author on reasonable request.
Abbreviations
- GO:
-
Gene ontology
- IUP:
-
Intrinsically unstructured protein
- MS:
-
Mass spectrometry
- PTM:
-
Post-translational modifications
- SILAC:
-
Stable isotope labeling of amino acids in cell culture
- TAM:
-
Transcription-associated mutagenesis
References
Andersen JS, Lam YW, Leung AK, Ong SE, Lyon CE, Lamond AI, Mann M. Nucleolar proteome dynamics. Nature. 2005;433(7021):77–83.
Belle A, Tanay A, Bitincka L, Shamir R, O'Shea EK. Quantification of protein half-lives in the budding yeast proteome. Proc Natl Acad Sci U S A. 2006;103(35):13004–9.
Cargile BJ, Bundy JL, Grunden AM, Stephenson JL Jr. Synthesis/degradation ratio mass spectrometry for measuring relative dynamic protein turnover. Anal Chem. 2004;76(1):86–97.
Doherty MK, Hammond DE, Clague MJ, Gaskell SJ, Beynon RJ. Turnover of the human proteome: determination of protein intracellular stability by dynamic SILAC. J Proteome Res. 2009;8(1):104–12.
Doherty MK, Whitehead C, McCormack H, Gaskell SJ, Beynon RJ. Proteome dynamics in complex organisms: using stable isotopes to monitor individual protein turnover rates. Proteomics. 2005;5(2):522–33.
Fern EB, Garlick PJ, McNurlan MA, Waterlow JC. The excretion of isotope in urea and ammonia for estimating protein turnover in man with [15N]glycine. Clin Sci (Lond). 1981;61(2):217–28.
Pratt JM, Petty J, Riba-Garcia I, Robertson DH, Gaskell SJ, Oliver SG, Beynon RJ. Dynamics of protein turnover, a missing dimension in proteomics. Mol Cell Proteomics. 2002;1(8):579–91.
Price JC, Guan S, Burlingame A, Prusiner SB, Ghaemmaghami S. Analysis of proteome dynamics in the mouse brain. Proc Natl Acad Sci U S A. 2010;107(32):14508–13.
Yen HC, Xu Q, Chou DM, Zhao Z, Elledge SJ. Global protein stability profiling in mammalian cells. Science (New York, NY). 2008;322(5903):918–23.
Friedel CC, Dolken L, Ruzsics Z, Koszinowski UH, Zimmer R. Conserved principles of mammalian transcriptional regulation revealed by RNA half-life. Nucleic Acids Res. 2009;37(17):e115.
Cambridge SB, Gnad F, Nguyen C, Bermejo JL, Kruger M, Mann M. Systems-wide proteomic analysis in mammalian cells reveals conserved, functional protein turnover. J Proteome Res. 2011;10(12):5275–84.
Shakhnovich BE, Deeds E, Delisi C, Shakhnovich E. Protein structure and evolutionary history determine sequence space topology. Genome Res. 2005;15(3):385–92.
Theis FJ, Latif N, Wong P, Frishman D. Complex principal component and correlation structure of 16 yeast genomic variables. Mol Biol Evol. 2011;28(9):2501–12.
Rattray DG, Foster LJ. Dynamics of protein complex components. Curr Opin Chem Biol. 2019;48:81–5.
Nicholson AL, Pasquinelli AE. Tales of detailed poly(a) tails. Trends Cell Biol. 2019;29(3):191–200.
Ong SE, Foster LJ, Mann M. Mass spectrometric-based approaches in quantitative proteomics. Methods. 2003;29(2):124–30.
Garlick PJ, Waterlow JC. Measurement of muscle protein turnover by constant intravenous infusion of [14C]glycine. Biochem J. 1969;113(2):7P–8P.
Dice JF, Goldberg AL. A statistical analysis of the relationship between degradative rates and molecular weights of proteins. Arch Biochem Biophys. 1975;170(1):213–9.
Capra JA, Williams AG, Pollard KS. ProteinHistorian: tools for the comparative analysis of eukaryote protein origin. PLoS Comput Biol. 2012;8(6):e1002567.
Datta RS, Meacham C, Samad B, Neyer C, Sjolander K. Berkeley PHOG: PhyloFacts orthology group prediction web server. Nucleic Acids Res. 2009;37(Web Server issue):W84–9.
Toll-Riera M, Bostick D, Alba MM, Plotkin JB. Structure and age jointly influence rates of protein evolution. PLoS Comput Biol. 2012;8(5):e1002542.
Chiti F, Dobson CM. Protein misfolding, amyloid formation, and human disease: a summary of progress over the last decade. Annu Rev Biochem. 2017;86:27–68.
De Baets G, Reumers J, Delgado Blanco J, Dopazo J, Schymkowitz J, Rousseau F. An evolutionary trade-off between protein turnover rate and protein aggregation favors a higher aggregation propensity in fast degrading proteins. PLoS Comput Biol. 2011;7(6):e1002090.
Gsponer J, Futschik ME, Teichmann SA, Babu MM. Tight regulation of unstructured proteins: from transcript synthesis to protein degradation. Science (New York, NY). 2008;322(5906):1365–8.
Schad E, Tompa P, Hegyi H. The relationship between proteome size, structural disorder and organism complexity. Genome Biol. 2011;12(12):R120.
Xue B, Dunker AK, Uversky VN. Orderly order in protein intrinsic disorder distribution: disorder in 3500 proteomes from viruses and the three domains of life. J Biomol Struct Dyn. 2012;30(2):137–49.
van der Lee R, Lang B, Kruse K, Gsponer J, Sanchez de Groot N, Huynen MA, Matouschek A, Fuxreiter M, Babu MM. Intrinsically disordered segments affect protein half-life in the cell and during evolution. Cell Rep. 2014;8(6):1832–44.
de Godoy LM, Olsen JV, Cox J, Nielsen ML, Hubner NC, Frohlich F, Walther TC, Mann M. Comprehensive mass-spectrometry-based proteome quantification of haploid versus diploid yeast. Nature. 2008;455(7217):1251–4.
Silva JC, Gorenstein MV, Li GZ, Vissers JP, Geromanos SJ. Absolute quantification of proteins by LCMSE: a virtue of parallel MS acquisition. Mol Cell Proteomics. 2006;5(1):144–56.
Wolf YI, Novichkov PS, Karev GP, Koonin EV, Lipman DJ. The universal distribution of evolutionary rates of genes and distinct characteristics of eukaryotic genes of different apparent ages. Proc Natl Acad Sci U S A. 2009;106(18):7273–80.
Kertesz M, Wan Y, Mazor E, Rinn JL, Nutter RC, Chang HY, Segal E. Genome-wide measurement of RNA secondary structure in yeast. Nature. 2010;467(7311):103–7.
Lian X, Guo J, Gu W, Cui Y, Zhong J, Jin J, He QY, Wang T, Zhang G. Genome-wide and experimental resolution of relative translation elongation speed at individual gene level in human cells. PLoS Genet. 2016;12(2):e1005901.
Hutz JE, Kraja AT, McLeod HL, Province MA. CANDID: a flexible method for prioritizing candidate genes for complex human traits. Genet Epidemiol. 2008;32(8):779–90.
Alba MM, Castresana J. Inverse relationship between evolutionary rate and age of mammalian genes. Mol Biol Evol. 2005;22(3):598–606.
Mann M, Kelleher NL. Precision proteomics: the case for high resolution and high mass accuracy. Proc Natl Acad Sci U S A. 2008;105(47):18132–8.
She X, Rohl CA, Castle JC, Kulkarni AV, Johnson JM, Chen R. Definition, conservation and epigenetics of housekeeping and tissue-enriched genes. BMC Genomics. 2009;10:269.
Chang CW, Cheng WC, Chen CR, Shu WY, Tsai ML, Huang CL, Hsu IC. Identification of human housekeeping genes and tissue-selective genes by microarray meta-analysis. PLoS One. 2011;6(7):e22859.
Vaquerizas JM, Kummerfeld SK, Teichmann SA, Luscombe NM. A census of human transcription factors: function, expression and evolution. Nat Rev Genet. 2009;10(4):252–63.
Liao BY, Zhang J. Null mutations in human and mouse orthologs frequently result in different phenotypes. Proc Natl Acad Sci U S A. 2008;105(19):6987–92.
Liao BY, Zhang J. Mouse duplicate genes are as essential as singletons. Trends Genet. 2007;23(8):378–81.
Lopez-Bigas N, De S, Teichmann SA. Functional protein divergence in the evolution of Homo sapiens. Genome Biol. 2008;9(2):R33.
Vishnoi A, Kryazhimskiy S, Bazykin GA, Hannenhalli S, Plotkin JB. Young proteins experience more variable selection pressures than old proteins. Genome Res. 2010;20(11):1574–81.
Milholland B, Dong X, Zhang L, Hao X, Suh Y, Vijg J. Differences between germline and somatic mutation rates in humans and mice. Nat Commun. 2017;8:15183.
Rutherford SL. Between genotype and phenotype: protein chaperones and evolvability. Nat Rev Genet. 2003;4(4):263–74.
Choudhary C, Kumar C, Gnad F, Nielsen ML, Rehman M, Walther TC, Olsen JV, Mann M. Lysine acetylation targets protein complexes and co-regulates major cellular functions. Science (New York, NY). 2009;325(5942):834–40.
Chen H, Xue Y, Huang N, Yao X, Sun Z. MeMo: a web tool for prediction of protein methylation modifications. Nucleic Acids Res. 2006;34(Web Server issue):W249–53.
Olsen JV, Vermeulen M, Santamaria A, Kumar C, Miller ML, Jensen LJ, Gnad F, Cox J, Jensen TS, Nigg EA, et al. Quantitative phosphoproteomics reveals widespread full phosphorylation site occupancy during mitosis. Sci Signal. 2010;3(104):ra3.
Shao J, Xu D, Tsai SN, Wang Y, Ngai SM. Computational identification of protein methylation sites through bi-profile Bayes feature extraction. PLoS One. 2009;4(3):e4920.
Macek B, Gnad F, Soufi B, Kumar C, Olsen JV, Mijakovic I, Mann M. Phosphoproteome analysis of E. coli reveals evolutionary conservation of bacterial Ser/Thr/Tyr phosphorylation. Mol Cell Proteomics. 2008;7(2):299–307.
Zhang J, Sprung R, Pei J, Tan X, Kim S, Zhu H, Liu CF, Grishin NV, Zhao Y. Lysine acetylation is a highly abundant and evolutionarily conserved modification in Escherichia coli. Mol Cell Proteomics. 2009;8(2):215–25.
Cao XJ, Dai J, Xu H, Nie S, Chang X, Hu BY, Sheng QH, Wang LS, Ning ZB, Li YX, et al. High-coverage proteome analysis reveals the first insight of protein modification systems in the pathogenic spirochete Leptospira interrogans. Cell Res. 2010;20(2):197–210.
Collins MO, Yu L, Campuzano I, Grant SG, Choudhary JS. Phosphoproteomic analysis of the mouse brain cytosol reveals a predominance of protein phosphorylation in regions of intrinsic sequence disorder. Mol Cell Proteomics. 2008;7(7):1331–48.
Schwanhausser B, Busse D, Li N, Dittmar G, Schuchhardt J, Wolf J, Chen W, Selbach M. Global quantification of mammalian gene expression control. Nature. 2011;473(7347):337–42.
Beck M, Claassen M, Aebersold R. Comprehensive proteomics. Curr Opin Biotechnol. 2011;22(1):3–8.
Hendriks G, Jansen JG, Mullenders LH, de Wind N. Transcription and replication: far relatives make uneasy bedfellows. Cell Cycle (Georgetown, Tex). 2010;9(12):2300–4.
Gaillard H, Aguilera A. Transcription as a threat to genome integrity. Annu Rev Biochem. 2016;85:291–317.
Zhu J, He F, Hu S, Yu J. On the nature of human housekeeping genes. Trends Genet. 2008;24(10):481–4.
Monsellier E, Ramazzotti M, Taddei N, Chiti F. Aggregation propensity of the human proteome. PLoS Comput Biol. 2008;4(10):e1000199.
Vogel C, Abreu Rde S, Ko D, Le SY, Shapiro BA, Burns SC, Sandhu D, Boutz DR, Marcotte EM, Penalva LO. Sequence signatures and mRNA concentration can explain two-thirds of protein abundance variation in a human cell line. Mol Syst Biol. 2010;6:400.
Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol. 2008;26(12):1367–72.
Schaab C, Geiger T, Stoehr G, Cox J, Mann M. Analysis of high accuracy, quantitative proteomics data in the MaxQB database. Mol Cell Proteomics. 2012;11(3):M111.014068.
Acknowledgements
Funding
This work was in part funded by a grant to S.C. from the Volkswagen-Stiftung. I acknowledge financial support by Deutsche Forschungsgemeinschaft within the funding program Open Access Publishing, by the Baden-Württemberg Ministry of Science, Research and the Arts and by Ruprecht-Karls-Universität Heidelberg.
Author information
Authors and Affiliations
Contributions
SBC conceived the hypothesis, performed all analyses, and wrote the manuscript. The author read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The author declares that he/she has no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Additional file 1: Figure S1.
Age-dependent changes of protein turnover. (a) Box plot of human protein turnover (Fig. 1a). (b) Median mouse protein turnover for taxonomic groups. Number of proteins: u org=751, u euk=1331, OBD=667, chor=273, mamm=46. (c) Median rat protein degradation rates for taxonomic groups. Number of proteins: u org=188, u euk=461, OBD=558, chor=439, mamm=139. (d) Mean Schizosaccharomyces pombe protein half-lives for taxonomic groups. Number of proteins: u org=409, u euk=915, Ophistokonta (Oph)=407, Ascomycota (Asc)=484, Sch.Po.=737.
Additional file 2: Figure S2.
Age-dependent changes of macromolecules. (a) Median level of protein disorder for taxonomic groups in yeast. Number of proteins: u org=1336, u euk=1693, OBD=562, chor=1192, mamm=1070. (b) Distribution of protein abundance of HeLa MS data set. (c) Median human mRNA abundance for taxonomic groups. u org=190, u euk=337, OBD=208, chor=163, mamm=51. (d) Median mouse mRNA half-lives for taxonomic groups. Number of different mRNA species: u org=1072, u euk=2553, OBD=1636, chor=1122, mamm=313. (e) Mean gene conservation score for taxonomic groups. This graph is based on the same data as Fig. 1h but without the grouping to more general taxonomies. The number of genes is given in parentheses (ANOVA, Bonferroni post-hoc analysis for Suppl. Fig. 1 a,c,d,e).
Additional file 3: Figure S3.
Age-dependent changes of human housekeeping genes. (a) Mean gene conservation score for taxonomic groups. Number of genes: u org=394, u euk=852, OBD=369, chor=237, mamm=39. (b) Median level of protein disorder for taxonomic groups. Number of proteins: u org=385, u euk=823, OBD=366, chor=225, mamm=38. (c) Median protein turnover for taxonomic groups. Number of proteins: u org=303, u euk=588, OBD=167, chor=85, mamm=10. (d) Median mRNA half-lives for taxonomic groups. Number of different mRNA species: u org=360, u euk=793, OBD=341, chor=220, mamm=40. (ANOVA, Bonferroni post-hoc analysis for all figures).
Additional file 4: Figure S4
. Age-dependent changes of human transcription factors and mouse essential genes. (a) Median protein turnover for taxonomic groups of transcription factors. Number of proteins: u org+ u euk=34, OBD=34, chor+ mamm=23. Groups with less than 10 genes were added to the neighboring group (also in b,c,d) (b) Median mRNA half-lives for taxonomic groups of transcription factors. Number of different mRNA species: u org+ u euk=51, OBD=172, chor=175, mamm=89.(c) Mean gene conservation score for taxonomic groups of transcription factors. Number of genes: u org+u euk=79, OBD=471, chor=509, mamm=241. (d) Median level of protein disorder for taxonomic groups of transcription factors. Number of proteins: prok+ u euk=82, OBD=531, chor=512, mamm=322. (e) Median protein turnover for taxonomic groups of mouse essential genes. Number of proteins: u org=92, u euk=168, OBD=129, chor+ mamm=73. (f) Median mRNA half-lives for taxonomic groups of mouse essential genes. Number of different mRNA species: u org=133, u euk=296, OBD=327, chor=222, mamm=13. (ANOVA, Bonferroni post-hoc analysis for all figures).
Additional file 5: Table S3
and Figure S5. Table S3: Median protein turnover values for GO categories based on mass spectrometry proteomics data from arrested human HeLa and differentiated mouse muscle C2C12 cells. Probability indicates difference to average turnover of all proteins. Figure S5: Evolutionary changes of protein turnover between human and mouse. (a) Graphical representation of Supplementary Table 3. (b) Median peptide turnover values of the ‘histone deacetylase and nucleosome remodeling activities complex’ proteins. (c) Grouping of proteins into equal size bins of 210 and quantification of homologous proteins with a turnover increase from human to mouse.
Additional file 6: Tables S1
and S2. Median protein turnover values for GO categories based on mass spectrometry proteomics data from arrested human HeLa and differentiated mouse muscle C2C12 cells. Probability indicates difference to average turnover of all proteins.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Cambridge, S.B. Hypothesis: protein and RNA attributes are continuously optimized over time. BMC Genomics 20, 1012 (2019). https://doi.org/10.1186/s12864-019-6371-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12864-019-6371-0