Original Research ARTICLE
D1FHS, the Type Strain of the Ammonia-Oxidizing Bacterium Nitrosococcus wardiae spec. nov.: Enrichment, Isolation, Phylogenetic, and Growth Physiological Characterization
- 1Evolutionary and Genomic Microbiology, Department of Biological Sciences, University of North Carolina at Charlotte, Charlotte, NC, USA
- 2State Key Laboratory of Marine Environmental Science and College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
- 3Joint Research Center for Carbon Sink: The Institute of Marine Microbes and Ecospheres, Xiamen University and the Qingdao Institute of BioEnergy and Bioprocess Technology, Chinese Academy of Sciences, Qingdao, China
- 4School of Marine Science and Policy, College of Earth, Ocean, and Environment, University of Delaware, Newark, DE, USA
- 5Delaware Biotechnology Institute, University of Delaware, Newark, DE, USA
- 6Evolutionary and Genomic Microbiology, Department of Biology and School of Earth and Environmental Sciences, Queens College, The City University of New York, Flushing, NY, USA
An ammonia-oxidizing bacterium, strain D1FHS, was enriched into pure culture from a sediment sample retrieved in Jiaozhou Bay, a hyper-eutrophic semi-closed water body hosting the metropolitan area of Qingdao, China. Based on initial 16S rRNA gene sequence analysis, strain D1FHS was classified in the genus Nitrosococcus, family Chromatiaceae, order Chromatiales, class Gammaproteobacteria; the 16S rRNA gene sequence with highest level of identity to that of D1FHS was obtained from Nitrosococcus halophilus Nc4T. The average nucleotide identity between the genomes of strain D1FHS and N. halophilus strain Nc4 is 89.5%. Known species in the genus Nitrosococcus are obligate aerobic chemolithotrophic ammonia-oxidizing bacteria adapted to and restricted to marine environments. The optimum growth (maximum nitrite production) conditions for D1FHS in a minimal salts medium are: 50 mM ammonium and 700 mM NaCl at pH of 7.5 to 8.0 and at 37°C in dark. Because pertinent conditions for other studied Nitrosococcus spp. are 100–200 mM ammonium and <700 mM NaCl at pH of 7.5 to 8.0 and at 28–32°C, D1FHS is physiologically distinct from other Nitrosococcus spp. in terms of substrate, salt, and thermal tolerance.
Jiaozhou Bay (36°7′24.44′′N 120°14′44.3′′E; Figure 1) is a semi-enclosed water body with a surface area of about 390 km2 and an average water depth of 7 m, located on the Southwestern coast of Shandong Peninsula, China. It connects to the Chinese Yellow Sea with a narrow mouth, about 2.5 km wide (Shen et al., 2006). The bay is bordered by Qingdao City, a large industrial and agricultural city in North China with a population of approximately 9 million residents. Freshwater is provided to Jiaozhou Bay by numerous tributaries (Yang River, Dagu River, Moshui River, Baishahe River, Wantou River, Loushan River, Banqiaofang River, Licun River, Haipo River, etc.). Most of these rivers have become the canals for industrial, agricultural, and urban wastes and are thus highly enriched in nutrients including concentrated dissolved inorganic nitrogen and inorganic phosphate sources (Shen, 2001). One of the frequently sampled locations is Station D1, (Figure 1 and Dang et al., 2008, 2009, 2010a). D1 is close to Huangdao industrial park, which includes one of the three largest crude oil storage facilities of China with the capacity to stockpile more than 3 million tons of crude oil as well as a large oil-refining facility with a designed processing capacity of 10 million tons per year (Shen et al., 2006; Dang et al., 2009, 2010a). Huangdao industrial park also contains one of the largest sea terminals for the shipment of crude and refined oil. Despite the advanced eutrophic state of the water column, the much less impacted sediments revealed a highly diverse community of ammonia-oxidizing bacteria (AOB) that facilitate aerobic (Dang et al., 2010b) and anaerobic (Dang et al., 2010a) removal of ammonium. This was the impetus for us to use sediments from this unique environment to search for unknown representatives of marine AOB.
Obligate aerobic chemolithotrophic AOB utilize ammonium as the sole source for energy and reducing power for growth by oxidizing its non-protonated form, ammonia, to nitrite with hydroxylamine as the intermediate. While initially defined primarily by differences in cell morphology and physiological characteristics (Watson and Mandel, 1971), the advent of 16S rRNA gene phylogeny led to the definition of three genera of AOB (Head et al., 1993), which was largely confirmed by the sequencing of available representatives in pure culture (Arp et al., 2007; Campbell et al., 2011). While the betaproteobacterial AOB (beta-AOB) in the family Nitrosomonadaceae represent two genera with numerous species (Urakawa et al., 2015), Nitrosococcus is the only genus of gammaproteobacterial AOB (gamma-AOB; Campbell et al., 2011). In the past, research of the ecology and phylogeny of AOB has almost entirely focused on the beta-AOB because this group is abundant in all oxic environments, accessible to molecular ecological analysis because the evolutionary histories of enzymes that are key to ammonia catabolism are monophyletic and congruent with 16S rRNA gene-based phylogeny, and beta-AOB usually grow faster than gamma-AOB (Kowalchuk and Stephen, 2001; Campbell et al., 2011). Gamma-AOB in the genus Nitrosococcus are adapted to and restricted to marine environments (Table 1; Ward and O’Mullan, 2002) and strains representing three species are being maintained in pure culture (Campbell et al., 2011).
As recently summarized (Campbell et al., 2011), the type strain of the genus Nitrosococcus was originally described by Winogradsky (1892), renamed ‘Micrococcus nitrosus’ by Migula (1900) and finally established as N. nitrosus (Buchanan, 1925). The binomial was approved as valid by Skerman et al. (1980); however, by then the culture was already lost thereby rendering the genus without a validly described type strain in culture (Campbell et al., 2011). Because Winogradsky enriched N. nitrosus reportedly from freshwater, marine strain D1FHS isolated into pure culture (from station D1 in Jiaozhou Bay, Fermented in a bioreactor and cultured in the presence of Hydroxylamine at 700 mM NaCl instead of the typical 500 mM in mineral Salts media) cannot be assigned to the valid type species of the Nitrosococcus genus. The characterization of its properties indicates that D1FHS is physiologically and taxonomically distinct from other Nitrosococcus spp. leading us to propose it as the type strain of N. wardiae spec. nov. Strain D1FHS has been submitted to the Japan Collection of Microorganisms (JCM) and the American Type Culture Collection (ATCC) for deposit.
Materials and Methods
Sample Collection, Cultivation, and Screening of the Enrichment
Sediment samples were collected from station D1 (120.16438°E, 36.00879°N) and from 13 other stations of Jiaozhou Bay (Figure 1) on October 26, 2008, as described previously by Dang et al. (2009). Briefly, surface sediment subcore samples were collected from 14 stations including D1 (Figure 1) in summer 2008 using a stainless steel 0.05-m2 Gray O’Hara box corer. Replicate samples were taken from inside the box core reaching to a depth of 5 cm, placed in sterile Ziploc bags and stored on dry ice and at 80°C after returning to the laboratory. Water was collected at each of the 14 sampling sites and temperatures, dissolved oxygen, salinity, and pH of all surface sediments were measured (Table 2). Only sediment samples from station D1 were used for this study.
Enrichments were started by mixing 0.5 g of frozen sediment into 100 mL of ammonia mineral salts media [AMS: 700 mM NaCl, 12.5 mM (NH4)2SO4, 30 mM MgSO4, 20 mM MgCl2, 10 mM CaCl2, 10 mM KCl, 0.2 mM NaCO3, 3.0 mM NaHCO3, 0.09 mM K2HPO4, 3 μM chelated iron, 0.4 μM NaMoO4, 1.0 μM MnCl2, 0.008 μM CoCl2, 0.35 μM ZnSO4, 0.08 μM CuSO4, phenol red 0.5% w/v] amended with hydroxylamine at 100 μM final concentration (AMS-H) in a 250 ml Erlenmeyer flask (Pyrex, USA) and maintained as standing batch cultures at 30°C, in the dark, without shaking (Campbell et al., 2011). Phenol-red was used as a pH indicator and pH was adjusted daily to 7.5∼8.0 using 0.25 M K2CO3. Batch cultures were propagated monthly by transferring 5 mL of culture into 100 mL of fresh AMS-H medium using a sterile pipette (Corning, USA). Hydroxylamine was added weekly to a final concentration of 200 μM to eliminate susceptible bacteria and archaea (Campbell et al., 2011). Nitrite (NO2-) was measured in spent medium samples (1.3 mL) as an indicator of catabolic activity of nitrifying bacteria by the sequential addition of 0.5 mL of 1% sulfanilamide and 0.5 mL of 0.02% N-(l-napthyl) ethylenediamine dihydrochloride using sterile medium as a control (Nicholas and Nason, 1957; Koops et al., 1990). Following incubation for 20 min at room temperature, absorbance at 543 nm was measured using a Nanodrop 200C spectrophotometer (Thermo scientific Inc.). Contamination with heterotrophs potentially introduced during handling and transfer of cultures was assessed by plating an aliquot on LB agar and incubation at 35°C.
Biomass was harvested from 200 mL of D1 enrichment in exponential phase (based on NO2- measurement) by centrifugation at 4,300 × g for 15 min using a Sorvall Evolution centrifuge and SS-34 rotor (Thermo Fisher Scientific, USA). The resulting cell pellet was resuspended in water to a volume of 200 μl for metagenomic DNA extraction using the MP FastDNA Spin Kit (MP Biomedicals, USA) and a MP FastPrep®-24 Instrument (MP Biomedicals, USA; Poret-Peterson et al., 2008). Genomic DNA, primers targeting universal and gammaproteobacteria-16S rRNA-specific sequences of the 16S rRNA gene (Campbell et al., 2011) as well as primers specific for detection of the haoA gene from gamma-AOB (Schmid et al., 2008) and the polymerase chain reaction (PCR) were used to generate inserts for clone libraries (Table 3). Amplicons (universal 16S ∼1470 bp, gamma-16S ∼1500 bp, and gamma-haoA ∼1220 bp) were generated using GoTaq Green Master® Mix (Promega, USA) on a gradient mastercycler (Eppendorf, Germany) following an optimized PCR protocol (30 cycles of 1 min at 95°C, 1 min at annealing temperature (Tm): 55°C for universal 16S, 56°C for gamma 16S and 60°C for gamma haoA, 1 min extension time at 72°C). For each 50-μl PCR reaction, 2 μl of extracted D1 metagenomic DNA was used as the template. The PCR using gamma-16S rRNA-gene primers generated products of the expected length; therefore, PCR products obtained with universal 16S rRNA and gamma-haoA gene primers were gel-purified (1.5% low-melt agarose) and ligated into pCR®2.1-TOPO TA vectors (Invitrogen, USA). The overnight ligation products were used to transform TOP10 competent cells for clone library construction. Forty randomly picked recombinants from each clone library were selected using X-Gal-LB plates with 25 μg/ml kanamycin (Dang et al., 2008, 2009, 2011; Campbell et al., 2011). Plasmids were extracted using the alkaline mini-prep method (Wizard Plus SV Minipreps DNA Purification Systems; Promega, USA) and inserts were sequenced with vector primers, M13F and M13R (Table 3; DNA Core Facility, University of Louisville). Insert sequences were trimmed, aligned using MUSCLE (Edgar, 2004) provided at the EMBL-EBI webserver1 and analyzed for closest-matches with sequences deposited in GenBank using the BLASTN program (Altschul et al., 1997).
Isolation of D1FHS into Pure Culture, DNA Extraction, Genome Sequencing, and Annotation
Sample aliquots of 0.1 ml were sequentially diluted, plated on sterile semi-solid (0.3% agar) AMS-H medium and incubated in the dark at 30°C (Tittsler and Sandholzer, 1936). After 2 months, visible single colonies were picked and transferred into a 6-well plate with 5 ml liquid AMS-H media. Subsamples (2.5 ml) of wells that turned yellow, indicating acidification due to nitrous acid production, were transferred into 250 ml Erlenmeyer flasks with 100 ml AMS-H media. Enriched culture (20 mL) was used to inoculate a benchtop bioreactor (BioFlo®/CelliGen® 115, New Brunswick/Eppendorf, Germany) with 1 L working volume of AMS-H media sterilized by autoclaving (121°C, 20 min, 100 Kpa, 15 psi). The chemostat, operated in batch culture mode, was incubated in the dark at 28°C with agitation of 50 rpm and 1.0 LPM flow of 0.22 μm filtered air; pH was maintained at 8.0 by automated addition of sterile 0.5 M K2CO3. Cell density (OD600) and nitrite concentration were measured daily by withdrawing a 5-ml sample from the reactor. Continuous culture for 7–10 days was alternated with batch culture (500 mL in 2-L Erlenmeyer flask) several times; each time an aliquot of the extracted genomic DNA was used to build a universal 16S rRNA cloning library for insert sequencing and analysis of sequence from 40 randomly picked clones (see text above) after each alternation. The alignments were inspected for single nucleotide polymorphisms. Once a single 16S rRNA sequence was observed across all clones in multiple subcultures, the culture was deemed pure and labeled “D1FHS.”
To inspect morphological homogeneity of the D1FHS culture, cells from 50 ml of a 10-day old culture were harvested by centrifugation (4,300 × g, 15 min). The cell pellet was washed twice with AMS medium and resuspended in 10 μl of AMS medium. A 2-μl aliquot of cell suspension was spread on a glass slide and covered with a cover slip to observe cell motility with a light microscope (Olympus, Japan) at 400× magnification. An 8-μl aliquot of cell suspension was heat-fixed on a slide for Gram staining (Sigma–Aldrich, USA) and observed using a light microscope (Olympus, Japan) at 1000× magnification.
For genome sequencing, ∼1.2 μg of gDNA were isolated from a 500-mL mid-exponential-phase culture. To obtain sufficient DNA for Pacific Biosciences (PacBio) single molecule sequencing, the purified genomic DNA was amplified by multiple displacement amplification (MDA) based on published protocols (Zhang et al., 2006). Briefly, 100 ng amounts of DNA were amplified in five MDA reactions (REPLI-g Mini, Qiagen). MDA products from each reaction (∼4 μg per reaction) were sequentially treated with 10 U Phi29 (Fermentas) and 100 U S1 nuclease (Thermo Fisher) followed by precipitation (SureClean, Bioline) and resuspension in sterile 18 Mω water. The identity of the amplified DNA was verified by PCR amplification with universal 16S rRNA primers. Half of the PCR product was digested with PstI (Fermentas) to verify that they contained a single site as predicted and the remainder was directly sequenced with the amplification primers to verify that the sequence matched that from D1FHS. Following confirmation, the amplified DNA (∼15 μg total) was submitted for PacBio sequencing at the University of Delaware Genotyping and Sequencing Core Facility, University of Delaware. Sequence reads were assembled using RS_HGAP_Assembly protocol v.3 on a computational cluster maintained by the University of Delaware Center for Bioinformatics and Computational Biology. Completeness of the assembled genome was assessed by BLASTN and TBLASTN searches (CLC Main Workbench, Qiagen) using the D1FHS 16S rRNA sequence and a collection of 31 single-copy genes (Wu and Eisen, 2008). Biosample (SAMN04324281) and Bioproject (PRJNA305330) information has been submitted to the NCBI.
HaoA Protein Sequence Analysis
HaoA protein sequences from Nitrosococcus species, including N. oceani strains C-107, C-27, AFC132, and AFC27, N. watsonii C-113, and N. halophilus Nc4, the deduced HaoA sequence from D1FHS and pertinent sequences of the “HAO” clade of octaheme cytochrome c proteins were aligned using MUSCLE and manually refined by comparison with previously published results from phylogenetic analyses including structural and protein sequence-analytical features (Klotz et al., 2008; Kern et al., 2011). N- and C-terminally extending sequences beyond the first and last heme-binding motif (CxxCH), respectively, were trimmed and the final alignment was subjected to a Bayesian inference of phylogeny using the BEAST package (v1.8.1 of BEAUti, BEAST and TreeAnnotator; FigTree v.1.4; Drummond et al., 2012). By utilizing unique sites, tree likelihoods (ignoring ambiguities) were determined for the alignment by creating a Monte–Carlo Markov Chain (10,000,000 generations) in three independent runs. The searches were conducted assuming an equal distribution of rates across sites, sampling every 1000th generation and using the WAG empirical amino acid substitution model (Whelan and Goldman, 2001). The resulting 10,000 trees (omitting the first 350 trees as burn-in) were used to construct a phylogenetic consensus tree.
Average Nucleotide Identity between Bacterial Genomes
The Average Nucleotide Identity (ANI) calculated for pairwise compared genomes is an advanced method developed to identify species boundaries for bacteria and archaea (Konstantinidis and Tiedje, 2005). The species delineation cut-off point is an ANI of 94% (Konstantinidis and Tiedje, 2005; Goris et al., 2007), which corresponds to the traditional species cut-off of 70% using DNA–DNA hybridization (Wayne et al., 1987). The high quality draft genome sequence of strain D1FHS (single chromosome, no plasmid) was compared with the genome sequences of N. halophilus Nc4 (CP001798; plasmid: CP001799), N. oceani ATCC19707 (CP000127; plasmid: CP000126), N. watsonii C-113 [CP002086; plasmids: CP002087 & CP002088] as well as those of other Chromatiaceae, including N. oceani strain AFC-27 [ABSG00000000], N. oceani strain AFC132 [JPFN00000000], N. oceani strain C-27 [JPGN00000000], Allochromatium vinosum strain DSM 180 [CP001896; plasmids: CP001897 & CP001898], Halorhodospira halophila strain SL1 [CP000544], Alkalimnicola ehrlichei strain MLHE-1 [CP000453], Thioalkalivibrio sp. strains HL-EbGR7 [CP001339] and K90mix [CP001905; plasmid: CP00196] and Halothiobacillus neopolitanus strain c2 [CP001801] as previously done for the identification of N. watsonii as a species (Campbell et al., 2011). The ANI values reported in this paper were calculated using the tool developed and provided online by Luis M. Rodriguez-R and Kostas T. Konstantinidis2 using default parameters.
Physiological Assays and Phenotypic Assessment
One-milliliter aliquots of mid-exponential-phase “D1FHS” pure culture were inoculated into 25 ml HEPES-buffered AMS-H media using 125 ml flasks and incubated in the dark at 30°C without shaking to determine growth as a function of time, salt and ammonium tolerance as well as optimum pH and temperature. A growth curve was estimated by measuring nitrite concentrations at eight times (days 0, 1, 3, 5, 6, 8, 10, and 12) after inoculation in AMS-H standard media [700 mM NaCl and 12.5 mM (NH4)2SO4 at a set pH of 8]. AMS-H media at eight different NaCl concentrations ranging from 100 to 1600 mM were used to determine the optimum NaCl concentration for the growth of D1FHS. Growth of the D1FHS culture was also tested at eight pH conditions (pH of 4, 5, 6,7, 7.5, 8, 9, and 10) and at five temperatures (4, 20, 28, 37, and 45C). Ammonium tolerance was examined by testing 10 ammonium concentrations [using (NH4)2 SO4] ranging from 0 to 600 mM. Cell cultures were monitored for nitrite production routinely. All tests were performed with three technical replicates and repeated with two or more biological replicates. Like N. halophilus Nc4 and unlike cultures of N. oceani C-107 and N. watsonii C-113 (Campbell et al., 2011), D1FHS did not grow on urea when provided as the sole source of ammonium (N, energy and reductant). Flagellar motility was observed under a light microscope at 1000× magnification.
Results and Discussion
Characterization of the Enrichment and Pure Isolate Cultures
PCR performed with metagenomic DNA extracted from D1 enrichment culture targeting at universal 16S rRNA, gamma 16S rRNA and gamma-haoA genes, yielded correctly sized fragments (universal 16S ∼1470 bp, gamma-16S ∼1500 bp, and gamma-haoA ∼1220 bp, respectively), indicating that the D1 enrichment culture contained gamma-AOB. Universal 16S rRNA gene sequences were analyzed based on BLASTN searches of the NCBI database; 29 of the 40 sequences were best hits to 16S rRNA gene sequences from Nitrosococcus species: 87–98% identity to N. halophilus Nc4; 93–98% identity to N. oceani C-107, and 93–98% identity to N. watsonii C-113. The rest of the sequences had best hits to sequences from Mesorhizobium, Janibacter sp., or uncultured bacteria.
Following the repeated enrichment procedure described in the experimental procedures that involved selection of colonies from semi-solid agar, alternating batch and continuous culture and the application of functional pressure by exposure to hydroxylamine and high sodium salt concentration, isolated genomic DNA from the final culture labeled “D1FHS” was, again, investigated using primers designed to detect universal 16S rRNA, gamma 16S rRNA and gamma-haoA genes and the PCR. Analysis of the obtained amplicon sequences using BLASTN consistently revealed highest identity with sequences from N. halophilus Nc4 and no clone had a best hit with sequences from bacteria outside of the Nitrosococcus genus.
Genome Sequencing and Characteristics
PacBio-based high-throughput sequencing of D1FHS DNA produced 2,267,036,463 bp of sequence in 191,850 reads (N50 read length = 16,723 bp, mean read length = 11,816 bp). These were assembled into a high quality draft sequence containing one contig of 4,022,640 base pairs; the genome has a GC content of 50.7%. Homologs were detected for each of 31 universal single copy genes (Wu and Eisen, 2008) by TBLASTN searches and two copies of the D1FHS 16S rRNA sequence by BLASTN indicating that the assembled contig contains a complete bacterial genome. The sequence of the D1FHS genome was analyzed using the RAST annotation server (Aziz et al., 2008) and found to contain 4,128 protein-coding DNA sequences, 51 RNA genes including two 16S-23S-5S rRNA operons (RAST-ID: 6666666.126356). The genome sequence of D1FHS included genetic markers known from other Nitrosococcus genomes: the genome of D1FHS encodes a complete set of protein inventories required for ammonia-dependent chemolithotrophic growth, including ammonia monooxygenase (EC 220.127.116.11), hydroxylamine dehydrogenase (EC 18.104.22.168), cytochromes c554 and cM552 as well as nitrosocyanin, although the functions of cytochromes c554 and nitrosocyanin are still elusive (Stein et al., 2013). Based on the identification of signal peptides, all of these proteins operate and have access to their reactants within the periplasm or at the periplasmic side of the plasma membrane (Arp et al., 2007). Consistent with the inability of D1FHS to grow on urea, the genome did not encode urea hydrolase.
Average Nucleotide Identity and HaoA Protein Phylogenetic Analysis
The ANI values for the comparison of the D1FHS genome with the genomes from Proteobacteria outside the genus Nitrosococcus were low as previously reported for the analysis of the N. watsonii C-113 genome (Campbell et al., 2011). ANI values calculated between the D1FHS genome and all available Nitrosococcus genomes including N. halophilus strain Nc4 (89.5%), N. watsonii C-113 (78.2%) and N. oceani (C-107, 78.4%; C-27, 78.2%; AFC27, 78.4%; AFC132, 78.2%) are significantly below the cut-off for the delineation of species (Konstantinidis et al., 2006), which suggests that strain D1FHS is representative of a species distinct from N. halophilus, N. watsonii, and N. oceani most closely related to N. halophilus. We propose to assign D1FHS as the type strain of N. wardiae spec. nov.
Phylogenetic analysis of protein sequences confirmed that the HaoA proteins from D1FHS and N. halophilus strain Nc4 are more closely related to one another than either is to HaoA proteins from strains of N. watsonii or N. oceani (Figure 2). Identical results were obtained from phylogenetic analysis of the other proteins implicated in ammonia-dependent chemolithotrophy (Arp et al., 2007; Stein et al., 2013): ammonia monooxygenase, cytochrome cM552 and nitrosocyanin (data not shown). The phylogram in Figure 2 supports the conclusion that the ancestor of the Nitrosococcus genus diverged into two phylotypes that gave rise to N. wardiae spec. nov. and N. halophilus as well as N. watsonii and N. oceani. It was previously hypothesized that the hao genes encoded on the chromosome (soil methane-oxidizing) or a plasmid (marine sulfur-oxidizing) of Alphaproteobacteria and the hao genes in the genomes of betaproteobacterial AOB were obtained by horizontal transfer from Gammaproteobacteria (Klotz et al., 2008; Kern et al., 2011; Stein et al., 2011). The position of the Hao protein sequence from the “comammox” bacterium Nitrospira inopinata ERN4 capable of oxidizing ammonia all the way to nitrate (Daims et al., 2015) included in this analysis suggests that the hao gene cluster has been horizontally transferred from Gammaproteobacteria into the genomes of nitrite-oxidizing Nitrospira bacteria ancestral to the comammox lineage as well. Preliminary analysis of the D1FHS genome suggests that its genome is likely most representative of the ancestral Nitrosococcus genome; a detailed analysis of the genome is forthcoming.
FIGURE 2. Consensus tree constructed after Bayesian inference of phylogeny from the MUSCLE alignments of sequences from octaheme cytochrome c proteins in the “Hao cluster” (Klotz et al., 2008; Kern et al., 2011) using the sequence from “reverse Hao” proteins as the out group (Hanson et al., 2013). The N- and C-terminal sequences outside the heme c binding motifs were eliminated from the final alignment that was subjected to Bayesian inference of phylogeny using the BEAST package (see Experimental procedures). Posterior probability values of all nodes were ≥0.99. Mean branch lengths are characterized by a scale bar indicating the evolutionary distance between the proteins (changes per amino acid position). The branches are annotated with labels indicating the source organism and the protein sequence accession number.
Growth-Physiological Characterization of D1FHS (Table 4)
TABLE 4. Classification and general features of Nitrosococcus wardiae D1FHS according to MIGS Recommendations (Chain et al., 2009).
In batch culture, cells of D1FHS routinely exhibited a lag time of 5 days before entering exponential growth phase. Concomitantly, cells started dividing and Gram-staining revealed that cells appear as short rods, arranged singly or in pairs but not in tetrads (Figure 3). The stain also confirmed the existence of only one cellular morphotype. In agreement with strains of other Nitrosococcus species, D1FHS contained intracellular membrane stacks, was lophotrichous and capable of flagellar motility (data not shown).
FIGURE 3. Gram Stain, 1000× magnification using a light microscope. Cell shape is short rods, arranged singly or in pairs.
The nitrite production rate was used to assess growth physiological characteristics of D1FHS (“optimum catabolism”), which was calculated based on a nitrite standard curve (calculated as μmol per L of media). The pH, temperature, NaCl and ammonium tolerances were determined by measuring the maximum nitrite production of log-phase cultures (day 5 to day 8; Figure 4A). Based on the observed maximum nitrite production rate, the optimum growth temperature for D1FHS is 37°C (Figure 4B), which is much higher than those reported (also based on nitrite production) for other strains of the genus Nitrosococcus (Table 5; Koops et al., 1990; Purkhold et al., 2000; Campbell et al., 2011). There is presently no published information for any Nitrosococcus strain about the need for maintenance energy at different temperatures. No growth (no changes in nitrite levels) was observed at 4 and 52°C, but D1FHS tolerated 45°C (Figure 4B). D1FHS was able to grow at a pH ranging from 5 to 9 with an optimum pH between 7.5 and 8.0 (Figure 4C), which is in the same range reported for other strains of the genus Nitrosococcus (Table 5). D1FHS was not able to tolerate ammonium concentrations higher than 300 mM, which represents a much lower tolerance than what is known for other Nitrosococcus strains (Table 5, Figure 4D; Koops et al., 1990; Purkhold et al., 2000; Campbell et al., 2011); consequently, the optimum ammonium concentration of 50 mM is also the lowest known optimum for all characterized strains of Nitrosococcus. D1FHS exhibited a NaCl tolerance similar to that of N. halophilus as it can grow at salt concentrations up to 1600 mM; the optimum salt concentration is 700 mM (Figure 4E). Taken together, Nitrosococcus wardiae D1FHS represents a new species that is characterized by a distinct growth physiology from other Nitrosococcus species.
FIGURE 4. (A) Growth of D1FHS (assessed by maximum nitrite production rate) in AMS-H standard medium with 700 mM NaCl and 12.5 mM (NH4)2SO4 at pH 8, grown at 28°C in the dark. (B,C) Physiological tests: growth of D1FHS in AMS-H in the dark (assessed by maximum nitrite production rate) was characterized as a function of (B) temperature and (C) pH conditions. Data points are mean values calculated from three biological replicates. (D) Physiological tests: growth of D1FHS in AMS-H in the dark (assessed by maximum nitrite production rate) was characterized as a function of ammonium concentration. Data points are mean values calculated from three biological replicates. (E) Physiological tests: growth of D1FHS in AMS-H in the dark (assessed by maximum nitrite production rate) was characterized as a function of NaCl concentration. Data points are mean values calculated from three biological replicates.
Emended Description of the Genus Nitrosococcus
Nitrosococcus (GenBank Taxonomy ID: 1227) synonym Nitrosococcus Winogradsky (Winogradsky, 1892); ex “Micrococcus nitrosus” Migula (Migula, 1900); ex “Nitrosococcus nitrosus” Buchanan (Buchanan, 1925; Editorial_Board, 1955; Commission, 1958). Name approved by Skerman et al. (1980), but the culture had been lost decades before. The etymology of Nitrosococcus is cataloged using digital optical identifiers as follows: Genus3, Family Chromatiaceae4 (Bavendamm, 1924), Order Chromatiales5 (Imhoff, 2005; List_Editor, 2005), Class Gammaproteobacteria6, phylum Proteobacteria7 (Garrity and Holt, 2001).
Description of Nitrosococcus wardiae spec. nov.
Nitrosus (Latin masculine adjective): nitrous; coccus: (Latin masculine adjective): sphere; wardiae (ward.i’ae. N.L. fem. gen. n. wardiae): of Ward, named after the American microbiologist Bess B. Ward for her pioneering work on the marine Nitrogen cycle including the study of AOB in the genus Nitrosococcus.
Locality as Well as Culture History
Collected on October 26, 2008, at station D1 in Jiaozhou Bay, a marginal sea bay in the Yellow Sea of China (36.067°N, 120.230°E) by Hongyue Dang and Martin G. Klotz. Maintained as enrichment culture D1 in the Klotz Lab (University of Louisville). Identified as a gammaproteobacterial AOB in an enrichment-culture in 2011 by Lin Wang and Martin G. Klotz and later purified using a hydroxylamine treatment regime by Lin Wang, Chee Kent Lim, and Martin G. Klotz (at UNC Charlotte).
Appearing as large cocci or very short rods. Cells contain a well-developed intra-cellular membrane system of an arrangement that appears as one stack of membrane vesicles packed mainly in the center of the cell. Polar flagella allow for motility, the genotype is chemotaxis-positive. Light sensitive; Strictly aerobic; Moderately alkaliphilic and mesophilic with an optimum growth temperature at 37°C; Grows in the presence of sodium salts between 200 to 800 mM with an optimum at 700 mM; Grows at ammonium concentrations below 300 mM with an optimum at 50 mM; The size of the genome is 4,022,640 base pairs with a GC content of 50.7%. No plasmid was identified. The genome of N. wardiae strain D1FHST (GenBank Taxonomy ID: 1814290) is presently being annotated and analyzed. MIGS information (Chain et al., 2009) to correlate the genotype with the environment of its isolation is provided in Table 4.
MGK conceived and designed the experiments; LW, CKL, HD, and TEH performed the experiments; LW, TEH, and MGK analyzed the data; and the manuscript was written by MGK with input from all authors.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Pertinent taxonomic information was accessed through the “Names for Life” online tool (http://www.namesforlife.com). Nitrogen cycle research in the Klotz lab has been supported during the duration of this project by the NSF (EF-0412129, EF-0541797, MCD1202648) the Gordon and Betty Moore Foundation, the China National Science Foundation (Grants 91328209 and 41076091), a distinguished fellowship at State Key Laboratory of Marine Environmental Science at Xiamen University, China, and institutional incentive funds from the University of Louisville (Office of the EVP for Research). LW and CKL were supported by incentive funds from UNC Charlotte (Office of the VC for Research and the College of Liberal Arts & Sciences). HD was supported by China MOST 973 grant 2013CB955700, NSFC grants 91328209, 41076091 and 91428308, and SOA grant GASI-03-01-02-05.
Altschul, S. F., Madden, T. L., Schaffer, A. A., Zhang, J., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402. doi: 10.1093/nar/25.17.3389
Arp, D. J., Chain, P. S., and Klotz, M. G. (2007). The impact of genome analyses on our understanding of ammonia-oxidizing bacteria. Annu. Rev. Microbiol. 61, 503–528. doi: 10.1146/annurev.micro.61.080706.093449
Aziz, R. K., Bartels, D., Best, A. A., DeJongh, M., Disz, T., Edwards, R. A., et al. (2008). The RAST server: rapid annotations using subsystems technology. BMC Genomics 9:75. doi: 10.1186/1471-2164-9-75
Campbell, M. A., Chain, P. S., Dang, H., El Sheikh, A. F., Norton, J. M., Ward, N. L., et al. (2011). Nitrosococcus watsonii sp. nov., a new species of marine obligate ammonia-oxidizing bacteria that is not omnipresent in the world’s oceans: calls to validate the names ‘Nitrosococcus halophilus’ and ‘Nitrosomonas mobilis’. FEMS Microbiol. Ecol. 76, 39–48. doi: 10.1111/j.1574-6941.2010.01027.x
Chain, P. S. G., Grafham, D. V., Fulton, R. S., FitzGerald, M. G., Hostetler, J., Muzny, D., et al. (2009). Genomics. Genome project standards in a new era of sequencing. Science 326, 236–237. doi: 10.1126/science.1180614
Commission (1958). Opinion 23. Rejection of the generic names Nitromonas Winogradsky 1890 and Nitromonas Orla-Jensen 1909, Conservation of the generic Names Nitrosomonas Winogradsky, 1892, Nitrosococcus Winogradsky, 1892, and the designation of the type species of these genera. Int. Bull. Bacteriol. Nomencl. Taxon 8, 169–170.
Dang, H., Chen, R., Wang, L., Guo, L., Chen, P., Tang, Z., et al. (2010a). Environmental factors shape sediment anammox bacterial communities in hypernutrified Jiaozhou Bay, China. Appl. Environ. Microbiol. 76, 7036–7047. doi: 10.1128/AEM.01264-10
Dang, H., Chen, R., Wang, L., Shao, S., Dai, L., Ye, Y., et al. (2011). Molecular characterization of putative biocorroding microbiota with a novel niche detection of Epsilon- and Zetaproteobacteria in Pacific Ocean coastal seawaters. Environ. Microbiol. 13, 3059–3074. doi: 10.1111/j.1462-2920.2011.02583.x
Dang, H., Luan, X. W., Chen, R., Zhang, X., Guo, L., and Klotz, M. G. (2010b). Diversity, abundance and distribution of amoA-encoding archaea in deep-sea methane seep sediments of the Okhotsk Sea. FEMS Microbiol. Ecol. 72, 370–385. doi: 10.1111/j.1574-6941.2010.00870.x
Dang, H., Wang, C., Li, J., Li, T., Tian, F., Jin, W., et al. (2009). Diversity and distribution of sediment NirS -encoding bacterial assemblages in response to environmental gradients in the eutrophied Jiaozhou Bay, China. Microb. Ecol. 58, 161–169. doi: 10.1007/s00248-008-9469-5
Dang, H., Zhang, X., Sun, J., Li, T., Zhang, Z., and Yang, G. (2008). Diversity and spatial distribution of sediment ammonia-oxidizing crenarchaeota in response to estuarine and environmental gradients in the Changjiang Estuary and East China Sea. Microbiology 154, 2084–2095. doi: 10.1099/mic.0.2007/013581-0
Garrity, G. M., and Holt, J. G. (2001). “The road map to the manual,” in Bergey’s Manual of Systematic Bacteriology, eds G. M. Garrity, D. R. Boone, and R. W. Castenholz (New York, NY: Springer), 119–169.
Goris, J., Konstantinidis, K. T., Klappenbach, J. A., Coenye, T., Vandamme, P., and Tiedje, J. M. (2007). DNA-DNA hybridization values and their relationship to whole-genome sequence similarities. Int. J. Syst. Evol. Microbiol. 57, 81–91. doi: 10.1099/ijs.0.64483-0
Head, I. M., Hiorns, W. D., Embley, T. M., McCarthy, A. J., and Saunders, J. R. (1993). The phylogeny of autotrophic ammonia-oxidizing bacteria as determined by analysis of 16S ribosomal RNA gene sequences. J. Gen. Microbiol. 139, 1147–1153. doi: 10.1099/00221287-139-6-1147
Kern, M., Klotz, M. G., and Simon, J. (2011). The Wolinella succinogenes mcc gene cluster encodes an unconventional respiratory sulphite reduction system. Mol. Microbiol. 82, 1515–1530. doi: 10.1111/j.1365-2958.2011.07906.x
Klotz, M. G., Schmid, M. C., Strous, M., op den Camp, H. J., Jetten, M. S., and Hooper, A. B. (2008). Evolution of an octahaem cytochrome c protein family that is key to aerobic and anaerobic ammonia oxidation by bacteria. Environ. Microbiol. 10, 3150–3163. doi: 10.1111/j.1462-2920.2008.01733.x
List_Editor (2005). Validation of publication of new names and new combinations previously effectively published outside the IJSEM. List no. 106. Int. J. Syst. Evol. Microbiol. 55, 2235–2238. doi: 10.1099/ijs.0.64108-0
Poret-Peterson, A. T., Graham, J. E., Gulledge, J., and Klotz, M. G. (2008). Transcription of nitrification genes by the methane-oxidizing bacterium, Methylococcus capsulatus strain Bath. ISME J. 2, 1213–1220. doi: 10.1038/ismej.2008.71
Purkhold, U., Pommerening-Roser, A., Juretschko, S., Schmid, M. C., Koops, H.-P., and Wagner, M. (2000). Phylogeny of all recognized species of ammonia oxidizers based on comparative 16S rRNA and amoA sequence analysis: implications for molecular diversity surveys. Appl. Environ. Microbiol. 66, 5368–5382. doi: 10.1128/AEM.66.12.5368-5382.2000
Schmid, M. C., Hooper, A. B., Klotz, M. G., Woebken, D., Lam, P., Kuypers, M. M., et al. (2008). Environmental detection of octahaem cytochrome c hydroxylamine/hydrazine oxidoreductase genes of aerobic and anaerobic ammonium-oxidizing bacteria. Environ. Microbiol. 10, 3140–3149. doi: 10.1111/j.1462-2920.2008.01732.x
Stein, L. Y., Bringel, F., DiSpirito, A. A., Han, S., Jetten, M. S., Kalyuzhnaya, M. G., et al. (2011). Genome sequence of the methanotrophic alphaproteobacterium Methylocystis sp. strain Rockwell (ATCC 49242). J. Bacteriol. 193, 2668–2669. doi: 10.1128/JB.00278-11
Stein, L. Y., Campbell, M. A., and Klotz, M. G. (2013). Energy-mediated vs. ammonium-regulated gene expression in the obligate ammonia-oxidizing bacterium, Nitrosococcus oceani. Front. Microbiol. 4:277. doi: 10.3389/fmicb.2013.00277
Turner, S., Pryer, K. M., Miao, V. P. W., and Palmer, J. D. (1999). Investigating deep phylogenetic relationships among cyanobacteria and plastids by small subunit rRNA sequence analysis. J. Eukaryot. Microbiol. 46, 327–338. doi: 10.1111/j.1550-7408.1999.tb04612.x
Urakawa, H., Garcia, J. C., Nielsen, J. L., Le, V. Q., Kozlowski, J. A., Stein, L. Y., et al. (2015). Nitrosospira lacus sp. nov., a psychrotolerant, ammonia-oxidizing bacterium from sandy lake sediment. Int. J. Syst. Evol. Microbiol. 65, 242–250. doi: 10.1099/ijs.0.070789-0
Ward, B. B., and O’Mullan, G. D. (2002). Worldwide distribution of Nitrosococcus oceani, a marine ammonia-oxidizing gamma-proteobacterium, detected by PCR and sequencing of 16S rRNA and amoA genes. Appl. Environ. Microbiol. 68, 4153–4157. doi: 10.1128/AEM.68.8.4153-4157.2002
Wayne, L. G., Brenner, D. J., Colwell, R. R., Grimont, P. A. D., Kandler, O., Krichevsky, M. I., et al. (1987). Report of the ad hoc committee on reconciliation of approaches to bacterial systematics. Int. J. Syst. Bacteriol. 37, 463–464. doi: 10.1099/00207713-37-4-463
Whelan, S., and Goldman, N. (2001). A general empirical model of protein evolution derived from multiple protein families using a maximum-likelihood approach. Mol. Biol. Evol. 18, 691–699. doi: 10.1093/oxfordjournals.molbev.a003851
Keywords: ammonia-oxidizing bacteria, enrichment, Jiaozhou Bay, Nitrosococcus wardiae
Citation: Wang L, Lim CK, Dang H, Hanson TE and Klotz MG (2016) D1FHS, the Type Strain of the Ammonia-Oxidizing Bacterium Nitrosococcus wardiae spec. nov.: Enrichment, Isolation, Phylogenetic, and Growth Physiological Characterization. Front. Microbiol. 7:512. doi: 10.3389/fmicb.2016.00512
Received: 31 December 2015; Accepted: 29 March 2016;
Published: 14 April 2016.
Edited by:Ludmila Chistoserdova, University of Washington, USA
Reviewed by:Sergey M. Stolyar, University of Idaho, USA
Jeanette M. Norton, Utah State University, USA
Copyright © 2016 Wang, Lim, Dang, Hanson and Klotz. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Martin G. Klotz, firstname.lastname@example.org