Previous Article | Next Article ![]()
Applied and Environmental Microbiology, August 2004, p. 4961-4970, Vol. 70, No. 8
0099-2240/04/$08.00+0 DOI: 10.1128/AEM.70.8.4961-4970.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
Center for Microbial Ecology,1 Department of Civil and Environmental Engineering, Michigan State University, East Lansing, Michigan 48824,5 Laboratory of Microbial Ecology and Technology, Ghent University, Ghent, Belgium,2 Chemical Engineering Department, University of Michigan, Ann Arbor, Michigan 48109,3 Xeotron Corporation, Houston, Texas 770544
Received 15 January 2004/ Accepted 22 April 2004
|
|
|---|
|
|
|---|
In collaboration with our laboratory, a draft sequence of LB400's genome was generated by The Joint Genome Institute. The latest annotated assembly (JGI/ORNL annotation December 2003) suggests that the total genome size is
9.7 Mbp (50 major contigs) and that there are 9,851 open reading frames (ORFs). Despite years of research on this PCB-degrading organism, very little is known about its general physiology, metabolism, biochemistry, and management of the metabolic warehouse of the large genome. To enhance our knowledge of this organism and, more specifically, how its genomic context allows it to be successful as a PCB degrader, a microarray containing one 45-mer oligonucleotide probe per ORF was constructed to explore genes that are potentially important in the biodegradation of PCBs. In collaboration with Xeotron Corporation (Houston, Tex.), we evaluated and used their new light-directed, in situ-synthesized microfluidic microarray platform (XeoChip) (14) for this purpose. We investigated the whole-genome expression patterns when LB400 was grown on biphenyl or benzoate compared to the patterns when the organism was grown on succinate. Furthermore, we investigated the expression patterns in the early stationary phase of biphenyl-grown cells. We validated our array results by using the data for a subset of genes measured by quantitative reverse transcription PCR (Q-RT-PCR).
|
|
|---|
Media and growth conditions.
LB400 was grown in liquid K1 mineral medium (46) supplemented with succinate (10 mM), benzoate (5 mM), or particulate biphenyl (5 mM; S = 6.99 x 103 g/liter) as the sole carbon source. Cells were grown at 29 ± 1°C in 200-ml batch cultures in 500-ml Erlenmeyer flasks on a rotary shaker at 250 rpm. At the start of each experiment, a glycerol stock culture was grown on R2A (Difco) agar plates, and biomass from pooled colonies was grown on a K1 agar plate containing biphenyl. K1 liquid medium containing succinate (reference conditions, transfer 1) was inoculated with pooled colonies from this plate and grown until the stationary phase. Two subsequent transfers were made in succinate medium prior to harvest (reference conditions, transfer 3). Cells from transfer 2 in K1 medium containing succinate were used to inoculate liquid K1 medium containing biphenyl and liquid K1 medium containing benzoate (treatment conditions, transfer 1), and cells were harvested after two more transfers in K1 medium containing biphenyl and liquid K1 medium containing benzoate (treatment conditions, transfer 3), respectively. Transfer 3 for each carbon source (reference and treatment) was done in duplicate to obtain biological replicates. The two biological replicates of succinate-, benzoate-, and biphenyl-grown cells were harvested at the maximum growth rates (0.40, 0.27, and 0.23 cell division/h, respectively) at optical densities of 0.45, 0.35, and 0.35, respectively. Additionally, two biological replicates of biphenyl-grown cells were harvested within 1-log-phase generation time after the start of the stationary phase, at an optical density of
0.75. We termed this sampling period the early stationary phase because the growth rate decreased to a much lower value after this point although the cells eventually reached an optical density of
1.2.
RNA extraction and labeling.
RNAlater (Ambion) was added at a 1:1 ratio to a culture to protect RNA against degradation. Bacterial cells were harvested by centrifugation at 5,000 x g for 10 min at 4°C. The RNA was extracted with an RNeasy RNA extraction kit (QIAGEN), and the remaining DNA was removed by 30 min of incubation at room temperature with 1.5 U of DNase I (10 U/µl; Roche) per µg of nucleic acid. The integrity of the RNA and the absence of DNA were verified by 1.2% agarose gel electrophoresis.
For direct incorporation (succinate versus early-stationary-phase biphenyl), 25 µg of total RNA was labeled by using the Xeotron cDNA labeling protocol (www.xeotron.com), except that 30 µg of random hexamers (Invitrogen) was used instead of poly(dT) primers. Reaction cleanup was performed by using the Qiaquick PCR purification system protocol (QIAGEN). Amino-allyl labeling (all other hybridizations) was performed by using a protocol adapted from a protocol of The Institute for Genomic Research (http://www.tigr.org/tdb/microarray/protocolsTIGR.shtml). Briefly, 2 µg of total RNA was labeled overnight at 45°C by using 6 µg of random primers (Invitrogen), amino-allyl-labeled dUTP (Sigma), and Superscript III reverse transcriptase (Invitrogen), and subsequently reactive Cy5 or Cy3 fluorophores (Amersham) were coupled to the amino-allyl groups. Purification after enzymatic incorporation and chemical coupling was performed by using QiaQuick PCR purification columns (QIAGEN) as described in The Institute for Genomic Research protocol. The quantity of labeled cDNA and the fluororophore incorporation efficiency were determined by using UV-visible spectrophotometry.
Microarray design and synthesis.
The LB400 genomic microarray consists of one unique 45-mer (specificity checked by BLAST versus the LB400 genome) for each of the 9,511 annotated ORFs (JGI/ORNL annotation December 2000b), designed by using OligoArray (30), with a melting temperature of 86 to 92°C and a G+C content of 50 to 60% and subsequently synthesized in situ by Xeotron by using their 4,000-chamber (spot) format. The genome is represented by a genomic set consisting of three chips, each containing
4,000 synthesized 45-mers. A total of 1,601 probes on chip 1 and chip 3 of the genomic set are identical and used as a chip-to-chip replication control in each experiment. For comparison tests of the XeoChips, we used a spotted 70-mer oligonucleotide microarray that contained triplicate spots of 507 probes (the probe for a gene did not necessarily overlap with the XeoChip probe for that gene) targeting a random set of the genes of LB400 (J. Park, T. V. Tsoi, S. A. Hashsham, J. R. Cole, S. Callister, and J. M. Tiedje, Abstr. NIEHS Symp. Bioremediat. Biodegradat., abstr. 23, 2002). Production was based on optimized protocols described previously (9).
Hybridization, scanning, and data analysis.
All buffers were passed through a 0.22-µm-pore-size filter to keep particulate matter out of the microfluidic channels. XeoChips were placed in Xeotron hybridization adaptors and prehybridized with a solution containing 25% formamide (Ambion), 6x SSPE (pH 6.6; Ambion), and 1 mg of acetylated bovine serum albumin per ml (Ambion) for 5 min at a flow rate of 100 µl/min (1x SSPE is 0.18 M NaCl, 10 mM NaH2PO4, and 1 mM EDTA). A 100-µl hybridization solution consisting of 25% formamide, 6x SSPE (pH 6.6), 1 mg of bovine serum albumin per ml, and labeled target containing 100 pmol of each dye was denatured at 95°C for 3 min, snap cooled on ice, filtered by using Spin-X 0.22-µm-pore-size spin column filters (Corning), and hybridized by circulating the hybridization mixture for 18 h at 32°C at a flow rate of 100 µl/min. Both solutions were stripped of dissolved gases by flushing them with nitrogen for 2 min prior to the denaturation step. After hybridization, each chip was washed at a flow rate of 100 µl/min with 6x SSPE (pH 6.6) for 5 min at 32°C and subsequently with 1x SSPE (pH 6.6) for 5 min at 32°C. The spotted 70-mer microarray was hybridized as described previously (9).
Fluorescence intensity data were acquired by using an Axon 4000A laser scanner (Axon Laboratories) and, for XeoChips, were subsequently imported in ArrayPro 4.0 (Media Cybernetics, L.P.) to extract the 635-nm (Cy5) and 532-nm (Cy3) signals. Based on an analysis by Fang et al. (13), signals without background subtraction were used for both platforms. Spotted microarray data were extracted by using Genepix Pro 3.0 (Axon Laboratories). Median signals for each channel (Cy5 and Cy3) were imported into GeneSpring 5.0 (Silicon Genetics) and normalized by using Lowess intensity-dependent normalization. P values were calculated by using GeneSpring's cross gene error model and the Student t test algorithm based on the variation between log2 ratio values of biological replicates. Plots of M [log2 (Cy5/Cy3)] versus A [log2 (
(Cy5 x Cy3)] were drawn as described by Dudoit el al (11) by using all data from one experiment grouped per chip. Chip-to-chip reproducibility was determined by calculating the coefficient of variation (CV) (CV = [standard deviation/mean] x 100) for the ratios for identical probes on separate microarrays (i.e., separate slides for spotted microarrays [1,536 probes] and chip 1 and chip 3 of a genomic set for XeoChips [1,601 probes]). The CV for biological reproducibility for XeoChip experiments was calculated separately for ratios of hybridizations with dye-swapped labeled biological replicate samples for all probes (9,511 probes), all nonreplicated probes (7,910 probes), and all replicated probes (1,601 probes).
Quantitative real-time PCR.
Triplicate Q-RT-PCR runs were performed for 25 genes by using the same RNA samples that were used for the microarray hybridizations. One microgram of total RNA was converted into cDNA, and 1/500 to 1/1,000 of the sample was utilized for a 40-cycle, two-step PCR with an ABI 7900HT (Applied Biosystems, Foster City, Calif.) by using 1x SYBR Green master mixture (Applied Biosystems) and each primer at a concentration of 125 nM (sequences available on request). Amplicon size (80 to 100 bp) and reaction specificity were confirmed by agarose gel electrophoresis and product dissociation curves. The number of target copies in each sample was interpolated from the detection threshold value by using a purified PCR product standard curve for bphA, which was constructed for each Q-RT-PCR run. l6S rRNA expression was measured as an internal control, and the measured internal control signal was used to normalize variations due to different reverse transcription efficiencies.
|
|
|---|
![]() View larger version (72K): [in a new window] |
FIG. 1. (a) Plots of M [log2 (Cy5/Cy3)] versus A [log2 ((Cy5 x Cy3)] per chip for a typical XeoChip hybridization, after Lowess normalization. (b) Image of a XeoChip (4,000-chamber [spot] platform) after hybridization. (c) Chip-to-chip reproducibility and biological reproducibility of XeoChip platform determined by using data from all hybridizations: frequency (Freq) and cumulative frequency (Cumm Freq) of the CVs of ratios between chips hybridized with same hybridization mixture (chip to chip), of ratios between nonreplicated probes of biological dye-swapped replicates (biological), and of ratios on spotted microarrays hybridized with the same hybridization mixture (spotted). The error bars indicate standard errors between experiments.
|
![]() View larger version (26K): [in a new window] |
FIG. 2. Evaluation of microarray results by Q-RT-PCR. (a) Log2 ratio of microarray versus log2 ratio of Q-RT-PCR for each experiment, averaged over biological replications. (b) Bias between microarray and Q-RT-PCR data. The ratios determined by Q-RT-PCR [log10 (ratio-QRTPCR)] are plotted against the ratio of ratios obtained by using microarrays and Q-RT-PCR [log10 (ratio-array)/(ratio-QRTPCR)]. The error bars indicate one standard deviation for biological replicates. Benzo, benzoate-grown cells; BiphML, mid-log-phase biphenyl-grown cells; BiphES, early-stationary-phase biphenyl-grown cells.
|
45% of the genes under either condition. Functional classification of differentially expressed genes based on clusters of orthologous groups of proteins (COGs) (Table 1) indicated resource reallocations depending on the growth substrate and phase. |
View this table: [in a new window] |
TABLE 1. Summary of differential expression based on COG protein classification
|
![]() View larger version (50K): [in a new window] |
FIG. 3. Biphenyl and benzoate metabolic pathways. Genes, whose annotations were manually verified by using BLAST, are presented as they are organized on the chromosome (position, orientation, and scaffold) for catabolic genes (dark grey arrows), putative regulatory genes (light grey arrows), and ORFs with unknown functions (open arrows). The ORF numbers in LB400 corresponding to the catabolic genes are indicated in the pathways. The graphs show the M values [log2 (test condition/succinate)] for benzoate-grown cells (open bars), log-phase biphenyl-grown cells (grey bars), and early-stationary-phase biphenyl-grown cells (black bars). The error bars indicate the standard errors (standard deviation/n) for the biological replications. Genes without data were not included in the sequence draft and are therefore not present on the microarray, while two neighboring ORFs were annotated as bphI due to replication during genome assembly. log, log phase; stat, stationary phase.
|
Growth phase-dependent metabolism.
The most striking observation for the expression patterns in the growth phase is the up-regulation of genes coding for a homolog of methanol dehydrogenase (xoxF), biosynthesis of its (and other quinoproteins') cofactor pyrroloquinoline quinone (PQQ), and a tetrahydromethanopterin-dependent oxidative pathway for conversion of its product, formaldehyde, to CO2 and further conversion to acetyl-CoA by carbon monoxide dehydrogenase (Fig. 4). Additionally, up-regulation of the chlorocatechol pathway (Table 2) was observed, which was not detected during exponential growth. When cells were grown on benzoate, however, the chlorocatechol pathway (as well as the 2-aminophenol pathway) was significantly (P < 0.05) up-regulated relative to the expression during growth with succinate (Table 2).
![]() View larger version (44K): [in a new window] |
FIG. 4. C1 metabolic pathways induced in early-stationary-phase biphenyl-grown cells and their organization on the chromosome, represented by squares left and right of the center based on their orientation. Squares representing genes that are significantly up-regulated in early-stationary-phase biphenyl-grown cells are black. The plot shows M values and standard errors for biological replicates. The solid lines in the pathway indicate up-regulated steps. The gene names (annotations were manually verified by BLAST) are those found in GenBank, except for the following names: gsh-faldh (glutathione-dependent formaldehyde dehydrogenase), fgsh-hydrol (formylglutathione hydrolase), M.e. (Methylobacterium extorquens), tr reg (transcriptional regulator), deHase (dehydrogenase), etoh deHase and etdh (ethanol dehydrogenase), and acdh II (type II acetaldehyde dehydrogenase). H4MPT, tetrahydromethanopterin.
|
|
View this table: [in a new window] |
TABLE 2. Expression ratios for 2-aminophenol and 3-chlorocatechol pathways (annotation manually verified by BLAST) as measured for benzoate-grown, mid-log-phase biphenyl-grown, and early-stationary-phase biphenyl-grown cells relative to succinate-grown cells
|
54, a
54 modulation protein, and 2 of 13
54-dependent transcriptional regulators were up-regulated (Table 4). |
View this table: [in a new window] |
TABLE 3. Summary of differential expression, expressed as ratios relative to succinate expression, for transport and regulation specific to biphenyl or benzoate (annotation manually verified by using BLAST)
|
|
View this table: [in a new window] |
TABLE 4. 54 transcription factor and its dependent transcriptional regulators, based on a GAFTGA amino acid motif searcha
|
|
|
|---|
30%) suggest that there is more to efficiently degrading biphenyl than expressing the biphenyl pathway genes.
By linking our data to previous reports, we can propose a refinement of the regulatory network of the biphenyl pathway in LB400. The observed equal expression of ORF0 on benzoate and on succinate and its up-regulation on biphenyl confirmed Erickson and Mondello's (12) report, based on S1 nuclease mapping, of a biphenyl-inducible promoter upstream of ORF0. These authors also suggested that there are two constitutive
70-dependent promoters upstream of bphA1. However, mutational studies and Northern blot analyses performed by Beltrametti et al. (2), who proposed that ORF0 is a positive regulator of bphA1A2, and our observation of increased expression ratios of the bph genes with increased ORF0 expression indicate that there is a positive regulatory function of ORF0 induced by biphenyl. The similarity in expression patterns among ORF0, bphJ, and bphI and possibly bphD suggests that ORF0 is an inducer of bphJ, bphI, and bphD. This is similar in part to what has been suggested for Pseudomonas pseudoalcaligenes KF707, whose bph genes exhibit 98 to 100% sequence similarity to LB400 bph genes (BLAST), except for ORF0 (86%), and are organized just like the LB400 bph genes (41). In KF707, bphR1 (formerly ORF0) was identified as a transcriptional regulator for itself and the lower biphenyl pathway genes, as well as bphD, while bphR2, which is not located near the bph cluster, is a necessary LysR-type regulator controlling all bph genes (40). Although the protein sequences are not significantly similar (15%), a LysR-type transcriptional regulator of LB400 (ORF 10103), located elsewhere in the genome, was up-regulated under both biphenyl growth conditions, possibly acting in the role of bphR2.
Based on sequence information, LB400 has two candidate benzoate degradation pathways, one of which, an aerobic degradation pathway with CoA activation (15), is present in two copies, which exhibit 70 to 86% nucleotide sequence similarity (ClustalW). The different pathway choices for benzoate degradation, depending on the carbon source added, and specifically the induction of at least one gene cluster encoding the CoA pathway during growth on biphenyl indicate the presence of a carbon source-specific regulatory network. Apart from pcaR, a positive regulator of the ß-ketoadipate pathway, the microarray failed to detect the Q-RT-PCR-measured up-regulation of the catechol-ß-ketoadipate pathway in biphenyl-grown cells compared to succinate-grown cells due to the pathway's background-level expression. Based on the Q-RT-PCR data, which indicated that there was
1,000-fold-lower expression of this pathway in biphenyl-grown cells than in benzoate-grown cells, and the knowledge that one of the pathway's intermediates, cis,cis-muconate, positively feeds back, inducing catA expression (19), we hypothesize that catechol-ß-ketoadipate pathway expression is positively correlated with the benzoate concentration in the medium, which can be assumed to be significantly lower when cells are growing on biphenyl than when cells are growing on the same concentration of benzoate. Instead of being induced by benzoate, CoA pathway regulation could be dependent on biphenyl, one of the metabolic intermediates, or a lower oxygen concentration due to upper Bph pathway dioxygenase (BphA/BphC) oxygen consumption. Different promoter regions and transcriptional regulators, which are neighbors of the two CoA pathway copies, offer a potential explanation for the differential expression of the two copies. The substrate specificities of the enzymes, especially for chlorinated benzoates, remain to be determined. Also, the LysR-type regulator (ORF 10103) up-regulated under both biphenyl growth conditions and an ArsR-type regulator (ORF 12464) that is specifically up-regulated in benzoate-grown cells are targets for further research on regulation of these pathways.
As expected due to the growth phase transition, the largest amount of differentially expressed genes was observed in the early-stationary-phase sample (36, 37, 39). Remarkably, a large number of the up-regulated genes were involved in C1 metabolism, specifically the oxidation of methanol (or possibly a methoxylated carbon source) to CO2 (38) and possibly further acetogenic assimilation to acetyl-CoA (10). Burkholderia species have been shown to assimilate C1 compounds (26, 28), and considering the ubiquitous presence of C1 compounds from microbial degradation of plant aromatic compounds present in the rhizosphere, one of the typical niches of Burkholderia species (7), the presence of these pathways in LB400 could provide an ecological advantage. As for regulation, the coexpression of
54, which is involved in regulation of a wide range of molecular processes (6), as well as two
54-dependent transcriptional regulators, which are located next to PQQ biosynthesis genes and tetrahydromethanopterin-dependent formaldehyde oxidation pathway genes, suggests that this metabolic shift to C1 metabolism is
54 dependent. The two transcriptional regulators were previously linked to mxaF (35), the homolog of which (xoxF) has several potential
54-dependent promoters in LB400 (data not shown). A last element in the regulation of this metabolic switch is the up-regulated two-component mxbM-mxbD regulatory system, which has been reported to regulate both PQQ biosynthesis and methanol dehydrogenase (16).
As described above, the most significant change in metabolism-related expression was the activation of a methylotrophic metabolic pathway. No involvement of the methylotrophic enzymes with biphenyl or PCB degradation has been reported thus far. The activation of C1 metabolic pathways, as well as the multitude of genes with transport annotation, could be an indication of LB400 switching to a scavenging mode and therefore turning to less energy-efficient, although more available, carbon sources. The wide diversity of up-regulated functions in the cells, the distribution of which is highly similar to the genome's contents, could be interpreted as another indication of this switch to scavenging. It is interesting that the C1 pathway coexpression with the core PCB degradation pathway fits well with the hypothesis that PCB degradation pathways evolved from (often methoxylated) plant aromatic degradation pathways (34).
In concert with this work we documented the performance of a new type of microarray especially suited to microbiological applications. Based on M-A plots, which have been shown to be useful for detecting artifacts and the signal dependency of ratios and for general normalization quality evaluation (11, 43), GeneSpring's Lowess normalization of data proved to be sufficient when the new microarrays were used. The average chip-to-chip CV compared favorably to the values reported for cDNA microarrays (10 to 25%) (44), as well as 70-mer oligonucleotide microarrays which we tested ourselves. A similar correlation was observed between significantly differentially expressed genes and ink jet technology arrays (r = 0.98) (23). High reproducibility and chip hybridization homogeneity due to forced flow of the hybridization mixture resulted in more reliable data, with less replication. The most prominent advantage of these chips lies in their design flexibility, which allows easy, fast, and low-cost alterations to the chip design simply by providing a new spreadsheet of probe sequences and locations to the in situ chip synthesizer (14). This should be useful when genome closure is attained and annotation updates are made, and it also offers opportunities for microarrays that target more complex biological problems, where unpredictability of bacterial diversity requires iterative probe design capability.
The qualitative correspondence of most of the evaluated microarray results with Q-RT-PCR results indicated the reliability of XeoChips for identifying differentially expressed genes. Sixty percent of the noncorrespondence could be explained by low signal detection at the probes, due to low transcript levels and the relative insensitivity of the microarray assay. Importantly, only one false positive (type I error) in 75 comparisons was observed, while the remainder of noncorresponding results were due to false negatives compared to the more sensitive Q-RT-PCR. These results make sense in the context of the Q-RT-PCR-versus-microarray bias (i.e., underestimation of ratios by the microarray or overestimation by the Q-RT-PCR). This observation is consistent with a previous report (45) which showed that there is a bias for cDNA microarray and Q-RT-PCR data that follows a consistent trend.
This study confirmed that the new and flexible in situ-synthesized XeoChip is a reliable DNA microarray platform for bacterial genomics research. Discoveries made by using this approach, such as LB400's benzoate metabolism and early-stationary-phase C1 metabolism, when linked to physiological data, could result in new hypotheses regarding the underlying molecular mechanisms for improved biodegradation of recalcitrant pollutants.
We acknowledge Lindsay Eltis for suggestions on regulation of benzoate metabolism, Christopher Marx for discussions on C1 metabolism, and Benli Chai for bioinformatic support.
|
|
|---|
54 (
N) transcription factor. J. Bacteriol. 182:4129-4136.
54-dependent transcriptional activators. J. Bacteriol. 185:1757-1767.
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»