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.     

A DNA polymerase V homologue encoded by TOL plasmid pWW0 confers evolutionary fitness on Pseudomonas putida under conditions of environmental stress

Plasmids in conjunction with other mobile elements such as transposons are major players in the genetic adaptation of bacteria in response to changes in environment. Here we show that a large catabolic TOL plasmid, pWW0, from Pseudomonas putida carries genes (rulAB genes) encoding an error-prone DNA polymerase Pol V homologue which increase the survival of bacteria under conditions of accumulation of DNA damage. A study of population dynamics in stationary phase revealed that the presence of pWW0-derived rulAB genes in the bacterial genome allows the expression of a strong growth advantage in stationary phase (GASP) phenotype of P. putida. When rulAB-carrying cells from an 8-day-old culture were mixed with Pol V-negative cells from a 1-day-old culture, cells derived from the aged culture out-competed cells from the nonaged culture and overtook the whole culture. At the same time, bacteria from an aged culture lacking the rulAB genes were only partially able to out-compete cells from a fresh overnight culture of the parental P. putida strain. Thus, in addition to conferring resistance to DNA damage, the plasmid-encoded Pol V genes significantly increase the evolutionary fitness of bacteria during prolonged nutritional starvation of a P. putida population. The results of our study indicate that RecA is involved in the control of expression of the pWW0-encoded Pol V.     

Loss of the toluene-xylene catabolic genes of TOL plasmid pWW0 during growth of Pseudomonas putida on benzoate is due to a selective growth advantage of 'cured' segregants

During growth on benzoate-minimal medium Pseudomonas putida mt-2 (PaW1) segregates derivative ('cured') strains which have lost the ability to use the pathway encoded by its resident catabolic plasmid pWW0. Experiments with two plasmids identical to pWW0 but each with an insert of Tn401, which confers resistance to carbenicillin, suggested that the 'benzoate curing' occurs far more frequently by the specific deletion of the 39 kbp region carrying the catabolic genes than by total plasmid loss. This effect was not pH-dependent, and was not produced during growth on other weak organic acids, such as succinate or propionate, or when benzoate was present in the medium with an alternative, preferentially used carbon source such as succinate. Growth on benzoate did not cause loss from strain PaW174 of the plasmid pWW0174, a derivative of pWW0 which has deleted the 39 kbp region but carries Tn401. Similarly the naphthalene-catabolic plasmid pWW60-1, of the same incompatibility group as pWW0, was not lost from PaW701 during growth on benzoate. Competition between wild-type PaW1 and PaW174, which has the 'cured' phenotype, showed that the latter has a distinct growth advantage on benzoate over the wild-type even when initially present as only 1% of the population: when PaW174 was seeded at lower cell ratios, spontaneously 'cured' derivatives of PaW1 took over the culture after 60-80 generations, indicating that they are present in PaW1 cultures at frequencies between 10(-2) and 10(-3). We conclude that the progressive takeover of populations of PaW1 only occurs when benzoate is present as the sole growth source and that neither benzoate, nor other weak acids, affect plasmid segregation or deletion events: a sufficient explanation is that the 'cured' segregants grow faster than the wild-type using the chromosomally determined beta-ketoadipate pathway.     

Nucleotide sequences of the Acinetobacter calcoaceticus benABC genes for benzoate 1,2-dioxygenase reveal evolutionary relationships among multicomponent oxygenases.

The nucleotide sequences of the Acinetobacter calcoaceticus benABC genes encoding a multicomponent oxygenase for the conversion of benzoate to a nonaromatic cis-diol were determined. The enzyme, benzoate 1,2-dioxygenase, is composed of a hydroxylase component, encoded by benAB, and an electron transfer component, encoded by benC. Comparison of the deduced amino acid sequences of BenABC with related sequences, including those for the multicomponent toluate, toluene, benzene, and naphthalene 1,2-dioxygenases, indicated that the similarly sized subunits of the hydroxylase components were derived from a common ancestor. Conserved cysteine and histidine residues may bind a [2Fe-2S] Rieske-type cluster to the alpha-subunits of all the hydroxylases. Conserved histidines and tyrosines may coordinate a mononuclear Fe(II) ion. The less conserved beta-subunits of the hydroxylases may be responsible for determining substrate specificity. Each dioxygenase had either one or two electron transfer proteins. The electron transfer component of benzoate dioxygenase, encoded by benC, and the corresponding protein of the toluate 1,2-dioxygenase, encoded by xylZ, were each found to have an N-terminal region which resembled chloroplast-type ferredoxins and a C-terminal region which resembled several oxidoreductases. These BenC and XylZ proteins had regions similar to certain monooxygenase components but did not appear to be evolutionarily related to the two-protein electron transfer systems of the benzene, toluene, and naphthalene 1,2-dioxygenases. Regions of possible NAD and flavin adenine dinucleotide binding were identified.     

TOL plasmid pWW0 in constructed halobenzoate-degrading Pseudomonas strains: enzyme regulation and DNA structure.

WR211 and WR216 are derivatives of halobenzoate-degrading Pseudomonas sp. strain B13 into which the 117-kilobase TOL degradative plasmid pWW0 has been transferred from Pseudomonas putida mt-2. WR211 has lost the ability to grow on the TOL-specific substrate m-xylene but retains the ability to grow on its metabolite, m-toluate. An analysis of the induction of enzymes was consistent with WR211 carrying a nonfunctional regulatory gene, xy1R, WR216 is a spontaneous derivative of WR211 which grows on one of the TOL substrates and yet expresses the nonspecific toluate oxidase, which enables it to grow on the novel substrate 4-chlorobenzoate. In addition to the xy1R lesion inherited from WR211, WR216 appears to carry a mutation in the structural gene for catechol 2,3-oxygenase, xy1E. The plasmids in both strains were analyzed by restriction endonuclease digestion. pWW0-1211 in WR211 has a large deletion (39 kilobases) compared with pWW0 and appears to be identical to a previously described plasmid (pWW0-8) which encodes none of the TOL degradative functions. pWW0-1216 in WR216 has undergone a major structural reorganization relative to its parent, pWW0-1211. This plasmid has a smaller deletion (19 kilobases), which is staggered relative to the deletion in pWW0-1211, and in addition it has two 3-kilobase insertions of unknown origin, one of which appears to cause the xylE mutation.     

DNA sequence determination of the TOL plasmid (pWW0) xylGFJ genes of Pseudomonas putida: implications for the evolution of aromatic catabolism

The meta operon of the Pseudomonas putida TOL plasmid (pWWO) encodes all enzymes of a meta-cleavage pathway for the metabolism of benzoic acids to Krebs-cycle intermediates. We have determined and analysed the nucleic acid sequence of a 3442 bp region of the meta operon containing the xyl-GFJ genes whose products are involved in the post meta-ring fission transformation of catechols. Homology analysis of the xylGFJ gene products revealed evidence of biochemical relatedness, suggested enzymatic mechanisms, and permitted us to propose evolutionary events which may have generated the current variety of aromatic degradative pathways. The xylG gene, which specifies 2-hydroxymuconic semialdehyde dehydrogenase (HMSD), was found to encode a protein of 51.7 kDa. The predicted protein sequence exhibits significant homology to eukaryotic aldehyde dehydrogenases (ADHs) and to the products of two other Pseudomonas catabolic genes, i.e. xylC and alkH. Expansion of the ADH superfamily to include these prokaryotic enzymes permitted a broader analysis of functionally critical ADH residues and phylogenetic relationships among superfamily members. The importance of three regions of these enzymes previously thought to be critical to ADH activity was reinforced by this analysis. However glutamine-487, also thought to be critical, is less well conserved. The revised ADH phylogeny proposed here suggests early catabolic ADH divergence with subsequent interkingdom gene exchange. The xylF gene, which specifies 2-hydroxymuconic semialdehyde hydrolase (HMSH), was delineated by N-terminal sequence analysis of the purified gene product and is shown to encode a protein of 30.6 kDa. Homology analysis revealed sequence similarity to a chromosomally encoded serine hydrolase, especially in the region of the previously identified active-site serine residue, suggesting that HMSH may also possess a serine hydrolytic enzymatic mechanism. Likewise, the xylJ gene, which specifies 2-hydroxy-pent-2,4-dienoate hydratase (HPH), was delineated by N-terminal sequence analysis of purified HPH, and was found to encode a 23.9 kDa protein. Sequence comparisons revealed that both HMSH and HPH have analogues in the tod gene cluster, which specifies a toluene/benzene degradative pathway. Although the newly identified todF and todJ genes had been at least partially sequenced (Zylstra and Gibson, 1989), the open reading frames had not been positively identified. The presence of todJ provides strong evidence that the reactions following ring fission in the tod pathway are identical to those of the TOL pathway.(ABSTRACT TRUNCATED AT 400 WORDS)     

Chromosomal location of TOL plasmid DNA in Pseudomonas putida.

The soil isolate Pseudomonas putida MW1000 can grow on toluene and other hydrocarbons; in this respect it is similar to strains of Pseudomonas which carry the TOL plasmid. By conjugation experiments, the genes conferring these growth abilities have been shown to be located on the bacterial chromosome, linked to vil and catB. A 56-kilobase segment of the bacterial chromosome of MW strains carrying the TOL genes can transpose to the IncP-1 plasmid R18-18. Physical analysis of these TOL R18-18 hybrids has shown that the TOL segment is almost identical to the same region found in the TOL plasmid pWW0.     

cis-diol dehydrogenases encoded by the TOL pWW0 plasmid xylL gene and the Acinetobacter calcoaceticus chromosomal benD gene are members of the short-chain alcohol dehydrogenase superfamily.

In the aerobic degradation of benzoate by bacteria, benzoate is first dihydroxylated by a ring-hydroxylating dioxygenase to form a cis-diol (1,2-dihydroxycyclohexa-3,4-diene carboxylate) which is subsequently transformed to a catechol by an NAD(+)-dependent cis-diol dehydrogenase. The structural gene for this dehydrogenase, encoded on TOL plasmid pWW0 of Pseudomonas putida (xylL) and that encoded on the chromosome of Acinetobacter calcoaceticus (benD), were sequenced. They encode polypeptides of about 28 kDa in size. These proteins are similar to each other, exhibiting 58% sequence identity. They are also similar to other proteins of at least 20 different functions, which are members of the short-chain alcohol dehydrogenase family. The alignment of these proteins suggest two amino acids, lysine and tyrosine, as catalytically important residues.     

Characterization of five genes in the upper-pathway operon of TOL plasmid pWW0 from Pseudomonas putida and identification of the gene products.

The upper operon of the TOL plasmid pWW0 of Pseudomonas putida encodes a set of enzymes which transform toluene and xylenes to benzoate and toluates. The genetic organization of the operon was characterized by cloning of the upper operon genes into an expression vector and identification of their products in Escherichia coli maxicells. This analysis showed that the upper operon contains at least five genes in the order of xylC-xylM-xylA-xylB-xylN. Between the promoter of the operon and xylC, there is a 1.7-kilobase-long space of DNA in which no gene function was identified. In contrast, most of the DNA between xylC and xylN consists of coding sequences. The xylC gene encodes the 57-kilodalton benzaldehyde dehydrogenase. The xylM and xylA genes encode 35- and 40-kilodalton polypeptides, respectively, which were shown by genetic complementation tests to be subunits of xylene oxygenase. The structural gene for benzyl alcohol dehydrogenase, xylB, encodes a 40-kilodalton polypeptide. The last gene of this operon is xylN, which synthesizes a 52-kilodalton polypeptide of unknown function.     

Substrate specificity of catechol 2,3-dioxygenase encoded by TOL plasmid pWW0 of Pseudomonas putida and its relationship to cell growth.

Catechol 2,3-dioxygenase encoded by TOL plasmid pWW0 of Pseudomonas putida consists of four identical subunits, each containing one ferrous ion. The enzyme catalyzes ring cleavage of catechol, 3-methylcatechol, and 4-methylcatechol but shows only weak activity toward 4-ethylcatechol. Two mutants of catechol 2,3-dioxygenases (4ECR1 and 4ECR6) able to oxidize 4-ethylcatechol, one mutant (3MCS) which exhibits only weak activity toward 3-methylcatechol but retained the ability to cleave catechol and 4-methylcatechol, and one phenotypic revertant of 3MCS (3MCR) which had regained the ability to oxidize 3-methylcatechol were characterized by determining their Km and partition ratio (the ratio of productive catalysis to suicide catalysis). The amino acid substitutions in the four mutant enzymes were also identified by sequencing their structural genes. Wild-type catechol 2,3-dioxygenase was inactivated during the catalysis of 4-ethylcatechol and thus had a low partition ratio for this substrate, whereas the two mutant enzymes, 4ECR1 and 4ECR6, had higher partition ratios for it. Similarly, mutant enzyme 3MCS had a lower partition ratio for 3-methylcatechol than that of 3MCR. Molecular oxygen was required for the inactivation of the wild-type enzyme by 4-ethylcatechol and of 3MCS by 3-methylcatechol, and the inactivated enzymes could be reactivated by incubation with FeSO4 plus ascorbic acid. The enzyme inactivation is thus most likely mechanism based and occurred principally by oxidation and/or removal of the ferrous ion in the catalytic center. In general, partition ratios for catechols lower than 18,000 did not support bacterial growth. A possible meaning of the critical value of the partition ratio is discussed.     

TOL plasmid pWW0 in constructed halobenzoate-degrading Pseudomonas strains: prevention of meta pathway.

The hybrid pathway for chlorobenzoate metabolism was studied in WR211 and WR216, which were derived from Pseudomonas sp. B13 by acquisition of TOL plasmid pWW0 from Pseudomonas putida mt-2. Chlorobenzoates are utilized readily by these strains when meta cleavage of chlorocatechols is suppressed. When WR211 utilizes 3-chlorobenzoate (3CB), the expression of catechol 2,3-dioxygenase (C23O) and the catabolic activities for chloroaromatics via the ortho pathway coexist as a consequence of inactivation of the meta cleavage activity by 3-chlorocatechol. Utilization of 4-chlorobenzoate (4CB) by WR216 presupposes the suppression of C23O by a spontaneous mutation in the structural gene, so that 4-chlorocatechol is not misrouted into the meta pathway. Such C23O- mutants were also selected when WR211 was grown continuously on 3CB. Our data explain why the phenotypic characters 3CB+ and Mtol+ (m-toluate) are compatible, whereas 4CB+ and Mtol+ are incompatible.     

Cloning and expression of Acinetobacter calcoaceticus catBCDE genes in Pseudomonas putida and Escherichia coli.

This report describes the isolation and preliminary characterization of a 5.0-kilobase-pair (kbp) EcoRI DNA restriction fragment carrying the catBCDE genes from Acinetobacter calcoaceticus. The respective genes encode enzymes that catalyze four consecutive reactions in the catechol branch of the beta-ketoadipate pathway: catB, muconate lactonizing enzyme (EC 5.5.1.1); catC, muconolactone isomerase (EC 5.3.3.4); catD, beta-ketoadipate enol-lactone hydrolase (EC 3.1.1.24); and catE, beta-ketoadipate succinyl-coenzyme A transferase (EC 2.8.3.6). In A. calcoaceticus, pcaDE genes encode products with the same enzyme activities as those encoded by the respective catDE genes. In Pseudomonas putida, the requirements for both catDE and pcaDE genes are met by a single set of genes, designated pcaDE. A P. putida mutant with a dysfunctional pcaE gene was used to select a recombinant pKT230 plasmid carrying the 5.0-kbp EcoRI restriction fragment containing the A. calcoaceticus catE structural gene. The recombinant plasmid, pAN1, complemented P. putida mutants with lesions in catB, catC, pcaD, and pcaE genes; the complemented activities were expressed constitutively in the recombinant P. putida strains. After introduction into Escherichia coli, the pAN1 plasmid expressed the activities constitutively but at much lower levels that those found in the P. putida transformants or in fully induced cultures of A. calcoaceticus or P. putida. When placed under the control of a lac promoter on a recombinant pUC13 plasmid in E. coli, the A. calcoaceticus restriction fragment expressed catBCDE activities at levels severalfold higher than those found in fully induced cultures of A. calcoaceticus. Thus there is no translational barrier to expression of the A. calcoaceticus genes at high levels in E. coli. The genetic origin of the cloned catBCDE genes was demonstrated by the fact that the 5.0-kbp EcoRI restriction fragment hybridized with a corresponding fragment from wild-type A. calcoaceticus DNA. This fragment was missing in DNA from an A. calcoaceticus mutant in which the cat genes had been removed by deletion. The properties of the cloned fragment demonstrate physical linkage of the catBCDE genes and suggest that they are coordinately transcribed.     

Genetic analysis of a relaxed substrate specificity aromatic ring dioxygenase, toluate 1,2-dioxygenase, encoded by TOL plasmid pWW0 of Pseudomonas putida

Summary Toluate 1,2-dioxygenase is the first enzyme of a meta-cleavage pathway for the oxidative catabolism of benzoate and substituted benzoates to Krebs cycle intermediates that is specified by TOL plasmid pWW0 of Pseudomonas putida. A collection of derivatives harbouring Tn1000 insertions and defective in toluate dioxygenase have been isolated from pPL392, a pBR322-based hybrid plasmid carrying the TOL plasmid meta-cleavage pathway operon. In parallel, a series of N-methyl-N-nitro-N-nitrosoguanidine-induced mutant plasmids defective in this enzyme activity were isolated from pNM72, a pKT231-based hybrid plasmid carrying the same operon. Pairs of mutant plasmids, consisting of one Tn1000 derivative and one nitrosoguanidine-induced derivative, were used for complementation analysis of toluate dioxygenase in Escherichia coli recA bacteria, in which the formation of 2-hydroxymuconic semialdehyde from benzoate was examined. Four cistrons for toluate 1,2-dioxygenase were thus identified. DNA fragments containing nitrosoguanidine-induced mutant cistrons plus the other meta-cleavage operon genes were cloned into pOT5, an R388-based vector, and complementation tests between different nitrosoguanidine-induced mutant cistrons were carried out in Pseudomonas putida cells, this time scoring for growth on p-toluate. This analysis also identified four cistrons. Examination of the products of these cistrons, by means of E. coli minicells containing pPL392 or its Tn1000 insertion derivatives, indicated that the first two cistrons of the operon comprise a single gene, xylX, which encodes a 57 kilodalton protein, and that the third cistron, xy/Y, encodes a 20 kilodalton protein.     

Purification and characterisation of the NADH:acceptor reductase component of xylene monooxygenase encoded by the TOL plasmid pWW0 of Pseudomonas putida mt-2.

The xylene monooxygenase system encoded by the TOL plasmid pWW0 of Pseudomonas putida catalyses the hydroxylation of a methyl side-chain of toluene and xylenes. Genetic studies have suggested that this monooxygenase consists of two different proteins, products of the xylA and xylM genes, which function as an electron-transfer protein and a terminal hydroxylase, respectively. In this study, the electron-transfer component of xylene monooxygenase, the product of xylA, was purified to homogeneity. Fractions containing the xylA gene product were identified by its NADH:cytochrome c reductase activity. The molecular mass of the enzyme was determined to be 40 kDa by SDS/PAGE, and 42 kDa by gel filtration. The enzyme was found to contain 1 mol/mol of tightly but not covalently bound FAD, as well as 2 mol/mol of non-haem iron and 2 mol/mol of acid-labile sulfide, suggesting the presence of two redox centers, one FAD and one [2Fe-2S] cluster/protein molecule. The oxidised form of the protein had absorbance maxima at 457 nm and 390 nm, with shoulders at 350 nm and 550 nm. These absorbance maxima disappeared upon reduction of the protein by NADH or dithionite. The NADH:acceptor reductase was capable of reducing either one- or two-electron acceptors, such as horse heart cytochrome c or 2,6-dichloroindophenol, at an optimal pH of 8.5. The reductase was found to have a Km value for NADH of 22 microM. The oxidation of NADH was determined to be stereospecific; the enzyme is pro-R (class A enzyme). The titration of the reductase with NADH or dithionite yielded three distinct reduced forms of the enzyme: the reduction of the [2Fe-2S] center occurred with a midpoint redox potential of -171 mV; and the reduction of FAD to FAD. (semiquinone form), with a calculated midpoint redox potential of -244 mV. The reduction of FAD. to FAD.. (dihydroquinone form), the last stage of the titration, occurred with a midpoint redox potential of -297 mV. The [2Fe-2S] center could be removed from the protein by treatment with an excess of mersalyl acid. The [2Fe-2S]-depleted protein was still reduced by NADH, giving rise to the formation of the anionic flavin semiquinone observed in the native enzyme, thus suggesting that the electron flow was NADH --> FAD --> [2Fe-2S] in this reductase. The resulting protein could no longer reduce cytochrome c, but could reduce 2,6-dichloroindophenol at a reduced rate.     

Genetic, functional and sequence analysis of the xylR and xylS regulatory genes of the TOL plasmid pWW0

Mutant derivatives of a plasmid, pCF20, which carries the XhoI-D fragment of the TOL plasmid pWW0 have been isolated using Tn5 transposon mutagenesis. Insertion mutations of the xylR and xylS regulatory genes of the catabolic pathway have been isolated and characterized and their ability to induce catechol 2,3-oxygenase activity determined. Analysis of the insertion mutants and also segments of the XhoI-D fragment cloned into plasmid pUC8 in maxicells has identified a 68 kDa polypeptide product encoded by the xylR gene. No clear candidate for the xylS polypeptide was observed. The nucleotide sequence of the xylS region, the intergenic region and part of the xylR region has been determined and open reading frames (ORFs) assigned for both genes. The ORF designated xylS appears capable of encoding a polypeptide of approximately 37 kDa.     

Evolutionary conservation of genes coding for meta pathway enzymes within TOL plasmids pWW0 and pWW53.

Pseudomonas putida MT53 contains a TOL plasmid, pWW53, that encodes toluene-xylene catabolism. pWW53 is nonconjugative, is about 105 to 110 kilobase pairs (kbp) in size, and differs significantly in its restriction endonuclease digestion pattern and incompatibility group from the archetypal TOL plasmid pWW0. An RP4::pWW53 cointegrate plasmid, pWW53-4, containing about 35 kbp of pWW53 DNA, including the entire catabolic pathway genes, was formed, and a restriction map for KpnI, HindIII, and BamHI was derived. The entire regulated meta pathway genes for the catabolism of m-toluate were cloned into pKT230 from pWW53 on a 17.5-kbp HindIII fragment. The recombinant plasmid supported growth on m-toluate when mobilized into plasmid-free P. putida PaW130. A restriction map of the insert for 10 restriction enzymes was derived, and the locations of xylD, xylL, xylE, xylG, and xylF were determined by subcloning and assaying for their gene products in both Escherichia coli and P. putida hosts. Good induction of the enzymes by m-toluate and m-methylbenzyl alcohol but not by m-xylene was measured in P. putida, but little or no regulation was found in E. coli. The restriction map and the gene order showed strong similarities with published maps of the DNA encoding both the entire meta pathway operon (xylDLEGFJIH) and the regulatory genes xylS and xylR on the archetype TOL plasmid pWW0, suggesting a high degree of conservation in DNA structure for the catabolic operon on the two different plasmids.     

Comparison of the nucleotide sequences of the meta-cleavage pathway genes of TOL plasmid pWW0 from Pseudomonas putida with other meta-cleavage genes suggests that both single and multiple nucleotide substitutions contribute to enzyme evolution

TOL plasmid pWW0 from Pseudomonas putida mt-2 encodes catabolic enzymes required for the oxidation of toluene and xylenes. The structural genes for these catabolic enzymes are clustered into two operons, the xylCMABN operon, which encodes a set of enzymes required for the transformation of toluene/xylenes to benzoate/toluates, and the xylXYZLTEGFJQKIH operon, which encodes a set of enzymes required for the transformation of benzoate/toluates to Krebs cycle intermediates. The latter operon can be divided physically and functionally into two parts, the xylXYZL cluster, which is involved in the transformation of benzoate/toluates to (methyl)catechols, and the xylTEGFJQKIH cluster, which is involved in the transformation of (methyl)catechols to Krebs cycle intermediates. Genes isofunctional to xylXYZL are present in Acinetobacter calcoaceticus, and constitute a benzoate-degradative pathway, while xylTEGFJQKIH homologous encoding enzymes of a methylphenol-degradative pathway and a naphthalene-degradative pathway are present on plasmid pVI150 from P. putida CF600, and on plasmid NAH7 from P. putida PpG7, respectively. Comparison of the nucleotide sequences of the xylXYZLTEGFJQKIH genes with other isofunctional genes suggested that the xylTEGFJQKIH genes on the TOL plasmid diverged from these homologues 20 to 50 million years ago, while the xylXYZL genes diverged from the A. calcoaceticus homologues 100 to 200 million years ago. In codons where amino acids are not conserved, the substitution rate in the third base was higher than that in synonymous codons. This result was interpreted as indicating that both single and multiple nucleotide substitutions contributed to the amino acid-substituting mutations, and hence to enzyme evolution. This observation seems to be general because mammalian globin genes exhibit the same tendency.