Skip Navigation


MBE Advance Access originally published online on November 12, 2008
Molecular Biology and Evolution 2009 26(2):399-406; doi:10.1093/molbev/msn261
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary Data
Right arrow All Versions of this Article:
26/2/399    most recent
msn261v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Lima, W. C.
Right arrow Articles by Menck, C. F.M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Lima, W. C.
Right arrow Articles by Menck, C. F.M.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2008. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oxfordjournals.org

Research Articles

NAD Biosynthesis Evolution in Bacteria: Lateral Gene Transfer of Kynurenine Pathway in Xanthomonadales and Flavobacteriales

Wanessa C. Lima*, Alessandro M. Varani{dagger} and Carlos F.M. Menck*

* Department of Microbiology, Institute of Biomedical Sciences, University of São Paulo, São Paulo, Brazil
{dagger} Department of Botany, Institute of Biosciences, University of São Paulo, São Paulo, Brazil

E-mail: cfmmenck{at}usp.br.


    Abstract
 TOP
 Abstract
 Background
 Methods
 Results and Discussion
 Conclusion
 Supplementary Material
 Acknowledgements
 References
 
The biosynthesis of quinolinate, the de novo precursor of nicotinamide adenine dinucleotide (NAD), may be performed by two distinct pathways, namely, the bacterial aspartate (aspartate-to-quinolinate) and the eukaryotic kynurenine (tryptophan-to-quinolinate). Even though the separation into eukaryotic and bacterial routes is long established, recent genomic surveys have challenged this view, because certain bacterial species also carry the genes for the kynurenine pathway. In this work, both quinolinate biosynthetic pathways were investigated in the Bacteria clade and with special attention to Xanthomonadales and Bacteroidetes, from an evolutionary viewpoint. Genomic screening has revealed that a small number of bacterial species possess some of the genes for the kynurenine pathway, which is complete in the genus Xanthomonas and in the order Flavobacteriales, where the aspartate pathway is absent. The opposite pattern (presence of the aspartate pathway and absence of the kynurenine pathway) in close relatives (Xylella ssp. and the order Bacteroidales, respectively) points to the idea of a recent acquisition of the kynurenine pathway through lateral gene transfer in these bacterial groups. In fact, sequence similarity comparison and phylogenetic reconstruction both suggest that at least part of the genes of the kynurenine pathway in Xanthomonas and Flavobacteriales is shared by eukaryotes. These results reinforce the idea of the role that lateral gene transfer plays in the configuration of bacterial genomes, thereby providing alternative metabolic pathways, even with the replacement of primary and essential cell functions, as exemplified by NAD biosynthesis.

Key Words: lateral gene transfer • Xanthomonas • Flavobacteriales • NAD metabolism • kynurenine pathway • tryptophan catabolism


    Background
 TOP
 Abstract
 Background
 Methods
 Results and Discussion
 Conclusion
 Supplementary Material
 Acknowledgements
 References
 
For several decades, it has been known that nicotinamide adenine dinucleotide (NAD) plays a major role as a coenzyme in numerous oxidation–reduction reactions. According to a widely accepted concept, a combination of de novo and salvage pathways contributes to the biosynthesis of NAD (Rongvaux et al. 2003Go). Depending on the organism, three different building blocks can be used for the biosynthesis of NAD, namely, quinolinic acid (or quinolinate) as the precursor in the so-called de novo pathway, nicotinic acid, and nicotinamide in the salvage pathway (Mattevi 2006Go).

Two different biosynthetic pathways are responsible for de novo quinolinate production (Begley et al. 2001Go) (fig. 1). In prokaryotes, quinolinate is usually formed from aspartate and dihydroxyacetone phosphate, in a two-step pathway involving the enzymes NadB (L-aspartate oxidase, EC 1.4.3.16 [EC] , encoded by nadB, which oxidizes the amino group of L-aspartate to the imino group, forming {alpha}-iminoaspartate) and NadA (quinolinate synthase, encoded by nadA, responsible for the production of quinolinic acid from {alpha}-iminoaspartate and dihydroxyacetone phosphate) (Katoh and Hashimoto 2004Go).


Figure 1
View larger version (21K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 1.— Schematic representation of the two biosynthetic pathways to quinolinic acid: the five-step "eukaryotic" route (from L-tryptophan) and the two-step "prokaryotic" route (from L-aspartate). On the left, the alternative anthranilate branch. The enzyme name and synonyms and the EC number are shown in each text box.

 
In eukaryotes, NAD biosynthesis is coupled to tryptophan catabolism through the so-called kynurenine pathway (or tryptophan-to-quinolinate pathway), whereas the common enzymatic substrate resides in {alpha}-amino-β-carboxymuconic acid {epsilon}-semialdehyde. This compound may either cyclize nonenzymatically to yield quinolinate or it can be enzymatically decarboxylated by a-amino-b-carboxy muconic acid e-semialdehyde decarboxylase (Fukuoka et al. 2002Go). In the first step of the kynurenine pathway, tryptophan is cleaved by tryptophan-2,3-dioxygenase (TDO, EC 1.13.11.11 [EC] , encoded by kynA) to yield N-formylkynurenine. Alternatively, in some organisms (as in Fungi), this step may be performed by indoleamine 2,3-dioxygenase (EC 1.13.11.52 [EC] ). Removal of the formyl group by N-formylkynurenine formamidase (or arylformamidase, KFA, EC 3.5.1.9 [EC] , encoded by kynB) results in L-kynurenine, which is then hydroxylated by kynurenine-3-monooxygenase (KMO, EC 1.14.13.9 [EC] , encoded by kmo), a FAD-dependent monooxygenase that incorporates molecular oxygen into 3-hydroxykynurenine. Cleavage of the amino acid side chain by the pyridoxal-5'-phosphate-dependent enzyme kynureninase (KYN, EC 3.7.1.3 [EC] , encoded by kynU) generates 3-hydroxyanthranilic acid, which is then converted by 3-hydroxyanthranilate-3,4-dioxygenase (HAD, EC 1.13.11.6 [EC] , encoded by haaO) to an unstable aliphatic compound, {alpha}-amino-β-carboxymuconic acid {epsilon}-semialdehyde, which spontaneously cyclizes to quinolinic acid (Kurnasov, Goral et al. 2003Go).

Pioneer studies dating from four to five decades ago have demonstrated that certain bacteria, such as Streptomyces antibioticus, Pseudomonas fluorescens, and Xanthomonas pruni, might utilize this pathway to synthesize quinolinic acid (Davis et al. 1951Go; Wilson and Henderson 1963Go; Tremblay et al. 1967Go) or at least present related enzymatic activities, as in Pseudomonas aureofaciens and Streptomyces parvulus (Salcher and Lingens 1980Go; Brown et al. 1986Go). In a keynote work undertaken by Brown and Wagner (1970)Go, enzymatic assays of the five steps in the kynurenine pathway were performed in X. pruni, whereby all the activities were found to be present.

However, until recently, very little was known regarding identification of the primary sequence or structural prediction of such bacterial enzymes (Kurnasov, Jablonski, et al. 2003Go), and no glimpse was given to clarify the evolutionary scenario of the quinolinate biosynthetic pathway. Recent extensive bacterial genome sequencing allowed us to address this evolutionary question. Analyses of these pathways help to understand how both de novo pathways are distributed and how they have evolved within bacterial clades, when taking into special consideration two groups, Xanthomonadales and Bacteroidetes.


    Methods
 TOP
 Abstract
 Background
 Methods
 Results and Discussion
 Conclusion
 Supplementary Material
 Acknowledgements
 References
 
Identification of the Quinolinate Biosynthetic Pathway
Sequence similarity analyses were performed (using BlastP and TBlastN programs) using a current data set of more than 600 complete or almost complete genomes (locally downloaded on February 7, 2008), from NCBI (Pruitt et al. 2005Go), Expasy Proteomics Server (Gasteiger et al. 2003Go), or CMR-TIGR (Peterson et al. 2001Go) databases. Metabolic pathways were analyzed through the KEGG web service (Kanehisa et al. 2008Go), BRENDA enzyme database (Schomburg et al. 2004Go), Genome Properties server (Selengut et al. 2007Go), and SEED Subsystems (Overbeek et al. 2005Go), up to April 2008. Genes coding for enzymes involved in the quinolinate biosynthesis were identified both by similarity (using as seed sequences homologs with a functional assignment [see supplementary table S1, Supplementary Material online]; threshold inclusion values were e-value <10–05 and similarity percent higher than 40%) or keyword searches and, wherever possible, by genome context analysis (for a review, see Osterman and Overbeek 2003Go). GenBank accession numbers for the most relevant genomes used in this study are presented in supplementary table S2, Supplementary Material online. Data concerning the other genomes here analyzed are available upon request to the authors. Protein sequences to all open reading frames (ORFs) referred in table 1 are provided as Supplementary Material online.


View this table:
[in this window]
[in a new window]

 
Table 1 Genes Related to Quinolinate Biosynthesis in Xanthomonadales and Bacteroidetes Genomes

 
Phylogenetic Analyses
Protein sequences were aligned by using ClustalX 2.0 program (Larkin et al. 2007Go), whereby regions of the alignments that were ambiguous, hypervariable, or containing gaps were excluded from subsequent analysis (GENEDOC program, Nicholas et al. 1997Go). Pairwise genetic distances were computed by using a PROTDIST program from the PHYLIP 3.67 package (Felsenstein 2005Go) or MEGA 4.0 (Tamura et al. 2007Go), with the Jones–Taylor–Thornton (JTT) model of amino acid replacement, gamma-distributed rate (four categories) and inferred {alpha} parameter. Distance-based phylogenetic trees were generated through the Neighbor-Joining (NJ) algorithm (NJ program). Bootstrap assessment of tree topology (1,000 replicates) was done with the SEQBOOT program. Maximum likelihood (ML) trees were set up with RAxML 7.0.4 (Stamatakis 2006Go), with the JTT model, and parameters for invariable sites and gamma-distributed rate heterogeneity (four categories). One-thousand bootstrap replicates were executed and bootstrap values drawn up on the best-scoring ML tree. Trees were visualized using MEGA 4.0 and TREEVIEW (Page 1996Go) programs and were arbitrarily rooted at midpoint (although trees should be fundamentally viewed as unrooted).


    Results and Discussion
 TOP
 Abstract
 Background
 Methods
 Results and Discussion
 Conclusion
 Supplementary Material
 Acknowledgements
 References
 
Identification of the Quinolinate Biosynthetic Pathway in Bacteria
In a review, Kurnasov, Goral, et al. (2003)Go described the presence of the kynurenine pathway in certain bacterial species (including Cytophaga hutchinsonii, Gemmata sp., Polaribacter filamentus, Ralstonia metallidurans, and Xanthomonas axonopodis), thus providing evidence that this pathway is found not only in Eukaryotes as formerly thought. However, since then hundreds of new genomes have been sequenced and published. Consequently, an up-to-date analysis is required to afford a broader evolutionary view of NAD biosynthesis in bacteria.

Attempts have been made to identify all the genes belonging to both the aspartate and kynurenine pathways for the more than 600 genomes publicly available, figure 2 representing a summary of all the collected data. According to current assumptions, a large majority of genomes bear only those genes responsible for the aspartate pathway (nadA and nadB). However, the presence of kynurenine pathway–coding genes in bacterial groups is not restricted as previously thought. In fact, six groups of Bacteria carry at least some of the genes for this pathway, although patchily distributed within each group: Actinobacteria (from the order Actinomycetales), Firmicutes (only species from the order Bacillales), Bacteroidetes (only Flavobacteriales), Alpha-proteobacteria (Rhizobiales and Rhodobacterales), Beta-proteobacteria (Burkholderiales), and Gamma-proteobacteria (Xanthomonadales and some Pseudomonadales). Most of these groups bear only kynA, kynB, and kynU genes (fig. 2), kmo and haaO having a very restricted phyletic distribution and being present only in the suborder Micromonosporineae and the orders Flavobacteriales and Xanthomonadales.


Figure 2
View larger version (30K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 2.— Schematic profile of the presence and absence of seven genes of the kynurenine (kynA, kynB, kmo, kynU, and haaO) and aspartate pathway (nadA and nadB). A consensus phylogenetic tree for the major bacterial groups is presented. Exceptions within each group are as follows: (a) species bearing kynA, kynB, and kynU: Herpetosiphon aurantiacus (Chlorofexi); Deinococcus ssp. (Deinococcus-Thermus); Kineococcus radiotolerans and Rhodococcus sp. RHA1 (Actinobacteria); Bradyrhizobium ssp., Mesorhizobium loti, Sphingopyxis alaskensis, and Erythrobacter litoralis (alpha-proteobacteria); Pseudomonas aeruginosa PAO1; Pseudomonas fluorescens Pf-5 and Photorhabdus luminescens TTO1 (Gamma-proteobacteria); (b) species bearing kynA or kynB and kynU: Renibacterium salmoninarum, Arthrobacter ssp., Thermobifida fusca, and Saccharopolyspora erythraea (Actinobacteria); Anaeromyxobacter and Bdellovibrio (Delta-proteobacteria); Acinetobacter baumannii, Idiomarina loihiensis, Nitrococcus mobilis, and Marinomonas sp. (Gamma-proteobacteria).

 
The presence of kynA, kynB, and kynU in few bacterial species has been linked to tryptophan catabolism by way of the anthranilate pathway (fig. 1) (Farrow and Pesci 2007Go). In Escherichia coli and many other bacteria, tryptophan degradation proceeds through pyridoxal phosphate–dependent tryptophanase TnaA (encoded by tnaA) (Vederas et al. 1978Go). In eukaryotes, catabolism of tryptophan is coupled to NAD biosynthesis through the kynurenine pathway, using as common intermediate 2-amino-3-carboxymuconate semialdehyde (Fukuoka et al. 2002Go). However, in some species of Pseudomonadaceae and Bacillaceae, the oxidative degradation of L-tryptophan occurs through the anthranilate pathway by using TDO, KFA, and KYN enzymes (encoded by kynA, kynB, and kynU genes, respectively), whereby tryptophan is converted to formylkynurenine and then to kynurenine by TDO and KFA, with the final conversion of kynurenine to anthranilic acid catalyzed by the kynureninase KYN (fig. 1) (D'Argenio et al. 2002Go; Kurnasov et al. 2003Go; Farrow and Pesci 2007Go). Interestingly, the anthranilate branch and the aspartate-to-quinolinate pathway always co-occur in the given species, playing different roles within the cell, such as quinolinate production and tryptophan degradation.

Nevertheless, only the order Flavobacteriales and Xanthomonas present the entire set of genes responsible for the complete kynurenine pathway (kynA, kynB, kmo, kynU, and haaO) and not those for the classical prokaryotic pathway (comprising nadA and nadB genes). The order Xanthomonadales is mainly composed of three genera (Xanthomonas, Stenotrophomonas, and Xylella), although only Xanthomonas and Stenotrophomonas bear genes for the kynurenine pathway, as Xylella shows the opposite profile by carrying only nadA and nadB (table 1). The same occurs for Bacteroidetes: Although species from the order Flavobacteriales reveal the same pattern seen in Xanthomonas, those from the order Bacteroidales bear only genes for the aspartate-to-quinolinate pathway (table 1).

Genomic Organization of the Kynurenine Pathway in Flavobacteriales and Xanthomonas
Genes responsible for the kynurenine pathway are normally found clustered along the genome. In Xanthomonas, the genetic arrangement is similar for all completely sequenced genomes (supplementary fig. S1, Supplementary Material online): Three genes are clustered together in a putative operon—kmo, kynU, and haaO—interrupted by one or two hypothetical genes, depending on the organism. The genes kynA and kynB do not belong to the same cluster and may be separated from the main operon by more than 1 Mb. This configuration gives support to the idea that these three genes were transferred en bloc, given their common evolutionary history (see below). Moreover, this pathway has been previously identified as part of a genomic island potentially originating from lateral gene transfer (LGT) in Xanthomonas (Lima et al. 2005Go, 2008Go).

Very few genomes display genes spread throughout the entire chromosome, this being the case of Flavobacteriales. Those related to the kynurenine pathway are not directly connected physically, kmo and kynU being the only ones closer together in some cases (even though separated by at least 5,000 bp). Moreover, there is no collinearity: each gene appearing in a different position, depending on the genome considered (supplementary fig. S1, Supplementary Material online, shows only the scheme for completely sequenced Flavobacteriales genomes). This fact may reflect gene dispersion throughout the genome by recombination, although this is difficult to prove, and other operons are easily found in the genome of these bacteria (data not shown).

Identification of the kynB gene has been a matter of debate, because it occurs in at least two nonorthologous forms, one of which presenting a cyclase or metal-dependent hydrolase domain, functionally characterized in some bacteria (as in Pseudomonas aeruginosa, R. metallidurans, and Bacillus cereus) (Kurnasov, Goral, et al. 2003Go; Kurnasov, Jablonski, et al. 2003Go) and the other possessing a conserved esterase/lipase domain (from the alpha/beta-hydrolase family), characterized as the kynurenine formamidase in eukaryotes (Pabarcus and Casida 2002Go, 2005Go). The distribution of these two different proteins differs significantly according to the bacteria analyzed, as shown in table 1. Xanthomonas homologs were identified based on their similarity to the eukaryotic kynB (reciprocal BlastP searches confirm these genes as being similar to KFA-coding genes in mice and humans). On the other hand, Flavobacteriales homologs were identified based on their similarity to the cyclase form (as identified by Kurnasov, Goral, et al. 2003Go; Kurnasov, Jablonski, et al. 2003Go). However, some species of this order also carry the eukaryotic esterase form (table 1). All genes identified as putative KFA possess the conserved catalytic residues characteristic of both esterase and cyclase-like genes (supplementary fig. S2, Supplementary Material online). Therefore, the functional characterization of these two homolog forms is still necessary to understand which protein is fulfilling the role of kynurenine formamidase in these organisms.

Evolution of the Kynurenine Pathway in Bacteria
The discovery of a typical eukaryotic metabolic pathway in bacterial species is uncommon and evokes an evolutionary explanation that accounts for such a scenario. In order to trace the evolution of each gene belonging to the kynurenine pathway, phylogenetic reconstructions were conducted (fig. 3 and supplementary fig. S3, Supplementary Material online). The general pattern observed reinforces the relationship between the eukaryotic homologs and the genes present in Xanthomonas and Flavobacteriales. In fact, the branching of both groups with Eukaryotes and the high level of sequence similarity to eukaryotic genes suggest a potential foreign origin for such genes.


Figure 3
View larger version (19K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
FIG. 3.— ML tree of kmo (A) and kynU (B) genes. Branch lengths are proportional to the number of amino acid substitutions per site. Numbers at the nodes represent the percentage of bootstrap support (upper values for the ML tree and lower values for the NJ tree. Only values >60% are shown).

 
As reflected by the profile of presence and absence (fig. 2), kmo (fig. 3A) and haaO (supplementary fig. S3, Supplementary Material online) genes are those that present the most restricted phyletic distribution. Both genes are present in eukaryotic species, mainly those belonging to Unikonts (comprising Metazoa and Fungi; for a review, see Cavalier-Smith 2002Go), as well as in Xanthomonas and Flavobacteriales. Few bacterial species (mostly proteobacterial) possess the kmo gene and a few more, including species belonging to Firmicutes and Proteobacteria, a copy of the haaO. However, these bacterial haaO homologs form a monophyletic group, independent of those found in Xanthomonas, Flavobacteriales, and Eukaryotes.

The kynU gene is present in a wider range of bacterial species of the Firmicutes and Proteobacteria (fig. 3B). However, the homologs of Xanthomonas and Flavobacteriales branch with Eukaryotes and only a few other bacteria species (those bearing kmo and haaO genes, such as Myxococcus xanthus, Salinispora ssp., and Acidobacteria bacterium). The separation of these homologs from the main group of Bacteria is supported by high bootstrap values, both for ML and NJ trees.

Phylogenetic reconstructions conducted for kynA gene (data not shown) place Flavobacteriales with Eukarya (although with low bootstrap support) and Xanthomonadales branches with gram-positive species (Actinobacteria and Firmicutes); Proteobacteria form a distinct and monophyletic group. Interestingly, the second copy of kynA gene from X. campestris clusters with Eukarya and Flavobacteriales, thus indicating an independent origin for this, probably recently acquired by an ancestor of both X. campestris pv. campestris strains (table 1).

As stated before, degradation of tryptophan in bacteria may be performed by the tryptophanase TnaA or by way of the anthranilate-degrading pathway, although production of quinolinic acid is mostly carried out by the aspartate pathway (nadA and nadB genes). Both nadA and nadB genes are absent in eukaryotes, with a few exceptions in the Viridiplantae kingdom. On the other hand, the production of quinolinate is fulfilled by the kynurenine pathway, similar to what occurs in Xanthomonas and Flavobacteriales. The question is, how did this profile arise during the evolution of these bacteria?

The simplest and most parsimonious explanation for the widespread presence of the aspartate route and the scattered distribution of the anthranilate pathway would be an ancient bacterial origin for genes belonging to both pathways. Genes related to the anthranilate pathway would have been lost in several bacterial groups, probably due to functional redundancy in other systems for amino acid degradation. However, another scenario should be evoked to explain the presence and distribution of an entire kynurenine pathway in Xanthomonas and Flavobacteriales. The odd phyletic distribution of kyn genes (especially kynU, kmo, and haaO) in Xanthomonadales and Bacteroidetes, along with the results of phylogenetic reconstructions and similarity searches, points to the lateral acquisition of kynurenine pathway genes in these groups. The phylogenetic placement of Xanthomonas, Flavobacteriales, and Eukaryotes, branching together even when homologs are present in other Bacteria (as in the case of kynU and haaO genes), supports the differential origin of such genes in Xanthomonas and Flavobacteriales. The polarity of LGT is difficult to assess because of uncertain rooting, but data indicate that these genes, necessary for the complete kynurenine pathway, would have been acquired later by gene transfer to an ancestor of both. They possibly originated from Eukaryotes, where homologs are widely distributed in Unikonts and then recruited to NAD biosynthesis.

Moreover, the apparent phylogenetic proximity of the homologs between Xanthomonas and Flavobacteriales indicates that transfer might have occurred to one of these bacteria, then spreading in one of the groups, before being transferred to the other. In fact, further exchange of genes between Bacteroidetes and Xanthomonadales may have also occurred. In a recent study, certain genes of the arginine biosynthetic pathway (argF, argG, and argH) were described as having been laterally transferred among Xanthomonadales, Bacteroidetes, and Eukaryotes (Lima and Menck 2008Go).

An alternative explanation, a bacterial origin for these genes followed by a gene transfer to an ancestor of the Unikonts, cannot be excluded. However, this is unlikely, as it would presuppose a massive gene loss in most of the bacterial groups. The loss would be particularly recent in Xylellas and Bacteroidales, thus arguing against the bacterial origin of such genes.

The loss of nadA and nadB genes in Xanthomonas and Flavobacteriales due to redundancy in the kynurenine pathway may be an example of xenologous pathway replacement. In fact, the success and maintenance of one quinolinate biosynthetic pathway over another in a given organism is probably governed by complex factors, and it is interesting to discuss the ecological reasons for this. In any scenario (of gene transfer from a bacterial or eukaryotic donor), the introgression of kynurenine genes may have coexisted with the aspartate pathway, which was then lost, due to duplicated function.

In both groups, there are free living and pathogenic species. Xylellas are obligatory plant parasites, living in an extremely dilute nutritional environment and present an intermediate insect vector, where it lives in the gut (Hopkins 1989Go; Monteiro-Vitorello et al. 2005Go; Moreira et al. 2005Go). Similarly, Bacteroidales are obligatory anaerobic parasites or symbionts of the mammalian intestinal tract (Xu et al. 2003Go, 2007Go). The presence of sufficient aspartate and/or tryptophan in the milieu of ancestral bacteria, during a specific period in evolution, may have dictated the fate of those pathways that were retained in these bacteria. However, these bacterial groups also present genes that potentially encode key proteins for the synthesis of aspartate from the citrate cycle (aspartate aminotransferase) or tryptophan (tryptophan synthase). Another explanation would be the aerobic or anaerobic status of each organism, as the kynurenine pathway requires three molecules of oxygen (in the reactions catalyzed by TDO, KMO, and HAD), whereas no oxygen is required for the aspartate pathway. Bacteroidales species are strict anaerobics and Xylella ssp., although growing in aerobic conditions in in vitro cultures, possess an unusual and least energy-efficient aerobic respiratory system, indicating that this is not the prevalent type of energy metabolism in these organisms (Bhattacharyya et al. 2002Go).

In any case, it is difficult to account for the benefit of such a change in the metabolic network. The aspartate pathway seems to be the most primitive. In an interesting discussion, Cleaves and Miller (2001)Go stated that NadA and NadB functions were even dispensable for the interconversion of aspartate and dihydroxyacetone phosphate into quinolinate in a hypothetically primitive world. Therefore, this pathway might have arisen very early during evolution, independent of any enzymatic chemical reaction. On the other hand, the kynurenine pathway is a multistep process with a rate-limiting substrate (as tryptophan is one of the rarest amino acids). The successful pathway replacement in both bacterial groups reveals the degree of complexity during evolution.


    Conclusion
 TOP
 Abstract
 Background
 Methods
 Results and Discussion
 Conclusion
 Supplementary Material
 Acknowledgements
 References
 
Lateral gene transfer is well accepted as changing bacterial gene content in prokaryotes by the increment of genes and/or pathways for secondary metabolism, thus likely increasing fitness (such as antibiotic resistance, pathogenesis, etc.). The transfer of genes involved in primary metabolism is considered to be a less frequent process (Nakamura et al. 2004Go; Pal et al. 2005Go; Merkl 2006Go; Homma et al. 2007Go), although these events have been reported (Boucher et al. 2003Go; Omelchenko et al. 2003Go). The results described in this work provide evidence that transfer may have occurred for pathways involved in the biosynthesis of quinolinate, an essential precursor of NAD, together with de novo acquisition of genes involved in the kynurenine pathway. In this case, the transfer of such a primary metabolic pathway could have triggered the replacement of a previously well-established pathway, which further supports the high relevance of these lateral gene transfer processes during the evolution of bacterial genomes.


    Supplementary Material
 TOP
 Abstract
 Background
 Methods
 Results and Discussion
 Conclusion
 Supplementary Material
 Acknowledgements
 References
 
Supplementary table S1 (List of seed sequences used to similarity search analysis) and table S2 (Genome's GenBank accession numbers); supplementary figure S1 (Chromosomal arrangement of the kynurenine pathway genes), figure S2 (Alignment showing the conserved residues in kynB orthologs), and figure S3 (phylogenetic tree of haaO gene); and protein sequences (in fasta format) for all ORFs referred in table 1 are available at Molecular Biology and Evolutiononline (http://www.mbe.oxfordjournals.org/).


    Acknowledgements
 TOP
 Abstract
 Background
 Methods
 Results and Discussion
 Conclusion
 Supplementary Material
 Acknowledgements
 References
 
This work was supported by Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP; São Paulo, SP, Brazil) and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq; Brasília, DF, Brazil). We also would like to thank the anonymous reviewers for their very interesting contributions to this manuscript.


    Footnotes
 
Hervé Philippe, Associate Editor


    References
 TOP
 Abstract
 Background
 Methods
 Results and Discussion
 Conclusion
 Supplementary Material
 Acknowledgements
 References
 

    Begley TP, Kinsland C, Mehl RA, Osterman A, Dorrestein P. The biosynthesis of nicotinamide adenine dinucleotides in bacteria. Vitam Horm (2001) 61:103–119.[Web of Science][Medline]

    Bhattacharyya A, Stilwagen S, Reznik G, et al, (22 co-authors). Draft sequencing and comparative genomics of Xylella fastidiosa strains reveal novel biological insights. Genome Res (2002) 12:1556–1563.[Abstract/Free Full Text]

    Boucher Y, Douady CJ, Papke RT, Walsh DA, Boudreau ME, Nesbo CL, Case RJ, Doolittle WF. Lateral gene transfer and the origins of prokaryotic groups. Annu Rev Genet (2003) 37:283–328.[CrossRef][Web of Science][Medline]

    Brown AT, Wagner C. Regulation of enzymes involved in the conversion of tryptophan to nicotinamide adenine dinucleotide in a colorless strain of Xanthomonas pruni. J Bacteriol (1970) 101:456–463.[Abstract/Free Full Text]

    Brown D, Hitchcock MJ, Katz E. Purification and characterization of kynurenine formamidase activities from Streptomyces parvulus. Can J Microbiol (1986) 32:465–472.[Web of Science][Medline]

    Cavalier-Smith T. The phagotrophic origin of eukaryotes and phylogenetic classification of Protozoa. Int J Syst Evol Microbiol (2002) 52:297–354.[Abstract]

    Cleaves HJ, Miller SL. The nicotinamide biosynthetic pathway is a by-product of the RNA world. J Mol Evol (2001) 52:73–77.[Web of Science][Medline]

    D'Argenio DA, Calfee MW, Rainey PB, Pesci EC. Autolysis and autoaggregation in Pseudomonas aeruginosa colony morphology mutants. J Bacteriol (2002) 184:6481–6489.[Abstract/Free Full Text]

    Davis D, Henderson LM, Powell D. The niacin–tryptophan relationship in the metabolism of Xanthomonas pruni. J Biol Chem (1951) 189:543–549.[Free Full Text]

    Farrow JM 3rd, Pesci EC. Two distinct pathways supply anthranilate as a precursor of the Pseudomonas quinolone signal. J Bacteriol (2007) 189:3425–3433.[Abstract/Free Full Text]

    Felsenstein J. PHYLIP (Phylogeny Inference Package) version 3.6. In: Distributed by the author (2005) Seattle (WA): Department of Genome Sciences, University of Washington.

    Fukuoka S, Ishiguro K, Yanagihara K, Tanabe A, Egashira Y, Sanada H, Shibata K. Identification and expression of a cDNA encoding human alpha-amino-beta-carboxymuconate-epsilon-semialdehyde decarboxylase (ACMSD). A key enzyme for the tryptophan-niacine pathway and "quinolinate hypothesis". J Biol Chem (2002) 277:35162–35167.[Abstract/Free Full Text]

    Gasteiger E, Gattiker A, Hoogland C, Ivanyi I, Appel RD, Bairoch A. ExPASy: the proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res (2003) 31:3784–3788.[Abstract/Free Full Text]

    Homma K, Fukuchi S, Nakamura Y, Gojobori T, Nishikawa K. Gene cluster analysis method identifies horizontally transferred genes with high reliability and indicates that they provide the main mechanism of operon gain in 8 species of gamma-Proteobacteria. Mol Biol Evol (2007) 24:805–813.[Abstract/Free Full Text]

    Hopkins DL. Xylella-fastidiosa—xylem-limited bacterial pathogen of plants. Annu Rev Phytopathol (1989) 27:271–290.[CrossRef][Web of Science]

    Kanehisa M, Araki M, Goto S, et al, (11 co-authors). KEGG for linking genomes to life and the environment. Nucleic Acids Res (2008) 36:D480–D484.[Abstract/Free Full Text]

    Katoh A, Hashimoto T. Molecular biology of pyridine nucleotide and nicotine biosynthesis. Front Biosci (2004) 9:1577–1586.[Web of Science][Medline]

    Kurnasov O, Goral V, Colabroy K, Gerdes S, Anantha S, Osterman A, Begley TP. NAD biosynthesis: identification of the tryptophan to quinolinate pathway in bacteria. Chem Biol (2003) 10:1195–1204.[CrossRef][Web of Science][Medline]

    Kurnasov O, Jablonski L, Polanuyer B, Dorrestein P, Begley T, Osterman A. Aerobic tryptophan degradation pathway in bacteria: novel kynurenine formamidase. FEMS Microbiol Lett (2003) 227:219–227.[CrossRef][Web of Science][Medline]

    Larkin MA, Blackshields G, Brown NP, et al, (13 co-authors). Clustal W and Clustal X version 2.0. Bioinformatics (2007) 23:2947–2948.[Abstract/Free Full Text]

    Lima WC, Menck CF. Replacement of the arginine biosynthesis operon in Xanthomonadales by lateral gene transfer. J Mol Evol (2008) 66:266–275.[CrossRef][Web of Science][Medline]

    Lima WC, Paquola AC, Varani AM, Van Sluys MA, Menck CF. Laterally transferred genomic islands in Xanthomonadales related to pathogenicity and primary metabolism. FEMS Microbiol Lett (2008) 281:87–97.[CrossRef][Web of Science][Medline]

    Lima WC, Van Sluys MA, Menck CF. Non-gamma-proteobacteria gene islands contribute to the Xanthomonas genome. OMICS (2005) 9:160–172.[CrossRef][Web of Science][Medline]

    Mattevi A. A close look at NAD biosynthesis. Nat Struct Mol Biol (2006) 13:563–564.[CrossRef][Web of Science][Medline]

    Merkl R. A comparative categorization of protein function encoded in bacterial or archeal genomic islands. J Mol Evol (2006) 62:1–14.[CrossRef][Web of Science][Medline]

    Monteiro-Vitorello CB, de Oliveira MC, Zerillo MM, Varani AM, Civerolo E, Van Sluys MA. Xylella and Xanthomonas Mobil'omics. OMICS (2005) 9:146–159.[CrossRef][Web of Science][Medline]

    Moreira LM, De Souza RF, Digiampietri LA, Da Silva AC, Setubal JC. Comparative analyses of Xanthomonas and Xylella complete genomes. OMICS (2005) 9:43–76.[CrossRef][Web of Science][Medline]

    Nakamura Y, Itoh T, Matsuda H, Gojobori T. Biased biological functions of horizontally transferred genes in prokaryotic genomes. Nat Genet (2004) 36:760–766.[CrossRef][Web of Science][Medline]

    Nicholas KB, Nicholas HBJ, Deerfield DWI. GeneDoc: analysis and visualization of genetic variation. EMBNEW News (1997) 4:14.

    Omelchenko MV, Makarova KS, Wolf YI, Rogozin IB, Koonin EV. Evolution of mosaic operons by horizontal gene transfer and gene displacement in situ. Genome Biol (2003) 4:R55.[CrossRef][Medline]

    Osterman A, Overbeek R. Missing genes in metabolic pathways: a comparative genomics approach. Curr Opin Chem Biol (2003) 7:238–251.[CrossRef][Web of Science][Medline]

    Overbeek R, Begley T, Butler RM, et al, (30 co-authors). The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes. Nucleic Acids Res (2005) 33:5691–5702.[Abstract/Free Full Text]

    Pabarcus MK, Casida JE. Kynurenine formamidase: determination of primary structure and modeling-based prediction of tertiary structure and catalytic triad. Biochim Biophys Acta (2002) 1596:201–211.[CrossRef][Medline]

    Pabarcus MK, Casida JE. Cloning, expression, and catalytic triad of recombinant arylformamidase. Protein Expr Purif (2005) 44:39–44.[CrossRef][Web of Science][Medline]

    Page RD. TreeView: an application to display phylogenetic trees on personal computers. Comput Appl Biosci (1996) 12:357–358.[Free Full Text]

    Pal C, Papp B, Lercher MJ. Adaptive evolution of bacterial metabolic networks by horizontal gene transfer. Nat Genet (2005) 37:1372–1375.[CrossRef][Web of Science][Medline]

    Peterson JD, Umayam LA, Dickinson T, Hickey EK, White O. The comprehensive microbial resource. Nucleic Acids Res (2001) 29:123–125.[Abstract/Free Full Text]

    Pruitt KD, Tatusova T, Maglott DR. NCBI Reference Sequence (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res (2005) 33:D501–D504.[Abstract/Free Full Text]

    Rongvaux A, Andris F, Van Gool F, Leo O. Reconstructing eukaryotic NAD metabolism. Bioessays (2003) 25:683–690.[CrossRef][Web of Science][Medline]

    Salcher O, Lingens F. Metabolism of tryptophan by Pseudomonas aureofaciens and its relationship to pyrrolnitrin biosynthesis. J Gen Microbiol (1980) 121:465–471.[Abstract/Free Full Text]

    Schomburg I, Chang A, Ebeling C, Gremse M, Heldt C, Huhn G, Schomburg D. BRENDA, the enzyme database: updates and major new developments. Nucleic Acids Res (2004) 32:D431–D433.[Abstract/Free Full Text]

    Selengut JD, Haft DH, Davidsen T, Ganapathy A, Gwinn-Giglio M, Nelson WC, Richter AR, White O. TIGRFAMs and genome properties: tools for the assignment of molecular function and biological process in prokaryotic genomes. Nucleic Acids Res (2007) 35:D260–D264.[Abstract/Free Full Text]

    Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics (2006) 22:2688–2690.[Abstract/Free Full Text]

    Tamura K, Dudley J, Nei M, Kumar S. MEGA4: molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol (2007) 24:1596–1599.[Abstract/Free Full Text]

    Tremblay GC, Gottlieb JA, Knox WE. Induction by L-tryptophan and an analogue, alpha-methyl-DL-tryptophan, of the enzymes catabolizing L-tryptophan in Pseudomonas. J Bacteriol (1967) 93:168–176.[Abstract/Free Full Text]

    Vederas JC, Schleicher E, Tsai MD, Floss HG. Stereochemistry and mechanism of reactions catalyzed by tryptophanase Escherichia coli. J Biol Chem (1978) 253:5350–5354.[Abstract/Free Full Text]

    Wilson RG, Henderson LM. Tryptophan–niacin relationship in Xanthomonas pruni. J Bacteriol (1963) 85:221–229.[Abstract/Free Full Text]

    Xu J, Bjursell MK, Himrod J, Deng S, Carmichael LK, Chiang HC, Hooper LV, Gordon JI. A genomic view of the human-Bacteroides thetaiotaomicron symbiosis. Science (2003) 299:2074–2076.[Abstract/Free Full Text]

    Xu J, Mahowald MA, Ley RE, et al, (19 co-authors). Evolution of symbiotic bacteria in the distal human intestine. PLoS Biol (2007) 5:e156.[CrossRef][Medline]

Accepted for publication November 7, 2008.


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?



This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary Data
Right arrow All Versions of this Article:
26/2/399    most recent
msn261v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Lima, W. C.
Right arrow Articles by Menck, C. F.M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Lima, W. C.
Right arrow Articles by Menck, C. F.M.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?