Previous Article | Next Article ![]()
Journal of Bacteriology, October 2008, p. 6697-6705, Vol. 190, No. 20
0021-9193/08/$08.00+0 doi:10.1128/JB.00543-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
and
Gary Stacey1,5,6*
National Center for Soybean Biotechnology, Division of Plant Sciences, University of Missouri, Columbia, Missouri 65211,1 Department of Statistics, University of Missouri, Columbia, Missouri 65211,2 Department of Soil, Water, and Climate, University of Minnesota, St. Paul, Minnesota 55108,3 Department of Biological Sciences, Wright State University, Dayton, Ohio,4 Division of Biochemistry,5 Department of Molecular Microbiology and Immunology, University of Missouri, Columbia, Missouri 652116
Received 21 April 2008/ Accepted 19 July 2008
|
|
|---|
|
|
|---|
Hydrogen oxidation in B. japonicum was extensively studied in strains USDA122 and USDA110 due to its potential for enhancing symbiotic nitrogen fixation efficiency. As a symbiont, B. japonicum USDA110 expresses an uptake hydrogenase that serves to recapture electrons lost via hydrogen evolution resulting from nitrogenase activity (15). Expression of the uptake hydrogenase system of B. japonicum also occurs during heterotrophic growth when an atmosphere containing 2% oxygen and 2% hydrogen is supplied (47). Under these conditions no ribulose 1,5-bisphosphate carboxylase/oxygenase (RuBisCO) activity is detected (47). In addition, the uptake hydrogenase is expressed when B. japonicum is cultured chemoautotrophically under microaerobic conditions (1% oxygen) when supplied with hydrogen (10%) and carbon dioxide (5%) in a minimal medium lacking an organic carbon source (33). RuBisCO activity is present under these conditions at levels sufficient to account for all assimilated carbon (33).
The recent availability of Affymetrix (25) and oligonucleotide microarray (5) platforms for B. japonicum USDA110 provides useful tools to analyze the environmental responses of the bacterium on a global scale. Transcriptomic studies of B. japonicum cultured in minimal and rich media (5), in the bacteroid state (5), and responding to osmotic stress (5) and desiccation (8) as well as iron limitation (41, 51) have all been reported. In addition, analysis of the NifA and
54 regulons from anaerobically cultured cells was recently completed (25). Here we extend the list of transcriptomic studies by generating mRNA profiles of B. japonicum USDA110 grown chemoautotrophically to better characterize the transcriptional and physiological responses to an alternative mode of growth based upon hydrogen oxidation and carbon dioxide fixation.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Bacterial strains and plasmids used in this study
|
Cells were grown for RNA extraction in 2-liter Fernbach flasks containing 500 ml of minimal medium. Inoculum was grown in 10 ml of PYE to late logarithmic phase and collected by centrifugation at 24°C and
9,000 x g. Cell pellets were gently resuspended with the same medium into which they would be inoculated (minimal medium with or without arabinose supplementation) to remove any residual PYE and centrifuged again. The cells were resuspended in 1 ml of growth medium and used to inoculate 500-ml cultures. Starting culture optical densities at 595 nm were typically 0.01 to 0.03. After inoculation, the cultures were sealed and flushed for 20 min with a gas mixture of 1% O2, 5% CO2, 10% H2, and 84% N2. Control flasks were sealed under an ambient air atmosphere. The flasks were statically incubated at 30°C and gently mixed twice daily until they were harvested after 7 days of incubation.
RNA isolation. Cells from 500-ml cultures were harvested by centrifugation for 10 min at 4°C and 6,000 x g in a fixed-angle rotor. The cell pellets were resuspended in 10 ml ice-cold medium, transferred to cold 15-ml conical tubes and centrifuged for 10 min at 4°C and 8,000 x g in a swinging bucket rotor. After discarding the supernatant, the cell pellets were immediately stored at –80°C. Total RNA was extracted using a hot phenol method (1). DNase treatment and RNA purification were carried out as previously described (5). RNA samples were electrophoresed in 0.8% agarose gels prior to use to ensure RNA integrity.
Microarray expression analysis. Transcriptional profiling of B. japonicum was performed using whole-genome 70-mer oligonucleotide microarrays containing probes specific to 8,453 predicted open reading frames printed by the Washington University Genome Sequencing Center (St. Louis, MO) and described by Chang et al. in 2007 (5). Synthesis of cDNA, labeling, and hybridization were completed as previously published (5). All hybridizations were performed with 2 µg of labeled cDNAs. Three independent biological replicates of cells grown heterotrophically, chemoautotrophically, and chemoautotrophically supplemented with arabinose were generated. Three independent comparisons were performed such that each culture condition was compared to the other two culture conditions. Three replicate hybridizations, including a dye swap for each replication, were performed for each comparison with the exception of the chemoautotrophy and heterotrophy comparison, for which two replications were performed.
Spot intensities were obtained from each slide using GenePix Pro 6.0 software (Molecular Devices Corp.). The base 2 logarithms of the mean foreground intensities from GenePix software were normalized by the joint lowess normalization method (7) using R/maanova software (0.98_7 version; Jackson Laboratory, Bar Harbor, ME; http://www.jax.org/staff/churchill/labsite/software/Rmaanova/index.html). No background subtraction was used since more than 20% of the spots would have negative values if background correction was applied. After normalization of each array, all 32 channels (2 channels per array) were median centered so that the median log-scale expression levels were all zero. The 1,280 empty and control spots were removed from the downstream analysis to reduce the dimension. A linear mixed effects model was applied to the normalized log-scale expression measures separately for each gene using R/maanova software. Each linear mixed effects model included array effect, dye effect, spot effect, biological replicate effect, and culture condition effect, where the array, spot, and biological replicate effects were treated as random effects. Since the spots are printed in duplicate on each array, spot effects were included to account for the technical replicates. t tests for three pairwise comparisons were conducted to identify significantly differentially expressed genes across any two culture conditions. P values from the t statistics were obtained via the Matest statement in the R/mmanova software, and a q-value (45) was computed for each P value to produce lists of differentially expressed genes with estimated false discovery rates of 5%.
qRT-PCR analysis. Microarray expression data were confirmed for a selected set of genes by quantitative reverse transcription-PCR (qRT-PCR). RNA was isolated, DNase treated, and purified as described above from cultures of two independent biological replicates. Three micrograms of RNA was reverse transcribed with 200 U of Moloney's murine leukemia virus reverse transcriptase (Promega, Madison, WI). Quantitative RT-PCR was performed with an ABI 7500 real-time PCR system with sequence detection system software version 1.3 (Applied Biosystems, Foster City, CA) utilizing Sybr green PCR master mix (Applied Biosystems, Foster City, CA). Gene-specific primers are listed in Table 2. PCR conditions were as follows: stage 1, 50°C for 2 min; stage 2, 95°C for 10 min; stage 3, 45 cycles of 95°C for 15 s and 60°C for 1 min; stage 4 (for dissociation), 95°C for 15 s, 60°C for 1 min, and 95°C for 15 s. Data obtained for two technical replicates of each biological replicate were normalized to the elongation factor fusA (bll5403), which is constitutively expressed under the given culture conditions.
|
View this table: [in a new window] |
TABLE 2. Gene-specific primers for quantitative RT-PCR
|
suicide vector (50), creating pKO
-ICL. Insert identity was confirmed by DNA sequence analysis. A kanamycin (Km) cassette from pHP45
-Km (16) was PCR amplified with the primers 5'-GCTGAACCGGTCCTAGGCCACTAACTAA-3' and 5'-CAGTAACCGGTTGATTGATTGAGCAAGC-3', containing AgeI sites (underlined) and subsequently inserted into an AgeI site in blr2455 to generate pKO
-ICL-Km. pKO
-ICL-Km was transferred from E. coli DH5
to B. japonicum USDA110 by triparental mating with the helper plasmid pRK2073 (32). Transconjugants were selected based upon kanamycin resistance and streptomycin sensitivity. Mutants were confirmed by colony PCR and Southern blot analysis (data not shown). Growth assays were performed with the mutant in minimal medium (5) buffered with 1.1 g liter–1 4-morpholinoethanesulfonic acid and 1.3 g liter–1 2-[4-(2-hydroxyethyl)-1-piperazine]ethanesulfonic acid and supplied with 20 mM acetate as a carbon source. The growth phenotype of the B. japonicum USDA110 aceA mutant, relative to the wild-type strain, was assessed using Biolog GN plates (Biolog, Inc., Hayward, CA). Wild-type and mutant strains were grown for 3 days at 30°C on R2A agar medium (Difco, Detroit, MI), and cells were suspended in Biolog GN/GP inoculation broth (Biolog, Inc.) to a final density of approximately 109 cells/ml. Cells (150 µl) were inoculated into the wells of Biolog GN plates and incubated for 24 h at 30°C, and color development in each well was determined at 540 nm using a microplate reader.
Microarray data accession numbers. All microarray data are available at the NCBI Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo) via GEO series accession numbers GSE10295, GSE10296, and GSE10298.
|
|
|---|
Comparison of the chemoautotrophic condition to the heterotrophic condition identified 1,485 genes as differentially expressed at a 1.5-fold cutoff (856 upregulated and 629 downregulated) (see Tables S1 and S2 in the supplemental material). Likewise, when the arabinose-supplemented chemoautotrophic condition was compared to the heterotrophic condition, 1,227 genes were identified as differentially expressed (641 upregulated and 586 downregulated) (see Tables S3 and S4 in the supplemental material). Substantial overlap was observed between the two comparisons; 601 upregulated genes and 528 downregulated genes were identical in the two data sets, suggesting that the chemoautotrophic and the arabinose-supplemented chemoautotrophic conditions are highly similar at the transcriptional level. In a direct comparison of the two chemoautotrophic conditions (i.e., the presence and absence of arabinose), no genes were identified as differentially expressed. These results indicate that the addition of arabinose to chemoautotrophic conditions does not influence gene expression under our growth conditions, and a mixotrophic mode of growth under microaerobic conditions in the presence of a carbon source was not observed.
Classification by functional category of the differentially expressed genes is summarized in Fig. 1. Genes annotated as hypothetical comprise about 50% of the differentially expressed genes in each data set. Differentially expressed genes belonging to the categories of cell envelope, cellular processes, central intermediary metabolism, and energy metabolism are more highly represented in cells cultured under chemoautotrophic conditions. Included in these categories are genes involved in expression of the uptake hydrogenase, the Calvin-Benson-Bassham (CBB) cycle, nitrogen fixation, ATP synthesis, and components of the respiratory chain and biosynthesis of pili and flagella, outer membrane proteins, and a subset of chaperones.
![]() View larger version (26K): [in a new window] |
FIG. 1. Functional classification of differentially expressed genes with transcript changes greater than 1.5-fold. (A) Chemoautotrophy and heterotrophy comparison. (B) Arabinose-supplemented chemoautotrophy and heterotrophy comparison. In both panels, black bars represent the chemoautotrophic condition and gray bars represent the heterotrophic condition. Functional classifications were derived from B. japonicum genome annotations available via Rhizobase (http://bacteria.kazusa.or.jp/rhizobase/) and the JCVI Comprehensive Microbial Resource (http://cmr.tigr.org/tigr-scripts/CMR/CmrHomePage.cgi).
|
Transcriptional data from chemoautotrophically cultured cells identified 18 of the 23 genes on the characterized uptake hydrogenase cluster as being upregulated (Table 3). Included are four of the most strongly upregulated genes, hupL and hupS, encoding the subunits of the uptake hydrogenase, hupC, a probable b-type cytochrome involved in electron transfer, and hupF, a hydrogenase accessory protein (20, 43). HoxA, a transcriptional activator of the NtrC-like family of two-component response regulators (13), was also upregulated greater than twofold. A deletion mutant of hoxA was demonstrated to have no free-living hydrogenase activity; however, bacteroids retained full hydrogenase activity, indicating differential regulation of free-living and symbiotic hydrogenase expression (13). Recently, Hauser et al. (25) identified a subset of hydrogenase genes, including hupK (bll6933), hupD (bll6939), hupC (bll6940), hupL (bll6941), and hupS (bll6942), as being NifA regulated under anaerobic conditions. Expression of the B. japonicum hydrogenase system is dependent on RpoN activity (2), and an RpoN binding site was identified upstream of hupS (2, 25). Both HoxA and NifA contain RpoN-interacting domains as identified by InterProScan (39). A possible interpretation of these results is that HoxA directly or indirectly induces expression of the B. japonicum uptake hydrogenase under free-living microaerobic conditions, while NifA may regulate the uptake hydrogenase under symbiotic conditions as previously suggested (13). Regulation of the uptake hydrogenase in Rhizobium leguminosarum bv. viciae was demonstrated to be under the control of NifA (3). R. leguminosarum bv. viciae does not contain a functional HoxA, and unlike B. japonicum, it does not express the uptake hydrogenase under free-living microaerobic conditions (3, 36). However, in B. japonicum no traditional NifA binding site was identified upstream of hupS (25), and binding of NifA to the uptake hydrogenase promoter region was not detected (12). Durmowicz and Maier provided evidence for a role of FixK2 in symbiotic hydrogenase expression (12). Mutants of fixJ, part of the FixLJ oxygen-responsive two-component system, and fixK2 lack hydrogenase activity under symbiotic conditions, but it is not clear if FixK2 acts directly to regulate hydrogenase expression (12). NifA and FixK2 are components of separate regulatory cascades (11), and it remains unclear how much cross talk exists between the two. At the present time, it is clear that hydrogenase regulation in free-living B. japonicum is subject to control by HoxA, although the exact means of symbiotic hydrogenase regulation remain unresolved.
|
View this table: [in a new window] |
TABLE 3. B. japonicum genes significantly upregulated during chemoautotrophic growtha
|
B. japonicum cultured heterotrophically derives energy from both substrate-level phosphorylation and oxidative phosphorylation. When grown chemoautotrophically, the primary means of energy generation is by oxidative phosphorylation requiring ATP synthase activity. Two ATP synthase complexes encoded by five genes each are found in the B. japonicum genome (bll0439 to bll0443 and bll1185 to bsl1189). As shown in Table 3, 8 of the 10 genes encoding ATP synthase components are induced during chemoautotrophic growth, demonstrating increased ATP generation by oxidative phosphorylation.
CO2 fixation and carbon metabolism. Facultative chemoautotrophs such as B. japonicum utilize the CBB reductive pentose phosphate cycle as the primary means of carbon assimilation during autotrophic growth (30). The B. japonicum genome contains a cluster of nine genes encoding the key components of the CBB cycle (bll2580 to blr2588) (27). Eight of the nine genes of this cluster are strongly upregulated during chemoautotrophic growth, as shown in Table 3. Included are two enzymes unique to this cycle (37, 44), RuBisCO and phosphoribulokinase, which are both strongly upregulated. The only gene of the cluster not upregulated is the divergently transcribed cbbR, encoding a LysR-type transcriptional regulator that is weakly constitutively expressed in another facultative autotroph, Ralstonia eutropha (44).
Carbon dioxide (or bicarbonate) is required not only for autotrophic growth but also heterotrophic growth of many bacteria (31). An enzymatic activity of importance to carbon dioxide metabolism is the interconversion of carbon dioxide and bicarbonate catalyzed by carbonate dehydratase/carbonic anhydrase. B. japonicum contains four putative carbonic anhydrases (bll1137, bll2065, bll4863, and bll4865) and one putative carbonate dehydratase (blr0500), of which two, bll1137 and bll4865, are induced under chemoautotrophic conditions (Table 3). The first, bll1137, is a putative
-type carbonic anhydrase with strong similarity (67% amino acid identity) to a characterized periplasmic
-type carbonic anhydrase designated acaP from Rhodopseudomonas palustris (38). R. palustris acaP is expressed in response to anaerobic conditions and functions to enhance the extracellular conversion of HCO3– to CO2 at neutral pH, where spontaneous conversion is slow (38). The second gene, bll4865, is a putative β-type carbonic anhydrase with similarity (44% amino acid identity) to a β-type carbonic anhydrase of Chlamydomonas reinhardtii. The β-type carbonic anhydrase of C. reinhardtii is expressed in mitochondria of cells grown under low CO2 and is believed to function in stabilization of the mitochondrial pH during photorespiration (14). Carbonic anhydrases also play a role in cyanobacteria by concentrating CO2 in carboxysomes to enhance RuBisCO activity (4). The role of carbonic anhydrases in chemoautotrophically grown B. japonicum is not known but it will be interesting to determine if they play a role in enhancing RuBisCO activity or CO2 incorporation.
Is isoctrate lyase essential for chemoautotrophic growth? Among the most highly upregulated genes found in chemoautotrophically cultured cells is a putative isocitrate lyase (blr2455), which is induced approximately 30-fold (Table 3). The putative B. japonicum isocitrate lyase has 75% amino acid identity to the characterized aceA gene of Pseudomonas aeruginosa PAO1 (10). Isocitrate lyase and malate synthase together comprise the glyoxylate cycle of plants and bacteria. The glyoxylate cycle serves as a bypass of the tricarboxylic acid (TCA) cycle and is used during growth on hexose sugars and glycolytic intermediates when an anaplerotic function is required to regenerate four-carbon compounds of the TCA cycle lost for biosynthetic purposes (29). It is also utilized by many bacteria for growth on acetate and other two-carbon substrates when avoidance of the two decarboxylation reactions of the TCA cycle is warranted to prevent a futile cycle (22). In B. japonicum, isocitrate lyase activity is readily detected when cells are grown on carbon sources that are metabolized via acetyl coenzyme A, including acetate and poly-β-hydroxybutyrate (22). Less activity is detected when grown on galactose or pyruvate, and no activity is supported by arabinose, malate, or glycerol (22).
Based upon the strong induction of isocitrate lyase, it was hypothesized that the isocitrate lyase activity of the glyoxylate cycle in B. japonicum would play an important role in partitioning carbon flow between energy generation and biosynthetic processes when cultured chemoautotrophically. To further address the role of isocitrate lyase in chemoautotrophic growth of B. japonicum, a deletion mutant in aceA (the putative isocitrate lyase, blr2455) was constructed. As indicated above, isocitrate lyase activity is most readily detected in B. japonicum when acetate is provided as a carbon source. The mutant and wild type were cultured aerobically in minimal medium with acetate provided as the sole carbon source. Consistent with isocitrate lyase mutations in other organisms, acetate supported growth of the wild-type strain but not the mutant (10, 40).
The metabolism of the aceA mutant was further characterized by using Biolog GN2 microplates. Results of this analysis indicated that the mutant failed to utilize leucine, Tween 40, Tween 80, or propionic acid, whereas growth on these substrates was observed in the wild-type strain. It was previously reported that aceA is also involved in the assimilation of acetate and leucine in P. aeruginosa (10). Similar to acetate, all of these compounds are metabolized through acetyl coenzyme A and enter metabolism via the TCA cycle or the glyoxylate cycle. These results support the identification of blr2455 as isocitrate lyase (aceA).
Chemoautotrophic growth of the isocitrate lyase mutant and wild type were similar, and no differences in growth rate or final culture density were observed (data not shown), indicating a nonessential role for this gene in chemoautotrophy. Therefore, irrespective of the strong expression of isocitrate lyase under chemoautotrophic growth and contrary to the initial hypothesis, the data indicate that a functioning TCA cycle is sufficient to meet both the energetic and biosynthetic needs of the organism in the absence of an isocitrate lyase-mediated glyoxylate bypass.
The second enzyme of the glyoxylate cycle, malate synthase, was not identified as differentially expressed under chemoautotrophic conditions. However, Green et al. (22) reported malate synthase activity in arabinose-grown B. japonicum, indicating transcripts may be present under all our growth conditions and as such the activity was not detected as differentially expressed. Another enzyme serving an anaplerotic function is phosphoenolpyruvate (PEP) carboxylase, which produces oxaloacetate via carboxylation of phosphoenolpyruvate. Transcripts for PEP carboxylase (blr2955) were upregulated 1.6- to 1.7-fold under chemoautotrophic conditions. Therefore, our data suggest that B. japonicum cultured chemoautotrophically employs the anaplerotic reactions of the glyoxylate cycle and PEP carboxylase to generate four-carbon intermediates for biosynthetic processes, while simultaneously limiting the two decarboxylation reactions of the TCA cycle.
Nitrogen metabolism. B. japonicum was previously reported to utilize nitrate and ammonia as nitrogen sources during chemoautotrophic growth (33). However, N2 as a sole nitrogen source would not support chemoautotrophic growth (33). Interestingly, when grown chemoautotrophically in the presence of ammonia, B. japonicum expresses the core components of the nitrogenase: nifH, nifD, and nifK, two ferrodoxins involved in nitrogen fixation, and fixC, and all were induced at levels 1.5- to 4.5-fold greater than under the heterotrophic condition (Table 3). B. japonicum nitrogenase expression in symbiotic bacteroids is primarily regulated by oxygen status, in contrast to free-living diazotrophs, in which oxygen and nitrogen status play a combined role (11, 17). We hypothesized that the unnecessary expression and synthesis of the nitrogenase complex under chemoautrophic conditions in the presence of adequate inorganic nitrogen would reduce fitness of the bacterium and as a result nitrogenase mutants would display enhanced chemoautotrophic growth. However, when tested for growth under chemoautotrophic conditions, the B. japonicum nifH and nifD mutants (23) grew the same as wild-type B. japonicum (data not shown), indicating there is no obvious fitness penalty associated with expression of the nitrogenase genes. Based upon the inability of B. japonicum to grow chemoautotrophically utilizing N as a sole nitrogen source, we concluded that the low-level expression of the nitrogenase components does not contribute significantly to nitrogen acquisition in chemoautotrophically grown cells but rather likely reflects a transcriptional artifact of growth under microaerobic conditions.
Chaperones. Molecular chaperones are required for proper folding and assembly of proteins and higher-order protein structures, often in response to environmental stresses. Fifteen chaperones were identified as upregulated and six downregulated in the two chemoautotrophic conditions, respectively (see Tables S1 and S2 in the supplemental material). Included are GroEL3, GroES4, GroEL4, GroES5, and GroEL5 (Table 3). Nitrogenase activity in B. japonicum bacteroids requires either GroEL3 or GroEL4 for full activity (19). GroEL3 was identified as a member of the NifA regulon and is expressed in bacteroids and anaerobic cultures (18, 25). groEL4 was previously reported to be expressed in aerobic and anaerobic cultures at similar levels and to a lesser extent in bacteroids (18). Moreover, our studies indicate that groEL4 was upregulated six- to sevenfold in chemoautotrophic cultures (Table 3). groEL5 was also reported to be expressed at low levels in anaerobic cultures (18), and we observed three- to fourfold-higher expression in chemoautotrophic cultures relative to heterotrophically grown cells. The apparent change in chaperone transcript profiles suggests that the microaerobic environment required for induction of chemoautotrophic growth also leads to expression of a different subset of chaperones that are optimized to function under low oxygen tension.
Quantitative RT-PCR. To validate and provide an indication of the sensitivity of the microarray data, we performed qRT-PCR for eight genes listed in Table 2. The genes were selected to represent the range of induction levels obtained for genes up- or downregulated in chemoautotrophically grown cells. The qRT-PCR analysis was performed on two biological replicates of chemoautotrophic and heterotrophic cells that were generated independently of those used in the microarray analysis. The comparison of qRT-PCR data to that obtained from the microarray analysis is presented in Fig. 2. The results demonstrate good agreement of expression values (R2 = 0.94) between the two methodologies with different RNA samples used for each approach, indicating that the transcriptional profiles are accurate and reproducible.
![]() View larger version (20K): [in a new window] |
FIG. 2. Comparison of log-transformed qRT-PCR data and microarray data of eight representative genes selected from Tables S1 to S4 in the supplemental material.
|
This research was supported by grant 2004-35064-14708 from the USDA/CSREES/NRI.
Published ahead of print on 8 August 2008. ![]()
Supplemental material for this article may be found at http://jb.asm.org/. ![]()
Present address: Applied Biology and Aerosol Technology, Battelle Memorial Institute, 505 King Avenue, Columbus, OH 43201. ![]()
|
|
|---|
54 regulon, and identification of a ferredoxin gene (fdxN) for symbiotic nitrogen fixation. Mol. Genet. Genomics 278:255-271.[CrossRef][Medline]
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»