Previous Article | Next Article ![]()
Applied and Environmental Microbiology, January 2008, p. 312-318, Vol. 74, No. 1
0099-2240/08/$08.00+0 doi:10.1128/AEM.01654-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Evolutionary and Genomic Microbiology Laboratory, Department of Biology,1 Department of Microbiology & Immunology, University of Louisville, Louisville, Kentucky 402922
Received 19 July 2007/ Accepted 30 October 2007
|
|
|---|
|
|
|---|
NH2OH + H2O (10). The subsequent oxidation of hydroxylamine to nitrite is facilitated by the soluble periplasmic enzyme hydroxylamine oxidoreductase (HAO): NH2OH + H2O
NO2– + 5H+ + 4e– (10). The oxidation of hydroxylamine to nitrite yields four electrons, of which two are returned to the upstream monooxygenase reaction and two are the sole source for generating useable energy and reductant. The mechanism of returning the electrons to AMO is unknown (2). AMO is encoded by at least three contiguous genes, amoCAB, arranged in a gene cluster that is conserved in all investigated genomes of AOB (3, 8, 12, 14, 15, 20-22). Prior work had identified a conserved open reading frame (ORF) following the terminator downstream of the single cluster of amoCAB genes in N. oceani (1, 20). A recent analysis of available genome sequences revealed that all amo gene clusters in betaproteobacterial AOB (beta-AOB) genomes, which contain two or three copies of nearly identical gene clusters, are actually succeeded by two conserved ORFs, orf4 and orf5, except for the amo-hao supercluster in Nitrosospira multiformis (2, 12, 22). In contrast, the terminator downstream of the amoB gene in gammaproteobacterial AOB (gamma-AOB) is succeeded only by the orf5 gene (1), which is conserved at the levels of DNA and protein sequence in all AOB. Our recent analysis of the N. oceani genome sequence revealed that the intergenic region between the amoB terminator and orf5 did not contain a promoter consensus sequence. Furthermore, examination of upstream flanking sequence of the N. oceani amoCAB cluster revealed an additional 213-bp ORF (orf1) that has no homologue in the nonredundant GenBank database including the published (8, 22) and unpublished (12) genomes of beta-AOB. The aforementioned gene structure differences between gamma- and beta-AOB indicate a possible divergent expression of amo genes in gamma- and beta-AOB, which may correspond to their respective niche adaptation in the environment (e.g., gamma-AOB are restricted to marine oligotrophic environments).
In this paper, we report that the orf1 and orf5 genes are cotranscribed with amoCAB in N. oceani, designate these ORFs as amoR and amoD, respectively, and propose that these five genes constitute the gamma-AOB-typical amo operon, amoRCABD.
In our experiments, N. oceani C-107 ATCC 19707 (18, 25) was grown in artificial seawater as 200- to 400-ml batch cultures in 2-liter Erlenmeyer flasks for 3 weeks at 30°C in the dark without shaking as described previously (1). Genomic DNA was isolated from cells in stationary growth phase using a Wizard genomic DNA (gDNA) purification kit (Promega, Madison, WI) according to the manufacturer's recommendations. For RNA preparations, N. oceani cells were harvested at mid-exponential to late exponential growth phase and resuspended for 24 h prior to RNA isolation in 200 ml of fresh marine medium. RNA was isolated using a Fast RNA Pro Blue kit (Q-Biogene, Solon, Ohio) according to the manufacturer's guidelines. Before cDNA synthesis, the RNA preparations were treated with RNase-free RQ1 DNase (Promega) according to the manufacturer's protocol.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Primers
|
![]() View larger version (32K): [in a new window] |
FIG. 3. Map of the amo gene cluster based on the genome sequence of N. oceani ATCC 19707 (A) and PCR amplification of intergenic regions in the amo gene cluster using cDNA (B) and gDNA (C) as templates. (A) The locations of primers used for amplification of intergenic sequence from cDNA are indicated above the map, and the locations of the sequence complementary to probes used for Northern analysis are indicated as horizontal bars within the arrows that indicate length and location of the genes in the amo gene cluster. Transcriptional start sites and terminators are indicated by flags and circles, respectively. Open circles indicate leaky terminators. (B and C) cDNA and gDNA, respectively, were PCR amplified with the primers indicated in panel A and listed in Table 1. Each reaction mixture also included primers for amplification of the 16S rRNA gene. Lane 1, ladder; lane 2, tspF and R-CR; lane 3, R-CF and R-CR; lane 4, C-AF and C-AR; lane 5, A-BF and A-BR; lane 6, B-DF and B-DR.
|
![]() View larger version (81K): [in a new window] |
FIG. 1. Northern hybridization analysis using probes based on the amoC, amoA, amoB, and orf5 genes. The sizes of the observed bands are based on a 9-kb RNA ladder. The blot with the probe targeting orf5 yielded only the larger (4.6-kb) band.
|
|
|
|---|
-32P]ATP (Perkin Elmer, Boston, MA; specific activity, 3,000 Ci mmol–1) using T4 polynucleotide kinase (Promega). X1 was complementary to positions 91 to 63 in the nucleotide sequence upstream of orf1 (Fig. 2A). Primer extension with total RNA was conducted using SuperScript III RT (Invitrogen) following the manufacturer's protocol. The extension products were electrophoresed in a 10% denaturing polyacrylamide 10- by 10-cm minigel, exposed to a PhosphorImager screen, and analyzed with molecular imager FX and Quantity One image analysis software (Bio-Rad, Hercules, CA). A 5'-end-labeled FX174 HinfI DNA marker (Promega) was used as a size estimate. The resulting cDNA product was 160 to 180 nucleotides (Fig. 2B) and identified the transcriptional start point between nucleotides 223 and 232 upstream of the translational start for orf1 (Fig. 2A). The –10 and –35 consensus sequences of the putative operon promoter were preceded by an A track (Fig. 2A), which is known to enhance transcription (9). In contrast to beta-AOB, where the distal amo promoter is located 166 nucleotides upstream of amoC (6), the operon promoter identified in N. oceani is much further upstream of the amoC gene, thereby generating a leader sequence of at least 600 bp. In addition to the identified orf1, this leader contained 220 to 230 nucleotides of untranslated RNA located upstream of orf1. While the function of orf1 and the significance of the 5' untranslated region were outside the scope of this study, the cotranscription of orf1 with the amoCAB genes, its size and predicted topology as a small cytoplasmic alpha-helical protein (using PSIPRED [7]), and its uniqueness to nitrosococci suggest a regulatory involvement in ammonia catabolism of gamma-AOB. For these reasons, we consider the orf1 gene a member of the amo operon in N. oceani and designate it as amoR.
![]() View larger version (40K): [in a new window] |
FIG. 2. Primer extension experiment performed to identify the TSPs of the two transcripts observed in Northern analysis (Fig. 1). (A) Annotated sequence upstream of the amoC gene, including orf1 (amoR). (B) Lane 1, single-stranded DNA ladder; lane 2, extension product obtained with primer X1 (see panel A and Table 1).
|
|
|
|---|
First-strand cDNA synthesis from N. oceani total RNA was carried out with SuperScript II (Invitrogen) and random nonamer primers according to the manufacturer's protocol. Second-strand synthesis was conducted in 100-µl reaction mixtures with the entire volume of the previous RT reaction and 15 U of Klenow fragment (Promega, Madison, WI) according to the manufacturer's protocol.
All final multiplex PCR assays were carried out using GoTaq Flexi DNA polymerase (Promega) with one of the following templates: (i) 8 to 10 ng of the cDNA, (ii) DNase-treated RNA preparations obtained before the reverse transcription for the negative controls, (iii) approximately 10 ng of gDNA as a positive control. Primers targeting the 16S ribosomal sequence were included as an internal standard in every multiplex PCR amplification reaction, including the positive and the negative controls. Negative-control PCRs did not result in any PCR products, which indicated the absence of DNA carryover in these preparations (data not shown). Positive-control PCRs resulted in bands of the expected sizes (Fig. 3C), indicating the specificity of the primers. Results of multiplex PCR assays conducted with cDNA as the template are presented in Fig. 3B. Amplification of the upstream leader sequence using tspF (Table 1; complementary to nucleotide positions 560 to 539 upstream of the amoC start codon) and the reverse primer R-CR (Table 1; complementary to nucleotide positions –12 to + 26 with respect to the amoC start codon) resulted in an approximately 585-bp band (Fig. 3B, lane 2). Amplification of the amoR-amoC intergenic region was conducted with forward primer R-CF (Table 1, complementary to nucleotide positions + 56 to + 82 within amoR), and the reverse primer R-CR yielded a PCR product of 347 bp (Fig. 3B, lane 3). These results indicate that amoC resides on a transcript with a leader of at least 600 bp including the amoR gene; therefore, amoR and amoC are transcriptionally linked. Amplification of the amoC-amoA and amoA-amoB intergenic spacers resulted in 290-bp and 450-bp bands, respectively (Fig. 3B, lanes 4 and 5). This confirms that, similarly to beta-AOB (19), the three amoCAB genes reside on a common transcript in N. oceani and that all four amplicons (Fig. 3B, lanes 2, 3, 4, and 5) originated from one or both transcripts observed in the Northern hybridization experiment (Fig. 1).
To investigate the difference between the two observed transcripts, the amoB-amoD intergenic spacer was amplified with forward primer B-DF, which targets the C terminus of amoB, and reverse primer B-DR, which is complementary to the N terminus of amoD. The observed band of
330 bp (Fig. 3B, lane 6) indicates that amoB and amoD are transcriptionally linked, in that the amoRCAB transcript extends beyond the amoB terminator, the amoB-amoD intergenic spacer, and the amoD gene and terminates at an in silico-identified rho-independent terminator downstream of amoD (
G = –35.3 kcal/mol; start at nucleotide 636,765 in the N. oceani genome sequence; GenBank no. CP000127 [15]), resulting in a 4.6-kb amoRCABD transcript as observed in the Northern hybridization experiment (Fig. 1).
|
|
|---|
None of the searched databases contained a sequence with significant similarity to the AmoR protein; hence, this protein appears to be unique to Nitrosococcus. The AmoR ORF is preceded by a Shine-Dalgarno sequence (Fig. 2A), and the protein is predicted to be cytoplasmic (19). Its sequence of 71 amino acids (Fig. 2A) contains three small helices of nine, four, and nine amino acid residues at the N terminus, whereas the C terminus contains a larger helical domain of 19 residues. Using a combination of sequence and predicted structure of the AmoR protein for a search of the PDB (http://www.rcsb.org/pdb/home/home.do), the N-terminal folding domain of frizzled-related protein 3 (PDB entry 1ijxA0) was found to be nearly identical in its fold and 23% identical to primary sequence of AmoR, extending from residue 25 to 105 (126 residues total). The deduced AmoR protein sequence contains three cysteine residues that are not part of a known coordination motif. By analogy to the structure of frizzled-related protein 3 (1ijxA0), cysteine 37 and cysteine 71 of AmoR could form a disulfide bond. Thus, AmoR could serve as a cytoplasmic redox sensor in that this disulfide bond would lock AmoR into a particular reactive secondary structure if the cytoplasm becomes less reducing. The N terminus of frizzled-related protein 3 is involved in protein-protein interactions. Based on this analogy and the fact that the amoR gene is a member of the amo operon, we speculate that the AmoR protein participates in the regulation of ammonia catabolism in Nitrosococcus. Experiments have been initiated to test this hypothesis.
A multiple sequence alignment was produced from a total of 25 (20 full-length and 5 partial) available protein sequences using the deduced full-length AmoD protein sequence from N. oceani ATCC 19707, respective BLAST hits retrieved from the GenBank/EMBL database, and unpublished genome sequences using ClustalX version 1.83 (24). Based on this alignment, a distance neighbor-joining tree was constructed with the BioNJ function in PAUP version 4.10b (23) and used as a guide tree for manual refinement of the ClustalX alignment. Sources for the protein sequences from AOB and other organisms used in the alignments are indicated in the presented phylogenetic tree (Fig. 4). The alignment was subjected to a Bayesian inference of phylogeny (MrBayes version 3.0b4; http://mrbayes.scs.fsu.edu) using four equally heated Markov chains over 1,000,000 generations in three independent runs. The searches were conducted assuming an equal or a gamma distribution of rates across sites, sampling every 100th generation and using the Whelan and Goldman empirical amino acid substitution model (26). A 50% majority rule consensus phylogram was constructed that displayed the mean branch lengths and posterior probability values of the observed clades.
![]() View larger version (37K): [in a new window] |
FIG. 4. Unrooted phylogenetic consensus trees constructed after Bayesian analysis of an alignment of available protein sequences homologous to AmoD (orf5) and AmoE (orf4). Posterior probability values smaller than 1.0 are given at the nodes; branch lengths reflect the evolutionary distance based on the standard provided for 20 changes over time. (A) Labels indicate the sequence source. amoD (pmoD) and amoE genes are clustered with the structural amoCAB and pmoCAB genes, whereas shaded boxes identify the singleton status of amoD and amoE in the source genome. , a homologue of the multicopper oxidase, MCA2129, is located upstream of the amo gene cluster in gamma-AOB; ¶, orf5 (amoD) is located between amoE and copC in beta-AOB. Structural analysis identified AmoD and AmoE as periplasmic membrane proteins. (B) Unrooted star tree presenting the relationships between major clades. BlastP top hits are provided in support of our hypothesis that amoE and singleton genes were derived from amoD by duplication (see the text). Abbreviations: Ns., Nitrosospira; Nm., Nitrosomonas; Nc., Nitrosococcus; Mc., Methylococcus (indicated with an asterisk). **, alpha-MOB.
|
Very recent work demonstrated that expression of the betaproteobacterial orf4 and orf5 gene tandems is significantly up-regulated during recovery from ammonia starvation (5). Interestingly, all orf4-orf5 gene tandems are flanked downstream by copper resistance genes (copCD) in the genome sequences of beta-AOB (2). Furthermore, our database search revealed that homologues of amoD (but not orf4) genes also reside in the vicinity of genes encoding copper enzymes in some genomes of MOB. In alpha-MOB, amoD homologues reside downstream of pmo genes (CAJ01564 to CAJ01562 and AAF37892 to AAF37894), which encode a homologue of AMO (1, 20). In gamma-MOB, an amoD homologue was found adjacent to genes encoding multicopper oxidases (MCA2129 and MCA2128) that are homologues of the copper oxidase gene copA upstream of the amoRCABD operon in N. oceani. Additionally, an amoD homologue, which had been missed in the annotation of the M. capsulatus genome, was identified downstream of a gene that encodes the copper resistance protein CopC (MCA2170), and both are likely coexpressed. Genes encoding proteins in the AmoD-AmoE family were found only in genomes that encode AMO or pMMO (Fig. 4).
In silico analysis of the deduced AmoD protein structure revealed that AmoD is likely exported to the periplasm but stays anchored with its C terminus in the inner membrane (19). Likewise, the AmoD protein is likely transported to and integrated into the extensive intracytoplasmic membrane system of AOB, where it may interact with electron transfer proteins and enzymes that facilitate the oxidation of ammonia or their maturation. This hypothesis is based on our findings that amoD is an expressed member of the amo operon in Nitrosococcus, that it is also expressed in beta-AOB (5), and that it is conserved in sequence and synteny in the amo gene clusters of all AOB, all of which suggests that the amo cluster genes encode proteins that interact physically (4, 11, 17).
In summary, we hypothesize that the amoD gene is distributed in AOB and MOB similar to other inventory that is involved in nitrification (i.e., amo [pmo] and hao [2, 4, 16]) and that AmoD is unique to organisms capable of ammonia oxidation and nitrification. Because the amoD gene is ancestral to orf4 (Fig. 4), which has likely arisen by complete gene duplication from amoD in the ancestor of all beta-AOB, we propose naming the orf4 gene of beta-AOB "amoE." By analogy to amoC and pmoC, we propose naming all amoD-homologous singleton genes identified in AOB and MOB "amoD" and "pmoD," respectively.
|
|
|---|
This project was supported, in part, by incentive funds provided by the University of Louisville-EVPR office, the KY Science and Engineering Foundation (KSEF-787-RDE-007), and the National Science Foundation (EF-0412129).
Published ahead of print on 9 November 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»