Comparative genomics of nucleotide metabolism: a tour to the past of the three cellular domains of life
BMC Genomics volume 15, Article number: 800 (2014)
Nucleotide metabolism is central to all biological systems, due to their essential role in genetic information and energy transfer, which in turn suggests its possible presence in the last common ancestor (LCA) of Bacteria, Archaea and Eukarya. In this context, elucidation of the contribution of the origin and diversification of de novo and salvage pathways of nucleotide metabolism will allow us to understand the links between the enzymatic steps associated with the LCA and the emergence of the first metabolic pathways.
In this work, the taxonomical distribution of the enzymes associated with nucleotide metabolism was evaluated in 1,606 complete genomes. 151 sequence profiles associated with 120 enzymatic reactions were used. The evaluation was based on profile comparisons, using RPS-Blast. Organisms were clustered based on their taxonomical classifications, in order to obtain a normalized measure of the taxonomical distribution of enzymes according to the average of presence/absence of enzymes per genus, which in turn was used for the second step, to calculate the average presence/absence of enzymes per Clade.
From these analyses, it was suggested that divergence at the enzymatic level correlates with environmental changes and related modifications of the cell wall and membranes that took place during cell evolution. Specifically, the divergence of the 5-(carboxyamino) imidazole ribonucleotide mutase to phosphoribosylaminoimidazole carboxylase could be related to the emergence of multicellularity in eukaryotic cells. In addition, segments of salvage and de novo pathways were probably complementary in the LCA to the synthesis of purines and pyrimidines. We also suggest that a large portion of the pathway to inosine 5’-monophosphate (IMP) in purines could have been involved in thiamine synthesis or its derivatives in early stages of cellular evolution, correlating with the fact that these molecules may have played an active role in the protein-RNA world. The analysis presented here provides general observations concerning the adaptation of the enzymatic steps in the early stages of the emergence of life and the LCA.
Metabolism represents an intricate ensemble of enzyme-catalyzed reactions that lead to synthesis and degradation of compounds within the cell. In recent years, an increasing amount of information on metabolism from different species has become available, allowing for comparative genomic-scale studies on the evolution of specific pathways or whole metabolic networks [1–4]. Metabolism can be considered one of the most ancient biological networks, where nodes represent substrates or enzymes and edges represent the relationships between them. From this perspective, the study of metabolic networks is focused on describing topological properties, such as the existence of functional modules, giving special relevance to clustering and motif formation, and showing the existence of similar attributes to the small-world and scale-free networks .
Therefore, two main hypotheses on the origin and evolution of enzyme-driven metabolism have been postulated based on the notion that gene duplication, followed by divergence, can lead to the origin of new metabolic reactions. The “stepwise hypothesis”  suggests that, in the case where a substrate tends to be depleted, gene duplication can provide an enzyme capable of supplying the exhausted substrate, giving rise to homologous enzymes that catalyze consecutive reactions. On the other hand, the “patchwork hypothesis”  proposes that duplication of genes encoding promiscuous enzymes (capable of catalyzing multiple reactions) allows each descendant enzyme to specialize in one of the ancestral reactions.
In this regard, it is plausible that a small number of enzymes with broad specificity existed in early stages of metabolic evolution. Genes encoding these enzymes would have been duplicated, generating enzymes that, through sequence divergence, became more specialized .
Collectively, these studies have highlighted the contribution of gene duplication in the evolution of metabolism . In recent works, the universal occurrence of some pathways and branches, such as diverse amino acid pathways, in modern species suggests that they existed in the last common ancestor (LCA) of Bacteria, Archaea and Eukarya[1, 8]. However, despite the importance of nucleotide metabolism in all organisms, few studies have addressed this issue by using genomic approaches [9, 10]. Nucleotide metabolism is central in all living systems, due to its role in transferring genetic information and energy. Indeed, it has been described as one of the ancient metabolisms in evolution. Specifically, the emergence of an ancestral folding or P-loop hydrolase appeared parallel to this metabolism , reinforcing its antiquity. In addition, many of the intermediates associated with this metabolic module have been intimately associated with prebiotic chemistry and the origin of life [10, 12, 13]. In this regard, we adopted a multigenomic strategy for the reconstruction and analysis of the metabolism of nucleotides, evaluating the contribution of the origin and diversification of de novo and salvage pathways for nucleotides in the evolution of organisms. In addition, these analyses allow the identification of a metabolic link between the LCA and the first steps in the structure of biological networks [14–16]. Our strategy reveals some general rules concerning the adaptation of the first predominant chemical reactions to enzymatic steps in the LCA and allows us to infer environmental issues in the early stages of the emergence of life. In addition, it was possible to determine the presence of de novo biosynthesis pathways of the ribonucleotides and deoxyribonucleotides of uracil and cytosine for pyrimidine metabolism associated with the LCA. Finally, we found differences in the enzymes for nucleotide metabolism that correlated with environmental changes and with associated cellular architecture adaptations; such is the case for the enzymes involved in the synthesis of 5-aminoimidazole ribonucleotide (AIR) to (4-carboxyaminoimidazole ribonucleotide) (CAIR). These findings suggest the presence of HCO3- in primitive seas and its use as one of the main carbon sources for the first organisms. Our main results derived from the taxonomical distribution of enzymatic families belonging to nucleotide metabolism, supported by experimental evidence, are described in this report.
Results and discussion
Taxonomic distribution of nucleotide metabolic enzymes
The taxonomic distribution of proteins provides clues concerning the relative occurrence of the enzymes and their branches and paths for the evolution of metabolism [1, 8, 11, 17, 18]. In this regard, the origin and evolution of nucleotide metabolism was traced in organisms belonging to the three cellular domains, Bacteria, Archaea and Eukarya, by an exhaustive evaluation of the taxonomic distribution of their enzymatic repertoires. Therefore, each enzymatic activity encoded by an EC number can exhibit different profiles, which are individual vectors of the presence/absence of homologous enzymes in all the genomes. In total, 151 profiles associated with 120 enzymatic reactions related to nucleotide metabolism were used to scan 1,606 genomes by using RPS-BLAST. (Additional file 1: Figure S1 and Additional file 2: Figure S2 and Additional file 3: Tables S1 and Additional file 4: Table S2). Based on these comparisons, evolutionary origins of nucleotide metabolism can be traced, even close to the LCA of all organisms. In this regard, enzymes widely distributed across the three cellular domains are proposed to have been present in the LCA [8, 11, 18–20]. Alternatively, enzymes constrained to specific clades or cellular domains would suggest adaptations of organisms or cellular domains to specific lifestyles. Based on these considerations, we discuss our most notable results in the following sections.
Evolution of purine metabolism
De novo purine biosynthesis
The de novo biosynthesis of purines, starting from d-ribose-1-phosphate to inosine 5’-monophosphate (IMP) production, the main intermediate in the synthesis of ribonucleotides and deoxyribonucleotides, guanine and adenine, follows a linear branch. The first step is associated with phosphoglucomutase (EC 126.96.36.199) or phosphopentomutase (188.8.131.52) and the second step is associated with ribose-phosphate diphosphokinase (184.108.40.206); both steps are necessary for the synthesis of 5-phospho-alpha-d-ribosy-1-pyrophosphate (PRPP), which starts from d-ribose-1-phosphate. Based on their taxonomical distribution, the enzymes associated with the 220.127.116.11 and 18.104.22.168 catalytic activities were identified as being widely distributed among Bacteria, Archaea and Eukarya, suggesting the probable existence of PRPP biosynthesis in the LCA. Indeed, PRPP is a key precursor for biosynthesis in the de novo and salvage pathways for purines and pyrimidines; however, this intermediary is unstable and susceptible to hydrolysis. Therefore, it is probable that its abiotic synthesis, if it occurred, was not enough to maintain the biosynthesis in the LCA .
Therefore, the first step for purine biosynthesis, the catalysis to ribose-5-phosphate starting from ribose 1-phosphate, is achieved by either of the two enzymes related to the enzymes EC 22.214.171.124 and 126.96.36.199. These two enzymes are analogous, since no homology at the sequence or structural level was detected. The enzyme 188.8.131.52 is partially distributed in Bacteria, mainly in free-living organisms associated with a host, such as Streptococcus pneumoniae and Lactobacillus rhamnosus; however, it was not found in archaeal and eukaryal organisms, suggesting that its emergence was posterior to the LCA divergence, probably as a secondary adaptation associated with the bacterial host.
Starting from the intermediary PRPP, in the linear branch towards IMP biosynthesis, we identified enzymes belonging to five catalytic steps (amidophosphoribosyltransferase, 184.108.40.206; phosphoribosylamine-glycine ligase, 220.127.116.11; phosphoribosylglycinamide formyltransferase, 18.104.22.168; phosphoribosylformylglycinamidine synthase, 22.214.171.124; phosphoribosylformylglycinamidine cyclo-ligase, 126.96.36.199) required for the transformation of PRPP into AIR . Most of these enzymes were identified as widely distributed in the three cellular domains, suggesting their presence in the LCA (Figure 1). The enzymatic step associated with EC 188.8.131.52 is responsible for the transformation of glycinamide ribotide (GAR) to formyglycinamide ribotide (FGAR) and could be carried out by two enzymes associated with different evolutionary families, PurN (Figure 1 Gold box) or phosphoribosylglycinamide formyltransferase and PurT or phosphoribosylglycinamide formyltransferase 2 . Proteins associated with the PurN family use derivatives from folate synthesis as substrates. This family was identified as widely distributed in Bacteria and Eukarya and partially distributed in Archaea. Alternatively, proteins from the PurT family were partially distributed in Archaea and Bacteria and sparsely in Eukarya. It is probable that the PurT enzymatic family could have been present in the LCA, with posterior loss events in Eukarya due to its requirement for formate as a substrate. In this regard, formate is described as one-carbon donor and one of the main molecules present in prebiotics conditions, prior to folate metabolism [22–24]. In a posterior phase, the emergence of folate biosynthesis might have facilitated the emergence of PurN (Figure 1, gold box), thereby achieving the co-occurrence of PurN and PurT in the LCA. Indeed, previous works have suggested that the emergence of PurT preceded the emergence of PurN, mainly because PurT utilizes a more primitive substrate prior to the folate-dependent pathway . One of the evolutionary pressures for the selection of PurN instead of PurT in eukaryotic organisms could be associated with the emergence of the mitochondrial respiratory chain. It has been shown that the PurT substrate, formate, is toxic and binds to cytochrome c oxidase-like [25, 26], uncoupling the redox reactions and favoring the selection of PurN in eukaryotic organisms.
Another enzymatic step identified in the purine biosynthetic pathway is achieved by the phosphoribosylformylglycinamidine synthase (184.108.40.206) that transforms FGAR to formylglycinamidine-ribonucleotide (FGAM). This enzyme is composed of two catalytic subunits, PurQ and PurL. In general, this enzyme may occur as a multidomain protein or as a multiprotein complex, where the subunits L and Q are encoded independently. Altogether, we identified a co-occurrence of these subunits in the three cellular domains. Finally, the last step towards AIR transformation is carried out by the cyclo-ligase phosphoribosylformylglycinamidine (EC 220.127.116.11), which transforms FGAM to AIR, and in turn this enzyme was found to be widely distributed in the three cellular domains, suggesting its presence in the LCA.
The product of the steps previously described, the AIR 5-amino-1-(5-phospho-d-ribosyl) imidazole-4-carboxylate, can be transformed by two pathways; in the first one, only one reaction, catalyzed by 18.104.22.168 (AIR carboxylase) (Figure 1, green box), is involved in CAIR synthesis; in the second, a two-step pathway, entails two reactions involving N5-carboxyaminoimidazole ribonucleotide synthase (22.214.171.124) and 5-(carboxyamino) imidazole ribonucleotide mutase (126.96.36.199). When the taxonomical distribution of these enzymes was evaluated, it was found that the N5-carboxyaminoimidazole ribonucleotide synthase (188.8.131.52) was widely distributed among all the organisms in the three cellular domains, whereas the enzyme 5-(carboxyamino) imidazole ribonucleotide mutase (184.108.40.206) exhibited a wide distribution in Bacteria and Archaea and in some eukaryal clades, suggesting their ancestry in all the organisms with posterior gene loss in Eukarya (Figure 1). Additionally, these enzymes were analyzed at the sequence level, where the 5-(carboxyamino) imidazole ribonucleotide mutase (220.127.116.11) was evolutionarily related to AIR carboxylase (18.104.22.168) . Indeed, at the structural level both enzymes share a common evolutionary origin, belonging to the N5-CAIR mutase (phosphoribosylaminoimidazole carboxylase superfamily). The observation that AIR carboxylase (22.214.171.124) is almost exclusively distributed in the Eukarya suggests its posterior subfunctionalization from 5-(carboxyamino) imidazole ribonucleotide mutase (126.96.36.199) .
The discovery of the emergence of these branches is interesting, because it allows us to infer the processes associated with environmental changes on Earth. In this context, Tribunskikh et al. suggested that the divergence of the phosphoribosylaminoimidazole carboxylase (188.8.131.52) to 5-(carboxyamino) imidazole ribonucleotide mutase (184.108.40.206) was a consequence of decreasing atmospheric CO2, which resulted in the addition of the 5-(carboxyamino) imidazole ribonucleotide synthase (220.127.116.11) and a change of specificity and, by consequence, in the two-step pathway. Those authors supported their proposal by the fact that the two-step pathways could work at low CO2 levels, under aerobic or anaerobic conditions. However, our data suggest an alternative scenario where, although the concentration of CO2 in the primitive atmosphere was high, acidic conditions in early oceans favored the formation of HCO3-. This is consistent with simulations showing that the concentration of HCO3- oscillated between levels 30 to 30,000 times higher in the early oceans than current levels . Interestingly, the enzyme 18.104.22.168, which is associated with the two-step pathway, uses HCO3- as a substrate, and HCO3- is considered a dominant form of CO2 in early oceans [30, 31], which together with our data of taxonomic distribution supports the notion that this pathway preceded the path where the enzymatic reaction catalyzed by 22.214.171.124 uses CO2 as the substrate. Indeed, it has been experimentally shown that in aqueous solutions with high concentrations of KHCO3, AIR is easily converted into CAIR in the absence of enzymes [29, 32, 33]. This mechanism could take place through the accumulation of the intermediate N-CAIR (N5-carboxy-AIR, or N5-carboxyaminoimidazole ribonucleotide), which then undergoes a rearrangement to CAIR. These reactions appear to be the template whereby the enzymatic activities adapted in the two-step pathway, widespread in the cellular domains and probably in the LCA. Probably, HCO3- levels had a steady decline due to the reduction of atmospheric CO2, as documented in the evolution of the terrestrial atmosphere . This reduction could have selected the enzymes 126.96.36.199 and 188.8.131.52, which transform AIR to N-CAIR and N-CAIR to CAIR, respectively.
Subsequently, once the atmosphere was provided with oxygen, the emergence of mitochondria and eukaryotic cells was possible. In this regard, HCO3- is one of the main products of mitochondrial respiration, which follows a pH-dependent conversion to CO2, converting an impermeable anion into a gas that can diffuse through membranes . The transformation of HCO3- accumulated in cells with mitochondrial activity towards CO2 may have resulted in a more efficient regulation of intracellular pH in eukaryotic cells, and in parallel, to the use of CO2 as a carbon source. All these processes might have favored multicellularity, because cells with nutrients availability and high mitochondrial activity as the consequence of oxidative respiration could provide a permeable carbon source to other cells with low nutrients availability. Altogether, the taxonomic distribution data, chemical synthesis information, and primitive ocean simulations data, support the divergence of the enzyme 184.108.40.206 to 220.127.116.11 by selecting a CO2 binding site. Although the selection of a CO2 binding site in metabolism is not exclusive to this protein family, the selection of a CO2 binding site in these enzymes could favor the development of multicellularity in eukaryotic cells.
The last two steps in the synthesis of IMP correspond to the transformation of 5-aminoimidazole-4-carboxamide ribonucleotide (AICAR) and 5-formamidoimidazole-4-carboxamide ribotide (FAICAR). In the first step (18.104.22.168), AICAR formyltransferase activity is required, whereas for the second step (22.214.171.124), IMP cyclohydrolase activity is required. The gene for the bifunctional folate-dependent enzyme PurH, which is widely distributed in Bacteria and Eukarya and partially in Archaea, encodes both activities . PurH exhibits two independent catalytic sites, in which each half of the enzyme catalyzes an independent reaction. The transformation to IMP by the bifunctional enzyme PurH may have been posterior to the emergence of the folate synthesis pathway (Figure 1, gold box). In addition, both catalytic reactions performed by PurH could have emerged in Archaea to be carried out by different proteins that achieve each step independently. In this regard, the first reaction is carried out by the formate-dependent PurP enzyme, whereas the second one is achieved by the IMP cyclohydrolase PurO (Figure 1, magenta). In both cases, the catalytic mechanism is very similar to that of its counterpart PurH, and they are exclusively distributed in Archaea; however, these two enzymes do not maintain relationships at the sequence or structural level with PurH [9, 35]. In summary, all these results suggest that PurP and PurO emerged posterior to archaeal divergence.
Alternatively, the archaeal cenancestor may have brought the formate (PurT) and folate-dependent (PurN-PurH) families, as seen in Halobacteria. The emergence of PurP in Archaea was probably related to the availability of formate in the environment, which agrees with the fact that most of the Archaea clades are methylotrophs. In turn, PurP may have duplicated from families associated with the LCA [carbamoyl-phosphate synthase, 126.96.36.199; phosphoribosylamine-glycine ligase, 188.8.131.52; 5-(carboxyamino) imidazole ribonucleotide synthase, 184.108.40.206; phosphoribosylglycinamide formyltransferase, 220.127.116.11] belonging to nucleotide metabolism, partially replacing PurH in some archaeal clades. For instance, PurH was identified in some archaeal clades, such as Halobacteria, Thermoplasmata, and Methanomicrobia among others, including co-occurrence with PurP, suggesting that their activities are complementary in these organisms.
In the case of PurO, we did not find homology relationships with any other enzymatic family associated with nucleotide biosynthesis. In addition, we found a correlation between the presence of PurP and PurO; however, in Archaea some clades do not exhibit either PurH or PurO, which suggests that alternative enzymes performing the PurO function are probably present and remain to be described.
Purine salvage pathway
Ribonucleotide adenine salvage pathway
In the purine salvage pathway, diverse pathways to generate adenine and guanine ribonucleotides and deoxynucleotides were identified. In this regard, the first pathway starts from hypoxanthine, via hypoxanthine-guanine phosphoribosyltransferase (PRTase) (18.104.22.168) (Figure 2), an enzyme with broad specificity and that is universally distributed in the three cellular domains. This enzyme achieves the biochemical transformation to IMP, starting from hypoxanthine and PRPP, for subsequent transformation to the ribonucleotide of adenine through the de novo adenyl-succinate pathway, which includes enzymes widely distributed in the three cellular domains (22.214.171.124, adenylosuccinate synthetase; 126.96.36.199, adenylosuccinate lyase; 188.8.131.52, adenylate kinase; 184.108.40.206, nucleoside diphosphate kinase; 220.127.116.11, pyruvate kinase), suggesting its presence in the LCA (Figure 2). For the transformation of alpha-deoxyribonucleotide, there are two routes, the pathway associated with the ribonucleoside-diphosphate reductase beta-subunit (small subunit) (18.104.22.168), the ribonucleoside-diphosphate reductase alpha-subunit (large chain family), and the 22.214.171.124 pathway, which contains adenosylcobalamin-dependent ribonucleoside-triphosphate reductase, and the anaerobic ribonucleoside-triphosphate reductase complex (NrdD and NrdG). The first of these pathways, for the formation of deoxyribonucleotides, is widely distributed in the three cellular domains, suggesting its presence in the LCA. Finally, the enzymes of the 126.96.36.199 pathway are partially distributed in Bacteria and Archaea but not in Eukarya, which makes it difficult to determine its presence in the LCA.
Alternatively, it is also possible to synthesize adenine deoxynucleotides and ribonucleotides. This pathway conserves the last three catalytic steps (188.8.131.52, 184.108.40.206 and 220.127.116.11) from the pathway previously described, i.e., only the first step has been added, which is catalyzed by adenine phosphoribosyltransferase (a PRTase) (18.104.22.168) (Figure 2), an enzyme widely distributed in Bacteria and Eukarya and that starts from the adenine salvage pathway.
This enzyme may be specialized for the adenine substrate, posterior to the divergence of the LCA, from a broad-specificity ancestor, similar to the HPRTases (22.214.171.124).
Route of salvage of guanine ribonucleotides
Starting from the first salvaging route, previously described, which starts with hypoxanthine-guanine phosphoribosyltransferase (126.96.36.199), it is also possible to synthesize guanine ribonucleotide, adding to the enzymatic families widely distributed in the three cellular domains: inosine-5’-monophosphate dehydrogenase 1 (188.8.131.52) and GMP synthase (glutamine-hydrolyzing) subunit A and GMP synthase (glutamine-hydrolyzing) subunit B (184.108.40.206) (Figure 2). Because the hypoxanthine-guanine phosphoribosyltransferase (220.127.116.11) also exhibits specificity for guanine, it is possible to synthesize (ribonucleotide monophosphate guanine) guanosine 5’-monophosphate (GMP) in one step, starting from the guanine salvage pathway with PRPP; however, although the subsequent step with 18.104.22.168, performed by the guanylate kinase family, for transforming GMP to GDP is widely distributed in Eukarya and Bacteria, in Archaea only the reaction but not the enzyme has been identified (Figure 2, star) . These data limit the possibility of making genomic and evolutionary comparisons between Bacteria and Eukarya to guarantee the possible presence of guanine deoxyribonucleotides in the LCA. However, subsequent steps for the synthesis of GDP to GTP are carried out with the nucleoside diphosphate kinase (22.214.171.124) and pyruvate kinase (126.96.36.199); both of these enzymes are widely distributed in the three cellular domains, and the step of GTP conversion to dGTP for the widely distributed enzyme ribonucleoside diphosphate reductase (188.8.131.52) suggests the presence of GTP and dGTP in the LCA.
Therefore, although the gene sequence associated with the guanylate kinase function in Archaea has not yet been identified, we suggest that this protein could exhibit a common origin with its counterpart in Bacteria and Eukarya, because most of the nucleotide kinases belong to the P-loop-containing nucleoside triphosphate hydrolases family. This family shows a broad specificity in recognition of ribonucleotides and deoxyribonucleotides . In this regard, it has been reported that only two mutations are enough to introduce the adenylate kinase activity into guanylate kinase , suggesting a masking of guanylate kinase assignment by adenylate kinase. In addition, it is possible that an enzyme with broad specificity similar to adenylate kinase (Figure 3, star) recognizes NMP and converts it into NDP; thus, its specialization toward the guanylate kinase function occurred posterior to the divergence of the three cellular domains.
At present, members of the adenylate kinase family and in particular from the AK6 subfamily have been described as proteins that exhibit the ability to transform different NMPs to their corresponding NDPs, specifically, the AK6 subfamily transforms AMP, dAMP, CMP, dCMP, IMP, and GMP to their corresponding NDPs and dNDPs [39, 40]. This finding also suggests that an ancestral enzyme with broad specificity for the synthesis of ribonucleotides and deoxynucleotides could have been utilized in the LCA. It is also interesting that this enzyme could have used IMP as a substrate, since this nucleotide could have had an active role in the evolution of nucleic acids and the genetic code . Interestingly, our results show that Archaea contain an adenylate kinase from the AK6-like subfamily, similar to Eukarya and Bacteria, suggesting that the function of guanylate kinase in archeal organisms, which has not yet been identified, could be carried out at least in part by enzymes of this subfamily.
Integrating the evolutionary analysis of the salvage and de novo purine pathways
The taxonomic distribution of enzymes associated with purine metabolism shows that segments of the de novo and salvage pathways were complementary and critical to the availability of nucleic acids before the divergence of the three cellular domains. This finding correlates with the chemical synthesis under prebiotic conditions for nitrogenous base precursors for the purine salvage pathway, as previously described [10, 13, 42, 43] (Figure 3, blue circles). The findings of those previous studies agree with the fact that these routes were dependent on PRPP biosynthesis, since it has been argued that PRPP, nucleosides, and nucleotides are susceptible to hydrolysis and thus are very unlikely prebiotic compounds. In addition, our data show a widely taxonomic distribution in the three cellular domains of the enzymes associated with the last two enzymatic steps (184.108.40.206 and 220.127.116.11) required for the transformation of ribose-1-phosphate to PRPP (Figure 3, yellow circle), supporting the biosynthesis of this molecule in the LCA.
Despite the possible absence of folate and folate-dependent enzymes PurH and PurN in early stages of life emergence (Figure 1, gold boxes), most of the enzymes associated with the de novo biosynthesis pathway for IMP may have occurred in the LCA (Figure 1), based on their wide taxonomical distribution in the three cellular domains. Posterior to early stages of life emergence and folate biosynthesis, the appearance of PurH could have completed the purine biosynthesis pathway, filling the gap between AICAR and IMP (Figure 1, PurH is shown in the gold boxes). In addition, the fact that PurP and PurO have been identified exclusively in Archaea, without homology with PurH, suggests that these proteins were generated posterior to the divergence of Archaea (Figure 1, magenta boxes).
Previous to the emergence of folate synthesis and by consequence of PurH, it is possible that in addition to the salvage routes, semienzymatic alternative routes, starting from the substrates 5-amino-4-imidazole carboxyamide AICA and PRPP, are transformed to AICAR by an ancestral phosphoribosyltransferase (PRTase), with subsequent nonenzymatic transformations to IMP, as has been previously proposed  (Figure 3, white boxes). Specifically, the adenine phosphoribosyltransferase (PRTases) (18.104.22.168) synthesize adenine ribonucleotides from the adenine salvage pathway and also transform AICA to AICAR. This enzyme may have become specialized posterior to the divergence of the LCA, from a broad-specificity ancestor, similar to the HGPRTase (PRTases) (22.214.171.124); indeed, this enzyme was found to be universally distributed among the three cellular domains and exhibits a similar catalytic mechanism. In evolutionary terms, the enzymes 126.96.36.199 and 188.8.131.52 belong to the phosphoribosyltransferases superfamily, suggesting a common evolutionary origin.
Based on data previously described, we suggest that an ancestral PRTase enzyme similar to the hypoxanthine-guanine phosphoribosyltransferase (184.108.40.206) (Figure 3, star) with broader specificity to diverse structurally related substrates, such as guanine, xanthine, hypoxanthine, adenine, and AICA, was present in the early stage of the emergence of life. This ancestral enzyme with prebiotic origins could have been involved in the transformation of AICA (Figure 3, light blue circle)  to AICAR, not only to feed semienzymatic IMP synthesis, as previously proposed  (Figure 3, white box), but also to a greater extent to thiamine (vitamin B1) biosynthesis. From this, it is interesting two pathways feed that thiamine metabolism as a bifurcation of the de novo IMP pathway and also with AICA in backflow as a seed substrate (Figure 3, following the direction of the dashed arrows). Additionally, the thiamine pathway may have been also fed from a third source by means of histidine metabolism (Figure 3, box with bars); the connection to purine biosynthesis results from an enzymatic step catalyzed by imidazole glycerol phosphate (IGP) synthase, which transforms N-(5-phosphoribosyl)-formimino-5-aminoimidazol-4-carboxamide ribonucleotide (PRFAR) into AICAR, which is then recycled into the de novo purine biosynthetic pathway, and imidazole-glycerol 3-phosphate, which leads to histidine. Interestingly, previous works have also suggested that the histidine biosynthetic route is ancient and related to the emergence of life . Our proposal of an ancestral branch related to thiamine synthesis is consistent with the catalytic functions that have been proposed for this molecule in the early evolution of life, as suggested by its essential catalytic role in most of organisms and its requirement at several central points of anabolic and catabolic intermediary metabolism , as well as in semienzymatic pathways that may precede the actuals . It is interesting this branch in the early emergence of life and before the constitution of the LCA could have been fed the thiamine metabolism in a semienzymatic way, as the AICAR-to-AIR transformation could occur in a facile, nonenzymatic chemical synthesis pathway .
In this regard, molecules of thiamine or its derivatives bind to the mRNA in the absence of cofactors or proteins in the three domains of life, forming a complex that sequesters the ribosome binding site and suggesting the existence of an ancestral form of riboswitches, which have been implicated in regulatory mechanisms . Additionally, thiamine could have interacted with the RNA, leading to catalytically versatile ribozymes related to the RNA world, due to its catalytic and RNA binding capabilities .
Finally, we suggest that the guanylate kinase function could have been carried out in the LCA by an ancestral enzyme with broad specificity to nucleoside monophosphates (NMPs), similar to those in the AK6 adenylate kinase subfamily. Its specialization to guanylate kinase could have been posterior to archaeal divergence. In parallel, this ancestral enzyme may have provided the nucleotide inosine (Figure 3, pink circle), which has been suggested to play an active role in a rudimentary stage of the genetic code . The ITP can compete with or replace the ATP and GTP binding sites in diverse proteins, such as RNA polymerase [48, 49], that may have maintained their affinity and specificity for ITP as a remnant of its ancestral role.
In addition, inosine maintains a strong structural similarity to guanine , and it may even have the same function for guanine pairing to cytosine in the codons. Currently, inosine has been found in the third position of anticodons, pairing with codons at bases U, C or A and thereby decreasing the need for the 61 tRNAs for each codon. In this regard, it has been suggested that inosine may have been produced by adenosine deamination or even in RNA-mediated catalysis through an early stage of emergence of the genetic code, and it was excluded in nucleic acids when the canonical Watson-Crick pairing evolved to avoid ambiguous rules in replication . Posterior to the divergence of Archaea from the LCA, the specialization and divergence to guanylate kinase of members with broad specificity of the ancestral AK6 subfamily type could have increased the availability of guanine, favoring the replacement of inosine for guanine.
Evolution of pyrimidine metabolism
De novo pyrimidine biosynthesis
Based on taxonomical distributions, we evaluated the enzymes associated with de novo pyrimidine biosynthesis: carbamoyl-phosphate synthase large chain, carbamoyl-phosphate synthase small chain, 220.127.116.11; aspartate carbamoyltransferase 1, 18.104.22.168; dihydroorotase, DHOase family, 22.214.171.124; dihydroorotate dehydrogenase B (NAD+), catalytic subunit, dihydroorotate oxidase B, electron transfer subunit, 126.96.36.199; dihydroorotate dehydrogenase (quinone), 188.8.131.52; dihydroorotate dehydrogenase A (fumarate), 184.108.40.206; orotate phosphoribosyltransferase 1, PyrE1, 220.127.116.11; orotidine 5’-phosphate decarboxylase, 18.104.22.168; cytidylate kinase, 22.214.171.124; nucleoside diphosphate kinase, 126.96.36.199 (Figure 4). These enzymes form the branch of UTP de novo biosynthesis starting from l-glutamine. This entire branch may have occurred in the universal ancestor, based on the wide taxonomical distribution of the enzymes that compose it (Figure 4, red boxes). One of the most interesting steps of this pathway is the conversion of dihydroorotate to orotate, which is carried out by enzymes classified into two families of dihydroorotate group dehydrogenases (DHODs), according to the terminal electron acceptor and relationships at the sequence level . In this regard, family 1 uses soluble electron acceptors. This family is widely distributed in gram-positive Bacteria, Archaea and in some unicellular eukaryotic organisms. In turn, this family is subdivided into DHODA (188.8.131.52) and DHODB (184.108.40.206), which are homodimeric and heterotetrameric, respectively. DHODA uses fumarate as the electron acceptor, whereas DHODB uses NAD+. Members of family 2 (220.127.116.11) are linked to the cell membrane and use quinones from the respiratory chain as electron acceptors . These enzymes are mainly found in most eukaryotic organisms and gram-negative bacteria, in agreement with our taxonomical distribution results.
In spite of these differences, both DHODs belong to the FMN-linked oxidoreductases superfamily, suggesting a common ancestor. Probably, this ancestral enzyme was similar in functional terms to members of the family 2 (Figure 4, star), according to its electron acceptor molecules, which have been previously described as one of the most abundant in extraterrestrial environments, and also that could have been delivered to the Earth around 4 billons years ago, making possible the prebiotic conditions needed for the emergence of life. Indeed, quinones have been found in meteorites in considerable amounts and also have been synthesized in a cloud chamber simulation with good yields [54, 55]. Additionally, quinones are spontaneously partitioned into model membrane systems, representing an evolutionary advantage to early organisms by providing some protection against UV radiation in the early Earth environment , and were posteriorly exploited for their capacity to pump protons across membrane bilayers . In a posterior step, changes associated with the cell membrane and cell wall could lead to divergence of the DHODs from family 2 to family 1, via incorporation of soluble electron acceptors. Surprisingly, our proposal for DHOD family divergence is consistent with previous reports describing a transition from gram-negative to gram-positive bacteria [58, 59] and to drastic changes associated with the chemical constituents of the cell membrane in Archaea, such as glycerol stereochemistry , posterior to the divergence of the LCA. Furthermore, subdivision of family 1 might have occurred in the direction of DHODB to DHODA, since most of these enzymes have been found in gram-positive bacteria that have adapted to parasitic or symbiotic relationships, as shown by our taxonomical distribution results. The DHODA uses fumarate as a terminal electron acceptor, which in turn is used instead of oxygen as a terminal electron acceptor for succinate production, one of the essential processes controlling redox homeostasis for many parasites living under anaerobic conditions .A wide taxonomic distribution of the enzyme CTP synthase (18.104.22.168), which converts UTP to CTP, was also found in the three cellular domains, in addition to the enzyme for the subsequent steps, nucleoside diphosphate kinase (22.214.171.124), thioredoxin reductase NTRA (126.96.36.199) and ribonucleotide reductase (188.8.131.52), giving rise to deoxy-CTP. Based on the distributions of these enzymes, the presence of UTP, CTP and deoxy-CTP in the LCA is suggested (Figure 5, red boxes).
Concerning the transformation to deoxyuridine, two main routes have been described. One pathway starts from deoxy-CDP, whereas the second one starts from deoxy-CMP and which can be extended to transforming deoxyuridine. The route of deoxy-CDP requires two enzymatic steps: the first one requires the nucleoside diphosphate kinase (184.108.40.206), which produces dUTP, which is subsequently transformed by the inosine triphosphate pyrophosphatase 1 (nucleoside-triphosphate pyrophosphatase 1; 220.127.116.11) to dUMP. These two enzymes were identified as universally distributed in the three cellular domains; therefore, they can be associated with the LCA (Figure 5, red boxes). The second pathway converts dCMP to dUMP via cytidine/deoxycytidylate deaminase, which is associated with the catalytic activity of 18.104.22.168, whose members are partially distributed among the three cellular domains. Finally, a third pathway, which starts from dCTP and is catalyzed by deoxycytidine triphosphate deaminase (22.214.171.124), is absent in eukaryotes and partially distributed in Bacteria and Archaea. The assignment of the second and third pathways toward deoxyuridine transformation was within the LCA, based on the complex evolutionary history of their enzymes.
Finally, the transformation of deoxy-UMP to deoxy-TMP can be carried out through two folate-dependent enzymes, thymidylate synthase ThyX (126.96.36.199) and thymidylate synthase ThyA (188.8.131.52)  (Figure 5 in gold boxes). These enzymes are not homologous, suggesting an independent evolutionary origin. Both enzymes promote methylation by using the 5,10-methylenetetrahydrofolate (CH2-H4 folate) as a carbon donor. ThyA also uses CH2-H4 folate to produce dihydrofolate (H2-folate). In contrast, ThyX uses flavin adenine dinucleotide (FAD) and NAD (P) H as cofactors to form reduced tetrahydrofolate (H4-folate) . ThyX is partially distributed in Bacteria and Archaea but is absent in eukaryotes; its counterpart, ThyA, is partially distributed in Bacteria, sparsely distributed in Archaea and widely distributed in Eukarya. In this context, it has been suggested that the catalytic differences between ThyA and ThyX influenced the evolution of bacterial genomes [61, 62]. In this regard, ThyA is 10 times more catalytically efficient than ThyX, and previous studies in more than 400 prokaryotic genomes have revealed that the catalytic capacity associated with ThyA correlates with its presence in organisms with large genome sizes [61, 62]. Our results related to its taxonomic distribution are consistent with those of earlier studies, which showed a pattern of anticorrelation for the presence/absence of these two enzyme families through bacterial clades, which has also been described as a mutual event of replacements between these two families , indicating that bacterial metabolism has modulated the size and composition of the bacterial genomes. Additionally, our data show the complete absence of ThyX and absolute presence of ThyA in eukaryotic genomes, suggesting that the influence of ThyA was a determinant in the eukaryotic genome size and in those organisms’ evolutionary potentials. In this context, it is difficult to determine whether some of these families were present in the LCA, based on their observed complex taxonomical distribution.
Pyrimidine salvage routes
In association with the salvage pathway for uracil ribonucleotide are two main routes that start from uracil. The first pathway comprises uridine kinase (184.108.40.206) and uridine phosphorylase (220.127.116.11). Uridine phosphorylase was found sparsely distributed in the three cellular domains and uses preformed uracil as a substrate. Uridine kinase is widely distributed in Bacteria and Eukarya and sparsely distributed in Archaea. Therefore, based on the taxonomical distribution patterns, it is difficult to determine whether all these enzymes were present in the LCA, due to their low distribution in Archaea genomes, which may also be a reflection of horizontal transfer events.
The enzymatic family of uridine phosphorylases uses uridine as a substrate, which can be provided by another pathway branch associated with cytidine deaminase (18.104.22.168) and pyrimidine-nucleoside phosphorylase (22.214.171.124). This branch uses cytosine as the starting substrate. Cytidine deaminase (126.96.36.199) maintains a similar taxonomical distribution as uridine kinase, whereas the pyrimidine-nucleoside phosphorylase (188.8.131.52) is only found partially distributed in Bacteria, suggesting that this pathway was not present in the LCA. In a third salvage pathway, a single enzymatic step, using uracil phosphoribosyltransferase (184.108.40.206), is required. This enzyme is widely distributed in the three cellular domains, suggesting its presence in the LCA (Figure 6B).
Regarding the salvage routes of ribonucleotides of cytosine, the main pathway occurs in two steps, involving uridine-cytidine kinase (220.127.116.11) and pyrimidine-nucleoside phosphorylase (18.104.22.168), and starting from cytosine as the substrate. Based on their taxonomical distribution, the first enzymatic family, with uridine-cytidine kinase, as previously described does not exhibit a clear distribution pattern to be considered in the LCA, while the second step is only sparsely distributed in Bacteria. In addition, there is another cytosine deoxyribonucleotide salvage pathway that starts from deoxycytidine, in which the deoxyadenosine/deoxycytidine kinase (22.214.171.124) participates (Figure 6B), and in turn it is only sparsely distributed in Eukarya.
In the salvage pathway of deoxythymine, two enzymatic steps are involved, one with thymidine kinase (126.96.36.199) and the other with thymidine phosphorylase (188.8.131.52), and they start with thymine as the substrate (Figure 6A). The enzyme associated with the first step is partially distributed in Bacteria and more sparsely distributed in Eukarya and Archaea, while the enzyme associated with the second step is sparsely distributed in the three cellular domains. It is difficult to determine whether this pathway was present in the LCA or if its distribution is a reflection of horizontal gene transfer or massive gene loss, due to the poor distribution of the enzymes belonging to this pathway in the three cellular domains.
Integrating the evolutionary analysis of the de novo and salvage pathways of pyrimidine metabolism
Based on the taxonomical distributions of enzymes associated with the de novo and salvage pathways for pyrimidine, it was possible to identify, in addition to the de novo pathway, a widespread salvage pathway (via 184.108.40.206) for the synthesis of ribonucleotide of uracil (Figure 6B), suggesting that the LCA had two pathways to carry out its synthesis. On the other hand, for the synthesis of uracil deoxyribonucleotides, two main routes, each with one step, have been identified; the first one starts from dCTP, which is converted by the enzyme 220.127.116.11, and the second route starts from dCMP, which is converted by the enzyme linked to 18.104.22.168. The wide taxonomic distribution in the three cellular domains of the 22.214.171.124 enzyme suggests that the second pathway was present in the LCA (Figure 7).
The taxonomical distribution of enzymatic families associated with DHDOS suggests that their divergence correlated with the transition from gram-negative to gram-positive bacteria and was affected by similar evolutionary pressures. The evolutionary pressures could be associated with a variety of environmental changes, such as increased atmospheric oxygen levels, temperature or changes from water to soil habitats. These changes could lead to modifications in the plasma membrane’s chemical properties, resulting in the divergence of family 2 enzymes (126.96.36.199) linked to the cell membrane, which uses quinones from the respiratory chain as electron acceptors, to the family 1 enzymes, which incorporate soluble electron acceptors. This transition from gram-negative to gram-positive was also suggested by Cavalier-Smith [59, 60]; altogether, with our data regarding the taxonomic distribution showing that family 1 is ubiquitous in gram-positive bacteria and archaea, we therefore suggest close evolutionary relationships between archaea and gram-positive bacteria.
Although gram-negative bacteria and eukarya contain the subfamily 2 of DHODS, it is not possible to deduce an evolutionary relationship, as we previously identified between gram-positive bacteria and archaea, because members of this subfamily in eukarya are linked to the mitochondrial membrane. In this regard, the mitochondrial acquisition in Eukarya has been described as a probable lateral gene transfer event, as described in the endosymbiont theory . Therefore, our results support the notion of phagocytosis of a gram-negative bacterium by a protoeukaryotic cell, with posterior specialization to mitochondria.
Additionally, de novo synthesis of cytosine ribonucleotides and deoxyribonucleotides was probably associated with the LCA. In de novo biosynthesis of thymine, the folate-dependent enzymes ThyA and ThyX can catalyze the transformation of this metabolite independently, suggesting that this pathway appeared posterior to the emergence of folate biosynthesis. For enzymes present in the thymine salvage pathway, it was not possible to determine their presence in the LCA.
The analysis presented here is based on multiple complete genomes belonging to organisms from the three cellular domains, along with current biochemical knowledge. These analyses allowed us to identify issues related to the origin and evolution of nucleotide metabolism. One of our main findings is that we were able to assess the ancestry of some segments of the purine salvage and the de novo pathways, which could be complementary and closely related to the LCA of the three cellular domains. Additionally, it was found that a large part of the de novo purine branch is widely distributed in Bacteria, Archaea and Eukarya, primarily towards de novo biosynthesis of IMP (a key precursor of purines). This branch may have been associated, in the early stages of cell evolution, with the metabolism of thiamine (vitamin B1) and posteriorly was complemented by the addition of two new enzymatic steps to complete the IMP biosynthesis pathway by means of the folate-dependent PurH enzyme, giving rise to the modern de novo synthesis of purines.
The ancestry and divergence of enzymes associated with these routes provide clues to the environmental changes in the early stages of the emergence of life. Such is the case with the divergence of the enzyme phosphoribosylaminoimidazole carboxylase (188.8.131.52) from N5-carboxyaminoimidazole ribonucleotide mutase (184.108.40.206); the divergence of these two enzymes supports the hypothesis of the origin of life in primitive seas with high levels of HCO3-. Once the atmosphere was provided with oxygen and the first eukaryotic organisms with mitochondria emerged, the enzyme 220.127.116.11 diverged to 18.104.22.168 through acquisition of a CO2 binding site. This divergence process, among other metabolic changes, may have facilitated the emergence of the first eukaryotic multicellular organisms.
In the case of pyrimidines, it was possible for us to infer that the LCA synthesized uracil ribonucleotides, by both the de novo and salvage pathways, suggesting that this ribonucleotide could have been involved in a great number of enzymatic functions and/or regulation as a remnant of the RNA world. Additionally, it was possible to associate the synthesis of cytosine and uracil deoxyribonucleotides in the LCA, and once folate biosynthesis was possible, the thymine deoxyribonucleotides emerged due to their enzymatic dependence on the folate precursor. Before the emergence of folate biosynthesis, some variants of these three bases, including cytosine or cytosine-methylated derivatives, based on its similarity to thymine, could have played a role similar to actual thymine in DNA. This inference is supported by the fact that several analogs of cytosine and uracil partially integrate into DNA as replacements for thymine [65–67]. Thus, thymine does not have an active role in the nucleotide transcript and it is more likely to be replaced. This argument is consistent with the evolution of bacterial strains , which shows that strains can be generated with the ability to incorporate derivatives of uracil (chloro-uracil) by replacing up to 98% of thymine and maintain cell viability.
In relation to purine synthesis, our results revealed that the biosynthesis of adenine could have been carried out in the LCA by the adenylate kinase. In the case of guanine biosynthesis, we found a complex evolutionary history; for instance, its synthesis has been detected in archaeal organisms, although its gene sequence has not been identified . The fact that we did not find the guanylate kinase bacterial type in Archaea suggests three possible scenarios: a) a high sequence divergence between these proteins, b) a masking function for some other enzyme, and/or c) a gene loss or gene nonorthologous displacement. As discussed above, it is likely that the guanylate kinase function has been masked by the adenylate kinase AK6 subfamily type (AK6). This enzyme is homologous to the canonical guanylate kinase of Eukarya and Bacteria and is probably closer to the possible ancestral enzyme present in the LCA, with promiscuity in the synthesis of NMPs and dNMPs to their respective NDPs and dNDPs. In addition, we found the presence of this multispecific adenylate kinase AK6 subfamily type in Archaea, suggesting that this could have been acquired from the LCA with no major posterior changes. In contrast, the guanylate kinase of Bacteria and Eukarya may have diverged to a greater extent from the adenylate kinase family subsequent to archaeal divergence. The structural similarity of inosine and guanine, together with the remnant affinity of inosine in several proteins such as RNA polymerase, in addition to the taxonomical distribution data that showed a possible LCA route for the synthesis of ITP (Figure 3, magenta circle) by means of a broad-specificity ancestral enzyme similar to the adenylate kinase Ak6 subfamily type, suggest that this base could have played an important role in early cell evolution, as has been previously proposed . The ancestral enzyme might have provided both IDP and GDP for subsequent processing to NTPs, with both playing similar informational process roles in the genetic code based on their structural similarities.
The subsequent specialization of guanylate kinase may have facilitated greater availability of guanine nucleotide to replace inosine, thus avoiding ambiguous rules in DNA replication that would have been achieved with this consolidation stage of the genetic code.
Interestingly, the divergence of the DHOD families in the LCA suggests transitions associated with changes in the cell wall and cell membrane, supporting an order of divergence from cell walls of gram-negative-like organisms and membranes similar to Eukarya-Bacteria, towards gram-positive cell wall, and/or membranes similar to Archaea. Since the plasmatic membrane is considered a matter of vertical inheritance, we suggest that the divergence of the family of DHODs can be associated with the divergence of Archaea and gram-positive bacteria.
In order to select the enzymatic numbers and their corresponding enzymes belonging to nucleotide metabolism, the KEGG  and MetaCyc  databases were exhaustively explored. In total, 120 enzymatic numbers and their corresponding enzymes were collected. (Additional file 5: Table S3 and Additional file 6: Table S4). In a second step, RPS-Blast profiles were used to search for the occurrence of members of enzyme families in complete genomes. These profiles were extracted from PRIAM, the specialized database for detection of enzymatic sequences . This database encompasses characteristics representative of alignments manually annotated, including members of a particular enzyme family, according to the Enzyme DB . In addition, enzymatic steps can be associated with more than one profile, such as protein complexes, or families of nonenzymatic homologous proteins (analogous enzymes) or steps performed by members of the same families (paralogs). Altogether, profiles were curated manually based on functional annotated domains according to Enzyme DB .
For the analysis of enzymatic functions, the best hit for each sequence, with an E-value of ≤10-10 and coverage of ≥55% in relation to the profile, was considered. Similar criteria have been previously described for enzymatic annotation of complete genomes .
In total, 151 sequence profiles associated with 120 enzymatic reactions related to nucleotide metabolism were evaluated in 2,044 complete genomes. The evaluation was based on profile comparisons, using RPS-Blast against the 2,044 complete genomes. Organisms classified as obligate parasites or those organisms with a reduced genome, i.e., less than 1,000 genes, were not considered in this study, with the aim of excluding a possible bias associated with massive gene loss, as previously described [1, 8], leaving us with a total of 1,606 of the 2,044 complete genomes. In order to exclude redundancy in the genomes analyzed, we clustered the organisms based on their taxonomical classifications, in order to obtain a normalized measure of the taxonomical distribution of enzymes according to the following steps. In the first step, we obtained the average presence/absence of enzymes per genus, which in turn was used for the second step, to calculate the average presence/absence of enzymes per Clade. Clades corresponded to the taxonomical categories from the Joint Genome Institute’s Integrated Microbial Genomics. Finally, we consider enzymes as widely distributed as those present in more than 50 percent of the clades of the same cellular domain. Additional file 1: Figure S1 and Additional file 2: Figure S2.
Identification of structural domains and families
Structural domains of proteins and families were identified using the models deposited in the Superfamily database . In brief, this database contains Hidden Markov models (HMMs) for each superfamily, which are classified according to their structural domains, which in turn are based on the classifications of SCOP ; all members of the same superfamily have a common evolutionary origin. Consensus sequences were derived from the alignments (from which profiles were constructed) and used to identify their corresponding superfamilies, using default parameters.
DA-M retrieved all enzymatic repertoires, performed the all sequence analysis and interpreted the results. LS interpreted the results; and EP-R interpreted the results. DA-M and EP-R drafted and wrote the manuscript. All authors read and approved the final manuscript.
Hernandez-Montes G, Diaz-Mejia JJ, Perez-Rueda E, Segovia L: The hidden universal distribution of amino acid biosynthetic networks: a genomic perspective on their origins and evolution. Genome Biol. 2008, 9 (6): R95-10.1186/gb-2008-9-6-r95.
Wixon J, Kell D: The Kyoto encyclopedia of genes and genomes–KEGG. Yeast. 2000, 17 (1): 48-55.
Barabasi AL, Oltvai ZN: Network biology: understanding the cell’s functional organization. Nat Rev Genet. 2004, 5 (2): 101-113. 10.1038/nrg1272.
Diaz-Mejia JJ, Perez-Rueda E, Segovia L: A network perspective on the evolution of metabolism by gene duplication. Genome Biol. 2007, 8 (2): R26-10.1186/gb-2007-8-2-r26.
Horowitz NH: On the evolution of biochemical syntheses. Proc Natl Acad Sci U S A. 1945, 31 (6): 153-10.1073/pnas.31.6.153.
Jensen RA: Enzyme recruitment in evolution of new function. Annu Rev Microbiol. 1976, 30: 409-425. 10.1146/annurev.mi.30.100176.002205.
Lopez de la Osa J, Bateman DA, Ho S, Gonzalez C, Chakrabartty A, Laurents DV: Getting specificity from simplicity in putative proteins from the prebiotic earth. Proc Natl Acad Sci U S A. 2007, 104 (38): 14941-14946. 10.1073/pnas.0706876104.
Delaye L, Becerra A, Lazcano A: The last common ancestor: what’s in a name?. Orig Life Evol Biosph. 2005, 35 (6): 537-554. 10.1007/s11084-005-5760-3.
Brown AM, Hoopes SL, White RH, Sarisky CA: Purine biosynthesis in archaea: variations on a theme. Biol Direct. 2011, 6: 63-10.1186/1745-6150-6-63.
Becerra A, Lazcano A: The role of gene duplication in the evolution of purine nucleotide salvage pathways. Orig Life Evol Biosph. 1998, 28 (4–6): 539-553.
Caetano-Anolles G, Yafremava LS, Gee H, Caetano-Anolles D, Kim HS, Mittenthal JE: The origin and evolution of modern metabolism. Int J Biochem Cell Biol. 2009, 41 (2): 285-297. 10.1016/j.biocel.2008.08.022.
Oro J, Yamrom T, Cortes S: The Prebiotic Synthesis of Imidazoles and Their Catalytic Role in the Synthesis of Oligonucleotides. Fed Proc. 1983, 42 (7): 1926-1926.
Basile B, Lazcano A, Oro J: Prebiotic syntheses of purines and pyrimidines. Adv Space Res. 1984, 4 (12): 125-131. 10.1016/0273-1177(84)90554-4.
Meléndez-Hevia E, Waddell TG, Montero F: Optimization of metabolism: the evolution of metabolic pathways toward simplicity through the game of the pentose phosphate cycle. J Theor Biol. 1994, 166 (2): 201-220. 10.1006/jtbi.1994.1018.
Morowitz H, Peterson E, Chang S: The Synthesis of Glutamic-Acid in the Absence of Enzymes - Implications for Biogenesis. Orig Life Evol Biosph. 1995, 25 (4): 395-399. 10.1007/BF01581777.
Waddell TG, Eilders LL, Patel BP, Sims M: Prebiotic methylation and the evolution of methyl transfer reactions in living cells. Orig Life Evol Biosph. 2000, 30 (6): 539-548. 10.1023/A:1026523222285.
Peregrín-Alvarez JM, Sanford C, Parkinson J: The conservation and evolutionary modularity of metabolism. Genome Biol. 2009, 10 (6): R63-10.1186/gb-2009-10-6-r63.
Fani R, Fondi M: Origin and evolution of metabolic pathways. Phys Life Rev. 2009, 6 (1): 23-52. 10.1016/j.plrev.2008.12.003.
Kyrpides N, Overbeek R, Ouzounis C: Universal protein families and the functional content of the last universal common ancestor. J Mol Evol. 1999, 49 (4): 413-423. 10.1007/PL00006564.
Koonin EV: Comparative genomics, minimal gene-sets and the last universal common ancestor. Nat Rev Microbiol. 2003, 1 (2): 127-136. 10.1038/nrmicro751.
Zrenner R, Stitt M, Sonnewald U, Boldt R: Pyrimidine and purine biosynthesis and degradation in plants. Annu Rev Plant Biol. 2006, 57: 805-836. 10.1146/annurev.arplant.57.032905.105421.
Vantrump JE, Miller SL: Carbon-Monoxide on Primitive Earth. Earth Planet Sc Lett. 1973, 20 (1): 145-150. 10.1016/0012-821X(73)90152-0.
Smirnov A, Schoonen MA: Evaluating experimental artifacts in hydrothermal prebiotic synthesis experiments. Orig Life Evol Biosph. 2003, 33 (2): 117-127. 10.1023/A:1024621500930.
Martin W, Russell MJ: On the origin of biochemistry at an alkaline hydrothermal vent. Philos Trans R Soc Lond B Biol Sci. 2007, 362 (1486): 1887-1925. 10.1098/rstb.2006.1881.
Wallace KB, Starkov AA: Mitochondrial targets of drug toxicity. Annu Rev Pharmacol Toxicol. 2000, 40: 353-388. 10.1146/annurev.pharmtox.40.1.353.
Haritos VS, Dojchinov G: Cytochrome c oxidase inhibition in the rice weevil Sitophilus oryzae (L.) by formate, the toxic metabolite of volatile alkyl formates. Comp Biochem Physiol C Toxicol Pharmacol. 2003, 136 (2): 135-143. 10.1016/S1532-0456(03)00173-X.
Kappock TJ, Ealick SE, Stubbe J: Modular evolution of the purine biosynthetic pathway. Curr Opin Chem Biol. 2000, 4 (5): 567-572. 10.1016/S1367-5931(00)00133-2.
Zhang Y, Morar M, Ealick SE: Structural biology of the purine biosynthetic pathway. Cell Mol Life Sci. 2008, 65 (23): 3699-3724. 10.1007/s00018-008-8295-8.
Tribunskikh IA, Alenin VV, Selivanov SI, Shavva AG, Inge-Vechtomov SG: Divergence of de novo biosynthesis of inosine-5’-triphosphate. Dokl Biochem Biophys. 2005, 400: 65-68. 10.1007/s10628-005-0034-6.
Dismukes G, Klimov V, Baranov S, Kozlov YN, DasGupta J, Tyryshkin A: The origin of atmospheric oxygen on Earth: the innovation of oxygenic photosynthesis. Proc Natl Acad Sci. 2001, 98 (5): 2170-2175. 10.1073/pnas.061514798.
Grotzinger JP, Kasting JF: New constraints on Precambrian ocean composition. J Geol. 1993, 101 (2): 235-243. 10.1086/648218.
Meyer E, Kappock TJ, Osuji C, Stubbe J: Evidence for the direct transfer of the carboxylate of N5-carboxyaminoimidazole ribonucleotide (N5-CAIR) to generate 4-carboxy-5-aminoimidazole ribonucleotide catalyzed by Escherichia coli PurE, an N5-CAIR mutase. Biochemistry. 1999, 38 (10): 3012-3018. 10.1021/bi9827159.
Cusack NJ, Shaw G, Litchfield GJ: Purines, pyrimidines, and imidazoles. XXXVI. Carboxylation of some 5-aminoimidazoles and related compounds, including nucleosides and nucleotides, with potassium hydrogen carbonate in aqueous solution. J Chem Soc Perkin 1. 1971, 8: 1501-1507.
Cordat E, Casey JR: Bicarbonate transport in cell physiology and disease. Biochem J. 2009, 417: 423-439. 10.1042/BJ20081634.
Kang YN, Tran A, White RH, Ealick SE: A novel function for the N-terminal nucleophile hydrolase fold demonstrated by the structure of an archaeal inosine monophosphate cyclohydrolase. Biochemistry. 2007, 46 (17): 5050-5062. 10.1021/bi061637j.
Grochowski LL, Censky K, Xu HM, White RH: A new class of adenylate kinase in methanogens is related to uridylate kinase. Arch Microbiol. 2012, 194 (2): 141-145. 10.1007/s00203-011-0759-9.
Yan H, Tsai MD: Nucleoside monophosphate kinases: structure, mechanism, and substrate specificity. Adv Enzymol Relat Areas Mol Biol. 1999, 73: 103-134.
Stolworthy TS, Black ME: The mouse guanylate kinase double mutant E72Q/D103N is a functional adenylate kinase. Protein Eng. 2001, 14 (11): 903-909. 10.1093/protein/14.11.903.
Feng X, Yang RN, Zheng XF, Zhang FY: Identification of a novel nuclear-localized adenylate kinase 6 from Arabidopsis thaliana as an essential stem growth factor. Plant Physiol Biochem. 2012, 61: 180-186.
Ren H, Wang L, Bennett M, Liang Y, Zheng X, Lu F, Li L, Nan J, Luo M, Eriksson S, Zhang C, Su XD: The crystal structure of human adenylate kinase 6: An adenylate kinase localized to the cell nucleus. Proc Natl Acad Sci U S A. 2005, 102 (2): 303-308. 10.1073/pnas.0407459102.
Keegan LP, Leroy A, Sproul D, O’Connell MA: Adenosine deaminases acting on RNA (ADARs): RNA-editing enzymes. Genome Biol. 2004, 5 (2): 209-10.1186/gb-2004-5-2-209.
Hitchings GH: The purine metablism of protozoa. Adv Enzyme Regul. 1982, 20: 375-386.
Zubay G: To what extent do biochemical pathways mimic prebiotic pathways. Chemtracts Biochem Mol Biol. 1993, 4: 317-323.
Frank RA, Leeper FJ, Luisi BF: Structure, mechanism and catalytic duality of thiamine-dependent enzymes. Cell Mol Life Sci. 2007, 64 (7–8): 892-905.
Lazcano A, Miller SL: On the origin of metabolic pathways. J Mol Evol. 1999, 49 (4): 424-431. 10.1007/PL00006565.
Groziak MP, Bhat B, Leonard NJ: Nonenzymatic synthesis of 5-aminoimidazole ribonucleoside and recognition of its facile rearrangement. Proc Natl Acad Sci U S A. 1988, 85 (19): 7174-7176. 10.1073/pnas.85.19.7174.
Winkler W, Nahvi A, Breaker RR: Thiamine derivatives bind messenger RNAs directly to regulate bacterial gene expression. Nature. 2002, 419 (6910): 952-956. 10.1038/nature01145.
Sakumi K, Abolhassani N, Behmanesh M, Iyama T, Tsuchimoto D, Nakabeppu Y: ITPA protein, an enzyme that eliminates deaminated purine nucleoside triphosphates in cells. Mutat Res. 2010, 703 (1): 43-50. 10.1016/j.mrgentox.2010.06.009.
Glasemacher J, Bock AK, Schmid R, Schonheit P: Purification and properties of acetyl-CoA synthetase (ADP-forming), an archaeal enzyme of acetate formation and ATP synthesis, from the hyperthermophile Pyrococcus furiosus. Eur J Biochem. 1997, 244 (2): 561-567. 10.1111/j.1432-1033.1997.00561.x.
Krepl M, Otyepka M, Banas P, Sponer J: Effect of Guanine to Inosine Substitution on Stability of Canonical DNA and RNA Duplexes: Molecular Dynamics Thermodynamics Integration Study. J Phys Chem B. 2013, 117 (6): 1872-1879. 10.1021/jp311180u.
Rowland P, Norager S, Jensen KF, Larsen S: Structure of dihydroorotate dehydrogenase B: electron transfer between two flavin groups bridged by an iron-sulphur cluster. Structure. 2000, 8 (12): 1227-1238. 10.1016/S0969-2126(00)00530-X.
Takashima E, Inaoka DK, Osanai A, Nara T, Odaka M, Aoki T, Inaka K, Harada S, Kita K: Characterization of the dihydroorotate dehydrogenase as a soluble fumarate reductase in Trypanosoma cruzi. Mol Biochem Parasitol. 2002, 122 (2): 189-200. 10.1016/S0166-6851(02)00100-7.
Norager S, Jensen KF, Bjornberg O, Larsen S: E-coli dihydroorotate dehydrogenase reveals structural and functional distinctions between different classes of dihydroorotate dehydrogenases. Structure. 2002, 10 (9): 1211-1223. 10.1016/S0969-2126(02)00831-6.
Bernstein MP, Sandford SA, Allamandola LJ, Gillette JS, Clemett SJ, Zare RN: UV irradiation of polycyclic aromatic hydrocarbons in ices: production of alcohols, quinones, and ethers. Science. 1999, 283 (5405): 1135-1138. 10.1126/science.283.5405.1135.
Deamer D, Dworkin JP, Sandford SA, Bernstein MP, Allamandola LJ: The first cell membranes. Astrobiology. 2002, 2 (4): 371-381. 10.1089/153110702762470482.
Wynn-Williams DDEH: Environmental UV radiation: biological strategies for protection and avoidance. Astrobiology The quest for the conditions of life. Edited by: Gerda Horneck CB-K. 2002, Berlin: Springer, 245-260.
Deamer DW: Boundary structures are formed by organic-components of the Murchison Sarbonaceous Chondrite. Nature. 1985, 317 (6040): 792-794. 10.1038/317792a0.
Cavalier-Smith T: The neomuran origin of archaebacteria, the negibacterial root of the universal tree and bacterial megaclassification. Int J Syst Evol Microbiol. 2002, 52 (Pt 1): 7-76.
Cavalier-Smith T: Cell evolution and Earth history: stasis and revolution. Philos Trans R Soc Lond B Biol Sci. 2006, 361 (1470): 969-1006. 10.1098/rstb.2006.1842.
Lombard J, Lopez-Garcia P, Moreira D: The early evolution of lipid membranes and the three domains of life. Nat Rev Microbiol. 2012, 10 (7): 507-515.
Myllykallio HSS, Grosjean H: Folate-Dependent Thymidylate-Forming Enzymes: Parallels between DNA and RNA Metabolic Enzymes and Evolutionary Implications. DNA and RNA Modification Enzymes: Structure, Mechanism, Function and Evolution. Edited by: Henri G. 2009, Austin (TX): Landes Bioscience, 275-288.
Escartin F, Skouloubris S, Liebl U, Myllykallio H: Flavin-dependent thymidylate synthase X limits chromosomal DNA replication. Proc Natl Acad Sci U S A. 2008, 105 (29): 9948-9952. 10.1073/pnas.0801356105.
Myllykallio H, Lipowski G, Leduc D, Filee J, Forterre P, Liebl U: An alternative flavin-dependent mechanism for thymidylate synthesis. Science. 2002, 297 (5578): 105-107. 10.1126/science.1072113.
Chan CX, Bhattacharya D, Reyes-Prieto A: Endosymbiotic and horizontal gene transfer in microbial eukaryotes: Impacts on cell evolution and the tree of life. Mob Genet Elements. 2012, 2 (2): 101-105. 10.4161/mge.20110.
Nygaard P: Utilisation of preformed purine bases and nucleosides. Metabolism of Nucleotides, Nucleosides and Nucleobases in Microorganisms. Edited by: Munch-Peterson A. 1983, London: Academic Press, 27-93.
Bick MD, Davidson RL: Total substitution of bromodeoxyuridine for thymidine in the DNA of a bromodeoxyuridine-dependent cell line. Proc Natl Acad Sci U S A. 1974, 71 (5): 2082-2086. 10.1073/pnas.71.5.2082.
Warren RA: Modified bases in bacteriophage DNAs. Annu Rev Microbiol. 1980, 34: 137-158. 10.1146/annurev.mi.34.100180.001033.
Marliere P, Patrouix J, Doring V, Herdewijn P, Tricot S, Cruveiller S, Bouzon M, Mutzel R: Chemical evolution of a bacterium’s genome. Angew Chem Int Ed Engl. 2011, 50 (31): 7109-7114. 10.1002/anie.201100535.
Caspi R, Foerster H, Fulcher CA, Hopkinson R, Ingraham J, Kaipa P, Krummenacker M, Paley S, Pick J, Rhee SY, Tissier C, Zhang P, Karp PD: MetaCyc: a multiorganism database of metabolic pathways and enzymes. Nucleic Acids Res. 2006, 34 (Database issue): D511-D516.
Claudel-Renard C, Chevalet C, Faraut T, Kahn D: Enzyme-specific profiles for genome annotation: PRIAM. Nucleic Acids Res. 2003, 31 (22): 6633-6639. 10.1093/nar/gkg847.
Bairoch A: The ENZYME database in 2000. Nucleic Acids Res. 2000, 28 (1): 304-305. 10.1093/nar/28.1.304.
Mavromatis K, Ivanova NN, Chen I-MA, Szeto E, Markowitz VM, Kyrpides NC: The DOE-JGI Standard operating procedure for the annotations of microbial genomes. Stand Genomic Sci. 2009, 1 (1): 63-10.4056/sigs.632.
Wilson D, Madera M, Vogel C, Chothia C, Gough J: The SUPERFAMILY database in 2007: families and functions. Nucleic Acids Res. 2007, 35 (Database issue): D308-D313.
Lo Conte L, Ailey B, Hubbard TJ, Brenner SE, Murzin AG, Chothia C: SCOP: a structural classification of proteins database. Nucleic Acids Res. 2000, 28 (1): 257-259. 10.1093/nar/28.1.257.
We thank Georgina Hernandez-Montes and Jessica Brambila-Tapia for their critical reading of the manuscript. DA-M acknowledges the Programa de Maestría y Doctorado en Ciencias Bioquímicas at Instituto de Biotecnologia-UNAM and his support by a PhD fellowship (35083) from CONACYT and (IN-204714) DGAPA. EP-R was supported by a grant (IN-204714) from DGAPA and (155116) from CONACYT. We also want to thank Juan Manuel Hurtado Ramirez by his technical support.
The authors declare that they have no competing interests.
Electronic supplementary material
Additional file 1: Figure S1: Average taxonomic distribution of purine metabolism enzymes distributed across the three domains of life. The taxonomic distribution for enzymes catalyzing the purine metabolism (vertical labels) was computed by searching for their ortholog distribution across diverse taxonomic groups (Archaea, Bacteria and Eukarya) (horizontal labels). Enzymes are sorted in terms of their E.C. number. Some identical E.C. numbers and superfamily assignation are complementary according to PRIAM db. The “S” notation defines protein subunits; “AN” defines ANalogous enzymes, those defined as enzymes with the same E.C. number and the different superfamily classification, suggesting different evolutionary origin. The colors are in rainbow scale, where dark blue are those enzymes with low average (0) whereas red denotes high average values (1). (TIFF 311 KB)
Additional file 2: Figure S2: Average taxonomic distribution of pyrimidine metabolism enzymes. Labels and colors are as in Additional file 1: Figure S1. In asterisks are indicated those enzymes common to purine and pyrimidine metabolisms. (TIFF 433 KB)
Additional file 3: Table S1: Matrix used to construct the average taxonomic distribution. Purine metabolism. Column 1 shows the E.C. numbers; columns 2 to 76 correspond to the taxonomical divisions analyzed here. Each cell includes normalized values as described in material and methods. (XLSX 35 kb). (XLSX 35 KB)
Additional file 4: Table S2: Matrix used to construct the average taxonomic distribution. Pyrimidine metabolism. Column 1 shows the E.C. numbers; columns 2 to 76 correspond to the taxonomical divisions analyzed here. Each cell includes normalized values as described in material and methods. (XLSX 34 KB)
Additional file 5: Table S3: Enzymes analyzed. Purine metabolism (Map 00230). Column 1 denotes the E.C. Number; column 2, the accepted name; column 3, reactions associated to each E. C. number in the purine metabolism; and column 4, the identification code on KEGG database. (DOCX 32 KB)
Additional file 6: Table S4: Enzymes analyzed. Pyrimidine metabolism (Map 00240). Column 1 denotes the E.C. Number; column 2, the accepted name; column 3, reactions associated to each E. C. number in the pyrimidine metabolism; and column 4, the identification code on KEGG database. (DOCX 35 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Armenta-Medina, D., Segovia, L. & Perez-Rueda, E. Comparative genomics of nucleotide metabolism: a tour to the past of the three cellular domains of life. BMC Genomics 15, 800 (2014). https://doi.org/10.1186/1471-2164-15-800
- Nucleotide metabolism
- Comparative genomics
- Sequence profiles