Ancient origin of the tryptophan operon and the dynamics of evolutionary change

Metabolic engineering for improving anthranilate synthesis from glucose in Escherichia coli

BACKGROUND

Anthranilate is an aromatic amine used industrially as an intermediate for the synthesis of dyes, perfumes, pharmaceuticals and other classes of products. Chemical synthesis of anthranilate is an unsustainable process since it implies the use of nonrenewable benzene and the generation of toxic by-products. In Escherichia coli anthranilate is synthesized from chorismate by anthranilate synthase (TrpED) and then converted to phosphoribosyl anthranilate by anthranilate phosphoribosyl transferase to continue the tryptophan biosynthetic pathway. With the purpose of generating a microbial strain for anthranilate production from glucose, E. coli W3110 trpD9923, a mutant in the trpD gene that displays low anthranilate producing capacity, was characterized and modified using metabolic engineering strategies.

RESULTS

Sequencing of the trpED genes from E. coli W3110 trpD9923 revealed a nonsense mutation in the trpD gene, causing the loss of anthranilate phosphoribosyl transferase activity, but maintaining anthranilate synthase activity, thus causing anthranilate accumulation. The effects of expressing genes encoding a feedback inhibition resistant version of the enzyme 3-deoxy-D-arabino-heptulosonate-7-phosphate synthase (aroGfbr), transketolase (tktA), glucokinase (glk) and galactose permease (galP), as well as phosphoenolpyruvate:sugar phosphotransferase system (PTS) inactivation on anthranilate production capacity, were evaluated. In shake flask experiments with minimal medium, strains W3110 trpD9923 PTS- and W3110 trpD9923/pJLBaroGfbrtktA displayed the best production parameters, accumulating 0.70-0.75 g/L of anthranilate, with glucose-yields corresponding to 28-46% of the theoretical maximum. To study the effects of extending the growth phase on anthranilate production a fed-batch fermentation process was developed using complex medium, where strain W3110 trpD9923/pJLBaroGfbrtktA produced 14 g/L of anthranilate in 34 hours.

CONCLUSION

This work constitutes the first example of a microbial system for the environmentally-compatible synthesis of anthranilate generated by metabolic engineering. The results presented here, including the characterization of mutation in the trpD gene from strain W3110 trpD9923 and the development of a fermentation strategy, establish a step forward towards the future improvement of a sustainable process for anthranilate production. In addition, the present work provides very useful data regarding the positive and negative consequences of the evaluated metabolic engineering strategies.

Origin and evolution of metabolic pathways

The emergence and evolution of metabolic pathways represented a crucial step in molecular and cellular evolution. In fact, the exhaustion of the prebiotic supply of amino acids and other compounds that were likely present in the ancestral environment, imposed an important selective pressure, favoring those primordial heterotrophic cells which became capable of synthesizing those molecules. Thus, the emergence of metabolic pathways allowed primitive organisms to become increasingly less-dependent on exogenous sources of organic compounds. Comparative analyses of genes and genomes from organisms belonging to Archaea, Bacteria and Eukarya revealed that, during evolution, different forces and molecular mechanisms might have driven the shaping of genomes and the arisal of new metabolic abilities. Among these gene elongations, gene and operon duplications undoubtedly played a major role since they can lead to the (immediate) appearance of new genetic material that, in turn, might undergo evolutionary divergence giving rise to new genes coding for new metabolic abilities. Gene duplication has been invoked in the different schemes proposed to explain why and how the extant metabolic pathways have arisen and shaped. Both the analysis of completely sequenced genomes and directed evolution experiments strongly support one of them, i.e. the patchwork hypothesis, according to which metabolic pathways have been assembled through the recruitment of primitive enzymes that could react with a wide range of chemically related substrates. However, the analysis of the structure and organization of genes belonging to ancient metabolic pathways, such as histidine biosynthesis and nitrogen fixation, suggested that other different hypothesis, i.e. the retrograde hypothesis or the semi-enzymatic theory, may account for the arisal of some metabolic routes.

Purification and characterization of anthranilate synthase component I (TrpE) from Mycobacterium tuberculosis H37Rv

The emergence of multi-drug resistant (MDR) strains of Mycobacterium tuberculosis is the main reason why tuberculosis (TB) continues to be a major health problem worldwide. It is urgent to discover novel anti-mycobacterial agents based on new drug targets for the treatment of TB, especially MDR-TB. Tryptophan biosynthetic pathway, which is essential for the survival of M. tuberculosis and meanwhile absent in mammals, provides potential anti-TB drug targets. One of the promising drug targets in this pathway is anthranilate synthase component I (TrpE), whose role is to catalyze the conversion of chorismate to anthranilate using ammonia as amino source. In order to get a deep understanding of TrpE, a study on purification and characteristic identification of TrpE is required. In this work, the putative trpE gene of M. tuberculosis H37Rv was expressed as a fusion protein with a 6x His-tag on the N-terminal (His-TrpE) in Escherichia coli. The recombinant TrpE protein was successfully purified and then its enzymatic characteristics were analyzed. The native TrpE without His-tag was obtained by removal of the N-terminal fusion partner of His-TrpE using enterokinase. It was found that N-terminal fusion partner had little influence on TrpE catalytic activity. In addition, the key residues related to enzyme catalytic activity and that involved in l-tryptophan inhibition were predicted in the structure of M. tuberculosis H37Rv TrpE. These results would be beneficial to the designing of novel anti-TB drugs with high potency and selectivity.

Evolution of bacterial trp operons and their regulation

Survival and replication of most bacteria require the ability to synthesize the amino acid L-tryptophan whenever it is not available from the environment. In this article we describe the genes, operons, proteins, and reactions involved in tryptophan biosynthesis in bacteria, and the mechanisms they use in regulating tryptophan formation. We show that although the reactions of tryptophan biosynthesis are essentially identical, gene organization varies among species--from whole-pathway operons to completely dispersed genes. We also show that the regulatory mechanisms used for these genes vary greatly. We address the question--what are some potential advantages of the gene organization and regulation variation associated with this conserved, important pathway?

The hidden universal distribution of amino acid biosynthetic networks: a genomic perspective on their origins and evolution

A core of widely distributed network branches biosynthesizing at least 16 out of the 20 standard amino acids is predicted using comparative genomics.

Background

Twenty amino acids comprise the universal building blocks of proteins. However, their biosynthetic routes do not appear to be universal from an Escherichia coli-centric perspective. Nevertheless, it is necessary to understand their origin and evolution in a global context, that is, to include more 'model' species and alternative routes in order to do so. We use a comparative genomics approach to assess the origins and evolution of alternative amino acid biosynthetic network branches.

Results

By tracking the taxonomic distribution of amino acid biosynthetic enzymes, we predicted a core of widely distributed network branches biosynthesizing at least 16 out of the 20 standard amino acids, suggesting that this core occurred in ancient cells, before the separation of the three cellular domains of life. Additionally, we detail the distribution of two types of alternative branches to this core: analogs, enzymes that catalyze the same reaction (using the same metabolites) and belong to different superfamilies; and 'alternologs', herein defined as branches that, proceeding via different metabolites, converge to the same end product. We suggest that the origin of alternative branches is closely related to different environmental metabolite sources and life-styles among species.

Conclusion

The multi-organismal seed strategy employed in this work improves the precision of dating and determining evolutionary relationships among amino acid biosynthetic branches. This strategy could be extended to diverse metabolic routes and even other biological processes. Additionally, we introduce the concept of 'alternolog', which not only plays an important role in the relationships between structure and function in biological networks, but also, as shown here, has strong implications for their evolution, almost equal to paralogy and analogy.

Cohesion group approach for evolutionary analysis of TyrA, a protein family with wide-ranging substrate specificities

Many enzymes and other proteins are difficult subjects for bioinformatic analysis because they exhibit variant catalytic, structural, regulatory, and fusion mode features within a protein family whose sequences are not highly conserved. However, such features reflect dynamic and interesting scenarios of evolutionary importance. The value of experimental data obtained from individual organisms is instantly magnified to the extent that given features of the experimental organism can be projected upon related organisms. But how can one decide how far along the similarity scale it is reasonable to go before such inferences become doubtful? How can a credible picture of evolutionary events be deduced within the vertical trace of inheritance in combination with intervening events of lateral gene transfer (LGT)? We present a comprehensive analysis of a dehydrogenase protein family (TyrA) as a prototype example of how these goals can be accomplished through the use of cohesion group analysis. With this approach, the full collection of homologs is sorted into groups by a method that eliminates bias caused by an uneven representation of sequences from organisms whose phylogenetic spacing is not optimal. Each sufficiently populated cohesion group is phylogenetically coherent and defined by an overall congruence with a distinct section of the 16S rRNA gene tree. Exceptions that occasionally are found implicate a clearly defined LGT scenario whereby the recipient lineage is apparent and the donor lineage of the gene transferred is localized to those organisms that define the cohesion group. Systematic procedures to manage and organize otherwise overwhelming amounts of data are demonstrated.

Cohesion group approach for evolutionary analysis of TyrA, a protein family with wide-ranging substrate specificities

Many enzymes and other proteins are difficult subjects for bioinformatic analysis because they exhibit variant catalytic, structural, regulatory, and fusion mode features within a protein family whose sequences are not highly conserved. However, such features reflect dynamic and interesting scenarios of evolutionary importance. The value of experimental data obtained from individual organisms is instantly magnified to the extent that given features of the experimental organism can be projected upon related organisms. But how can one decide how far along the similarity scale it is reasonable to go before such inferences become doubtful? How can a credible picture of evolutionary events be deduced within the vertical trace of inheritance in combination with intervening events of lateral gene transfer (LGT)? We present a comprehensive analysis of a dehydrogenase protein family (TyrA) as a prototype example of how these goals can be accomplished through the use of cohesion group analysis. With this approach, the full collection of homologs is sorted into groups by a method that eliminates bias caused by an uneven representation of sequences from organisms whose phylogenetic spacing is not optimal. Each sufficiently populated cohesion group is phylogenetically coherent and defined by an overall congruence with a distinct section of the 16S rRNA gene tree. Exceptions that occasionally are found implicate a clearly defined LGT scenario whereby the recipient lineage is apparent and the donor lineage of the gene transferred is localized to those organisms that define the cohesion group. Systematic procedures to manage and organize otherwise overwhelming amounts of data are demonstrated.

Predicted transcription factor binding sites as predictors of operons in Escherichia coli and Streptomyces coelicolor

BACKGROUND

As a polycistronic transcriptional unit of one or more adjacent genes, operons play a key role in regulation and function in prokaryotic biology, and a better understanding of how they are constituted and controlled is needed. Recent efforts have attempted to predict operonic status in sequenced genomes using a variety of techniques and data sources. To date, non-homology based operon prediction strategies have mainly used predicted promoters and terminators present at the extremities of transcriptional unit as predictors, with reasonable success. However, transcription factor binding sites (TFBSs), typically found upstream of the first gene in an operon, have not yet been evaluated.

RESULTS

Here we apply a method originally developed for the prediction of TFBSs in Escherichia coli that minimises the need for prior knowledge and tests its ability to predict operons in E. coli and the 'more complex', pharmaceutically important, Streptomyces coelicolor. We demonstrate that through building genome specific TFBS position-specific-weight-matrices (PSWMs) it is possible to predict operons in E. coli and S. coelicolor with 83% and 93% accuracy respectively, using only TFBS as delimiters of operons. Additionally, the 'palindromicity' of TFBS footprint data of E. coli is characterised.

CONCLUSION

TFBS are proposed as novel independent features for use in prokaryotic operon prediction (whether alone or as part of a set of features) given their efficacy as operon predictors in E. coli and S. coelicolor. We also show that TFBS footprint data in E. coli generally contains inverted repeats with significantly (p < 0.05) greater palindromicity than random sequences. Consequently, the palindromicity of putative TFBSs predicted can also enhance operon predictions.

TrpY regulation of trpB2 transcription in Methanothermobacter thermautotrophicus

TrpY binds specifically to TRP box sequences upstream of trpB2, but the repression of trpB2 transcription requires additional TrpY assembly that is stimulated by but not dependent on the presence of tryptophan. Inhibitory complex formation is prevented by insertions within the regulatory region and by a G149R substitution in TrpY, even though TrpY(G149R) retains both TRP box DNA- and tryptophan-binding abilities.

The tryptophan pathway genes of the Sargasso Sea metagenome: new operon structures and the prevalence of non-operon organization

BACKGROUND

The enormous database of microbial DNA generated from the Sargasso Sea metagenome provides a unique opportunity to locate genes participating in different biosynthetic pathways and to attempt to understand the relationship and evolution of those genes. In this article, an analysis of the Sargasso Sea metagenome is made with respect to the seven genes of the tryptophan pathway.

RESULTS

At least 5% of all the genes that are related to amino acid biosynthesis are tryptophan (trp) genes. Many contigs and scaffolds contain whole or split operons that are similar to previously analyzed trp gene organizations. Only two scaffolds discovered in this analysis possess a different operon organization of tryptophan pathway genes than those previously known. Many marine organisms lack an operon-type organization of these genes or have mini-operons containing only two trp genes. In addition, the trpB genes from this search reveal that the dichotomous division between trpB_1 and trpB_2 also occurs in organisms from the Sargasso Sea. One cluster was found to contain trpB sequences that were closely related to each other but distinct from most known trpB sequences.

CONCLUSION

The data show that trp genes are widely dispersed within this metagenome. The novel organization of these genes and an unusual group of trpB_1 sequences that were found among some of these Sargasso Sea bacteria indicate that there is much to be discovered about both the reason for certain gene orders and the regulation of tryptophan biosynthesis in marine bacteria.

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

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.

Comparison of tryptophan biosynthetic operon regulation in different Gram-positive bacterial species

The tryptophan biosynthetic operon has been widely used as a model system for studying transcription regulation. In Bacillus subtilis, the trp operon is primarily regulated by a tryptophan-activated RNA-binding protein, TRAP. Here we show that in many other Gram-positive species the trp operon is regulated differently, by tRNA(Trp) sensing by the RNA-based T-box mechanism, with T-boxes arranged in tandem. Our analyses reveal an apparent relationship between trp operon organization and the specific regulatory mechanism(s) used.

Spontaneous trpY mutants and mutational analysis of the TrpY archaeal transcription regulator

Over 90% of Methanothermobacter thermautotrophicus mutants isolated as spontaneously resistant to 5-methyl tryptophan had mutations in trpY. Most were single-base-pair substitutions that identified separate DNA- and tryptophan-binding regions in TrpY. In vivo and in vitro studies revealed that DNA binding was sufficient for TrpY repression of trpY transcription but that TrpY must bind DNA and tryptophan to assemble a complex that represses trpEGCFBAD.

Symbiosis as an adaptive process and source of phenotypic complexity

Genomics has revealed that inheritance systems of separate species are often not well segregated: genes and capabilities that evolve in one lineage are often stably acquired by another lineage. Although direct gene transfer between species has occurred at some level in all major groups, it appears to be far more frequent in prokaryotes than in multicellular eukaryotes. An alternative to incorporating novel genes into a recipient genome is acquiring a stable, possibly heritable, symbiotic association and thus enjoying benefits of complementary metabolic capabilities. These kinds of symbioses have arisen frequently in animals; for example, many insect groups have diversified on the basis of symbiotic associations acquired early in their evolutionary histories. The resulting associations are highly complex, often involving specialized cell types and organs, developmental mechanisms that ensure transfer of symbionts between generations, and mechanisms for controlling symbiont proliferation and location. The genomes of long-term obligate symbionts often undergo irreversible gene loss and deterioration even as hosts evolve dependence on them. In some cases, animal genomes may have acquired genes from symbionts, mirroring the gene uptake from mitochondrial and plastid genomes. Multiple symbionts often coexist in the same host, resulting in coadaptation among several phylogenetically distant genomes.

Environmental selection of the feed-forward loop circuit in gene-regulation networks

Gene-regulation networks contain recurring elementary circuits termed network motifs. It is of interest to understand under which environmental conditions each motif might be selected. To address this, we study one of the most significant network motifs, a three-gene circuit called the coherent feed-forward loop (FFL). The FFL has been demonstrated theoretically and experimentally to perform a basic information-processing function: it shows a delay following ON steps of an input inducer, but not after OFF steps. Here, we ask under what environmental conditions might the FFL be selected over simpler gene circuits, based on this function. We employ a theoretical cost-benefit analysis for the selection of gene circuits in a given environment. We find conditions that the environment must satisfy in order for the FFL to be selected over simpler circuits: the FFL is selected in environments where the distribution of the input pulse duration is sufficiently broad and contains both long and short pulses. Optimal values of the biochemical parameters of the FFL circuit are determined as a function of the environment such that the delay in the FFL blocks deleterious short pulses of induction. This approach can be generally used to study the evolutionary selection of other network motifs.

Modelling the evolution of the archeal tryptophan synthase

Background

Microorganisms and plants are able to produce tryptophan. Enzymes catalysing the last seven steps of tryptophan biosynthesis are encoded in the canonical trp operon. Among the trp genes are most frequently trpA and trpB, which code for the alpha and beta subunit of tryptophan synthase. In several prokaryotic genomes, two variants of trpB (named trpB1 or trpB2) occur in different combinations. The evolutionary history of these trpB genes is under debate.

Results

In order to study the evolution of trp genes, completely sequenced archeal and bacterial genomes containing trpB were analysed. Phylogenetic trees indicated that TrpB sequences constitute four distinct groups; their composition is in agreement with the location of respective genes. The first group consisted exclusively of trpB1 genes most of which belonged to trp operons. Groups two to four contained trpB2 genes. The largest group (trpB2_o) contained trpB2 genes all located outside of operons. Most of these genes originated from species possessing an operon-based trpB1 in addition. Groups three and four pertain to trpB2 genes of those genomes containing exclusively one or two trpB2 genes, but no trpB1. One group (trpB2_i) consisted of trpB2 genes located inside, the other (trpB2_a) of trpB2 genes located outside the trp operon. TrpA and TrpB form a heterodimer and cooperate biochemically. In order to characterise trpB variants and stages of TrpA/TrpB cooperation in silico, several approaches were combined. Phylogenetic trees were constructed for all trp genes; their structure was assessed via bootstrapping. Alternative models of trpB evolution were evaluated with parsimony arguments. The four groups of trpB variants were correlated with archeal speciation. Several stages of TrpA/TrpB cooperation were identified and trpB variants were characterised. Most plausibly, trpB2 represents the predecessor of the modern trpB gene, and trpB1 evolved in an ancestral bacterium.

Conclusion

In archeal genomes, several stages of trpB evolution, TrpA/TrpB cooperation, and operon formation can be observed. Thus, archeal trp genes may serve as a model system for studying the evolution of protein-protein interactions and operon formation.

On the origin and evolution of biosynthetic pathways: integrating microarray data with structure and organization of the Common Pathway genes

BACKGROUND

The lysine, threonine, and methionine biosynthetic pathways share the three initial enzymatic steps, which are referred to as the Common Pathway (CP). In Escherichia coli three different aspartokinases (AKI, AKII, AKIII, the products of thrA, metL and lysC, respectively) can perform the first step of the CP. Moreover, two of them (AKI and AKII) are bifunctional, carrying also homoserine dehydrogenasic activity (hom product). The second step of the CP is catalyzed by a single aspartate semialdehyde dehydrogenase (ASDH, the product of asd). Thus, in the CP of E. coli while a single copy of ASDH performs the same reaction for three different metabolic routes, three different AKs perfom a unique step. Why and how such a situation did emerge and maintain? How is it correlated to the different regulatory mechanisms acting on these genes? The aim of this work was to trace the evolutionary pathway leading to the extant scenario in proteobacteria.

RESULTS

The analysis of the structure, organization, phylogeny, and distribution of ask and hom genes revealed that the presence of multiple copies of these genes and their fusion events are restricted to the gamma-subdivision of proteobacteria. This allowed us to depict a model to explain the evolution of ask and hom according to which the fused genes are the outcome of a cascade of gene duplication and fusion events that can be traced in the ancestor of gamma-proteobacteria. Moreover, the appearance of fused genes paralleled the assembly of operons of different sizes, suggesting a strong correlation between the structure and organization of these genes. A statistic analysis of microarray data retrieved from experiments carried out on E. coli and Pseudomonas aeruginosa was also performed.

CONCLUSION

The integration of data concerning gene structure, organization, phylogeny, distribution, and microarray experiments allowed us to depict a model for the evolution of ask and hom genes in proteobacteria and to suggest a biological significance for the extant scenario.

Evolution of proteins and proteomes: a phylogenetics approach

The study of evolutionary relationships among protein sequences was one of the first applications of bioinformatics. Since then, and accompanying the wealth of biological data produced by genome sequencing and other high-throughput techniques, the use of bioinformatics in general and phylogenetics in particular has been gaining ground in the study of protein and proteome evolution. Nowadays, the use of phylogenetics is instrumental not only to infer the evolutionary relationships among species and their genome sequences, but also to reconstruct ancestral states of proteins and proteomes and hence trace the paths followed by evolution. Here I survey recent progress in the elucidation of mechanisms of protein and proteome evolution in which phylogenetics has played a determinant role.

Complete genome sequence of the marine, chemolithoautotrophic, ammonia-oxidizing bacterium Nitrosococcus oceani ATCC 19707

The gammaproteobacterium Nitrosococcus oceani (ATCC 19707) is a gram-negative obligate chemolithoautotroph capable of extracting energy and reducing power from the oxidation of ammonia to nitrite. Sequencing and annotation of the genome revealed a single circular chromosome (3,481,691 bp; G+C content of 50.4%) and a plasmid (40,420 bp) that contain 3,052 and 41 candidate protein-encoding genes, respectively. The genes encoding proteins necessary for the function of known modes of lithotrophy and autotrophy were identified. Contrary to betaproteobacterial nitrifier genomes, the N. oceani genome contained two complete rrn operons. In contrast, only one copy of the genes needed to synthesize functional ammonia monooxygenase and hydroxylamine oxidoreductase, as well as the proteins that relay the extracted electrons to a terminal electron acceptor, were identified. The N. oceani genome contained genes for 13 complete two-component systems. The genome also contained all the genes needed to reconstruct complete central pathways, the tricarboxylic acid cycle, and the Embden-Meyerhof-Parnass and pentose phosphate pathways. The N. oceani genome contains the genes required to store and utilize energy from glycogen inclusion bodies and sucrose. Polyphosphate and pyrophosphate appear to be integrated in this bacterium's energy metabolism, stress tolerance, and ability to assimilate carbon via gluconeogenesis. One set of genes for type I ribulose-1,5-bisphosphate carboxylase/oxygenase was identified, while genes necessary for methanotrophy and for carboxysome formation were not identified. The N. oceani genome contains two copies each of the genes or operons necessary to assemble functional complexes I and IV as well as ATP synthase (one H(+)-dependent F(0)F(1) type, one Na(+)-dependent V type).

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

Methanococcus maripaludis is a strictly anaerobic, methane-producing archaeon and facultative autotroph capable of biosynthesizing all the amino acids and vitamins required for growth. In this work, the novel 6-deoxy-5-ketofructose-1-phosphate (DKFP) pathway for the biosynthesis of aromatic amino acids (AroAAs) and p-aminobenzoic acid (PABA) was demonstrated in M. maripaludis. Moreover, PABA was shown to be derived from an early intermediate in AroAA biosynthesis and not from chorismate. Following metabolic labelling with [U-(13)C]-acetate, the expected enrichments for phenylalanine and arylamine derived from PABA were observed. DKFP pathway activity was reduced following growth with aryl acids, an alternative source of the AroAAs. Lastly, a deletion mutant of aroA', which encodes the first step in the DKFP pathway, required AroAAs and PABA for growth. Complementation of the mutants by an aroA' expression vector restored the wild-type phenotype. In contrast, a deletion of aroB', which encodes the second step in the DKFP pathway, did not require AroAAs or PABA for growth. Presumably, methanococci contain an alternative activity for this step. These results identify the initial reactions of a new pathway for the biosynthesis of PABA in methanococci.

Analysis of gene expression in operons of Streptomyces coelicolor

Analysis of the relative transcript levels of intra-operonic genes in Streptomyces coelicolor suggests significant levels of internal regulation.

Background

Recent studies have shown that microarray-derived gene-expression data are useful for operon prediction. However, it is apparent that genes within an operon do not conform to the simple notion that they have equal levels of expression.

Results

To investigate the relative transcript levels of intra-operonic genes, we have used a Z-score approach to normalize the expression levels of all genes within an operon to expression of the first gene of that operon. Here we demonstrate that there is a general downward trend in expression from the first to the last gene in Streptomyces coelicolor operons, in contrast to what we observe in Escherichia coli. Combining transcription-factor binding-site prediction with the identification of operonic genes that exhibited higher transcript levels than the first gene of the same operon enabled the discovery of putative internal promoters. The presence of transcription terminators and abundance of putative transcriptional control sequences in S. coelicolor operons are also described.

Conclusion

Here we have demonstrated a polarity of expression in operons of S. coelicolor not seen in E. coli, bringing caution to those that apply operon prediction strategies based on E. coli 'equal-expression' to divergent species. We speculate that this general difference in transcription behavior could reflect the contrasting lifestyles of the two organisms and, in the case of Streptomyces, might also be influenced by its high G+C content genome. Identification of putative internal promoters, previously thought to cause problems in operon prediction strategies, has also been enabled.

Rate and Polarity of Gene Fusion and Fission in Oryza sativa and Arabidopsis thaliana

Eukaryotic gene fusion and fission events are mechanistically more complicated than in prokaryotes, and their quantitative contributions to genome evolution are still poorly understood. We have identified all differentially composite or split genes in 2 fully sequenced plant genomes, Oryza sativa and Arabidopsis thaliana. Out of 10,172 orthologous gene pairs, 60 (0.6% of the total) revealed a verified fusion or fission event in either lineage after the divergence of O. sativa and A. thaliana. Polarizing these events by outgroup comparison revealed differences in the rate of gene fission but not of gene fusion in the rice and Arabidopsis lineages. Gene fission occurred at a higher rate than gene fusion in the O. sativa lineage and was furthermore more common in rice than in Arabidopsis. Nucleotide insertion bias has promoted gene fission in the O. sativa lineage, consistent with its generally longer nucleotide sequences than A. thaliana in selectively neutral regions, and with the abundance of transposable elements in rice. The divergence time of monocots and dicots (140-200 Myr) indicates that gene fusion/fission events occur at an average rate of 1x10(-11) to 2x10(-11) events per gene per year, approximately 100-fold slower than the average per site nuclear nucleotide substitution rate in these lineages. Gene fusion and fission are thus rare and slow processes in higher plant genomes; they should be of utility to address deeper evolutionary relationships among plants--and the relationship of plants to other eukaryotic lineages--where sequence-based phylogenies provide equivocal or conflicting results.

Optimal gene partition into operons correlates with gene functional order

Gene arrangement into operons varies between bacterial species. Genes in a given system can be on one operon in some organisms and on several operons in other organisms. Existing theories explain why genes that work together should be on the same operon, since this allows for advantageous lateral gene transfer and accurate stoichiometry. But what causes the frequent separation into multiple operons of co-regulated genes that act together in a pathway? Here we suggest that separation is due to benefits made possible by differential regulation of each operon. We present a simple mathematical model for the optimal distribution of genes into operons based on a balance of the cost of operons and the benefit of regulation that provides 'just-when-needed' temporal order. The analysis predicts that genes are arranged such that genes on the same operon do not skip functional steps in the pathway. This prediction is supported by genomic data from 137 bacterial genomes. Our work suggests that gene arrangement is not only the result of random historical drift, genome re-arrangement and gene transfer, but has elements that are solutions of an evolutionary optimization problem. Thus gene functional order may be inferred by analyzing the operon structure across different genomes.

Analytical study of the multiplicity of regulatory mechanisms in the tryptophan operon

In this paper we study the stability of a previously introduced model for the tryptophan operon regulatory pathway. For this, we make use of the second Lyapunov's method. The results obtained for the wild-type and for a couple ofin-silico mutant bacterial strains allow a deeper understanding of the multiplicity of regulatory mechanisms in this operon. In particular, we confirm that enzyme inhibition and transcription attenuation strengthen the system stability, the effect of transcription attenuation being much shorter than that of enzyme inhibition. Furthermore, the analysis here presented provides some insights about how enzyme inhibition affects the system stability.

Screening a wide host-range, waste-water metagenomic library in tryptophan auxotrophs of Rhizobium leguminosarum and of Escherichia coli reveals different classes of cloned trp genes

A metagenomic cosmid library was constructed, in which the insert DNA was derived from bacteria in a waste-water treatment plant and the vector was the wide host-range cosmid pLAFR3. The library was screened for clones that could correct defined tryptophan auxotrophs of the alpha-proteobacterium Rhizobium leguminosarum and of Escherichia coli. A total of 26 different cosmids that corrected at least one trp mutant in one or both of these species were obtained. Several cosmids corrected the auxotrophy of one or more R. leguminosarum trp mutants, but not the corresponding mutants in E. coli. Conversely, one cosmid corrected trpA, B, C, D and E mutants of E. coli but none of the trp mutants of R. leguminosarum. Two of the Trp+ cosmids were examined in more detail. One contained a trp operon that resembled that of the pathogen Chlamydophila caviae, containing the unusual kynU gene, which specifies kynureninase. The other, whose trp genes functioned in R. leguminosarum but not in E. coli, contained trpDCFBA in an operon that is likely co-transcribed with five other genes, most of which had no known link with tryptophan synthesis. The sequences of these TRP proteins, and the products of nine other genes encoded by this cosmid, failed to affiliate them with any known bacterial lineage. For one metagenomic cosmid, lac reporter fusions confirmed that its cloned trp genes were transcribed in R. leguminosarum, but not in E. coli. Thus, rhizobia, with their many sigma-factors, may be well-suited hosts for metagenomic libraries, cloned in wide host-range vectors.

Directed evolution of metabolic pathways

The modification of cellular metabolism is of biotechnological and commercial significance because naturally occurring metabolic pathways are the source of diverse compounds used in fields ranging from medicine to bioremediation. Directed evolution is the experimental improvement of biocatalysts or cellular properties through iterative genetic diversification and selection procedures. The creation of novel metabolic functions without disrupting the balanced intracellular pool of metabolites is the primary challenge of pathway manipulation. The introduction of coordinated changes across multiple genetic elements, in conjunction with functional selection, presents an integrated approach for the modification of metabolism with benign physiological consequences. Directed evolution formats take advantage of the dynamic structures of genomes and genomic sub-structures and their ability to evolve in multiple directions in response to external stimuli. The elucidation, design and application of genome-restructuring mechanisms are key elements in the directed evolution of cellular metabolic pathways.

Evolutionary plasticity of methionine biosynthesis

Methionine is an essential cellular constituent, the initiator of protein synthesis and a precursor in many metabolic activities, such as methylation and formylation. Here we investigate the genomic distribution of the methionine biosynthetic pathway and analyze its evolutionary history by reconstructing the phylogeny of its enzymatic components. We demonstrate the evolutionary complexity of methionine synthesis and describe the various mechanisms that have shaped this biosynthetic pathway: gene duplication, functional reassignment, lateral acquisition and gene loss. Lateral gene transfer within and between domains and gene recruitment have played an important role in the evolution of this pathway, especially in its first and third enzymatic steps--homoserine activation and homocysteine methylation. These analyses are also the basis of predictions regarding methionine synthesis in Archaea, where the pathway is yet to be characterized. This study illustrates how diverse molecular solutions can fulfill a conserved function in living beings.

The subsystems approach to genome annotation and its use in the project to annotate 1000 genomes

The release of the 1000th complete microbial genome will occur in the next two to three years. In anticipation of this milestone, the Fellowship for Interpretation of Genomes (FIG) launched the Project to Annotate 1000 Genomes. The project is built around the principle that the key to improved accuracy in high-throughput annotation technology is to have experts annotate single subsystems over the complete collection of genomes, rather than having an annotation expert attempt to annotate all of the genes in a single genome. Using the subsystems approach, all of the genes implementing the subsystem are analyzed by an expert in that subsystem. An annotation environment was created where populated subsystems are curated and projected to new genomes. A portable notion of a populated subsystem was defined, and tools developed for exchanging and curating these objects. Tools were also developed to resolve conflicts between populated subsystems. The SEED is the first annotation environment that supports this model of annotation. Here, we describe the subsystem approach, and offer the first release of our growing library of populated subsystems. The initial release of data includes 180 177 distinct proteins with 2133 distinct functional roles. This data comes from 173 subsystems and 383 different organisms.

New insights into regulation of the tryptophan biosynthetic operon in Gram-positive bacteria

The tryptophan operon of Bacillus subtilis serves as an excellent model for investigating transcription regulation in Gram-positive bacteria. In this article, we extend this knowledge by analyzing the predicted regulatory regions in the trp operons of other fully sequenced Gram-positive bacteria. Interestingly, it appears that in eight of the organisms examined, transcription of the trp operon appears to be regulated by tandem T-box elements. These regulatory elements have recently been described in the trp operons of two bacterial species. Single T-box elements are commonly found in Gram-positive bacteria in operons encoding aminoacyl tRNA synthetases and proteins performing other functions. Different regulatory mechanisms appear to be associated with variations of trp gene organization within the trp operon.

Regulation of tryptophan operon expression in the archaeon Methanothermobacter thermautotrophicus

Conserved trp genes encode enzymes that catalyze tryptophan biosynthesis in all three biological domains, and studies of their expression in Bacteria and eukaryotes have revealed a variety of different regulatory mechanisms. The results reported here provide the first detailed description of an archaeal trp gene regulatory system. We have established that the trpEGCFBAD operon in Methanothermobacter thermautotrophicus is transcribed divergently from a gene (designated trpY) that encodes a tryptophan-sensitive transcription regulator. TrpY binds to TRP box sequences (consensus, TGTACA) located in the overlapping promoter regions between trpY and trpE, inhibiting trpY transcription in the absence of tryptophan and both trpY and trpEGCFBAD transcription in the presence of tryptophan. TrpY apparently inhibits trpY transcription by blocking RNA polymerase access to the site of trpY transcription initiation and represses trpEGCFBAD transcription by preventing TATA box binding protein (TBP) binding to the TATA box sequence. Given that residue 2 (W2) is the only tryptophan in TrpY and in TrpY homologues in other Euryarchaea and that there is only one tryptophan codon in the entire trpEGCFBAD operon (trpB encodes W175), expression of the trp operon may also be regulated in vivo by the supply of charged tRNA(Trp) available to translate the second codon of the trpY mRNA.

The origin and evolution of operons: the piecewise building of the proteobacterial histidine operon

The structure and organization of 470 histidine biosynthetic genes from 47 different proteobacteria were combined with phylogenetic inference to investigate the mechanisms responsible for assembly of the his pathway and the origin of his operons. Data obtained in this work showed that a wide variety of different organization strategies of his gene arrays exist and that some his genes or entire his operons are likely to have been horizontally transferred between bacteria of the same or different proteobacterial branches. We propose a "piecewise" model for the origin and evolution of proteobacterial his operons, according to which the initially scattered his genes of the ancestor of proteobacteria coded for monofunctional enzymes (except possibly for hisD) and underwent a stepwise compacting process that reached its culmination in some gamma-proteobacteria. The initial step of operon buildup was the formation of the his "core," a cluster consisting of four genes (hisBHAF) whose products interconnect histidine biosynthesis to both de novo synthesis of purine metabolism and that occurred in the common ancestor of the alpha/beta/gamma branches, possibly after its separation from the epsilon one. The following step was the formation of three mini-operons (hisGDC, hisBHAF, hisIE) transcribed from independent promoters, that very likely occurred in the ancestor of the beta/gamma-branch, after its separation from the alpha one. Then the three mini-operons joined together to give a compact operon. In most gamma-proteobacteria the two fusions involving the gene pairs hisN-B and hisI-E occurred. Finally the gamma-proteobacterial his operon was horizontally transferred to other proteobacteria, such as Campylobacter jejuni. The biological significance of clustering of his genes is also discussed.

Regulation of transcription in a reduced bacterial genome: nutrient-provisioning genes of the obligate symbiont Buchnera aphidicola

Buchnera aphidicola, the obligate symbiont of aphids, has an extremely reduced genome, of which about 10% is devoted to the biosynthesis of essential amino acids needed by its hosts. Most regulatory genes for these pathways are absent, raising the question of whether and how transcription of these genes responds to the major shifts in dietary amino acid content encountered by aphids. Using full-genome microarrays for B. aphidicola of the host Schizaphis graminum, we examined transcriptome responses to changes in dietary amino acid content and then verified behavior of individual transcripts using quantitative reverse transcriptase PCR. The only gene showing a consistent and substantial (>twofold) response was metE, which underlies methionine biosynthesis and which is the only amino acid biosynthetic gene retaining its ancestral regulator (metR). In another aphid host, Acyrthosiphon pisum, B. aphidicola has no functional metR and shows no response in metE transcript levels to changes in amino acid concentrations. Thus, the only substantial transcriptional response involves the one gene for which an ancestral regulator is retained. This result parallels that from a previous study on heat stress, in which only the few genes retaining the global heat shock promoter showed responses in transcript abundance. The irreversible losses of transcriptional regulators constrain ability to alter gene expression in the context of environmental fluctuations affecting the symbiotic partners.

The TyrA family of aromatic-pathway dehydrogenases in phylogenetic context

BACKGROUND

The TyrA protein family includes members that catalyze two dehydrogenase reactions in distinct pathways leading to L-tyrosine and a third reaction that is not part of tyrosine biosynthesis. Family members share a catalytic core region of about 30 kDa, where inhibitors operate competitively by acting as substrate mimics. This protein family typifies many that are challenging for bioinformatic analysis because of relatively modest sequence conservation and small size.

RESULTS

Phylogenetic relationships of TyrA domains were evaluated in the context of combinatorial patterns of specificity for the two substrates, as well as the presence or absence of a variety of fusions. An interactive tool is provided for prediction of substrate specificity. Interactive alignments for a suite of catalytic-core TyrA domains of differing specificity are also provided to facilitate phylogenetic analysis. tyrA membership in apparent operons (or supraoperons) was examined, and patterns of conserved synteny in relationship to organismal positions on the 16S rRNA tree were ascertained for members of the domain Bacteria. A number of aromatic-pathway genes (hisHb, aroF, aroQ) have fused with tyrA, and it must be more than coincidental that the free-standing counterparts of all of the latter fused genes exhibit a distinct trace of syntenic association.

CONCLUSION

We propose that the ancestral TyrA dehydrogenase had broad specificity for both the cyclohexadienyl and pyridine nucleotide substrates. Indeed, TyrA proteins of this type persist today, but it is also common to find instances of narrowed substrate specificities, as well as of acquisition via gene fusion of additional catalytic domains or regulatory domains. In some clades a qualitative change associated with either narrowed substrate specificity or gene fusion has produced an evolutionary "jump" in the vertical genealogy of TyrA homologs. The evolutionary history of gene organizations that include tyrA can be deduced in genome assemblages of sufficiently close relatives, the most fruitful opportunities currently being in the Proteobacteria. The evolution of TyrA proteins within the broader context of how their regulation evolved and to what extent TyrA co-evolved with other genes as common members of aromatic-pathway regulons is now feasible as an emerging topic of ongoing inquiry.

A core catalytic domain of the TyrA protein family: arogenate dehydrogenase from Synechocystis

The TyrA protein family includes prephenate dehydrogenases, cyclohexadienyl dehydrogenases and TyrA(a)s (arogenate dehydrogenases). tyrA(a) from Synechocystis sp. PCC 6803, encoding a 30 kDa TyrA(a) protein, was cloned into an overexpression vector in Escherichia coli. TyrA(a) was then purified to apparent homogeneity and characterized. This protein is a model structure for a catalytic core domain in the TyrA superfamily, uncomplicated by allosteric or fused domains. Competitive inhibitors acting at the catalytic core of TyrA proteins are analogues of any accepted cyclohexadienyl substrate. The homodimeric enzyme was specific for L-arogenate (K(m)=331 microM) and NADP+ (K(m)=38 microM), being unable to substitute prephenate or NAD+ respectively. L-Tyrosine was a potent inhibitor of the enzyme (K(i)=70 microM). NADPH had no detectable ability to inhibit the reaction. Although the mechanism is probably steady-state random order, properties of 2',5'-ADP as an inhibitor suggest a high preference for L-arogenate binding first. Comparative enzymology established that both of the arogenate-pathway enzymes, prephenate aminotransferase and TyrA(a), were present in many diverse cyanobacteria and in a variety of eukaryotic red and green algae.

Variation and evolution of biomolecular systems: Searching for functional relevance

The availability of genome sequences and functional genomics data from multiple species enables us to compare the composition of biomolecular systems like biochemical pathways and protein complexes between species. Here, we review small- and large-scale, "genomics-based" approaches to biomolecular systems variation. In general, caution is required when comparing the results of bioinformatics analyses of genomes or of functional genomics data between species. Limitations to the sensitivity of sequence analysis tools and the noisy nature of genomics data tend to lead to systematic overestimates of the amount of variation. Nevertheless, the results from detailed manual analyses, and of large-scale analyses that filter out systematic biases, point to a large amount of variation in the composition of biomolecular systems. Such observations challenge our understanding of the function of the systems and their individual components and can potentially facilitate the identification and functional characterization of sub-systems within a system. Mapping the inter-species variation of complex biomolecular systems on a phylogenetic species tree allows one to reconstruct their evolution.

Selection for gene clustering by tandem duplication

In prokaryotic genomes, related genes are frequently clustered in operons and higher-order arrangements that reflect functional context. Organization emerges despite rearrangements that constantly shuffle gene and operon order. Evidence is presented that the tandem duplication of related genes acts as a driving evolutionary force in the origin and maintenance of clusters. Gene amplification can be viewed as a dynamic and reversible regulatory mechanism that facilitates adaptation to variable environments. Clustered genes confer selective benefits via their ability to be coamplified. During evolution, rearrangements that bring together related genes can be selected if they increase the fitness of the organism in which they reside. Similarly, the benefits of gene amplification can prevent the dispersal of existing clusters. Examples of frequent and spontaneous amplification of large genomic fragments are provided. The possibility is raised that tandem gene duplication works in concert with horizontal gene transfer as interrelated evolutionary forces for gene clustering.

The origin and evolution of eucaryal HIS7 genes: from metabolon to bifunctional proteins?

The fifth step of histidine biosynthesis is catalysed by an imidazole glycerol-phosphate (IGP) synthase. In Archaea and Bacteria, the active form of IGP synthase is a stable 1:1 dimeric complex constituted by a glutamine amidotransferase (GAT) and a cyclase, the products of hisH and hisF. In Eucarya, the two activities are associated with a single bifunctional polypeptide encoded by HIS7. In this work, we report a comparative analysis of the amino acid sequence of all the available HisH, HisF and HIS7 proteins, which allowed depicting a likely evolutionary pathway leading to the present-day bifunctional HIS7 genes. According to the model that we propose, the bifunctional HIS7 gene is the outcome of a gene fusion event between two independent ancestral cistrons encoding an amidotransferase and a cyclase, respectively. The phylogenetic distribution of the eucaryal HIS7 genes and the analysis of all the available prokaryotic counterparts (hisH and hisF) revealed the absence of such fusions in prokaryotes, suggesting that the fusion event very likely occurred in an early stage of eucaryal evolution and was fixed in the nucleated cells. The biological significance of this gene fusion is also discussed.

Dynamic influence of feedback enzyme inhibition and transcription attenuation on the tryptophan operon response to nutritional shifts

A mathematical model of the tryptophan operon is developed. This model considers all of the system known regulatory mechanisms: repression, transcription attenuation, and feedback enzyme inhibition. Special attention is paid to the estimation of all the model parameters from reported experimental data. The model equations are numerically solved. An analysis of these solutions reveals that transcription attenuation helps to speed up the operon response to nutritional shifts, while enzyme inhibition increases the operon stability.

Two biosynthetic pathways for aromatic amino acids in the archaeon Methanococcus maripaludis

Methanococcus maripaludis is a strictly anaerobic, methane-producing archaeon. Aromatic amino acids (AroAAs) are biosynthesized in this autotroph either by the de novo pathway, with chorismate as an intermediate, or by the incorporation of exogenous aryl acids via indolepyruvate oxidoreductase (IOR). In order to evaluate the roles of these pathways, the gene that encodes the third step in the de novo pathway, 3-dehydroquinate dehydratase (DHQ), was deleted. This mutant required all three AroAAs for growth, and no DHQ activity was detectible in cell extracts, compared to 6.0 +/- 0.2 mU mg(-1) in the wild-type extract. The growth requirement for the AroAAs could be fulfilled by the corresponding aryl acids phenylacetate, indoleacetate, and p-hydroxyphenylacetate. The specific incorporation of phenylacetate into phenylalanine by the IOR pathway was demonstrated in vivo by labeling with [1-(13)C]phenylacetate. M. maripaludis has two IOR homologs. A deletion mutant for one of these homologs contained 76, 74, and 42% lower activity for phenylpyruvate, p-hydoxyphenylpyruvate, and indolepyruvate oxidation, respectively, than the wild type. Growth of this mutant in minimal medium was inhibited by the aryl acids, but the AroAAs partially restored growth. Genetic complementation of the IOR mutant also restored much of the wild-type phenotype. Thus, aryl acids appear to regulate the expression or activity of the de novo pathway. The aryl acids did not significantly inhibit the activity of the biosynthetic enzymes chorismate mutase, prephenate dehydratase, and prephenate dehydrogenase in cell extracts, so the inhibition of growth was probably not due to an effect on these enzymes.

Inter-genomic displacement via lateral gene transfer of bacterial trp operons in an overall context of vertical genealogy

Background

The growing conviction that lateral gene transfer plays a significant role in prokaryote genealogy opens up a need for comprehensive evaluations of gene-enzyme systems on a case-by-case basis. Genes of tryptophan biosynthesis are frequently organized as whole-pathway operons, an attribute that is expected to facilitate multi-gene transfer in a single step. We have asked whether events of lateral gene transfer are sufficient to have obscured our ability to track the vertical genealogy that underpins tryptophan biosynthesis.

Results

In 47 complete-genome Bacteria, the genes encoding the seven catalytic domains that participate in primary tryptophan biosynthesis were distinguished from any paralogs or xenologs engaged in other specialized functions. A reliable list of orthologs with carefully ascertained functional roles has thus been assembled and should be valuable as an annotation resource. The protein domains associated with primary tryptophan biosynthesis were then concatenated, yielding single amino-acid sequence strings that represent the entire tryptophan pathway. Lateral gene transfer of several whole-pathway trp operons was demonstrated by use of phylogenetic analysis. Lateral gene transfer of partial-pathway trp operons was also shown, with newly recruited genes functioning either in primary biosynthesis (rarely) or specialized metabolism (more frequently).

Conclusions

(i) Concatenated tryptophan protein trees are congruent with 16S rRNA subtrees provided that the genomes represented are of sufficiently close phylogenetic spacing. There are currently seven tryptophan congruency groups in the Bacteria. Recognition of a succession of others can be expected in the near future, but ultimately these should coalesce to a single grouping that parallels the 16S rRNA tree (except for cases of lateral gene transfer). (ii) The vertical trace of evolution for tryptophan biosynthesis can be deduced. The daunting complexities engendered by paralogy, xenology, and idiosyncrasies of nomenclature at this point in time have necessitated an expert-assisted manual effort to achieve a correct analysis. Once recognized and sorted out, paralogy and xenology can be viewed as features that enrich evolutionary histories.