Marine oxygen minimum zones (OMZs) support diverse microbial communities with roles in major elemental cycles. It is unclear how the taxonomic composition and metabolism of OMZ microorganisms vary between particle-associated and free-living size fractions. We used amplicon (16S rRNA gene) and shotgun metagenome sequencing to compare microbial communities from large (41.6 lm) and small (0.2-1.6 lm) filter size fractions along a depth gradient in the OMZ off Chile. Despite steep vertical redox gradients, size fraction was a significantly stronger predictor of community composition compared to depth. Phylogenetic diversity showed contrasting patterns, decreasing towards the anoxic OMZ core in the small size fraction, but exhibiting maximal values at these depths within the larger size fraction. Fraction-specific distributions were evident for key OMZ taxa, including anammox planctomycetes, whose coding sequences were enriched up to threefold in the 0.2-1.6 lm community. Functional gene composition also differed between fractions, with the 41.6 lm community significantly enriched in genes mediating social interactions, including motility, adhesion, cell-to-cell transfer, antibiotic resistance and mobile element activity. Prokaryotic transposase genes were three to six fold more abundant in this fraction, comprising up to 2% of protein-coding sequences, suggesting that particle surfaces may act as hotbeds for transposition-based genome changes in marine microbes. Genes for nitric and nitrous oxide reduction were also more abundant (three to seven fold) in the larger size fraction, suggesting microniche partitioning of key denitrification steps. These results highlight an important role for surface attachment in shaping community metabolic potential and genome content in OMZ microorganisms.

The genetic composition of marine microbial communities varies at the microscale between particleassociated (PA; 41.6 μm) and free-living (FL; 0.2-1.6 μm) niches. It remains unclear, however, how metabolic activities differ between PA and FL fractions. We combined rate measurements with metatranscriptomics to quantify PA and FL microbial activity in the oxygen minimum zone (OMZ) of the Eastern Tropical North Pacific, focusing on dissimilatory processes of the nitrogen (N) cycle. Bacterial gene counts were 8-to 15-fold higher in the FL compared with the PA fraction. However, rates of all measured N cycle processes, excluding ammonia oxidation, declined significantly following particle (41.6 μm) removal. Without particles, rates of nitrate reduction to nitrite (1.5-9.4 nM N d − 1 ) fell to zero and N 2 production by denitrification (0.5-1.7 nM N d − 1 ) and anammox (0.3-1.9 nM N d − 1 ) declined by 53-85%. The proportional representation of major microbial taxa and N cycle gene transcripts in metatranscriptomes followed fraction-specific trends. Transcripts encoding nitrate reductase were uniform among PA and FL fractions, whereas anammox-associated transcripts were proportionately enriched up to 15-fold in the FL fraction. In contrast, transcripts encoding enzymes for N 2 O and N 2 production by denitrification were enriched up to 28-fold in PA samples. These patterns suggest that the majority of N cycle activity, excluding N 2 O and N 2 production by denitrification, is confined to a FL majority that is critically dependent on access to particles, likely as a source of organic carbon and inorganic N. Variable particle distributions may drive heterogeneity in N cycle activity and gene expression in OMZs.

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.

Iron (Fe) and copper (Cu) are essential cofactors for microbial metalloenzymes, but little is known about the metalloenyzme inventory of anaerobic marine microbial communities despite their importance to the nitrogen cycle. We compared dissolved O2, NO3−, NO2−, Fe and Cu concentrations with nucleic acid sequences encoding Fe and Cu-binding proteins in 21 metagenomes and 9 metatranscriptomes from Eastern Tropical North and South Pacific oxygen minimum zones and 7 metagenomes from the Bermuda Atlantic Time-series Station. Dissolved Fe concentrations increased sharply at upper oxic-anoxic transition zones, with the highest Fe:Cu molar ratio (1.8) occurring at the anoxic core of the Eastern Tropical North Pacific oxygen minimum zone and matching the predicted maximum ratio based on data from diverse ocean sites. The relative abundance of genes encoding Fe-binding proteins was negatively correlated with O2, driven by significant increases in genes encoding Fe-proteins involved in dissimilatory nitrogen metabolisms under anoxia. Transcripts encoding cytochrome c oxidase, the Fe- and Cu-containing terminal reductase in aerobic respiration, were positively correlated with O2 content. A comparison of the taxonomy of genes encoding Fe- and Cu-binding vs. bulk proteins in OMZs revealed that Planctomycetes represented a higher percentage of Fe genes while Thaumarchaeota represented a higher percentage of Cu genes, particularly at oxyclines. These results are broadly consistent with higher relative abundance of genes encoding Fe-proteins in the genome of a marine planctomycete vs. higher relative abundance of genes encoding Cu-proteins in the genome of a marine thaumarchaeote. These findings highlight the importance of metalloenzymes for microbial processes in oxygen minimum zones and suggest preferential Cu use in oxic habitats with Cu > Fe vs. preferential Fe use in anoxic niches with Fe > Cu.

Molecular surveys are revealing diverse eukaryotic assemblages in oxygen-limited ocean waters. These communities may play pivotal ecological roles through autotrophy, feeding, and a wide range of symbiotic associations with prokaryotes. We used 18S rRNA gene sequencing to provide the first snapshot of pelagic microeukaryotic community structure in two cellular size fractions (0.2–1.6 μm, >1.6 μm) from seven depths through the anoxic oxygen minimum zone (OMZ) off northern Chile. Sequencing of >154,000 amplicons revealed contrasting patterns of phylogenetic diversity across size fractions and depths. Protist and total eukaryote diversity in the >1.6 μm fraction peaked at the chlorophyll maximum in the upper photic zone before declining by ~50% in the OMZ. In contrast, diversity in the 0.2–1.6 μm fraction, though also elevated in the upper photic zone, increased four-fold from the lower oxycline to a maximum at the anoxic OMZ core. Dinoflagellates of the Dinophyceae and endosymbiotic Syndiniales clades dominated the protist assemblage at all depths (~40–70% of sequences). Other protist groups varied with depth, with the anoxic zone community of the larger size fraction enriched in euglenozoan flagellates and acantharean radiolarians (up to 18 and 40% of all sequences, respectively). The OMZ 0.2–1.6 μm fraction was dominated (11–99%) by Syndiniales, which exhibited depth-specific variation in composition and total richness despite uniform oxygen conditions. Metazoan sequences, though confined primarily to the 1.6 μm fraction above the OMZ, were also detected within the anoxic zone where groups such as copepods increased in abundance relative to the oxycline and upper OMZ. These data, compared to those from other low-oxygen sites, reveal variation in OMZ microeukaryote composition, helping to identify clades with potential adaptations to oxygen-depletion.

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.

Endosymbioses between animals and chemoautotrophic bacteria are ubiquitous at hydrothermal vents. These environments are distinguished by high physico-chemical variability, yet we know little about how these symbioses respond to environmental fluctuations. We therefore examined how the γ-proteobacterial symbionts of the vent snail Ifremeria nautilei respond to changes in sulfur geochemistry. Via shipboard high-pressure incubations, we subjected snails to 105 μM hydrogen sulfide (LS), 350 μM hydrogen sulfide (HS), 300 μM thiosulfate (TS) and seawater without any added inorganic electron donor (ND). While transcript levels of sulfur oxidation genes were largely consistent across treatments, HS and TS treatments stimulated genes for denitrification, nitrogen assimilation, and CO2 fixation, coincident with previously reported enhanced rates of inorganic carbon incorporation and sulfur oxidation in these treatments. Transcripts for genes mediating oxidative damage were enriched in the ND and LS treatments, potentially due to a reduction in O2 scavenging when electron donors were scarce. Oxidative TCA cycle gene transcripts were also more abundant in ND and LS treatments, suggesting that I. nautilei symbionts may be mixotrophic when inorganic electron donors are limiting. These data reveal the extent to which I. nautilei symbionts respond to changes in sulfur concentration and species, and, interpreted alongside coupled biochemical metabolic rates, identify gene targets whose expression patterns may be predictive of holobiont physiology in environmental samples.

Bacterial communities colonized on submerged substrata are recognized as a key factor in the formation of complex biofouling phenomenon in the marine environment. Despite massive maritime activities and a large industrial sector in the nearshore of the Laccadive Sea, studies describing pioneer bacterial colonizers and community succession during the early-stage biofilm are scarce. We investigated the biofilm-forming bacterial community succession on three substrata viz. stainless steel, high-density polyethylene, and titanium over 15 days of immersion in the seawater intake area of a power plant, located in the southern coastal region of India. The bacterial community composition of biofilms and peripheral seawater were analyzed by Illumina MiSeq sequenced 16S rRNA gene amplicons. The obtained metataxonomic results indicated a profound influence of temporal succession over substrate type on the early-stage biofilm-forming microbiota. Bacterial communities showed vivid temporal dynamics that involved variations in abundant bacterial groups. The proportion of dominant phyla viz. Proteobacteria decreased over biofilm succession days, while Bacteroidetes increased, suggesting their role as initial and late colonizers, respectively. A rapid fluctuation in the proportion of two bacterial orders viz. Alteromonadales and Vibrionales were observed throughout the successional stages. LEfSe analysis identified specific bacterial groups at all stages of biofilm development, whereas no substrata type-specific groups were observed. Furthermore, the results of PCoA and UPGMA hierarchical clustering demonstrated that the biofilm-forming community varied considerably from the planktonic community. Phylum Proteobacteria preponderated the biofilm-forming community, while the Bacteroidetes, Cyanobacteria, and Actinobacteria dominated the planktonic community. Overall, our results refute the common assumption that substrate material has a decisive impact on biofilm formation; rather, it portrayed that the temporal succession overshadowed the influence of the substrate material. Our findings provide a scientific understanding of the factors shaping initial biofilm development in the marine environment and will help in designing efficient site-specific anti-biofouling strategies.