Size-fraction partitioning of community gene transcription and nitrogen metabolism in a marine oxygen minimum zone

Nitrite-oxidizing bacteria adapted to low-oxygen conditions dominate nitrite oxidation in marine oxygen minimum zones

Nitrite is a central molecule in the nitrogen cycle because nitrite oxidation to nitrate (an aerobic process) retains fixed nitrogen in a system and its reduction to dinitrogen gas (anaerobic) reduces the fixed nitrogen inventory. Despite its acknowledged requirement for oxygen, nitrite oxidation is observed in oxygen-depleted layers of the ocean's oxygen minimum zones (OMZs), challenging the current understanding of OMZ nitrogen cycling. Previous attempts to determine whether nitrite-oxidizing bacteria in the anoxic layer differ from known nitrite oxidizers in the open ocean were limited by cultivation difficulties and sequencing depth. Here, we construct 31 draft genomes of nitrite-oxidizing bacteria from global OMZs. The distribution of nitrite oxidation rates, abundance and expression of nitrite oxidoreductase genes, and relative abundance of nitrite-oxidizing bacterial draft genomes from the same samples all show peaks in the core of the oxygen-depleted zone (ODZ) and are all highly correlated in depth profiles within the major ocean oxygen minimum zones. The ODZ nitrite oxidizers are not found in the Tara Oceans global dataset (the most complete oxic ocean dataset), and the major nitrite oxidizers found in the oxygenated ocean do not occur in ODZ waters. A pangenomic analysis shows the ODZ nitrite oxidizers have distinct gene clusters compared to oxic nitrite oxidizers and are microaerophilic. These findings all indicate the existence of nitrite oxidizers whose niche is oxygen-deficient seawater. Thus, specialist nitrite-oxidizing bacteria are responsible for fixed nitrogen retention in marine oxygen minimum zones, with implications for control of the ocean's fixed nitrogen inventory.

Age, metabolisms, and potential origin of dominant anammox bacteria in the global oxygen-deficient zones

Anammox bacteria inhabiting oxygen-deficient zones (ODZs) are a major functional group mediating fixed nitrogen loss in the global ocean. However, many basic questions regarding the diversity, broad metabolisms, origin, and adaptive mechanisms of ODZ anammox bacteria remain unaddressed. Here we report two novel metagenome-assembled genomes of anammox bacteria affiliated with the Scalindua genus, which represent most, if not all, of the anammox bacteria in the global ODZs. Metagenomic read-recruiting and comparison with historical data show that they are ubiquitously present in all three major ODZs. Beyond the core anammox metabolism, both organisms contain cyanase, and the more dominant one encodes a urease, indicating most ODZ anammox bacteria can utilize cyanate and urea in addition to ammonium. Molecular clock analysis suggests that the evolutionary radiation of these bacteria into ODZs occurred no earlier than 310 million years ago, ~1 billion years after the emergence of the earliest modern-type ODZs. Different strains of the ODZ Scalindua species are also found in benthic sediments, and the first ODZ Scalindua is likely derived from the benthos. Compared to benthic strains of the same clade, ODZ Scalindua uniquely encodes genes for urea utilization but has lost genes related to growth arrest, flagellum synthesis, and chemotaxis, presumably for adaptation to thrive in the global ODZ waters. Our findings expand the known metabolisms and evolutionary history of the bacteria controlling the global nitrogen budget.

Particle-associated denitrification is the primary source of N2O in oxic coastal waters

The heavily human-perturbed coastal oceans are hotspots of nitrous oxide (N2O) emission to the atmosphere. The processes underpinning the N2O flux, however, remain poorly understood, leading to large uncertainties in assessing global N2O budgets. Using a suite of nitrogen isotope labeling experiments, we show that multiple processes contribute to N2O production throughout the estuarine-coastal gradient, sustaining intensive N2O flux to the atmosphere. Unexpectedly, denitrification, rather than ammonia oxidation as previously assumed, constitutes the major source of N2O in well-oxygenated coastal waters. Size-fractionated manipulation experiments with gene analysis further reveal niche partitioning of ammonia oxidizers and denitrifiers across the particle size spectrum; denitrification dominated on large particles and ammonia oxidizers on small particles. Total N2O production rate increases with substrate and particle concentrations, suggesting a crucial interplay between nutrients and particles in controlling N2O production. The controlling factors identified here may help understand climate feedback mechanisms between human activity and coastal oceans.

The diversity, community dynamics, and interactions of the microbiome in the world's deepest blue hole: insights into extreme environmental response patterns and tolerance of marine microorganisms

IMPORTANCE

This study comprehensively examined the community dynamics, functional profiles, and interactions of the microbiome in the world's deepest blue hole. The findings revealed a positive correlation between the α-diversities of Symbiodiniaceae and archaea, indicating the potential reliance of Symbiodiniaceae on archaea in an extreme environment resulting from a partial niche overlap. The negative association between the α-diversity and β-diversity of the bacterial community suggested that the change rule of the bacterial community was consistent with the Anna Karenina effects. The core microbiome comprised nine microbial taxa, highlighting their remarkable tolerance and adaptability to sharp environmental gradient variations. Bacteria and archaea played significant roles in carbon, nitrogen, and sulfur cycles, while fungi contributed to carbon metabolism. This study advanced our understanding of the community dynamics, response patterns, and resilience of microorganisms populating the world's deepest blue hole, thereby facilitating further ecological and evolutional exploration of microbiomes in diverse extreme environments.

Microscale dynamics promote segregated denitrification in diatom aggregates sinking slowly in bulk oxygenated seawater

Sinking marine particles drive the biological pump that naturally sequesters carbon from the atmosphere. Despite their small size, the compartmentalized nature of particles promotes intense localized metabolic activity by their bacterial colonizers. Yet the mechanisms promoting the onset of denitrification, a metabolism that arises once oxygen is limiting, remain to be established. Here we show experimentally that slow sinking aggregates composed of marine diatoms-important primary producers for global carbon export-support active denitrification even among bulk oxygenated water typically thought to exclude anaerobic metabolisms. Denitrification occurs at anoxic microsites distributed throughout a particle and within microns of a particle's boundary, and fluorescence-reporting bacteria show nitrite can be released into the water column due to segregated dissimilatory reduction of nitrate and nitrite. Examining intact and broken diatoms as organic sources, we show slowly leaking cells promote more bacterial growth, allow particles to have lower oxygen, and generally support greater denitrification.

Partitioning of the denitrification pathway and other nitrite metabolisms within global oxygen deficient zones

Oxygen deficient zones (ODZs) account for about 30% of total oceanic fixed nitrogen loss via processes including denitrification, a microbially mediated pathway proceeding stepwise from NO3- to N2. This process may be performed entirely by complete denitrifiers capable of all four enzymatic steps, but many organisms possess only partial denitrification pathways, either producing or consuming key intermediates such as the greenhouse gas N2O. Metagenomics and marker gene surveys have revealed a diversity of denitrification genes within ODZs, but whether these genes co-occur within complete or partial denitrifiers and the identities of denitrifying taxa remain open questions. We assemble genomes from metagenomes spanning the ETNP and Arabian Sea, and map these metagenome-assembled genomes (MAGs) to 56 metagenomes from all three major ODZs to reveal the predominance of partial denitrifiers, particularly single-step denitrifiers. We find niche differentiation among nitrogen-cycling organisms, with communities performing each nitrogen transformation distinct in taxonomic identity and motility traits. Our collection of 962 MAGs presents the largest collection of pelagic ODZ microorganisms and reveals a clearer picture of the nitrogen cycling community within this environment.

A compendium of bacterial and archaeal single-cell amplified genomes from oxygen deficient marine waters

Oxygen-deficient marine waters referred to as oxygen minimum zones (OMZs) or anoxic marine zones (AMZs) are common oceanographic features. They host both cosmopolitan and endemic microorganisms adapted to low oxygen conditions. Microbial metabolic interactions within OMZs and AMZs drive coupled biogeochemical cycles resulting in nitrogen loss and climate active trace gas production and consumption. Global warming is causing oxygen-deficient waters to expand and intensify. Therefore, studies focused on microbial communities inhabiting oxygen-deficient regions are necessary to both monitor and model the impacts of climate change on marine ecosystem functions and services. Here we present a compendium of 5,129 single-cell amplified genomes (SAGs) from marine environments encompassing representative OMZ and AMZ geochemical profiles. Of these, 3,570 SAGs have been sequenced to different levels of completion, providing a strain-resolved perspective on the genomic content and potential metabolic interactions within OMZ and AMZ microbiomes. Hierarchical clustering confirmed that samples from similar oxygen concentrations and geographic regions also had analogous taxonomic compositions, providing a coherent framework for comparative community analysis.

Chemical Links Between Redox Conditions and Estimated Community Proteomes from 16S rRNA and Reference Protein Sequences

Environmental influences on community structure are often assessed through multivariate analyses in order to relate microbial abundances to separately measured physicochemical variables. However, genes and proteins are themselves chemical entities; in combination with genome databases, differences in microbial abundances directly encode for chemical variability. We predicted that the carbon oxidation state of estimated community proteomes, obtained by combining taxonomic abundances from published 16S rRNA gene sequencing datasets with reference microbial proteomes from the NCBI Reference Sequence (RefSeq) database, would reflect environmental oxidation-reduction conditions. Analysis of multiple datasets confirms the geobiochemical predictions for environmental redox gradients in hydrothermal systems, stratified lakes and marine environments, and shale gas wells. The geobiochemical signal is largest for the steep redox gradients associated with hydrothermal systems and between injected water and produced fluids from shale gas wells, demonstrating that microbial community composition can be a chemical proxy for environmental redox gradients. Although estimates of oxidation state from 16S amplicon and metagenomic sequences are correlated, the 16S-based estimates show stronger associations with redox gradients in some environments.

Diverse secondary metabolites are expressed in particle-associated and free-living microorganisms of the permanently anoxic Cariaco Basin

Secondary metabolites play essential roles in ecological interactions and nutrient acquisition, and are of interest for their potential uses in medicine and biotechnology. Genome mining for biosynthetic gene clusters (BGCs) can be used for the discovery of new compounds. Here, we use metagenomics and metatranscriptomics to analyze BGCs in free-living and particle-associated microbial communities through the stratified water column of the Cariaco Basin, Venezuela. We recovered 565 bacterial and archaeal metagenome-assembled genomes (MAGs) and identified 1154 diverse BGCs. We show that differences in water redox potential and microbial lifestyle (particle-associated vs. free-living) are associated with variations in the predicted composition and production of secondary metabolites. Our results indicate that microbes, including understudied clades such as Planctomycetota, potentially produce a wide range of secondary metabolites in these anoxic/euxinic waters.

Composition and assembly of the bacterial community in the overlying waters of the coral reef of China's Xisha Islands

Coral reef ecosystems are one of the most diverse and productive habitats on Earth. Microbes in the reef-overlying waters are key players in maintaining this ecosystem through regulating biogeochemical and ecological processes. However, the composition structure and assembly mechanism of microbial community in the reef-overlying waters remain largely unknown. In the present study, the bacterial communities from the overlying waters of atolls and fringing reefs as well as the surface waters of the adjacent open ocean of the Xisha Islands in the South China Sea were investigated using 16S rRNA high-throughput sequencing combined with a size-fractionation strategy. The results showed that environments of all sampling stations were similar, characterized by an almost complete lack of inorganic nutrients such as nitrogen and phosphorus. Proteobacteria, Cyanobacteria and Bacteroidetes were the dominant phyla, and Synechococcus was most abundant at the genus level in both large fraction (LF; 1.6-200 μm) and small fraction (SF; 0.2-1.6 μm) communities. Only a slight difference in community composition between LF and SF samples was observed. The bacterial communities among the three habitat types showed noticeable differences, and the bacterial composition among the atoll reefs was more varied than that among the fringing reefs. The similarity of bacterial communities significantly declined with the increasing geographic distance, and stochastic processes were more important than deterministic processes in bacterial community assembly. This study sheds lights on the bacterial biodiversity of coral reefs and the importance of stochastic process in structuring bacterial communities.

Habitat-dependent prokaryotic microbial community, potential keystone species, and network complexity in a subtropical estuary

Microbes (e.g., bacteria and archaea) are indispensable components for the key biological processes of estuarine ecosystems and three main habitats (sediment, particle, and water) are harboring diverse estuarine microbes. However, we still know little about how the microbial community structures, potential keystone species, and network properties change among these three habitats in estuarine ecosystems. In this study, we collected size-fractioned water and sediment samples from the Pearl River Estuary to reveal their microbial diversity, community structures, network properties, and potential keystone taxa. We found that the sediment microbial community was remarkably more diverse than particle-attached (PA) and free-living (FL) communities, whereas its ecological network was less complex in terms of node distance and connectivity. TOC was determined as the main driver of sediment community, while the PA and FL communities were predominantly shaped by NO₂^-, non-ionic ammonia (NH) and pH. Among the bulk water, there were no significant differences between PA and FL communities in diversity, community structure, and network complexity. However, the PA community was more susceptible to metal elements, suggesting their higher level of involvement in physiological metabolism. Potential keystone taxa among community networks were taxonomically divergent in three habitats. Specifically, Synechococcales (Cyanobacteria) and Actinomarinales (Actinobacteria) exclusively served as the module-hubs in FL network, while members from phylum Proteobacteria and Bacteroidetes were the module-hubs and connectors in PA network. Potential keystone taxa in sediment network were more diverse and covered 9 phyla, including the only archaeal lineage Bathyarchaeia (Crenarchaeota). Overall, our study provided more detailed information about estuarine microbial communities in three habitats, especially the potential keystone species, which provided new perspectives on evaluating further effects of anthropogenic disturbances on estuarine microbes and facilitated the environment monitoring based on microbial community.

Insights into Prokaryotic Community and Its Potential Functions in Nitrogen Metabolism in the Bay of Bengal, a Pronounced Oxygen Minimum Zone

Ocean oxygen minimum zones (OMZs) around the global ocean are expanding both horizontally and vertically. Multiple studies have identified the significant influence of anoxic conditions (≤1 μM O₂) on marine prokaryotic communities and biogeochemical cycling of elements. However, little attention has been paid to the expanding low-oxygen zones where the oxygen level is still above the anoxic level. Here, we studied the abundance and taxonomic and functional profiles of prokaryotic communities in the Bay of Bengal (BoB), where the oxygen concentration is barely above suboxic level (5 μM O₂). We found the sinking of Trichodesmium into deep water was far more efficient than that of Prochlorococcus, suggesting Trichodesmium blooms might be an essential carbon and nitrogen source for the maintenance of the BoB OMZ. In addition to the shift in the prokaryotic community composition, the abundance of some functional genes also changed with the change of oxygen concentration. Compared to oxic (>60 μM O₂) Tara Ocean and high-hypoxic (>20 to ≤60 μM O₂) BoB samples, we found more SAR11-nar sequences (responsible for reducing nitrate to nitrite) in low-hypoxic (>5 to ≤20 μM O₂) BoB waters. This suggested SAR11-nar genes would be more widespread due to the expansion of OMZs. It seems that the nitrite-N was not further reduced to nitrogen through denitrification but likely oxidized to nitrate by Nitrospinae in the BoB OMZ and then accumulated in the form of nitrate-N. However, the lack of N₂ production in the BoB would change if the BoB OMZ became anoxic. Together, these results suggested that reduction of oxygen concentration and OMZ expansion may increase the use of nitrate by SAR11 and N₂ production in the BoB.

IMPORTANCE Recognizing the prokaryotic community and its functions in hypoxic (>5 to ≤60 μM O₂) environments before further expansion of OMZs is critical. We demonstrate the prokaryotic community and its potential functions in nitrogen metabolism in the Bay of Bengal (BoB), where oxygen concentration is barely above suboxic level. This study highlighted that Trichodesmium might be an essential carbon and nitrogen source in the maintenance of the BoB OMZ. Additionally, we suggest that the lack of N₂ production in the BoB would change if the BoB OMZ became anoxic, and the expansion of OMZs in the global ocean may potentially increase the use of nitrate by SAR11.

Anaerobic methane oxidation in a coastal oxygen minimum zone: spatial and temporal dynamics

Coastal waters are a major source of marine methane to the atmosphere. Particularly high concentrations of this potent greenhouse gas are found in anoxic waters, but it remains unclear if and to what extent anaerobic methanotrophs mitigate the methane flux. Here we investigate the long-term dynamics in methanotrophic activity and the methanotroph community in the coastal oxygen minimum zone (OMZ) of Golfo Dulce, Costa Rica, combining biogeochemical analyses, experimental incubations and 16S rRNA gene sequencing over 3 consecutive years. Our results demonstrate a stable redox zonation across the years with high concentrations of methane (up to 1.7 μmol L-1 ) in anoxic bottom waters. However, we also measured high activities of anaerobic methane oxidation in the OMZ core (rate constant, k, averaging 30 yr-1 in 2018 and 8 yr-1 in 2019-2020). The OPU3 and Deep Sea-1 clades of the Methylococcales were implicated as conveyors of the activity, peaking in relative abundance 5-25 m below the oxic-anoxic interface and in the deep anoxic water respectively. Although their genetic capacity for anaerobic methane oxidation remains unexplored, their sustained high relative abundance indicates an adaptation of these clades to the anoxic, methane-rich OMZ environment, allowing them to play major roles in mitigating methane fluxes.

Viral community analysis in a marine oxygen minimum zone indicates increased potential for viral manipulation of microbial physiological state

Microbial communities in oxygen minimum zones (OMZs) are known to have significant impacts on global biogeochemical cycles, but viral influence on microbial processes in these regions are much less studied. Here we provide baseline ecological patterns using microscopy and viral metagenomics from the Eastern Tropical North Pacific (ETNP) OMZ region that enhance our understanding of viruses in these climate-critical systems. While extracellular viral abundance decreased below the oxycline, viral diversity and lytic infection frequency remained high within the OMZ, demonstrating that viral influences on microbial communities were still substantial without the detectable presence of oxygen. Viral community composition was strongly related to oxygen concentration, with viral populations in low-oxygen portions of the water column being distinct from their surface layer counterparts. However, this divergence was not accompanied by the expected differences in viral-encoded auxiliary metabolic genes (AMGs) relating to nitrogen and sulfur metabolisms that are known to be performed by microbial communities in these low-oxygen and anoxic regions. Instead, several abundant AMGs were identified in the oxycline and OMZ that may modulate host responses to low-oxygen stress. We hypothesize that this is due to selection for viral-encoded genes that influence host survivability rather than modulating host metabolic reactions within the ETNP OMZ. Together, this study shows that viruses are not only diverse throughout the water column in the ETNP, including the OMZ, but their infection of microorganisms has the potential to alter host physiological state within these biogeochemically important regions of the ocean.

Tracking the impacts of nutrient inputs on estuary ecosystem function

Estuaries are one of the most impacted coastal environments globally, subjected to multiple stressors from urban, industry and coastal development. With increasing anthropogenic activity surrounding estuarine systems, sewage inputs have become a common concern. Stable isotope analysis provides a well-established tool to investigate the incorporation of nitrogen into marine organisms and identify major nutrient sources. Benthic macroinvertebrate communities are often used as bioindicators in ecological studies as they typically display predictable responses to anthropogenic pressures, however have a suite of limitations and costs associated with their use. 16S rDNA amplicon sequencing techniques allow for investigation of the microbial communities inhabiting complex environmental samples, with potential as a tool in the ecological assessment of pollution. These communities have not yet been adequately considered for ecological studies and biomonitoring, with a need to better understand interactions with environmental stressors and implications for ecosystem function. This study used a combination of stable isotope analysis to trace the uptake of anthropogenic nitrogen in biota, traditional assessment of benthic macroinvertebrate communities, and 16S rDNA genotyping of benthic microbial communities. Stable isotope analysis of seagrass and epiphytes identified multiple treated and untreated sewage inputs, ranges of 5.2-7.2‰ and 1.9-4.0‰ for δ¹⁵N respectively, as the dominant nitrogen source at specific locations. The benthic macroinvertebrate community reflected these inputs with shifts in dominant taxa and high abundances of polychaetes at some sites. Microbial communities provided a sensitive indication of impact with a breadth of information not available using traditional techniques. Composition and predicted function reflected sewage inputs, particularly within sediments, with the relative abundance of specific taxa and putative pathogens linked to these inputs. This research supports the growing body of evidence that benthic microbial communities respond rapidly to anthropogenic stressors and have potential as a monitoring tool in urban estuarine systems.

The marine nitrogen cycle: new developments and global change

The ocean is home to a diverse and metabolically versatile microbial community that performs the complex biochemical transformations that drive the nitrogen cycle, including nitrogen fixation, assimilation, nitrification and nitrogen loss processes. In this Review, we discuss the wealth of new ocean nitrogen cycle research in disciplines from metaproteomics to global biogeochemical modelling and in environments from productive estuaries to the abyssal deep sea. Influential recent discoveries include new microbial functional groups, novel metabolic pathways, original conceptual perspectives and ground-breaking analytical capabilities. These emerging research directions are already contributing to urgent efforts to address the primary challenge facing marine microbiologists today: the unprecedented onslaught of anthropogenic environmental change on marine ecosystems. Ocean warming, acidification, nutrient enrichment and seawater stratification have major effects on the microbial nitrogen cycle, but widespread ocean deoxygenation is perhaps the most consequential for the microorganisms involved in both aerobic and anaerobic nitrogen transformation pathways. In turn, these changes feed back to the global cycles of greenhouse gases such as carbon dioxide and nitrous oxide. At a time when our species casts a lengthening shadow across all marine ecosystems, timely new advances offer us unique opportunities to understand and better predict human impacts on nitrogen biogeochemistry in the changing ocean of the Anthropocene.

Sediments alleviate the inhibition effects of antibiotics on denitrification: Functional gene, microbial community, and antibiotic resistance gene analysis

Both antibiotics and sediments can affect the denitrification in aquatic systems. However, little is known how antibiotics influence the denitrification in the presence of sediments. Here, the effects of antibiotics (sulfamethoxazole, tetracycline and ofloxacin) on denitrification in the absence and presence of sediments were investigated. The influencing mechanisms were revealed by quantifying the denitrification functional genes (DNGs), 16S-seq of bacteria, and antibiotic resistance genes (ARGs). The results showed that the presence of antibiotics inhibited NO₃-N reduction by decreasing the abundances of narG, nirK, nosZ, total DNGs, and denitrifying bacteria. However, the inhibition effect was alleviated by sediments, which promoted the growth of bacteria and decreased the selective pressure of antibiotics as the vector of bacteria and antibiotics, thus increasing the abundances of denitrifying bacteria and all the DNGs. Partial least-squares path model disclosed that antibiotics had negative effects on bacteria, ARGs and DNGs, while sediments had negative effects on ARGs but positive effects on bacteria and DNGs. The network analysis further revealed the close relation of the genera Bacillus, Acinetobacter, and Enterobacter with the ARGs and DNGs. The findings are helpful to understand the denitrification in antibiotic-polluted natural waters.

Substantial oxygen consumption by aerobic nitrite oxidation in oceanic oxygen minimum zones

Oceanic oxygen minimum zones (OMZs) are globally significant sites of biogeochemical cycling where microorganisms deplete dissolved oxygen (DO) to concentrations <20 µM. Amid intense competition for DO in these metabolically challenging environments, aerobic nitrite oxidation may consume significant amounts of DO and help maintain low DO concentrations, but this remains unquantified. Using parallel measurements of oxygen consumption rates and ¹⁵N-nitrite oxidation rates applied to both water column profiles and oxygen manipulation experiments, we show that the contribution of nitrite oxidation to overall DO consumption systematically increases as DO declines below 2 µM. Nitrite oxidation can account for all DO consumption only under DO concentrations <393 nM found in and below the secondary chlorophyll maximum. These patterns are consistent across sampling stations and experiments, reflecting coupling between nitrate reduction and nitrite-oxidizing Nitrospina with high oxygen affinity (based on isotopic and omic data). Collectively our results demonstrate that nitrite oxidation plays a pivotal role in the maintenance and biogeochemical dynamics of OMZs.

Microbial Ecology of Oxygen Minimum Zones Amidst Ocean Deoxygenation

Oxygen minimum zones (OMZs) have substantial effects on the global ecology and biogeochemical processes of marine microbes. However, the diversity and activity of OMZ microbes and their trophic interactions are only starting to be documented, especially in regard to the potential roles of viruses and protists. OMZs have expanded over the past 60 years and are predicted to expand due to anthropogenic climate change, furthering the need to understand these regions. This review summarizes the current knowledge of OMZ formation, the biotic and abiotic factors involved in OMZ expansion, and the microbial ecology of OMZs, emphasizing the importance of bacteria, archaea, viruses, and protists. We describe the recognized roles of OMZ microbes in carbon, nitrogen, and sulfur cycling, the potential of viruses in altering host metabolisms involved in these cycles, and the control of microbial populations by grazers and viruses. Further, we highlight the microbial community composition and roles of these organisms in oxic and anoxic depths within the water column and how these differences potentially inform how microbial communities will respond to deoxygenation. Additionally, the current literature on the alteration of microbial communities by other key climate change parameters such as temperature and pH are considered regarding how OMZ microbes might respond to these pressures. Finally, we discuss what knowledge gaps are present in understanding OMZ microbial communities and propose directions that will begin to close these gaps.

Insights into the cyanosphere: capturing the respective metabolisms of cyanobacteria and chemotrophic bacteria in natural conditions?

Specific interactions have been highlighted between cyanobacteria and chemotrophic bacteria within the cyanosphere, suggesting that nutrients recycling could be optimized by cyanobacteria/bacteria exchanges. In order to determine the respective metabolic roles of the cyanobacterial and bacterial consortia (microbiome), a day-night metatranscriptomic analysis was performed on Dolichospermum sp. (N₂ -fixer) and Microcystis sp. (non N₂ -fixer) natural blooms occurring successively within a French peri-urban lake. The taxonomical and functional analysis of the metatranscriptoms have highlighted specific association of bacteria within the cyanosphere, driven by the cyanobacteria identity, without strongly modifying the functional composition of the microbiomes, suggesting functional redundancy within the cyanosphere. Moreover, the functional composition of these active communities was driven by the living mode. During the two successive bloom events, it appeared that NH₄ ⁺ (newly fixed and/or allochthonous) was preferentially transformed into amino acids for the both the microbiome and the cyanobacteria, while phosphate metabolism was enhanced, suggesting that due to a high cellular growth, P limitation might take place within the cyanosphere consortium.

Small sinking particles control anammox rates in the Peruvian oxygen minimum zone

Anaerobic oxidation of ammonium (anammox) in oxygen minimum zones (OMZs) is a major pathway of oceanic nitrogen loss. Ammonium released from sinking particles has been suggested to fuel this process. During cruises to the Peruvian OMZ in April–June 2017 we found that anammox rates are strongly correlated with the volume of small particles (128–512 µm), even though anammox bacteria were not directly associated with particles. This suggests that the relationship between anammox rates and particles is related to the ammonium released from particles by remineralization. To investigate this, ammonium release from particles was modelled and theoretical encounters of free-living anammox bacteria with ammonium in the particle boundary layer were calculated. These results indicated that small sinking particles could be responsible for ~75% of ammonium release in anoxic waters and that free-living anammox bacteria frequently encounter ammonium in the vicinity of smaller particles. This indicates a so far underestimated role of abundant, slow-sinking small particles in controlling oceanic nutrient budgets, and furthermore implies that observations of the volume of small particles could be used to estimate N-loss across large areas.

Up to 40% of the ocean’s fixed nitrogen is lost in oxygen minimum zones (OMZs) by anammox, but despite the importance of this process, nitrogen loss patterns in OMZs are difficult to predict. Here the authors show that ammonium release from small particles is a major control of anammox in the Peruvian OMZ.

Deep ocean metagenomes provide insight into the metabolic architecture of bathypelagic microbial communities

The deep sea, the largest ocean's compartment, drives planetary-scale biogeochemical cycling. Yet, the functional exploration of its microbial communities lags far behind other environments. Here we analyze 58 metagenomes from tropical and subtropical deep oceans to generate the Malaspina Gene Database. Free-living or particle-attached lifestyles drive functional differences in bathypelagic prokaryotic communities, regardless of their biogeography. Ammonia and CO oxidation pathways are enriched in the free-living microbial communities and dissimilatory nitrate reduction to ammonium and H2 oxidation pathways in the particle-attached, while the Calvin Benson-Bassham cycle is the most prevalent inorganic carbon fixation pathway in both size fractions. Reconstruction of the Malaspina Deep Metagenome-Assembled Genomes reveals unique non-cyanobacterial diazotrophic bacteria and chemolithoautotrophic prokaryotes. The widespread potential to grow both autotrophically and heterotrophically suggests that mixotrophy is an ecologically relevant trait in the deep ocean. These results expand our understanding of the functional microbial structure and metabolic capabilities of the largest Earth aquatic ecosystem.

Denitrifying bacteria respond to and shape microscale gradients within particulate matrices

Heterotrophic denitrification enables facultative anaerobes to continue growing even when limited by oxygen (O₂) availability. Particles in particular provide physical matrices characterized by reduced O₂ permeability even in well-oxygenated bulk conditions, creating microenvironments where microbial denitrifiers may proliferate. Whereas numerical particle models generally describe denitrification as a function of radius, here we provide evidence for heterogeneity of intraparticle denitrification activity due to local interactions within and among microcolonies. Pseudomonas aeruginosa cells and microcolonies act to metabolically shade each other, fostering anaerobic processes just microns from O₂-saturated bulk water. Even within well-oxygenated fluid, suboxia and denitrification reproducibly developed and migrated along sharp 10 to 100 µm gradients, driven by the balance of oxidant diffusion and local respiration. Moreover, metabolic differentiation among densely packed cells is dictated by the diffusional supply of O₂, leading to distinct bimodality in the distribution of nitrate and nitrite reductase expression. The initial seeding density controls the speed at which anoxia develops, and even particles seeded with few bacteria remain capable of becoming anoxic. Our empirical results capture the dynamics of denitrifier gene expression in direct association with O₂ concentrations over microscale physical matrices, providing observations of the co-occurrence and spatial arrangement of aerobic and anaerobic processes.

Smriga et al. investigate heterotrophic denitrification in particulate matrices. Their model shows that Pseudomonas aeruginosa cells and microcolonies metabolically shade each other, allowing denitrification to occur microns away from O₂-saturated bulk water, and simulates microscale effects of a bulk system.

Microbial niche differentiation explains nitrite oxidation in marine oxygen minimum zones

Nitrite is a pivotal component of the marine nitrogen cycle. The fate of nitrite determines the loss or retention of fixed nitrogen, an essential nutrient for all organisms. Loss occurs via anaerobic nitrite reduction to gases during denitrification and anammox, while retention occurs via nitrite oxidation to nitrate. Nitrite oxidation is usually represented in biogeochemical models by one kinetic parameter and one oxygen threshold, below which nitrite oxidation is set to zero. Here we find that the responses of nitrite oxidation to nitrite and oxygen concentrations vary along a redox gradient in a Pacific Ocean oxygen minimum zone, indicating niche differentiation of nitrite-oxidizing assemblages. Notably, we observe the full inhibition of nitrite oxidation by oxygen addition and nitrite oxidation coupled with nitrogen loss in the absence of oxygen consumption in samples collected from anoxic waters. Nitrite-oxidizing bacteria, including novel clades with high relative abundance in anoxic depths, were also detected in the same samples. Mechanisms corresponding to niche differentiation of nitrite-oxidizing bacteria across the redox gradient are considered. Implementing these mechanisms in biogeochemical models has a significant effect on the estimated fixed nitrogen budget.

Potential virus-mediated nitrogen cycling in oxygen-depleted oceanic waters

Viruses play an important role in the ecology and biogeochemistry of marine ecosystems. Beyond mortality and gene transfer, viruses can reprogram microbial metabolism during infection by expressing auxiliary metabolic genes (AMGs) involved in photosynthesis, central carbon metabolism, and nutrient cycling. While previous studies have focused on AMG diversity in the sunlit and dark ocean, less is known about the role of viruses in shaping metabolic networks along redox gradients associated with marine oxygen minimum zones (OMZs). Here, we analyzed relatively quantitative viral metagenomic datasets that profiled the oxygen gradient across Eastern Tropical South Pacific (ETSP) OMZ waters, assessing whether OMZ viruses might impact nitrogen (N) cycling via AMGs. Identified viral genomes encoded six N-cycle AMGs associated with denitrification, nitrification, assimilatory nitrate reduction, and nitrite transport. The majority of these AMGs (80%) were identified in T4-like Myoviridae phages, predicted to infect Cyanobacteria and Proteobacteria, or in unclassified archaeal viruses predicted to infect Thaumarchaeota. Four AMGs were exclusive to anoxic waters and had distributions that paralleled homologous microbial genes. Together, these findings suggest viruses modulate N-cycling processes within the ETSP OMZ and may contribute to nitrogen loss throughout the global oceans thus providing a baseline for their inclusion in the ecosystem and geochemical models.

Description of Candidatus Mesopelagibacter carboxydoxydans and Candidatus Anoxipelagibacter denitrificans: Nitrate-reducing SAR11 genera that dominate mesopelagic and anoxic marine zones

The diverse and ubiquitous members of the SAR11 lineage (Alphaproteobacteria) represent up to 30-40% of the surface and mesopelagic oceanic microbial communities. However, the molecular and ecological mechanisms that differentiate closely related, yet distinct, SAR11 members that often co-occur under similar environmental conditions remain speculative. Recently, two mesopelagic and oxygen minimum zone (OMZ)-associated subclades of SAR11 (Ic and IIa.A) were described using single-cell amplified genomes (SAGs) linked to nitrate reduction in OMZs. In this current study, the collection of genomes belonging to these two subclades was expanded with thirteen new metagenome-assembled genomes (MAGs), thus providing a more detailed phylogenetic and functional characterization of these subclades. Gene content-based predictions of metabolic functions revealed similarities in central carbon metabolism between subclades Ic and IIa.A and surface SAR11 clades, with small variations in central pathways. These variations included more versatile sulfur assimilation pathways, as well as a previously predicted capacity for nitrate reduction that conferred unique versatility on mesopelagic-adapted clades compared to their surface counterparts. Finally, consistent with previously reported abundances of carbon monoxide (CO) in surface and mesopelagic waters, subclades Ia (surface) and Ic (mesopelagic) have the genetic potential to oxidize carbon monoxide (CO), presumably taking advantage of this abundant compound as an electron donor. Based on genomic analyses, environmental distribution and metabolic reconstruction, we propose two new SAR11 genera, Ca. Mesopelagibacter carboxydoxydans (subclade Ic) and Ca. Anoxipelagibacter denitrificans (subclade IIa.A), which represent members of the mesopelagic and OMZ-adapted SAR11 clades.

Depth-Dependent Variables Shape Community Structure and Functionality in the Prince Edward Islands

Physicochemical variables limit and control the distribution of microbial communities in all environments. In the oceans, this may significantly influence functional processes such the consumption of dissolved organic material and nutrient sequestration. Yet, the relative contributions of physical factors, such as water mass variability and depth, on functional processes are underexplored. We assessed microbial community structure and functionality in the Prince Edward Islands (PEIs) using 16S rRNA gene amplicon analysis and extracellular enzymatic activity assays, respectively. We found that depth and nutrients substantially drive the structural patterns of bacteria and archaea in this region. Shifts from epipelagic to bathypelagic zones were linked to decreases in the activities of several extracellular enzymes. These extracellular enzymatic activities were positively correlated with several phyla including several Alphaproteobacteria (including members of the SAR 11 clade and order Rhodospirillales) and Cyanobacteria. We show that depth-dependent variables may be essential drivers of community structure and functionality in the PEIs.

Microbial N2O consumption in and above marine N2O production hotspots

The ocean is a net source of N₂O, a potent greenhouse gas and ozone-depleting agent. However, the removal of N₂O via microbial N₂O consumption is poorly constrained and rate measurements have been restricted to anoxic waters. Here we expand N₂O consumption measurements from anoxic zones to the sharp oxygen gradient above them, and experimentally determine kinetic parameters in both oxic and anoxic seawater for the first time. We find that the substrate affinity, O₂ tolerance, and community composition of N₂O-consuming microbes in oxic waters differ from those in the underlying anoxic layers. Kinetic parameters determined here are used to model in situ N₂O production and consumption rates. Estimated in situ rates differ from measured rates, confirming the necessity to consider kinetics when predicting N₂O cycling. Microbes from the oxic layer consume N₂O under anoxic conditions at a much faster rate than microbes from anoxic zones. These experimental results are in keeping with model results which indicate that N₂O consumption likely takes place above the oxygen deficient zone (ODZ). Thus, the dynamic layer with steep O₂ and N₂O gradients right above the ODZ is a previously ignored potential gatekeeper of N₂O and should be accounted for in the marine N₂O budget.

Non-denitrifier nitrous oxide reductases dominate marine biomes

Microbial enzymes often occur as distinct variants that share the same substrate but differ in substrate affinity, sensitivity to environmental conditions, or phylogenetic ancestry. Determining where variants occur in the environment helps identify thresholds that constrain microbial cycling of key chemicals, including the greenhouse gas nitrous oxide (N₂O). To understand the enzymatic basis of N₂O cycling in the ocean, we mined metagenomes to characterize genes encoding bacterial nitrous oxide reductase (NosZ) catalyzing N₂O reduction to N₂. We examined data sets from diverse biomes but focused primarily on those from oxygen minimum zones where N₂O levels are often elevated. With few exceptions, marine nosZ data sets were dominated by 'atypical' clade II gene variants. Atypical nosZ has been associated with low oxygen, enhanced N₂O affinity, and organisms lacking enzymes for complete denitrification, i.e., non-denitrifiers. Atypical nosZ often occurred in metagenome-assembled genomes (MAGs) with nitrate or nitrite respiration genes, although MAGs with genes for complete denitrification were rare. We identified atypical nosZ in several taxa not previously associated with N₂O consumption, in addition to known N₂O-associated groups. The data suggest that marine environments generally select for high N₂O-scavenging ability across diverse taxa and have implications for how N₂O concentration may affect N₂O removal rates.

Spatial Distribution Patterns of Bacterioplankton in the Oxygen Minimum Zone of the Tropical Mexican Pacific

Microbial communities within oxygen minimum zones (OMZs) are crucial drivers of marine biogeochemical cycles; however, we still lack an understanding of how these communities are distributed across an OMZ. We explored vertical (from 5 to 500 m depth) and horizontal (coast to open ocean) distribution of bacterioplankton and its relationships with the main oceanographic conditions in three transects of the tropical Mexican Pacific OMZ. The distribution of the microbial diversity and the main clades changed along the transition from oxygen-rich surface water to the OMZ core, demonstrating the sensitivity of key bacterial groups to deoxygenation. The euphotic zone was dominated by Synechococcales, followed by Flavobacteriales, Verrucomicrobiales, Rhodobacterales, SAR86, and Cellvibrionales, whereas the OMZ core was dominated by SAR11, followed by SAR406, SAR324, SAR202, UBA10353 marine group, Thiomicrospirales and Nitrospinales. The marked environmental gradients along the water column also supported a high potential for niche partitioning among OMZ microorganisms. Additionally, in the OMZ core, bacterial assemblages from the same water mass were more similar to each other than those from another water mass. There were also important differences between coastal and open-ocean communities: Flavobacteriales, Verrucomicrobiales, Rhodobacterales, SAR86, and Cellvibrionales were more abundant in coastal areas, while Synechococcales, SAR406, SAR324, SAR202, UBA10353 marine group, and Thiomicrospirales were more abundant in the open ocean. Our results suggest a biogeographic structure of the bacterioplankton in this OMZ region, with limited community mixing across water masses, except in upwelling events, and little dispersion of the community by currents in the euphotic zone.

Using modern low-oxygen marine ecosystems to understand the nitrogen cycle of the Paleo- and Mesoproterozoic oceans

During the productive Paleoproterozoic (2.4-1.8 Ga) and less productive Mesoproterozoic (1.8-1.0 Ga), the ocean was suboxic to anoxic and multicellular organisms had not yet evolved. Here, we link geologic information about the Proterozoic ocean to microbial processes in modern low-oxygen systems. High iron concentrations and rates of Fe cycling in the Proterozoic are the largest differences from modern oxygen-deficient zones. In anoxic waters, which composed most of the Paleoproterozoic and ~40% of the Mesoproterozoic ocean, nitrogen cycling dominated. Rates of N₂ production by denitrification and anammox were likely linked to sinking organic matter fluxes and in situ primary productivity under anoxic conditions. Additionally autotrophic denitrifiers could have used reduced iron or methane. 50% of the Mesoproterozoic ocean may have been suboxic, promoting nitrification and metal oxidation in the suboxic water and N₂ O and N₂ production by partial and complete denitrification in anoxic zones in organic aggregates. Sulfidic conditions may have composed ~10% of the Mesoproterozoic ocean focused along continental margins. Due to low nitrate concentrations in offshore regions, anammox bacteria likely dominated N₂ production immediately above sulfidic zones, but in coastal regions, higher nitrate concentrations probably promoted complete S-oxidizing autotrophic denitrification at the sulfide interface.

Seasonality and Community Separation of Fungi in a Municipal Wastewater Treatment Plant

Fungi are known to play important roles in pollutant transformation in activated sludge-based wastewater treatment plants (WWTPs). However, the seasonality and distributions of fungal populations in different-sized flocs have still remained largely unknown. In this study, seasonal population dynamics and community separation of fungi in a municipal WWTP across a 1-year period were investigated. We classified all taxa into six categories based on abundances to assess their roles and contributions to the whole community. The results showed that the rare taxa (<0.01%) contributed greatly to species richness (95.27%). Conversely, although low in species diversity, abundant taxa (≥1%) accounted for the majority (89.45%) of the total relative abundance, which suggested that a few core abundant fungi existed in the activated sludge ecosystem. The abundant, conditionally rare, and rare taxa contributed 30.14%, 31.11%, and 38.75%, respectively, to temporal shifts in community structure, and their abundances responded differently to environmental variables, suggesting that these three subcommunities exhibited a large difference in environmental sensitivity. Importantly, the results revealed seasonal dynamics of the whole fungal community and the subcommunities of all the microbial taxon categories, resulting in significant differences in community structures between warm and cold seasons. Furthermore, fungal diversity and the compositions of the whole community and subcommunities differed significantly among flocs of different sizes, which underlined the size-based fungal community separation in activated sludge of WWTPs. The findings of this work improved our understanding of fungal population dynamics and community separation in WWTPs.IMPORTANCE Fungi are important contributors to the various functions of activated sludge in wastewater treatment plants (WWTPs). Unlike previous studies, this work demonstrated the seasonality of the fungal community over a longer time span while it also systematically assessed the contributions of abundant, conditionally rare, and rare taxa to the whole community. Importantly, in the present study, we considered sludge flocs of a certain size range rather than the whole sludge flocs as a community. Our results revealed significant differences in fungal community structure among different-sized flocs, which supported the idea that size-based fungal community segregation is occurring in activated sludge ecosystems. The findings provide new insights into the dynamic changes or distribution of fungi in the bioaggregates of sludge flocs in WWTPs.

Diverse nitrogen cycling pathways across a marine oxygen gradient indicate nitrogen loss coupled to chemoautotrophic activity

Genetic markers and geochemical assays of microbial nitrogen cycling processes, including autotrophic and heterotrophic denitrification, anammox, ammonia oxidation, and nitrite oxidation, were examined across the oxycline, suboxic, and anoxic zones of the Cariaco Basin, Venezuela. Ammonia and nitrite oxidation genes were expressed through the entire gradient. Transcripts associated with autotrophic and heterotrophic denitrifiers were mostly confined to the suboxic zone and below but were also present in particles in the oxycline. Anammox genes and transcripts were detected over a narrow depth range near the bottom of the suboxic zone and coincided with secondary NO₂ ^- maxima and available NH₄ ⁺ . Dissolved inorganic nitrogen (DIN) amendment incubations and comparisons between our sampling campaigns suggested that denitrifier activity may be closely coupled with NO₃ ^- availability. Expression of denitrification genes at depths of high rates of chemoautotrophic carbon fixation and phylogenetic analyses of nitrogen cycling genes and transcripts indicated a diverse array of denitrifiers, including chemoautotrophs capable of using NO₃ ^- to oxidize reduced sulfur species. Thus, results suggest that the Cariaco Basin nitrogen cycle is influenced by autotrophic carbon cycling in addition to organic matter oxidation and anammox.

Colonization kinetics and implantation follow-up of the sewage microbiome in an urban wastewater treatment plant

The Seine-Morée wastewater treatment plant (SM_WWTP), with a capacity of 100,000 population-equivalents, was fed with raw domestic wastewater during all of its start-up phase. Its microbiome resulted from the spontaneous evolution of wastewater-borne microorganisms. This rare opportunity allowed us to analyze the sequential microbiota colonization and implantation follow up during the start-up phase of this WWTP by means of regular sampling carried out over 8 months until the establishment of a stable and functional ecosystem. During the study, biological nitrification–denitrification and dephosphatation occurred 68 days after the start-up of the WWTP, followed by flocs decantation 91 days later. High throughput sequencing of 18S and 16S rRNA genes was performed using Illumina's MiSeq and PGM Ion Torrent platforms respectively, generating 584,647 16S and 521,031 18S high-quality sequence rDNA reads. Analyses of 16S and 18S rDNA datasets show three colonization phases occurring concomitantly with nitrification, dephosphatation and floc development processes. Thus, we could define three microbiota profiles that sequentially colonized the SM_WWTP: the early colonizers, the late colonizers and the continuous spectrum population. Shannon and inverse Simpson diversity indices indicate that the highest microbiota diversity was reached at days 133 and 82 for prokaryotes and eukaryotes respectively; after that, the structure and complexity of the wastewater microbiome reached its functional stability. This study demonstrates that physicochemical parameters and microbial metabolic interactions are the main forces shaping microbial community structure, gradually building up and maintaining a functionally stable microbial ecosystem.

Aggregate size-dependence of bacterial community separation in a full-scale membrane bioreactor plant

The size of bacterial aggregates can determine both nutrient removal and sludge/water separation in activated sludge processes. In this study, the bacterial community structures and network associations of different sized aggregates obtained from a full-scale membrane bioreactor plant over a one-year period were investigated. Our results showed that biodiversity of larger sized aggregates was significantly higher than that of smaller ones and that the bacterial compositions of different sized aggregates differed significantly from each other. Bacteria related to nutrient removal (e.g. denitrification, hydrolysis and fermentation) were found to be significantly more abundant in larger aggregates than smaller ones. Network analysis revealed significant difference in species-species interactions, topological roles of individual OTU and keystone populations among different sized aggregates. Moreover, the occurrence of keystone OTUs affiliated with denitrifiers (Thermomonas) in networks of large and medium aggregates may suggest that denitrification influences bacterial interactions in large and medium aggregates. Overall, our results indicate the aggregates size-dependence of bacterial community separation in activated sludge. The current findings not only can provide guidance for process design and operation optimization, but also highlight the necessity for paying more attentions to the aggregate-based community, metabolic function and gene expression of activated sludge in future studies.

Microbial Diversity and Metabolic Potential in the Stratified Sansha Yongle Blue Hole in the South China Sea

The Sansha Yongle Blue Hole is the world’s deepest (301 m) underwater cave and has a sharp redox gradient, with oligotrophic, anoxic, and sulfidic bottom seawater. In order to discover the microbial communities and their special biogeochemical pathways in the blue hole, we analyzed the 16S ribosomal RNA amplicons and metagenomes of microbials from seawater depths with prominent physical, chemical, and biological features. Redundancy analysis showed that dissolved oxygen was the most important factor affecting the microbial assemblages of the blue hole and surrounding open sea waters, and significantly explained 44.7% of the total variation, followed by silicate, temperature, sulfide, ammonium, methane, nitrous oxide, nitrate, dissolved organic carbon, salinity, particulate organic carbon, and chlorophyll a. We identified a bloom of Alteromonas (34.9%) at the primary nitrite maximum occurring in close proximity to the chlorophyll a peak in the blue hole. Genomic potential for nitrate reduction of Alteromonas might contribute to this maximum under oxygen decrease. Genes that would allow for aerobic ammonium oxidation, complete denitrification, and sulfur-oxidization were enriched at nitrate/nitrite-sulfide transition zone (90 and 100 m) of the blue hole, but not anammox pathways. Moreover, γ-Proteobacterial clade SUP05, ε-Proteobacterial genera Sulfurimonas and Arcobacter, and Chlorobi harbored genes for sulfur-driven denitrification process that mediated nitrogen loss and sulfide removal. In the anoxic bottom seawater (100-300 m), high levels of sulfate reducers and dissimilatory sulfite reductase gene (dsrA) potentially created a sulfidic zone of ~200 m thickness. Our findings suggest that in the oligotrophic Sansha Yongle Blue Hole, O₂ deficiency promotes nitrogen- and sulfur-cycling processes mediated by metabolically versatile microbials.

Nitrifier adaptation to low energy flux controls inventory of reduced nitrogen in the dark ocean

Ammonia oxidation to nitrite and its subsequent oxidation to nitrate provides energy to the two populations of nitrifying chemoautotrophs in the energy-starved dark ocean, driving a coupling between reduced inorganic nitrogen (N) pools and production of new organic carbon (C) in the dark ocean. However, the relationship between the flux of new C production and the fluxes of N of the two steps of oxidation remains unclear. Here, we show that, despite orders-of-magnitude difference in cell abundances between ammonia oxidizers and nitrite oxidizers, the two populations sustain similar bulk N-oxidation rates throughout the deep waters with similarly high affinities for ammonia and nitrite under increasing substrate limitation, thus maintaining overall homeostasis in the oceanic nitrification pathway. Our observations confirm the theoretical predictions of a redox-informed ecosystem model. Using balances from this model, we suggest that consistently low ammonia and nitrite concentrations are maintained when the two populations have similarly high substrate affinities and their loss rates are proportional to their maximum growth rates. The stoichiometric relations between the fluxes of C and N indicate a threefold to fourfold higher C-fixation efficiency per mole of N oxidized by ammonia oxidizers compared to nitrite oxidizers due to nearly identical apparent energetic requirements for C fixation of the two populations. We estimate that the rate of chemoautotrophic C fixation amounts to ∼1 × 1013 to ∼2 × 1013 mol of C per year globally through the flux of ∼1 × 1014 to ∼2 × 1014 mol of N per year of the two steps of oxidation throughout the dark ocean.

Uncultured Nitrospina-like species are major nitrite oxidizing bacteria in oxygen minimum zones

Oxygen minimum zones (OMZs) are marine regions where O2 is undetectable at intermediate depths. Within OMZs, the oxygen-depleted zone (ODZ) induces anaerobic microbial processes that lead to fixed nitrogen loss via denitrification and anammox. Surprisingly, nitrite oxidation is also detected in ODZs, although all known marine nitrite oxidizers (mainly Nitrospina) are aerobes. We used metagenomic binning to construct metagenome-assembled genomes (MAGs) of nitrite oxidizers from OMZs. These MAGs represent two novel Nitrospina-like species, both of which differed from all known Nitrospina species, including cultured species and published MAGs. Relative abundances of different Nitrospina genotypes in OMZ and non-OMZ seawaters were estimated by mapping metagenomic reads to newly constructed MAGs and published high-quality genomes of members from the Nitrospinae phylum. The two novel species were present in all major OMZs and were more abundant inside ODZs, which is consistent with the detection of higher nitrite oxidation rates in ODZs than in oxic seawaters and suggests novel adaptations to anoxic environments. The detection of a large number of unclassified nitrite oxidoreductase genes in the dataset implies that the phylogenetic diversity of nitrite oxidizers is greater than previously thought.

Niche Differentiation of Aerobic and Anaerobic Ammonia Oxidizers in a High Latitude Deep Oxygen Minimum Zone

To elucidate the potential for nitrification and denitrification processes in a high latitude deep oxygen minimum zone (OMZ) we determined the abundance and community composition of the main microbial players in the aerobic and anaerobic (anammox) ammonium oxidation and denitrification processes in the Gulf of Alaska throughout the water column. Within the dominant bacterial groups, Flavobacterales, Rhodobacterales, Actinomarinales, and SAR86 were more abundant in epipelagic waters and decreased with depth, whereas SAR11, SAR324, Marinimicrobia, and Thiomicrospirales increased their contribution to the bacterial community with depth. Nitrosopumilaceae also increased with depth and dominated the OMZ and bathypelagic archaeal communities. Euryarchaeota Marine Group II exhibited an opposite depth pattern to Nitrosopumilaceae, whereas Marine Group III and Woesearchaeota were more abundant in the bathypelagic realm. Candidatus Brocadia contributed 70-100% of the anammox bacterial community throughout the water column. Archaeal ammonia oxidizers (AOA) dominated the microbial community involved in the nitrogen cycle. Two AOA ecotypes, the high ammonia (HAC) and low ammonia (LAC)-AOA, characterized by distinct genes for aerobic ammonia oxidation (amoA) and for denitrification (nirK), exhibited a distinct distribution pattern related to depth and ammonia concentrations. HAC-AOA dominated in epipelagic (80.5 ± 28.3% of total AOA) oxygenated and ammonia-rich waters, and LAC-AOA dominated in the OMZ (90.9 ± 5.1%) and bathypelagic waters (85.5 ± 13.5%), characterized by lower oxygen and ammonia concentrations. Bacterial denitrifiers (3.7 ± 6.9 bacterial nirK gene mL-1) and anaerobic ammonia oxidizers (78 ± 322 anammox 16S rRNA genes L-1) were low in abundance under the oxygen conditions in the Gulf of Alaska throughout the water column. The widespread distribution of bacterial denitrifiers and anaerobic ammonia oxidizers in low abundances reveals a reservoir of genetic and metabolic potential ready to colonize the environment under the predicted increase of OMZs in the ocean. Taken together, our results reinforce the niche partitioning of archaeal ammonia oxidizers based on their distinct metabolic characteristics resulting in the dominance of LAC-AOA in a high latitude deep OMZ. Considering the different ecological roles and functions of the two archaeal ecotypes, the expansion of the zones dominated by the LAC-ecotype might have implications for the nitrogen cycle in the future ocean.

Microbial niches in marine oxygen minimum zones

In the ocean's major oxygen minimum zones (OMZs), oxygen is effectively absent from sea water and life is dominated by microorganisms that use chemicals other than oxygen for respiration. Recent studies that combine advanced genomic and chemical detection methods are delineating the different metabolic niches that microorganisms can occupy in OMZs. Understanding these niches, the microorganisms that inhabit them, and their influence on marine biogeochemical cycles is crucial as OMZs expand with increasing seawater temperatures.

Geomicrobiology of the carbon, nitrogen and sulphur cycles in Powell Lake: a permanently stratified water column containing ancient seawater

We present the first geomicrobiological characterization of the meromictic water column of Powell Lake (British Columbia, Canada), a former fjord, which has been stably stratified since the last glacial period. Its deepest layers (300-350 m) retain isolated, relict seawater from that period. Fine-scale vertical profiling of the water chemistry and microbial communities allowed subdivision of the water column into distinct geomicrobiological zones. These zones were further characterized by phylogenetic and functional marker genes from amplicon and shotgun metagenome sequencing. Binning of metagenomic reads allowed the linkage of function to specific taxonomic groups. Statistical analyses (analysis of similarities, Bray-Curtis similarity) confirmed that the microbial community structure followed closely the geochemical zonation. Yet, our characterization of the genetic potential relevant to carbon, nitrogen and sulphur cycling of each zone revealed unexpected features, including potential for facultative anaerobic methylotrophy, nitrogen fixation despite high ammonium concentrations and potential micro-aerobic nitrifiers within the chemocline. At the oxic-suboxic interface, facultative anaerobic potential was found in the widespread freshwater lineage acI (Actinobacteria), suggesting intriguing ecophysiological similarities to the marine SAR11. Evolutionary divergent lineages among diverse phyla were identified in the ancient seawater zone and may indicate novel adaptations to this unusual environment.

Utilization of urea and cyanate in waters overlying and within the eastern tropical north Pacific oxygen deficient zone

In marine oxygen deficient zones (ODZs), which contribute up to half of marine N loss, microbes use nitrogen (N) for assimilatory and dissimilatory processes. Here, we examine N utilization above and within the ODZ of the Eastern Tropical North Pacific Ocean, focusing on distribution, uptake and genes for the utilization of two simple organic N compounds, urea and cyanate. Ammonium, urea and cyanate concentrations generally peaked in the oxycline while uptake rates were highest in the surface. Within the ODZ, concentrations were lower, but urea N and C and cyanate C were taken up. All identified autotrophs had an N assimilation pathway that did not require external ammonium: ODZ Prochlorococcus possessed genes to assimilate nitrate, nitrite and urea; nitrite oxidizers (Nitrospina) possessed genes to assimilate nitrite, urea and cyanate; anammox bacteria (Scalindua) possessed genes to utilize cyanate; and ammonia-oxidizing Thaumarchaeota possessed genes to utilize urea. Urease genes were present in 20% of microbes, including SAR11, suggesting the urea utilization capacity was widespread. In the ODZ core, cyanate genes were largely (∼95%) associated with Scalindua, suggesting that, within this ODZ, cyanate N is primarily used for N loss via anammox (cyanammox), and that anammox does not require ammonium for N loss.

The water depth-dependent co-occurrence patterns of marine bacteria in shallow and dynamic Southern Coast, Korea

To investigate the interactions between bacterial species in relation to the biotic and abiotic environmental fluctuations, free-living (FL), nanoparticle-associated (NP), and microparticle-associated (MP) bacterial community compositions (BCCs) were analyzed. A total of 267 samples were collected from July to December 2016 in the dynamic and shallow southern coastal water of Korea. The variations in BCC mostly depended on planktonic size fraction. Network analysis revealed water depth-dependent co-occurrence patterns of coastal bacterial communities. Higher interspecies connectivity was observed within FL bacteria than NP/MP bacteria, suggesting that FL bacteria with a streamlined genome may need other bacterial metabolites for survival, while the NP/MP copiotrophs may have the self-supporting capacity to produce the vital nutrients. The analysis of topological roles of individual OTUs in the network revealed that several groups of metabolically versatile bacteria (the marine Roseobacters, Flavobacteriales, Desulfobacterales, and SAR406 clade) acted as module hubs in different water depth. In conclusion, interspecies interactions dominated in FL bacteria, compared to NP and MP bacteria; modular structures of bacterial communities and keystone species strongly depended on the water depth-derived environmental factors. Furthermore, the multifunctional, versatile FL bacteria could play pivotal roles in dynamic shallow coastal ecosystems.

Complete arsenic-based respiratory cycle in the marine microbial communities of pelagic oxygen-deficient zones

Microbial capacity to metabolize arsenic is ancient, arising in response to its pervasive presence in the environment, which was largely in the form of As(III) in the early anoxic ocean. Many biological arsenic transformations are aimed at mitigating toxicity; however, some microorganisms can respire compounds of this redox-sensitive element to reap energetic gains. In several modern anoxic marine systems concentrations of As(V) are higher relative to As(III) than what would be expected from the thermodynamic equilibrium, but the mechanism for this discrepancy has remained unknown. Here we present evidence of a complete respiratory arsenic cycle, consisting of dissimilatory As(V) reduction and chemoautotrophic As(III) oxidation, in the pelagic ocean. We identified the presence of genes encoding both subunits of the respiratory arsenite oxidase AioA and the dissimilatory arsenate reductase ArrA in the Eastern Tropical North Pacific (ETNP) oxygen-deficient zone (ODZ). The presence of the dissimilatory arsenate reductase gene arrA was enriched on large particles (>30 um), similar to the forward bacterial dsrA gene of sulfate-reducing bacteria, which is involved in the cryptic cycling of sulfur in ODZs. Arsenic respiratory genes were expressed in metatranscriptomic libraries from the ETNP and the Eastern Tropical South Pacific (ETSP) ODZ, indicating arsenotrophy is a metabolic pathway actively utilized in anoxic marine water columns. Together these results suggest arsenic-based metabolisms support organic matter production and impact nitrogen biogeochemical cycling in modern oceans. In early anoxic oceans, especially during periods of high marine arsenic concentrations, they may have played a much larger role.

Changes in Carbon Oxidation State of Metagenomes Along Geochemical Redox Gradients

There is widespread interest in how geochemistry affects the genomic makeup of microbial communities, but the possible impacts of oxidation-reduction (redox) conditions on the chemical composition of biomacromolecules remain largely unexplored. Here we document systematic changes in the carbon oxidation state, a metric derived from the chemical formulas of biomacromolecular sequences, using published metagenomic and metatranscriptomic datasets from 18 studies representing different marine and terrestrial environments. We find that the carbon oxidation states of DNA, as well as proteins inferred from coding sequences, follow geochemical redox gradients associated with mixing and cooling of hot spring fluids in Yellowstone National Park (USA) and submarine hydrothermal fluids. Thermodynamic calculations provide independent predictions for the environmental shaping of the gene and protein composition of microbial communities in these systems. On the other hand, the carbon oxidation state of DNA is negatively correlated with oxygen concentration in marine oxygen minimum zones. In this case, a thermodynamic model is not viable, but the low carbon oxidation state of DNA near the ocean surface reflects a low GC content, which can be attributed to genome reduction in organisms adapted to low-nutrient conditions. We also present evidence for a depth-dependent increase of oxidation state at the species level, which might be associated with alteration of DNA through horizontal gene transfer and/or selective degradation of relatively reduced (AT-rich) extracellular DNA by heterotrophic bacteria. Sediments exhibit even more complex behavior, where carbon oxidation state minimizes near the sulfate-methane transition zone and rises again at depth; markedly higher oxidation states are also associated with older freshwater-dominated sediments in the Baltic Sea that are enriched in iron oxides and have low organic carbon. This geobiochemical study of carbon oxidation state reveals a new aspect of environmental information in metagenomic sequences, and provides a reference frame for future studies that may use ancient DNA sequences as a paleoredox indicator.

Seaweed Loads Cause Stronger Bacterial Community Shifts in Coastal Lagoon Sediments Than Nutrient Loads

The input of nutrients from anthropogenic sources is the leading cause of coastal eutrophication and is usually coupled with algal/seaweed blooms. Effects may be magnified in semi-enclosed systems, such as highly productive coastal lagoon ecosystems. Eutrophication and seaweed blooms can lead to ecosystem disruption. Previous studies have considered only one of these factors, disregarding possible interactive effects and the effect of the blooming species' identity on sediment bacterial communities. We tested the effect of experimental nutrient loading and two common blooming seaweeds (Ulva rigida and Gracilaria vermiculophylla) in coastal lagoon sediments, on the structure of bacterial communities (using 16S rRNA amplicon sequencing) and corresponding putative functional potential (using PiCRUSt). At the Operational Taxonomic Unit (OTU) level, the addition of nutrients reduced bacterial community α-diversity and decreased the abundance of sulfate reducers (Desulfobacterales) compared to sulfur oxidizers/denitrifiers (Chromatiales and Campylobacterales), whereas this was not the case at the order level. Seaweed addition did not change bacterial α-diversity and the effect on community structure depended on the taxonomic level considered. The addition of Gracilaria increased the abundance of orders and OTUs involved in sulfate reduction and organic matter decomposition (Desulfobacterales, Bacteroidales, and Clostridiales, respectively), an effect which was also detected when only Ulva was added. Nutrients and the seaweeds combined only interacted for Ulva and nutrients, which increased known sulfide oxidizers and denitrifiers (order Campylobacterales). Seaweed enrichment affected putative functional profiles; a stronger increase of sulfur cycling KEGG pathways was assigned to nutrient-disturbed sediments, particularly with the seaweeds and especially Ulva. In contrast, nitrogen and sulfur cycle pathways showed a higher abundance of genes related to dissimilatory nitrate reduction to ammonium (DNRA) in Ulva+nutrients treatments. However, the other seaweed treatments increased the nitrogen fixation genes. Thiosulfate reduction, performed by sulfate-reducing bacteria, increased in seaweed treatments except when Ulva was combined with nutrients. In conclusion, the in situ addition of nutrients and the seaweeds to intertidal sediments affected the bacterial communities differently and independently. The predicted functional profile suggests a shift in relative abundances of putative pathways for nitrogen and sulfur cycles, in line with the taxonomic changes of the bacterial communities.

Depth-Dependent Environmental Drivers of Microbial Plankton Community Structure in the Northern Gulf of Mexico

The Gulf of Mexico (GoM) is a dynamic marine ecosystem influenced by multiple natural and anthropogenic processes and inputs, such as the intrusion of warm oligotrophic water via the Loop Current, freshwater and nutrient input by the Mississippi River, and hydrocarbon inputs via natural seeps and industrial spills. Microbial plankton communities are important to pelagic food webs including in the GoM but understanding the drivers of the natural dynamics of these passively distributed microorganisms can be challenging in such a large and heterogeneous system. As part of the DEEPEND consortium, we applied high throughput 16S rRNA sequencing to investigate the spatial and temporal dynamics of pelagic microbial plankton related to several environmental conditions during two offshore cruises in 2015. Our results show dramatic community shifts across depths, especially between photic and aphotic zones. Though we only have two time points within a single year, observed temporal shifts in microbial plankton communities were restricted to the seasonally influenced epipelagic zone (0-200 m), and appear mainly driven by changes in temperature. Environmental selection in microbial plankton communities was depth-specific, with variables such as turbidity, salinity, and abundance of photosynthetic taxa strongly correlating with community structure in the epipelagic zone, while variables such as oxygen and specific nutrient concentrations were correlated with community structure at deeper depths.

Particle-Attached and Free-Living Archaeal Communities in the Benthic Boundary Layer of the Mariana Trench

The benthic boundary layer (BBL) is the part of the water column that is situated near to the sediment surface, where active oceanic biogeochemical cycling occurs. Archaea play an important role in mediating this cycling, however, their composition and diversity in the BBL remain largely unknown. We investigated the community composition and abundance of both particle-attached (PA) and free-living (FL) archaea in the BBL on the slopes of the Mariana Trench using Illumina sequencing and quantitative PCR (qPCR), at both the DNA and RNA levels. Our results showed that Thaumarchaeota (>90%) and Woesearchaeota (1–10%) dominated in all the BBL samples, and that the former was composed mainly of Marine Group I (MGI). A clear separation of PA and FL samples was observed, and they showed a high level of similarity to the subsurface sediments and the water column, respectively. No significant differences were detected in the archaeal communities located in the southern and northern slopes of the Mariana Trench, or between the levels of DNA and RNA. However, lower RNA/DNA ratios (estimated by qPCR) were found in the PA samples than in the FL samples, indicating higher transcriptional activities in the FL fractions. A distinct archaeal community structure was found in the middle of the trench when compared with samples collected at the same depth at other stations along the trench slopes. This indicates that a dynamic deep current might affect the distribution of organic matter on the slopes. Our study provides direct information regarding the archaeal communities in the BBL of the Mariana Trench. We suggest that this might promote further exploration of the ecological roles and microbial processes of such communities located in deep-sea ecosystems.