Previous Article | Next Article ![]()
Journal of Bacteriology, March 2007, p. 1931-1945, Vol. 189, No. 5
0021-9193/07/$08.00+0 doi:10.1128/JB.01259-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
,
*
Anu Y. Chakicherla,1,
Patrick S. G. Chain,1,4
Radomir Schmidt,2
Maria W. Shin,1
Tina C. Legler,1
Kate M. Scow,2
Frank W. Larimer,3,4
Susan M. Lucas,4
Paul M. Richardson,4 and
Krassimira R. Hristova2
Lawrence Livermore National Laboratory, Livermore, California,1 Department of Land Air and Water Resources, University of California, Davis, California,2 Genome Analysis Group, Oak Ridge National Laboratory, Oak Ridge, Tennessee,3 Joint Genome Institute Production Genomics Facility, Walnut Creek, California4
Received 10 August 2006/ Accepted 29 November 2006
|
|
|---|
4-Mb circular chromosome and an
600-kb megaplasmid, containing 3,831 and 646 genes, respectively. Aromatic hydrocarbon and alkane degradation, metal resistance, and methylotrophy are encoded on the chromosome. The megaplasmid contains an unusual t-RNA island, numerous insertion sequences, and large repeated elements, including a 40-kb region also present on the chromosome and a 29-kb tandem repeat encoding phosphonate transport and cobalamin biosynthesis. The megaplasmid also codes for alkane degradation and was shown to play an essential role in MTBE degradation through plasmid-curing experiments. Discrepancies between the insertion sequence element distribution patterns, the distributions of best BLASTP hits among major phylogenetic groups, and the G+C contents of the chromosome (69.2%) and plasmid (66%), together with comparative genome hybridization experiments, suggest that the plasmid was recently acquired and apparently carries the genetic information responsible for PM1's ability to degrade MTBE. Comparative genomic hybridization analysis with two PM1-like MTBE-degrading environmental isolates (
99% identical 16S rRNA gene sequences) showed that the plasmid was highly conserved (ca. 99% identical), whereas the chromosomes were too diverse to conduct resequencing analysis. PM1's genome sequence provides a foundation for investigating MTBE biodegradation and exploring the genetic regulation of multiple biodegradation pathways in M. petroleiphilum and other MTBE-degrading beta-proteobacteria. |
|
|---|
99% similarity to PM1, based on the 16S rRNA gene) (38, 67, 82). PM1-like bacteria occur naturally in a number of MTBE-contaminated aquifers in the United States (46, 47, 82), Mexico (21), and Europe (55, 61), and their presence has been correlated with MTBE degradation activity in numerous sites (47, 67, 82), using real-time PCR analysis (37). These results suggest that PM1-like organisms may play a major role in MTBE biodegradation under aerobic conditions in contaminated aquifers. The genetic basis for MTBE metabolism is not currently understood, although there is general agreement that the initial enzymatic steps are similar to cometabolic degradation pathways (27, 68, 74), and recent reports have described genes involved in degradation of the MTBE downstream metabolites 2-methyl-1,2-propanediol (52) and 2-hydroxyisobutyrate (61). The complex regulation of the metabolism of fuel hydrocarbons and MTBE, often occurring in mixtures, is relatively unknown, while limited studies showed that MTBE degradation could be inhibited in mixtures with BTEX (benzene, toluene, ethylbenzene, xylene) compounds (20, 47).
In this paper, we present our analysis of the whole genome sequence of M. petroleiphilum PM1. We present comparative sequence analysis results between PM1 and bacteria with homologous individual genes and operons as well as comparative whole-genome hybridization analysis between PM1 and PM1-like MTBE-degrading isolates (
99% identical 16S rRNA gene sequences) from gasoline-contaminated sites. General genome features are discussed, including interesting repeated elements as well as genes and operons involved in methylotrophy, degradation of aromatic hydrocarbons, degradation of cyclic and straight-chain alkanes, cofactor biosynthesis, motility, secretion, and heavy metal resistance and transport. A noteworthy finding was the presence of a large,
600-bp plasmid in PM1 that was highly conserved among PM1-like bacteria. Furthermore, plasmid-curing experiments showed that the plasmid was essential for MTBE and tert-butyl alcohol (TBA) biodegradation in PM1. The PM1 genome sequence has provided a foundation for understanding novel pathways and interactions in this important subsurface bacterium as well as those in phylogenetically similar MTBE-degrading bacteria.
|
|
|---|
Sequencing, gene prediction, and annotation. Genomic DNA was isolated and purified from M. petroleiphilum PM1, and whole-genome shotgun libraries (3-kb, 8-kb, and 40-kb DNA inserts) were constructed and sequenced as previously described (14). After quality control of the 90,327 total initial reads of draft sequence, 83,180 sequences were assembled, producing an average of 10.7-fold coverage across the genome. The whole genome sequence was assembled using the Phred/Phrap/Consed package (P. Green, University of Washington) (25, 26, 32). The reads were assembled into 24 high-quality draft sequence contigs, which were linked into three larger scaffolds by paired-end sequence information. Gaps in the sequence were closed either by walking on gap-spanning clones or with PCR products generated from genomic DNA. Physical (uncaptured) gaps were closed by combinatorial (multiplex) PCR. Sequence finishing and polishing added 308 reads, and the final assessment of the genome assembly was completed as described previously (14). The final genome assembly quality of PM1 adheres to conventional standards of less than one error per 10,000 bp. Each base is covered by at least two quality sequences, with an average of 10.7-fold coverage. Proper assembly was verified by fosmid coverage coupled with PCR data. Gene modeling and genome annotation were performed as previously described (14) to identify open reading frames likely encoding proteins (coding sequences [CDSs]).
Comparative genomics. Orthologs and CDSs unique to M. petroleiphilum PM1 were identified using the Integrated Microbial Genomes system from the Joint Genome Institute. Results were based on BLASTP analysis, with cutoff values of <105 (E value) and 30% identity for orthologs and <102 (E value) and 20% identity for unique CDSs.
Phylogenetic tree analysis. Homologs of M. petroleiphilum PM1 translated CDSs were identified using BLASTP searches against the nonredundant GenBank database from the National Center for Biotechnology Information. Sequences were aligned and alignments were refined using ClustalX, version 1.8 (41), along with manual adjustments. The Protdist program and the Neighbor program of the Phylip package (28) were used to generate the phylogenetic tree for MpeA3393. The pairwise parameters included a gap opening of 35 and a gap extension of 0.75. The multiple alignment parameters included a gap opening of 15, a gap extension of 0.3, delay divergence of 30%, and use of the Gonnet series for the protein weight matrix.
Comparative genomic hybridization and comparative genomic sequencing analyses.
Comparative genomic hybridization was conducted in order to analyze the conservation of genes from MTBE-degrading isolates MG4 and 312 with those from PM1 across the entire genome. High-density arrays (
400,000 oligomers) were designed and produced by NimbleGen Systems, Inc. (Madison, WI), using 29-mer probes specific for every 26 bp for both strands of the entire genome and for every 7 bases for both strands of the plasmid. Arrays were hybridized with labeled genomic preparations of MG4, 312, and PM1. Genomic DNA was isolated (5) and digested (5 µg per array) with 0.005 U DNase I (Amersham) in 1x One-Phor-All buffer (Amersham, Piscataway, NJ) at 37°C for 5 min, with subsequent inactivation (95°C for 15 min). To the DNA digest were added 4 µl 5x terminal transferase buffer, 1 nmol biotin-N6-ddATP, and 25 U terminal transferase. The sample was incubated at 37°C for 90 min, followed by inactivation at 95°C for 15 min. Hybridization of arrays was conducted in 1x hybridization buffer for 16 h at 45°C, using a Hybriwheel device (NimbleGen). PM1 was used as a reference in the analysis and was hybridized to separate arrays. Duplicate arrays were processed per strain. Arrays were washed with nonstringent wash buffer (6x SSPE [1x SSPE is 0.18 M NaCl, 10 mM NaH2PO4, and 1 mM EDTA {pH 7.7}], 0.01% [vol/vol] Tween 20), followed by two 5-min washes with stringent buffer (100 mM morpholineethanesulfonic acid [MES], 0.1 M NaCl, 0.01% [vol/vol] Tween 20) at 47.5°C. Arrays were stained with Cy3-streptavidin conjugate (Amersham Piscataway, NJ) for 10 min, followed by washing in nonstringent buffer. Signal amplification was achieved by secondary labeling with biotinylated goat anti-streptavidin (Vector Laboratories, Burlingame, CA), washing in nonstringent buffer, and restaining with Cy3-streptavidin. Finally, arrays were washed in nonstringent wash buffer, in 0.5x SSC (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate) two times for 30 s each, and in 70% ethanol for 15 s. Arrays were spun dry by centrifugation. Scanning was conducted at 5-µm resolution with a Genepix 4000b scanner (Axon Instruments, Union City, CA), and NimbleScan software (NimbleGen) was used to obtain pixel intensities. For higher-resolution resequencing of the MG4 and 312 plasmids, arrays were synthesized and hybridized with genomic DNA from each strain and scanned as described above. Single-nucleotide polymorphism (SNP) positions were tested for uniqueness in the genome, using custom algorithms (NimbleGen). The PM1 annotation (http://genome.ornl.gov/microbial/rgel/) was used to generate the output file in SignalMap analysis software (NimbleGen). The predicted SNPs were confirmed by producing and sequencing amplicons, using PCR primers located external to the SNP locations.
Random mutagenesis, mutant characterization, and plasmid curing of PM1.
In order to label the megaplasmid with a selectable marker, random transposon mutagenesis was employed using the mini-transposon derivative pTnMod-SmO, containing the streptomycin (Sm)/spectinomycin adenylyltransferase gene (aadA) and an oriR origin of replication between the inverted repeats (22). Electrocompetent PM1 cells were prepared by being cultured in 0.5x tryptic soy broth at 27°C, with shaking, to log phase. Cells were collected by centrifugation, washed in 10% glycerol four times, and suspended in 10% glycerol to a final volume of 100 µl. A mixture of 50 µl cells and 2 µl pTnMod-SmO DNA (1 µg/µl) was electroporated in 0.1-mm-gap cuvettes, using 1.8 kV, 200
, and 25 µF capacitance settings (22) on a Bio-Rad Gene Pulser electroporator (Bio-Rad, Hercules, CA). Following a 4-h recovery in 0.5x tryptic soy broth at 27°C with shaking, transposon mutants were selected on 0.5x TSA plates with 50 µg/ml Sm. Sm-resistant colonies were present after incubation for 5 to 6 days at 27°C, and stable transposon integration was confirmed by PCR analysis of genomic DNA, using pTnMod-SmO-specific primers.
Using the rapid cloning strategy outlined by Dennis and Zylstra (22), the <SmO> insert location was mapped in several PM1 subclones containing oriR within the transposon. Briefly genomic DNA was extracted, digested with the AvaII restriction endonuclease, self-ligated, and transformed into Escherichia coli TOP10 cells (Invitrogen, Carlsbad, CA). The resulting transformants were selected on LB agar containing 50 µg/ml Sm. Sequencing with primers against the ends of the <SmO> insert was used to determine the exact insert location. One transposon mutant, MP0005, was shown to have the <SmO> insert on the megaplasmid (in MpeB636). MP0005 was subjected to plasmid curing by heat stress as described by Trevors (78). Specifically, strain MP0005 was incubated at 37°C for 6 to 8 h before being plated on 0.5x TSA. Following replica plating on 0.5x TSA, with and without 50 µg/ml Sm, Sm-sensitive colonies were selected, and megaplasmid loss was confirmed by PCR analysis. MTBE and TBA degradation activities by strain MP0005 and a megaplasmid-free strain, MP0007, were determined by gas chromatography analysis as previously described (34).
Nucleotide sequence accession number. The nucleotide sequence annotation is available at the Joint Genome Institute website (http://genome.ornl.gov/microbial/rgel/) and has been deposited in the GenBank/EMBL database under accession numbers CP000555 for the chromosome sequence and CP000556 for the plasmid sequence.
|
|
|---|
![]() View larger version (25K): [in a new window] |
FIG. 1. Circular maps of the chromosome (A) and megaplasmid pPM1 (B) of M. petroleiphilum PM1. Outer rings 1 and 2 show all CDSs, colored by functional category (dark gray, hypothetical proteins; light gray, conserved hypothetical and unknown function; brown, general function prediction; red, replication and repair; green, energy metabolism; blue, carbon and carbohydrate metabolism; cyan, lipid metabolism; magenta, transcription; yellow, translation; orange, amino acid metabolism; pink, metabolism of cofactors and vitamins; light red, purine and pyrimidine metabolism; lavender, signal transduction; sky blue, cellular processes; and pale green, structural RNAs). Genes coding for major metabolic features (green, methylotrophy; red, aromatic hydrocarbon degradation; and light blue, alkane degradation) are shown on rings 3 and 4. Large repeat regions are indicated on rings 5 and 6. The larger repeated regions colored light gray (29-kb repeat) and dark gray (40-kb repeat) are described in the text, while the other colors represent repeated IS elements, ISmp1-8. Ring 7 shows the deviation from the average G+C, and the innermost ring, ring 8, shows the GC skew (GC)/(G+C). The plasmid and chromosome are not drawn to scale.
|
|
View this table: [in a new window] |
TABLE 1. General features of the Methylibium petroleiphilum PM1 genomea
|
Cell motility, secretion, and transport systems. M. petroleiphilum PM1 possesses the genes necessary for flagellar biosynthesis (for one polar flagellum), chemotactic responses, type IV pilus synthesis, and the type II secretion pathway as well as several genes related to the Agrobacterium tumefaciens type IV secretion pathway (Fig. 2; see Table S1 in the supplemental material). Type IV secretion mechanisms are often involved in pathogenesis. However, homologs to only three of the five core type IV secretion proteins (VirB9, -10, and -11, not VirB4 or -7) (6) were identified, so it is unclear at present if PM1 possesses a functional type IV secretion pathway. PM1 likely moves by both flagellum-facilitated swimming and pilus-facilitated twitching motility. Three copies of tra genes on pPM1 suggest that PM1 may be capable of conjugative transfer, a possibility currently under investigation. Thirteen chromosomally and one plasmid-encoded methyl-accepting chemotaxis protein (MCP) allow PM1 to respond to a range of environmental stimuli (see Table S1 in the supplemental material). As in other organisms (83), MCP genes in PM1 are found scattered throughout the genome. Only three MCP genes are found within taxis operons, namely, the pilG-L operon, the cheYA(MCP)W operon, and the flg/flh gene cluster. The apparent mobility of MCPs, together with the fact that six PM1 MCPs appear to be paralogs, complicates functional assignation. Nevertheless, there are five other MCP genes in addition to those already mentioned whose gene environments may offer insight into their possible functions, as follows: two paralogous MCP genes are located immediately downstream of and may be part of the same operon as the genes for the two toluene/benzene monooxygenase pathways; one of the aer-like MCPs is immediately downstream of, and in the same orientation as, the LysR-type activator gene of the ribulose 1,5-bisphophate carboxylase/oxygenase (RuBisCO) operon, which is upstream of the two regulatory genes; one MCP gene may be cotranscribed with a gene showing similarity to the direct oxygen sensor dos of E. coli; and one MCP gene may be cotranscribed with a gene showing low percent similarity to those for bacteriophytochromes and motility sensors. Neighbor-joining analysis of the putative PM1 MCPs against their homologs showed that eight MCPs cluster close to MCP1-4 of E. coli and MCPA-D of Salmonella enterica serovar Typhimurium, two appear to be related to the aerotaxis and energy sensor AER of E. coli, and one is very similar to the twitching motility protein PilJ of Pseudomonas aeruginosa.
![]() View larger version (41K): [in a new window] |
FIG. 2. Schematic diagram of an M. petroleiphilum PM1 cell showing structural features and cellular processes, including predicted methylotrophy and fuel hydrocarbon degradation pathways. Two adjacent arrows imply multiple steps. Abbreviations: MTBE, methyl tert-butyl ether; HMTBE, hydroxymethyl tert-butyl ether; TBF, tert-butyl formate; TBA, tert-butyl alcohol; HIBA, hydroxyisobutyric acid. Systems involved in metal resistance and transport, secretion, motility, chemotaxis, and electron transport are shown within or associated with the cell membrane. Cofactors are labeled green, and electron donors are labeled red in the cytosol. The figure is not drawn to scale.
|
Repeated elements. The genome has a number of complex repetitive elements, including eight families of insertion sequences (ISmp1-8) (up to 12 copies) and two large genomic segments (29 and 40 kb) flanked by IS elements that have undergone what appear to be recent duplications (Fig. 1). The two replicons do not equally share the repeated insertion elements; five of the eight families are located only on the chromosome and one family is found strictly on the plasmid. The distribution patterns of the IS elements lend support to the dissimilar phylogenetic distributions of best KEGG hits among sequenced genomes and strengthen the notion of the megaplasmid's recent acquisition.
Parallel copies of ISmp8 flank two tandem copies of a 29-kb repeat, each consisting of two operons involved in phosphonate and cobalamin metabolism. The phosphonate operons (phnFDC-htxFGHIJKLMN) include genes for putative C-P lyase subunits which are 54 to 83% similar to those of Pseudomonas stutzeri WM88 (81) (see Table S2 in the supplemental material). The Htx and Phn C-P lyases support growth on methylphosphonate and additional alkylphosphonates, respectively; growth on these substrates is not yet known for PM1. Also contained in the repeat are cobalamin (vitamin B12) synthesis genes encoding the conversion of uroporphyrinogen III to cobinamide and the synthesis of dimethylbenzimidazole in the aerobic pathway for cobalamin biosynthesis (MpeB437-453 and MpeB472-488) (see Table S1 in the supplemental material). Downstream of the tandem repeat are the remaining genes (MpeB509-522) for the covalent linkage of dimethylbenzimidazole, cobinamide, and a phosphoryl group to complete the cobalamin synthesis pathway. Genes coding for the anaerobic pathway of cobalamin biosynthesis (cbi genes) are also present in the cob clusters, but a complete pathway is lacking. PM1 also lacks a cobG gene, encoding the monooxygenase that converts precorrin-3 to precorrin-4 in the aerobic pathway, but one or more of the multiple copies of cbiG (MpeB479, -480, -444, and -445) may code for a functional enzyme that performs this reaction without oxygen (60). Cobalt and cobalamin (vitamin B12) have been shown to enhance PM1's ability to grow on and degrade MTBE and its primary metabolite, TBA (Hristova, unpublished data), so it is not surprising that multiple copies of genes involved in cobalamin synthesis are present in PM1, including tandem repeats of cob and cbi genes. Recently, Rohwerder et al. (61) showed that cobalamin synthesis affected the growth rate on the MTBE metabolites TBA and 2-hydroxyisobutyrate (2-HIBA) for a beta-proteobacterial MTBE-degrading strain that was phylogenetically similar to PM1 (95.6% identical based on the 16S rRNA gene sequence). In this strain, cobalt or cobalamin was necessary for the activity of an enzyme, isobutyryl-CoA mutase, involved in metabolism of 2-HIBA (61), and 99% identical homologs to this two-component mutase are present in PM1 (MpeB538/541). As mentioned, a relatively large percentage of predicted proteins whose genes are carried on the plasmid (7.0%) belong to the COG category for coenzyme transport and metabolism. A cluster of ethanolamine utilization (eut) genes (MpeB0499-502) found between the tandem repeats and the third cobalamin cluster on the plasmid encode putative proteins which are 48 to 85% similar to EutJEMN from the cobalamin-dependent ethanolamine utilization pathway of S. enterica serovar Typhimurium (49). The latter two proteins are homologs of carboxysome shell proteins (CcmKL). While S. enterica serovar Typhimurium also contains the ethanolamine lyase subunits and regulator (EutBC and EutR) in its eut operon, in PM1 these genes (MpeA2417-8 and MpeA2415) are on the chromosome.
The largest (40 kb) repeated element is present on both the plasmid and the chromosome and encodes a putative PinR-like site-specific recombinase, a replicative DNA helicase, two putative SpoJ-like transcriptional regulators or plasmid partitioning proteins, a SpoIIIE-like DNA translocase, a tellurite resistance protein, and many hypothetical products. The presence of the repeat on both replicons suggests a recent duplicative transfer event. Although the types of genes found in this region suggest a plasmid origin, both copies interrupt similar but nonidentical copies of dcd genes (for dCTP deaminase), and the direction of duplicative transfer remains to be proven. Since this duplication, the 40-kb repeat on pPM1 has been interrupted by a transposase between genes MpeB0184 and MpeB0187.
Heavy metal tolerance and metal homeostasis. One interesting outcome of the genome analysis is evidence of PM1's potential resistance to heavy metals, suggesting promise in using the organism to treat sites containing mixed wastes. Arsenic extrusion in PM1 is probably mediated by arsRBC, present in two copies on the chromosome that are 58 to 81% similar to each other (MpeA1581-1584, arsHBCR; MpeA2343-4 and MpeA2347, arsCB and arsR, respectively). While some other bacteria have five-gene operons (arsRDABC) and use the ArsAB pump, PM1 probably extrudes arsenite by a carrier protein, ArsB, energized by the membrane potential (62), although resistance to arsenic oxyanions needs to be evaluated. ArsC is an arsenic reductase responsible for the transformation of As(V) into As(III), and ArsR is a transcriptional repressor that responds to As(III) and Sb(III) (62). The function of the fourth gene product, ArsH, found in several bacteria (Yersinia enterocolitica, Acidothiobacillus ferrooxidans, Pseudomonas aeruginosa, and Pseudomonas putida KT2440), still remains unclear (10, 63). Three chromosomal copies of chrA (MpeA2204, MpeA2205, and MpeA2526), belonging to the CHR family of transporter genes, may mediate chromate resistance in PM1. One copy of the ChrA gene (MpeA2204) is 63 and 61% similar to its homologs in Dechloromonas aromatica RCB and P. putida KT2440, respectively, although a homologous chromate reductase (ChrR) (42) was not evident in PM1.
Genome analysis revealed 15 copper resistance genes in a large cluster (
14.4 kb, at positions 1760297 to 1775675) on the PM1 chromosome with structural analogy to the plasmid-mediated (pMOL30) copper resistance cluster copOAIPRSFG in Ralstonia metallidurans (54) (see Table S2 in the supplemental material). The copOAIP-copRS cluster in PM1 is likely involved in the efflux of periplasmic copper (analogous to the cop systems of Pseudomonas syringae, R. metallidurans, and R. solanacearum and the cos system of E. coli) (13), whereas the efflux of cytoplasmic copper is mediated by a P1-ATPase, CopF1F2 (analogous to R. metallidurans CopF) (54). The genome of PM1 also encodes a putative chemiosmotic antiporter efflux system similar to CzcCBA of R. metallidurans, conferring resistance to Cd, Zn, and Co (54). In addition to copF, there are two other genes encoding putative metal-transporting P1-type ATPases, namely, MpeA2479 and MpeA3535. Additional proteins potentially involved in metal transport include NikBCDE for nickel (MpeA3117-3120), CbiOQMK for cobalt and nickel (MpeA2799-2802), and ModABC (MpeA3707, MpeA3714, and MpeA3715) for molybdenum uptake (see Table S2 in the supplemental material).
Ferric iron has also been shown to enhance PM1's ability to grow on and degrade MTBE and TBA (Hristova, unpublished data), so it is likely that active transport of iron is of particular importance. As with other gram-negative bacteria, PM1 acquires its iron supply via Fe3+ siderophores. fep genes, which function in the synthesis of polypeptides required for uptake of ferric enterobactin, were identified inside each of the four cob operons of PM1 (see Table S1 in the supplemental material). Polaromonas sp., Rhodoferax ferrireducens, and Methylobacillus flagellatus all have iron transport genes either within or in close proximity to their cob operons (data not shown). Minimally, the FepABC proteins are required for ferric enterobactin uptake. MpeA2292 and MpeA2605 have been annotated as fepA, coding for an outer membrane receptor for an iron siderophore, but it is possible that btuB genes located near the febBDC genes are also involved in iron assimilation. The TonB-dependent energy transduction complex (tonB, exbB, or tolQ and exbB) (see Table S1 in the supplemental material) encoded on the chromosome likely provides the mechanism for active transport of iron siderophores and cobalamin across the outer membrane (8, 37). The PM1 genome encodes about 39 putative proteins involved in iron transport and homeostasis, which implies the importance of iron in its physiology.
Methylotrophy. Methylotrophic metabolism of PM1 is of great interest because formaldehyde and formate are common intermediates of both methanol and MTBE or TBA oxidation by PM1 and other degraders (9, 59). M. petroleiphilum PM1 is capable of aerobic growth on methanol, formate, and succinate. Unlike other methylotrophic beta-proteobacteria, PM1 grows on MTBE, toluene, benzene, ethylbenzene, and dihydroxybenzoates as sole carbon sources (57, 59). PM1 possesses genes for the serine cycle and methylotrophy scattered in several different clusters on its chromosome (Fig. 1; Table 2). The strain does not grow on methylamine (Hristova, unpublished data), lacks a gene encoding the methylamine dehydrogenase large subunit, and likely lacks methylamine dehydrogenase activity. Despite the ability of PM1 to grow on methanol, its genome lacks true homologs of mxaF and mxaI, known genes coding for the methanol dehydrogenase (MeDH) large and small subunits, respectively, present in several methylotrophs known to date. PM1 contains an MeDH-like cluster, xoxF-J (MpeA3393-5), that is present in Methylobacterium extorquens AM1 (17), which also contains the true mxaF cluster. Comparative sequence analysis of the product of the MpeA3393 gene revealed high similarity to MxaF (74% to that of M. extorquens AM1) and the XoxF homologs present in several nonmethylotrophs (77% to that of Burkholderia fungorum). Based on phylogenetic analysis, MpeA3393 clusters with the MxaF homologs of unknown function from other methylotrophic and nonmethylotrophic Rhizobium and Burkholderia spp., while true MeDH large subunits (MxaF) cluster together and are distinct from the MxaF homologs (Fig. 3).
|
View this table: [in a new window] |
TABLE 2. M. petroleiphilum PM1 genes putatively involved in metabolism of methanol
|
![]() View larger version (26K): [in a new window] |
FIG. 3. Phylogenetic tree of MpeA3393, which putatively codes for the large subunit of methanol dehydrogenase. The protein product from MpeA3393 is most similar to glucose dehydrogenase (GluDH) from B. fungorum (81% similar; GenBank accession no. ZP_00283396) and the methanol dehydrogenase large subunit (Mdh large) from M. capsulatus Bath (75% similar; accession no. AAU90462). Protein names are shown in square brackets.
|
Three different formate dehydrogenase genes are present in the PM1 genome, with homologs to genes in M. extorquens and M. capsulatus Bath. The functions of the tungsten-dependent formate dehydrogenase gene fdh1 (MpeA0337-339), the NAD-linked formate dehydrogenase gene fdh2 (MpeA3708-12), and the cytochrome-linked formate dehydrogenase gene fdh3 (MpeA1170-71 and -1173) for energy generation during growth on C1 substrates or for MTBE oxidation need to be explored further. MpeA3377, coding for a putative ABC-type tungstate transport system permease, links gene clusters fdh1 and fdh2 (Table 2). The fdh2 genes in PM1 have the same gene arrangement as and significant sequence identity (52 to 81%) to the NAD-dependent formate dehydrogenase gene cluster, fdsGBACD, of Ralstonia eutropha (58). Pathways involved in metabolism and detoxification of formaldehyde, a central intermediate of both methanol and MTBE degradation, may also function in MTBE metabolism. PM1 has two pathways for formaldehyde oxidation to CO2, a tetrahydromethenopterin (H4MPT)-linked metabolic module that includes an archaea-like gene cluster and an H4F-linked metabolic module. Recently, phylogenetic analysis of a subset of bacterial and archaeal H4MPT-linked C1 transfer genes placed PM1 sequences with those of other described beta-proteobacteria (44).
Fuel hydrocarbon degradation pathways. PM1 contains an operon (MpeA0814-0821) likely encoding enzymes for conversion of benzene to phenol (and catechol) and toluene to methylphenol (and methylcatechol) that is 62 to 74% similar to the benzene monooxygenase pathway in P. aeruginosa JI104 (BmoA-D1) (48) and 50 to 71% similar to the toluene para-monooxygenase (TpMO) pathway in Ralstonia pickettii PKO1 (TbuA1UBVA2CX) (76) (Table 3). A second operon (MpeA2539-2547) is 55 to 74% identical to the first operon, but it likely does not yield a functional monooxygenase since the TbuA1 homolog is interrupted by a transposon insertion and the TbuC homolog is a pseudogene. Both operons have two-component response regulator-sensor histidine kinases which are located upstream and divergently transcribed (MpeA0811-812 and MpeA2536-2537), although MpeA2537 may be truncated due to the transposon insertion. MpeA0821 encodes a putative TbuX protein (65% similar to that in PKO1), an outer membrane protein regulated by TbuT and involved in toluene uptake (43). The benzene monooxygenase pathway has been implicated in benzene and toluene degradation (48), as has the TpMO pathway (76). In addition to benzene and toluene, PM1 has been shown to degrade o-xylene (20), although the biochemical pathway used has not been elucidated to date. It is likely that m- and p-xylene can also be metabolized via the toluene monooxygenase (TMO) pathway of PM1, as described for PKO1 and other bacteria (29).
|
View this table: [in a new window] |
TABLE 3. M. petroleiphilum PM1 genes putatively involved in metabolism of aromatic hydrocarbons denoting the predicted role of the gene product and percent similarity with well-characterized homologs
|
PM1 has nine CDSs encoding putative proteins with various similarities to cyclohexanone monooxygenases (CHMOs), sometimes referred to as Baeyer-Villiger-type MOs (MpeB0579, -B607, -B610, -A393, -A898, -A1038, -A1351, -A2885, and -A2915). Their protein products may play a role in hydroxylation of either alicyclic, aliphatic, or aryl ketones to form a corresponding ester, which can easily be hydrolyzed. Alicyclic hydrocarbons represent up to 12% (wt/wt) of total hydrocarbons in petroleum mixtures (according to the American Petroleum Institute). Aryl ketones such as acetophenone can be produced directly from atmospheric breakdown of ethylbenzene (a major petroleum component) or following abiotic conversion of ethylbenzene to ethylphenol (3) and subsequent biological conversion to the ketone. The putative CHMO genes are scattered across the genome and are not present in operons with other genes coding for subsequent metabolism after the MO reaction (i.e., esterases and alcohol and aldehyde dehydrogenases). The CHMOs have a narrow substrate range, possibly explaining the number of different flavoprotein MOs in PM1 with various levels of similarity to representatives from this class; the putative CHMOs in PM1 are more similar to phenylacetone MO (46 to 67%) (53) than to 4-hydroxyacetophenone MO (43 to 53%) (45). In PM1, the nine Baeyer-Villiger MOs have a putative NADP+-binding site that is 72 to 88% similar to the proposed site in a CHMO from Acinetobacter sp. strain NCIMB 9871 (16).
An alkane monooxygenase pathway on pPM1 may facilitate PM1's growth on n-alkanes. In addition, alkane monooxygenase (hydroxylase) has been proposed to play a role in cometabolic MTBE hydroxylation, since acetylene, an inactivator of short-chain alkane monooxygenase, was shown to inhibit MTBE degradation (68). In PM1, the hydroxylase subunit, AlkB (MpeB0606), is 69% and 66% similar to that of Alcanivorax borkumensis AP1 (72) and P. putida PGo1 (contained on the OCT plasmid) (79), respectively, and contains all eight of the conserved His residues observed in other integral membrane binuclear-iron hydrocarbon monooxygenases (33). Also present are two rubredoxin (Rd) genes (MpeB0603 and MpeB0602), whose products are 76 and 78% similar to rubredoxin 3 and 4, respectively, in Gordonia sp. strain TF6 (31). The Rd encoded by MpeB0603 is an AlkG1-type Rd, whereas that encoded by MpeB0602 is an AlkG2-type Rd, based on the CXXCG motif criterion described by van Beilen et al. (80). Only AlkG2-type Rds were shown to be functional in electron transfer from the rubredoxin reductase to alkane hydroxylase. In addition, MpeB0601 codes for an ATP-dependent transcriptional regulator which is 38% similar to AlkS from A. borkumensis SK2 (35). Separated from the putative alkS gene by three hypothetical genes is a rubredoxin reductase gene, alkT (MpeB0597), whose protein product is 49% similar to that in Gordonia sp. strain TF6. PM1 does not appear to possess any long-chain alkane (>C13) oxidation pathways, such as an alkane dioxygenase (64), P-450 monooxygenase (2), or two alkane hydroxylase complexes (AlkMa and AlkMb), similar to Acinetobacter sp. strain M-1 (75), although PM1's single AlkB protein is 54% similar to both AlkMa and AlkMb. The gene organization of the alk operon in PM1 is somewhat similar to that of Gordonia sp. strain TF6 (alkB2G1G2T), except that a transposase gene (MpeB605) and putative esterase gene (MpeB604) are between alkB and alkG1G2, and as mentioned, a putative alkS gene and three hypothetical genes (MpeB600-598) are located between alkG1G2 and alkT in PM1. Homologs to alkHJKL from P. putida GPo1, coding for aldehyde dehydrogenase, alcohol dehydrogenase, acyl-CoA synthetase, and an outer membrane protein (79), respectively, were not present on the plasmid, although the PM1 chromosome contains homologs to the AlkH (MpeA2324, 47% similar), AlkJ (MpeA3803, 58% similar), AlkK (MpeA1769, 71% similar) and AlkL (MpeA3010, 49% similar) genes.
In addition, the PM1 chromosome contains a putative propane monooxygenase pathway (MpeA950-953) whose predicted proteins are 41 to 64% identical to PrmABCD in Gordonia sp. strain TY-5, including the large hydroxylase subunit, the NADH-dependent acceptor oxidoreductase, the small hydroxylase subunit, and the coupling protein, respectively. The Prm complex in strain TY-5 was shown to catalyze the subterminal oxidation of propane, yielding 2-propanol (50), while propane oxidation in PM1 is currently under investigation. As for PrmA in Gordonia TY-5, a pair of conserved Glu-X-X-His sequences are present in the putative PrmA protein of PM1, at residues 138 to 141 and 237 to 240. The presence of these sequences is consistent with the case for other monooxygenases in the binuclear-iron oxygenase family, including soluble methane monooxygenases (23, 71), suggesting that PrmA in PM1 may catalyze the hydroxylation of propane. Like the operon in strain TY-5, a chaperone similar to GroEL is encoded adjacent to the prm cluster in PM1 (MpeA0954). Finally, PM1 has homologs to strain TY-5 alcohol dehydrogenase genes, namely, adh-1 (MpeA936) and adh-3 (MpeA599), that are 72 and 83% similar, respectively, and may facilitate 2-propanol degradation. The putative monooxygenase genes in PM1 are summarized in Table S3 in the supplemental material, including the methanesulfonate monooxygenase gene msmA and the alkanesulfonate monooxygenase gene ssuD, which are part of the msmABDCEFGHG and ssuAADCB operons, respectively. PM1 may not utilize methanesulfonate, since its msmA gene lacks the sequence encoding CXH-X26-CXXH unique to methanesulfonate utilizers (7). In general, PM1 possesses several homologous genes with other soil bacteria, including Gordonia, Alcinovorax, and Pseudomonas spp., capable of biodegradation of petroleum hydrocarbons as well as xenobiotic and recalcitrant compounds such as phthalates.
MTBE biodegradation. Although MTBE is a recent anthropogenic contaminant (released within the last 15 years), various microorganisms can utilize the compound for carbon and energy under aerobic conditions (30, 61, 65, 74). M. petroleiphilum PM1 is the best characterized of the few bacterial pure cultures reported to grow on and completely degrade MTBE and its daughter product, TBA (20, 30, 36, 65). The genetic basis for MTBE and TBA conversion is not known, although different classes of monooxygenases have been proposed to play a role in metabolism or cometabolism of these compounds (30, 36, 51, 68, 74), including P-450 monooxygenase and alkane monooxygenase (hydroxylase) systems, the latter of which was shown to play a role in cometabolic degradation of MTBE by P. putida GPo1 (69) and possibly also by Pseudomonas mendocina KR-1 (70). A known inducer of alkane hydroxylase, dicyclopropylketone, was also shown to induce MTBE conversion to TBA in GPo1 (69). As reported above, an alkane MO (AlkB) system was detected in the PM1 genome on the megaplasmid. The AlkB system in PM1 is likely involved in MTBE hydroxylation, based on its similarity to other AlkB proteins in organisms shown to be involved in MTBE degradation. Whereas the Ks values for MTBE in n-alkane-grown GPo1 were reported to be quite high (20 to 40 mM), the apparent half-saturation constant for MTBE by PM1 was 88 µM, which is in the range of Ks values for MTBE by butane-degrading bacteria (51). Unlike GPo1 and KR-1, PM1 further degrades TBA, ultimately producing CO2 and biomass. The putative AlkB protein in PM1 is proposed to oxidize only MTBE, not TBA, based on kinetic experiments with MTBE- and TBA-grown cells (19). Two separate enzyme systems were also suggested for MTBE and TBA degradation by Hydrogenophaga flava ENV735 (36). Because of its potential role in MTBE metabolism, the coding region for AlkB is the focus of current gene knockout studies. Biodegradation of a similar molecule, ethyl-tert butyl ether (ETBE), occurs via a cytochrome P-450 pathway in Rhodococcus ruber IFP2001 (15). However, homologs of protein complexes involved in ETBE degradation from R. ruber were not found in PM1; like GPo1 (69), PM1 has not been shown to degrade ETBE.
Many pollutant degradation genes are located on bacterial catabolic plasmids. Significantly, the two strains, MG4 and 312, that are capable of MTBE degradation have nearly identical plasmids to that of PM1 (ca. 99% identical), as determined by comparative genome sequencing analysis. The MG4 and 312 plasmids showed only five and four SNPs, respectively, relative to PM1 (Table 4). The MG4 and 312 plasmids also lacked transposase genes (three copies of Tra5 and transposase 8 genes and one copy of a DDE domain transposase gene) that are present on the PM1 plasmid and a 1.2-kb deletion putatively containing an esterase/lipase gene (MpeB0604) and a DDE domain transposase gene (MpeB0605) (Table 3). The promoter and coding regions for alkB (MpeB0606) did not appear to be affected by this deletion, since it is significantly upstream, although there was a SNP mapped within alkB of MG4 and 312 resulting in a putative amino acid change. As mentioned, two other PM1 plasmid-carried genes, MpeB541 and MpeB0538, code for putative large and small subunits of isobutyryl-CoA mutase, respectively. The plasmids of strains MG4 and 312 also contained identical copies of MpeB541 and MpeB0538. These gene products were shown to have 99% identical homologs in Ideonella sp. strain L108 which were predicted to allow conversion of 2-HIBA to 3-hydroxybutyrate in the presence of CoA and ATP (61). It is not known whether these mutase genes are contained on a megaplasmid in L108, although horizontal gene transfer is often evoked when such high similarities in gene sequences between bacteria are observed. In addition to alkB, the PM1 plasmid (as well as the MG4 and 312 plasmids) contains a gene coding for 3-hydroxybutyryl-CoA dehydrogenase (MpeB0547), putatively involved in conversion of 3-hydroxybutyryl-CoA, a proposed MTBE metabolite (61), to acetoacetyl-CoA.
|
View this table: [in a new window] |
TABLE 4. Genomic differences between the plasmid of PM1 and those from isolates MG4 and 312
|
![]() View larger version (16K): [in a new window] |
FIG. 4. MTBE and TBA degradation by M. petroleiphilum PM1. The graph shows data for the parent strain MP0005, which carries the <SmO> marker on the megaplasmid (, MTBE; , TBA), and the megaplasmid-free mutant MP0007 ( , MTBE; , TBA).
|
We acknowledge the University of California Office of the President for funding through the Campus-Laboratory Collaboration Program as well as the LLNL Laboratory-Directed Research and Development Program. This research was also supported by grant number 5 P42 ES04699-16 from the National Institute of Environmental Health Sciences (NIEHS), NIH, with funding provided by the U.S. EPA. This work was performed under the auspices of the U.S. Department of Energy by the University of California, Lawrence Livermore National Laboratory, under contract no. W-7405-Eng-48.
The contents of this study are solely the responsibility of the authors and do not necessarily represent the official views of the NIEHS, NIH, or EPA.
Published ahead of print on 8 December 2006. ![]()
Supplemental material for this article may be found at http://jb.asm.org/. ![]()
S.R.K. and A.Y.C. gave equal contributions to this study. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»