Journal of Bacteriology, June 2007, p. 3935-3944, Vol. 189, No. 11
0021-9193/07/$08.00+0 doi:10.1128/JB.01861-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

Department of Microbiology, University of Washington,1 Department of Civil and Environmental Engineering, University of Washington, Seattle, Washington2
Received 11 December 2006/ Accepted 14 March 2007
|
|
|---|
70 type amoCAB promoter was constitutively active in the presence of ammonia, but the proximal promoter was only active when cells were recovering from ammonia starvation. The lone, divergent copy of amoC (amoC3) was expressed only during recovery. Both the proximal amoC1,2 promoter and the amoC3 promoter are similar to gram-negative
E promoters, thus implicating
E in the regulation of the recovery response. Although modeling of subunit interactions suggested that a nonconservative proline substitution in AmoC3 may modify the activity of the holoenzyme, characterization of a
amoC3 strain showed no significant difference in starvation recovery under conditions evaluated. In contrast to the amo transcripts, a delayed appearance of transcripts for a gene required for CO2 fixation (cbbL) suggested that its transcription is retarded until sufficient energy is available. Overall, these data revealed a programmed exit from starvation likely involving regulation by
E and the coordinated regulation of catabolic and anabolic genes. |
|
|---|
Ammonia-oxidizing bacteria convert ammonia to nitrite (NO2) in two enzymatic reactions. Ammonia is first oxidized to hydroxylamine (NH2OH) in an energy-consuming step by the membrane-bound ammonia monooxygenase (AMO), a three-subunit holoenzyme encoded by the amoCAB operon. Electrons for energy generation and biosynthesis are derived from the oxidation of hydroxylamine to nitrite by the periplasmic hydroxylamine oxidoreductase (52). All betaproteobacterial ammonia oxidizers characterized to date have multiple copies of amoCAB and an additional copy of amoC (4, 25, 33, 37, 43). In N. europaea, the sequence of the protein encoded by the monocistronic amoC (AmoC3) diverges from the two nearly identical copies of AmoC (AmoC1,2) encoded by duplicate amoCAB operons (67.5% identity, 81.4% similarity). Although cotranscribed with amoA and amoB (43), an abundant monocistronic transcript of amoC has been shown to be remarkably stable in N. europaea during starvation (43).
Very little is known about the function of AmoC, other than a likely localization to the cell membrane as a subunit of ammonia monooxygenase (25), nor is there direct information about the possible role for the divergent copy, even though the conservation of amoC3 among ammonia-oxidizing bacteria suggests that it serves an important function. Our current analyses demonstrate that transcription of amoC3 is specific to recovery from ammonia starvation, thus implicating this subunit of the holoenzyme in the recovery of ammonia-oxidizing bacteria from starvation. Comparative promoter analysis showed that amoC3 transcription is correlated with elevated transcription of the operon copies of amoC as part of the poststarvation recovery response. Additional transcript analyses revealed a programmed exit from starvation in which transcription of genes encoding energy demanding anabolic functions are delayed relative to the synthesis/repair of central energy-generating systems.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Bacterial strains and plasmids
|
|
View this table: [in a new window] |
TABLE 2. Oligonucleotide primers and probes
|
Primer extensions.
Plasmids are described in Table 1, and oligonucleotide primers are described in Table 2. Primer extensions were performed as previously described (39). Optikinase (USB) was used to end-label primers C23, NEUamoC3_54r, NEUcbbL40r, NEUglyA32r, and NEUrpsJ95r with [
-32P]dATP (3,000 Ci/mmol; PerkinElmer). For primer extensions using RNA from the 500-h starvation of N. europaea ATCC 19718, 2.5 µg total RNA was used with primer C23 and 4.5 µg total RNA was used with primers NEUamoC3_54r and NEUcbbL40r. For primer extensions using RNA from the 144-h starvation of N. europaea NCIMB 11850, 3 µg of RNA was used in reactions with all of the above primers. Superscript III reverse transcriptase (Invitrogen) was used with N. europaea ATCC 19718 RNA, and Superscript II reverse transcriptase (Invitrogen) was used with NCIMB 11850 RNA. Sequencing ladders were prepared by cycle sequencing using the SequiTherm EXCEL II DNA sequencing kit (Epicenter). Plasmids pPMB8, pPMB10, and pPMB11 were used as templates with the respective labeled primers C23, NEUamoC3_54r, and NEUcbbL40r. Primer extension reactions were fractionated on a 6% polyacrylamide gel and visualized by autoradiography of dried gels.
S1 nuclease digestion.
S1 analysis was preformed as previously described (39) using RNA from N. europaea ATCC 19718 starved for 2 weeks in spent medium before the addition of 12.5 mM (NH4)2SO4. Each 80-mer probe (amoC1-140r and amoC1-66r) was purified by polyacrylamide gel electrophoresis. Optikinase (USB) was used to end-label amoC1-140r (targeting amoCp1) and amoC1-66r (targeting amoCp2) with [
-32P]dATP (3,000 Ci/mmol; PerkinElmer). Each probe was hybridized with 5 µg of total RNA before digestion with 50, 250, or 500 U/ml S1 nuclease (Invitrogen). The digestion products were fractionated on a 6% polyacrylamide gel and visualized as described above.
mRNA secondary structure prediction.
The web interface for the mfold version 3.0 algorithm (http://www.bioinfo.rpi.edu/applications/mfold/) was used for all RNA structure predictions (49). As this version allows predictions at different temperatures, initial foldings were run using temperatures between 65°C and 85°C to predict the most stable secondary structures for each RNA molecule. To obtain
G values for each structure, mfold version 3.1 was employed at a fixed temperature of 37°C (32, 54).
Construction of a
amoC3 strain.
The allelic exchange vector, pCM351, was used to replace the amoC3 gene with the gentamicin resistance determinant, aacC1, in the same orientation. The upstream flanking region of amoC3 was amplified by PCR using primers C3upEcoR1F1 and C3upNco1R1. The resulting PCR product was digested with EcoRI and NcoI and cloned upstream of aacC1 in pCM351 to create pPMB25. The downstream flanking region of amoC3 was amplified using primers C3downApa1F2 and C3downSac1R2. The resulting PCR product was digested with ApaI and SacI and cloned downstream of aacC1 to create pPMB34. pPMB34 was linearized using NheI prior to subsequent electroporation of N. europaea ATCC 19718, as previously described (15). Recombinant clones were selected by plating on solid media containing 5 µg/ml gentamicin (USB). Replacement of amoC3 was confirmed by PCR using primers specific to amoC3, the gentamicin cassette, and upstream and downstream of the recombination sites. In addition, primer extension analysis using RNA from recovering cells was used to confirm that the mutant strain did not express amoC3. The mutant was designated N. europaea PMB1.
Assay of the recovery phenotype for the
amoC3 strain.
N. europaea ATCC 19718 and N. europaea PMB1 were grown in triplicate 500-ml batch cultures. At late-exponential phase, cells were harvested by filtration and resuspended in ammonia-free medium. Cell densities were normalized using optical density at 600 nm and starved for up to 4.5 months. Periodically during the starvation, an aliquot of cells was spiked with 0.020 mM or 25 mM NH3. In a separate pulse-feed experiment, 7-day stationary-phase cultures of N. europaea ATCC 19718 and N. europaea PMB1 were transferred in triplicate to ammonia-free medium. NH3 was periodically added to the cultures to induce several cycles of starvation and recovery. The recovery response in all experiments was monitored by measuring nitrite production as previously described (23).
AMO structure modeling. Structural modeling of AMO was performed with the comparative modeling protocol on the PROTINFO structure prediction server (http://protinfo.compbio.washington.edu) using the particulate methane monooxygenase crystal structure as a template (28). This comparative modeling protocol has been shown to work well in the CASP protein structure prediction experiments (18, 19). Initial models were constructed using a minimum perturbation approach that aims to preserve as much information as possible from the template structure solved by X-ray crystallography. Variable side chains and main chains were then built using a graph theory clique-finding approach that explores a variety of possible conformations for the respective side chains and main chains and finds the optimal combination using an all-atom scoring function (42). These approaches are described in further detail in references 40 and 41. PDB files of all AMO structure models are available from the Computational Biology Research Group in the Department of Microbiology at the University of Washington (http://data.compbio.washington.edu/misc/downloads/amoc/).
|
|
|---|
![]() View larger version (71K): [in a new window] |
FIG. 1. Primer extension analysis of amoC1,2, amoC3, and cbbL expression during exponential growth and recovery of N. europaea ATCC 19718 from 1, 6, and 21 days of starvation. (A) Nitrite accumulation and ratios of amoCp2 transcripts to amoCp1 transcripts (determined by densitometry) during recovery. (B) amoC1,2 expression levels, as determined by densitometry, showing the contribution of each promoter to total amoC1,2 abundance. (C, D, E) Primer extension analysis of amoC1,2, amoC3, and cbbL, respectively. Transcripts derived from each amoC1,2 promoter are indicated at left.
|
![]() View larger version (109K): [in a new window] |
FIG. 2. (A) Primer extension analysis of amoC1,2 following a 3-h recovery of a 1-month stationary-phase culture of N. europaea NCIMB 11850. Transcription start sites associated with each promoter are indicated at left. (B) amoC1,2 expression levels, as determined by densitometry.
|
Transcription of cbbL, glyA, and rpsJ. As a control, the expression of several additional genes of various functional role categories was also examined. The expression of cbbL (large subunit of ribulose 1,5-bisphosphate carboxylase/oxygenase), glyA (glycine hydroxymethyltransferase), and rpsJ (S10 ribosomal protein) were monitored in actively growing, starved, and recovering cells of N. europaea. Since cbbL encodes a subunit of an anabolic enzyme involved in CO2 assimilation (50), it was hypothesized that its transcription would be regulated differently than the genes encoding the catabolic ammonia monooxygenase. The cbbL transcript was detected in exponential-phase cells but not in starved cells (Fig. 1). Expression of cbbL during the recovery period was delayed several hours with respect to the expression of amoC1,2 and amoC3 (Fig. 1). Longer periods of starvation resulted in a longer lag between the initiation of recovery and the expression of cbbL (Fig. 1). The "housekeeping" genes, glyA and rpsJ, exhibited similar expression profiles in N. europaea NCIMB 11850. Both transcripts were detected in actively growing cells but not during the starvation period. Expression of glyA and rpsJ was induced within the first 3 h of recovery following 6 days of starvation at levels similar to that of exponential phase cells (data not shown).
Identification of transcription start sites and promoter consensus sequences.
The transcription start sites for amoC1,2, amoC3, cbbL, glyA, and rpsJ were first determined by primer extension. S1 nuclease protection was then used to confirm the transcription start sites for amoC1,2. The 35 and 10 promoter sequences were identified by visual inspection of sequences upstream of the identified transcription start sites (53). The start sites for the promoters amoCp1 and amoCp2 are, respectively, 168 to 169 and 105 to 107 nucleotides upstream of the start codon for amoC1,2 (Fig. 3). These positions differ by approximately 2 nucleotides (nt) from those reported by Hommes et al. for amoC1,2 (16), likely reflecting the use of a heterologous sequence for positional identification in the prior study. As shown previously, the amoCp1 promoter most closely matches the eubacterial
70 consensus sequence (16). Although the amoCp2 promoter also shows homology to the eubacterial
70 consensus sequence, as reported by Hommes et al. (16), further inspection of this region revealed features highly similar to the consensus sequence for the eubacterial extracytoplasmic function (ECF)
factor, encoded by the rpoE gene in N. europaea (Fig. 3). This promoter assignment also places the putative promoter closer to the 35 and 10 regions. The putative amoC3 promoter for the major primer extension product, 227 nt upstream of the start codon, shows some homology to
70 and
32 consensus sequences, but also most closely matches the consensus sequence for ECF
factors (Fig. 3). No obvious promoter sequence could be found upstream of the putative start site corresponding to the minor primer extension product for amoC3 (209 nt upstream of the amoC3 start codon). The transcription start sites for cbbL, glyA, and rpsJ are 83 to 84, 59 to 60, and 23 nt upstream of their respective start codons. These genes, all primarily up-regulated during exponential growth, appear to be controlled by
70 promoters (Fig. 3).
![]() View larger version (31K): [in a new window] |
FIG. 3. Transcription start sites, promoters, and predicted 5' untranslated region secondary structures associated with amoC1,2, amoC3, glyA, rpsJ, and cbbL. Transcription start sites determined by primer extension and S1 analysis are indicated by open circles and diamonds, respectively. Putative 70 promoters are underlined and putative ECF promoters are boxed. G values for each structure were determined by mfold version 3.1 at a fixed temperature of 37°C. Start codons (AUG) and stop codons (UAA) are indicated as needed. The positions of the transcription start sites relative to the start codons are indicated next to each start site. Bacterial promoter consensus sequences were obtained from the review by Wösten (53).
|
G values of 12.9 and 11.3 kcal/mol, respectively (Fig. 3). In addition, a highly stable (
G = 19.1 kcal/mol) stem-loop is predicted immediately downstream of the stop codon for amoC1,2 (Fig. 3). The 5' termini for both primer extension products for amoC3 are located at the base of a hairpin structure with a
G value of 7.3 kcal/mol (Fig. 3). In contrast, the cbbL transcript lacks a stable stem-loop structure associated with its transcription start site but does have a predicted structure located approximately 32 nt downstream of the start site with a
G value of 10.9 kcal/mol (Fig. 3).
Characterization of a
amoC3 mutant.
An amoC3 deletion mutant of N. europaea ATCC 19718 (
amoC3::aacC1 strain PMB1) was constructed to further evaluate possible differences in the function of AmoC1,2 and AmoC3. The deletion was confirmed by PCR and the absence of an amoC3 transcript during recovery from starvation (Fig. 4). The deletion mutant was initially characterized by examining its recovery from starvation relative to the wild type. No significant differences in nitrite production between the wild type and mutant were observed over a 10-day recovery period following the addition of 25 mM NH3 to cultures starved of ammonia for 4.5 months (data not shown). The recovery phenotype was also assayed after addition of 20 µM NH3 to a culture starved for 47 days, a concentration below the estimated KS of 50 µM NH3 for N. europaea (26) (data not shown). As with the higher concentration of NH3, no significant differences in nitrite production were observed. Pulse-feeding of ammonia to 7-day-old stationary-phase cultures of wild-type and mutant strains was then examined, since our transcription data suggest that this feeding regimen would selectively enrich for AmoC3 (Fig. 3). A slight but replicable difference (standard deviation) in ammonia oxidation activity was observed when the cultures were initially pulse-fed 400 µM ammonia. However, nitrite production converged with subsequent additions of ammonia (data not shown).
![]() View larger version (54K): [in a new window] |
FIG. 4. Verification of the amoC3 deletion. (A) PCR analysis of the recombination site. (B) Primers used (Table 1) and PCR products corresponding to each lane. (C) Primer extension analysis of amoC3 expression in wild-type (WT) and amoC3 cells following 20 min of recovery in 25 mM NH3 after 24 h of starvation.
|
A comparative alignment of PmoC and AmoC sequences revealed an amino acid substitution in AmoC3 near a conserved metal binding center. As revealed by the crystal structure of PmoC, this tetrahedral metal binding center involves 3 residues (D156, H160, and H173) from PmoC and 1 residue (E195) from PmoA. All 3 residues from PmoC are conserved in both AmoC1,2 and AmoC3 (Fig. 5). However, there is a nonconservative proline substitution at position 157 in AmoC3 (V155 in AmoC1 and I146 in PmoC1) (Fig. 5). PmoC3 does not have a substitution at this position compared to the operon copies of PmoC (Fig. 5). The structural consequences of this mutation were then examined using a knowledge-based all-atom scoring function developed for protein structure prediction. The all-atom function computes a stability score by summing up the individual scores of all the atomic interactions in a protein structure. The individual scores are derived from atomic preferences in a database of experimentally determined structures (42).
|
View larger version (33K): [in a new window] |
FIG. 5. ClustalW alignment of AmoC and PmoC amino acid sequences. Identities are highlighted in black, and similarities are highlighted in gray. Conserved residues involved in coordination of the tetrahedral metal binding site are indicated by asterisks. The proline substitutions in the AmoC3 proteins from N. europaea ATCC 19718 and Nitrosomonas eutropha C71 are boxed.
|
|
|
|---|
Our studies demonstrate a major role of the AmoC subunit in recovery from starvation. The long-lived amoC1,2 message presumably contributes to recovery following relatively short-term starvation periods of several days. Recovery from longer-term starvation is associated with a high level of transcription from both promoters for the amoC1,2 operons and a recovery-specific activation of transcription of amoC3. While these analyses confirmed the presence of two promoters for amoC1,2, as previously reported by Hommes et al. (16), the earlier study showed, in contrast to our results, that transcripts originating from both promoters are present in starved cells and that only the distal promoter is active during outgrowth from starvation (16). Differences in starvation and recovery conditions may account for some differences in results between these two studies. In the Hommes et al. study, N. europaea ATCC 19718 was grown for 3 days to late exponential phase, harvested by centrifugation, washed three times with buffer (0.2 mM MgCl2 and 50 mM NaH2PO4, pH 7.5), and stored as a pellet at 4°C for at least 24 h to allow the endogenous mRNA to decay before transcription was induced by the addition of 50 mM NH3 (16). In contrast, our protocol limited the manipulation of the cells to more physiologically relevant conditions; cells were starved as a suspension in fresh growth medium lacking ammonia at room temperature (approximately 21°C). To address these differences, the effects of starvation time and temperature on transcription of amoC1,2 were investigated. When we used our protocol to examine cells starved for either 1 or 6 days at 4°C or 20°C in ammonia-free medium, the addition of ammonia induced transcription from both promoters (data not shown). These studies also showed that mRNA was more stable when cells were stored at 4°C; the longer transcript could still be detected after 6 days of starvation at 4°C but not after 6 days of starvation at 20°C (data not shown). Thus, variation in conditions antecedent to harvest may have influenced the results earlier reported (16).
Control of amoC1,2 and amoC3 transcription by an ECF sigma factor implies a role of AmoC in stress response. ECF sigma factors are often responsive to environmental stimuli and membrane stress conditions that result in misfolded proteins in the periplasm (1, 53). The fact that AMO is a membrane protein is consistent with this function of ECF sigma factors. Inactivation of the membrane-bound AMO during starvation may serve as one signal for increased transcription of the amoCAB operons and recovery-specific transcription of amoC3 during the subsequent recovery period.
As initiation is typically the rate-limiting step of transcription (17), our data indicate that the tandem promoters upstream of the amoCAB operon are utilized to enhance transcription during the recovery period. There are a number of examples of tandem promoters in bacteria, many of which are regulated by different sigma factors (5, 10, 12, 13, 36). For instance, in Shewanella violacea, the tandem
70 and
54 promoters that control expression of glnA are differentially regulated with respect to hydrostatic pressure (20). Also, the tandem rrn promoters in E. coli are differentially regulated with respect to growth rate such that the proximal promoter is preferentially utilized during the outgrowth of cells from stationary phase (35). Utilization of amoCp2 presumably enables cells to increase expression levels beyond that observed in actively growing cells (Fig. 1).
Transcript stability can determine protein stoichiometry in some systems and has also been suggested to conserve energy at the expense of tight regulatory control (14, 27). The stability of the bulk mRNA pool increases during nutrient limitation in many bacteria and is hypothesized to enhance survival by decreasing the energetic requirements for transcription (2, 3, 24). The proximity of a stable secondary structure to the 5' end of a transcript (within 4 to 7 bp) has been shown to be the primary determinant of messenger stability (7, 9, 11, 45). The 5' stem-loop structures associated with amoC1,2 transcripts are within 4 bp of the termini and should be sufficient to impart significant stability to the message (Fig. 3). In addition, there is a stable structure located downstream of the amoC1,2 stop codon that likely protects these transcripts from 3'-5' degradation (9). Although this structure lacks a downstream poly(U) sequence characteristic of many intrinsic transcription terminators in E. coli, there is evidence that a poly(U) sequence is not necessary for efficient termination in some species (46). While further studies are warranted, these predicted secondary structures could explain the high abundance of the 1.1-kb monocistronic amoC message present in N. europaea cells (16, 43).
In contrast to the elevated expression of genes encoding enzymes responsible for catabolic processes and energy generation, strict control over CO2 fixation is supported by a delay in the expression of cbbL encoding the large subunit of RuBisCO (Fig. 1). These data are consistent with a previous study by Johnstone and Jones which demonstrated that CO2 uptake was delayed approximately 4 h before it began to steadily increase in recovering cells of N. cryotolerans that had been starved for 5 weeks (22). Unlike glyA and rpsJ which also have putative
70 promoters, cbbL is not expressed immediately upon addition of ammonia to starved cells. This result indicates that another regulatory factor is involved in the control of cbbL expression during recovery. Since there is no apparent stable stem-loop structure associated with the 5' terminus of the cbbL transcript, a lower stability compared to amoC may impart greater transcriptional control of cbbL messenger abundance.
The third, divergent, lone copy of amoC (amoC3) apparently has a specific role in the recovery of cells from starvation. Similar to the operon copies of amoC, the transcription start site of amoC3 is associated with a short stem-loop structure that may influence the stability of the amoC3 message (Fig. 3). Since no promoter could be assigned to the minor primer extension product, it is most likely an artifact caused by premature termination of the reverse transcriptase at this stem-loop. A poly(U) sequence immediately downstream of this stem-loop structure implicates another regulatory control, since probable transcription termination in this region would need to be suppressed during recovery. However, the amoC3 deletion mutant showed only modest differences from the wild type in its recovery phenotype under the conditions so far evaluated, suggesting a subtle role in recovery physiology that may in part be masked by the operon variants.
The conservation of the metal binding site in PmoC and AmoC implies a function beyond that of a chaperone or membrane anchor as suggested previously (25). Lieberman and Rosenzweig have hypothesized that the metal binding site of PmoC may function in catalysis or as an electron transfer center due to the proximity of a putative quinone binding site (29). The close proximity of the proline substitution to the metal center in AmoC3 could affect these possible functions by altering the local structure of this metal binding site. In contrast, our observations of the increased stability of the AMO holoenzyme containing AmoC3 and that an ECF sigma factor appears to control the increased transcription of all three amoC copies during recovery indicate that AmoC3 may serve a role as a specialized chaperone as suggested previously. As such, AmoC3 may facilitate association of AMO subunits in the membrane during recovery but not function optimally, as illustrated by the lack of amoC3 transcription during active growth (Fig. 1) and the likely disruptive effects of a proline substitution near the metal binding site of AmoC3 (Fig. 5). We are now examining the stoichiometry of the AMO subunits under different starvation and recovery scenarios to identify conditions that may enhance phenotypic differences between the wild type and the amoC3 deletion mutant. Many proteins involved in ammonia oxidation are stable and highly abundant (38, 51). It may not be possible to detect subtle differences between the wild type and the amoC3 deletion mutant if the level of AmoC1,2 is significantly greater that the quantities of AmoC3 produced during the recovery of wild-type cells.
Our results demonstrate that cellular resources are primarily directed toward the regeneration of ammonia oxidation activity during the recovery of ammonia-oxidizing bacteria from starvation. In particular, the high abundance and putative stability of amoC transcripts suggest that the physiological importance of AmoC has been underestimated, during both active growth and recovery from starvation. The presence of a conserved, divergent copy of amoC that is primarily expressed during recovery underscores this point and presents a unique opportunity to make further advances in understanding the structure and function of AMO, which has been recalcitrant to purification and in vitro enzymatic studies. Although technically challenging, the development of strains in which the stoichiometry of amoC copies is controlled would facilitate in vivo studies aimed at understanding the functional nature of AMO. In addition, our results indicate that the RpoE stress response regulon has a role in the recovery physiology of N. europaea. Characterizing other members of this regulon is likely to further our understanding of the unique starvation and recovery response of ammonia-oxidizing bacteria.
Published ahead of print on 23 March 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»