Ancient Origin of the Tryptophan Operon and the Dynamics of Evolutionary Change

The seven conserved enzymatic domains required for tryptophan (Trp) biosynthesis are encoded in seven genetic regions that are organized differently (whole-pathway operons, multiple partial-pathway operons, and dispersed genes) in prokaryotes. A comparative bioinformatics evaluation of the conservation and organization of the genes of Trp biosynthesis in prokaryotic operons should serve as an excellent model for assessing the feasibility of predicting the evolutionary histories of genes and operons associated with other biochemical pathways. These comparisons should provide a better understanding of possible explanations for differences in operon organization in different organisms at a genomics level. These analyses may also permit identification of some of the prevailing forces that dictated specific gene rearrangements during the course of evolution. Operons concerned with Trp biosynthesis in prokaryotes have been in a dynamic state of flux. Analysis of closely related organisms among the Bacteria at various phylogenetic nodes reveals many examples of operon scission, gene dispersal, gene fusion, gene scrambling, and gene loss from which the direction of evolutionary events can be deduced. Two milestone evolutionary events have been mapped to the 16S rRNA tree of Bacteria, one splitting the operon in two, and the other rejoining it by gene fusion. The Archaea, though less resolved due to a lesser genome representation, appear to exhibit more gene scrambling than the Bacteria. The trp operon appears to have been an ancient innovation; it was already present in the common ancestor of Bacteria and Archaea. Although the operon has been subjected, even in recent times, to dynamic changes in gene rearrangement, the ancestral gene order can be deduced with confidence. The evolutionary history of the genes of the pathway is discernible in rough outline as a vertical line of descent, with events of lateral gene transfer or paralogy enriching the analysis as interesting features that can be distinguished. As additional genomes are thoroughly analyzed, an increasingly refined resolution of the sequential evolutionary steps is clearly possible. These comparisons suggest that present-day trp operons that possess finely tuned regulatory features are under strong positive selection and are able to resist the disruptive evolutionary events that may be experienced by simpler, poorly regulated operons.     

Multiple signaling systems target a core set of transition metal homeostasis genes using similar binding motifs

Bacterial response to metals can require complex regulation. We report an overlapping regulation for copper and zinc resistance genes in the denitrifying bacterium, Pseudomonas stutzeri RCH2, by three two‐component regulatory proteins CopR1, CopR2 and CzcR. We conducted genome‐wide evaluations to identify gene targets of two paralogous regulators, CopR1 and CopR2, annotated for copper signaling, and compared the results with the gene targets for CzcR, implicated in zinc signaling. We discovered that the CopRs and CzcR have largely common targets, and crossregulate a core set of P. stutzeri copper and zinc responsive genes. We established that this crossregulation is enabled by a conserved binding motif in the upstream regulatory regions of the target genes. The crossregulation is physiologically relevant as these regulators synergistically and antagonistically target multicopper oxidases, metal efflux and sequestration systems. CopR1 and CopR2 upregulate two cop operons encoding copper tolerance genes, while all three regulators downregulate a putative copper chaperone, Psest_1595. CzcR also upregulated the oprD gene and the CzcIABC Zn²⁺ efflux system, while CopR1 and CopR2 downregulated these genes. Our study suggests that crossregulation of copper and zinc homeostasis can be advantageous, and in P. stutzeri this is enabled by shared binding motifs for multiple response regulators.     

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

The seven conserved enzymatic domains required for tryptophan (Trp) biosynthesis are encoded in seven genetic regions that are organized differently (whole-pathway operons, multiple partial-pathway operons, and dispersed genes) in prokaryotes. A comparative bioinformatics evaluation of the conservation and organization of the genes of Trp biosynthesis in prokaryotic operons should serve as an excellent model for assessing the feasibility of predicting the evolutionary histories of genes and operons associated with other biochemical pathways. These comparisons should provide a better understanding of possible explanations for differences in operon organization in different organisms at a genomics level. These analyses may also permit identification of some of the prevailing forces that dictated specific gene rearrangements during the course of evolution. Operons concerned with Trp biosynthesis in prokaryotes have been in a dynamic state of flux. Analysis of closely related organisms among the Bacteria at various phylogenetic nodes reveals many examples of operon scission, gene dispersal, gene fusion, gene scrambling, and gene loss from which the direction of evolutionary events can be deduced. Two milestone evolutionary events have been mapped to the 16S rRNA tree of Bacteria, one splitting the operon in two, and the other rejoining it by gene fusion. The Archaea, though less resolved due to a lesser genome representation, appear to exhibit more gene scrambling than the Bacteria. The trp operon appears to have been an ancient innovation; it was already present in the common ancestor of Bacteria and Archaea. Although the operon has been subjected, even in recent times, to dynamic changes in gene rearrangement, the ancestral gene order can be deduced with confidence. The evolutionary history of the genes of the pathway is discernible in rough outline as a vertical line of descent, with events of lateral gene transfer or paralogy enriching the analysis as interesting features that can be distinguished. As additional genomes are thoroughly analyzed, an increasingly refined resolution of the sequential evolutionary steps is clearly possible. These comparisons suggest that present-day trp operons that possess finely tuned regulatory features are under strong positive selection and are able to resist the disruptive evolutionary events that may be experienced by simpler, poorly regulated operons.     

Characterization of the CopR regulon of Lactococcus lactis IL1403

To identify components of the copper homeostatic mechanism of Lactococcus lactis, we employed two-dimensional gel electrophoresis to detect changes in the proteome in response to copper. Three proteins upregulated by copper were identified: glyoxylase I (YaiA), a nitroreductase (YtjD), and lactate oxidase (LctO). The promoter regions of these genes feature cop boxes of consensus TACAnnTGTA, which are the binding site of CopY-type copper-responsive repressors. A genome-wide search for cop boxes revealed 28 such sequence motifs. They were tested by electrophoretic mobility shift assays for the interaction with purified CopR, the CopY-type repressor of L. lactis. Seven of the cop boxes interacted with CopR in a copper-sensitive manner. They were present in the promoter region of five genes, lctO, ytjD, copB, ydiD, and yahC; and two polycistronic operons, yahCD-yaiAB and copRZA. Induction of these genes by copper was confirmed by real-time quantitative PCR. The copRZA operon encodes the CopR repressor of the regulon; a copper chaperone, CopZ; and a putative copper ATPase, CopA. When expressed in Escherichia coli, the copRZA operon conferred copper resistance, suggesting that it functions in copper export from the cytoplasm. Other member genes of the CopR regulon may similarly be involved in copper metabolism.     

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.     

The life-cycle of operons

Operons are a major feature of all prokaryotic genomes, but how and why operon structures vary is not well understood. To elucidate the life-cycle of operons, we compared gene order between Escherichia coli K12 and its relatives and identified the recently formed and destroyed operons in E. coli. This allowed us to determine how operons form, how they become closely spaced, and how they die. Our findings suggest that operon evolution may be driven by selection on gene expression patterns. First, both operon creation and operon destruction lead to large changes in gene expression patterns. For example, the removal of lysA and ruvA from ancestral operons that contained essential genes allowed their expression to respond to lysine levels and DNA damage, respectively. Second, some operons have undergone accelerated evolution, with multiple new genes being added during a brief period. Third, although genes within operons are usually closely spaced because of a neutral bias toward deletion and because of selection against large overlaps, genes in highly expressed operons tend to be widely spaced because of regulatory fine-tuning by intervening sequences. Although operon evolution may be adaptive, it need not be optimal: new operons often comprise functionally unrelated genes that were already in proximity before the operon formed.     

Evolution of mosaic operons by horizontal gene transfer and gene displacement <Emphasis Type="Italic">in situ</Emphasis>

Background

Shuffling and disruption of operons and horizontal gene transfer are major contributions to the new, dynamic view of prokaryotic evolution. Under the 'selfish operon' hypothesis, operons are viewed as mobile genetic entities that are constantly disseminated via horizontal gene transfer, although their retention could be favored by the advantage of coregulation of functionally linked genes. Here we apply comparative genomics and phylogenetic analysis to examine horizontal transfer of entire operons versus displacement of individual genes within operons by horizontally acquired orthologs and independent assembly of the same or similar operons from genes with different phylogenetic affinities.

Results

Since a substantial number of operons have been identified experimentally in only a few model bacteria, evolutionarily conserved gene strings were analyzed as surrogates of operons. The phylogenetic affinities within these predicted operons were assessed first by sequence similarity analysis and then by phylogenetic analysis, including statistical tests of tree topology. Numerous cases of apparent horizontal transfer of entire operons were detected. However, it was shown that apparent horizontal transfer of individual genes or arrays of genes within operons is not uncommon either and results in xenologous gene displacement in situ, that is, displacement of an ancestral gene by a horizontally transferred ortholog from a taxonomically distant organism without change of the local gene organization. On rarer occasions, operons might have evolved via independent assembly, in part from horizontally acquired genes.

Conclusions

The discovery of in situ gene displacement shows that combination of rampant horizontal gene transfer with selection for preservation of operon structure provides for events in prokaryotic evolution that, a priori, seem improbable. These findings also emphasize that not all aspects of operon evolution are selfish, with operon integrity maintained by purifying selection at the organism level.

     

Comparative genomics reveals the acquisition of mobile genetic elements by the plant growth-promoting Pantoea eucrina OB49 in polluted environments

Heavy metal-tolerant plant growth-promoting bacteria (PGPB) have gained popularity in bioremediation in recent years. A genome-assisted study of a heavy metal-tolerant PGPB Pantoea eucrina OB49 isolated from the rhizosphere of wheat grown on a heavy metal-contaminated site is presented. Comparative pan-genome analysis indicated that OB49 acquired heavy metal resistance genes through horizontal gene transfer. On contigs S10 and S12, OB49 has two arsRBCH operons that give arsenic resistance. On the S12 contig, an arsRBCH operon was discovered in conjunction with the merRTPCADE operon, which provides mercury resistance. P. eucrina OB49 may be involved in an ecological alternative for heavy metal remediation and growth promotion of wheat grown in metal-polluted soils. Our results suggested the detection of mobile genetic elements that harbour the ars operon and the fluoride resistance genes adjacent to the mer operon.     

Germline Expression Influences Operon Organization in the Caenorhabditis elegans Genome

Operons are found across multiple kingdoms and phyla, from prokaryotes to chordates. In the nematode Caenorhabditis elegans, the genome contains >1000 operons that compose ~15% of the protein-coding genes. However, determination of the force(s) promoting the origin and maintenance of operons in C. elegans has proved elusive. Compared to bacterial operons, genes within a C. elegans operon often show poor coexpression and only sometimes encode proteins with related functions. Using analysis of microarray and large-scale in situ hybridization data, we demonstrate that almost all operon-encoded genes are expressed in germline tissue. However, genes expressed during spermatogenesis are excluded from operons. Operons group together along chromosomes in local clusters that also contain monocistronic germline-expressed genes. Additionally, germline expression of genes in operons is largely independent of the molecular function of the encoded proteins. These analyses demonstrate that mechanisms governing germline gene expression influence operon origination and/or maintenance. Thus, gene expression in a specific tissue can have profound effects on the evolution of genome organization.     

Theczc operon ofAlcaligenes eutrophus CH34: from resistance mechanism to the removal of heavy metals

Summary The plasmid-borneczc operon ensures for resistance to Cd²⁺, Zn²⁺ and Co²⁺ ions through a tricomponent export pathway and is associated to various conjugative plasmids ofA. eutrophus strains isolated from metal-contaminated industrial areas. Theczc region of pMOL30 was reassessed especially for the segments located upstream and downstream the structural genesczc CBA. In cultures grown with high concentrations of heavy metals,czc-mediated efflux of cations is followed by a process of metal bioprecipitation. These observations led to the development of bioreactors designed for the removal of heavy metals from polluted effluents.     

Operons Are a Conserved Feature of Nematode Genomes

The organization of genes into operons, clusters of genes that are co-transcribed to produce polycistronic pre-mRNAs, is a trait found in a wide range of eukaryotic groups, including multiple animal phyla. Operons are present in the class Chromadorea, one of the two main nematode classes, but their distribution in the other class, the Enoplea, is not known. We have surveyed the genomes of Trichinella spiralis, Trichuris muris, and Romanomermis culicivorax and identified the first putative operons in members of the Enoplea. Consistent with the mechanism of polycistronic RNA resolution in other nematodes, the mRNAs produced by genes downstream of the first gene in the T. spiralis and T. muris operons are trans-spliced to spliced leader RNAs, and we are able to detect polycistronic RNAs derived from these operons. Importantly, a putative intercistronic region from one of these potential enoplean operons confers polycistronic processing activity when expressed as part of a chimeric operon in Caenorhabditis elegans. We find that T. spiralis genes located in operons have an increased likelihood of having operonic C. elegans homologs. However, operon structure in terms of synteny and gene content is not tightly conserved between the two taxa, consistent with models of operon evolution. We have nevertheless identified putative operons conserved between Enoplea and Chromadorea. Our data suggest that operons and “spliced leader” (SL) trans-splicing predate the radiation of the nematode phylum, an inference which is supported by the phylogenetic profile of proteins known to be involved in nematode SL trans-splicing.     

Divergence and redundancy of 16S rRNA sequences in genomes with multiple rrn operons

The level of sequence heterogeneity among rrn operons within genomes determines the accuracy of diversity estimation by 16S rRNA-based methods. Furthermore, the occurrence of widespread horizontal gene transfer (HGT) between distantly related rrn operons casts doubt on reconstructions of phylogenetic relationships. For this study, patterns of distribution of rrn copy numbers, interoperonic divergence, and redundancy of 16S rRNA sequences were evaluated. Bacterial genomes display up to 15 operons and operon numbers up to 7 are commonly found, but ~40% of the organisms analyzed have either one or two operons. Among the Archaea, a single operon appears to dominate and the highest number of operons is five. About 40% of sequences among 380 operons in 76 bacterial genomes with multiple operons were identical to at least one other 16S rRNA sequence in the same genome, and in 38% of the genomes all 16S rRNAs were invariant. For Archaea, the number of identical operons was only 25%, but only five genomes with 21 operons are currently available. These considerations suggest an upper bound of roughly threefold overestimation of bacterial diversity resulting from cloning and sequencing of 16S rRNA genes from the environment; however, the inclusion of genomes with a single rrn operon may lower this correction factor to ~2.5. Divergence among operons appears to be small overall for both Bacteria and Archaea, with the vast majority of 16S rRNA sequences showing <1% nucleotide differences. Only five genomes with operons with a higher level of nucleotide divergence were detected, and Thermoanaerobacter tengcongensis exhibited the highest level of divergence (11.6%) noted to date. Overall, four of the five extreme cases of operon differences occurred among thermophilic bacteria, suggesting a much higher incidence of HGT in these bacteria than in other groups.     

Identification of a copper-responsive two-component system on the chromosome of Escherichia coli K-12

Using a genetic screen we have identified two chromosomal genes, cusRS (ylcA ybcZ), from Escherichia coli K-12 that encode a two-component, signal transduction system that is responsive to copper ions. This regulatory system is required for copper-induced expression of pcoE, a plasmid-borne gene from the E. coli copper resistance operon pco. The closest homologs of CusR and CusS are plasmid-borne two-component systems that are also involved in metal responsive gene regulation: PcoR and PcoS from the pco operon of E. coli; CopR and CopS from the cop operon, which provides copper resistance to Pseudomonas syringae; and SilR and SilS from the sil locus, which provides silver ion resistance to Salmonella enterica serovar Typhimurium. The genes cusRS are also required for the copper-dependent expression of at least one chromosomal gene, designated cusC (ylcB), which is allelic to the recently identified virulence gene ibeB in E. coli K1. The cus locus may comprise a copper ion efflux system, because the expression of cusC is induced by high concentrations of copper ions. Furthermore, the translation products of cusC and additional downstream genes are homologous to known metal ion antiporters.