PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of bmcgenoBioMed Centralsearchsubmit a manuscriptregisterthis articleBMC Genomics
 
BMC Genomics. 2007; 8: 245.
Published online Jul 23, 2007. doi:  10.1186/1471-2164-8-245
PMCID: PMC1971073

Comparative genomics of bacterial and plant folate synthesis and salvage: predictions and validations

Abstract

Background

Folate synthesis and salvage pathways are relatively well known from classical biochemistry and genetics but they have not been subjected to comparative genomic analysis. The availability of genome sequences from hundreds of diverse bacteria, and from Arabidopsis thaliana, enabled such an analysis using the SEED database and its tools. This study reports the results of the analysis and integrates them with new and existing experimental data.

Results

Based on sequence similarity and the clustering, fusion, and phylogenetic distribution of genes, several functional predictions emerged from this analysis. For bacteria, these included the existence of novel GTP cyclohydrolase I and folylpolyglutamate synthase gene families, and of a trifunctional p-aminobenzoate synthesis gene. For plants and bacteria, the predictions comprised the identities of a 'missing' folate synthesis gene (folQ) and of a folate transporter, and the absence from plants of a folate salvage enzyme. Genetic and biochemical tests bore out these predictions.

Conclusion

For bacteria, these results demonstrate that much can be learnt from comparative genomics, even for well-explored primary metabolic pathways. For plants, the findings particularly illustrate the potential for rapid functional assignment of unknown genes that have prokaryotic homologs, by analyzing which genes are associated with the latter. More generally, our data indicate how combined genomic analysis of both plants and prokaryotes can be more powerful than isolated examination of either group alone.

Background

Folates are tripartite molecules comprising pterin, p-aminobenzoate (pABA), and glutamate moieties to which one-carbon units at various oxidation levels can be attached at the N5 and N10 positions (Figure (Figure1).1). In natural folates the pterin ring is in the dihydro or tetrahydro state, and a short, γ-linked polyglutamyl tail of up to about eight residues is usually attached to the first glutamate.

Figure 1
The structure of tetrahydrofolate. In natural folates, the pterin ring exists in tetrahydro form (as shown) or in 7,8-dihydro form (as in DHF). The ring is fully oxidized in folic acid, which is not a natural folate. Folates usually have a γ-linked ...

Tetrahydrofolates serve as cofactors in one-carbon transfer reactions during the synthesis of purines, formylmethionyl-tRNA, thymidylate, pantothenate, glycine, serine, and methionine [1] (Figure (Figure2).2). Most folate-dependent enzymes strongly prefer polyglutamates to monoglutamates, but the opposite is usually true of folate transporters so that polyglutamylation is generally considered to favor folate retention within cells and subcellular compartments [2,3].

Figure 2
Major folate-dependent reactions of one-carbon metabolism. The gene names are for E. coli (except for sarcosine dehydrogenase). Note that the formation of 5-formyl-THF from 5,10-methenyl-THF occurs via a second catalytic activity of serine hydroxymethyltransferase ...

Plants, fungi, certain protists, and most bacteria make folates de novo, starting from GTP and chorismate, but higher animals lack key enzymes of the synthetic pathway and so require dietary folate [4-7]. Folates are crucial to human nutrition and health [3], and antifolate drugs are widely used in cancer chemotherapy and as antimicrobials [3,7,8]. For these reasons, folate synthesis and salvage pathways have been extensively characterized in model organisms, and the folate synthesis pathway in both bacteria and plants has been engineered in order to boost the folate content of foods [9-11].

The de novo folate synthesis pathway has the same steps in bacteria and plants, and consists of a pterin branch and a pABA branch (Figure (Figure3,3, rose and blue color, respectively). The first enzyme of the pterin branch is GTP cyclohydrolase I (GCHY-I, EC 3.5.4.16), which catalyzes a complex reaction in which the five-membered imidazole ring of GTP is opened, C8 is expelled as formate, and a six-membered dihydropyrazine ring is formed using C1 and C2 of the ribose moiety of GTP [5]. The resulting 7,8-dihydroneopterin triphosphate is then converted to the corresponding monophosphate by a specific pyrophosphatase [5,12]. Removal of the last phosphate is believed to be mediated by a non-specific phosphatase [5]. Dihydroneopterin aldolase (DHNA, EC 4.1.2.25) then releases glycolaldehyde to produce 6-hydroxymethyl-7,8-dihydropterin, which is then pyrophosphorylated by hydroxymethyldihydropterin pyrophosphokinase (HPPK, EC 2.7.6.3). DHNA also interconverts 7,8-dihydroneopterin and 7,8-dihydromonapterin, and cleaves the latter to 6-hydroxymethyl-7,8-dihydropterin. A paralog of DHNA, FolX, interconverts the triphosphates of 7,8-dihydroneopterin and 7,8-dihydromonapterin, and also catalyzes the same reactions as DHNA at very slow rates [13].

Figure 3
Folate synthesis and salvage pathways. Gene names are white-on-gray; all except folQ and PTR1 are from E. coli. The folQ gene has been identified only in Lactococcus lactis and plants, and PTR1 only in Leishmania and other trypanosomatids. Note that DHN ...

In the pABA branch of the pathway, chorismate is aminated to aminodeoxychorismate (ADC) by ADC synthase (EC 6.3.5.8) using the amide group of glutamine as amino donor [5]. ADC is then converted to pABA by ADC lyase (EC 4.1.3.38) [5].

6-Hydroxymethyl-7,8-dihydropterin pyrophosphate and pABA moieties are condensed by dihydropteroate synthase (DHPS, EC 2.5.1.15). The resulting dihydropteroate is glutamylated by dihydrofolate synthase (DHFS, EC 6.3.2.12) giving dihydrofolate (DHF), which is reduced by dihydrofolate reductase (DHFR, EC 1.5.1.3) to tetrahydrofolate (THF). Folylpolyglutamate synthase (FPGS, EC 6.3.2.17) then adds a γ-glutamyl tail. In Escherichia coli, it has been reported that there can also be α linkages in the distal part of the polyglutamyl tail [14].

Although the biosynthetic steps are the same in plants and bacteria, the plant pathway is split between three subcellular compartments, with pterin synthesis in the cytosol, pABA synthesis in chloroplasts, and the other steps in mitochondria (Figure (Figure4)4) [6]. FPGS isoforms are present in all three of these compartments, as are folates themselves [15,16]. Folates – both poly- and monoglutamates – are also found in plant vacuoles [16]. The highly compartmented nature of folate synthesis in plants implies the existence of pterin and folate transporters that are integral components of the pathway.

Figure 4
Compartmentation of the folate synthesis pathway in plants. The steps in the pterin branch of the pathway are in blue, those in the pABA branch are in green, and the others (condensation, glutamylation, and reduction) are in red. Note that the compartmentation ...

Folate-related salvage pathways are of three kinds. The first ('intact folate salvage') (Figure (Figure3,3, green color) enables utilization of supplied folic acid and DHF, and relies on a DHFR activity to reduce these oxidized folates to THF, and on an FPGS activity [7]. DHFR activity is also required to recycle the DHF produced in the reaction catalyzed by thymidylate synthase (TS, EC 2.1.1.45). The second kind of salvage ('pterin salvage') (Figure (Figure3,3, yellow color), known in Leishmania and other trypanosomatid parasites, involves the reduction of fully oxidized pterins to the dihydro and tetrahydro levels by pteridine reductase 1 (PTR1, EC 1.5.1.33) [17]. This enables oxidized pterins to be used (after reduction to dihydro forms) for folate synthesis, and (after reduction to tetrahydro forms) as cofactors for aromatic hydroxylases and other pterin-dependent enzymes. Finally, some bacteria, plants, and protists probably carry out a more radical kind of salvage, in which the pterin and pABA-glutamate fragments produced by folate breakdown are recycled for folate synthesis [18]. This type of salvage has been little studied and will not be considered further in this article.

Genes for all the enzymes of folate synthesis have been identified in model organisms such as Escherichia coli, Saccharomyces cerevisiae, and Arabidopsis thaliana [4-6]. Likewise, the intact folate salvage pathway has been well characterized in mammals, the malaria parasite Plasmodium, and Lactobacillus casei [7,19,20], and pterin salvage in Leishmania [17]. However, analysis of the distribution of known folate synthesis and salvage genes in hundreds of bacterial genomes using the SEED platform [21] reveals that much remains to be learnt about both synthesis and salvage.

The SEED is a freely available, open-source database that provides efficient ways to discover new genes or pathways, to generate predictions about gene function, and to improve annotations, based on a 'functional subsystem approach' [21]. This approach has much in common with metabolic reconstruction [22,23]. A functional subsystem may be defined as a set of functional roles (usually ten to twenty) jointly involved in a biological process. A typical subsystem is a group of enzymes, transporters, and regulatory components that participate in a metabolic pathway such as folate synthesis or salvage. Subsystem analysis examines which components are actually present in a genome and which should be present but cannot be identified, and so provides a picture of what is actually missing. This sets the stage to pursue the 'missing genes', also termed 'pathway holes' [24-26]. Homology-based searches alone are usually unable to locate missing genes that have not been previously identified in any genome ('globally missing genes') or those that are missing due to non-orthologous gene replacement ('locally missing genes') [27].

In this study, we first predicted the pathways (de novo folate synthesis, intact folate salvage, and pterin salvage) present in around four hundred sequenced bacteria and identified cases of missing genes for almost every step of the synthesis pathway. Candidates for such missing genes in bacteria and plants were then predicted using comparative genomic tools and representative candidates were tested experimentally.

Results and Discussion

Are folates essential in all bacteria?

As folate-dependent formylation of the initiator tRNA is a hallmark of bacterial translation and bacteria cannot import formylmethionyl-tRNA [28], we investigated the distribution of the fmt gene encoding methionyl-tRNA formyltransferase (EC 2.1.2.9) as a signature gene for a folate requirement. Homologs of fmt are found in all sequenced genomes except Mycoplasma hyopneumoniae and Onion yellows phytoplasma OY-M (Table (Table1).1). We confirmed the observation [29] that M. hyopneumoniae lacks all the enzymes of folate-mediated one-carbon metabolism except for glycine hydroxymethyltransferase (GlyA), which has aldolase activities that do not require folate [30]. Another widespread folate-dependent metabolic step is the conversion of dTMP to dUMP, catalyzed by thymidylate synthase (ThyA, EC 2.1.1.45). This step can also be performed by a folate- and flavin-dependent thymidylate synthase (ThyX) [31]. As first observed by Myllykallio et al. [32], most bacteria have a thyA or a thyX homolog, some have both, and the few that have neither – such as M. hyopneumoniae or Ureaplasma parvum – contain the tdk gene encoding the thymidine (dT) salvage enzyme thymidine kinase. Our genomic analysis suggests that M. hyopneumoniae strains are the only sequenced bacteria that do not require folate for initiator tRNA formylation or thymidylate synthesis. The situation in the phytoplasma that lacks the fmt gene (Table (Table1)1) is different; it contains a thyA homolog like most Mycoplasma species and therefore presumably requires intact folates.

Table 1
Examples of bacteria dependent on folate salvage

Intact folate transport and salvage

As just discussed, folate is most probably essential for all sequenced bacteria except M. hyopneumoniae. However, not all bacteria synthesize folate de novo but instead rely on an external supply [see Additional File 1, variant 001; see "Methods" for an explanation of the variant code]. To predict the absence of the de novo synthesis pathway, the HPPK (FolK) and DHPS (FolP) proteins were used as signature proteins (for reasons described below). Many bacteria lack homologs of both these genes (Table (Table1)1) and so almost certainly rely on reducing and glutamylating intact folates taken up from the environment. These are mainly host-associated bacteria such as Mycoplasma or Treponema or organisms that live in folate-rich environments such as Lactobacilli. Chloroplasts and vacuoles must likewise take up folates from the cytoplasm (Figure (Figure4),4), and there is also evidence for folate uptake by intact plant cells [6].

(i) Transport

Systems that mediate folate uptake in auxotrophs such as Lactobacillus casei and L. salivarius have been partially biochemically characterized [33,34], but the corresponding genes remain unknown. Whatever they are, they are most likely unrelated to mammalian folate carriers (i.e., the reduced folate carrier, the folate receptor, the intestinal folate carrier, and the mitochondrial folate carrier) since these lack close homologs among bacteria and plants. However, cyanobacteria, which are folate prototrophs, have a protein with significant similarity to a folate carrier from Leishmania species (FT1), and the cyanobacterial protein has a close homolog in plants (52% amino acid identity), as well as several more distant relatives in plants and in alpha-, beta-, and gamma-proteobacteria. We showed first that the cyanobacterial protein (Synechocystis slr0642) conferred the ability to transport folates and folate analogs when expressed in E. coli, and then that its plant homolog (Arabidopsis At2g32040) did the same [35]. We further showed that the Arabidopsis At2g32040 protein is located in the chloroplast envelope [35]. The weak slr0642 homolog in some alpha-proteobacteria (Silicibacter, Roseobacter) clusters with the folate-dependent enzyme sarcosine dehydrogenase, suggesting that this protein may also be a folate transporter.

Thus, despite progress in identifying folate transporters in cyanobacteria and in the chloroplast envelope, there are as yet no candidates for the folate carriers in many folate-requiring bacterial taxa, or in plant mitochondrial, vacuolar, and plasma membranes. These still-missing genes are future prospects for discovery by comparative genomics methods [36].

(ii) Reduction

As noted above, DHFR is essential in both de novo and salvage pathways. Most bacteria have a folA gene (DHFR0), but two other bacterial enzymes able to reduce DHF are now known: FolM (DHFR1) belonging to the short-chain dehydrogenase/reductase (SDR) family [37], and a flavin-dependent dihydropteroate reductase that is fused to dihydropteroate synthase (DHFR2). [38]. The trypanosomatid enzyme PTR1 can also reduce DHF and folic acid [17]. As folM occurs in E. coli and other bacteria that also have a folA gene, its normal function is most probably not folate reduction, as discussed in a later section. The annotation of DHFR0 family members is complicated by their similarity to pyrimidine dehydrogenase family members (Pfam01872), which are numerous in Actinomycetes like Streptomyces coelicolor. At this stage we named them all DHFR0 but further genetic or biochemical analysis is needed to check these assignments.

Analysis of the distribution of DHFR genes in bacterial genomes reinforced the conclusions [32] that many bacteria such as Prochlorococcus marinus lack any recognizable DHFR proteins, and that most of these organisms use ThyX and not ThyA. Even if a high capacity for DHF reduction is not needed in ThyX-dependent organisms [39], these do require some DHFR activity to complete the de novo or salvage pathways so the corresponding gene(s) have yet to be identified in these organisms [32] (see Additional File 1, variants 106, 116, 006).

(iii) Glutamylation

FolC-like proteins can have FPGS activity alone [20] or both DHFS and FPGS activities [40], which complicates annotation. Although the bifunctional type has a unique dihydropteroate binding site [41], it overlaps the rest of the substrate binding site and we could not derive a motif to distinguish mono- and bifunctional enzymes. We therefore annotated them all as bifunctional. By analogy with the Lactococcus. lactis situation, we predict that organisms reliant on the salvage pathway (see Additional File 1, variants 001 and 011) will have a monofunctional FPGS. The folC gene is missing in the Mycoplasma species that contain an fmt, a thyA and a folA gene and must therefore rely on a salvage pathway (Table (Table1).1). This absence points to three possibilities for these species: (a) they import folate polyglutamates; (b) they have a novel type of FPGS gene; or (c) they import monoglutamyl folates and polyglutamylation is not needed. We favor the last hypothesis as there is evidence for monoglutamyl folate uptake in Mycoplasma mycoides [42]. A similar situation must exist in bacteria such as Borrelia burgdorferi that lack all folate synthesis genes but contain THF-dependent enzymes such as Fmt (Table (Table11).

De novo folate biosynthesis

The majority of sequenced bacteria (250 out of 400) contain all genes of the pathway and are therefore predicted to be prototrophic for folate (see examples in Table Table22 and Additional File 1, variant 111). However, a substantial minority lack just one or a few genes of the pterin or pABA branches, and detailed analysis of these cases reveals several biologically significant points.

Table 2
Examples of bacteria capable of de novo folate synthesis and of genes that are still missing

(i) The pterin branch

The first enzyme of this branch, GCHY-I, is encoded in E. coli by the folE gene. A recent analysis of the distribution of folE genes among bacterial genomes showed the folE gene to be locally missing in one-third of them [43]. Another protein family, COG1469, was found to responsible for 7,8-dihydroneopterin triphosphate formation in these organisms. This protein was named GCHY-IB and the corresponding gene folE2 [43] (Table (Table2).2). Further analysis revealed that a few bacteria such as Wolbachia, Chlamydia, and Chlamydophila species lack both folE and folE2 homologs whereas they contain the signature genes of the pathway folKP (see Table Table22 and additional File 1, variants 701, 702), suggesting that another family of GCHY-I enzymes has yet to be identified. For instance, at least certain Chlamydia species are known to synthesize folates de novo [44], but lack folE and folE2. A candidate for the missing GCYH-I enzyme was the Chlamydia trachomatis protein CT610 and its homologs, which cluster with the folABKP folate genes in Chlamydia and Wolbachia species (Figure (Figure5A).5A). The protein is homologous to the pyrroloquinoline quinone (PQQ) biosynthesis protein PqqC that catalyzes an overall eight-electron oxidation, leading to a pyrrole and pyridine ring, but their active sites are not conserved, consistent with a different enzymatic activity [45]. The CT610 gene was cloned in pBAD24 but failed to complement the dT auxotrophy of the E. coli folE mutant. The strong linkage of CT610 homologs with folate genes certainly points to a function in folate metabolism as other de novo folate genes than folE are missing in chlamydiae such as folQ or pabAabc (Table (Table2),2), but further studies are needed to determine its functional role.

Figure 5
Clustering of predicted folate-related genes with known folate synthesis genes. Gene names are as described in the text or given below. [For full gene and genome names, see Additional File 1.] Matching colors correspond to orthologous genes. Pale grey ...

The second step of folate synthesis is the removal of pyrophosphate. Although an enzyme mediating this step had been demonstrated in E. coli [46], no gene was known from any organism. We identified a DHNTP pyrophosphatase (FolQ) candidate in L. lactis as part of the folKEPQC gene cluster (Table (Table2)2) [12]. FolQ belongs to the Nudix (Nucleoside diphosphate X) hydrolase family [47]. Biochemical and genetic tests confirmed DHNTP pyrophosphatase activity [12]. Furthermore, the closest Arabidopsis homolog of L. lactis FolQ was also shown to have this activity [12].

Since the Nudix family is large and functionally heterogeneous it is not very amenable to projection of annotations just by homology. FolQ homologs with a high homology score occur in rather few bacteria, so that the DHNTP pyrophosphatase gene is still missing in most genomes, including E. coli. Other putative phosphohydrolases unrelated to FolQ, FolQ2 members of the HDIG superfamily, are found in some folate-related gene clusters (Figure (Figure5B),5B), such as CPE1020 in Clostridium perfringens; these genes are good candidates for alternatives to FolQ but again have a limited phylogenetic distribution leaving the problem open in most bacterial species (Table (Table22).

The third specific enzyme of the pathway, DHNA, is encoded in E. coli by the folB gene. This gene and its paralog folX [13] appear to be missing in many phylogenetically diverse bacteria such as Geobacter metallireducens. Genome and functional context analysis allows the prediction that the DHNA role is played by members of the transaldolase (EC 2.2.1.2) family (e.g. DVU1658 in Desulfovibrio vulgaris). Specifically, about half the bacteria that lack DHNA have a transaldolase encoding gene that clusters with folK genes in several organisms (Figure (Figure5C).5C). This prediction awaits experimental validation as this transaldolase family is broad and only some of its members might encode a DHNA aldolase. Some genomes such as Rickettsia felis lack both FolB and transaldolase homologs while containing all the other de novo enzymes (see Table Table22 and additional File 1, variant 401), again suggesting that another family of FolB enzymes has yet to be identified unless the pathway is on its way to elimination in these organisms specifically.

HPPK (FolK) and DHPS (FolP) are distinctive proteins found in all organisms that make folate de novo and so, as noted above, these were used as pathway signature genes. A few sporadic organisms apparently lack one of the two genes, but further analysis shows that this is usually because of a gene-calling problem (a homolog can be found using the tblastn algorithm) or because the corresponding genome is still incomplete. Some organisms, however, have two folP genes or two folK genes (Table (Table2).2). Are these functionally redundant or catalyzing different reactions? In most cases one paralog is clustered with folate genes and the other clusters with genes involved in different pathways (see Table Table22 and additional File 1). For instance, in the high-GC gram-positive group the second folP (folP2) clusters with cell wall synthesis genes. In Mycobacterium leprae the folP2 gene does not complement an E. coli folP mutant whereas the copy that clusters with the folate genes (folP1) does, suggesting that folP2 is involved in another pathway [49].

FolK is duplicated in many organisms. In most cases such as Shewanella denitrificans (Table (Table2),2), one copy is in a folate operon and the other in a pantothenate operon but there are several cases where both genes are close to other folate biosynthesis genes (see also Additional File 1). Only experimental testing will show whether both copies are active. It is of note that an internal duplication of FolK and fusion with FolB is found in Bifidobacterium longum.

The sequenced chlamydiae all lack homologs of folC (DHFS/FPGS) but have folPK homologs (see Table Table22 and additional File 1), making folC a locally missing gene in this group. Inspection revealed that a member of gene family COG1478 is clustered in chlamydiae with folate biosynthesis genes (Figure (Figure5A,5A, folC2). This COG1478 family contains the F420:γ-glutamyl ligase CofE of Archaea and Mycobacteria [50]. CofE catalyzes the GTP-dependent successive addition of two γ-linked L-glutamates to the L-lactyl phosphodiester of 7,8-didemethyl-8-hydroxy-5-deazariboflavin (F420), a reaction analogous to that mediated by FolC. Chlamydiae almost certainly do not make F420 since they lack all the other known cof genes [50]. We accordingly predicted that the CofE homolog in chlamydiae has FolC activity. A cofE homolog (CT611) was shown to complement the methionine and glycine requirements of the E. coli folC mutant SF4 [40] indicating that CT611 can indeed functionally replace FolC (Figure (Figure6).6). The E. coli folC gene from the ASKA collection [51] was used as a positive control.

Figure 6
Complementation of folC function by Chlamydia trachomatis CT611. Complementation of E. coli folC mutant SF4 by a pBAD24 plasmid harboring CT611 from Chlamydia trachomatis on MS minimal medium with or without glycine plus methionine. E. coli folC was included ...

(ii) The pABA branch

We adopted the nomenclature of Xie et al. [52] for the pABA branch genes. These genes are hard to annotate for several reasons. In the first place, they can be fused in various combinations. A fusion between the subunits of ADC synthase (PabAa and PabAb) is a common arrangement, as is fusion between PabAa and ADC lyase (PabAc). In one genome, Corynebacterium diphtheriae, our analysis indicated a triple fusion. The functions of this PabAa-PabAb-PabAc fusion gene (DIP1790) were tested experimentally. The gene was cloned into an expression vector and introduced into an E. coli pabAa pabAb mutant (strain BN1163), which cannot grow on minimal medium unless it expresses a recombinant enzyme with ADC synthase activity. A bifunctional PabAa-PabAb ADC synthase protein from Arabidopsis served as a positive control. Like the positive control, expression of the DIP1790 protein restored pABA prototrophy (Figure (Figure7).7). This result shows that the DIP1790 protein has ADC synthase activity but does not demonstrate ADC lyase activity because the BN1163 strain has endogenous ADC lyase (PabAc). Enzyme assays were therefore used to test DIP1790 for ADC lyase activity. BN1163 cultures harboring plasmids encoding DIP1790, Arabidopsis ADC synthase, and E. coli PabAc were grown and induced, and proteins were extracted. Extracts of cells expressing DIP1790 were incubated with chorismate and glutamine, without or with E. coli PabAc; pABA was formed in the absence of PabAc whereas, as expected, Arabidopsis PabAa-PabAb formed pABA only if PabAc was added. Reaction rates (nmol pABA h-1 mg-1 protein) were: DIP1790 – PabAc, 7.0; DIP170 + PabAc, 6.0; Arabidopsis ADCS – PabAc, <0.01; Arabidopsis ADCS + PabAc, 4.0. These data establish that DIP1790 has ADC lyase as well as ADC synthase activity.

Figure 7
Complementation of pabAa pabAb function by Corynebacterium diphtheriae DIP1790. Complementation of E. coli pabAa pabAb double mutant BN1163 by a pLOI707HE plasmid harboring DIP1790 from Corynebacterium diphtheriae on M9 minimal medium with or without ...

Another difficulty in annotating the pabAabc genes is that most organisms contain paralogs of pabAa and pabAb (trpAa and trpAb, respectively) that participate in tryptophan biosynthesis [52], and in some cases the PabAb (amidotransferase) subunit is shared between the pABA and tryptophan pathways [53]. Finally, PabAc belongs to the large branched-chain amino acid aminotransferase family (EC 2.6.1.42) and is hard to distinguish from these enzymes. These problems mean that the current SEED annotation of the pABA branch of folate synthesis should be taken as tentative. That said, analysis of the distribution of these genes reveals that most bacteria make pABA from chorismate. As expected, many intracellular bacteria lack all pabA genes. In cases where the organisms have the pterin branch but lack all enzymes of the pABA branch, annotation problems cannot be ruled out but an alternative pathway for the biosynthesis of pABA, starting for example with dehydroquinate instead of chorismate, could also be the answer [54].

Pterin salvage

The Leishmania pterin reductase PTR1 is a member of the SDR family, but has a highly characteristic motif TGX3RXG (in place of the TGX3GXG motif that is typical of this family) [55]. This motif is shared with E. coli FolM and similar SDR family proteins in a variety of bacterial taxa. Several of the folM-like genes are clustered with genes of the pterin branch of folate synthesis (Figure (Figure8),8), suggesting a function in folate or pterin synthesis Since E. coli [56] and other bacteria [57] are known to contain tetrahydromonapterin or other tetrahydropterins that could serve as cofactors for pterin-dependent enzymes, we predict that folM-like genes are not primarily involved in folate synthesis but rather are pteridine reductases that, like PTR1, produce and/or reduce 7,8-dihydropterins. (Note that such reductases are distinct from 6,7-dihydropterin reductases [also termed quinonoid pteridine reductases], of which E. coli has two [58,59].)

Figure 8
Pterin utilization by an E. coli folE deletant harboring Xylella fastidiosa PD0677 or Leishmania major PTR1. Deletant cells transformed with pBluescript (pBS) alone, or pBS harboring Xylella fastidiosa Temecula1 PD0677 or Leishmania major PTR1, were streaked ...

Consistent with this prediction, the recombinant FolM protein catalyzes reduction of dihydrobiopterin to the tetrahydro form; unlike PTR1, however, it does not mediate reduction of fully oxidized biopterin to the dihydro form [37]. Supporting the latter observation, we found that an E. coli GCHY-I mutant (which is unable to make pterins) can use the dihydro but not the oxidized forms of neopterin, monapterin, or 6-hydroxymethylpterin to support folate synthesis [60]. Futhermore, expression of a typical folM-like gene (Xylella fastidiosa PD0677, Figure Figure5D,5D, Table Table2)2) from a plasmid did not enable this mutant to use oxidized pterins, indicating that – like FolM – the PD0677 gene product does not act on oxidized pterins (Figure (Figure8).8). In control experiments in which Leishmania PTR1 was expressed from a plasmid, the mutant was able to use oxidized pterins, confirming that it is oxidized pterin reduction (and not uptake) that is lacking in E. coli (Figure (Figure8)8) [60].

Searching the Arabidopsis genome revealed some 86 members of the SDR family, of which none had the TGX3RXG motif. This led to the prediction that Arabidopsis would be unable to salvage oxidized pterins, which was verified by showing that 6-hydroxymethylpterin was not reduced in vivo or in vitro, and was not incorporated into folates [60].

Conclusion

This analysis and integration study demonstrates that simple phylogenomic analysis of a biochemical pathway – even a well-known one – can unearth globally missing (e.g., folQ) or locally missing (e.g., folE2 or folC2) genes in bacteria and plants and reveals that many open questions remain (such as the missing folQ, folB, folE cases listed in Table Table2).2). It can also identify, or suggest functions for, additional genes related to the pathway (e.g., folM). Such analysis can thus lead to discovery of potential new drug or herbicide targets such as GCHY-IB, which occurs in many pathogenic bacteria but not in mammals, or the chloroplast folate carrier that is likewise absent from mammals.

It should be noted that content of the current SEED folate subsystem captures the present status of an ongoing annotation effort, that the content will be refined and improved as more bacterial and plant genomes are added, and that further predictions are expected to emerge. Finally, we emphasize that the predictions herein are offered with the hope that others will find them useful in their own research.

Methods

Bioinformatics

Analysis of the folate subsystem was performed in the SEED database [61]. Results are available in the 'Folate biosynthesis sub-system' on the public SEED server at [62]. The snapshot of this analysis on the SEED database is given in the additional file. Phylogenetic pattern searches were made on the NMPDR SEED server at [63] to find candidates for the missing folE and folC genes. We also used the Blast tools and resources at NCBI [64] and the comparative genomics platforms STRING [65] for additional gene clustering analysis tools.

Annotations for paralog families were made using physical clustering on the chromosome when possible or by building phylogenetic trees using the ClustalW tool [66] integrated in SEED or deriving specific protein motifs. Pseudogenes (i.e., those encoding clearly aberrant proteins) were ignored; these are not uncommon in the folate pathways of intracellular parasites undergoing genome reduction [67].

The 'variant code' is used in SEED to schematize the type of pathways found in a given organism [21]. A three-digit code was used. Digit one describes the pterin branch of the pathway: 1 = complete, 0 = HPPK and DHPS missing, 4 = DHNA missing, 7 = GCHY-I missing. Digit two describes the pABA branch: 1 = two or three of the pabAabc genes present, 0 = all pabAabc genes missing or just one present. Digit three describes the salvage pathway: 1 = complete; 0 = FPGS and DHFR missing; 2 = FPGS missing; 6 = DHFR missing. Variant -1 represents genomes with no pathway genes but no need for them because no folate-dependent enzymes are present. Particular care was given to annotation of fused proteins, which are common in both branches of the pathway; SEED has annotation tools to deal with fusion proteins [21].

Strains, growth conditions, and cloning

Bacteria were routinely grown at 37°C in LB medium (BD Diagnostic Systems), in minimal medium [68] supplemented with 0.2% (v/v) glycerol, or in M9 medium [35]. Agar (BD Diagnostic Systems) concentration in plates was 15 g l-1. Transformations were by standard procedures [69,70]. Thymidine (dT, 300 μM), ampicillin (100 μg ml-1), kanamycin (50 or 100 μg ml-1), tetracycline (10 μg ml-1), isopropyl-β-D-thiogalactopyranoside (IPTG, 0.5 or 1 mM), methionine (100 μg ml-1), glycine (100 μg ml-1), pABA (0.5 μg ml-1) and L-arabinose (0.02–0.2%, w/v) were added as required. Strains Topo10 (Invitrogen), BL21-CodonPlus (DE3)-RIL (Stratagene), DH10B, or DH5α were used for cloning and expression. SF4 (F-strA recA folC srlC::Tn10) [40], BN1163 (pabA1, pab-B::Kan, rpsL704, ilvG-, rfb-50, rph-1) (B. Nichols, University of Chicago), and MG1655 (ΔfolE::KanR) [35] were used for complementation tests.

The Chlamydia trachomatis CT610 and CT611 genes were cloned in pBAD24 [71] using the following primers: CT610, 5'-AATACCATGGTGGAGGTGTTTATGAA-3' and 5'-AATAAAGCTTTTAATAAGATTGATGACAACTAC-3'; CT611, 5'-AATACCATGGAAATAACTCCGATCAAAACAC-3' and 5'-AATAAAGCTTTCATTTCTTTTCTTGACTCCAC-3'. Genomic DNA from C. trachomatis, LGV-II, strain 434 was obtained from ABI (Maryland). PCR products were obtained and purified as described [43], then digested with NcoI/HindIII before ligation into plasmid pBAD24 digested with the same enzymes and transformation into Topo10 cells (Invitrogen). The respective plasmids named pBY149.9 (expressing CT610) and pBY143.1 (expressing CT611) were checked by sequencing.

The Corynebacterium diphtheriae DIP1790 gene was cloned into pGEM-T Easy (Promega), after amplification from genomic DNA (obtained from the American Type Culture Collection) using the primers 5'-GCGGCCGCCACAGGAAACAGCTATGGTTATGCAACGCGCGCA-3' and 5'-GAGCTCTCACACTTGGGCGATATTCT-3'. The SstI site in the gene was ablated by PCR using the internal primers 5'-TCATCACCGAaCtTGAAGGCA-3' and 5'-TTTGCCTTCAaGtTCGGTGATG-3' (changed nucleotides in lower case). The modified gene was ligated into pGEM-T Easy and verified by sequencing. It was then excised with NotI and SstI and ligated into pLOI707HE [72]. This construct was used to transform E. coli BN1163. Complementation tests were made using minimal medium, appropriately supplemented as above.

The Xylella fastidiosa Temecula1 PD0677 amplicon preceded by a Shine-Dalgarno sequence and a stop codon in frame with LacZα was cloned between the EcoRI and KpnI sites of pBluescript SK-. The PCR template was genomic DNA from the American Type Culture Collection; primers were 5'-AGTCAGAATTCGTGAAGGAAACAGCTATGTCAGATCCCTCTAAAGTC-3' and 5'-AGTAGGTACCTCATGTCAGCGTGCGGCC-3'; amplification was with KOD HiFi polymerase. The deduced amino acid sequence differed from that published in having serine not isoleucine at position 57. The PTR1 construct was as described [60]. The constructs were introduced into E. coli folE deletant cells [35]. Transformants were grown on LB plates supplemented appropriately as above.

ADCS/ADCL enzyme assays

Protein extracts were prepared from IPTG-induced cultures as described [73] from strain BN1163 harboring three pLOI707HE constructs – DIP1790, Arabidopsis ADCS [73], or vector alone – and from BL21-CodonPlus (DE3)-RIL harboring E. coli PabC cloned in pJMG30 [74]. Glutamine-dependent pABA synthesis activity was assayed as described [73]. PabC extract (3.8 μg protein) was added when indicated. Assays were incubated for 1 h at 37°C, stopped with 20 μl of 75% (v/v) acetic acid, held on ice for 1 h, then stored at -80°C until analysis. pABA was estimated by HPLC with fluorescence detection [73].

Authors' contributions

BElY carried out the complementation studies on the Chlamydia cofE homolog. RDdelaG made the complementation and biochemical assays on the Corynebacterium pABA synthesis protein. AN carried out the studies on the Xylella folM homolog. VdeC-L and ADH conceived the study, carried out bioinformatic work, and drafted the manuscript. All authors read and approved the final manuscript.

Supplementary Material

Additional File 1:

Spreadsheet summarizing the distribution, clustering, and fusions of bacterial and plant folate synthesis genes. This table represents a snapshot for the record of the "BMCgenomics 2007" table that can be found in the Folate Biosynthesis sub-system on the public SEED website http://theseed.uchicago.edu/FIG/index.cgi. Clustering is shown by similar color backgrounds. Genome and protein IDs are from the SEED database. Abbreviations for the functional roles are given in the first page of the spreadsheet, gene distribution in all analyzed genomes in the second page. Note that the SEED table is the primary source to which the reader is directed; it is not static but develops with time as new genomes become available and predictions are tested and validated.

Acknowledgements

We thank B. Shane for the SF4 strain, B. Nichols for the BN1163 strain, H. Mori for the pfolCEc plasmid, and G. Basset for help with enzyme assays. We thank A. Osterman for insightful suggestions and help with figure design. This work was supported in part by National Institutes of Health grants R01 GM70641-01 (to V. de C.-L.) and R01 GM071382 (to A.D.H.), by the National Research Initiative of the USDA Cooperative State Research, Education, and Extension Service, grant number 2005-35318-15228 (to A.D.H.), and by an endowment from the C.V. Griffin, Sr. Foundation.

References

  • Matthews RG. One-carbon metabolism. In: Neidhardt FC, Curtiss R 3rd, Ingraham JL, Lin ECC, Low KB, Magasanik B, Reznikoff WS, Riley M, Schaechter M, Umbarger HE, editor. Escherichia coli and Salmonella: Cellular and Molecular Biology. Second. Vol. 1. Washington DC, ASM Press; 1996. pp. 600–611.
  • Huennekens FM, Vitols KS, Pope LE, Fan J. Membrane transport of folate compounds. J Nutr Sci Vitaminol (Tokyo) 1992:52–57.
  • Lucock M. Folic acid: nutritional biochemistry, molecular biology, and role in disease processes. Mol Genet Metab. 2000;71:121–138. doi: 10.1006/mgme.2000.3027. [PubMed] [Cross Ref]
  • Cossins EA, Chen L. Folates and one-carbon metabolism in plants and fungi. Phytochemistry. 1997;45:437–452. doi: 10.1016/S0031-9422(96)00833-3. [Cross Ref]
  • Green JC, Nichols BP, Matthews RG. Folate biosynthesis, reduction, and polyglutamylation. In: Neidhardt FC, Curtiss R 3rd, Ingraham JL, Lin ECC, Low KB, Magasanik B, Reznikoff WS, Riley M, Schaechter M, Umbarger HE, editor. Escherichia coli and Salmonella: Cellular and Molecular Biology. Second. Vol. 1. Washington DC, ASM Press; 1996. pp. 665–673.
  • Hanson AD, Gregory JF., 3rd Synthesis and turnover of folates in plants. Curr Opin Plant Biol. 2002;5:244–249. doi: 10.1016/S1369-5266(02)00249-2. [PubMed] [Cross Ref]
  • Hyde JE. Exploring the folate pathway in Plasmodium falciparum. Acta Trop. 2005;94:191–206. [PMC free article] [PubMed]
  • Huovinen P, Sundström L, Swedberg G, Sköld O. Trimethoprim and sulfonamide resistance. Antimicrob Agents Chemother. 1995;39:279–289. [PMC free article] [PubMed]
  • Sybesma W, Starrenburg M, Kleerebezem M, Mierau I, de Vos WM, Hugenholtz J. Increased production of folate by metabolic engineering of Lactococcus lactis. Appl Environ Microbiol. 2003;69:3069–3076. doi: 10.1128/AEM.69.6.3069-3076.2003. [PMC free article] [PubMed] [Cross Ref]
  • Hossain T, Rosenberg I, Selhub J, Kishore G, Beachy R, Schubert K. Enhancement of folates in plants through metabolic engineering. Proc Natl Acad Sci USA. 2004;101:5158–5163. doi: 10.1073/pnas.0401342101. [PubMed] [Cross Ref]
  • Díaz de la Garza R, Quinlivan EP, Klaus SM, Basset GJ, Gregory JF, 3rd, Hanson AD. Folate biofortification in tomatoes by engineering the pteridine branch of folate synthesis. Proc Natl Acad Sci USA. 101:13720–13725. doi: 10.1073/pnas.0404208101. [Cross Ref]
  • Klaus SM, Wegkamp A, Sybesma W, Hugenholtz J, Gregory JF, 3rd, Hanson AD. A nudix enzyme removes pyrophosphate from dihydroneopterin triphosphate in the folate synthesis pathway of bacteria and plants. J Biol Chem. 2005;280:5274–5280. doi: 10.1074/jbc.M413759200. [PubMed] [Cross Ref]
  • Haussmann C, Rohdich F, Schmidt E, Bacher A, Richter G. Biosynthesis of pteridines in Escherichia coli. Structural and mechanistic similarity of dihydroneopterin-triphosphate epimerase and dihydroneopterin aldolase. J Biol Chem. 1998;273:17418–17424. doi: 10.1074/jbc.273.28.17418. [Cross Ref]
  • Ferone R, Singer SC, Hunt DF. In vitro synthesis of alpha-carboxyl-linked folylpolyglutamates by an enzyme preparation from Escherichia coli. J Biol Chem. 1986;261:16363–16371. [PubMed]
  • Ravanel S, Cherest H, Jabrin S, Grunwald D, Surdin-Kerjan Y, Douce R, Rebeille F. Tetrahydrofolate biosynthesis in plants: molecular and functional characterization of dihydrofolate synthetase and three isoforms of folylpolyglutamate synthetase in Arabidopsis thaliana. Proc Natl Acad Sci USA. 2001;98:15360–15365. doi: 10.1073/pnas.261585098. [Cross Ref]
  • Orsomando G, de la Garza RD, Green BJ, Peng M, Rea PA, Ryan TJ, Gregory JF, 3rd, Hanson AD. Plant gamma-glutamyl hydrolases and folate polyglutamates: characterization, compartmentation, and co-occurrence in vacuoles. J Biol Chem. 2005;280:28877–28884. doi: 10.1074/jbc.M504306200. [PubMed] [Cross Ref]
  • Nare B, Hardy LW, Beverley SM. The roles of pteridine reductase 1 and dihydrofolate reductase-thymidylate synthase in pteridine metabolism in the protozoan parasite Leishmania major. J Biol Chem. 1997;272:13883–13891. doi: 10.1074/jbc.272.21.13883. [PubMed] [Cross Ref]
  • Orsomando G, Bozzo GG, de la Garza RD, Basset GJ, Quinlivan EP, Naponelli V, Rébeillé F, Ravanel S, Gregory JF, 3rd, Hanson AD. Evidence for folate-salvage reactions in plants. Plant J. 2006;46:426–435. doi: 10.1111/j.1365-313X.2006.02685.x. [Cross Ref]
  • Green JM, Ballou DP, Matthews RG. Examination of the role of methylenetetrahydrofolate reductase in incorporation of methyltetrahydrofolate into cellular metabolism. FASEB J. 1988;2:42–47. [PubMed]
  • Toy J, Bognar AL. Cloning and expression of the gene encoding Lactobacillus casei folylpoly-gamma-glutamate synthetase in Escherichia coli and determination of its primary structure. J Biol Chem. 1990;265:2492–2499. [PubMed]
  • Overbeek R, Begley T, Butler RM, Choudhuri JV, Chuang HY, Cohoon M, de Crécy-Lagard V, Diaz N, Disz T, Edwards R, et al. The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes. Nucleic Acids Res. 2005;33:5691–5702. doi: 10.1093/nar/gki866. [PMC free article] [PubMed] [Cross Ref]
  • Galperin MY, Brenner SE. Using metabolic pathway databases for functional annotation. Trends Genet. 1998;14:332–333. doi: 10.1016/S0168-9525(98)01523-6. [PubMed] [Cross Ref]
  • Selkov E, Maltsev N, Olsen GJ, Overbeek R, Whitman WB. A reconstruction of the metabolism of Methanococcus jannaschii from sequence data. Gene. 1997;197:GC11–26. doi: 10.1016/S0378-1119(97)00307-7. [Cross Ref]
  • Gerlt JA, Babbitt PC. Can sequence determine function? Genome Biol. 2000;1:REVIEWS0005. doi: 10.1186/gb-2000-1-5-reviews0005. [PMC free article] [PubMed] [Cross Ref]
  • Karp PD. Call for an enzyme genomics initiative. Genome Biol. 2004;5:401. doi: 10.1186/gb-2004-5-8-401. [PMC free article] [PubMed] [Cross Ref]
  • Osterman A, Overbeek R. Missing genes in metabolic pathways: a comparative genomics approach. Curr Opin Chem Biol. 2003;7:238–251. doi: 10.1016/S1367-5931(03)00027-9. [PubMed] [Cross Ref]
  • Koonin EV, Mushegian AR, Bork P. Non-orthologous gene displacement. Trends Genet. 1996;12:334–336. doi: 10.1016/0168-9525(96)20010-1. [PubMed] [Cross Ref]
  • Clark BF, Marcker KA. The role of N-formyl-methionyl-sRNA in protein biosynthesis. J Mol Biol. 1966;17:394–406. [PubMed]
  • Vasconcelos AT, Ferreira HB, Bizarro CV, Boniato SL, Carvalho MO, Pinto PM, Almeida DF, Almeida LG, Almeida R, Alves-Filho L, et al. Swine and poultry pathogens: the complete genome sequences of two strains of Mycoplasma hyopneumoniae and a strain of Mycoplasma synoviae. J Bacteriol. 2005;187:5568–5577. doi: 10.1128/JB.187.16.5568-5577.2005. [PMC free article] [PubMed] [Cross Ref]
  • Schirch L, Gross T. Serine transhydroxymethylase. Identification as the threonine and allothreonine aldolases. J Biol Chem. 1968;243:5651–5655. [PubMed]
  • Myllykallio H, Lipowski G, Leduc D, Filee J, Forterre P, Liebl U. An alternative flavin-dependent mechanism for thymidylate synthesis. Science. 2002;297:105–107. doi: 10.1126/science.1072113. [Cross Ref]
  • Myllykallio H, Leduc D, Filee J, Liebl U. Life without dihydrofolate reductase FolA. Trends Microbiol. 2003;11:220–223. [PubMed]
  • Henderson GB, Zevely EM, Huennekens FM. Purification and properties of a membrane-associated, folate-binding protein from Lactobacillus casei. J Biol Chem. 1977;252:3760–3765.
  • Kumar HP, Tsuji JM, Henderson GB. Folate transport in Lactobacillus salivarius. Characterization of the transport mechanism and purification and properties of the binding component. J Biol Chem. 1987;262:7171–7179. [PubMed]
  • Klaus SM, Kunji ER, Bozzo GG, Noiriel A, de la Garza RD, Basset GJ, Ravanel S, Rébeillé F, Gregory JF, 3rd, Hanson AD. Higher plant plastids and cyanobacteria have folate carriers related to those of trypanosomatids. J Biol Chem. 2005;280:38457–38463. doi: 10.1074/jbc.M507432200. [Cross Ref]
  • Rodionov DA, Hebbeln P, Gelfand MS, Eitinger T. Comparative and functional genomic analysis of prokaryotic nickel and cobalt uptake transporters: evidence for a novel group of ATP-binding cassette transporters. J Bacteriol. 2006;188:317–327. doi: 10.1128/JB.188.1.317-327.2006. [PMC free article] [PubMed] [Cross Ref]
  • Giladi M, Altman-Price N, Levin I, Levy L, Mevarech M. FolM, a new chromosomally encoded dihydrofolate reductase in Escherichia coli. J Bacteriol. 2003;185:7015–7018. doi: 10.1128/JB.185.23.7015-7018.2003. [PMC free article] [PubMed] [Cross Ref]
  • Levin I, Mevarech M, Palfey BA. Characterization of a novel bifunctional dihydropteroate synthase/dihydropteroate reductase enzyme from Helicobacter pylori. J Bacteriol. 2007;189:4062–4069. doi: 10.1128/JB.01878-06. [PMC free article] [PubMed] [Cross Ref]
  • Graziani S, Xia Y, Gurnon JR, Van Etten JL, Leduc D, Skouloubris S, Myllykallio H, Liebl U. Functional analysis of FAD-dependent thymidylate synthase ThyX from Paramecium bursaria chlorella virus-1. J Biol Chem. 2004;279:54340–54347. doi: 10.1074/jbc.M409121200. [PubMed] [Cross Ref]
  • Bognar AL, Osborne C, Shane B, Singer SC, Ferone R. Folylpoly-gamma-glutamate synthetase-dihydrofolate synthetase. Cloning and high expression of the Escherichia coli folC gene and purification and properties of the gene product. J Biol Chem. 1985;260:5625–5630.
  • Mathieu M, Debousker G, Vincent S, Viviani F, Bamas-Jacques N, Mikol V. Escherichia coli FolC structure reveals an unexpected dihydrofolate binding site providing an attractive target for anti-microbial therapy. J Biol Chem. 2005;280:18916–18922. doi: 10.1074/jbc.M413799200. [PubMed] [Cross Ref]
  • Neale GA, Mitchell A, Finch LR. Formylation of methionyl-transfer ribonucleic acid in Mycoplasma mycoides subsp. mycoides. J Bacteriol. 1981;146:816–818.
  • El Yacoubi B, Bonnett S, Anderson JN, Swairjo MA, Iwata-Reuyl D, de Crécy-Lagard V. Discovery of a new prokaryotic type I GTP cyclohydrolase family. J Biol Chem. 2006 [PubMed]
  • Fan H, Brunham RC, McClarty G. Acquisition and synthesis of folates by obligate intracellular bacteria of the genus Chlamydia. J Clin Invest. 1992;90:1803–1811. [PMC free article] [PubMed]
  • Schwarzenbacher R, Stenner-Liewen F, Liewen H, Robinson H, Yuan H, Bossy-Wetzel E, Reed JC, Liddington RC. Structure of the Chlamydia protein CADD reveals a redox enzyme that modulates host cell apoptosis. J Biol Chem. 2004;279:29320–29324. doi: 10.1074/jbc.M401268200. [PubMed] [Cross Ref]
  • Suzuki Y, Brown GM. The biosynthesis of folic acid. XII. Purification and properties of dihydroneopterin triphosphate pyrophosphohydrolase. J Biol Chem. 1974;249:2405–2410.
  • Bessman MJ, Frick DN, O'Handley SF. The MutT proteins or "Nudix" hydrolases, a family of versatile, widely distributed, "housecleaning" enzymes. J Biol Chem. 1996;271:25059–25062. doi: 10.1074/jbc.271.41.25059. [PubMed] [Cross Ref]
  • Tucker AM, Winkler HH, Driskell LO, Wood DO. S-Adenosylmethionine transport in Rickettsia prowazekii. J Bacteriol. 2003;185:3031–3035. doi: 10.1128/JB.185.10.3031-3035.2003. [Cross Ref]
  • Williams DL, Spring L, Harris E, Roche P, Gillis TP. Dihydropteroate synthase of Mycobacterium leprae and dapsone resistance. Antimicrob Agents Chemother. 2000;44:1530–1537. doi: 10.1128/AAC.44.6.1530-1537.2000. [PMC free article] [PubMed] [Cross Ref]
  • Li H, Graupner M, Xu H, White RH. CofE catalyzes the addition of two glutamates to F420-0 in F420 coenzyme biosynthesis in Methanococcus jannaschii. Biochemistry. 2003;42:9771–9778. doi: 10.1021/bi034779b. [PubMed] [Cross Ref]
  • Kitagawa M, Ara T, Arifuzzaman M, Ioka-Nakamichi T, Inamoto E, Toyonaga H, Mori H. Complete set of ORF clones of Escherichia coli ASKA library (a complete set of E. coli K-12 ORF archive): unique resources for biological research. DNA Res. 2005;12:291–299. doi: 10.1093/dnares/dsi012. [PubMed] [Cross Ref]
  • Xie G, Keyhani NO, Bonner CA, Jensen RA. Ancient origin of the tryptophan operon and the dynamics of evolutionary change. Microbiol Mol Biol Rev. 2003;67:303–342. doi: 10.1128/MMBR.67.3.303-342.2003. [PMC free article] [PubMed] [Cross Ref]
  • Yanofsky C. Advancing our knowledge in biochemistry, genetics, and microbiology through studies on tryptophan metabolism. Annu Rev Biochem. 2001;70:1–37. doi: 10.1146/annurev.biochem.70.1.1. [PubMed] [Cross Ref]
  • Porat I, Sieprawaska-Lupa M, Teng Q, Bohanon FJ, White RH, Whitman WB. Biochemical and genetic characterization of an early step in a novel pathway for the biosynthesis of aromatic amino acids and p-aminobenzoic acid in the archaeon Methanococcus maripaludis. Mol Microbiol. 2006;62:1117–1131. doi: 10.1111/j.1365-2958.2006.05426.x. [Cross Ref]
  • Gourley DG, Schuttelkopf AW, Leonard GA, Luba J, Hardy LW, Beverley SM, Hunter WN. Pteridine reductase mechanism correlates pterin metabolism with drug resistance in trypanosomatid parasites. Nat Struct Biol. 2001;8:521–525. doi: 10.1038/88584. [PubMed] [Cross Ref]
  • Ikemoto K, Sugimoto T, Murata S, Tazawa M, Nomura T, Ichinose H, Nagatsu T. (6R)-5,6,7,8-tetrahydro-L-monapterin from Escherichia coli, a novel natural unconjugated tetrahydropterin. Biol Chem. 2002;383:325–330. doi: 10.1515/BC.2002.035. [PubMed] [Cross Ref]
  • Lee HW, Oh CH, Geyer A, Pfleiderer W, Park YS. Characterization of a novel unconjugated pteridine glycoside, cyanopterin, in Synechocystis sp. PCC 6803. Biochim Biophys Acta. 1999;1410:61–70. doi: 10.1016/S0005-2728(98)00175-3. [PubMed] [Cross Ref]
  • Vasudevan SG, Shaw DC, Armarego WL. Dihydropteridine reductase from Escherichia coli. Biochem J. 1988;255:581–588. [PubMed]
  • Vasudevan SG, Armarego WL, Shaw DC, Lilley PE, Dixon NE, Poole RK. Isolation and nucleotide sequence of the hmp gene that encodes a haemoglobin-like protein in Escherichia coli K-12. Mol Gen Genet. 1991;226:49–58. doi: 10.1007/BF00273586. [PubMed] [Cross Ref]
  • Noiriel A, Naponelli V, Gregory JF, 3rd, Hanson AD. Pterin and folate salvage: plants and Escherichia coli lack capacity to reduce oxidized pterins. Plant Physiol. 2007;143:1101–1109. doi: 10.1104/pp.106.093633. [PubMed] [Cross Ref]
  • http://anno-3.nmpdr.org/anno/FIG/subsys.cgi
  • http://theseed.uchicago.edu/FIG/index.cgi
  • http://www.nmpdr.org/FIG/sigs.cgi?SPROUT=1
  • Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–3402. doi: 10.1093/nar/25.17.3389. [PMC free article] [PubMed] [Cross Ref]
  • von Mering C, Jensen LJ, Snel B, Hooper SD, Krupp M, Foglierini M, Jouffre N, Huynen MA, Bork P. STRING: known and predicted protein-protein associations, integrated and transferred across organisms. Nucleic Acids Res. 2005;33:D433–D437. doi: 10.1093/nar/gki005. [PMC free article] [PubMed] [Cross Ref]
  • Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, Higgins DG, Thompson JD. Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res. 2003;31:3497–3500. doi: 10.1093/nar/gkg500. [PMC free article] [PubMed] [Cross Ref]
  • Davis RE, Jomantiene R, Zhao Y. Lineage-specific decay of folate biosynthesis genes suggests ongoing host adaptation in phytoplasmas. DNA Cell Biol. 2005;24:832–840. doi: 10.1089/dna.2005.24.832. [PubMed] [Cross Ref]
  • Richaud C, Mengin-Lecreulx D, Pochet S, Johnson EJ, Cohen GN, Marlière P. Directed evolution of biosynthetic pathways. Recruitment of cysteine thioethers for constructing the cell wall of Escherichia coli. J Biol Chem. 1993;268:26827–26835.
  • Miller JH. Experiments in Molecular Genetics. Cold Spring Harbor, Cold Spring Harbor Press; 1972.
  • Sambrook J, Fitsch EF, Maniatis T. Molecular Cloning: A Laboratory Manual. Cold Spring Harbor, Cold Spring Harbor Press; 1989.
  • Guzman LM, Belin D, Carson MJ, Beckwith J. Tight regulation, modulation, and high-level expression by vectors containing the arabinose PBAD promoter. J Bacteriol. 1995;177:4121–4130. [PMC free article] [PubMed]
  • Arfman N, Worrell V, Ingram LO. Use of the tac promoter and lacIq for the controlled expression of Zymomonas mobilis fermentative genes in Escherichia coli and Zymomonas mobilis. J Bacteriol. 1992;174:7370–7378. [PMC free article] [PubMed]
  • Basset GJ, Ravanel , Quinlivan EP, White R, Giovannoni JJ, Rébeillé F, Nichols BP, Shinozaki K, Seki M, Gregory JF, 3rd, Hanson AD. Folate synthesis in plants: the last step of the p-aminobenzoate branch is catalyzed by a plastidial aminodeoxychorismate lyase. Plant J. 2004;40:453–461. doi: 10.1111/j.1365-313X.2004.02231.x. [PubMed] [Cross Ref]
  • Green JM, Merkel WK, Nichols BP. Characterization and sequence of Escherichia coli pabC, the gene encoding aminodeoxychorismate lyase, a pyridoxal phosphate-containing enzyme. J Bacteriol. 1992;174:5317–5323. [PMC free article] [PubMed]

Articles from BMC Genomics are provided here courtesy of BioMed Central