Previous Article | Next Article ![]()
Applied and Environmental Microbiology, September 2006, p. 6299-6315, Vol. 72, No. 9
0099-2240/06/$08.00+0 doi:10.1128/AEM.00463-06
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
University of Louisville, Louisville, Kentucky 40292,1 Oregon State University, Corvallis, Oregon 97331,2 Lawrence Livermore National Laboratory, Livermore, California 94550,3 Joint Genome Institute, Walnut Creek, California 94598,4 Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831,5 Utah State University, Logan, Utah 84322,6 Princeton University, Princeton, New Jersey 085447
Received 26 February 2006/ Accepted 16 June 2006
|
|
|---|
|
|
|---|
![]() View larger version (159K): [in a new window] |
FIG. 1. Transmission electron micrograph of Nitrosococcus oceani ATCC 19707. Cultures of N. oceani were grown until slowing readjustment of the pH with potassium carbonate indicated the beginning of transition of the culture into stationary growth phase. Cells were harvested by centrifugation and sent to the University of WisconsinMadison Electron Microscopy facility for further processing and electron microscopy. The scale bar at the lower right (500 nm) indicates the average cell size of 1.5 µm in diameter.
|
The susceptibility of ammonia oxidizers to inhibition by sunlight (due to the light sensitivity of ammonia monooxygenase [AMO]) is probably responsible for the characteristic distribution of nitrification in the water column; maximal rates occur in surface waters near the bottom of the euphotic zone (74). Nitrification rates decrease with increasing depth as the rate of organic matter decomposition (and thus ammonium supply) decreases with depth. As a consequence of both nitrifier and denitrifier activities, the oceans emit large amounts of the greenhouse gas nitrous oxide (57). The total oceanic N2O inventory is about two-thirds the size of the total atmospheric inventory, and the oceanic N2O flux to the atmosphere is estimated to be 4 Tg N/year (57). It has been shown that nitrous oxide and nitric oxide can be produced by aerobic AOB either through the reduction of nitrite (NO2 [17, 66]) or the oxidation of hydroxylamine (NH2OH [35]). Goreau et al. (25) reported on N2O production by a marine Nitrosomonas isolate and a culture of N. oceani that had been isolated from the western Atlantic by Stan Watson. N2O production was found to be much higher under subatmospheric O2 conditions than under atmospheric O2 conditions, indicating that nitrifiers may produce significant amounts of N2O in the interface between the oxic and anoxic zones. Recent geochemical evidence indicates that most of the N2O in the ocean is derived from nitrification (56); therefore, the detailed understanding of the marine nitrification process is crucial for any global management strategy of greenhouse gas production. This significance of nitrification in the oceans and its indirect influence on the oceanic carbon budget led to the selection of N. oceani ATTC 19707 as a target for genome sequencing in the Genomes to Life microbial sequencing program of the U.S. Department of Energy.
|
|
|---|
Library construction, sequencing, and sequence assembly.
Genomic DNA was isolated from late-exponential-phase cultures of N. oceani strain ATCC 19707 as described by McTavish et al. (49), modified following the recommendations by the Department of Energy's Joint Genome Institute (DOE-JGI; Walnut Creek, CA). The genome was sequenced using the whole-genome shotgun method as previously described (16, 23). Briefly, random 3- and 8-kb DNA fragments were isolated and cloned into pUC18 and pMCL200 vectors, respectively, for amplification in Escherichia coli. A larger fosmid library was constructed containing approximately 40-kb inserts of sheared genomic DNA cloned into the pCC1Fos cloning vector. Double-ended plasmid sequencing reactions were performed by the DOE-JGI using ABI 3730xl DNA Analyzers and MegaBACE 4500 Genetic Analyzers as described on the JGI website (http://www.jgi.doe.gov/).
After quality control of the 60,402 total initial reads of draft sequence, 51,334 were used for the final assembly, producing an average of 9.3-fold coverage across the genome. Processing of sequence traces, base calling, and assessment of data quality were performed with PHRED and PHRAP. Assembled sequences were visualized with CONSED. The initial assembly consisted of 36 contigs. Gaps in the sequence were primarily closed by primer walking on gap-spanning library clones or with PCR products from genomic DNA. True physical gaps were closed by combinatorial (multiplex) PCR. Sequence finishing and polishing added 249 reads, and assessment of final assembly quality was completed as described previously (16).
Sequence analysis and annotation.
Automated gene modeling was completed by combining results from Critica, Generation, and Glimmer modeling packages and comparing the translations to GenBank's nonredundant (NR) database using the basic local alignment search tool for proteins (BLASTP). The protein set was also searched against KEGG genes, InterPro, TIGRFams, PROSITE, and Clusters of Orthologous Groups of proteins (COGs) databases to further assess function. Manual functional assignments were assessed on an individual gene-by-gene basis as needed.
Nucleotide sequence accession numbers.
The sequence and annotation of the complete N. oceani strain ATCC 19707 genome is available at GenBank/EMBL/DDBJ accession numbers CP000127 (chromosome) and NC_007483 (plasmid).
|
|
|---|
![]() View larger version (23K): [in a new window] |
FIG. 2. (A) Chromosome and (B) circular plasmid of Nitrosococcus oceani ATCC 19707. The outer two circles depict predicted
protein-encoding and structural RNA genes on the plus and minus strand, respectively (green, energy metabolism; red, DNA replication; magenta, transcription; yellow, translation; orange, amino acid metabolism; dark blue, carbohydrate metabolism; pale red, nucleotide metabolism; black, coenzyme metabolism; cyan, lipid metabolism; light blue, cellular processes; brown, general function; gray, hypothetical and conserved hypothetical genes; pale green, structural RNAs). Circles 3 and 4 (in panel A only) indicate on the plus and minus strand, respectively, the
locations of the two rrn operons (black), loci involved in ammonia and urea catabolism (blue), loci involved in electron transfer and the generation of reducing equivalents (red), and predicted terminal electron acceptors (green). The two inner circles indicate GC
bias and GC skew. Cyt, cytochrome.
|
|
View this table: [in a new window] |
TABLE 1. General features of N. oceani ATCC 19707a
|
The plasmid comprises mostly hypothetical and conserved hypothetical proteins. A transposase (Noc_A0021), phage integrase (Noc_A0015), and a small number of other phage-related genes along with restriction modification systems are also found on this replicon. A possible replication protein (Noc_A0039) whose only putative homologue (56% similar over 80% of its length) is the RepA protein encoded by plasmid pRA2 from Pseudomonas alcaligenes (44), together with a possible partitioning system (Noc_A0013 to Noc_A0014), may help this plasmid be maintained within the N. oceani population.
Families and clusters of foreign and repeat sequences.
Surprisingly, we identified several large blocks of genes that were identified as putatively phage related, indicating that N. oceani has been a frequent target for bacteriophages in the open ocean (Table 2). Ten regions, ranging in size from 6.4 to 44.3 kb (totaling >175 kb), were found to harbor bacteriophage remains, such as phage integrase, terminase, primase, and tail genes. These regions were found to be associated with tRNAs (flanking), restriction modification systems, transposases, virulence, and many hypothetical genes. Other phage hits are not associated with any cluster of genes but are distributed throughout the genome.
|
View this table: [in a new window] |
TABLE 2. Mobile elements in the genome of N. oceani ATCC 19707
|
83% identity and represent an IS superfamily. Interestingly, the genome encodes many predicted transposases (127), some of which cluster in large groups, the largest of which comprises 15 consecutive ORFs (Noc_0003 to Noc_0017) encompassing >7 kb. Also, of the 14 predicted phage integrase genes or pseudogenes, two were nearly identical frameshifted integrase genes (Noc_0080 and Noc_1095).
Two identical copies of the Tu translation elongation factor (Noc_2326 and Noc_2338) and two gene copies encoding a fatty acid desaturase and a metalloprotein were identified. Additionally, there are two copies of a gene annotated as ammonia permease (Noc_2700 and Noc_2701); however, both were much shorter than AmtB proteins from other organisms and aligned only with their C termini. Furthermore, the region containing the membrane-spanning domains lacked amino acid residues that are conserved in other AmtB proteins. Hence, operation of AmtB-facilitated ammonia uptake by one or both products of the Noc_2700 and Noc_2701 genes needs to be experimentally verified. In addition to these examples, we have found several lines of evidence suggesting that genes or gene families are undergoing duplication and diversification. Examples include Noc_1310, which appears to be a truncated C-terminal version of Noc_1552, a full-length phosphoenolpyruvate-protein phosphotransferase, and the gene Noc_0725, which is only
70% identical to one of three nearly identical copies (Noc_0343, Noc_0724, and Noc_0973) of a predicted hypothetical protein of 492 amino acids.
Information processing and modification systems.
The genome of N. oceani contains the complete sets of genes necessary to encode DNA-directed DNA polymerases I (Noc_0554) and III (Noc_0002, Noc_0288, Noc_0846, Noc_1659, Noc_2593, Noc_2663, and Noc_2814) as well as multiple copies of the gene encoding the epsilon subunit of polymerase III. Genes encoding polymerase II were not identified. Nitrosococcus oceani is equipped with a full complement of genes to carry out repair of DNA lesions (uvrABCD, recFRO, mutSHLT, radC, and recN) and recombination (recA, recD, and recG; ruvABC), which may be necessary as a consequence of exposure to mutagens and uptake of foreign DNA. The genome of N. oceani contains a large number of open reading frames (a total of 24) that code for type I and type III site-specific restriction endonucleases, as concluded from a comparison with available complete genome sequences of ammonia-oxidizing Betaproteobacteria and closely related Gammaproteobacteria, whereof the next highest number of ORFs was found in the genome of Xylella fastidiosa Ann-1 (a total of 13). The presumed function of these restriction modification systems is protection against phage infection and foreign DNA by recognizing specific methylation patterns and distinguishing between host and foreign DNA (11, 15). The reason for this high number of restriction endonucleases in the genome of N. oceani is presently unclear; however, a similarly high number of ORFs encoding site-specific endonucleases have been identified in the unfinished genome of Nitrosococcus halophilus (A. F. El-Sheikh and M. G. Klotz, unpublished results).
The genes encoding the subunits of DNA-directed core RNA polymerase (EC 2.7.7.6) did not reside in a single gene cluster in the N. oceani genome. Whereas the genes encoding the beta and beta' subunits were arranged in tandem (Noc_2331 and Noc_2330, respectively), the gene encoding the alpha subunit resided upstream of this tandem as gene Noc_2300. A gene encoding an omega subunit (COG1758) was also found (Noc_1213). The genome also contained a variety of genes that encode alternative sigma factors, some of which were found in multiple copies (see below).
Metabolism and transport: genomic basis of ammonia lithotrophy. (i) Energy metabolism: acquisition of reductant from the environment.
A cluster of three contiguous genes encoding the subunits of ammonia monooxygenase (Noc_2503 to Noc_2501) were found to be organized in overlapping operons as described previously (1, 59). No additional functional amo genes or amo pseudogenes were found in the genome. As reported previously, the terminator of the amo operon was succeeded by a transcriptional unit containing the orf5 gene (Noc_2500) (59). An additional orf5-like gene was found as an orphan in the genome (Noc_3006). The deduced Orf5 protein sequence revealed a signal peptide and an additional membrane-spanning domain at its C terminus. Such orf5 genes with a high degree of sequence identity and in conservation of synteny have also been identified in betaproteobacterial nitrifier genomes downstream of the amo operon (U92432 and AF016003); however, they have been identified as one of two sequence-related genes in an orf45 transcriptional unit (16). Interestingly, an orf5 homologue was also found in the whole-genome sequence of the gammaproteobacterial methanotroph Methylococcus capsulatus (Bath) (MCA2130), where it was not in proximity of either of the two gene clusters that encode particulate methane monooxygenase (pMMO), a homologue of AMO (29, 59). Instead, the orf5 gene was resident in a transcriptional unit together with a gene that encodes a pan1-type multicopper oxidase (MCO; MCA2129), which is likely a homologue to the pan1-type MCO-encoding gene found upstream of the amo operon in the N. oceani genome (Noc_2506). Because both AMO and pMMO can oxidize ammonia to hydroxylamine and because of the sequence and genome organizational similarities between the associated mco and orf5 genes in the N. oceani and M. capsulatus genomes, a functional role of their membrane-associated expression products in transfer of electrons or intermediates related to ammonia oxidation is proposed.
The oxidation of hydroxylamine is the core of ammonia catabolism, as it provides electrons for redox-dependent proton pumping (cytochrome bc1 complex, terminal cytochrome c oxidase). Cytochrome P460 (Noc_0890) is a likely ancient hydroxylamine (and nitric oxide) dehydrogenase suited to detoxify the mutagenic hydroxylamine at low turnover rates (9, 10). However, the electrons gained during the oxidation process cannot be transferred directly to the cytochrome bc1 complex in the membrane and are likely relayed via the periplasmic soluble cytochrome c552 (Noc_0751) to the terminal oxidase. While this pathway contributes to the proton-motive force, it is an inefficient use of reductant. A more efficient pathway occurs when hydroxylamine is catalytically oxidized by hydroxylamine oxidoreductase (HAO; Noc_0892) and electrons are funneled via cytochromes c554 (Noc_0894) and cm552 (Noc_0895) to the ubiquinone pool (Q/QH2; Noc_1248 to Noc_1252). The structure and sequence conservation of the HAO-c554-cm552 pathway in N. oceani was recently evaluated (7), based in part on the whole-genome sequence reported in this paper. The reduced quinone pool provides reductant to the AMO complex, to the cytochrome bc1 proton-pumping complex (Noc_0297 to Noc_0299) in the electron transport chain, and to the NUO complex I responsible for reverse electron transport needed to generate NADH (see below). In addition, numerous uncharacterized MCO-encoding genes were found in the genome (Noc_0889 and Noc_2605, type 1 MCOs; Noc_1542, exported MCO; Noc_1741, copper resistance protein), which need further characterization to assess their involvement in catalytic and electron transfer processes of catabolic pathways in N. oceani. In addition, the finding of a gene encoding the red copper protein nitrosocyanin (Noc_1090) in the N. oceani genome suggests that this protein is unique to and important for ammonia catabolism.
The hydrolysis of urea to ammonia and carbon dioxide can be carried out in N. oceani by the ATP-independent heteromultimeric nickel enzyme urea-amidohydrolase (urease) (40). Ureolysis could be beneficial to the cell because it produces N. oceani's sole sources for energy and reductant (ammonia) as well as carbon (CO2). Given the low concentration of dissolved urea in the oceans, the role of urease in this organism is unclear.
Some betaproteobacterial nitrifiers can utilize H2 (12). However, the genome of N. oceani did not contain any genes that encode subunits of a hydrogenase, which is in agreement with the fact that Nitrosococcus cannot grow on H2 as the sole source of energy and reductant.
(ii) Energy metabolism: electron flow, generation of universal reductant, generation of the proton gradient, and ATP production.
The gene profile for N. oceani reveals complete sets of genes for electron transfer from NADH to O2 via NADH quinone oxidoreductase (complex I), cytochrome bc1 complex (complex III), and a Cu-aa3-type cytochrome c oxidase (complex IV). Coupled with genes for a complete tricarboxylic acid (TCA) cycle and glycolytic pathway, it seems that N. oceani has the potential to gain energy through the oxidation of organic compounds. A complete pentose phosphate pathway provides an alternative mechanism for oxidizing sugars and generating NAD(P)H. It appears that the obstacle to an organotrophic mode of catabolism is N. oceani's inability to import suitable organic substrates. When using ammonia as the energy source (lithotrophy), there is a need to generate NAD(P)H through reverse electron flow, and a complete TCA cycle is not needed. The complexes involved in electron flow, generation of reductant and the proton gradient, and production of ATP are described.
Two complete yet different sets of genes encoding complex I (NDH-1/NADH quinone oxidoreductase) are present in the genome of N. oceani. Genes Noc_1115 to Noc_1127, which encode one copy of complex I, are most similar to complex I genes found in other Gammaproteobacteria. This operon contains only 13 genes; the c and d subunits are fused into a single gene. The second set of genes encoding complex I (Noc_2552 to Noc_2565) includes genes with top BLAST hits to Nitrosomonas europaea (6 genes), to other Betaproteobacteria (3 genes), and to Gammaproteobacteria (5 genes). The role of these distinct complexes in N. oceani is unknown. It may be that the complex with the strongest similarity to that found in Nitrosomonas europaea plays a role in reverse electron flow when ammonia is the sole electron donor, whereas the gammaproteobacterial complex may be important in forward electron flow associated with NADH oxidation. The N. oceani genome also contains the genes needed to encode a Na+-transporting NADH:ubiquinone oxidoreductase (Noc_0970 and Noc_1170 to Noc_1174).
The presence of candidate genes for a Na+-dependent complex I (Noc_0970 and Noc_1170 to Noc_1174), a Na+-dependent V-type ATPase (Noc_2081 to Noc_2089), and several Na+/H+ antiporters (Noc_0159, Noc_0521, Noc_1282, Noc_2134, and Noc_2952) raises the possibility of a sodium circuit in addition to the proton circuit in N. oceani, which is likely an adaptation to N. oceani's high-salt environment. Under chemolithotrophic conditions, reverse operation of the sodium-dependent complex I in the plasma membrane could generate additional NADH; however, the sodium-dependent ATPase would be needed to remove excess sodium from the cytoplasm. Under mixotrophic conditions, whether lithoheterotroph or organoautotroph, a sodium gradient might help provide at least some of the organic needs of the cell through import (4). Our finding of several sodium-dependent transporters in the genome (Noc_0779, Noc_1365, Noc_1575, Noc_1600, Noc_2446, and Noc_2711) support this possibility; however, operation of the sodium-dependent ATPase would, again, be needed to remove excess sodium from the cytoplasm. If import of certain organics enabled N. oceani of organotrophy, additional NADH could be generated. The additional NADH produced by the sodium circuit (Na+-dependent complex I and transporters) on the expense of ATP generated by the proton circuit could serve the "forward" complex I as an extension of the bacterium's electron transport chain, thereby effectively converting a sodium-motive force into a proton-motive force. This could be of utility for chemotaxis and proton-dependent transport.
The genes encoding a ubiquinol-cytochrome c reductase (the cytochrome bc1 complex) are located in an operon (Noc_0297 to Noc_0299). Nitrosococcus oceani has genes encoding two complete terminal cytochrome c oxidases (COX), both of the Cu-aa3 type. One complex (Noc_3044 to Noc_3047) appears to be of gamma- and betaproteobacterial descent (Pseudomonas/Nitrosomonas), and the other (Noc_1244 to Noc_1247) of Bacteroidetes/Chlorobi decent. There are two additional copies of the genes for subunit 1 and subunit 2 but not for subunit 3. In addition, there are eight copies of genes encoding class I c-type cytochromes.
Additional electron sinks are also encoded in the genome, i.e., there are three di-heme cytochrome c peroxidases (Noc_0488, Noc_1263, and Noc_2697). A partial denitrification pathway is present, including a Cu-type nitrite reductase (nirK; Noc_0089) and nitric oxide reductase (nor; Noc_1847 to Noc_1851). However, genes for nitrate reductase (Nar) and nitrous oxide reductase (Nos) were not identified.
Under lithotrophic conditions, the proton gradient is presumably generated by the action of the quinol-cytochrome c oxidoreductase and cytochrome c oxidase. Under organotrophic conditions, we assume that one or both of the NADH-ubiquinone oxidoreductases could also contribute to the generation of the proton gradient. In addition, gene Noc_1901 encodes a proton-translocating inorganic pyrophosphatase which, in the presence of a flux of pyrophosphate, could contribute to the generation of a proton gradient (Fig. 3). Alternatively, the enzyme could use the proton gradient to generate pyrophosphate necessary for glycolysis and other processes. A potential source of pyrophosphate could be polyphosphate, since a polyphosphate kinase is present (Noc_2388). ATP formed from degradation of polyphosphate could release pyrophosphate through the action of one of the many nucleoside diphosphate hydrolases (NUDIX hydrolases) encoded in the genome (Noc_0193, Noc_0306, Noc_2018, Noc_2420, Noc_2512, Noc_2643, and Noc_2749).
![]() View larger version (21K): [in a new window] |
FIG. 3. Proposed roles of polyphosphate and pyrophosphate in the Nitrosococcus oceani cell. The figure illustrates the mechanisms for phosphate uptake and the proposed flow of phosphate, pyrophosphate, and polyphosphate and their involvement in energy metabolism (ATP, pmf), central pathways (PP-Pfk), and stress tolerance (copper detox) in the N. oceani cell. The proposed phosphate-proton exchanger to remove polyphosphate-copper complexes has yet to be experimentally identified.
|
(iii) Energy storage strategies: sucrose synthase and sucrose-phosphate synthase.
Nitrosococcus oceani has genes encoding a sucrose synthase (SuSy; Noc_3068) and sucrose phosphate synthase (SPS; Noc_3069). Sucrose synthase in plants functions predominantly in the degradation of sucrose (47). Sucrose phosphate synthase functions in plants to synthesize sucrose phosphate, and sucrose phosphate phosphatase completes the synthesis of sucrose in plants. A separate gene for sucrose phosphate phosphatase is not present in the genome of N. oceani; however, gene Noc_3069 is a fusion of the genes for sucrose phosphate synthase and sucrose phosphate phosphatase. The haloacid dehalogenase triad conserved in this superfamily of glycohydrolases is encoded by gene Noc_3069. UDP-glucose is the glucosyl donor for both enzymes, and gene Noc_2280 encodes a UTP glucose-1-phosphate uridylyltransferase. Gene Noc_3067 encodes a fructokinase that could make the fructose-6-phosphate required by sucrose phosphate synthase.
Sucrose synthesis in bacteria is not well documented or understood (47). Sucrose and sucrose-synthesizing activities were detected in two species of halotolerant methanotrophs, Methylobacter alcaliphilus 20Z and Methylobacter modestohalophilus 10S (37). Genes for sucrose synthesis have been identified, so far, in just a few Proteobacteria, including Acidithiobacillus ferrooxidans and Nitrosomonas europaea. The genes are also present in cyanobacteria, where their function in sucrose synthesis has been demonstrated. Their role in N. oceani is unknown. Perhaps sucrose can serve as an osmoprotectant, as has been suggested for other prokaryotes.
(iv) Energy storage strategies: carbon storage products.
The genome of N. oceani also contains genes encoding the five enzymes necessary for the synthesis of glycogen from fructose-6-phosphate. The genes encoding glucose-1-phosphate adenyl transferase (Noc_0905) and 1,4-alpha glucan branching enzyme (Noc_0904) are contiguous. The gene for phosphoglucomutase was identified in N. oceani (Noc_1719), but genes necessary for the synthesis of poly-beta-hydroxybutyrate do not appear to be present in the genome of N. oceani.
(v) Central carbon metabolism.
The gene profiles are consistent with complete pathways for glycolysis (from phosphorylated sugars or glucose) and gluconeogenesis. However, the mechanism for the interconversion of fructose-6-phosphate and fructose-1,6-bisphosphate is not clear. Gene Noc_0021 encodes fructose-1,6-bisphosphatase. A candidate for an ATP-dependent phosphofructokinase is not present, but gene Noc_2846 shows some similarity to genes encoding diphosphate-fructose-6-phosphate 1-phosphotransferase (EC 2.7.1.90; pyrophosphate-dependent phosphofructokinase). Both a membrane-bound, proton-translocating pyrophosphatase (Noc_1901) and a soluble pyrophosphatase (Noc_1134) are encoded, providing a mechanism to hydrolyze pyrophosphate generated during gluconeogenesis. The proton-translocating enzyme could also use energy from the proton gradient to synthesize pyrophosphate necessary for glycolysis.
The genes encoding a complete pentose phosphate pathway are also present. We cannot determine if glucose-6-phosphate dehydrogenase would couple to NADP+ (as is most often the case) or NAD+ (which occurs less frequently). A transhydrogenase, encoded by genes Noc_0261 and Noc_0262, would allow transfer of reductant from NADH to NADP+ and from NADPH to NAD+.
All the genes necessary for a complete TCA cycle are present. Of particular interest was the presence of genes encoding alpha-ketoglutarate dehydrogenase, because this enzyme activity was missing in several obligate chemolithotrophs, including Nitrosomonas europaea (30). The absence of alpha-ketoglutarate dehydrogenase activity would create a break in the TCA cycle that would preclude organotrophy; however, if present, the role of alpha-ketoglutarate dehydrogenase in these exclusively or predominantly lithotrophic organisms is not known. As was the case for Nitrosomonas europaea (16), the genes encoding the three subunits (E1 to E3) of alpha-ketoglutarate dehydrogenase were identified in the N. oceani genome, and they were contiguous (Noc_0111 to Noc_0113). Because pyruvate dehydrogenase catalyzes a reaction mechanistically similar to that of alpha-ketoglutarate dehydrogenase and has a similar subunit structure, one gene (lpd) encodes the E3 subunit for both enzymes in many organisms. However, in N. oceani, pyruvate dehydrogenase is encoded by separate genes, Noc_1254 to Noc_1256, which include the gene for subunit E3.
Amino acid and nucleotide metabolism.
The N. oceani genome contains genes encoding the biosynthesis of the 20 amino acids required for the synthesis of proteins. Unless otherwise indicated, all amino acids mentioned were the L form. The genes identified indicate most pathways are similar to previously identified synthesis pathways (5, 26, 68, 69, 77), and the genes encoding the enzymes have values highly similar to those found in other Proteobacteria, most of them with best matches to Methylococcus capsulatus or within the pseudomonads. While most biosynthetic pathway elements have been identified, missing enzymatic steps are typically involved with dual-function enzymes that may be difficult to identify based on sequence alone (for example, Noc_0176 encodes 3-phosphoshikimate 1-carboxyvinyltransferase with dual functions [EC 1.3.1.1 and EC 2.5.1.19]). In contrast, few amino acid degradative enzymes (catabolic enzymes) or specific transport genes were identified. Scavenging from leucine, valine, and isoleucine may be possible. Several transaminases were identified, but most could not be assigned to specific amino acids. Special attention was paid to amino acid biosynthetic functions found in the last common ancestor, as indicated by their universal distribution in the three domains of life (45).
Several large amino acid biosynthetic operons were identified, including a mixed-function supraoperon similar to that found in Pseudomonas encoding aromatic amino acid biosynthesis (76) and the his operon (22). The aromatic supraoperon is found encoded by Noc_0172 to Noc_0177. The dual-function chorismate mutase/prephenate dehydratase P protein is encoded by gene Noc_0174. The histidine biosynthetic genes are not in a single cluster and are located in the genome as genes Noc_2778 and Noc_2779 (hisDG) and Noc_3051 to Noc_3057 (hisCBHAFIE and hitA). Split organizations of the his operon are relatively common and were also found in Nitrosomonas europaea (16, 22). The genes hisI and hisE overlap by 8 bp but do not appear to be fused. The gene encoding histidinol phosphatase (EC 3.1.3.15; Noc_0374) was found outside of the operon.
In the genome of N. oceani, 20 aminoacyl-tRNA synthetases (AARS) were identified, including two forms of LysRS (class I and II) and two distinct forms of GlxRS (EC 6.1.1.17), but AsnRS (EC 6.1.1.22) and GlnRS (EC 6.1.1.18) were missing. The most common organism for the top match for the AARS was Methylococcus capsulatus, as is the case for the overall taxonomic distribution of top matches. The class II aminoacyl-tRNA synthetases for Phe and Gly have two nonidentical subunits, and genes for both the alpha and beta subunits have been identified adjacent to each other in the genome. While no AARS was identified as the specific GlnRS type (EC 6.1.1.18), it is likely that this function is mediated by the product of gene Noc_0264 or Noc_2250, the nondiscriminatory-type GlxRS. Genes encoding a possible glutamyl-tRNA-Gln amidotransferase (gatCAB) were found in the genome, although gatB (Noc_2014) is separated from gatCA (Noc_2635 and Noc_2636). So an indirect route for synthesis of glutaminyl-tRNA is probable. The sequence data alone are insufficient for specifically designating either GlxRS gene a discriminatory role. The presence of two genes does not appear to be a recent duplication event, as the peptides are more similar to their homologues in Methylococcus than to each other (only 44% identical to each other over 315 aligned residues versus 55 and 67% identical to the Methylococcus peptides).
No candidate for AsnRS was identified in the genome. The AspRS encoded by gene Noc_0302 does contain the GAD domain typically found in AspRS involved in the indirect transamination route to Asn-tRNAASN synthesis. The gatCAB-encoded transamination function has a likely role in this pathway. Interestingly, the genome contains several distinct asparagine synthetase genes of the glutamine-dependent asnB type (Noc_0777, Noc_1965, Noc_1975, and Noc_2478). It is somewhat surprising but not without precedent (51) that multiple asnB genes are present in an organism without an identified asnRS gene.
Although typically individual organisms contain only one class of LysRS, genes encoding both classes of LysRS have been identified in the genome of N. oceani. This case of LysRS existing in both class I and class II forms in the genome of N. oceani is the only known bacterial exception to the AARS "class rule" (60). The only known examples of the presence of both classes of LysRS in a single genome are in the archaea Methanosarcina barkeri and M. acetivorans. In N. oceani, the class I LysRS is encoded by gene Noc_1618, with the best match to M. barkeri LysRS. Both selective retention and horizontal gene transfer have played roles in the distribution of class I LysRS in bacteria (2). The gene Noc_2625, encoding 356 amino acids, represents the C-terminal region of class II LysRS with an intact core domain (the best match is to Pseudoalteromonas haloplanktis TAC125 [CAI85572]). Several other Gammaproteobacteria have a similar protein, as described for COG2269. In N. oceani this combination may be an example of gene displacement of the LysRS class II by a functional LysRS class I gene (33). Alternatively, in M. barkeri both forms are involved in the incorporation of the rare amino acid pyrrolysine into the enzyme monomethylamine methyltransferase. This is the first report of the presence of class I and class II LysRS genes together in a bacterial genome (62).
Genes for the synthesis of all five purine and pyrimidine nucleotides are present. In contrast, genes for degradation are very limited. There appears to be no capacity to catabolize nucleotides with the exception of uridine, which can only be converted to pseudouridine.
(vi) Transport.
Approximately 263 ORFs in N. oceani, about 9% of the total, are involved with transport. Included are P-P-bond-hydrolysis-driven transporters, electrochemical-potential-driven transporters, and channels/porins from a large number of protein families (Table 3). Predominant among these are 23 ATP binding cassette (ABC)-type transporters (85 genes) for a variety of organic and inorganic substrates, the resistance-nodulation-cell division family (26 genes) of H+ antiport-driven efflux transporters, and a large number of genes (22) involved with iron transport.
|
View this table: [in a new window] |
TABLE 3. Transporters in the genome of N. oceani ATCC 19707
|
In addition to iron transport, a number of uptake systems (67 genes) for other inorganic ions were identified. Inorganic N may be imported either as nitrite via a formate-nitrite (FNT family) transporter (Noc_0109) or via a putative ammonia/ammonium permease (Amount) (Noc_2700 and Noc_2701). Sulfate could be imported via a sulfate permease (SulP family) (Noc_1626), which may function as an SO42:H+ symport or an SO42:HCO3 antiport, or by a sodium:sulfate symporter (Noc_1175). Phosphate transporters included two ABC transporters (Noc_2396 to Noc_2399 and Noc_0581 to Noc_0584) and two phosphate-selective porins (Noc_2417 and Noc_2418) (see below). A number of metal ion uptake and efflux systems were identified, including an ABC transporter system dedicated to Mn2+/Zn2+ transport (Noc_2421 to Noc_2423), a CorABC-type Mg2+/Co2+ ion channel of the mitochondrial tricarboxylate carrier family (Noc_0240, Noc_1416, and Noc_2263), three divalent cation transporters for Mg2+/Co2+/Ni2+ of the MgtE family of magnesium transporters (Noc_1840, Noc_1785, and Noc_2801), and three divalent heavy-metal cation transporters (Noc_0092, Noc_04234, and Noc_1342). Efflux systems include two small multidrug resistance (SMR) family transporters of cations and cationic drugs (Noc_0601 and Noc_2235), CopCD copper export proteins (Noc_1741 and Noc_1742), and four cation diffusion facilitator-type cation efflux proteins (Co2+/Zn2+/Cd2+) (Noc_0595, Noc_1534, Noc_1782, and Noc_2871). Other ion transporters include chloride channel protein EriC (Noc_0358) and an major facilitator superfamily (MFS)-type cyanate transporter (Noc_1456). Two P-type cation-transporting ATPases were found (Noc_1406 and Noc_2130).
Likely of importance to N. oceani's salty habitat in the oceans are a number of monovalent cation transporters. Several Na+/H+ antiporter systems were found, including an NhaD type (Noc_1492), an NhaC type (Noc_2134), and two NhaP types (Noc_0159 and Noc_0521). These play important roles in maintaining intracellular pH and conferring salt tolerance. Five genes of the small conductance mechanosensitive ion channel (MscS) family (Noc_0108, Noc_0602, Noc_1030, Noc_1853, and Noc_1914) were found which play a role in fast osmoregulatory responses. Other transport systems identified included the genes encoding both the NAD-binding component (Noc_3020 and Noc_3064) and the membrane component (Noc_0960, Noc_2194, and Noc_2952) of the KefB-type K+ transport systems, a Trk-type K+ uptake system (Noc_0242, Noc_1639, and Noc_1640), and two DASS (divalent anion/sodium ion symporter) family members (Noc_1175 and Noc_2446).
The genome of N. oceani revealed only a few transporters for importing organic compounds. ABC-type transporters may import dipeptides/oligopeptides (Noc_1344, Noc_1768, and Noc_2770 to Noc_2773), proline/glycine betaine (Noc_0539 and Noc_0540), and unspecified sugars (Noc_0279 and Noc_0282). Other transporters included an amino acid transporter (Noc_3063), a lactate permease (Noc_1578), and a urea transporter (Noc_2884). A number of potential transport systems for dicarboxylates were found, including a TRAP-type C4-dicarboxylate permease (Noc_0598, Noc_0709, and Noc_0710), a tellurite-resistance/dicarboxylate transporter (TDT) (Noc_0077 and Noc_0542) which may be involved in transporting dicarboxylic acid intermediates, and a DAACS (di-/tricarboxylate/amino acid:cation symporter) family members (Noc_1175 and Noc_2446).
Export systems for organic and toxic compounds included ABC transporters for organic solvents (Noc_1746, Noc_2096, and Noc_2782 to Noc_2785), multidrug resistance (Noc_1779 to Noc_1781, Noc_2141, Noc_2142, Noc_1835, Noc_1836, Noc_2644 to Noc_2646, Noc_2832, and Noc_2833), heme export (Noc_0946 to Noc_0948), dipeptides (Noc_0933 to Noc_0936, Noc_2184, and Noc_2185), polar amino acids (Noc_1558 and Noc_1559), MFS-type arabinose efflux transporters (Noc_1547, Noc_1759, Noc_2803, and Noc_3022), colicin export (Noc_0144, Noc_0145, and Noc_2616), polysaccharide/polyol-phosphate (Noc_1227, Noc_1228, Noc_2181, and Noc_2182), lipoprotein export (Noc_1272, Noc_1273, Noc_2214 to Noc_2216, Noc_2669, and Noc_2670), as well as a transporter involved in lysophospholipase L1 biosynthesis (Noc_0427 and Noc_0428). Unlike the uptake ABC transporters, these ABC transporters generally lacked any associated periplasmic component.
Nitrosococcus oceani possesses several protein export and secretion systems, including a preprotein translocase (Noc_2305), genes for exporting folded redox proteins via the sec-independent protein secretion system (TatABC; Noc_3058 to Noc_3060), genes encoding the type II general secretion/pilus synthesis pathway, and genes encoding the type IV conjugal DNA/protein transfer system.
Nitrosococcus oceani has a cluster of genes encoding parts of a phosphotransferase (PTS)-type sugar transport system. Specifically, genes for E1 (Noc_2800), HPr (Noc_2799), an Ntr-type IIa (Noc_2795), a mannose/fructose type IIa (Noc_2313), as well as an HPr kinase/phosphorylase (Noc_2796) were found. However, genes encoding components IIb, IIc, and IId were absent, making it unlikely that N. oceani expresses a functional sugar-transporting PTS system. Instead, as has been seen in other Proteobacteria, these genes may form part of a regulatory cascade involving RpoN (13).
Growth, sensing, responses, and stress tolerance. (i) Carbon fixation.
The N. oceani genome encodes a form I ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO). The form I genes (cbbL and cbbS) occur in an operon with cbbX and a hypothetical gene (Noc_0330 to Noc_0333). The operon is preceded by a divergently transcribed cbbR gene (Noc_0334) encoding a LysR-type transcriptional regulator. The cbbX gene is required for efficient autotrophic growth in Rhodobacter sphaeroides (24) and is predicted to be an AAA family ATPase (which can be involved in chaperonin-like functions).
Both prokaryotic/plant-type (Noc_1341) and eukaryotic-type (Noc_1132) carbonic anhydrases are encoded by the genome. There is no obvious candidate for a bicarbonate acquisition system, however. Furthermore, the genome of N. oceani lacks genes for carboxysome formation.
Genes for all enzymes to complete the Calvin-Benson-Bassham cycle are present. Transketolase, NAD-dependent glyeraldehyde-3-phosphate dehydrogenase (EC 1.2.1.12), phosphoglycerate kinase, pyuvate kinase, and fructose-1,6-bisphosphate aldolase are encoded by an operon (Noc_2804 to Noc_2808), whereas fructose-1,6-/sedoheptulose-1,7-bisphosphatase (Noc_0021), ribose-5-phosphate isomerase (Noc_2667), and phosphoribulokinase (Noc_2826) are encoded by isolated genes. Genes encoding ribulose-5-phosphate 3-epimerase and phosphoglycolate phosphatase are grouped in an operon (Noc_2492 to Noc_1493).
Due to the aforementioned homology of ammonia monooxygenase and particular methane monooxygenases as well as the hypothesized lateral transfer of genes, the genome of N. oceani was investigated for the presence of C1 metabolic pathways (see Table 4). Whereas the N. oceani genome does not contain the genes needed to fix carbon directly from methane or methanol via the ribulose monophosphate or serine pathways, it contains the genetic inventory to funnel C1 carbon into the Calvin-Benson-Basham cycle (Table 4). Remarkably, two independent pathways for formaldehyde oxidation to formate (Noc_1394, Noc_1440, and Noc_2006) and two gene clusters encoding the capacity for formate dehydrogenation (Noc_1122 to Noc_1124 and Noc_2559 to Noc_2561) were identified. The identification of this inventory provides a theoretical explanation of earlier findings that labeled C1-carbon was assimilated into the biomass of N. oceani cultures growing in an ammonia-dependent fashion (36, 71).
|
View this table: [in a new window] |
TABLE 4. C-1 carbon metabolism in the genome of N. oceani ATCC 19707 in comparison with the genomes of the methanotroph M. capsulatus (Bath) and the beta-AOB N. europaea
|
As a nitrifying bacterium, N. oceani affects the nitrogen cycle in its marine environment by assimilatory and dissimilatory activities (Fig. 4). While the ammonia oxidation capacity (AMO and HAO; see above) of AOB links the oceanic pools of reduced (ammonia) and oxidized (nitrite) nitrogen, their classical denitrification capacity through dissimilatory nitrite reductase (NirK; EC 1.7.2.1; Noc_0089) and nitric oxide reductase (Nor; Noc_1847 to Noc_1851) is likely a major source of nitrous oxide emitted from the oceans (20, 21, 57). Additionally, there are various genes that encode putative cytoplasmic (Noc_0889 and Noc_2605) and exported (Noc_1542) multicopper oxidases (MCO), as well as other MCOs that have been implicated in oxidation of NOx, such as the pan1-type MCO (Noc_0889), an alternative NOx reductase (6). In contrast to its organization in the genome of Nitrosomonas europaea, where it is clustered in a four-gene operon with the NirK-type nitrite reductase (6, 16), the pan1-type MCO gene Noc_0889 is clustered, but not necessarily in the same transcriptional unit, with a gene encoding cytochrome P460 (Noc_0890) directly upstream of the hao gene cluster (7). Cytochrome P460 has been implicated in hydroxylamine detoxification in several bacteria, including N. oceani (9, 10, 31, 79). Because cytochrome P460 can be reduced by NO and because of the physical proximity of the two genes in the N. oceani genome, P460 and pan1-type MCO may be components of an alternative N-oxidation pathway that, in contrast to HAO, produces nitrite with NO as an intermediate. Considering the toxicity of NO and the fact that AOB have a high O2 consumption rate but tend to live near the oxic/anoxic interface, additional NO detoxification mechanisms that are operational under low-O2 conditions are likely necessary for survival of AOB (75). It is thus not surprising to find in the genome of N. oceani a gene, cycP, that encodes a beta-sheet-structured cytochrome c' ("c'-beta"; Noc_2696). Alpha-helical c' cytochromes have been implicated in microaerobic NO sequestration and dehydrogenation by strains of the betaproteobacterial pathogen genus Neisseria (18, 19, 52, 64, 67). It has been proposed only recently that cytochrome c'-beta, which is evolutionarily related to and a putative redox partner of cytochrome P460 (8), has evolved from an alpha-helical monoheme cytochrome c that is ancestral to both c' and c'-beta cytochromes (M. G. Klotz and A. B. Hooper, unpublished results). Both enzymes, c'-beta and P460, are also likely redox partners of the soluble periplasmic cytochrome c552 (Noc_0751), as are the three periplasmic di-heme cytochrome c peroxidases. The P460-heme coordination site, found in HAO and cytochrome P460, is highly sensitive to hydrogen peroxide (32); therefore, it is interesting that the c'-beta-encoding gene is clustered with one of the three di-heme cytochrome c peroxidase genes (Noc_2697), suggesting a dedicated protective function of key periplasmic enzymes by this peroxidase. The genome of N. oceani also contains a gene encoding the red copper protein, nitrosocyanin (Noc_1090), a putative enzyme with a cupredoxin fold (3, 55). Because this gene and the encoded protein have so far been uniquely found in a beta-AOB (3, 16) and now also in a gamma-AOB, its putative catalytic function is likely involved in and specific to ammonia-oxidizing catabolism. It appears that the catabolic dependence on ammonia oxidation and the ultimate production of reactive and toxic NOx intermediates have imparted selective pressure on all AOB to maintain this suite of periplasmic enzymes, because a similar complement of genes, albeit with different genomic organization, was also identified in the genome of Nitrosomonas europaea (16). It appears that the strategy of maintaining a complement of NOx-detoxifying enzymes in the periplasm is similar to the strategy of active oxygen defense and is designed to avoid the formation and presence of nitrosating agents such as nitrous anhydride in the cytoplasm, where they have mutagenic activity (75).
![]() View larger version (24K): [in a new window] |
FIG. 4. Schematic representation of nitrogen metabolism in Nitrosococcus oceani. The figure illustrates transport and metabolism of nitrogenous compounds and provides involved proteins identified by their gene numbers. Nitrogen is either acquired as an inorganic compound in the form of ammonia/ammonium or as an organic compound in the form of urea, which is hydrolyzed by urea hydrolase (UreABC) to ammonia and carbon dioxide. For synthesis, nitrogen can be assimilated from ammonia into glutamate via GDH (GdhA) or the GS-GOGAT (GlnA, GltBD) pathway. For catabolism, ammonia is oxidized by AMO and HAO and the extracted electrons are used to reduce the quinone pool via cytochromes c554 and cm552. Some of these electrons are employed to generate pmf by the cytochrome bc1 complex, and these electrons are then relayed to a terminal electron acceptor (COXaa3, Cu-NIR, cytochrome [Cyt] c', NORbb3, etc.) via cytochrome c552. The question marks indicate that the active site of ammonia monooxidation by AMO needs to be elucidated and that the annotation of AmtB-type ammonia transport genes needs experimental verification.
|
As a member of the Chromatiaceae, N. oceani should participate in the sulfur cycle beyond acquiring sulfur for biosynthesis. The identified sulfate uptake capacity allows N. oceani to acquire and process sulfate. Sulfate reduction may proceed via sulfate adenylyltransferase (cysND; Noc_2288 and Noc_2289; EC 2.7.7.4) to adenosine phosphosulfate via adenylylsulfate kinase (cysC; Noc_2482; EC 2.7.1.25) to phosphoadenosine phosphosulfate and via phosphoadenosine phosphosulfate reductase (cysH; Noc_2290; EC:1.8.4.8) to sulfite, which may be further reduced to H2S by an NADPH-dependent sulfite reductase, EC 1.8.1.2 (alpha-subunit CysI, Noc_1305; beta-subunit CysI, Noc_1306). Dihydrogen sulfide is required for cysteine biosynthesis, and genes encoding a thioredoxin-disulfide reductase (Noc_0345; EC 1.8.1.9) and a thiol:disulfide interchange protein (Noc_0551) were identified, as was the gene encoding thiosulfate sulfurtransferase (Noc_0593; EC 2.8.1.1). The genome of N. oceani also contains a gene cluster that encodes a putative polysulfide reductase (psr; Noc_1238 to Noc_1240), a monoheme cytochrome (cccA, Noc_1241), a transporter (Noc_1242), a cytochrome c oxidase (Noc_1244 to Noc_1247), and the five genes encoding the ubiquinone complex (Noc_1248 to Noc_1252). This cluster of genes is absent from the genomes of Nitrosomonas europaea,Nitrosomonas eutropha, Nitrospira multiformis, and Nitrobacter winogradskyi but is conserved (sequence and synteny) in the genomes of Nitrobacter hamburgensis, Rhodopirelulla baltica, Cytophaga hutchinsonii ATCC 3406, and Cupriavidus necator (Ralstonia eutropha JMP134). Functional polysulfite reductase is a molybdopterin oxidoreductase complex that has been experimentally described to act as a quinole oxidase in Wolinella succinogenes (42, 43). Analysis of the N. oceani genome did not reveal a molybdopterin guanine dinucleotide-binding protein subunit-encoding gene in the vicinity of the cluster. This putative polysulfide reductase activity awaits experimental verification in N. oceani, which would indicate the residence of an alternative catalytic center for polysulfide reduction. If present, this molybdopterin oxidoreductase could theoretically also be involved in the anaerobic reduction of nitrate, chlorate, selenate, or other highly oxidized minerals, thereby accommodating electron disposal in a microaerophilic environment near the oxic/anoxic interface.
In contrast to many other purple sulfur bacteria, the genome of N. oceani lacked genes for the formation of internal or external granules of sulfur compounds.
(iii) Cellular growth and motility.
The genome of N. oceani contains almost all of the typical complement of genes with an identified role in cell cycle and division of other Gammaproteobacteria, such as Escherichia coli. The genome clearly lacks the genes encoding FtsEX (involved in localization and stabilization of the septal ring), FtsN, and SulA. On the other hand, the genome contains three genes (Noc_0272, Noc_1903, and Noc_2569) encoding proteins with domains (conserved zinc-binding motif HEXXH, ATPase domain, and peptidase domain) matching the cell division metalloprotease FtsH (COG0465). An alignment of all three FtsH proteins showed that they differ in sequence mostly at their N and C termini, which flank the ATPase and peptidase domains.
Genes encoding flagellation and motility are in the N. oceani genome in two large clusters (Noc_2354 to Noc_2378 and Noc_2155 to Noc_2166) and several smaller clusters (Noc_0833, Noc_0834, Noc_0124 to Noc_0131, Noc_2052, Noc_2053, and Noc_2683 to Noc_2685). The master switch operon flhCD was not identified and is likely absent, as it is from the genomes of other AOB. As a likely adaptation to life in the open ocean, N. oceani appears to have only limited chemotactic capacity, because just one methyl-accepting chemotaxis protein (MCP) (Noc_0128) of the PilJ type was identified. In contrast, the genome of Nitrosomonas europaea contained three MCPs (16). The presence of a sodium-driven polar flagellar motor protein (MotA; Noc_0833), which can assemble with the product of the adjacent pomB gene (Noc_0834) in addition to the usual pmf-dependent flagellar rotation mechanism, may be an adaptation to N. oceani's marine lifestyle.
(iv) Two-component systems.
Considering the reductive evolution of AOB as concluded from analysis of the Nitrosomonas europaea genome (16), the genome of N. oceani harbored an impressive complement of complete two-component systems. The genome contained 13 genes encoding histidine protein kinases (HPK), of which 12 were paired (mostly succeeded) by a response regulator (RR). In addition, 1 HPK- and 11 RR-encoding orphaned genes were identified (Table 5 and Table S3 in the supplemental material). One of the HPK genes (Noc_1756) was succeeded by a tandem of two RR genes (Noc_1757 and Noc_1758).
|
View this table: [in a new window] |
TABLE 5. Two-component systems in the genome of N. oceani ATCC 19707a
|
(v) Stress tolerance.
In comparison with other AOB, the N. oceani genome contains only a limited inventory that contributes to stress tolerance in general and oxidative stress tolerance in particular. The genome encodes a heme-containing monofunctional large-subunit catalase (KatE; Noc_1165) and an iron-containing superoxide dismutase (Fe-SOD; Noc_2428), both of which are supplied with iron and heme by bacterioferritin (Bfr; Noc_1411). The genome also contains genes encoding glutaredoxin (Noc_2427), thioredoxin (TRX; Noc_0603 and Noc_2583), and a thioredoxin-dependent peroxide reductase (AhpC; peroxiredoxin; Noc_1307), but it lacks genes for bacterioferritin-comigratory protein, NADH-peroxiredoxin reductase (AhpF), glutathione oxidoreductase, and other isozymes of hydroperoxidases (KatG, KatA, and Mn-Cat) and SOD (Mn-SOD and Cu/Zn-SOD). Like Nitrosomonas europaea, the genome also lacks genes for OxyR redox-autoregulatory protein, which regulates oxidative stress tolerance (KatG, AhpC, and SOD), iron and zinc transport proteins (Fur and Zur), and the stationary-phase-specific sigma factor RpoS in many bacteria. In contrast to Nitrosomonas europaea, which lacks an RpoS gene, the genome of N. oceani contains two genes encoding RpoS (sigma 38; Noc_0183 and Noc_1702). RpoS is known to regulate hydroperoxidase (KatE) and the cell shape protein BolA (Noc_2387). Nitrosococcus oceani seems minimally prepared to respond to other stresses. In addition to two genes that encode the minimal growth sigma factor RpoD (sigma 70; Noc_0045 and Noc_2066), the genome contains genes that encode alternative sigma factors involved in heat (RpoH; sigma 32; Noc_1935) and extreme heat (RpoE; sigma 24; Noc_2463) stresses, nitrogen starvation (RpoN; sigma 54; Noc_2793), and the need to move by flagellar motility (FliA; sigma 28; Noc_2155). Despite the absence of the flagellar master operon (fhlCD), FliA-dependent regulation of a complete complement of flagella synthesis and chemotaxis gene clusters is aided by the presence of only one MCP of the pseudomonad PilJ type (with MA and HAMP domains; Noc_0128), whereas other MCPs in the Tar (CheM), Tsr (CheD), Tap, and Aer categories were absent from the genome. In addition, the N. oceani genome did not contain genes with significant similarity to the two-component regulatory systems of LasRI/RhlR involved in homoserine lactone autoinducer synthesis (quorum sensing) as well as the regulation of motility, virulence, starvation response, and iron homeostasis in several Gammaproteobacteria. The ferric uptake regulation protein (Fur; Noc_1194) regulates, for instance, ferric citrate (FecIR) and ferrichrome (fhu operon) transport, exotoxin synthesis, and the expression of hydroperoxidases in many proteobacteria. Interestingly, the zinc uptake regulation protein (Zur; Noc_2424) was found adjacent to a gene cluster that encodes a binding protein-dependent zinc ABC transporter system (Noc_2421 to Noc_2423) in the genome.
Conclusions.
Nitrosococcus oceani is one of only two known ammonia-oxidizing bacteria classified as Gammaproteobacteria, while the large majority of isolated ammonia-oxidizing bacteria are classified as Betaproteobacteria. The genome sequence of the betaproteobacterium Nitrosomonas europaea is available and facilitates a comparison of the genes most similar between these two bacteria. Of the 224 genes in N. oceani that were most similar to genes from Nitrosomonas europaea, 76 were classified as hypothetical or proteins of unknown function. These hypothetical and unknown proteins are slightly overrepresented (34%) in this subset of genes relative to genes without function prediction in the complete N. oceani (31.5%) and Nitrosomonas europaea (29%) genomes. Whereas the number of genes without function and without similarity to other known genes (hypothetical ORFs) is 10 times higher in the Nitrosomonas europaea (4.7%) versus the N. oceani (0.57%) genome, the number of genes without function but with similarity to genes in other genomes (conserved hypothetical ORFs) was significantly higher in the N. oceani genome (30.9%) than in Nitrosomonas europaea (24.4%). For those genes assigned a function, the functions included central carbon and nitrogen metabolism, electron transport, gene regulation, and transport. The Nitrosomonas europaea-like genes were distributed throughout the N. oceani genome. Most of these genes were not flanked by additional Nitrosomonas europaea-like genes, though in some cases two or three such genes were contiguous. A cluster of Nitrosomonas europaea-like genes was present from Noc_1955 to Noc_1986, where 20 of these 31 genes were most similar to Nitrosomonas europaea genes. As additional genome sequences become available, it will be of interest to carry out more detailed comparisons of the genes in common among the ammonia-oxidizing bacteria. Of particular interest will be the uncharacterized genes, some of which may encode functions unique to the use of ammonia as a growth substrate, whereas others may be responsible for the difference between marine and freshwater sediment soil AOB in their abilities to form nitrification consortia (e.g., ammonia- and nitrate-oxidizing bacteria). Progress in ongoing and future research with nitrosococci will be made available at the Nitrosococcus Project website (http://nitrosococcus.org).
Computational annotation was carried out at the Oak Ridge National Laboratory. Our thanks go to high school student Ariella Barhen of Oak Ridge, TN, for help with the annotation. Special thanks to A. B. Hooper (UMN) and D. J. Bergmann (BHSU) and an anonymous reviewer for helpful comments.
Supplemental material for this article may be found at http://aem.asm.org/. ![]()
|
|
|---|
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»