Single cell analyses reveal contrasting life strategies of the two main nitrifiers in the ocean.

Nitrification, the oxidation of ammonia via nitrite to nitrate, is a key process in marine nitrogen (N) cycling. Although oceanic ammonia and nitrite oxidation are balanced, ammonia-oxidizing archaea (AOA) vastly outnumber the main nitrite oxidizers, the bacterial Nitrospinae. The ecophysiological reasons for this discrepancy in abundance are unclear. Here, we compare substrate utilization and growth of Nitrospinae to AOA in the Gulf of Mexico. Based on our results, more than half of the Nitrospinae cellular N-demand is met by the organic-N compounds urea and cyanate, while AOA mainly assimilate ammonium. Nitrospinae have, under in situ conditions, around four-times higher biomass yield and five-times higher growth rates than AOA, despite their ten-fold lower abundance. Our combined results indicate that differences in mortality between Nitrospinae and AOA, rather than thermodynamics, biomass yield and cell size, determine the abundances of these main marine nitrifiers. Furthermore, there is no need to invoke yet undiscovered, abundant nitrite oxidizers to explain nitrification rates in the ocean.

Phylogenetic and structural diversity of aromatically dense pili from environmental metagenomes.

Electroactive type IV pili, or e-pili, are used by some microbial species for extracellular electron transfer. Recent studies suggest that e-pili may be more phylogenetically and structurally diverse than previously assumed. Here, we used updated aromatic density thresholds (≥9.8% aromatic amino acids, ≤22-aa aromatic gaps and aromatic amino acids at residues 1, 24, 27, 50 and/or 51, and 32 and/or 57) to search for putative e-pilin genes in metagenomes from diverse ecosystems with active microbial metal cycling. Environmental putative e-pilins were diverse in length and phylogeny, and included truncated e-pilins in Geobacter spp., as well as longer putative e-pilins in Fe(II)-oxidizing Betaproteobacteria and Zetaproteobacteria.

Cyanate and urea are substrates for nitrification by Thaumarchaeota in the marine environment.

Ammonia-oxidizing archaea of the phylum Thaumarchaeota are among the most abundant marine microorganisms1. These organisms thrive in the oceans despite ammonium being present at low nanomolar concentrations2,3. Some Thaumarchaeota isolates have been shown to utilize urea and cyanate as energy and N sources through intracellular conversion to ammonium4-6. Yet, it is unclear whether patterns observed in culture extend to marine Thaumarchaeota, and whether Thaumarchaeota in the ocean directly utilize urea and cyanate or rely on co-occurring microorganisms to break these substrates down to ammonium. Urea utilization has been reported for marine ammonia-oxidizing communities7-10, but no evidence of cyanate utilization exists for marine ammonia oxidizers. Here, we demonstrate that in the Gulf of Mexico, Thaumarchaeota use urea and cyanate both directly and indirectly as energy and N sources. We observed substantial and linear rates of nitrite production from urea and cyanate additions, which often persisted even when ammonium was added to micromolar concentrations. Furthermore, single-cell analysis revealed that the Thaumarchaeota incorporated ammonium-, urea- and cyanate-derived N at significantly higher rates than most other microorganisms. Yet, no cyanases were detected in thaumarchaeal genomic data from the Gulf of Mexico. Therefore, we tested cyanate utilization in Nitrosopumilus maritimus, which also lacks a canonical cyanase, and showed that cyanate was oxidized to nitrite. Our findings demonstrate that marine Thaumarchaeota can use urea and cyanate as both an energy and N source. On the basis of these results, we hypothesize that urea and cyanate are substrates for ammonia-oxidizing Thaumarchaeota throughout the ocean.

Broad Phylogenetic Diversity Associated with Nitrogen Loss through Sulfur Oxidation in a Large Public Marine Aquarium.

Denitrification by sulfur-oxidizing bacteria is an effective nitrate removal strategy in engineered aquatic systems. However, the community taxonomic and metabolic diversity of sulfur-driven denitrification (SDN) systems, as well as the relationship between nitrate removal and SDN community structure, remains underexplored. This is particularly true for SDN reactors applied to marine aquaria, despite the increasing use of this technology to supplement filtration. We applied 16S rRNA gene, metagenomic, and metatranscriptomic analyses to explore the microbial basis of SDN reactors operating on Georgia Aquarium's Ocean Voyager, the largest indoor closed-system seawater exhibit in the United States. The exhibit's two SDN systems vary in water retention time and nitrate removal efficiency. The systems also support significantly different microbial communities. These communities contain canonical SDN bacteria, including a strain related to Thiobacillus thioparus that dominates the system with the higher water retention time and nitrate removal but is effectively absent from the other system. Both systems contain a wide diversity of other microbes whose metagenome-assembled genomes contain genes of SDN metabolism. These include hundreds of strains of the epsilonproteobacterium Sulfurimonas, as well as gammaproteobacterial sulfur oxidizers of the Thiotrichales and Chromatiales, and a relative of Sedimenticolathiotaurini with complete denitrification potential. The SDN genes are transcribed and the taxonomic richness of the transcript pool varies markedly among the enzymatic steps, with some steps dominated by transcripts from noncanonical SDN taxa. These results indicate complex and variable SDN communities that may involve chemical dependencies among taxa as well as the potential for altering community structure to optimize nitrate removal.IMPORTANCE Engineered aquatic systems such as aquaria and aquaculture facilities have large societal value. Ensuring the health of animals in these systems requires understanding how microorganisms contribute to chemical cycling and waste removal. Focusing on the largest seawater aquarium in the United States, we explore the microbial communities in specialized reactors designed to remove excess nitrogen through the metabolic activity of sulfur-consuming microbes. We show that the diversity of microbes in these reactors is both high and highly variable, with distinct community types associated with significant differences in nitrogen removal rate. We also show that the genes encoding the metabolic steps of nitrogen removal are distributed broadly throughout community members, suggesting that the chemical transformations in this system are likely a result of microbes relying on other microbes. These results provide a framework for future studies exploring the contributions of different community members, both in waste removal and in structuring microbial biodiversity.

Single cell genomic and transcriptomic evidence for the use of alternative nitrogen substrates by anammox bacteria.

Anaerobic ammonium oxidation (anammox) contributes substantially to ocean nitrogen loss, particularly in anoxic marine zones (AMZs). Ammonium is scarce in AMZs, raising the hypothesis that organic nitrogen compounds may be ammonium sources for anammox. Biochemical measurements suggest that the organic compounds urea and cyanate can support anammox in AMZs. However, it is unclear if anammox bacteria degrade these compounds to ammonium themselves, or rely on other organisms for this process. Genes for urea degradation have not been found in anammox bacteria, and genomic evidence for cyanate use for anammox is limited to a cyanase gene recovered from the sediment bacterium Candidatus Scalindua profunda. Here, analysis of Ca. Scalindua single amplified genomes from the Eastern Tropical North Pacific AMZ revealed genes for urea degradation and transport, as well as for cyanate degradation. Urease and cyanase genes were transcribed, along with anammox genes, in the AMZ core where anammox rates peaked. Homologs of these genes were also detected in meta-omic datasets from major AMZs in the Eastern Tropical South Pacific and Arabian Sea. These results suggest that anammox bacteria from different ocean regions can directly access organic nitrogen substrates. Future studies should assess if and under what environmental conditions these substrates contribute to the ammonium budget for anammox.

Metabolic potential and in situ activity of marine Marinimicrobia bacteria in an anoxic water column.

Marinimicrobia bacteria are widespread in subeuphotic areas of the oceans and particularly abundant in oxygen minimum zones (OMZs). Information on Marinimicrobia metabolism is sparse, making the biogeochemical influence of this group challenging to predict. Here, metagenome-assembled genomes representing Marinimicrobia subgroups PN262000N21 and ARCTIC96B-7 were retrieved to near completion (97% and 94%) from OMZ metagenomes, with contamination (14.1%) observed only in ARCTIC96B-7. Genes for aerobic carbon monoxide (CO) oxidation, polysulfide metabolism and hydrogen utilization were identified only in PN262000N21, while genes for partial denitrification occurred in both genomes. Transcripts mapping to these genomes increased from <0.3% of total mRNA from the oxic zone to a max of 22% under anoxia. ARCTIC96B-7 transcript representation decreased an order of magnitude from non-sulfidic to sulfidic depths. In contrast, PN262000N21 representation was relatively constant throughout the OMZ, although transcripts encoding sulfur-utilizing proteins, including sulfur transferases, were enriched at sulfidic depths. PN262000N21 transcripts encoding a protein with fibronectin domains similar to those in cellulosome-producing bacteria were also abundant, suggesting a potential for high molecular weight carbon cycling. These data provide omic-level descriptions of metabolic potential and activity in OMZ-associated Marinimicrobia, suggesting differentiation between subgroups with roles in carbon and dissimilatory inorganic nitrogen and sulfur cycling.

Cryptic oxygen cycling in anoxic marine zones.

Oxygen availability drives changes in microbial diversity and biogeochemical cycling between the aerobic surface layer and the anaerobic core in nitrite-rich anoxic marine zones (AMZs), which constitute huge oxygen-depleted regions in the tropical oceans. The current paradigm is that primary production and nitrification within the oxic surface layer fuel anaerobic processes in the anoxic core of AMZs, where 30-50% of global marine nitrogen loss takes place. Here we demonstrate that oxygenic photosynthesis in the secondary chlorophyll maximum (SCM) releases significant amounts of O2 to the otherwise anoxic environment. The SCM, commonly found within AMZs, was dominated by the picocyanobacteria Prochlorococcus spp. Free O2 levels in this layer were, however, undetectable by conventional techniques, reflecting a tight coupling between O2 production and consumption by aerobic processes under apparent anoxic conditions. Transcriptomic analysis of the microbial community in the seemingly anoxic SCM revealed the enhanced expression of genes for aerobic processes, such as nitrite oxidation. The rates of gross O2 production and carbon fixation in the SCM were found to be similar to those reported for nitrite oxidation, as well as for anaerobic dissimilatory nitrate reduction and sulfate reduction, suggesting a significant effect of local oxygenic photosynthesis on Pacific AMZ biogeochemical cycling.

Putative archaeal viruses from the mesopelagic ocean.

Oceanic viruses that infect bacteria, or phages, are known to modulate host diversity, metabolisms, and biogeochemical cycling, while the viruses that infect marine Archaea remain understudied despite the critical ecosystem roles played by their hosts. Here we introduce "MArVD", for Metagenomic Archaeal Virus Detector, an annotation tool designed to identify putative archaeal virus contigs in metagenomic datasets. MArVD is made publicly available through the online iVirus analytical platform. Benchmarking analysis of MArVD showed it to be >99% accurate and 100% sensitive in identifying the 127 known archaeal viruses among the 12,499 viruses in the VirSorter curated dataset. Application of MArVD to 10 viral metagenomes from two depth profiles in the Eastern Tropical North Pacific (ETNP) oxygen minimum zone revealed 43 new putative archaeal virus genomes and large genome fragments ranging in size from 10 to 31 kb. Network-based classifications, which were consistent with marker gene phylogenies where available, suggested that these putative archaeal virus contigs represented six novel candidate genera. Ecological analyses, via fragment recruitment and ordination, revealed that the diversity and relative abundances of these putative archaeal viruses were correlated with oxygen concentration and temperature along two OMZ-spanning depth profiles, presumably due to structuring of the host Archaea community. Peak viral diversity and abundances were found in surface waters, where Thermoplasmata 16S rRNA genes are prevalent, suggesting these archaea as hosts in the surface habitats. Together these findings provide a baseline for identifying archaeal viruses in sequence datasets, and an initial picture of the ecology of such viruses in non-extreme environments.

Piecewise Structural Equation Model (SEM) Disentangles the Environmental Conditions Favoring Diatom Diazotroph Associations (DDAs) in the Western Tropical North Atlantic (WTNA).

Diatom diazotroph associations (DDAs) are important components in the world's oceans, especially in the western tropical north Atlantic (WTNA), where blooms have a significant impact on carbon and nitrogen cycling. However, drivers of their abundances and distribution patterns remain unknown. Here, we examined abundance and distribution patterns for two DDA populations in relation to the Amazon River (AR) plume in the WTNA. Quantitative PCR assays, targeting two DDAs (het-1 and het-2) by their symbiont's nifH gene, served as input in a piecewise structural equation model (SEM). Collections were made during high (spring 2010) and low (fall 2011) flow discharges of the AR. The distributions of dissolved nutrients, chlorophyll-a, and DDAs showed coherent patterns indicative of areas influenced by the AR. A symbiotic Hemiaulus hauckii-Richelia (het-2) bloom (>106 cells L-1) occurred during higher discharge of the AR and was coincident with mesohaline to oceanic (30-35) sea surface salinities (SSS), and regions devoid of dissolved inorganic nitrogen (DIN), low concentrations of both DIP (>0.1 µmol L-1) and Si (>1.0 µmol L-1). The Richelia (het-1) associated with Rhizosolenia was only present in 2010 and at lower densities (10-1.76 × 105nifH copies L-1) than het-2 and limited to regions of oceanic SSS (>36). The het-2 symbiont detected in 2011 was associated with H. membranaceus (>103nifH copies L-1) and were restricted to regions with mesohaline SSS (31.8-34.3), immeasurable DIN, moderate DIP (0.1-0.60 µmol L-1) and higher Si (4.19-22.1 µmol L-1). The piecewise SEM identified a profound direct negative effect of turbidity on the het-2 abundance in spring 2010, while DIP and water turbidity had a more positive influence in fall 2011, corroborating our observations of DDAs at subsurface maximas. We also found a striking difference in the influence of salinity on DDA symbionts suggesting a niche differentiation and preferences in oceanic and mesohaline salinities by het-1 and het-2, respectively. The use of the piecewise SEM to disentangle the complex and concomitant hydrography of the WTNA acting on two biogeochemically relevant populations was novel and underscores its use to predict conditions favoring abundance and distributions of microbial populations.

SAR11 bacteria linked to ocean anoxia and nitrogen loss.

Bacteria of the SAR11 clade constitute up to one half of all microbial cells in the oxygen-rich surface ocean. SAR11 bacteria are also abundant in oxygen minimum zones (OMZs), where oxygen falls below detection and anaerobic microbes have vital roles in converting bioavailable nitrogen to N2 gas. Anaerobic metabolism has not yet been observed in SAR11, and it remains unknown how these bacteria contribute to OMZ biogeochemical cycling. Here, genomic analysis of single cells from the world's largest OMZ revealed previously uncharacterized SAR11 lineages with adaptations for life without oxygen, including genes for respiratory nitrate reductases (Nar). SAR11 nar genes were experimentally verified to encode proteins catalysing the nitrite-producing first step of denitrification and constituted ~40% of OMZ nar transcripts, with transcription peaking in the anoxic zone of maximum nitrate reduction activity. These results link SAR11 to pathways of ocean nitrogen loss, redefining the ecological niche of Earth's most abundant organismal group.

NC10 bacteria in marine oxygen minimum zones.

Bacteria of the NC10 phylum link anaerobic methane oxidation to nitrite denitrification through a unique O2-producing intra-aerobic methanotrophy pathway. A niche for NC10 in the pelagic ocean has not been confirmed. We show that NC10 bacteria are present and transcriptionally active in oceanic oxygen minimum zones (OMZs) off northern Mexico and Costa Rica. NC10 16S rRNA genes were detected at all sites, peaking in abundance in the anoxic zone with elevated nitrite and methane concentrations. Phylogenetic analysis of particulate methane monooxygenase genes further confirmed the presence of NC10. rRNA and mRNA transcripts assignable to NC10 peaked within the OMZ and included genes of the putative nitrite-dependent intra-aerobic pathway, with high representation of transcripts containing the unique motif structure of the nitric oxide (NO) reductase of NC10 bacteria, hypothesized to participate in O2-producing NO dismutation. These findings confirm pelagic OMZs as a niche for NC10, suggesting a role for this group in OMZ nitrogen, methane and oxygen cycling.

Standard filtration practices may significantly distort planktonic microbial diversity estimates.

Fractionation of biomass by filtration is a standard method for sampling planktonic microbes. It is unclear how the taxonomic composition of filtered biomass changes depending on sample volume. Using seawater from a marine oxygen minimum zone, we quantified the 16S rRNA gene composition of biomass on a prefilter (1.6 µm pore-size) and a downstream 0.2 µm filter over sample volumes from 0.05 to 5 L. Significant community shifts occurred in both filter fractions, and were most dramatic in the prefilter community. Sequences matching Vibrionales decreased from ~40 to 60% of prefilter datasets at low volumes (0.05-0.5 L) to less than 5% at higher volumes, while groups such at the Chromatiales and Thiohalorhabdales followed opposite trends, increasing from minor representation to become the dominant taxa at higher volumes. Groups often associated with marine particles, including members of the Deltaproteobacteria, Planctomycetes, and Bacteroidetes, were among those showing the greatest increase with volume (4 to 27-fold). Taxon richness (97% similarity clusters) also varied significantly with volume, and in opposing directions depending on filter fraction, highlighting potential biases in community complexity estimates. These data raise concerns for studies using filter fractionation for quantitative comparisons of aquatic microbial diversity, for example between free-living and particle-associated communities.