Codon usage patterns suggest independent evolution of two catabolic operons on toluene-degradative plasmid TOL pWW0 of Pseudomonas putida

Life Within a Contaminated Niche: Comparative Genomic Analyses of an Integrative Conjugative Element ICEnahCSV86 and Two Genomic Islands From Pseudomonas bharatica CSV86T Suggest Probable Role in Colonization and Adaptation

Comparative genomic and functional analyses revealed the presence of three genomic islands (GIs, >50 Kb size): ICEnahCSV86, Pseudomonas bharatica genomic island-1 (PBGI-1), and PBGI-2 in the preferentially aromatic-degrading soil bacterium, Pseudomonas bharatica CSV86^T. Site-specific genomic integration at or near specific transfer RNAs (tRNAs), near-syntenic structural modules, and phylogenetic relatedness indicated their evolutionary lineage to the type-4 secretion system (T4SS) ICEclc family, thus predicting these elements to be integrative conjugative elements (ICEs). These GIs were found to be present as a single copy in the genome and the encoded phenotypic traits were found to be stable, even in the absence of selection pressure. ICEnahCSV86 harbors naphthalene catabolic (nah-sal) cluster, while PBGI-1 harbors Co-Zn-Cd (czc) efflux genes as cargo modules, whereas PBGI-2 was attributed to as a mixed-function element. The ICEnahCSV86 has been reported to be conjugatively transferred (frequency of 7 × 10^–8/donor cell) to Stenotrophomonas maltophilia CSV89. Genome-wide comparative analyses of aromatic-degrading bacteria revealed nah-sal clusters from several Pseudomonas spp. as part of probable ICEs, syntenic to conjugatively transferable ICEnahCSV86 of strain CSV86^T, suggesting it to be a prototypical element for naphthalene degradation. It was observed that the plasmids harboring nah-sal clusters were phylogenetically incongruent with predicted ICEs, suggesting genetic divergence of naphthalene metabolic clusters in the Pseudomonas population. Gene synteny, divergence estimates, and codon-based Z-test indicated that ICEnahCSV86 is probably derived from PBGI-2, while multiple recombination events masked the ancestral lineage of PBGI-1. Diversifying selection pressure (dN-dS = 2.27–4.31) imposed by aromatics and heavy metals implied the modular exchange-fusion of various cargo clusters through events like recombination, rearrangement, domain reshuffling, and active site optimization, thus allowing the strain to evolve, adapt, and maximize the metabolic efficiency in a contaminated niche. The promoters (Pnah and Psal) of naphthalene cargo modules (nah, sal) on ICEnahCSV86 were proved to be efficient for heterologous protein expression in Escherichia coli. GI-based genomic plasticity expands the metabolic spectrum and versatility of CSV86^T, rendering efficient adaptation to the contaminated niche. Such isolate(s) are of utmost importance for their application in bioremediation and are the probable ideal host(s) for metabolic engineering.

A computational system for identifying operons based on RNA-seq data

An operon is a set of neighboring genes in a genome that is transcribed as a single polycistronic message. Genes that are part of the same operon often have related functional roles or participate in the same metabolic pathways. The majority of all bacterial genes are co-transcribed with one or more other genes as part of a multi-gene operon. Thus, accurate identification of operons is important in understanding co-regulation of genes and their functional relationships. Here, we present a computational system that uses RNA-seq data to determine operons throughout a genome. The system takes the name of a genome and one or more files of RNA-seq data as input. Our method combines primary genomic sequence information with expression data from the RNA-seq files in a unified probabilistic model in order to identify operons. We assess our method's ability to accurately identify operons in a range of species through comparison to external databases of operons, both experimentally confirmed and computationally predicted, and through focused experiments that confirm new operons identified by our method. Our system is freely available at https://cs.wellesley.edu/~btjaden/Rockhopper/.

Degradation of toluene by ortho cleavage enzymes in Burkholderia fungorum FLU100

Burkholderia fungorum FLU100 simultaneously oxidized any mixture of toluene, benzene and mono-halogen benzenes to (3-substituted) catechols with a selectivity of nearly 100%. Further metabolism occurred via enzymes of ortho cleavage pathways with complete mineralization. During the transformation of 3-methylcatechol, 4-carboxymethyl-2-methylbut-2-en-4-olide (2-methyl-2-enelactone, 2-ML) accumulated transiently, being further mineralized only after a lag phase of 2 h in case of cells pre-grown on benzene or mono-halogen benzenes. No lag phase, however, occurred after growth on toluene. Cultures inhibited by chloramphenicol after growth on benzene or mono-halogen benzenes were unable to metabolize 2-ML supplied externally, even after prolonged incubation. A control culture grown with toluene did not show any lag phase and used 2-ML as a substrate. This means that 2-ML is an intermediate of toluene degradation and converted by specific enzymes. The conversion of 4-methylcatechol as a very minor by-product of toluene degradation in strain FLU100 resulted in the accumulation of 4-carboxymethyl-4-methylbut-2-en-4-olide (4-methyl-2-enelactone, 4-ML) as a dead-end product, excluding its nature as a possible intermediate. Thus, 3-methylcyclohexa-3,5-diene-1,2-diol, 3-methylcatechol, 2-methyl muconate and 2-ML were identified as central intermediates of productive ortho cleavage pathways for toluene metabolism in B. fungorum FLU100.

Transcriptome dynamics-based operon prediction in prokaryotes

Background

Inferring operon maps is crucial to understanding the regulatory networks of prokaryotic genomes. Recently, RNA-seq based transcriptome studies revealed that in many bacterial species the operon structure vary with the change of environmental conditions. Therefore, new computational solutions that use both static and dynamic data are necessary to create condition specific operon predictions.

Results

In this work, we propose a novel classification method that integrates RNA-seq based transcriptome profiles with genomic sequence features to accurately identify the operons that are expressed under a measured condition. The classifiers are trained on a small set of confirmed operons and then used to classify the remaining gene pairs of the organism studied. Finally, by linking consecutive gene pairs classified as operons, our computational approach produces condition-dependent operon maps. We evaluated our approach on various RNA-seq expression profiles of the bacteria Haemophilus somni, Porphyromonas gingivalis, Escherichia coli and Salmonella enterica. Our results demonstrate that, using features depending on both transcriptome dynamics and genome sequence characteristics, we can identify operon pairs with high accuracy. Moreover, the combination of DNA sequence and expression data results in more accurate predictions than each one alone.

Conclusion

We present a computational strategy for the comprehensive analysis of condition-dependent operon maps in prokaryotes. Our method can be used to generate condition specific operon maps of many bacterial organisms for which high-resolution transcriptome data is available.

Complete Sequence Determination Combined with Analysis of Transposition/Site-specific Recombination Events to Explain Genetic Organization of IncP-7 TOL Plasmid pWW53 and Related Mobile Genetic Elements

Recent studies have indicated that the evolutionarily common catabolic gene clusters are loaded on structurally diverse toluene-catabolic (TOL) plasmids and their residing transposons. To elucidate the mechanisms supporting the diversification of catabolic plasmids and transposons, we determined here the complete 107,929 bp sequence of pWW53, a TOL plasmid from Pseudomonas putida MT53. pWW53 was found to belong to the IncP-7 incompatibility group that play important roles in the catabolism of several xenobiotics. pWW53 carried two distinct transposase-resolvase gene clusters (tnpAR modules), five short terminal inverted repeats (IRs), and three site-specific resolution (res) sites that are all typical of class II transposons. This organization of pWW53 suggested the four possible transposable regions, Tn4657 to Tn4660. The largest 86 kb region (Tn4657) spanned the three other regions, and Tn4657 and Tn4660 (62 kb) covered all of the 36 xyl genes for toluene catabolism. Our subsequent transposition experiments clarified that the three transposons, Tn4657 to Tn4659, indeed exhibit their transposability, and that pWW53 also generated another 37 kb toluene-catabolic transposon, Tn4656, which carried the two separated and inversely oriented segments of pWW53: the tnpRA-IR module of Tn4658 and a part of xyl gene clusters on Tn4657. The Tn4658 transposase was able to mediate the transposition of Tn4658, Tn4657, and Tn4656, while the Tn4659 transposase catalyzed only the transposition of Tn4659. Tn4656 was formed by the Tn4658 resolvase-mediated site-specific inversion between the two inversely oriented res sites on pWW53. These findings and comparison with other catabolic plasmids clearly indicate multiple copies of transposition-related genes and sites on one plasmid and their recombination activities contribute greatly to the diversification of plasmid structures as well as wide dissemination of the evolutionary common gene clusters in various plasmids.

The evolution of IncP catabolic plasmids

The recent adoption of whole plasmid genome sequencing as a routine analytical technique has provided the basis for cataloging the historical events through which plasmids are assembled from the available families of modular plasmid components. Horizontal gene transfer mediated by plasmids plays an important role in the adaptation of bacteria to the presence of specific metabolizable compounds, including man-made chemicals, in the surrounding environment. Bacterial plasmid genome sequence comparisons indicate that plasmids have complex genetic histories resulting from transposition, homologous recombination, and illegitimate recombinational events. Evidence from IncP plasmid genome sequences indicates that cryptic plasmid backbones acquire diverse catabolic pathways through gene capture and horizontal gene transfer.

Molecular characterization of an inducible gentisate 1,2-dioxygenase gene, xlnE, from Pseudomonas alcaligenes NCIMB 9867

Pseudomonas alcaligenes NCIMB 9867 (strain P25X) produces isofunctional enzymes of the gentisate pathway that enables the degradation of xylenols and cresols via gentisate. Previous reports had indicated that one set of enzymes is constitutively expressed whereas the other set is strictly inducible by aromatic hydrocarbon substrates. The gene encoding gentisate 1,2-dioxygenase (GDO), the enzyme that catalyzes the cleavage of the gentisate aromatic ring, was cloned from strain P25X. The GDO gene, designated xlnE, is 1,044 bp, and is part of a 5.4 kb operon which consists of six genes, xlnEFGHID. Transcription of this operon was driven by a sigma 70-type promoter, PxlnE, located 123 bp upstream of the xlnE start codon. Primer extension analysis showed that the xlnE transcription start point is located at the -87 adenine residue. In a P25X xlnE knockout mutant, GDO activity could only be detected when cells were grown in the presence of aromatic substrates, suggesting that xlnE encodes for the constitutive copy of GDO. This was verified by constructing a P25X strain with xlnE transcriptionally fused to a promoterless catechol 2,3-dioxygenase gene. In this strain, catechol 2,3-dioxygenase activity was detected in cells that were grown in the absence of aromatic inducers. However, catechol 2,3-dioxygenase activity increased up to four fold when these cells were grown in the presence of aromatic substrates, in particular 3-hydroxybenzoate. Thus, xlnE is in fact, inducible and the constitutive activity observed under non-inducing conditions was due to its relatively high basal levels of expression.

Complete sequence of the IncP-9 TOL plasmid pWW0 from Pseudomonas putida

The TOL plasmid pWW0 (117 kb) is the best studied catabolic plasmid and the archetype of the IncP-9 plasmid incompatibility group from Pseudomonas. It carries the degradative (xyl) genes for toluenes and xylenes within catabolic transposons Tn4651 and Tn4653. Analysis of the complete pWW0 nucleotide sequence revealed 148 putative open reading frames. Of these, 77 showed similarity to published sequences in the available databases predicting functions for: plasmid replication, stable maintenance and transfer; phenotypic determinants; gene regulation and expression; and transposition. All identifiable transposition functions lay within the boundaries of the 70 kb transposon Tn4653, leaving a 46 kb sector containing all the IncP-9 core functions. The replicon and stable inheritance region was very similar to the mini-replicon from IncP-9 antibiotic resistance plasmid pM3, with their Rep proteins forming a novel group of initiation proteins. pWW0 transfer functions exist as two blocks encoding putative DNA processing and mating pair formation genes, with organizational and sequence similarity to IncW plasmids. In addition to the known Tn4651 and IS1246 elements, two additional transposable elements were identified as well as several putative transposition functions, which are probably genetic remnants from previous transposition events. Genes likely to be responsible for known resistance to ultraviolet light and free radicals were identified. Other putative phenotypic functions identified included resistance to mercury and other metal ions, as well as to quaternary ammonium compounds. The complexity and size of pWW0 is largely the result of the mosaic organization of the transposable elements that it carries, rather than the backbone functions of IncP-9 plasmids.

Marinobacter strain NCE312 has a Pseudomonas-like naphthalene dioxygenase

One strain of bacteria, designated NCE312, was isolated from a naphthalene-digesting chemostat culture that was inoculated with creosote-contaminated marine sediment. The strain was isolated based on its ability to grow using naphthalene as a sole carbon source. In addition, the strain degraded 2-methylnaphthalene and 1-methylnaphthalene. Analysis of a 16S rRNA gene sequence from NCE312 placed the isolate in the genus Marinobacter. Degenerate PCR primers were used to amplify a fragment of a naphthalene 1,2-dioxygenase large subunit gene. A phylogenetic analysis indicated the Marinobacter naphthalene dioxygenase is similar to those from Pseudomonas and Burkholderia strains suggesting that the dioxygenase gene may have been transferred horizontally between these lineages of bacteria.

Molecular diversity of plasmids bearing genes that encode toluene and xylene metabolism in Pseudomonas strains isolated from different contaminated sites in Belarus

Twenty different Pseudomonas strains utilizing m-toluate were isolated from oil-contaminated soil samples near Minsk, Belarus. Seventeen of these isolates carried plasmids ranging in size from 78 to about 200 kb (assigned pSVS plasmids) and encoding the meta cleavage pathway for toluene metabolism. Most plasmids were conjugative but of unknown incompatibility groups, except for one, which belonged to the IncP9 group. The organization of the genes for toluene catabolism was determined by restriction analysis and hybridization with xyl gene probes of pWW0. The majority of the plasmids carried xyl-type genes highly homologous to those of pWW53 and organized in a similar manner (M. T. Gallegos, P. A. Williams, and J. L. Ramos, J. Bacteriol. 179:5024-5029, 1997), with two distinguishable meta pathway operons, one upper pathway operon, and three xylS-homologous regions. All of these plasmids also possessed large areas of homologous DNA outside the catabolic genes, suggesting a common ancestry. Two other pSVS plasmids carried only one meta pathway operon, one upper pathway operon, and one copy each of xylS and xylR. The backbones of these two plasmids differed greatly from those of the others. Whereas these parts of the plasmids, carrying the xyl genes, were mostly conserved between plasmids of each group, the noncatabolic parts had undergone intensive DNA rearrangements. DNA sequencing of specific regions near and within the xylTE and xylA genes of the pSVS plasmids confirmed the strong homologies to the xyl genes of pWW53 and pWW0. However, several recombinations were discovered within the upper pathway operons of the pSVS plasmids and pWW0. The main genetic mechanisms which are thought to have resulted in the present-day configuration of the xyl operons are discussed in light of the diversity analysis carried out on the pSVS plasmids.

Cloning and expression of ntnD, encoding a novel NAD(P)(+)-independent 4-nitrobenzyl alcohol dehydrogenase from Pseudomonas sp. Strain TW3

Pseudomonas sp. strain TW3 is able to metabolize 4-nitrotoluene to 4-nitrobenzoate and toluene to benzoate aerobically via a route analogous to the upper pathway of the TOL plasmids. We report the cloning and characterization of a benzyl alcohol dehydrogenase gene (ntnD) which encodes the enzyme for the catabolism of 4-nitrobenzyl alcohol and benzyl alcohol to 4-nitrobenzaldehyde and benzaldehyde, respectively. The gene is located downstream of the previously reported ntn gene cluster. NtnD bears no similarity to the analogous TOL plasmid XylB (benzyl alcohol dehydrogenase) protein either in its biochemistry, being NAD(P)(+) independent and requiring assay via dye-linked electron transfer, or in its deduced amino acid sequence. It does, however, have significant similarity in its amino acid sequence to other NAD(P)(+)-independent alcohol dehydrogenases and contains signature patterns characteristic of type III flavin adenine dinucleotide-dependent alcohol oxidases. Reverse transcription-PCR demonstrated that ntnD is transcribed during growth on 4-nitrotoluene, although apparently not as part of the same transcript as the other ntn genes. The substrate specificity of the enzyme expressed from the cloned and overexpressed gene was similar to the activity expressed from strain TW3 grown on 4-nitrotoluene, providing evidence that ntnD is the previously unidentified gene in the pathway of 4-nitrotoluene catabolism. Examination of the 14.8-kb region around the ntn genes suggests that one or more recombination events have been involved in the formation of their current organization.

Complete nucleotide sequence and evolutionary significance of a chromosomally encoded naphthalene-degradation lower pathway from Pseudomonas stutzeri AN10

Pseudomonas stutzeri strain AN10 is a naphthalene-degrading strain whose dissimilatory genes are chromosomally encoded. We sequenced the entire naphthalene-degradation lower pathway of P. stutzeri AN10, this being, together with the upper-pathway reported previously (Bosch R. et al., 1999a. Gene 236, 149-157) the first complete DNA sequence for an entire naphthalene-catabolic pathway. Eleven open reading frames were identified. The nahGTHINLOMKJ genes encode enzymes for the metabolism of salicylate to pyruvate and acetyl-CoA, and nahR encodes the NahR regulatory protein. Our findings suggest that catabolic modules were recruited through transposition events and recombination among tnpA-like genes, and subsequent rearrangements and deletions of non-essential DNA fragments allowed the formation of the actual catabolic pathway. Our results also suggest that the genes encoding the xylene/toluene-degradation enzymes of P. putida mt-2 (pWW0) have coexisted with the nah genes of the P. stutzeri AN10 ancestral genome. This could allow the selection, via recombination events among homologous genes, for a combination of genes enabling the metabolism of a given aromatic compound in the ancestral host strain. Such events accelerate the evolution of modern catabolic pathways and provide new genetic material to the environment, ultimately resulting in improved, natural, bioremediation potential.

Genetic characterization and evolutionary implications of a chromosomally encoded naphthalene-degradation upper pathway from Pseudomonas stutzeri AN10

Pseudomonas stutzeri strain AN10 is a naphthalene-degrading strain whose dissimilatory genes are chromosomally encoded. We sequenced a total of 11514bp including the entire naphthalene-degradation upper pathway (nah) of P. stutzeri AN10. Nine open reading frames, nahAaAbAcAdBFCED, encoding the enzymes for the degradation of naphthalene to salicylate, were identified. The nah genes of P. stutzeri AN10 have been compared with genes encoding isofunctional proteins from other Pseudomonas naphthalene-degradation upper pathways. The implications of the sequence homologies to the evolution of aromatic catabolic pathways are discussed. Our findings indicate that this entire catabolic module of P. stutzeri AN10 was recruited from other microorganisms and a short period of time has elapsed after its incorporation within the P. stutzeri AN10 genome. Comparisons also suggest the coexistence of two entire nah upper pathways in a host strain, and further recombination between them. These events could accelerate the evolution of modern catabolic pathways.

ntn genes determining the early steps in the divergent catabolism of 4-nitrotoluene and toluene in Pseudomonas sp. strain TW3

Pseudomonas sp. strain TW3 is able to oxidatively metabolize 4-nitrotoluene and toluene via a route analogous to the upper pathway of the TOL plasmids. We report the sequence and organization of five genes, ntnWCMAB*, which are very similar to and in the same order as the xyl operon of TOL plasmid pWW0 and present evidence that they encode enzymes which are expressed during growth on both 4-nitrotoluene and toluene and are responsible for their oxidation to 4-nitrobenzoate and benzoate, respectively. These genes encode an alcohol dehydrogenase homolog (ntnW), an NAD+-linked benzaldehyde dehydrogenase (ntnC), a two-gene toluene monooxygenase (ntnMA), and part of a benzyl alcohol dehydrogenase (ntnB*), which have 84 to 99% identity at the nucleotide and amino acid levels with the corresponding xylWCMAB genes. The xylB homolog on the TW3 genome (ntnB*) appears to be a pseudogene and is interrupted by a piece of DNA which destroys its functional open reading frame, implicating an additional and as-yet-unidentified benzyl alcohol dehydrogenase gene in this pathway. This conforms with the observation that the benzyl alcohol dehydrogenase expressed during growth on 4-nitrotoluene and toluene differs significantly from the XylB protein, requiring assay via dye-linked electron transfer rather than through a nicotinamide cofactor. The further catabolism of 4-nitrobenzoate and benzoate diverges in that the former enters the hydroxylaminobenzoate pathway as previously reported, while the latter is further metabolized via the beta-ketoadipate pathway.