Previous Article | Next Article ![]()
Applied and Environmental Microbiology, November 2007, p. 6930-6938, Vol. 73, No. 21
0099-2240/07/$08.00+0 doi:10.1128/AEM.01697-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Department of Civil and Environmental Engineering, University of California, Berkeley, California 94720-1710,1 Department of Microbiology and Immunology, Life Sciences Institute, University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z3,2 Departments of Chemical Engineering, Biology, and Civil Engineering, Texas A&M University, College Station, Texas 77843-3122,3 Earth Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 947204
Received 24 July 2007/ Accepted 6 September 2007
|
|
|---|
|
|
|---|
Microorganisms grown on substrates such as propane, methane, and toluene have been shown to rapidly oxidize NDMA in the laboratory (7, 25). In these cases, evidence from inhibition and induction experiments along with observations of requisite oxygen consumption suggests that propane monooxygenases (PrMO), soluble methane monooxygenases (sMMO), and toluene monooxygenases (TMO) are most likely involved in these transformations. In addition, experiments with Escherichia coli clones expressing TMO inserts confirmed the role of toluene 4-monooxygenase (T4MO) in NDMA oxidation, while cupric selection for soluble rather than particulate MMO confirmed the role of sMMO (25). The involvement of PrMO is less understood, as the traditional boundary can blur between enzymes oxidizing gaseous and liquid n-alkanes. Liquid alkanes are typically oxidized by alkane monooxygenases (AlkMO), but AlkMO can be induced by propane in some bacteria but not in others (10, 16). Regardless of the class of monooxygenase involved, NDMA is transformed with little observed benefit to the cells and no evidence of cellular growth, despite the production of oxidized products, including formaldehyde, methylamine, and methanol, that can be incorporated into primary metabolic pathways (7, 24, 33). Limited evidence for metabolism suggests that non-energy-generating transformations, such as cometabolic oxidation reactions, play an important role in the biological attenuation of NDMA.
Rhodococci are soil heterotrophs of the order Actinomycetales with a noted diversity of functional enzymatic activities (9, 14). Collectively, this order is biologically and economically significant for the production of a diverse array of enzymes involved in the production of commercial secondary metabolites, antibiotics, and the metabolism of xenobiotic compounds for environmental and industrial applications (21). Global genomic, transcriptomic, and functional analyses of Rhodococcus sp. strain RHA1 reveal tremendous enzymatic diversity with the potential to grow on a wide variety of aromatic compounds, carbohydrates, nitriles, and steroids as the sole carbon and energy sources (see reference 17 and references therein). Indeed, the genome of Rhodococcus sp. strain RHA1 is predicted to encode over 200 oxygenases, including both PrMO and AlkMO gene clusters. The genetic blueprint provided by the annotated genome and the development of a corresponding global microarray facilitate the identification of genes responsible for physiological traits of RHA1, especially for growth and enzymatic activity. For this reason, strain RHA1 was selected as a model organism to better understand the genetics and biochemistry of NDMA transformation.
To gain insight into NDMA degradation by Rhodococcus sp. strain RHA1, we analyzed the effects of propane on gene expression and NDMA removal. Here, we report the first experimental evidence for NDMA degradation by a PrMO. First, the kinetics associated with NDMA degradation in strain RHA1 are characterized. Then, the candidate genes associated with propane and NDMA oxidation are identified and quantified through differential expression, as assayed by global transcriptional microarray analysis and reverse transcriptase-quantitative PCR (RT-qPCR). Finally, targeted gene disruption is used to confirm the role of PrMO in NDMA degradation and exclude the role of AlkMO in the observed degradation.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Bacterial strains and plasmids used in this study
|
Culture growth phase was determined by monitoring optical density at 600 nm (OD600). Cells were harvested from culture medium in the late exponential phase of growth at OD600 of
0.7 and 1.5 for cells grown, respectively, on propane or liquid organics. Cells were harvested with a Beckman Avanti J-301 centrifuge (Palo Alto, CA) at 15,000 x g for 5 min. Cells exposed to propane were centrifuged and transferred to a new container to remove residual propane. The resulting pellet was then suspended in 0.1 M phosphate buffer (pH 7) solution. For cells grown with liquid substrates, this washing process was repeated two more times. Washed cells were suspended with buffer to a target density to optimize measurement of transformation rates (adjusted OD600 of between 0.1 and 7.0), and a fraction of the cells were frozen for future protein analysis.
Quantification of NMDA removal.
Cell suspensions to evaluate and quantify the biodegradation of NDMA were incubated in 125-ml bottles sealed with Teflon-lined Miniert valves (Altech, Deerfield, IL). These flasks contained washed cells suspended in phosphate buffer as described above. NMDA (99+%) was purchased from Acros Organics (Geel, Belgium), and additions to experimental cultures, controls, and standards have been described previously (25). NDMA extractions were performed by removing 2-ml samples at each time point followed by equilibration with an equal volume of high-purity methylene chloride (EM Science, Darmstadt, Germany). Methylene chloride extracts containing NDMA were analyzed by previously described methods (18, 25) involving tandem mass spectrometry. The detection limit for the liquid-liquid extraction as determined by standard curve was approximately 5 µg NDMA liter–1.
NDMA degradation rates as a function of protein density were obtained by monitoring initial loss during cellular incubations in phosphate buffer (24). Each rate consisted of an average of four linearly spaced time points run over 2 h. Biomass was quantified as mass of cellular protein, and there was no significant change in cell density during the course of these incubations. Cellular pellets from frozen 1.5-ml samples were suspended in 210 µl of 48 mM NaOH. The mixture was sheared by bead beating for 5 min and boiled at 100°C for 20 min. The digest was then centrifuged in an Eppendorf 5417C (Hamburg, Germany) for 10 min at 10,000 x g to remove cellular debris from the supernatant. Protein mass per volume was quantified from 50 µl or appropriate dilutions of the supernatant using the Coomassie Plus protein assay reagent kit with bovine serum albumin as the standard (Pierce Biotechnology, Rockford, IL). Degradation rates were graphed as a function of substrate concentration to determine the interdependence of these variables. An iterative best fit for nonlinear regression with 95% confidence intervals was applied to the Monod kinetic model for a constant cell density (equation 1):
![]() | (1) |
Analysis of global gene expression using spotted microarrays.
For transcriptomic analysis, 65-ml or 80-ml liquid cultures of strain RHA1 were grown in sealed 1-liter flasks containing minimal medium amended, respectively, with 23 mM pyruvate or atmospheric air containing 20% gaseous bulk propane (+99.5% purity) as the sole electron donor and carbon source. Triplicate cultures for each condition were harvested in late-exponential growth. OD600 of approximately 0.7 and 1.3, respectively, were selected to correspond to 70% of the maximal OD600 reached in these incubations as determined by prior growth curves. Upon achieving the target density, 1/10 volumes of "stop solution" (5% phenol [pH 5] in ethanol) were added (1). Cells were collected by centrifugation at 4,900 x g for 10 min at 4°C, suspended in 1.0 ml of the supernatant plus 2.0 ml RNAprotect (QIAGEN), and incubated for 5 min at room temperature. Cells were then centrifuged at 13,000 x g for 2 min at room temperature. Pellets representing 40 ml of culture were each frozen on dry ice and stored at –80°C. RNA extraction was performed on the harvested pellets by adapting previous methods (8). Total RNA isolation involved vortexing with glass beads, hot phenol, and sodium dodecyl sulfate at final concentrations of 14.3% and 0.9% (vol/vol), respectively. Debris was precipitated with acetate followed by the addition of 4.0 ml phenol-chloroform (1:1 [vol/vol]). Nucleic acids were precipitated with acetate plus isopropanol, treated with DNase, and purified with an RNeasy mini column (QIAGEN).
Synthesis of cDNA from the extracted RNA, indirect Cy labeling, and microarray hybridizations were performed as described previously (8), with the following modifications. The cDNA synthesis mixture included 1.5 µg random hexamer primers (Invitrogen) per 6.0 µg RNA, which was brought to 15.3 µl with diethyl pyrocarbonate-treated water. After RNA denaturation for 10 min at 70°C followed by cooling for 5 min on ice, cDNA synthesis components were added to final concentrations of 0.46 mM each dATP, dCTP, and dGTP; 0.19 mM dTTP; 0.28 mM aminoallyl-dUTP (Ambion); 0.01 M dithiothreitol; 10 U RNaseOUT (Invitrogen); and other ingredients as described previously (8). Equal amounts of differentially labeled cDNA, consisting of 50 million pixels measured by ImageQuant 5.2 (Molecular Dynamics), from propane- and pyruvate-grown cells were hybridized at 42°C for 17 h. After the automated washes, the slides were dipped in 0.2x SSC (1xSSC is 0.15 M NaCl plus 0.015 M sodium citrate) and dried by centrifugation at 225 x g for 5 min at room temperature. For one of the three hybridizations, the Cy3 and Cy5 dyes were reversed (i.e., cDNA from the propane treatment was labeled with Cy5 rather than Cy3) to control for dye bias (29). The microarray contained duplicate 70-mer oligonucleotide probes for 8,213 RHA1 genes, representing 89% of the predicted genes (17). The probes were designed and synthesized by Operon Biotechnologies, Inc. (Huntsville, AL).
Microarray spot intensities were quantified using Imagene 6.0 (BioDiscovery, Inc.). Expression ratios were normalized using GeneSpring version 7.2 (Silicon Genetics) by the intensity-dependent Lowess method, with 20% of the data used for smoothing. Average normalized expression ratios were calculated for each gene. Significant differential expression on propane versus pyruvate was defined as absolute ratios of
4.0 and Student's t test P < 0.05.
Details of the microarray design, transcriptomic experimental design, and transcriptomic data have been deposited in the NCBI Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo/) and are accessible through GEO Series accession no. GSE8480.
Quantification of gene expression by RT-qPCR.
The above extracted RNA was also used for RT-qPCR analysis. While trace genomic DNA contamination present after DNase/RNeasy mini column cleaning was acceptable for the cDNA-specific Cy labeling used in the microarray study, further removal of contaminating DNA was conducted prior to reverse transcription. This was accomplished through two more rounds of DNase I treatments using the DNA-free kit (Ambion) followed by an additional cleanup step in a RNeasy MinElute Cleanup kit (QIAGEN) to remove any other impurities. All treated RNA was stored at –80°C prior to further use. The transcripts of two genes from the PrMO operon (prmA and prmB) and one gene from the AlkMO operon (alkB) were selected for quantification. TaqMan primer-probe sets labeled with 6-carboxyfluorescein (FAM) using 6-carboxytetramethylrhodamine (TAMRA) as a quencher were purchased from PE Applied Biosystems (Foster City, CA). The primer-probe set for prmA, prmB, alkB, and DNA polymerase IV (DNA pol IV) genes are listed in Table 2. The DNA pol IV gene was used as an internal reference (housekeeping gene) for quantification. Primers and probes were designed for quantitative PCR using ABI Prism Primer Express Software (Applied Biosystems). For each design, sequence specificity was confirmed using the NCBI BLAST algorithm optimized for short nucleotide sequences on the GenBank database (www.ncbi.nlm.nih.gov).
|
View this table: [in a new window] |
TABLE 2. RT-qPCR primer and probe set used in this study
|
0.001 ng total RNA) and 0.5 µM of each reverse primer. For the second step of the RT-qPCR method, the reverse-transcribed samples were amplified on an ABI Prism 7000 sequence detection system (Applied Biosystems). The target and housekeeping genes were quantified in triplicate. Each 25-µl qPCR volume contained 2 µl of the reverse-transcribed RNA samples, 12.5 µl of 2xTaqMan universal PCR master mix (Applied Biosystems), 0.2 µM of probe, and 0.7 µM of each primer (forward and reverse). Thermocycling conditions were as follows: 2 min at 50°C, 10 min at 95°C, and 40 cycles of 15 s at 95°C and 1 min at 60°C. Differential expression was calculated by the Pfaffl method (22), which takes into account the amplification efficiency of qPCR for each target gene. The mass of DNA per volume was quantified using a NanoDrop ND-1000 spectrophotometer (Wilmington, DE), according to the manufacturer's instructions.
Concentrated plasmid DNA standards were synthesized by cloning separately the gene cluster for PrMO and AlkMO into E. coli TG1/pBS(Kan) (4). This resulted in the creation of E. coli TG1/pBS(Kan)PrMO.RHA1 and E. coli TG1/pBS(Kan)AlkMO.RHA1 (Table 1). The 4.3-kb PrMO gene cluster containing prmA was PCR amplified from chromosomal DNA of Rhodococcus sp. strain RHA1 using front primer R.PrMO.f.pBS and rear primer R.PrMO.r.pBS (Table 3). The front primer introduced a restriction site, a new ribosome binding site, a stop codon for the upstream lacZ
gene, and an altered start codon for the first gene in the RHA1 PrMO gene cluster; the rear primer introduced an alternate restriction site downstream of the last gene in the cluster. PCR products were gel extracted prior to restriction digestion, and plasmid pBS(Kan) was dephosphorylated using Antarctic phosphatase (New England Biolabs) after digestion. Analogous methods and design were used for constructing pBS(Kan)AlkMO-RHA1, with the 2.9-kb AlkMO cluster containing alkB, except front primer R.AlkMO.f.pBS and rear primer R.AlkMO.r.pBS were used (Table 3). The plasmid inserts were both confirmed by DNA sequencing.
|
View this table: [in a new window] |
TABLE 3. Cloning, gene deletion, and knockout confirmation screening PCR primers used in this study
|
by electroporation, verified by PCR, and then transformed into S17 E. coli competent donor cells (27) maintained in 25 µg ml–1 kanamycin (Table 1). Conjunctive plasmid transfer was achieved by coculturing the donor and Rhodococcus sp. strain RHA1 on selective LB peptone plates amended with 30 µg ml–1 nalidixic acid and 50 µg ml–1 kanamycin followed by sacB counterselection. Final confirmation of the removal of the target gene in kanamycin-sensitive, sucrose-resistant colonies was verified by colony PCR using a pair of primers that matched sequences flanking the target gene (Table 3). |
|
|---|
![]() View larger version (12K): [in a new window] |
FIG. 1. Constitutive removal of NDMA occurs at a fraction of the propane-induced rate. A Monod kinetic model (line) is fit to NDMA removal rates measured after RHA1 was grown on propane ( ). The constitutive NDMA degradation rate at 200 µg/liter ( ) represents average removal after independent growth on three noninducing substrates (pyruvate, LB medium, or soy broth).
|
Effect of propane on gene expression.
The increased NDMA degradation rates observed after growth on propane were explored by investigating the effect of propane on gene expression. Specifically, we employed a microarray with probes for 8,213 of 9,225 predicted genes of Rhodococcus sp. strain RHA1 to identify global transcriptional differences between triplicate batches grown on propane versus those grown on pyruvate. Table 4 lists genes with significant differential expression defined as absolute expression ratios no less than 4.0 and with 95% significance (P < 0.05) according to Student's t test.
|
View this table: [in a new window] |
TABLE 4. Genes up-regulated or down-regulated for growth on propane relative to pyruvate
|
![]() View larger version (13K): [in a new window] |
FIG. 2. Putative operon containing prm genes. Gene annotations are available in Table 4. Hatched arrows represent genes that had significant expression ratios, as determined by microarray and/or RT-qPCR, indicating up-regulation on propane.
|
Finally, significant differential expression was not demonstrated for a cluster of genes encoding a putative AlkMO (ro02534 to ro02538), despite its predicted functional similarity to PrMO (16).
RT-qPCR analysis of the expression of PrMO and AlkMO genes.
To better quantify the differential expression of genes encoding PrMO and AlkMO, we employed RT-qPCR. Specifically, three genes were targeted: prmA, encoding the large hydroxylase subunit of PrMO; prmB, encoding a reductase component of PrMO not represented by a probe on the microarray; and alkB, encoding the large subunit of AlkMO. RNA extracts used for the prior microarray experiment were also used for RT-qPCR. The relative amount of each gene transcript was normalized to that of a housekeeping gene coding for polymerase IV (DNA pol IV).
As shown in Fig. 3, the RT-qPCR results are generally consistent with those from the microarray. Using the Pfaffl method (22) of relative quantification, the prmA and prmB genes of the PrMO had propane/pyruvate expression ratios of 2,450 and 3,020, respectively. Conversely, the alkB gene had a much lower expression ratio of 3.2. The levels of expression of each of these genes were significantly different on the two substrates, as determined by Student's t test (
= 0.05 and n = 9). For prmA, the RT-qPCR expression ratio is more than 10-fold greater than the microarray value. This is not unusual for such highly expressed genes and probably indicates that the change in expression exceeded the dynamic range of the microarray analysis (8). The results for prmB confirm that, like the other genes in the putative PrMO operon, it is also up-regulated on propane (Fig. 2).
![]() View larger version (9K): [in a new window] |
FIG. 3. Effect of propane on transcription of three aliphatic monooxygenase components as quantified by RT-qPCR ( ; Pfaffl method) and spotted microarray ( ). The microarray did not code for prmB (NP), preventing its quantification by that method. Asterisks denote that the propane-grown values are statistically different from those of the pyruvate-grown controls based upon Student's t test (P < 0.05; n = 9 analytical replicates for RT-qPCR and n = 6 for microarray).
|
The constitutive removal of NDMA (Fig. 1) enabled screening of these knockout mutants. Parallel batch cultures of wild-type RHA1 and the two mutants were grown in liquid LB medium and harvested in the late exponential phase of growth. The cells were washed, suspended in phosphate buffer containing 200 µg NDMA liter–1, and assayed for NDMA removal (Fig. 4). In less than 4 h, both the wild-type strain and the alkB mutant removed NDMA to below detection limits. In contrast, NDMA removal by the prmA mutant was indistinguishable from that of the abiotic control. After 19 h, an additional sample was analyzed and still no significant biological NDMA removal by the prmA mutant was detected (not shown).
![]() View larger version (10K): [in a new window] |
FIG. 4. Excision of the prmA component of Rhodococcus sp. strain RHA1 eliminates this strain's ability to transform NDMA. , wild-type RHA1; , alkB knockout (RHA028); , prmA knockout (RHA027); and , no-cell control. LB-grown cells were suspended in phosphate buffer (amended with 200 µg NDMA liter–1) to cellular densities of 510, 530, 730, and 0 mg protein liter–1, respectively. Error bars portray the mean deviation of biological replicates.
|
|
|
|---|
The prm genes of RHA1 are part of a gene cluster (Fig. 2) which appears to be conserved in other actinobacteria. The first eight genes of this cluster have homologs in Gordonia sp. strain TY-5 (13) and in Mycobacterium smegmatis (GenBank accession no. NC008596). The orders of the genes are identical in the three organisms, and the encoded proteins of RHA1 are 64% to 93% identical to their homologs in the other two organisms. The prmABCD genes in TY-5 are part of a cotranscript induced by propane. Knockout mutagenesis in TY-5 showed that prmB and the seventh gene in the cluster, which we named prmE, are involved in propane catabolism. The prmE gene appears to encode a secondary alcohol dehydrogenase that catalyzes the second step of the catabolic pathway.
Our results provide both transcriptomic and phenotypic evidence for the involvement of the annotated PrMO in NDMA biotransformation. Combined oligonucleotide microarray and RT-qPCR demonstrate that transcripts from the prm gene cluster increased by orders of magnitude following growth on propane (Table 4 and Fig. 3). Next, partial excision of prmA prevented an otherwise genetically identical bacterium from growing on propane and eliminated both its constitutive and induced capability to degrade NDMA (Fig. 4). Despite the presumptive functional similarity between PrMO and AlkMO, deletion of alkB, encoding the large catalytic subunit of the latter enzyme, had no appreciable effect on either growth on propane or removal of NDMA (Fig. 4). Accordingly, the alk operon was not significantly up-regulated during growth (Table 4 and Fig. 3). These findings are consistent with observations of Nocardioides in which degradation of C2-to-C16 n-alkanes was the result of two distinct systems which included one alkane hydroxylase with a homolog to alkB (65% nucleotide identity by pairwise alignment) that was active on alkanes larger than C6 (10). However, a survey of alkB expression in three strains of Mycobacterium austroafricanum has shown that expression of this gene can correlate with the transformation of smaller gaseous alkanes, including propane (16). NDMA removal after growth on propane and the presumed expression of AlkMO was observed in one of these strains, Mycobacterium austroafricanum JOB5 (J. O. Sharp, C. M. Sales, and L. Alvarez-Cohen, unpublished data); however, the strain did not share RHA1's ability to constitutively remove NDMA. Though environmental NDMA degradation through induction on propane appears to extend beyond homologues of prmA, degradation and gene expression in systems not exposed to propane are less well understood.
Since the phenotype of the prmA deletion strain indicates that PrMO is solely responsible for both propane and NDMA oxidation in RHA1, the 20- to 30-fold up-regulation of etb genes on propane was surprising. Interestingly, EtbA has been implicated in the transformation of polychlorinated biphenyls, biphenyl, and ethylbenzene (8). Thus it is possible that propane, a comparatively inexpensive, benign, and mobile carbon source could serve as an effective alternative to ethylbenzene or biphenyl to induce expression of EtbA in environmental settings. It has been suggested (17) that large genomes with multiple broad-specificity catabolic enzymes such as those reported in strain RHA1 could have a competitive advantage in constantly changing soil environments. Such metabolic diversity could result in bacteria that can sustain growth by simultaneously metabolizing an array of compounds present in trace quantities. The coactivation of multiple oxygenase enzymes, while a surprising allocation of biochemical resources, could contribute to such a strategy.
Rhodococci and other members of the Actinomycetales are common soil bacteria. Given the involvement of PrMO in NDMA degradation and the previously discussed identification of similar activity in related strains, quantification of genes such as prmA in uncharacterized communities could provide a proxy for NDMA transformation potential. Furthermore, PrMOs have been reported to degrade a diverse array of organic compounds, including chlorinated C1-to-C6 alkanes, vinyl chloride, chlorinated ethylenes, methyl and ethyl tert-butyl ether, and tert-amyl methyl ether (26, 28), suggesting broader applications.
Interestingly, a correlation between desiccation-induced cell stress and induction of the prm operon in RHA1 was previously observed (15). The reason for up-regulation of prmA under these conditions is not obvious. However, other genes in the operon, such as groEL encoding a chaperone protein, may be part of a general stress response. Due to this response, stressed RHA1 cells, such as those likely to occur in a subsurface vadose zone experiencing alternating wet and dry cycles, varying oxygen content, or periods of growth and starvation, could have increased activity toward low-concentration environmental contaminates such as NDMA. While it is unclear how common this feature of prmA regulation is, it is possible that attenuation strategies involving stressed biomass could hold promise for remediating aquifers containing analogous micropollutants without the introduction of exogenous inducers.
|
|
|---|
Plasmid vector pK18mobsacB was provided by the van der Geize laboratory (Haren, The Netherlands). Aid in the form of microbial and analytical assistance was supplied by Gordon Stewart, Shaily Mahendra, Rebecca Davis, Kristin Robrock, and Wenjin Liu. Two anonymous reviewers provided helpful comments and criticisms during the submission process.
Published ahead of print on 14 September 2007. ![]()
|
|
|---|
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»