Genomic analysis of the aromatic catabolic pathways from Pseudomonas putida KT2440

Analysis of the catabolic potential of Pseudomonas putida KT2440 against a wide range of natural aromatic compounds and sequence comparisons with the entire genome of this microorganism predicted the existence of at least four main pathways for the catabolism of central aromatic intermediates, that is, the protocatechuate (pca genes) and catechol (cat genes) branches of the beta-ketoadipate pathway, the homogentisate pathway (hmg/fah/mai genes) and the phenylacetate pathway (pha genes). Two additional gene clusters that might be involved in the catabolism of N-heterocyclic aromatic compounds (nic cluster) and in a central meta-cleavage pathway (pcm genes) were also identified. Furthermore, the genes encoding the peripheral pathways for the catabolism of p-hydroxybenzoate (pob), benzoate (ben), quinate (qui), phenylpropenoid compounds (fcs, ech, vdh, cal, van, acd and acs), phenylalanine and tyrosine (phh, hpd) and n-phenylalkanoic acids (fad) were mapped in the chromosome of P. putida KT2440. Although a repetitive extragenic palindromic (REP) element is usually associated with the gene clusters, a supraoperonic clustering of catabolic genes that channel different aromatic compounds into a common central pathway (catabolic island) was not observed in P. putida KT2440. The global view on the mineralization of aromatic compounds by P. putida KT2440 will facilitate the rational manipulation of this strain for improving biodegradation/biotransformation processes, and reveals this bacterium as a useful model system for studying biochemical, genetic, evolutionary and ecological aspects of the catabolism of aromatic compounds.     

Global features of the Pseudomonas putida KT2440 genome sequence

The compositional bias of the G+C, di- and tetranucleotide contents in the 6 181 862 bp Pseudomonas putida KT2440 genome was analysed in sliding windows of 4000 bp in steps of 1000 bp. The genome has a low GC skew (mean 0.066) between the leading and lagging strand. The values of GC contents (mean 61.6%) and of dinucleotide relative abundance exhibit skewed Gaussian distributions. The variance of tetranucleotide frequencies, which increases linearly with increasing GC content, shows two overlapping Gaussian distributions of genome sections with low (minor fraction) or high variance (major fraction). Eighty per cent of the chromosome shares similar GC contents and oligonucleotide bias, but 105 islands of 4000 bp or more show atypical GC contents and/or oligonucleotide signature. Almost all islands provide added value to the metabolic proficiency of P. putida as a saprophytic omnivore. Major features are the uptake and degradation of organic chemicals, ion transport and the synthesis and secretion of secondary metabolites. Other islands endow P. putida with determinants of resistance and defenceor with constituents and appendages of the cell wall. A total of 29 islands carry the signature of mobile elements such as phage, transposons, insertion sequence (IS) elements and group II introns, indicating recent acquisition by horizontal gene transfer. The largest gene carries the most unusual sequence that encodes a multirepeat threonine-rich surface adhesion protein. Among the housekeeping genes, only genes of the translational apparatus were located in segments with an atypical signature, suggesting that the synthesis of ribosomal proteins is uncoupled from the rapidly changing translational demands of the cell by the separate utilization of tRNA pools.     

Molecular characterization of fprB (ferredoxin-NADP+ reductase) in Pseudomonas putida KT2440

The fpr gene, which encodes a ferredoxin-NADP+ reductase, is known to participate in the reversible redox reactions between NADP+/NADPH and electron carriers, such as ferredoxin or flavodoxin. The role of Fpr and its regulatory protein, FinR, in Pseudomonas putida KT2440 on the oxidative and osmotic stress responses has already been characterized [Lee at al. (2006). Biochem. Biophys. Res. Commun. 339, 1246-1254]. In the genome of P. putida KT2440, another Fpr homolog (FprB) has a 35.3% amino acid identity with Fpr. The fprB gene was cloned and expressed in Escherichia coli. The diaphorase activity assay was conducted using purified FprB to identify the function of FprB. In contrast to the fpr gene, the induction of fprB was not affected by oxidative stress agents, such as paraquat, menadione, H2O2 and t-butyl hydroperoxide. However, a higher level of fprB induction was observed under osmotic stress. Targeted disruption of fprB by homologous recombination resulted in a growth defect under high osmotic conditions. Recovery of oxidatively damaged aconitase activity was faster for the fprB mutant than for the fpr mutant, yet still slower than that for the wild type. Therefore, these data suggest that the catalytic function of FprB may have evolved to augment the function of Fpr in P. putida KT2440.     

Complete genome sequence and comparative analysis of the metabolically versatile Pseudomonas putida KT2440

Pseudomonas putida is a metabolically versatile saprophytic soil bacterium that has been certified as a biosafety host for the cloning of foreign genes. The bacterium also has considerable potential for biotechnological applications. Sequence analysis of the 6.18 Mb genome of strain KT2440 reveals diverse transport and metabolic systems. Although there is a high level of genome conservation with the pathogenic Pseudomonad Pseudomonas aeruginosa (85% of the predicted coding regions are shared), key virulence factors including exotoxin A and type III secretion systems are absent. Analysis of the genome gives insight into the non-pathogenic nature of P. putida and points to potential new applications in agriculture, biocatalysis, bioremediation and bioplastic production.     

Combined Physical and Genetic Map of the Pseudomonas putida KT2440 Chromosome

A combined physical and genetic map of the Pseudomonas putida KT2440 genome was constructed from data obtained by pulsed-field gel electrophoresis techniques (PFGE) and Southern hybridization. Circular genome size was estimated at 6.0 Mb by adding the sizes of 19 SwaI, 9 PmeI, 6 PacI, and 6 I-CeuI fragments. A complete physical map was achieved by combining the results of (i) analysis of PFGE of the DNA fragments resulting from digestion of the whole genome with PmeI, SwaI, I-CeuI, and PacI as well as double digestion with combinations of these enzymes and (ii) Southern hybridization analysis of the whole wild-type genome digested with different enzymes and hybridized against a series of probes obtained as cloned genes from different pseudomonads of rRNA group I and Escherichia coli, as P. putida DNA obtained by PCR amplification based on sequences deposited at the GenBank database, and by labeling of macrorestriction fragments of the P. putida genome eluted from agarose gels. As an alternative, 10 random mini-Tn5-Km mutants of P. putida KT2440 were used as a source of DNA, and the band carrying the mini-Tn5 in each mutant was identified after PFGE of a series of complete chromosomal digestions and hybridization with the kanamycin resistance gene of the mini-Tn5 as a probe. We established a circular genome map with an average resolution of 160 kb. Among the 63 genes located on the genetic map were key markers such as oriC, 6 rrn loci (rnnA to -F), recA, ftsZ, rpoS, rpoD, rpoN, and gyrB; auxotrophic markers; and catabolic genes for the metabolism of aromatic compounds. The genetic map of P. putida KT2440 was compared to those of Pseudomonas aeruginosa PAO1 and Pseudomonas fluorescens SBW25. The chromosomal backbone revealed some similarity in gene clustering among the three pseudomonads but differences in physical organization, probably as a result of intraspecific rearrangements.     

Multiple and interconnected pathways for L-lysine catabolism in Pseudomonas putida KT2440

L-lysine catabolism in Pseudomonas putida KT2440 was generally thought to occur via the aminovalerate pathway. In this study we demonstrate the operation of the alternative aminoadipate pathway with the intermediates D-lysine, L-pipecolate, and aminoadipate. The simultaneous operation of both pathways for the use of L-lysine as the sole carbon and nitrogen source was confirmed genetically. Mutants with mutations in either pathway failed to use L-lysine as the sole carbon and nitrogen source, although they still used L-lysine as the nitrogen source, albeit at reduced growth rates. New genes were identified in both pathways, including the davB and davA genes that encode the enzymes involved in the oxidation of L-lysine to delta-aminovaleramide and the hydrolysis of the latter to delta-aminovalerate, respectively. The amaA, dkpA, and amaB genes, in contrast, encode proteins involved in the transformation of Delta1-piperidine-2-carboxylate into aminoadipate. Based on L-[U-13C, U-15N]lysine experiments, we quantified the relative use of pathways in the wild type and its isogenic mutants. The fate of 13C label of L-lysine indicates that in addition to the existing connection between the D- and L-lysine pathways at the early steps of the catabolism of L-lysine mediated by a lysine racemase, there is yet another interconnection at the lower end of the pathways in which aminoadipate is channeled to yield glutarate. This study establishes an unequivocal relationship between gene and pathway enzymes in the metabolism of L-lysine, which is of crucial importance for the successful colonization of the rhizosphere of plants by this microorganism.     

Genome-scale reconstruction and analysis of the Pseudomonas putida KT2440 metabolic network facilitates applications in biotechnology

A cornerstone of biotechnology is the use of microorganisms for the efficient production of chemicals and the elimination of harmful waste. Pseudomonas putida is an archetype of such microbes due to its metabolic versatility, stress resistance, amenability to genetic modifications, and vast potential for environmental and industrial applications. To address both the elucidation of the metabolic wiring in P. putida and its uses in biocatalysis, in particular for the production of non-growth-related biochemicals, we developed and present here a genome-scale constraint-based model of the metabolism of P. putida KT2440. Network reconstruction and flux balance analysis (FBA) enabled definition of the structure of the metabolic network, identification of knowledge gaps, and pin-pointing of essential metabolic functions, facilitating thereby the refinement of gene annotations. FBA and flux variability analysis were used to analyze the properties, potential, and limits of the model. These analyses allowed identification, under various conditions, of key features of metabolism such as growth yield, resource distribution, network robustness, and gene essentiality. The model was validated with data from continuous cell cultures, high-throughput phenotyping data, (13)C-measurement of internal flux distributions, and specifically generated knock-out mutants. Auxotrophy was correctly predicted in 75% of the cases. These systematic analyses revealed that the metabolic network structure is the main factor determining the accuracy of predictions, whereas biomass composition has negligible influence. Finally, we drew on the model to devise metabolic engineering strategies to improve production of polyhydroxyalkanoates, a class of biotechnologically useful compounds whose synthesis is not coupled to cell survival. The solidly validated model yields valuable insights into genotype-phenotype relationships and provides a sound framework to explore this versatile bacterium and to capitalize on its vast biotechnological potential.     

Engineering Pseudomonas putida KT2440 for simultaneous degradation of carbofuran and chlorpyrifos

Currently, chlorpyrifos (CP) and carbofuran are often applied together to control major agricultural pests in many developing countries, in most cases, they are simultaneously detected in agricultural soils. Some cost‐effective techniques are required for the remediation of combined pollution caused by multiple pesticides. In this work, we aim at constructing a detectable recombinant microorganism with the capacity to simultaneously degrade CP and carbofuran. To achieve this purpose, CP/carbofuran hydrolase genes and gfp were integrated into the chromosome of a biosafety strain Pseudomonas putida KT2440 using a chromosomal scarless modification strategy with upp as a counter‐selectable marker. The toxicity of the hydrolysis products was significantly lower compared with the parent compounds. The recombinant strain could utilize CP or carbofuran as the sole source of carbon for growth. The inoculation of the recombinant strain to soils treated with carbofuran and CP resulted in a higher degradation rate than in noninoculated soils. Introduced green fluorescent protein can be employed as a biomarker to track the recombinant strain during bioremediation. Therefore, the recombinant strain has potential to be applied for in situ bioremediation of soil co‐contaminated with carbofuran and CP.