Previous Article | Next Article ![]()
Journal of Bacteriology, June 2005, p. 4050-4063, Vol. 187, No. 12
0021-9193/05/$08.00+0 doi:10.1128/JB.187.12.4050-4063.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Department of Microbiology and Immunology, University of British Columbia, Vancouver, BC, V6T 1Z3, Canada
Received 24 September 2004/ Accepted 15 March 2005
|
|
|---|
|
|
|---|
Rhodococcus sp. strain RHA1 was isolated from lindane-contaminated soil (69) for its exceptional ability to aerobically degrade polychlorinated biphenyls (PCBs), a class of toxic and persistent pollutants. As in other aerobic PCB-degrading bacteria, these pollutants are cometabolized by the bph pathway, which is responsible for the aerobic degradation of biphenyl. The bph pathway consists of four enzymatic activities which act sequentially to transform biphenyl to benzoate and 2-hydroxypenta-2,4-dienoate as recently reviewed by Furukawa (30). For each of these four steps, RHA1 appears to possess multiple isozymes, including at least three bph-type ring-hydroxylating dioxygenases (51) and at least seven different bph-type ring cleavage enzymes (67). While most of the genes of the upper bph pathway are located on two of three large linear plasmids (72), pRHL1 (1,100 kb) and pRHL2 (450 kb), genes encoding related isozymes are distributed throughout the 9.7-Mb genome. It is unclear which of these isozymes is involved in the catabolism of biphenyl or related compounds, how these different activities are regulated, and whether this apparent redundancy is a general characteristic of catabolic pathways in rhodococci.
Rhodococcus sp. strain RHA1 utilizes benzoate and phthalate as sole sources of carbon and energy. The catabolism of these compounds is initiated by ring-hydroxylating oxygenases encoded by the ben (50) and pad (Fukuda, personal communication) genes, respectively. In gram-negative bacteria, the initial phthalate dioxygenase mediates 4,5-dihydroxylation (5, 14). In contrast, the phthalate dioxygenase of gram-positive bacteria mediates 3,4-dihydroxylation (26). The involvement of the ben and pad genes implies that benzoate and phthalate are catabolized via catechol and protocatechuate, respectively, in RHA1.
In other bacteria, both protocatechuate and catechol can be further degraded via either meta- or ortho-cleavage pathways. The ortho-cleavage pathway, commonly known as the ß-ketoadipate pathway, is widely conserved among diverse soil bacteria, and separate branches catabolize catechol and protocatechuate (37). Permutations of the pathway occur in different bacterial groups with respect to enzyme distribution (isozymes and points of convergence), regulation, and gene organization (37). For example, at least three arrangements of the branches have been reported in different genera: in Ralstonia eutropha (formerly Alcaligenes eutrophus), the two branches converge at ß-ketoadipate (44); in pseudomonads, they converge at an upstream enol-lactone intermediate (42); and in Acinetobacter sp. strain ADP1 (formerly Acinetobacter calcoaceticus), the branches do not converge at all (23).
As part of an effort to better understand the metabolism and physiology of rhodococci and related strains, we are sequencing and annotating the genome of Rhodococcus sp. strain RHA1 and studying this organism using a number of functional genomic approaches. One important methodology that has not been well developed for rhodococci is targeted gene deletion. The development of reliable methodologies has been hampered in part by the genetic instability and nonhomologous recombination typical of actinomycetes (53). Insertion mutagenesis has been used to disrupt genes by single crossover followed by genetic complementation (50, 57, 68, 78). However, this approach often generates polar effects on downstream genes. Using a counterselectable marker, Van der Geize et al. (82) constructed unmarked gene deletions in Rhodococcus erythropolis SQ1. Recently, Gust et al. described a strategy using
Red-mediated double-crossover recombination to create in-frame, nonpolar gene deletions in Streptomyces coelicolor A3(2) (35). In this approach, the gene is first replaced in a cosmid carrying the genomic DNA of interest. The mutagenized cosmid is then introduced into the streptomycete to effect allelic exchange. The method enables the deletion of entire gene clusters as well as single and multiple deletions. In principle, a similar approach could be used in Rhodococcus.
The current study describes an investigation of the catabolism of phthalate and benzoate in Rhodococcus sp. strain RHA1. Analyses of the nearly completed genome assembly revealed the presence of several putative operons that are involved in the catabolism of these aromatic carboxylic acids. Proteomic analyses of cells grown on pyruvate, phthalate, and benzoate identified pathway genes and enzymes. The involvement of key enzymes was confirmed by activity assays and gene disruption. The latter was accomplished by adapting a gene replacement methodology developed for Streptomyces and should be of general utility to study Rhodococcus and other actinomycetes. To investigate aspects of pathway and gene organization that might be unique to actinomycetes, the deduced pathways in RHA1 were compared with those reported for other bacteria and what could be deduced from genomic sequence data. These studies provide insights into the origin of the catabolic diversity of rhodococci.
|
|
|---|
Strains, media, and plasmids. The bacterial strains and plasmids used in this study are listed in Table 1. Rhodococcus sp. strains RHA1 and RHA1_005 were grown at 30°C on W medium (57) supplemented with 20 mM of an appropriate carbon source (benzoate, phthalate, or pyruvate). For proteomic studies and enzyme assays, rhodococcal strains were typically grown as 500-ml cultures in 2-liter Erlenmeyer flasks shaken at 200 rpm to mid-log phase as determined by optical density at 600 nm. Cells were harvested by centrifugation for 10 min at 16,887 x g at 25°C. The cell pellets were flash frozen in liquid nitrogen and stored at 80°C.
|
View this table: [in a new window] |
TABLE 1. Strains and plasmids used in this study
|
|
View this table: [in a new window] |
TABLE 2. PCR primers
|
For enzyme assays, cell-free protein extracts were prepared in essentially the same manner except that the cell pellets were washed twice using 50 mM Tris-HCl, pH 7.0, and the extracts were used immediately. Protein concentration was determined using the 2D Quant kit (Amersham Biosciences).
Two-dimensional gel electrophoresis. Protein two-dimensional gel electrophoresis was performed as previously described (31, 32), with the following modifications. The first-dimensional separation was carried out using nonlinear immobilized pH gradient (IPG) strips (24 cm, pH 3 to 7). The strips were rehydrated with the protein sample (90 µg of protein extract) in 400 µl rehydration solution (10 M urea, 2 M thiourea, 30 mM dithiothreitol, 3% CHAPS, Pharmalyte pH 3 to 10). Isoelectric focusing in the IPG strips was carried out for a total of 73.5 kVh at 20°C under mineral oil using ETTAN IPGphor (Amersham Biosciences). The IPG strips were then equilibrated and run into 12% sodium dodecyl sulfate (SDS)-polyacrylamide gel electrophoresis (PAGE) gels (24 by 20 cm) using the ETTAN DALTtwelve System (Amersham Biosciences). Broad-range molecular mass markers (Invitrogen) were run on each side of the gel. Protein was detected using silver nitrate (for screening purposes) or Sypro Ruby (for quantitative analysis). Silver- and Sypro Ruby-stained gels were imaged using a flat-bed image scanner and a variable mode imager Typhoon 9400 (excitation 488 nm, emission 610 nm; Amersham Biosciences), respectively.
Analysis of two-dimensional gels. Two-dimensional gels were differentially analyzed using Progenesis Workstation software (Nonlinear Dynamics, Durham, NC). Accordingly, each image was processed as follows: (i) sharp spikes were removed from the image as noise; (ii) background values were calculated using a mathematical surface model and subtracted; (iii) spots were detected; (iv) a signal intensity was assigned to each spot; and (v) the signal intensity of each spot was normalized against the total signal intensity of the gel. Processed gel images were matched using the combined warping and matching algorithm. Finally, the signal intensity of each spot was averaged over gels obtained from three biological replicates. In the current experiments, only spots with a minimum normalized volume of 0.002 or greater were analyzed. Molecular mass and isoelectric point values were assigned using the calibration standards. Protein spots whose intensities increased or decreased at least twofold versus the control (pyruvate-grown cells) were recorded as more or less abundant, respectively.
Protein identification. Proteins were identified based on peptide mass and/or peptide fragment mass fingerprint analyses (mass spectrometry [MS] and/or MS/MS). Spots of interest were excised from Sypro ruby-stained gels and digested in-gel using trypsin (49). Mass spectrometry analyses were performed using either a Voyager DESTR matrix-assisted laser desorption ionization-time of flight (MALDI-TOF) or a Sciex linear ion trap quadrupole liquid chromatography (LC)-MS/MS (Applied Biosystems). Proteins were identified using the MASCOT search engine (www.matrixscience.com) and a database generated by in silico digestion of the total RHA1 proteome predicted from build 34.03 of the genome assembly (http://www.bcgsc.bc.ca/cgi-bin/rhodococcus/blast_rha1.pl).
Searches were performed without constraining protein molecular mass or isoelectric point and allowing for the following modifications: carbamidomethylation of cysteine, partial oxidation of methionine residues, and up to one missed trypsin cleavage. A protein was considered identified if the hit fulfilled four criteria: the hit was statistically significant (a MASCOT search score above 55 for the RHA1 database), the number of matched peptides was four or higher, the protein sequence coverage was above 20%, and predicted molecular mass and isoelectric point values were consistent with the experimentally determined ones. Theoretical molecular masses and isoelectric points of the proteins of interest were calculated using EXPASY tools (8).
Enzyme assays.
Enzyme assays were performed using a Varian Cary IE UV-visible spectrophotometer equipped with a thermostatted cuvette holder. Protocatechuate 3,4-dioxygenase activity was determined by monitoring the transformation of protocatechuate to ß-carboxy-cis,cis-muconic acid at 290 nm (
= 2.3 mM1 cm1 (41)) in an assay mixture containing 50 mM Tris-HCl buffer (pH 8.8) and 160 µM protocatechuate. Catechol 1,2-dioxygenase activity was measured by monitoring the formation of cis,cis-muconate at 260 nm (
= 16.8 mM1 cm1 (22)) in an assay mixture containing 50 mM Tris-HCl buffer (pH 8.8) and 200 µM catechol. The activity of ß-ketoadipate:succinyl-coenzyme A (CoA) transferase was measured by monitoring the increase in absorbance at 305 nm (
= 16.3 mM1 cm1) using an assay mixture containing 35 µM Tris-HCl buffer (pH 8.0), 25 µM MgCl2, 3.5 µM ß-ketoadipate, and 0.15 µM succinyl-CoA (46). In all cases, 1 unit of enzyme activity was defined as the amount of enzyme required to produce 1 µmol of product per minute at 25°C.
Gene replacement on fosmids.
Fosmid RF00111bAO4 (Fig. 1b ), created as part of a fosmid library of the RHA1 genome (84), contains 46 kb of RHA1 genomic DNA, including pcaL. The parent vector, EpiFOS (Epicentre), carries a Cmr gene, cat, that is ineffective in RHA1. Accordingly, the
RED-based methodology was employed to replace the cat gene of RF00111bAO4 with an Hmr gene, hyg, and pcaL with an Aprar cassette.
![]() View larger version (13K): [in a new window] |
FIG. 1. Outline of the targeted gene disruption strategy. The method used to disrupt pcaL was adapted from the -Red-based system (35) and involves three major steps: (1) Replacement of the resistance gene in RF00111bAO4. The fosmid RF00111bAO4 (b) carries the chloramphenicol resistance gene (cat) and a 46-kb insert of RHA1 genomic DNA containing pcaL. An amplicon was generated containing the hygromycin resistance gene (hyg) flanked by the same 39-nucleotide sequences that flank the cat gene. The cat gene was replaced by electroporating the hyg amplicon into E. coli BW25113 carrying pKD46 and RF00111bAO4 and inducing the -Red system.(2) Replacement of pcaL on the fosmid. An amplicon (c) containing the apramycin resistance gene [aac(3)IV] flanked by the same 39-nucleotide sequences that flank pcaL. Allelic replacement in the fosmid was achieved as described in step 1, yielding RFMD2 (d). (3) Replacement of the pcaL gene on RHA1 chromosome. RFMD2, carrying the hyg gene and the disrupted pca cluster, was conjugated into RHA1 cells. Allelic exchange between the fosmid and the chromosome (e) resulted in pcaL replacement with Aprar cassette (f). The latter was selected by screening for apramycin resistance and hygromycin sensitivity and verified by PCR. The Aprar cassette contains oriT and FRT sites. oriT allows conjugal transfer into RHA1. The FRT sites, allowing FLP recombinase-mediated elimination of the disruption cassette, were not used in the current study. Additional details are provided in Materials and Methods.
|
Amplification was performed in a 50-µl reaction with 100 ng of template, 5% dimethyl sulfoxide, 50 pmol each primer, and 200 mM deoxynucleoside triphosphates as described (35). Electrocompetent cells of E. coli BW25113, containing the
RED recombination plasmid pKD46 (Ampr), were transformed with the fosmid and selected on ampicillin and chloramphenicol at 30°C to prevent the loss of pKD46. Transformants were used to prepare electrocompetent cells grown at 30°C in SOB (68a) containing ampicillin, chloramphenicol, 20 mM MgSO4 and 10 mM L-arabinose. The latter induces the RED genes. Competent cells were electrotransformed with the PCR-extended Hmr cassette and selected on LB containing hygromycin and ampicillin at 37°C to induce the loss of pKD46. Successful transformants (containing a copy of RFMD1, the Hmr fosmid) were checked by PCR (reactions with 4% dimethyl sulfoxide for 30 cycles; 45 s at 94°C, 1 min at 60°C, and 90 s at 72°C in a 50-µl reaction mixture) using two sets of primers: HYGfor2 and HYGrev2 and HYGfor3 and HYGrev3 (Table 2).
The same approach was used to replace the pcaL gene on RFMD1 with an Aprar gene flanked by FRT sequences (Fig. 1c). The template for amplification was a gel-purified 1,384-bp EcoRI/HindIII fragment of pIJ773. The 3' ends of the PCR primers (PCALfor1 and PCALrev1) matched 19-nucleotide and 20-nucleotide extensions, respectively, of the pUJ773 sequence flanking the FRT sequences. The 5' end of PCALfor1 contained 39 nucleotides corresponding to the sense strand upstream of pcaL ending in the start codon. The 5' end of PCALrev1 contained 39 nucleotides corresponding to the antisense strand ending in the stop codon. Transformants of E. coli BW25113 containing pKD46 and RFMD1 were selected on ampicillin and hygromycin. Cells transformed with the PCR-extended Aprar cassette were selected on LB containing hygromycin plus apramycin at 37°C. Successful transformants containing a copy of the Hmr- and Aprar-mutagenized fosmid, RFMD2, were checked by PCR using three sets of oligonucleotides: PCALfor2 and PCALrev2; APRAfor and APRArev; and PCALfor2 and APRArev (Table 2).
Allelic exchange in RHA1.
Fosmid RFMD2 was isolated, purified, and introduced into DH10b/pUZ8002 by electroporation. The fosmid was then transferred to Rhodococcus sp. strain RHA1 by intergeneric conjugation as described in (47). Aprar Hms exconjugants were selected as double-crossover mutants. The mutants were analyzed by PCR with three sets of primers: PCALfor2 and PCALrev2; PCALfor2 and APRArev; and PCALfor3 and PCALrev3 (Table 2). The frequency of double-crossover mutants among the exconjugants was
12%.
Bioinformatic analyses. Gene prediction and annotation were performed using integrated automated and manual approaches developed at Oak Ridge National Laboratories. Briefly, the automated step made use of three gene finders: Critica v.1.05 (3), Glimmer (20), and Generation. The annotation and locations of predicted open reading frames (ORFs) of interest were then evaluated using a variety of tools. ORF function and position were confirmed with BLASTP sequence alignments (1) to NCBI-nr and PFAM, TIGRfam, COGS, KEGG protein databases. Interproscan (2) was used with the ProfileScan, BLASTProDom, HMMPfam, HMMSMART, and ScanRegExp databases to search for conserved domains and motifs and to validate predicted gene function. PSI-BLAST alignments were used to identify ORF function that was not predicted by the above-mentioned searches. When BLAST alignments were performed, the global percent identity (over the full sequence length) was recorded. Finally, for the genes whose protein products were identified, the sequence information of signature peptides was used to verify gene coordinates.
|
|
|---|
A search of the current RHA1 genome assembly (www.rhodococcus.ca) revealed four clusters of genes which together could encode the pathways responsible for the catabolism of benzoate and phthalate in RHA1 (Fig. 2). The chromosomal ben genes encode a dioxygenase (benABC), dihydrodiol dehydrogenase (benD), and transporter (benK) as previously reported (50) (Fig. 3a). However, the sequence of benA reported here and supported by proteomic data (see below) differed to the published sequence (gi 16506124). It is unclear whether these differences represent sequencing errors or spontaneous mutations.
![]() View larger version (15K): [in a new window] |
FIG. 2. Proposed pathways of benzoate and phthalate degradation in Rhodococcus sp. strain RHA1. Catabolites are numbered as follows: 1, phthalate; 2, phthalate 3,4-dihydrodiol; 3, 3,4-dihydroxyphthalate; 4, protocatechuate; 5, -carboxymuconate; 6, -carboxymuconolactone; 7, benzoate; 8, cis-1,6 dihydroxy-2,4-cyclohexadiene-1-carboxylic acid; 9, catechol; 10, cis, cis-muconate; 11, muconolactone; 12, ß-ketoadipate enol-lactone; 13, ß-ketoadipate. The molecular mass-corrected normalized volumes of each protein observed on benzoate and phthalate are indicated in parentheses. In this figure, a value of 0 indicates that the protein was not detected on the gel.
|
![]() View larger version (19K): [in a new window] |
FIG. 3. Benzoate and phthalate catabolic genes in Rhodococcus sp. strain RHA1. (a) the chromosomal cat-ben cluster, (b) the pad-pat-tpa cluster found on pRHL1 and pRHL2, (c) the chromosomal pca cluster, and (d) a cluster of unknown function that contains a catD homolog. Arrows represent different types of genes as follows: diagonally striped, pad cluster; vertically striped, pat cluster; horizontally striped, tpa cluster; diamonds, cat cluster; hatched diamonds, ben cluster; squares, pca cluster; grey, transcriptional regulator; black, DNA recombination; white, miscellaneous; and waved, hypothetical. Numbers above the arrows indicate the ORF numbers (Table 1) from the RHA1 genome assembly (www.rhodococcus.ca). Solid lines below the pad-pat-tpa cluster indicate the two sections that are present in pRHL2 but not in pRHL1. The products of the depicted genes identified in the study are listed in Table 4.
|
|
View this table: [in a new window] |
TABLE 3. Gene annotationa
|
30% sequence identity with BenR from Pseudomonas putida (17) and may regulate the transcription of the ben and cyp genes, respectively. However, we could not assign specific functions based on the current data. Inspection of the 183 bp between fpr254A1 and benA revealed a 14-bp repeat centered at 94 bp with respect to the mapped transcription start (50). Two benK homologs and two additional benR homologs were found in the RHA1 genome (http://www.rhodococcus.ca/publications/supplementary/JBact05B.pdf). Identical copies of a pad cluster were found on plasmids pRHL1 and pRHL2 (Fig. 3b). The cluster contains seven genes predicted to encode a regulatory protein (padR), and the enzymes that transform phthalate to protocatechuate: a 3,4-dioxygenase (padAaAbAcAd; in gene names, lowercase letters are used to designate subunits and numbers are used to designate isozymes), a dehydrogenase (padB), and a decarboxylase (padC) (Fig. 2, Table 3). These ORFs share 99% sequence identity with putative ORFs carried by pDK3, the 750-kb plasmid of Rhodococcus sp. DK17 (accession number AY502076).
Further analysis of the nucleotide sequence surrounding the pad genes revealed that they are part of a 32.1-kb duplication which contains 34 predicted ORFs organized in four major clusters flanked by transposase-encoding genes (Fig. 3b; Table 3). The pad cluster is at the 3' end of this duplication. Two apparently related upstream clusters were named pat and tpa based on their predicted catabolic functions. The pat genes were predicted to encode four subunits of an ABC-type phthalate transporter system (patDACB) and a phthalate ester hydrolase (patE). The predicted gene products are 62 to 71% identical to the corresponding ptr-encoded proteins from Arthrobacter keyseri 12B (25). In this strain, the ptr cluster is flanked by operons involved in phthalate and protocatechuate catabolism.
The tpa genes were predicted to be involved in terephthalate catabolism. These include tpaAaAb, whose products share 68% sequence identity to the
and ß subunits of terephthalate 1,2-dioxygenase from Delftia tsuruhatensis T7 (71), tpaB, which likely encodes the cognate reductase, tpaC, predicted to encode a dehydrogenase, and tpaK, predicted to encode an aromatic acid permease. The cluster also includes a divergently transcribed gene whose product shares 39% identity with PcaR, an IclR-type transcriptional regulator. The fourth cluster in the 32.1-kb duplication lies several kb upstream of the tpa cluster. Although several of the predicted gene products show high similarity to catabolic enzymes such as cyclohexanone monooxygenase, the precise biological function of these genes is unknown.
The 32.1-kb duplication also includes genes predicted to encode an integrase (rha10163, 4 transposases (rha10188 rha10164 rha10165 and rha10167, and a reverse transcriptase (rha10166. Three of these genes are missing from the corresponding gene cluster on pRHL1 (solid bars in Fig. 3b). Moreover, the entire 32.1-kb region found on pRHL1 shares 99% sequence identity with a similar region of pDK3 of Rhodococcus sp. DK17. Indeed pDK3 and pRHL1 share two regions of 99% sequence identity: a 73-kb region that contains the tpa, pat, and pad clusters and a 22-kb region (data not shown). The similarity of the copy of the tpa-pat-pad island present on pRHL1 to that of pDK3 as well as the presence of additional transposase genes in the copy found on pRHL2 suggest that the copy on pRHL1 is the original in RHA1.
The RHA1 genome contains several genes potentially involved in catechol and protocatechuate catabolism. The chromosome contains genes predicted to encode two branches of the ß-ketoadipate pathway: the cat and pca genes specify the catabolism of catechol and protocatechuate, respectively. The putative catRABC operon is located within 7 kb of the ben genes (Fig. 3a) and shares over 96% sequence identity and the same organization as the corresponding genes in Rhodococcus opacus 1CP (28) (Table 3). The three encoded enzymes, catechol 1,2-dioxygenase (catA), cis,cis-muconate lactonizing enzyme (catB), and muconolactone isomerase (catC), are predicted to transform catechol to ß-ketoadipate enol-lactone but not to trichloroacetic acid (TCA) cycle intermediates. In contrast, the pca cluster (Fig. 3c) includes genes predicted to encode the enzymes (PcaJIHGBLF) required to convert protocatechuate to the TCA cycle intermediates (Fig. 2). The predicted products of these genes share high 97 to 99% sequence identity with characterized homologs from related actinomycetes (Table 3).
The organization of the pca genes in two putative divergently transcribed operons (pcaJI and pcaHGBLRF) is similar to their organization in R. opacus 1CP (29). Of particular note, pcaL in RHA1, R. opacus sp. 1CP, and Streptomyces sp. 2065 (41) appears to encode a bifunctional ß-ketoadipate enol-lactone-hydrolyzing enzyme. The gene appears to have arisen from a fusion of pcaD and pcaC, which encode ß-ketoadipate enol-lactone hydrolase and
-carboxymuconolactone decarboxylase, respectively, in Acinetobacter sp. ADP1 and pseudomonads (reviewed in reference 43).
Paralogs of the cat and pca genes found on the RHA1 genome include two homologs of catA, two of pcaIJ, two of pcaC encoding decarboxylases, and five of pcaF (http://www.rhodococcus.ca/publications/supplementary/JBact05B.pdf). The homologs appear to be functional as key catalytic residues are conserved, although their physiological roles in RHA1 are for the most part unclear. One of the pcaF homologs is involved in phenylacetate catabolism and was annotated as paaE (60). The annotation of the genes mentioned in this section is also available at www.rhodococcus.ca.
Identification of phthalate- and benzoate-catabolizing enzymes. To identify the pathways responsible for the catabolism of benzoate and phthalate, respectively, in Rhodococcus sp. strain RHA1, the cytosolic proteomes of cells grown on phthalate, benzoate or pyruvate as sole carbon source were quantitatively compared by two-dimensional gel electrophoresis. Approximately 1,500 protein spots were resolved per gel (see Fig. 4 for representative gel sections). Some of the proteins appeared on the gel as a horizontal series of spots suggesting that they represent multiple species differing in charge (e.g., PadAa in Fig. 3c). This may be caused by either posttranslational modification within the cell or chemical modification during sample preparation. Mass spectrometric analyses supports the latter for the proteins reported here as changes in the masses of peptides originating from "multiple" protein spots were consistent with carbamylation. Optimization of sample preparation, particularly with respect to temperature and buffer composition, minimized but did not completely eliminate carbamylation. For carbamylated proteins, the isoelectric point and molecular mass of the major spot in a series were recorded, and the expression difference was calculated based on the summed signal intensities of all the spots in the series.
![]() View larger version (72K): [in a new window] |
FIG. 4. Sections of two-dimensional gels showing parts of the Rhodococcus sp. strain RHA1 cellular proteome. Corresponding gel sections are shown from cells grown on pyruvate (a), benzoate (b), or phthalate (c). The protein spots of interest are shown with arrows. Panel 1 shows the proteins (BenA and PadAaAb) that were only observed in benzoate- or phthalate-grown cells. Panel 2 shows proteins (PcaLFB) that appeared in all the conditions tested but were more abundant in phthalate and benzoate samples.
|
|
View this table: [in a new window] |
TABLE 4. Identification of the proteins involved in benzoate and phthalate catabolism.
|
Three of the four ben gene products (BenABD) and two of the four cat gene products (CatAB) were identified in benzoate-grown cells (Fig. 4b, Table 4). Each of these five proteins was only observed in benzoate-grown cells (Fig. 2). Failure to identify BenC, CatC, and CatR does not imply their absence in benzoate-grown cells, as
50% of the protein spots analyzed by mass spectrometry did not yield usable data. Moreover, CatC is unlikely to be detected using the current two-dimensional gel electrophoresis protocols due to its low molecular mass (10 kDa). The components of the type V cytochrome P450 were not identified on the gels. Moreover, spectrophotometric analyses failed to detect the presence of a cytochrome P450 in benzoate-grown cells (data not shown). Finally, transcriptomic data analysis confirmed that cyp254A1 and frp254A1 are not up-regulated in benzoate-grown RHA1 cells (H. Hara and W. W. Mohn, in preparation).
Seven of the nine predicted pca-encoded enzymes (PcaIJHGBLF) were identified. These were more abundant in phthalate- and benzoate-grown cells than in the pyruvate-grown controls (Fig. 4, Table 4). Protocatechuate 3,4-dioxygenase (PcaHG) was only detected in the phthalate and benzoate proteomes. The other pca products, including the regulatory protein PcaR, were present in low but detectable levels in pyruvate-grown cells (Fig. 2, Table 4).
Other than the enzymes that are involved in benzoate and phthalate catabolic pathways, we identified 31 other proteins that were more abundant in phthalate- and/or benzoate-grown cell samples (http://www.rhodococcus.ca/publications/supplementary/JBact05C.pdf). Importantly, these included none of the other ben, pca, and cat homologues identified in the RHA1 genome. Identified proteins of known physiological functions included TCA cycle enzymes as well as amino acid, fatty acid, and nucleotide biosynthetic enzymes. None of the products of the tpa genes, thought to be involved in terephthalate catabolism, were detected. Transcriptomic data indicate that these genes, but not the pad genes, are strongly up-regulated in terephthalate-grown cells (H. Hara and W. W. Mohn, in preparation). The more abundant TCA enzymes included citrate synthase, isocitrate dehydrogenase, succinate dehydrogenase, and malate dehydrogenase as well as the
-subunit of succinyl-CoA ligase, an enzyme that couples the hydrolysis of succinyl-CoA to the synthesis of ATP.
Two enzymes of the central metabolic pathway were more abundant in pyruvate-grown cells: isocitrate lyase, an enzyme of the glyoxylate shunt; and acetate-CoA ligase, one of the two enzymes of the pyruvate oxidation pathway. The more abundant TCA cycle enzymes during growth on benzoate and phthalate is consistent with end products of the ß-ketoadipate pathway, succinyl-CoA and acetyl-CoA, feeding directly into the cycle. In contrast, increased abundance of the glyoxylate shunt enzymes during growth on pyruvate is consistent with the general observation that these enzymes are required for net assimilation of carbon when the carbon source enters the TCA cycle solely at the level of acetyl-CoA (52).
Relative protein abundance. To better estimate the relative levels of protein abundance, the normalized volumes (NV) of the protein spots were corrected for their molecular mass in kDa (MWc): (MWc-NV = [NV/MW] x 1,000). For three of the four two-subunit enzymes of the studied pathways, the MWc-NV of the two subunits are in good agreement (Fig. 2). The lack of agreement for the fourth enzyme, PcaIJ, likely reflects the different intensities with which Sypro Ruby interacts with these two denatured polypeptides. Nevertheless, the data revealed several interesting features. First, the Pca proteins were consistently more abundant in benzoate-grown cells compared to phthalate-grown cells (Fig. 2, Table 4). This may reflect the higher growth rate of RHA1 on benzoate. Second, a subset of Pca proteins, PcaBLRF, were detected in pyruvate-grown cells, but PcaHG were not. This indicates that pcaBLRF expression may be regulated independently from pcaHG. Finally, the data indicate that the most abundant proteins in each pathway branch are the oxygenases BenAB, PadAaAb, CatA, and PcaHG. This may reflect the relatively low cytoplasmic concentration of O2 in exponentially growing cells.
Enzyme activity.
The RHA1 genome sequence and proteomic data suggested that the final steps of phthalate and benzoate catabolism in RHA1 are accomplished by the same set of pca-encoded enzymes (Fig. 2). To verify this conclusion, the activities of three enzymes representing different branches of the ß-ketoadipate pathway were determined in whole-cell extracts of phthalate-, benzoate-, and pyruvate-grown cells. The specific activity of protocatechuate 3,4-dioxygenase (pcaHG), which cleaves protocatechuate to
-carboxymuconate, was 10- and 6-fold higher in RHA1 cells grown in benzoate and phthalate, respectively, than in pyruvate (Table 5). Catechol 1,2-dioxygenase (catA) activity was 133-fold greater in benzoate- and only 1.3-fold higher in phthalate-grown cells than in pyruvate-grown cells. No protocatechuate or catechol meta-cleavage activities were detected in extracts of cells grown on phthalate or benzoate. However, catechol 2,3-dioxygenase activity was detected in pyruvate-grown cells (data not shown). ß-Ketoadipate:succinyl-CoA transferase (pcaIJ) catalyzes the penultimate step in the ß-ketoadipate pathway. This activity was detected in all the samples and was twofold higher in cells grown on an aromatic carboxylic acid.
|
View this table: [in a new window] |
TABLE 5. Enzyme activities in cell lysates of RHA1 cells grown on different substratesa
|
Analysis of pcaL gene knockout.
The first shared step in the catabolism of phthalate and benzoate in RHA1 is catalyzed by the bifunctional, pcaL-encoded ß-ketoadipate enol-lactone hydrolase. The essential role of this enzyme in benzoate and phthalate catabolism was tested by investigating the phenotype of RHA1_005, a mutant of RHA1 in which pcaL was replaced with an Aprar cassette using a
RED-based methodology developed for Streptomyces spp. (35) (Fig. 1). RHA1_005 was resistant to apramycin but sensitive to hygromycin, consistent with allelic exchange and loss of the donor fosmid. PCR analysis of RHA1_005 genomic DNA using two sets of primers confirmed the presence of the Aprar cassette and loss of pcaL (Fig. 5). More particularly, primers flanking the pcaL gene yielded a larger amplicon from RHA1_005 genomic DNA (lane 2) than from RHA1 genomic DNA (lane 3). PCR performed with primers internal to pcaL yielded the expected amplicon from RHA1 genomic DNA (lane 7) and no product from RHA1_005 DNA (lane 6).
![]() View larger version (43K): [in a new window] |
FIG. 5. PCR analysis of the pcaL replacement mutant RHA1_005. The results and expected fragment sizes are shown in panels A and B, respectively. Reactions were performed with primers that are external (ext., lanes 2 to 5, PCALfor2/PCALrev2) or internal to pcaL (int, lanes 6 to 9, PCALfor3/PCALrev3). Reactions contained RHA1_005 genomic DNA (lanes 2 and 6), RHA1 genomic DNA (lanes 3 and 7), fosmid RF00111bAO4 (lanes 4 and 8), or no template DNA (lanes 5 and 9). Lanes 1 and 10 were loaded with molecular mass markers (1-kb DNA ladder, Invitrogen).
|
We performed proteomic analyses to investigate the catabolism of benzoate in RHA1_005. Two-dimensional gel electrophoresis analysis confirmed the presence of all ben-, cat-, and pca-encoded proteins detected in benzoate-grown wild-type cells except for PcaL. In addition, 24 proteins were at least fivefold more abundant in benzoate-grown RHA1_005, of which 20 were not detected in benzoate-grown RHA1. Half of these 24 proteins were identified by peptide fingerprint analysis. One of the identified proteins, encoded by rha05436 shares greatest sequence identity (18%) with CatD from Acinetobacter sp. ADP1 (70) and 19% sequence identity with the N-terminal hydrolase domain of PcaL (Tables 3 and 4). It is therefore likely that this enzyme functions as a ß-ketoadipate enol-lactone hydrolase, compensating for the loss of PcaL in RHA1_005.
Rha05436 resides in a putative operon containing five other genes (Fig. 3D), three of whose products were also more abundant in RHA1_005 (Table 4). These show sequence similarity to an acetyl-CoA acetyltransferase (CatF homolog), an O-succinylbenzoate-CoA ligase, and an acyl-CoA dehydrogenase, respectively (Table 3). None of these proteins were detected in wild-type RHA1 cells under any of the growth conditions studied. Four of the identified proteins appear to be involved in iron acquisition and storage, including two that appear to be involved in the biosynthesis of a catecholic siderophore, possibly an enterobactin-type molecule. Other proteins were identified in benzoate-grown RHA1_005 (http://www.rhodococcus.ca/publications/supplementary/JBact05B.pdf and http://www.rhodococcus.ca/publications/supplementary/JBact05C.pdf), and their current annotation is available at www.rhodococcus.ca.
|
|
|---|
-carboxy-muconolactone decarboxylase and an enol-lactone hydrolase in separate domains. In other bacteria, including pseudomonads, these enzymes are encoded by pcaC and pcaD/catD, respectively (43). In RHA1, only the protocatechuate branch utilizes the decarboxylase activity of PcaL. The pcaL gene is not essential for growth of RHA1 on benzoate, apparently because the CatD homolog is expressed in the pcaL mutant grown on benzoate, compensating for loss of PcaL. The current data also identify an ABC-type transporter that is likely involved in the uptake of phthalates. Genes encoding a related transporter were identified in A. keyseri 12B (25), but no functional data were obtained. The ATPase of the transporter (PatA) and a probable phthalate ester hydrolase (PatE) were highly abundant in the cytoplasm of phthalate-grown cells. Consistent with its cytoplasmic localization, PatE lacked a predicted secretion signal sequence. The localization of this enzyme in the cytoplasm suggests that phthalate esters are hydrolyzed in the cytoplasm after their uptake, either by the pat-encoded ABC transporter or by another route. In this respect, symporter-type phthalate permeases are inhibited by substituted phthalates but not structurally related compounds, such as 2-Cl benzoate, that lack one of the carboxylates, suggesting that vicinal carboxylate is an important substrate-binding determinant (13). We are currently investigating the role of PatDACB in the transport of phthalates, terephthalates, and related compounds.
Our comparison of the ß-ketoadipate pathway of RHA1 to that of other bacteria in which it occurs revealed several features of the pathway that appear to be unique to rhodococci and other actinomycetes with respect to component enzymes, gene organization, and regulation. Thus, analyses of genome sequences indicate that the same branched ß-ketoadipate pathway also occurs in Streptomyces avermitilis MA-4680, Streptomyces coelicolor A3(2), and Corynebacterium glutamicum (45). Moreover a bifunctional PcaL has been reported in R. opacus sp. strain 1CP (29) and Streptomyces sp. strain 2065 (41). In contrast, the pathway branches do not converge at all in Acinetobacter sp. strain ADP1 (23) and converge at a different point, ß-ketoadipate, in Ralstonia eutropha (44).
The pathway in pseudomonads appears to be similar to that in actinomycetes in that the branches converge at the enol-lactone (42, 61). However, the pseudomonads do not use a bifunctional enzyme at the point of convergence: the
-carboxymuconolactone decarboxylase and the enol-lactone hydrolase are encoded by pcaC and pcaD, respectively (43). Nonetheless, the bifunctional PcaL is not unique to actinomycetes: it also occurs in Ralstonia metallidurans and Caulobacter crescentus (43). Finally, the ß-ketoadipate pathways of pseudomonads and actinomycetes appear to differ by the presence of a ß-ketoadipate transporter (PcaT) in the former (43).
An analysis of 22 complete or partial bacterial genomes containing the pca and cat genes revealed that the latter appear to be organized in a fashion that is unique to rhodococci and most similar to that of the closely related corynebacteria. The pca genes are organized in a single cluster in all actinomycetes in which they have been found as well as in Caulobacter crescentus (62), an
-proteobacterium, and Acinetobacter sp. ADP1 (11), a
-proteobacterium. Actinomycetes containing the pca genes include R. opacus 1CP (29), C. glutamicum sp. ATCC 13032 (BX927155) (45), Streptomyces sp. 2065 (41), S. coelicolor A3(2) (SC0939128) (7), and S. avermitilis MA-4680 (AP005027.1) (63). In contrast, the pca genes can be arranged in up to three clusters in pseudomonads (42). Multiple pca clusters also occur in ß-proteobacteria such as Burkholderia pseudomallei and R. metallidurans (43).
Nevertheless, the organization of the pca genes in a single cluster of two divergently transcribed operons with the gene order of RHA1 appears to be unique to rhodococci. In C. glutamicum the gene order is different, and in streptomycetes, the genes appear to be arranged in a single operon. In all bacteria, the cat genes are usually organized in a single cluster (reviewed in reference 43). However, permutations occur with respect to gene order (e.g., catRBAC in the streptomycetes sequenced to date) and the presence of additional genes in the transcriptional unit (e.g., in Acinetobacter, there are six genes in the cat operon). The order of the genes in RHA1, catRABC, is seen in R. opacus 1CP (29), R. erythropolis AN-13 (D83237), and R. erythropolis CCM2595 (AJ605581). In C. glutamicum ATCC 13032, catR is not adjacent to the other genes.
There are insufficient data to compare the regulation of the ß-ketoadipate pathway genes in rhodococci to that in other organisms. It is nonetheless interesting to note the proteomic data which suggest that pcaBLRF are independently regulated to pcaHG despite their organization in an apparent operon (Fig. 3c). Consistent with this possibility, the gap between pcaG and pcaB, 28 bp, is larger than that between the other genes, which typically overlap 1 or 4 bp. This same spacing is seen in the pcaHGBLRF cluster of R. opacus 1CP (29). Further related to regulation, we note that PcaR, BenR, and PadR of RHA1 all belong to the same family of regulators as the identically named regulators in other bacteria. In contrast, the CatRs of RHA1 and R. opacus 1CP (28) belong to the IclR family, not to the LysR family to which the CatRs of other strains often belong (79). Finally, not all ß-ketoadipate pathway enzymes that are expressed in RHA1 during growth on benzoate are utilized (e.g., PcaHG and PcaB), as observed in pseudomonads (37). The physiological relevance of this apparent inefficiency is unclear, particularly given the apparent differential regulation of pcaHG and pcaBLRF: it is possible that these bacteria only encounter mixtures of compounds degraded via benzoate and protocatechuate in their natural environments.
The emerging data, including those from the current study, suggest that the catabolism of aromatic compounds in rhodococci is organized in a fashion similar to that found in the better-studied pseudomonads (42): a large number of "peripheral" pathways funnel a range of natural and xenobiotic compounds into a restricted number of "central" pathways. The latter, exemplified by the ß-ketoadipate pathway, complete the transformation of these compounds to TCA cycle intermediates. Analyses of the genomic sequences of four pseudomonad strains, P. putida KT2440, P. fluorescens Pf0-1, P. aeruginosa PAO1, and P. syringae pv. tomato DC3000 (42, 61), together with functional studies have identified at least 38 peripheral pathways, some of which are strain specific, and five conserved central pathways. The current genomic and proteomic data suggest that in RHA1, three peripheral pathways funnel phthalate, terephthalate, and 4-hydroxybenzoate to the ß-ketoadipate pathway via protocatechuate, while three others funnel benzoate, phenol, and 2-ethoxyphenol to the ß-ketoadipate pathway via catechol (M. A. Patrauchan, H. Hara, and L. D. Eltis, unpublished data).
The duplication of the peripheral pathways responsible for phthalate and terephthalate catabolism (tpa-pat-pad) establishes that catabolic redundancy in RHA1 is not confined to the bph, etb, and ebd genes, which specify pathways involved in biphenyl and ethylbenzene catabolism (51, 55, 67). While redundancy in catabolic genes has been cited as a trait of Rhodococcus (53, 81), the cited examples appear to involve paralogs involved in distinct, nonredundant physiological processes (40, 77). The duplication of catabolic pathways in RHA1 may increase the organism's potential to adapt to new carbon sources. This hypothesis is supported by the growth of the pcaL deletion mutant on benzoate: in this case, adaptation apparently involved recruitment of a catD ortholog. Nevertheless, it is unclear whether such duplication is shared by divergent rhodococcal strains and whether this redundancy is associated primarily with catabolic genes.
The physiological role of the gene cluster containing the catD ortholog, rha05436 remains unclear. Considering the likely enzymatic activities of the encoded proteins (Table 3), as well as the simultaneous expression of a probable catecholic siderophore biosynthetic pathway (Tables 3 and 4), it is tempting to speculate that the cluster containing the catD homolog is involved in the degradation of the catecholic siderophore and that the two gene clusters are under the control of a common regulatory mechanism. Interestingly, the RHA1 genome contains 21 genes whose products share as much or greater sequence identity with the enol-lactone hydrolase domain of PcaL than the rha0543-encoded ortholog does (results not shown). Nevertheless, it is unclear whether all of these enzymes could function in the ß-ketoadipate pathway. Moreover, although the pcaL mutant did not adapt to phthalate, the RHA1 genome also contains two pcaC homologs whose products share 32 and 53% sequence identity with the decarboxylase domain of PcaL. Finally, the growth of the pcaL mutant on benzoate illustrates the limitation of assigning gene function based on gene deletion studies.
In conclusion, this study establishes the particular configuration of the ß-ketoadipate pathway in RHA1 and indicates that the organization of the genes encoding this pathway is characteristic of rhodococci and related actinomycetes. Moreover, our analyses indicate that the overall organization of the catabolism of aromatic compounds in rhodococci is similar to that described in pseudomonads, with multiple peripheral pathways feeding into a limited number of central pathways. The redundancy of the peripheral phthalate and terephthalate pathways in RHA1 stands in marked contrast to the convergent nature of the ß-ketoadipate pathway. It is possible that this configuration augments the ability of RHA1 to adapt to newly encountered carbon sources, much as the pcaL mutant adapted to benzoate by recruiting an ortholog of the deleted gene. Finally, several lines of evidence provided by the current study, including the arrangement of the catabolic islands and the expression of pcaHG during growth on benzoate, indicate that RHA1 might normally grow on mixtures of aromatic compounds in the environment. We are currently investigating the regulation of these different pathways on single and multiple carbon sources.
Masao Fukuda is thanked for helpful discussions. Robert Olafson, Derek Smith, and other members of the Proteomics Centre, University of Victoria, are thanked for their assistance with the mass spectrometry. Ritesh Patel, Matthew J. Myhre, Clinton Fernandez, and Michael McLeod are thanked for their assistance in bioinformatic analyses. Hirofumi Hara is thanked for valuable assistance with the enzyme assays and for sharing microarray data.
|
|
|---|
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»