- 1University of Maryland Center for Environmental Science Horn Point Laboratory, Cambridge, MD, United States
- 2Center for Dark Energy Biosphere Investigations, University of Southern California, Los Angeles, CA, United States
- 3Branchpoint Sciences, Los Angeles, CA, United States
- 4Department of Oceanography, Texas A&M University, College Station, TX, United States
Introduction: Marine particles form in the ocean surface sink through the water column into the deep ocean, sequestering carbon. Microorganisms inhabit and consume carbon in these particles. The East Pacific Rise (EPR) harbors both an Oxygen Deficient Zone (ODZ) and a non-buoyant plume region formed from hydrothermal vents located on the ocean floor, allowing us to explore relationships between microbial community and particle size between a range of environments.
Methods: In this study, we quantified microbial diversity using a fractionation method that separated particles into seven fine scale fractions (0.2–1.2, 1.2–5, 5–20, 20–53, 53–180,180–500, >500 μm), and included a spike-in standard for sequencing the 16S rRNA gene. Size fractionated organic carbon into the same fractions enabled the calculation of bacterial 16S rRNA copies per μg C and per liter.
Results: There was a large increase in the bacterial 16S rRNA copies/ug C and copies/L on particles >180 μm between the upper water column and the deep water column. Though the total concentration of organic C in particles decreased in the deep water column, the density of bacteria on large particles increased at depth. The microbial community varied statistically significantly as a function of particle size and depth. Quantitative abundance estimates found that ostensibly obligate free-living microbes, such as SAR11 and Thaumarcheota, were more abundant in the free-living fraction but also common and abundant in the particulate size fractions. Conversely, ostensibly obligate particle attached bacteria such as members of Bacteroidetes and Planctomycetes, while most abundant on particles, were also present in the free living fraction. Total bacterial abundance, and the abundance of many taxonomic groups, increased in the ODZ region, particularly in the free-living fraction. Contrastingly, in the non-buoyant plume, there were highly abundant bacteria in the 5–20 and 20–53 μm fractions but reduced bacteria present in the 53–180 and 180–500 μm fractions.
Conclusion: Quantitative examination of microbial communities highlights the distribution of microbial taxa unburdened by compositional effects. These data are congruent with existing models which suggest high levels of exchange between particle-attached and free-living communities.
Introduction
Photosynthesis in the upper ocean is responsible for about 50% of the primary production on earth (Guidi et al., 2009). A fraction of the CO2 incorporated into this organic matter sinks to the deep ocean, removing it from the atmosphere for several centuries, and a small fraction of that sinking carbon is entrained in the sea floor removing that carbon from the atmosphere for millions of years (Le Moigne, 2019). The sinking of organic matter into deep waters constitutes a process known as the gravitational biological carbon pump (Le Moigne, 2019). Ocean regions are connected vertically by these marine particles that transport microorganisms from the surface towards the deep (Mestre et al., 2018).
Particulate Organic Carbon (POC) forms through several mechanical and biological processes including zooplankton ingestion/excretion and molting and aggregation of phytoplankton or of smaller particles into marine snow (Simon et al., 2002). Different members of the phytoplankton group contribute differently to the POC export to deep ocean; for example, modeling of carbon export has shown that diatoms could be contributing 40% of carbon export whereas coccolithophores contribute 10% (Jin et al., 2006). Small eukaryotic algae sink in fecal pellets while large eukaryotic algae can sink as single cells (Durkin et al., 2022). Picocyanobacteria aggregates can also contribute to export, but by forming small particles that are aggregated into larger particles below the mixed layer (Cavan et al., 2018; Fuchsman et al., 2019a). The size, and chemical composition of marine particles will depend on the combination of factors that contributed to their formation.
However, not all particles sink; many particles smaller than <53 μm are suspended (Clegg and Whitfield, 1990). Chemical data indicates that suspended and sinking sediment trap material have different isotopic, lipid, pigment and amino acid compositions (Abramson et al., 2010; Altabet, 1988; Wakeham and Canuel, 1988). Older microbial community fingerprinting techniques indicated that communities on sinking particles from sediment traps were significantly different from communities on suspended particles at depth, with trap material more reflective of surface waters (Tamelander, 2013). Recently, amplicon sequencing also indicated differences in microbial communities between suspended and sinking material (Stephens et al., 2024). However, there is connectivity between sinking and suspended particles with disaggregation of sinking particles forming suspended particles (Alldredge, 1998) and aggregation of suspended particles to form sinking particles (Cavan et al., 2018). Additionally, some small particles sink. At the Bermuda Ocean Time series station, in the upper mesopelagic, 39 ± 20% of carbon sinking flux was in 11–63 μm particles (Durkin et al., 2015). This is because density is actually more important to sinking rate than size.
Marine particles can be colonized by microbial communities that use the organic matter for microbial respiration. Metagenomic data from the Hawaii Ocean Time Series indicate that sinking particles are enriched in bacteria that have extracellular peptidases, which degrade proteins, and carbohydrate degrading enzymes as well as large numbers of ABC transporters (Leu et al., 2022). In the absence of oxygen, organic matter can be respired via anaerobic respiration such as denitrification, nitrate reduction or sulfate reduction (Fuchsman et al., 2017; Lam et al., 2009; Raven et al., 2021). The presence or absence of high oxygen concentrations in the water appears to impact the microbial diversity found in particles (Fuchsman et al., 2011; Kellogg and Deming, 2009), and particles themselves have a diffusive boundary layer, and thus can have zones of anoxia despite being in oxic waters (Bianchi et al., 2018; Ploug et al., 1997). Oxygen consumption due to the degradation of organic matter in particles leads to the development of microenvironments within the particles that can support metabolisms that would not be able to take place in the water column, such as denitrification in hypoxic waters (Bianchi et al., 2018; Fuchsman et al., 2019b) and sulfate, and metal reduction in nitrogenous anoxic zones (Raven et al., 2021; Saunders et al., 2019). The formation of anoxic zones inside particles is related to particle size (Ploug et al., 1997).
Most studies tend to assign microbial diversity to two main categories; free-living (FL) or particle-attached (PA). Groups consistently enriched on particles include the phylum Planctomycetes, the order Rhodobacterales (from the Alphaproteobacteria phylum), the phylum Cytophaga-Flavobacteria-Bacteroidetes (CFB) and some types of Gammaproteobacteria such as the order Alteromonadales (Eloe et al., 2011; Fuchsman et al., 2012, 2011; Ganesh et al., 2014; LeCleir et al., 2014; Leu et al., 2022). While some taxa can be classified as strictly FL or PA, other taxa are present in both states (Mestre et al., 2017). In the ocean, marine particles exist in a continuum of sizes through the water column (Nguyen and Harvey, 1997). A simple free-living or particle-attached characterization misses the fact that diffusion into particles is directly related to particle size and surface area to volume ratio, and therefore the interior of large particles may have different redox and nutrient concentrations than found in smaller particles (Ploug et al., 1997). The microbial diversity in individual particles increases with particle size (Stephens et al., 2024). Additionally, particle size directly affects sinking rates and by extension is thereby related to the origin and age of the particle.
A recent approach to fractionating particles, which divided them into 6 size fractions, was used both in the Mediterranean and at 8 open ocean stations during the Malaspina expedition, and showed that different sized particles host different bacterial communities (Mestre et al., 2018, 2017). The diversity was highest in the larger particles (20–200 μm), presumably because larger particles can provide more niches for microorganisms (Mestre et al., 2018, 2017). Mestre et al. (2018) found that communities on 3–20 μm and 20–200 μm size fractions at depth were distinct with the 20–200 μm size fraction more reflecting surface communities. Observations in the lab (Bižić-Ionescu et al., 2018; Datta et al., 2016) and models (Nguyen et al., 2022) suggest that microbial communities likely vary from particle to particle even within a size fraction.
The analysis of microbial amplicon sequences has been limited by the use of relative abundance as an indicator of microbial diversity; relative abundance makes it difficult to directly compare the changes in abundance of microbial groups across different samples. Indeed, relative abundance community data, when analyzed with conventional statistical approaches can produce misleading results, with changes in the actual abundances of one organism creating apparent changes in the opposite direction of other organisms (Gloor et al., 2016). Recently, several studies in microbial ecology have added quantitative elements, such as collecting microscopy data in concert with metagenomic or amplicon based measurements of relative abundance (Fu et al., 2023; Yeh and Fuhrman, 2022), combining quantitative PCR with amplicon sequencing (Shelton et al., 2023; Tettamanti Boshier et al., 2020), and using internal standard spikes (Cram et al., 2024; Jones-Kellett et al., 2024; Tourlousse et al., 2017; Wang et al., 2021). The internal spike approach produces quantitative information by incorporating an artificially designed 16S rRNA gene sequence that amplifies like other 16S rRNA genes, but that is distinct from known 16S rRNA genes of microorganisms (Tourlousse et al., 2017). The use of spike-in standards allows the quantitative estimation of numbers of bacterial 16S rRNA copies per unit of volume and mass (Cram et al., 2024). Additionally, fine scale size fractionation into seven fractions allows us to compare free-living microbes to those inhabiting on small (0.2–1.2, 1.2–5, 5–20 μm), intermediate (20–53, 53–180 μm) and large marine particles (180–500, >500 μm). Compared to the previous work by Mestre et al. (2017, 2018), we have removed the 3 μm size bin and added two larger size fractions (180–500 μm and >500 μm), increasing our chances of capturing larger sinking particles. Thirdly, we size fractionated organic matter into the same size fractions as the fractions analyzed for DNA, so the carbon content of particles in each size fraction were known. The combination of these three techniques has been used in the Chesapeake Bay Estuary (Cram et al., 2024). Here we use this approach for the first time in a depth profile in the open ocean. We specifically address how microbial composition changes among different particle sizes in a water column depth profile above the 9°50’N East Pacific Rise hydrothermal vent field.
Study site
The 9°50’N, 104°18 W East Pacific Rise (EPR) study site is of particular interest due to the presence of both an Oxygen Deficient Zone in the mesopelagic region and a hydrothermal plume in the bathypelagic region of the same water column. Oxygen Deficient Zones (ODZs) are regions of the ocean that have undetectable oxygen concentration (Revsbech et al., 2009; Tiano et al., 2014). In ODZs, bacteria and archaea use nitrate and nitrite as an electron acceptor instead of oxygen (Lam et al., 2009), leading to N2 production by heterotrophic denitrification and the autotrophic anammox process (Babbin et al., 2014; Dalsgaard et al., 2012). ODZs host 30–50% of marine N2 production (DeVries et al., 2013). The oxygen content of the Pacific Ocean has been decreasing since the 1980s (Ito et al., 2017), in consequence ODZs are expanding (Evans et al., 2023; Stramma et al., 2008). ODZs have a large anoxic core (O2 < 10 nM) surrounded by shells of suboxic (O2 < 5 μM), severely hypoxic (O2 < 20 μM) and then hypoxic waters (O2 < 60 μM). In the geographic core of the ODZ, the vertical width of anoxic waters can be >700 m. However, the vertical width of the ODZ is much smaller at the edges of the ODZ region.
At the ocean floor of 9°50’N EPR is a hydrothermal vent field that underlies the Southern edge of the Eastern Tropical North Pacific (ETNP) ODZ. It is unclear if the Oxygen Deficient Zone in this part of the EPR is truly anoxic or if oxygen is present but not measurable with Seabird sensors typically found on CTD rosettes. In oceanic regions and datasets where determination of oxygen concentration with a STOX sensor (Revsbech et al., 2009) or equivalent is not available, the presence of nitrite at concentrations of 0.5 μM or greater indicates functional anoxia (Banse et al., 2017). Thus, we measured nitrite concentrations here.
Deep sea hydrothermal vents are regions in the ocean floor where ocean water percolates to the subsurface where geothermal processes heat and enrich the water with chemicals and gas, forming a fluid that is released through orifices in the seafloor to the cold deep ocean waters (Dick, 2019). The hydrothermal vent plume at the bottom of the ocean floor at the EPR affects the chemistry and particle composition of the area surrounding the vent and the water column at depth (Dick, 2019). There are two types of plumes: buoyant plumes that rise for hundreds of meters above vents and non-buoyant plumes that spread horizontally away from the vent field. Previous analysis of buoyant and non-buoyant hydrothermal plumes using sediment traps in the EPR, revealed the presence of several microbial groups: Firmicutes, Actinobacteria, Planctomyces, and epsilon proteobacteria, which are now referred to as Campylobacterota (Sylvan et al., 2012). Non-Buoyant plumes evolve as they move away from the hydrothermal vent, with fewer vent related Campylobacterota and more SAR11, S-oxidizing Thioglobus Gammaproteobacteria and S-oxidizing Sulfurimonas Campylobacterota farther from the vent (Li et al., 2020). In our depth profile, we examined one sample designated as a non-buoyant plume sample that was 163 m from the seafloor above the vent field.
Methods
Hydrographic data
Samples were collected as part of the Hot2Cold Vents expedition aboard the R/V Atlantis in April 2019, cruise AT42-09. The expedition examined a station in the EPR located at 9°50’N, 104°18.144’W (Fig. S1) which had a bottom depth of ~2,500 m (Fig. S2). A CTD-Rosette was used to collect water samples, and attached sensors recorded oxygen (Seabird 43 Dissolved oxygen sensor), chlorophyll fluorescence (Fluorescence, WET Labs ECO-AFL/FL), beam transmission (Beam Transmission, WET Labs C-Star), temperature, depth, conductivity and salinity. Chlorophyll fluorescence was not calibrated with discrete samples during the cruise. The hydrographic data from the cruise can be found at the R2R rolling deck to repository: https://www.rvdata.us/search/fileset/144040. The EPR station was sampled with CTD casts that targeted 6 different regions: Surface at 4 m, Deep Chlorophyll Maximum (DCM) at 92 m, the mesopelagic (MESO) with two depths (275 m and 316 m), the ODZ with two depths (490 m and 527 m), bathypelagic at 1001 m and Non-buoyant Plume (NBP) at 2369 m (Table 1). The Non-buoyant Plume (NBP) was found using beam transmission data from CTD cast AT420912 (Fig. S2a) and sampled for DNA during that same cast. The temperature of this plume was not elevated compared to the surrounding seawater (Fig. S2b). Due to the volume of water and time needed to size-fractionate samples, only 2–3 depths could be sampled per cast, and samples were obtained over a 6 day period (Table 1). Samples from all the different CTD casts were combined as they were part of the same station.

Table 1. Date (local time), coordinates, CTD cast number, Sample category, Depth and Volume of water collected (L) for each water sample.
Sample preparation
Between 80 and 160 L of seawater was collected at each depth by the shipboard rosette (Table 1). Samples were obtained from Niskin bottles by opening the bottom of the Niskin bottle into an acid cleaned bucket. This ensured that particles that had sunk below the spigot of the Niskin bottle (Suter et al., 2017) were included. The water was gravity filtered through a stacked series of 142 mm nylon filters fixed between PVC connectors. The nylon filters were organized and placed one on top of each other in a gradient from the largest pore size at the top towards the smallest size pore at the bottom (500, 180, 53, 20, 5 μm). By pouring the water through the stack of filters, the particles were separated into size fractions on the nylon filters. After filtration, each nylon filter was back rinsed thoroughly using an air powered power water pistol filled with prefiltered (0.2 μm Pall AcroPak 1,500 Capsule with a Supor Polyethersulfone membrane) seawater. The fractionated marine particles were removed from the nylon mesh and collected into a clean bottle (Table 2). All the sampled fractions were divided into two parts: one part was used for DNA analysis and the second part was used to analyze particulate organic C. For DNA analysis, the 5 μm filtrate that passed through all the nylon mesh sizes was filtered onto a Polyethersulfone 1.2 μm filter and a Polyethersulfone 0.2 μm filter using a vacuum manifold. This final 0.2 μm filter collected the free-living microbial fraction. The filters for DNA analysis were placed inside 2 mL sterile screw cap cryogenic tubes, frozen in liquid nitrogen and stored at −80°C until DNA extraction in the laboratory.

Table 2. Filter sizes, particles range, filter type and treatment used to recover particles from filters.
For organic matter analysis, the fractionated resuspended particulate organic matter from each sample was filtered using vacuum via a 1.2 μm nominal pore size, 25 mm diameter, GF/C glass fiber filter (Whatman WHA1822025), that had been pre-combusted at 400°C and pre-weighed using a Sartorius microbalance (precision of 1 μg). Samples for POM analysis were stored in foil packs and also stored at −80°C until analysis in the laboratory.
Filtering may disrupt particles. Due to low organic matter concentrations in the EPR region, to carry out this study, ~120 L of water per depth (Table 1) was filtered in order to collect enough bacteria and particles to be measured. This means that filter clogging was not an issue in this study, but that the large amount of water passing through each filter could theoretically cause disaggregation. However, by passing particles through different filter sizes, and comparing their size distribution using optical scattering before and after filtration, this disaggregation was previously been shown to be minimal (Dougherty et al., 2025).
Nutrient analysis
Prefiltered (0.2 μm) dissolved nitrate, nitrite, and phosphate from the water samples were analyzed at the Nutrient Analytical Services Laboratory at the Chesapeake Biological Laboratory (CBL) where nitrite plus nitrate were measured by using the Method EPA 353.2. Phosphate by using method EPA 365.1, and Nitrate with the method ASTM D-7781-14. Nitrate concentrations were calculated by subtracting nitrate concentrations from the nitrite plus nitrate numbers. The combined nitrate + nitrite and phosphate concentrations were used to calculate N deficit as described in Chang et al. (2012), a variation on N* that uses N:P ratios from nearby oxic waters, which are ~14, rather than assuming a ratio of 16.
Particulate organic matter analysis
POM filters were dried at 40°C for >24 h and then were fumed with concentrated HCl for 24 h to remove carbonates. After this treatment the POM filters were packed into silver capsules and sent to the UC Davis Stable Isotope Facility (Davis, CA) to obtain organic C and N concentrations and the natural stable isotopes composition of δ15N and δ13C. Both carbon and nitrogen were quantified in an Elementar Vario EL Cube coupled with an Isoprim VisION isotope ratio mass spectrometer. A more detailed examination of organic C and N concentrations, ratios and stable isotope data are shown elsewhere (Fuchsman and Cram, 2024).
DNA extraction and spike-in preparation
DNA was extracted following the Cram et al. (2024) protocol. Frozen samples were extracted from DNA filters using hot SDS lysis followed by phenol-chloroform purification. The DNA concentration of the samples was measured using a Qubit 4.0 (ThermoFisher). Samples were then preserved at −80°C until further analysis.
A synthetic 16S rRNA gene sequence LC140931.1 Synthetic 16S Spike-In (Tourlousse et al., 2017) was synthesized as gBlocks Gene Fragments by Integrated DNA Technologies, Inc. The gBlock was resuspended in a solution of 0.1 mg/mL of PolyA (Roche Diagnostics Deutschland GmbH) in PCR grade water (DEPC-Treated Water Ambion™) and then serially diluted to reach a final concentration of 106 copies/μl following the manufacturer’s instruction.
DNA concentration per sample was diluted using PCR-grade water to reach a concentration of 2 ng/μl DNA (in the case of samples with higher DNA concentrations); then the spike resuspended in PolyA was added to each sample to reach a final concentration of 104 copies/μl or for samples with low DNA concentrations, 5*103 copies/μl.
DNA sequencing and amplicon bioinformatics
The extracted DNA was sequenced with an Illumina MiSeq by a third-party sequencing company, Molecular Research (MRDNA) Shallowater, Texas, USA. In brief, at MRDNA, 1 ng of sample DNA plus the 104 or 5 * 103 copies of spike were amplified with the HotStarTaq Plus Master Mix Kit (Quiagen, USA) using universal 515FY-926R primers (Parada et al., 2016) with a barcode added to the forward primer. Amplification conditions were 94°C for 3 min followed by 30 cycles of 94°C for 30 s, 53°C for 40 s and 72°C for 60 s. There was a final elongation step of 72°C for 5 min. Following amplification, PCR products were checked for length and intensity on a 2% agarose gel. Samples were then pooled together in equal proportions and purified using Ampure XP beads. These pooled cleaned products were sequenced on a MiSeq instrument (2×250 base pair reads) following manufacturers guidelines. Samples were returned to our group as FASTQ files.
The EPR sequenced amplicons were processed following the same pipeline described in Cram et al. (2024), which will be briefly summarized here. The sequenced amplicons were processed using the DADA2 pipeline (Callahan et al., 2016). The primers were removed from the sequences using Cutadapt (Martin, 2011), with the following modification, the filterandtrim() algorithm was used to retain sequences with < three errors in the forward read and fewer than five errors on the reverse read. Then DADA2 was used to dereplicate, obtain error rates, and merge forward and reverse reads. The taxonomy of each ASV was assigned using the Silva database version 132 released Dec 2017 (Quast et al., 2013). This database was modified by adding the spike-in sequences (Cram et al., 2024). Data were analyzed using the R software for statistical analysis (V4.0.2) using the packages Vegan, Phyloseq for data analysis, Dplyr for data curation and ggplot2 to visualize data.
Calculating taxon abundance
The taxon abundance (16S + 18S rRNA gene copies/L) of each microbial group in each sample was estimated by normalizing with the spike in read numbers, DNA extracted and volumes filtered, following Equation 1. In two cases, no spikes were observed in our sample. Our lowest sequence depth was 7,052 sequences in the 0.2 μm 1,000 m sample, which corresponded to no observed spikes, presumably because 104 copies was not enough spike to detect in every sample. In the other sample without spiked standards at 5 μm 490 m, there were 25,162 read copies, which is large enough that we would have expected to see spiked standards in this sample. We fear that spiked standards were not added to this sample. For these two cases, we used a spike count of one, rather than zero, in order to show the taxonomic information for these samples.
The quantitative examination of microbial amplicons in particles using spike-in standards allowed for calculating the number of bacterial and archaeal gene copy numbers. Using these same primers, picocyanobacteria, which are known to have only one copy of the 16S rRNA gene, closely corresponded to cell counts from flow cytometry data (Jones-Kellett et al., 2024). When examining whole microbial communities, microbial abundance is directly related to 16S rRNA copy numbers, but the copy numbers are often somewhat higher than microscopic count data (Cram et al., 2024). Therefore, we have been careful herein to report 16S rRNA gene copy abundances, rather than cellular abundances.
Calculating particle carbon content
To estimate the amount of organic carbon present in the original volume of water sampled, total observed μg C was normalized to the volume of the water that went through the stacked set of nylon filters and the fraction of water that was rinsed through the GF/F filter, using Equation 2.
Microbial abundance was normalized to carbon weight by dividing the taxon abundance by the carbon composition of marine particles following Equation 3.
We conducted a redundancy analysis to see if there was a statistical effect of particle size (μm) or depth (m) on community structure. In order to analyze the relative similarities in community structure between all samples, microbial abundance was log-transformed and depth was treated as continuous data. Samples were color and shape coded based on size class and the depth from which they were collected. In particular, we used an RDA of form ‘log(OTU_Abundances) ~ log(Depth) + log(Depth)^2 + log(Size) + log(Size)^2 + log(Depth) * log(Size).’ This model looks for non linear effects of depth and size on community structure as well as interactions between depth and size effects. We tested for statistical significance using an ANOVA-like Permutation Test for Redundancy Analysis.
Results
Hydrographic and biochemical context for the EPR
First we examine the hydrographic and biogeochemical context for this station. Salinity (33.5), temperature (29.5°C) and oxygen concentrations (198 μM) were constant for the top 75 m of the water column. Then values changed sharply in the 75–140 m region with salinity increasing to 34.8 and temperature and oxygen concentrations decreasing to 14.5°C and 6 μM, respectively, (Figure 1). The concentration of chlorophyll a was low in the top 10 m (~0.1 mg/m3), and there was a deep primary chlorophyll maximum at 60–80 m of 1.3 mg/m3 (Figure 2A). Between 140 and 400 m depth, oxygen concentrations fluctuated between 5 and 10 μM and temperature slowly decreased to 10°C. In the area between 400 and 700 m, the oxygen concentration was constant at 2.6 μM (Figure 1). This value may have been the sensor detection limit. While other Seabird sensors can reach 1 μM, detection limits can depend on the individual instrument (Kwiecinski and Babbin, 2021). Thus, it is not possible to tell if the oxygen depleted region was anoxic. Below this region, oxygen concentrations increased again, reaching 17 μM at 1000 m while temperatures continued to decrease, reaching 4.7°C. In the deepest waters, oxygen concentrations reached 100 μM and temperatures decreased to 1.9°C (Figure 1).

Figure 1. Hydrographic parameters for the EPR station. (A) Oxygen concentrations for the full water column with an inset for the top 300 m. (B) Salinity and temperature for the full water column with an inset for the top 300 m. (C) Nitrate and phosphate concentrations. Dashed lines indicate the boundaries of the ODZ.

Figure 2. Biochemical profiles at the EPR station. (A) Oxygen (μM) and Chlorophyll concentrations (mg/m3). (B) Nitrite concentration and nitrogen deficit, through the top 1,000 m of the water column. Dashed lines indicate the boundaries of the ODZ.
The presence of nitrite was observed in the water column in the EPR (Figure 2B). Nitrite accumulation was observed at two peaks: a primary nitrite maximum at 103 m, where nitrite reached 0.71 μM, and a secondary nitrite maximum at 0.32 μM at 530 m (Figure 2B). Banse et al. (2017) used 0.5 μM nitrite as a cut off for using nitrite to indicate true anoxia. The secondary nitrite maximum found here is close but less than this cut off, but higher values could have missed due to the sampling scheme. Nitrate and phosphate were also measured in the EPR water column. Nitrate was ~1.5 μM in the top 60 m and then increased to 5 μM at 100 m and to 40 μM at 1000 m depth (Figure 1C). Phosphate generally had a depth profile pattern similar to nitrate with 0.3 μM phosphate in the surface, 0.5 μM at 100 m and a maximum of 3.6 μM at 1000 m (Figure 1C). In addition, a nitrogen deficit of 10.9 μM at 530 m depth in the ODZ region was calculated from this nutrient data (Figure 2B). Outside the ODZ, the nitrogen deficit was <2 μM. These N deficit values are typical for those measured in the ETNP ODZ (Fuchsman et al., 2018), and could reflect local denitrification under anoxic conditions, however the nitrogen could also have been removed elsewhere and this low N water advected into the local system (Peters et al., 2018). Thus, the suboxic versus anoxic condition of the water in the EPR ODZ cannot be determined by this chemical data alone.
Concentrations of carbon and microbial 16S rRNA in particles
In the surface sample (4 m), the total combined POC concentration was 16 μg/L. Two samples were obtained in the DCM (92 m) three days apart (Table 1), and their total combined POC concentrations varied between 25 and 35 μg/L. The highest total combined POC concentration (54 μg/L org C) was at 275 m (Figure 3A), and it was still elevated at 316 m (30 μg/L). In the ODZ and at 1000 m, the total combined POC concentrations were 14 μg/L or less (Figure 3A). The non-buoyant plume (2,369 m) contained the smallest amount of total organic C (8 μg/L). At depths above the ODZ, particulate organic C (POC) concentrations were highest in the 5–20 μm fractions but there were still significant concentrations in the 20–53 μm and the 53–180 μm fractions (Figure 3A). However, in the ODZ (490 and 527 m) and below, the POC concentrations in the 5–20 μm fraction were particularly reduced (Figure 3A). In fact, at 1000 m, POC concentrations were highest in the 20–53 μm fraction (Figure 3A). Thus in the top 350 m of the water column, the POC concentrations in small particles (<20 μm) were larger than the concentrations in intermediate particles (20–180 μm), which in turn were larger than the concentration in large particles (>180 μm). However, at depths >350 m, intermediate sized particles often had the largest concentrations of organic C (Figure 3A).

Figure 3. Organic C concentrations and bacterial 16S rRNA copies across size fractions in the EPR. (A) Particulate organic carbon concentrations in marine particles. The >500 μm organic matter fraction for 275 m is missing. (B) Bacterial 16S rRNA copies/μg organic C across seven size fractions. Inset: Bacterial 16S rRNA copies/μg organic C zoomed in on 4, 92 and 316 m. (C) Bacterial 16S rRNA copies/L across seven size fractions. Inset: Bacterial 16S rRNA copies/L zoomed in on 92 and 316 m. (D) Archaeal 16S rRNA copies/L across seven size fractions. For 16S rRNA copies, the 275 m and one 92 m sample were not included because too many size fractions were missing. Additionally, 1.2–5 μm size fraction is missing from the 527 m sample, 53–180 μm fraction is missing from the 4 m sample, and the >500 μm fraction is missing from the 316 m sample. The 5–20 μm size fraction for 490 m was not included. DCM, deep chlorophyll maxima, ODZ, oxygen deficient zone, MESO, mesopelagic, BATHY, bathypelagic; NBP, non-buoyant plume.
We examined the abundance of archaea and bacteria on particles in two ways: in 16S rRNA copies per L, and copies/μg org C. When the size fractions were examined together, one can see different patterns in the upper and lower water column. In the upper water column, the numbers of bacterial 16S rRNA copies/μg org C on >180–500 μm and >500 μm particles was much smaller (108 copies/μg org C) than found in the deeper water column (109 copies/μg org C) (Figure 3B). However, for the smaller 5–20 μm particle size, the abundance of 16S rRNA gene copies/μg org C decreased as depth increased; bacterial copies/μg org C abundance was higher in surface (108 copies/μg org C), decreasing to 107 copies/μg org C in the mesopelagic and ODZ and decreasing again to 105 copies/μg org C at 1000 m (Figure 4). However, values then increased to ~109 copies/μg org C in the Non-buoyant plume at 2,369 m depth (Figure 4). The Non-buoyant plume had a fundamentally different structure of bacteria on particles, with extremely high numbers of 16S rRNA copies/μg org C in the 5–20, 20–53, and >500 μm fractions (>109 copies/μg org C), but extremely low numbers in the 53–180 and 180–500 μm fractions (105 copies/μg org C) (Figures 3, 4). For these copies/μg org C numbers, it is difficult to separate changes in 16S rRNA copies from changes in organic C concentrations. Therefore, we also examine 16S rRNA copies/L for each size fraction, which indicates whether the numbers of microbes are increasing or decreasing irrespective of changes in organic C.

Figure 4. Depth profile of 16S rRNA gene copy numbers per μg organic C, aggregated to domain level. Each panel represents particle size (in μm). The Y axis corresponds to depth (m), and the X axis corresponds to 16S rRNA gene copies/μg organic C. Dashed lines indicate the boundaries of the ODZ.
In the upper water column, bacterial 16S rRNA copies per L were higher in the free-living fraction and smaller particles and lower in large particles. However, in the ODZ and below, bacterial 16S rRNA copies per L were higher in intermediate and large particles (Figure 3C). The mesopelagic samples (316 m) had particularly low numbers of bacterial 16S rRNA copies overall, which was surprising since the organic C concentrations were particularly high at this depth (Figure 3A). Archaeal 16S rRNA copies per L, on the other hand, were primarily found in the free-living and 1.2–5 μm fractions (Figure 3D). Archaea abundance was always lower than bacteria abundance usually by one or two orders of magnitude, but the depth profiles for bacteria and archaea had similar shapes (Figure 5).

Figure 5. Depth profile of 16S rRNA gene copy numbers per L, aggregated to domain level. The Y axis corresponds to depth (m), and the X axis corresponds to 16S rRNA gene copies/L. Dashed lines indicate the boundaries of the ODZ.
The number of bacterial 16S rRNA copies per L in the different size fractions varied with depth. In the free-living size fraction, there were high 16S rRNA copies per L in the surface (~109 bacterial 16S rRNA copies/L), which decreased in the mesopelagic oxycline (107 copies/L), but increased again in the ODZ (~109 bacterial 16S rRNA copies/L) (Figure 5). Similar to the copies/μg org C data, the numbers of bacterial 16S rRNA copies/L in the 5–20 μm fraction decreased from the surface to the oxycline and decreased further in the bathypelagic. Also similar to the copies/μg org C data, in the non-buoyant plume, the highest copy numbers of 16S rRNA were found in the 5–20 and 20–53 μm fractions and the >500 μm fractions (~109 copies/L); however, 16S rRNA copy numbers were quite low (~106 copies/L) in the 53–180 and 180–500 μm fractions. Thus, the general trends were similar between copies per L and copies/μg org C data.
Abundance of bacterial taxa on particles
The most abundant phyla on particles were Bacteroidetes, Cyanobacteria, Planctomycetes and Proteobacteria, but Actinobacteria, Marinimicrobia, Chloroflexi, and Thaumarcheota were also key players in the system that were primarily free-living (Figure 6). We calculated relative microbial abundance (%) and quantitative microbial abundance (in 16S Copies/L). Proteobacteria dominate in the quantitative data at each particle size and across all depths (Figure 6B). The abundance of 16S rRNA gene copy numbers/L in Proteobacteria was up to 109 16S rRNA gene copies in the larger (>500 μm) particle size, which was similar to values for the free-living fractions (Figure 6B). In relative abundance units, Proteobacteria were still the dominant group (reaching almost 100% in the particles ≥500 μm) and we can see similar changes in Proteobacteria relative abundance as in the quantitative data (Supplementary Figures S3, S4). However, as seen previously (Gloor et al., 2016), the variation of smaller groups appear to be determined by variations in the determined by variations in the dominant group. For example, when examined quantitatively, we can see the changes in the less abundant Planctomycetes phylum across depth in the free-living fraction, including a maximum in the ODZ (Figure 6B). With relative abundance, we lost the information about the presence and variation of Planctomycetes, compared to the most abundant group (Proteobacteria) (Supplementary Figures S3, S4).

Figure 6. Depth profile of volume normalized 16S rRNA gene copy abundance of top 10 most abundant Phyla in the EPR. The ten most abundant Phyla in the entire dataset are shown. Divided into parts (A,B) for easy visualization. Each panel represents different particle sizes (in μm). The Y axis represents depth in meters, and the x axis represents 16S rRNA gene copies per L. The area between the dashed lines corresponds to the Oxygen Deficient Zone.
The majority of cyanobacteria were from the order Synechococcales, otherwise known as picocyanobacteria (Synechococcus and Prochlorococcus) (Supplementary Figure S5). The picocyanobacteria are generally free-living, but known to attach to particles under specific conditions (Capovilla et al., 2023). Abundances of Cyanobacteria, which use light for photosynthesis, changed sharply with depth. In the euphotic zone, Cyanobacteria were the second most abundant phyla on particles for almost all size fractions (Figure 6A). For the free-living fractions 0.2 and 1.2–5 μm particles, Cyanobacteria decreased from 109 16S rRNA gene copies/L in the surface to 106 16S rRNA gene copies/L in the mesopelagic region. For particles of size 5–20 μm, the decrease was almost linear, changing from 108 16S rRNA gene copies/L in the surface to 103 in the mesopelagic. In the size fractions 20–53 μm and 53–180 μm, Cyanobacteria were 108 16S rRNA gene copies/L in the euphotic zone while for particles >500 μm Cyanobacteria were 107 16S rRNA gene copies/L. Thus, though cyanobacteria were abundant on particles, the cyanobacteria were one order of magnitude less abundant on small particles than in the free-living fraction and two orders of magnitude less abundant on large particles.
Among Archaea, the phylum with the highest 16S rRNA gene copies/L corresponded to Thaumarchaeota (Figure 6B). Euryarcheota were more abundant in the <5 μm fractions in the euphotic zone, but otherwise Thaumarchaeota had the highest abundances (Figure 6A). Thaumarchaeota abundance changed with depth and across particle sizes. For example, in the free-living size fraction (0.2–1.2 μm) Thaumarchaeota increased from <106 16S rRNA gene copies/L in the surface to almost 109 16S rRNA gene copies/L in the ODZ region and remained high at depth, and a similar pattern was observed in particle size 1.2–5 μm. However, in particle 5–20 μm, the abundance of Thaumarchaeota decreased with depth from ~106 16S rRNA gene copies/L close to the surface to ~105 16S copies/L in the ODZ region to 102 16S rRNA gene copies/L in the bathypelagic region (1,001 m depth). In particle sizes ≥20 μm the Thaumarchaeota abundance mostly remained between 106 and 107 16S rRNA gene copies/L, including in the bathypelagic region (1,001 m). The abundance of Thaumarchaeota in the non-buoyant plume depth was different at each particle size; it was high 107–108 copies/L in the small size fractions (from 0.2 to 20 μm), but in the 53–180 μm and 180–500 μm size fractions, Thaumarchaeota was <102 copies/L. Finally, on the ≥500 μm size fraction in the non-buoyant plume, there were ~107 16S rRNA gene copies/L. Thus, while Thaumarchaeota were most abundant in the free-living fraction, they were also present in the majority of particles, but at several orders of magnitude lower abundance.
Among the Proteobacteria, the most abundant class in the particle fractions was consistently the Gammaproteobacteria, reaching >109 copies/L in each of the size fractions >20 μm followed by Alphaproteobacteria, reaching 109 copies/L (Figure 7). The most abundant Gammaproteobacterial orders on large particles were Alteromonas (~108 copies/L) and Oceanospiralles (Figure 8) (~108 copies/L). In the free-living fraction Alphaproteobacteria (>109 copies/L) were more abundant than Gammaproteobacteria (108 copies/L) (Figure 7). In the free-living and small particles 1.2–5 μm, the abundance of Alphaproteobacteria varied from 107 to 109.5 copies/L in the water column. In the 5–20 μm size class Alphaproteobacterial abundance decreased almost linearly from the surface ~108 to ~105 copies/L in the Bathypelagic region. In particles in sizes >53 μm, Alphaproteobacteria abundance varied between ~107 to ~109 copies/L, but Alphaproteobacteria abundance increased from the surface to the ODZ. Additionally, Deltaproteobacteria were present throughout the water column, and in all size fractions though at lower abundances than Gammaproteobacteria and Alphaproteobacteria (Figure 7). The abundance of Deltaproteobacteria decreased as particles increased in size; in the free-living and 1.2–5 μm fractions abundance was 106–108 copies/L, but in each of the particle size fractions, abundance was mostly 106 copies/L or less. Members of Zetaproteobacteria were only found in the ODZ and had low numbers (104 copies/L) as free living (0.2–1.2 μm) fraction but were associated with particles 105–106 copies/L on 20–53 and 53–180 μm in the ODZ region (Figure 7).

Figure 7. Depth profile of volume normalized 16S rRNA gene copy abundance of Proteobacteria, aggregated to class level in the EPR. The five most abundant classes in the entire dataset are shown. Each panel represents different particle sizes (in μm). The Y axis represents depth in meters, and the x axis represents 16S rRNA gene copies per L. The area between the dashed lines corresponds to the Oxygen Deficient Zone.

Figure 8. Depth profile of volume and bin size normalized 16S rRNA gene copy abundance of Gammaproteobacteria, aggregated to order level. The five most abundant orders in the entire dataset are shown. Each panel represents different particle sizes (in μm). The Y axis represents depth in meters, and the x axis represents 16S rRNA gene copies per L. Dashed lines indicate the boundaries of the ODZ.
SAR11 and Rhodobacterales were the most abundant Alphaproteobacteria orders present on all sizes of particles (Figure 9). Rhodobacterales clearly increased in abundance in the free-living fraction in the ODZ as did Rhodospirillales and SAR11. Rhodobacterales were the most abundant Alphaproteobacteria on >5 μm particles in the ODZ and below while SAR11 was the most abundant Alphaproteobacteria on particles in the top 300 m. SAR11 bacteria are an order of Alphaproteobacteria that are thought to be free-living (Giovannoni, 2017). In general, the total SAR11 abundance was higher in the free-living size class (0.2–1.2 μm) where the abundances of SAR11 decreased from ~109 16S rRNA gene copies/L in the surface to 106 16S rRNA gene copies in the mesopelagic oxycline above the ODZ, but in the ODZ SAR11 increased to ~108 16S rRNA gene copies/L, and remained relatively constant at deeper depths. On larger particles, the abundance of SAR11 varied between ~107 to 108 16S rRNA gene copies/L (Figure 9). Within SAR11, we detected four clades (Clades I, II, III, and IV) as well as Unclassified SAR11 that did not fall into any clade (Supplementary Figure S6). The clade that had the highest abundances both in the free-living and particle fractions was Clade I, closely followed by Clade II. Both clades were found at almost every analyzed depth. The abundance of all SAR11 clades decreased as particle size increased (Supplementary Figure S6). Thus, there was no significant shift in clades between the water and particles. All these clades break down into smaller ecotypes, which is the more appropriate level to examine for environmental distributions (Haro-Moreno et al., 2020), but which we cannot resolve with these methods.

Figure 9. Depth profile of volume normalized 16S rRNA gene copy abundance of Alphaproteobacteria, aggregated to order level in the EPR. The five most abundant orders in the entire dataset are shown. Each panel represents different particle sizes (in μm) The Y axis represents depth in meters, and the x axis represents 16S rRNA gene copies per L. Dashed lines indicate the boundaries of the ODZ.
The majority of taxa from the Planctomycetes phylum are known particle attached bacteria except for anammox bacteria which are free-living in marine ODZs (Fuchsman et al., 2012). Examination of the Planctomycetes phylum as a whole, indicates that the samples with the highest Planctomycetes were the free-living fractions in the ODZ region ~108 copies/L (Figure 6A). In the surface waters (0.2–1.2 μm), when examined in comparison to Proteobacteria >109 copies/L, there were few Planctomycetes 106–107 copies/L. At non-ODZ depths, Planctomycetes remained at 106–107 copies/L in each of the particle size fractions 20 μm and above while Proteobacteria decreased to ~108 copies/L (Figure 6A). Thus, in terms of relative abundance, Planctomycetes were enriched in particles (Fig. S3 and S4), but in actual abundances, they were similar between free-living and suspended particles (Figure 6). We looked into the abundance of classes inside Planctomycetes. Planctomycetes of the group Brocadia, which mediate the anammox reaction, were highly abundant (~108 copies/L) in the free-living fraction in the ODZ (Figure 10), which is consistent with previous data (Fuchsman et al., 2017, 2012; Ganesh et al., 2014). The class Planctomycetacia, which includes Pirellulaceae and Planctomycetales orders among others, had the highest abundance across particle sizes among the Planctomycetes, with 106–107 16S rRNA gene copies/L in each fraction (Figure 10). However, in the ≥500 μm fraction the abundance of Planctomycetacia decreased to ~106 16S rRNA gene copies/L (Figure 10).

Figure 10. Depth profile of volume normalized 16S rRNA gene copy abundance of Planctomycetes, aggregated to Class level. The five most abundant classes are shown for the entire dataset. Each panel represents different particle sizes (in μm). The Y axis represents depth in meters, and the x axis represents 16S rRNA gene copies per L. The area between the dashed lines corresponds to the Oxygen Deficient Zone.
Bacteroidetes, particularly Flavobacteria, is a second group that is generally considered particle attached (Delong et al., 1993; Eloe et al., 2011; Fuchsman et al., 2011; LeCleir et al., 2014; Leu et al., 2022). In the EPR, Bacteroidetes was the second most abundant phyla on large particles below the euphotic zone, after Proteobacteria. Our data indicate that in the free-living fraction the Bacteroidetes group was more prevalent in the euphotic zone than at depth. However, in the large particles, Bacteroidetes copies were similar in the euphotic zone and at depth. We observed that Bacteriodetes had a higher number in small particle fractions with 108 copies/L in both free-living and < 20 μm particles but 107–108 copies/L in particle fractions >20 μm (Figure 6A). Thus, the pattern was fairly similar between Planctomycetes and Bacteroidetes, both of which are generally considered particle attached.
We also examined the abundance of nitrite-oxidizing Nitrospina in our samples due to its ecological relevance in the nitrogen cycle. Nitrospina abundance was higher in the free-living fraction (0.2–1.2 μm) specially inside the ODZ region, where Nitrospina abundance increased to ~108 copies/L (Supplementary Figure S7). In the surface and in the non-buoyant plume the abundance was one order of magnitude lower ~106 copies/L. The abundance of Nitrospina decreased one order of magnitude in the larger particle sizes, at 1.2–5 μm the highest abundance was in the ODZ with ~106 copies/L, and in the larger particle size (≥20 μm) the Nitrospina abundance decreased to ~105–106 16S rRNA gene copies/L. Thus, Nitrospina was primarily free-living.
Redundancy Analysis (RDA) comparing taxonomic diversity against depth and size fraction generally suggested that large particles (≥180 μm) appear to be more similar to each other, regardless of location, and that smaller particles appeared to be more dissimilar between depth regions (Figure 11). We observed significant nonlinear (quadratic) relationships of both depth (p = 0.005) and size fraction (p = 0.005), as well as a significant interaction between these variables (p = 0.001), in shaping community structure.

Figure 11. Redundancy analysis plot (RDA) of the relationship between microbial community structure, particle size “fraction” and depth. Points that are closer together in space have more similar community structure. The RDA model was fitted using log-transformed values of depth and fraction, including quadratic and interaction terms, to capture nonlinearities and synergistic effects. For clarity, however, the biplot displays arrows labeled simply as depth and fraction, which represent the main-effect coefficients derived from the full model. Arrow lengths indicate the strength of the relationships, while their directions show the gradients along which community composition shifts. Although the plotted coefficients are not log-transformed, all analyses were conducted on log-transformed data, and the higher-order terms were statistically significant (p < 0.01). The first two axes explain 18 and 11% of the total sample variance.
Discussion
The analysis of microbial amplicon sequences has been limited by the use of relative abundance as an indicator of microbial diversity. Here we use internal spike-in standards to quantitatively address how microbial composition and organic C change among seven particle size fractions in a depth profile of the open ocean water column located above the 9°50’N East Pacific Rise hydrothermal vent field. The 9°50’N EPR region examined here had a defined low oxygen zone where oxygen was below detection (400–700 m depth), and in the deepest region sampled reached a non-buoyant plume (2,369 m depth) sourced from a hydrothermal vent.
Ecology of sinking particles in the EPR
Our data showed a statistically significant, though non-linear relationship between depth and size fraction and the microbial community structure (Figure 11), reflecting previous findings using related methods in the global ocean (Mestre et al., 2018) and in the Chesapeake Bay (Cram et al., 2024). Together these data suggest that microbial communities have biogeographic ranges that are defined not only at the macro-scale across depth and space, but also at the microscale between different sizes of particles. There was a large increase in the bacterial 16S rRNA copies/ug C and copies/L on particles 180–500 μm and >500 μm between the upper water column and the deep water column (Figure 3). This implies that, even though the total concentration of organic C in particles decreased in the deep water column, the density of bacteria on large particles increased at depth.
In our dataset, the most abundant bacterial phyla on particles were Bacteroidetes, Cyanobacteria, Planctomycetes, Gammaproteobacteria and Alphaproteobacteria (Figures 6, 7). These results are generally comparable to those found previously, with some key exceptions. Campylobacterota have mismatches with the primers used (McNichol et al., 2021), and were barely present in this dataset. This was unfortunate since Campylobacterota (Arcobacter), and Alteromonadales Gammaproteobacteria dominated metagenomes from sinking particles at 4000 m depth at HOT (Boeuf et al., 2019). Additionally, Campylobacterota (Sulfurovumaceae and Sulfuricurvaceae) have been found to be particularly enriched in sinking particles from hydrothermal vent plumes sampled 30 m and 115 m from the closest vent chimneys (Sylvan et al., 2012). In our 2,369 m depth particles, Campylobacterota, were <104 copies/L, but Alteromonadales were present in high numbers at >109 copies/L as expected (Figure 8). Thus we illustrate that, even when using well established primer sets regularly used in the marine environment, (e.g., Reintjes et al., 2023; Yeh and Fuhrman, 2022), PCR bias causes environmentally important groups to be missed.
Considering traditional free-living and particle attached bacteria quantitatively
Marine particles were inhabited by both traditionally particle attached bacteria like Planctomycetes and Bacteroidetes, as well as by free-living SAR11, Thaumarchaeota and picocyanobacteria. In surface waters, where Thaumarchaeota are typically not abundant (Santoro et al., 2017), numbers were similar between Thaumarchaeota on particles and in the water column. At depth, Thaumarchaeota were present in >5 μm particles, but they were an order of magnitude more abundant in the free-living and 1.2–5 μm fractions. 16S rRNA copies of SAR11 bacteria were abundant in the free-living fraction (0.2–1.2 μm), and were also present in larger particle sizes, but in lower abundances (Figure 9). SAR11 has been previously described as a free-living organism (Giovannoni, 2017), and it dominates free-living size fractions, (e.g., Reintjes et al., 2023). This is in concordance with the high abundances in the free-living fraction we observed in our samples (Figure 9). However, in the top 300 m, SAR11 was also often the most abundant Alphaproteobacteria on particles (Figure 9). Despite SAR11’s presence on particles, it was still two orders of magnitude lower in abundance than in the free-living fraction. Several other papers with size fractionated amplicon data have found SAR11 on particles (Cram et al., 2024; Mestre et al., 2020; Yeh and Fuhrman, 2022). However, it is possible that abundant free-living microbes are entrained in particles as a result of viral infection (Fuchsman et al., 2019a), grazing and defecation, or physical aggregation processes. Previous work indicated that the surface Thaumarcheota ecotypes were the main ecotype found in deep particles (Ijichi et al., 2019), which is consistent with entrainment of free-living microbes into particles in the euphotic zone. The equal numbers of Thaumarcheota on particles and free-living in surface waters is harder to explain. It is also possible that free-living microbes were entrained into the particle fractions while sampling, but as seen in the methods, these samples were gravity filtered and had no clogging. SAR11 bacteria have evolved a particularly non-sticky less hydrophobic cell surface to avoid predation (Dadon-Pilosof et al., 2017). Thus, SAR11 bacteria should be less likely to stick to filters than other bacteria, so we suspect that their presence on particles is not an artifact.
Planctomycetes are considered a particle attached group (Delong et al., 1993; Eloe et al., 2011; Fuchsman et al., 2012; Ganesh et al., 2014; LeCleir et al., 2014; Leu et al., 2022). The abundance of Planctomycetes in particles, could be explained by their ecological role as they respire complex organic matter such as polysaccharides (Sun et al., 2021). The highest abundance of Planctomycetes 16S rRNA was actually found in the free-living fractions in the ODZ (~108 copies/L) due to autotrophic anammox bacteria (Brocadia; Figure 10). At non-ODZ depths, there were similar numbers of Planctomycetes 16S rRNA copies in the free-living fractions and in each particle fraction (Figure 6A). However, if the particle fractions are added together, more Planctomycetes 16S rRNA was present on particles than in the water column. So, in our analysis, Planctomycetes are still a particle attached group, but are also found in the water column in higher abundances than previously thought.
We can also examine Bacteroidetes, a second group that is generally considered particle attached, (Delong et al., 1993; Eloe et al., 2011; Fuchsman et al., 2011; LeCleir et al., 2014; Leu et al., 2022). We observed that Bacteroidetes 16S rRNA had a higher abundance in small particle fractions (< 20 μm particles) and in the free-living fraction, with ~108 copies/L in both, but <107 copies/L in particle fractions >20 μm. Thus, the pattern was fairly similar between Planctomycetes and Bacteroidetes, both of which are generally considered particle attached. These data, in which particle associated organisms are actually similarly abundant in the free-living phase, cause us to rethink the conventional understanding of the habitats of “particle attached” organisms.
We can see that, integrated across all particles, Planctomycetes and Bacteroidetes are more abundant in particles than in the water column. However, their abundance in the free-living phase has previously been overlooked. The significant presence of particle attached bacteria in the free living phase is consistent with the published models that assume that particle attached bacteria attach to a particle and then detach again as the particle sinks (Nguyen et al., 2022). In this model, microbes are constantly colonizing particles as they sink rather than having colonization only at the surface, allowing the microbes to operate in the pressure zone where they are adapted (Nguyen et al., 2022). Microbial lifestyles complicate, but are generally consistent with, this picture. For example, many Planctomycetes have a sessile attached life phase that produces swimming swarmer daughter cells (Fuerst and Sagulenko, 2011). In this case, the Planctomycetes in the free-living and particle attached fractions may be in different life stages. We do know that some particle attached microbes, such as Bacteroidetes and Planctomycetes, are transported to the bathypelagic ocean during high flux pulses (Poff et al., 2021), but these microbes may not thrive there. About 10% of the bathypelagic microbes (1000–4,000 m) have reduced activity under high pressure, including Bacteroidetes and Alteromonas (Amano et al., 2022). It is the free-living reservoir of particle attached microbes that allows attaching to particles to be a legitimate evolutionary strategy.
Microbial communities in the ODZ
In the EPR, the Seabird oxygen sensor on the CTD reported oxygen concentrations of 2.6 μM from 450 to 850 m. Standard CTD oxygen sensors are unable to differentiate very low oxygen concentrations and oxygen concentration that are at a minimum value but unchanging over depth have previously been used to identify anoxic water (Kwiecinski and Babbin, 2021), suggesting that oxygen concentrations were likely less than the detection threshold of the CTD sensor. Additionally, this same region harbored a N deficit and a secondary nitrite maximum, suggesting the presence of denitrification (Figure 2). The secondary nitrite maximum was only 0.3 μM, which is close to but below the threshold of 0.5 μM suggested to confirm anoxia (Banse et al., 2017). However, our nutrient sampling was sparse with 1 point in the ODZ, so depths with higher concentrations could easily have been missed. The N deficit was in the same concentration range as previously seen in the core of the ETNP ODZ (Fuchsman et al., 2018), and we cannot prove that the this signal was produced in situ in the EPR region rather than advected from the ETNP core. With this chemical dataset, we cannot determine the difference between suboxia and anoxia. However, the microbial community can give us insight in this manner. The high abundance of Thaumarchaeota (ammonia oxidizing archaea) and the presence of nitrite oxidizing bacteria Nitrospina in the ODZ (Figure 6 and Supplementary Figure S7), implies that small amounts of oxygen were likely still present. While Nitrospina can use oxygen at <4 nM, Thaumarcheota archaea have a higher oxygen threshold (Km = 333 nM O2) (Bristow et al., 2016). In the core of the truly anoxic ETNP ODZ, Thaumarcheota and Nitrospina were not present (Fuchsman et al., 2017). On the other hand, anammox bacteria, which are 50% inhibited at 880 nm O2 (Dalsgaard et al., 2014), were also present at this EPR station. The ecosystem in the EPR ODZ may be similar to that of the extremely hypoxic Bay of Bengal, with oxidative N processes dominating in the water column but anammox still present (Bristow et al., 2017).
Bacterial and archaeal 16S rRNA abundance increased in the ODZ region in the free-living fraction, and in fractions >20 μm (Figure 5). The 16S rRNA abundance in the 5–20 μm fraction did not increase in the ODZ. Rather, abundances had a linear decrease from ~108 copies/μg org C at the surface towards 106 copies/μg org C for bacteria at 1000 m depth (Figure 4). However, on larger particle sizes (≥20 μm) we observed a decrease in abundance from the surface to the low oxygen mesopelagic water above the ODZ, but in the ODZ, (400–700 m depth), the abundance increased in all fractions ≥20 μm, reaching ~109 copies/μg org C. Assuming the EPR is representative of ODZs, reduced flux attenuation in ODZs (Keil et al., 2016; Van Mooy et al., 2002) does not appear to be due to the presence of fewer microbes on particles.
We can examine which microbes are causing this increase in 16S rRNA abundance In the ODZ region. In the free-living fraction, Thaumarcheota, Marinimicrobia, Actinobacteria, Nitrospina and Chloroflexi were all abundant (Supplementary Figure S9). One of the groups that only appears in the ODZ were Zetaproteobacteria (Figure 7). Although present in low numbers in the ODZ, Zetaproteobacteria were almost absent in the oxic part of the water column. Known Zetaproteobacteria are all autotrophic iron oxidizers, generally found at hydrothermal vents and in the sediments (McAllister et al., 2019). Reduced iron is often associated with ODZs (Kondo and Moffett, 2013; Vedamati et al., 2014). Additionally, some other groups like Planctomycetes and Nitrospina had unusually high abundances in the ODZ water column (Figure 6 and Supplementary Figure S7 respectively). Brocadia, or anammox bacteria, were the most abundant Planctomycetes in the ODZ though Planctomycetacia and Phycisphaera were also abundant (Figure 10). Anammox bacteria are known to primarily reside in suboxic and anoxic waters (Dalsgaard et al., 2014). Inside the ODZ region, the nitrite oxidizer Nitrospina had high abundance in the free-living fraction with ~107 -108 copies/L (Supplementary Figure S7). This increase in Nitrospina abundance can be related to the presence of nitrite in the ODZ (Figure 2B). Additionally, more ubiquitous bacteria such as SAR11, which have specific ODZ ecotypes (Tsementzi et al., 2016), also had a maximum in the ODZ (Figure 9).
Previous work has also found increases in particle number (Cram et al., 2022), particle flux (Fuchsman et al., 2019a) and cell counts (Ulloa et al., 2012) in ODZ regions. Potential reasons hypothesized for these increases include higher rates of chemosynthesis, photosynthesis when ODZs are at shallower depths than seen here and active transport by zooplankton who migrate into ODZs to avoid predation (Cram et al., 2022; Fuchsman et al., 2019a; Ulloa et al., 2012).
The non-buoyant plume (NBP) had an unique microbial community structure
Although our sampling did not reach the ocean floor, where the hydrothermal vents are located in the EPR, we analyzed a sample at 2369 m depth, 163 m from the bottom, which was characterized using beam transmission as containing particles from the hydrothermal vent’s non-buoyant plume (Supplementary Figure S2). This sample had a unique structure, where the 16S rRNA abundance did not follow the trends or patterns we observed through the water column at shallower depths. For example for the free living fraction the 16S rRNA abundance of bacteria in the Non-Buoyant plume region was similar to the numbers in the ODZ (~109 copies/L; Figure 5), but for the 53–180 and 180–500 μm size fractions the abundance was ~109 copies/L in the ODZ, but 106 copies/L at the Non-Buoyant plume. In the 53–180 and 180–500 μm size fractions, many of the typical bacteria, such as Planctomycetes, were missing and the number of Proteobacteria 16S rRNA copies were greatly decreased (Figure 6B, Supplementary Figure S9). In these size fractions, Pseudomonodales, rather than Alteromonadales, were the dominant Gammaproteobacteria (Figure 8). Non-Buoyant plumes evolve as they move away from the hydrothermal vent, with fewer vent related Campylobacterota and more SAR11, S-oxidizing Thioglobus Gammaproteobacteria and S-oxidizing Sulfurimonas Campylobacterota farther from the vent (Li et al., 2020). In our samples, in the Non-Buoyant plume, we found Thioglobus Gammaproteobacteria (Supplementary Figure S8), but Campylobacterota were almost absent in all of our samples, likely due to PCR biases. We observed SAR11 in our free-living non-buoyant plume sample (2,369 m depth) in high abundance (108 copies/L), but it was much less abundant in in 53–180 and 180–500 μm plume particle samples (103 copies/L) (Figure 9). The phylum Acidobacteria was also abundant in the free-living plume samples, but was not abundant at other depths (Supplementary Figure S9).
Methodological advances and caveats
This work includes three main technological advances. First, the use of 16S rRNA gene sequence spike-in allowed us to look quantitatively into microbial abundance across different particle sizes at different ocean depth ranges. This approach allowed us to compare microbial abundance of less abundant taxa directly across particle samples, rather than using relative abundance which is greatly affected by the abundance of the most abundant microbes. Secondly, we describe the microbial diversity and abundance in marine particles separated according to seven size-fractions from 0.2 to > 500 μm. This is many more particle sizes compared to other studies in ODZs that have only focused into 2 size fractions: 0.2–1.6 μm and >1.6 μm (Ganesh et al., 2014) or 0.2–30 μm and >30 μm (Fuchsman et al., 2017). This work has allowed us to examine how microbial communities change across the full particle size spectrum, and consider small, intermediate and large particles separately. Thirdly, we also examine microbial abundance quantitatively in relation to organic matter concentrations.
However, there are several caveats to this work. We learned, based on the non-detection of spike in some samples, and the low detection of spike in others that 104 spikes per ng of DNA is not enough spike copies. We therefore recommend using at least 105 spikes per ng of environmental DNA if using an illumina MiSeq or similar instrumentation. Greater sequence depth may allow use of similar or lower spike abundance. All PCR based environmental sequencing approaches including this one are subject to amplification biases which under-represent some species and over-represent others (Elbrecht and Leese, 2015). For example, the primers used here are known to not amplify Campylobacterota (McNichol et al., 2021), which are important in particle systems (Boeuf et al., 2019; Leu et al., 2022). Surprisingly, while McNichol et al. (2021) identified a primer mismatch to the order Brocadiales, which include anammox bacteria, we identified many Brocadiales bacteria in this dataset (Figure 10), suggesting that the mismatch may still allow amplification of this group. However, such a mismatch might lead to under amplification – it is possible that Brocadiales are even more abundant than indicated herein. Additionally, 16S rRNA gene corresponds to identity not function; in previous particle metagenomes, there were found several genes involved in denitrification (nosZ, nirK, qnorB) enriched in particles (Fuchsman et al., 2017; Ganesh et al., 2014). However, denitrifiers are phylogenetically diverse and it is very difficult to identify them based on 16S rRNA gene sequences (Bertagnolli et al., 2020; Zhang et al., 2023; Zumft, 1997). Another confounding factor is that the 16S rRNA gene is not a single copy core gene. Many bacteria contain multiple copies of the 16S rRNA gene (1–15) (Espejo and Plaza, 2018; Venter et al., 2004). Multiple gene copies per cell would cause our gene copy counts to be higher than true microbial cellular abundances. Indeed in the Chesapeake Bay, the amplicon-based estimates of microbial abundance were generally higher than microscopy counts from the particles (Cram et al., 2024). We advocate similar quantitative and size fractionated work with metagenomes, which would eliminate PCR amplification biases, allow the examination of functional genes for denitrification and sulfate reduction, and would allow the use of known single copy core genes for microbial enumeration. Despite these limitations, the patterns in 16S rRNA gene abundances here should reflect the patterns in their host organisms, allowing us to examine patterns in microbial abundance across a full size spectrum and in relation to organic C in the open ocean for the first time.
Concluding remarks
Combining 16S rRNA gene standard spike in amplicon sequencing and fractionation of particles at a finer scale allowed us to look into the diversity and abundance of free-living and particle-attached microbes from a new perspective. When we compared the quantitative abundance of different microbial taxa across depths and particle sizes, we observed increases in abundance of many groups in the ODZ region and in the non-buoyant hydrothermal vent plume, both regions that promote chemosynthesis.
There was a large increase in the bacterial 16S rRNA copies/ug C and copies/L on particles >180 μm between the upper water column and the deep water column. We expect that the density of bacteria on large particles was increasing with depth. Quantitative abundance estimates found that ostensibly particle attached bacteria such as Bacteroidetes and Planctomycetes were most abundant on particles, but still had many members in the free living state. Such patterns would not have been evident if we had looked at relative abundance data only, where the changes in total numbers of microbes between the water and particles are not taken into account. The significant presence of particle attached bacteria in the free living phase is consistent with the published models that assume that microbes are constantly colonizing and detaching from particles as they sink, allowing the microbes to operate in the pressure zone where they are adapted (Nguyen et al., 2022).
Indeed, models of microorganisms on particles have provided a valuable description of the microbial attachment, detachment and growth dynamics on particles (Mislan et al., 2014; Nguyen et al., 2022; Stephens et al., 2024). However, these models are limited by existing data in two ways. First, they have either not included observations of microbial diversity (Mislan et al., 2014; Nguyen et al., 2022) or else compared their model results to relative abundances of observed marine microbes on sinking particles (Stephens et al., 2024). Secondly, considering the continuum of particle sizes provides valuable insights into microbial processes (Cram et al., 2022, 2018; Leung et al., 2021; Weber and Bianchi, 2020) but microbially resolved models have not to date considered the full size spectrum. In fact, this work is the first time microbial abundance and particle organic C concentrations have been measured together across the size spectrum in the ocean. These data pave the way for more quantitative models of microbial processes on particles of a range of sizes.
Data availability statement
The sequence data are available in NCBI’s Short read archive under BioProject PRJNA1191024. CTD data are in the R2R Rolling deck repository (https://www.rvdata.us/search/cruise/AT42-09). Nutrient data have been deposited in BCO-DMO (doi: 10.26008/1912/bco-dmo.948718.1). Organic C concentrations have been deposited in BCO-DMO (doi: 10.26008/1912/bco-dmo.948709.1). A table with the abundances of each ASV has been deposited at FigShare (doi: 10.6084/m9.figshare.28887038.v1).
Author contributions
PH-V: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review & editing. CF: Conceptualization, Resources, Methodology, Supervision, Writing – original draft, Writing – review & editing. BT: Conceptualization, Investigation, Methodology, Writing – review & editing. JS: Funding acquisition, Resources, Writing – review & editing. JC: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Resources, Software, Supervision, Validation, Visualization, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. Funding for publication of this manuscript was provided by NSF OCE-2421095 to JAC. The cruise was funded by NSF OCE-1756339 to JS. BT was supported by the Center for Dark Energy Biosphere Investigations (C-DEBI) through NSF OCE-0939654.
Acknowledgments
We thank Ashley Hollins for help with DNA extractions. We thank Captain and Crew from the R/V Atlantis as well as members of the scientific crew who aided in the sampling effort. This paper was a chapter of PHV’s thesis. We appreciate the comments from her committee on the original thesis chapter. We acknowledge funding from Horn Point Laboratory startup and student funds.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2025.1552305/full#supplementary-material
References
Abramson, L., Lee, C., Liu, Z., Wakeham, S. G., and Szlosek, J. (2010). Exchange between suspended and sinking particles in the Northwest Mediterranean as inferred from the organic composition of in situ pump and sediment trap samples. Limnol. Oceanogr. 55, 725–739. doi: 10.4319/lo.2009.55.2.0725
Alldredge, A. (1998). The carbon, nitrogen and mass content of marine snow as a function of aggregate size. Deep-Sea Res. Part Oceanogr. Res. Pap. 45, 529–541. doi: 10.1016/S0967-0637(97)00048-4
Altabet, M. A. (1988). Variations in nitrogen isotopic composition between sinking and suspended particles: implications for nitrogen cycling and particle transformation in the open ocean. Deep Sea res. Part Oceanogr. Res. Pap. 35, 535–554. doi: 10.1016/0198-0149(88)90130-6
Amano, C., Zhao, Z., Sintes, E., Reinthaler, T., Stefanschitz, J., Kisadur, M., et al. (2022). Limited carbon cycling due to high-pressure effects on the deep-sea microbiome. Nat. Geosci. 15, 1041–1047. doi: 10.1038/s41561-022-01081-3
Babbin, A. R., Keil, R. G., Devol, A. H., and Ward, B. B. (2014). Organic matter stoichiometry, flux, and oxygen control nitrogen loss in the ocean. Science 344, 406–408. doi: 10.1126/science.1248364
Banse, K., Naqvi, S. W. A., and Postel, J. R. (2017). A zona incognita surrounds the secondary nitrite maximum in open-ocean oxygen minimum zones. Deep-Sea Res. Part Oceanogr. Res. Pap. 127, 111–113. doi: 10.1016/j.dsr.2017.07.004
Bertagnolli, A. D., Konstantinidis, K. T., and Stewart, F. J. (2020). Non-denitrifier nitrous oxide reductases dominate marine biomes. Environ. Microbiol. Rep. 12, 681–692. doi: 10.1111/1758-2229.12879
Bianchi, D., Weber, T. S., Kiko, R., and Deutsch, C. (2018). Global niche of marine anaerobic metabolisms expanded by particle microenvironments. Nat. Geosci. 11, 263–268. doi: 10.1038/s41561-018-0081-0
Bižić-Ionescu, M., Ionescu, D., and Grossart, H.-P. (2018). Organic particles: heterogeneous hubs for microbial interactions in aquatic ecosystems. Front. Microbiol. 9:2569. doi: 10.3389/fmicb.2018.02569
Boeuf, D., Edwards, B. R., Eppley, J. M., Hu, S. K., Poff, K. E., Romano, A. E., et al. (2019). Biological composition and microbial dynamics of sinking particulate organic matter at abyssal depths in the oligotrophic open ocean. Proc. Natl. Acad. Sci. USA 116, 11824–11832. doi: 10.1073/pnas.1903080116
Bristow, L. A., Callbeck, C. M., Larsen, M., Altabet, M. A., Dekaezemacker, J., Forth, M., et al. (2017). N2 production rates limited by nitrite availability in the bay of Bengal oxygen minimum zone. Nat. Geosci. 10, 24–29. doi: 10.1038/NGEO2847
Bristow, L. A., Dalsgaard, T., Tiano, L., Mills, D. B., Bertagnolli, A. D., Wright, J. J., et al. (2016). Ammonium and nitrite oxidation at nanomolar oxygen concentrations in oxygen minimum zone waters. Proc. Natl. Acad. Sci. USA 113, 10601–10606. doi: 10.1073/pnas.1600359113
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869
Capovilla, G., Braakman, R., Fournier, G. P., Hackl, T., Schwartzman, J., Lu, X., et al. (2023). Chitin utilization by marine picocyanobacteria and the evolution of a planktonic lifestyle. Proc. Natl. Acad. Sci. 120, –e2213271120. doi: 10.1073/pnas
Cavan, E. L., Giering, S. L. C., Wolff, G. A., Trimmer, M., and Sanders, R. (2018). Alternative particle formation pathways in the eastern tropical North Pacific’s biological carbon pump. J. Geophys. Res. Biogeosciences 123, 2198–2211. doi: 10.1029/2018JG004392
Chang, B. X., Devol, A. H., and Emerson, S. R. (2012). Fixed nitrogen loss from the eastern tropical North Pacific and Arabian Sea oxygen deficient zones determined from measurements of N2:Ar. Glob. Biogeochem. Cycles 26:GB004207. doi: 10.1029/2011GB004207
Clegg, S. L., and Whitfield, M. (1990). A generalized model for the scavenging of trace metals in the open ocean—I. Particle cycling. Deep Sea res. Part Oceanogr. Res. Pap. 37, 809–832. doi: 10.1016/0198-0149(90)90008-J
Cram, J. A., Fuchsman, C. A., Duffy, M. E., Pretty, J. L., Lekanoff, R. M., Neibauer, J. A., et al. (2022). Slow particle remineralization, rather than suppressed disaggregation, drives efficient flux transfer through the eastern tropical North Pacific oxygen deficient zone. Glob. Biogeochem. Cycles 36:e2021GB007080. doi: 10.1029/2021GB007080
Cram, J. A., Hollins, A., McCarty, A. J., Martinez, G., Cui, M., Gomes, M., et al. (2024). Microbial diversity and abundance vary along salinity, oxygen and particle size gradients in the Chesapeake Bay. Environ. Microbiol. 26:e16557. doi: 10.1111/1462-2920.16557
Cram, J. A., Weber, T., Leung, S. W., McDonnell, A. M. P., Liang, J. H., and Deutsch, C. (2018). The role of particle size, ballast, temperature, and oxygen in the sinking flux to the Deep Sea. Glob. Biogeochem. Cycles 32, 858–876. doi: 10.1029/2017GB005710
Dadon-Pilosof, A., Conley, K. R., Jacobi, Y., Haber, M., Lombard, F., Sutherland, K. R., et al. (2017). Surface properties of SAR11 bacteria facilitate grazing avoidance. Nat. Microbiol. 2, 1608–1615. doi: 10.1038/s41564-017-0030-5
Dalsgaard, T., Stewart, F. J., Thamdrup, B., Brabandere, L. D., Revsbech, P., and Ulloa, O. (2014). Oxygen at Nanomolar levels reversibly suppresses process rates and gene expression in Anammox and denitrification in the oxygen minimum zone off northern Chile. MBio 5, e01966–e01914. doi: 10.1128/mBio.01966-14
Dalsgaard, T., Thamdrup, B., Farías, L., and Peter Revsbech, N. (2012). Anammox and denitrification in the oxygen minimum zone of the eastern South Pacific. Limnol. Oceanogr. 57, 1331–1346. doi: 10.4319/lo.2012.57.5.1331
Datta, M. S., Sliwerska, E., Gore, J., Polz, M. F., and Cordero, O. X. (2016). Microbial interactions lead to rapid micro-scale successions on model marine particles. Nat. Commun. 7:11965. doi: 10.1038/ncomms11965
Delong, E. F., Franks, D. G., and Alldredge, A. L. (1993). Phylogenetic diversity of aggregate-attached marine bacterial assemblages. Limnol. Oceanogr. 38, 924–934. doi: 10.4319/lo.1993.38.5.0924
DeVries, T., Deutsch, C., Rafter, P. A., and Primeau, F. (2013). Marine denitrification rates determined from a global 3-D inverse model. Biogeosciences 10, 2481–2496. doi: 10.5194/bg-10-2481-2013
Dick, G. J. (2019). The microbiomes of deep-sea hydrothermal vents: distributed globally, shaped locally. Nat. Rev. Microbiol. 17, 271–283. doi: 10.1038/s41579-019-0160-2
Dougherty, E. C., Hollins, A., and Cram, J. A. (2025). Declines of suspended particle abundances in the oxygen-deficient zone of the Chesapeake Bay. Estuar. Coasts 48:29. doi: 10.1007/s12237-024-01464-2
Durkin, C. A., Cetinić, I., Estapa, M., Ljubešić, Z., Mucko, M., Neeley, A., et al. (2022). Tracing the path of carbon export in the ocean through DNA sequencing of individual sinking particles. ISME J. 16, 1896–1906. doi: 10.1038/s41396-022-01239-2
Durkin, C. A., Estapa, M. L., and Buesseler, K. O. (2015). Observations of carbon export by small sinking particles in the upper mesopelagic. Mar. Chem. 175, 72–81. doi: 10.1016/j.marchem.2015.02.011
Elbrecht, V., and Leese, F. (2015). Can DNA-based ecosystem assessments quantify species abundance? Testing primer bias and biomass-sequence relationships with an innovative metabarcoding protocol. PLoS One 10, –e0130324. doi: 10.1371/journal.pone.0130324
Eloe, E. A., Shulse, C. N., Fadrosh, D. W., Williamson, S. J., Allen, E. E., and Bartlett, D. H. (2011). Compositional differences in particle-associated and free-living microbial assemblages from an extreme deep-ocean environment. Environ. Microbiol. Rep. 3, 449–458. doi: 10.1111/j.1758-2229.2010.00223.x
Espejo, R. T., and Plaza, N. (2018). Multiple ribosomal RNA operons in Bacteria; their concerted evolution and potential consequences on the rate of evolution of their 16S rRNA. Front. Microbiol. 9:1232. doi: 10.3389/fmicb.2018.01232
Evans, N., Tichota, J., Ruef, W., Moffett, J., and Devol, A. (2023). Rapid expansion of fixed nitrogen deficit in the eastern Pacific Ocean revealed by 50-year time series. Glob. Biogeochem. Cycles 37:e2022GB007575. doi: 10.1029/2022GB007575
Fu, S., Wang, R., Xu, Z., Zhou, H., Qiu, Z., Shen, L., et al. (2023). Metagenomic sequencing combined with flow cytometry facilitated a novel microbial risk assessment framework for bacterial pathogens in municipal wastewater without cultivation. iMeta 2:e77. doi: 10.1002/imt2.77
Fuchsman, C. A., and Cram, J. A. (2024). Size fractionated suspended organic carbon and nitrogen from the offshore eastern tropical North Pacific oxygen deficient zone suggest contributions of picocyanobacteria and vertically migrating metazoans to organic matter. ESS Open Arch. doi: 10.22541/essoar.173046855.50289201/v1
Fuchsman, C. A., Devol, A. H., Casciotti, K. L., Buchwald, C., Chang, B. X., and Horak, R. E. A. (2018). An N isotopic mass balance of the eastern tropical North Pacific oxygen deficient zone. Deep-Sea res. Part II Top. Stud. Oceanogr. 156, 137–147. doi: 10.1016/j.dsr2.2017.12.013
Fuchsman, C. A., Devol, A. H., Saunders, J. K., McKay, C., and Rocap, G. (2017). Niche partitioning of the N cycling microbial community of an offshore oxygen deficient zone. Front. Microbiol. 8:2384. doi: 10.3389/fmicb.2017.02384
Fuchsman, C. A., Kirkpatrick, J. B., Brazelton, W. J., Murray, J. W., and Staley, J. T. (2011). Metabolic strategies of free-living and aggregate-associated bacterial communities inferred from biologic and chemical profiles in the Black Sea suboxic zone. FEMS Microbiol. Ecol. 78, 586–603. doi: 10.1111/j.1574-6941.2011.01189.x
Fuchsman, C. A., Palevsky, H. I., Widner, B., Duffy, M., Carlson, M. C. G., Neibauer, J. A., et al. (2019a). Cyanobacteria and cyanophage contributions to carbon and nitrogen cycling in an oligotrophic oxygen-deficient zone. ISME J. 13, 2714–2726. doi: 10.1038/s41396-019-0452-6
Fuchsman, C. A., Paul, B., Staley, J. T., Yakushev, E. V., and Murray, J. W. (2019b). Detection of transient denitrification during a high organic matter event in the Black Sea. Glob. Biogeochem. Cycles 33, 143–162. doi: 10.1029/2018GB006032
Fuchsman, C. A., Staley, J. T., Oakley, B. B., Kirkpatrick, J. B., and Murray, J. W. (2012). Free-living and aggregate-associated Planctomycetes in the Black Sea. FEMS Microbiol. Ecol. 80, 402–416. doi: 10.1111/j.1574-6941.2012.01306.x
Fuerst, J. A., and Sagulenko, E. (2011). Beyond the bacterium: Planctomycetes challenge our concepts of microbial structure and function. Nat. Rev. Microbiol. 9, 403–413. doi: 10.1038/nrmicro2578
Ganesh, S., Parris, D. J., DeLong, E. F., and Stewart, F. J. (2014). Metagenomic analysis of size-fractionated picoplankton in a marine oxygen minimum zone. ISME J. 8, 187–211. doi: 10.1038/ismej.2013.144
Giovannoni, S. J. (2017). SAR11 Bacteria: the Most abundant plankton in the oceans. Annu. Rev. Mar. Sci. 9, 231–255. doi: 10.1146/annurev-marine-010814-015934
Gloor, G. B., Wu, J. R., Pawlowsky-Glahn, V., and Egozcue, J. J. (2016). It’s all relative: analyzing microbiome data as compositions. Ann. Epidemiol. 26, 322–329. doi: 10.1016/j.annepidem.2016.03.003
Guidi, L., Stemmann, L., Jackson, G. A., Ibanez, F., Claustre, H., Legendre, L., et al. (2009). Effects of phytoplankton community on production, size and export of large aggregates: a world-ocean analysis. Limnol. Oceanogr. 54, 1951–1963. doi: 10.4319/lo.2009.54.6.1951
Haro-Moreno, J. M., Rodriguez-Valera, F., Rosselli, R., Martinez-Hernandez, F., Roda-Garcia, J. J., Lluesma Gomez, M., et al. (2020). Ecogenomics of the SAR11 clade. Environ. Microbiol. 22, 1748–1763. doi: 10.1111/1462-2920.14896
Ijichi, M., Itoh, H., and Hamasaki, K. (2019). Vertical distribution of particle-associated and free-living ammonia-oxidizing archaea in Suruga bay, a deep coastal embayment of Japan. Arch. Microbiol. 201, 1141–1146. doi: 10.1007/s00203-019-01680-6
Ito, T., Minobe, S., Long, M. C., and Deutsch, C. (2017). Upper ocean O2 trends: 1958–2015. Geophys. Res. Lett. 44, 4214–4223. doi: 10.1002/2017GL073613
Jin, X., Gruber, N., Dunne, J. P., Sarmiento, J. L., and Armstrong, R. A. (2006). Diagnosing the contribution of phytoplankton functional groups to the production and export of particulate organic carbon, CaCO3, and opal from global nutrient and alkalinity distributions. Glob. Biogeochem. Cycles 20:2005GB002532. doi: 10.1029/2005GB002532
Jones-Kellett, A. E., McNichol, J. C., Raut, Y., Cain, K. R., Ribalet, F., Armbrust, E. V., et al. (2024). Amplicon sequencing with internal standards yields accurate picocyanobacteria cell abundances as validated with flow cytometry. ISME Commun. 4:ycae115. doi: 10.1093/ismeco/ycae115
Keil, R. G., Neibauer, J. A., Biladeau, C., Van Der Elst, K., and Devol, A. H. (2016). A multiproxy approach to understanding the “enhanced” flux of organic matter through the oxygen-deficient waters of the Arabian Sea. Biogeosciences 13, 2077–2092. doi: 10.5194/bg-13-2077-2016
Kellogg, C. T. E., and Deming, J. W. (2009). Comparison of free-living, suspended particle, and aggregate-associated bacterial and archaeal communities in the Laptev Sea. Aquat. Microb. Ecol. 57, 1–18. doi: 10.3354/ame01317
Kondo, Y., and Moffett, J. W. (2013). Dissolved Fe (II) in the Arabian Sea oxygen minimum zone and western tropical Indian Ocean during the inter-monsoon period. Deep-Sea Res. Part I 73, 73–83. doi: 10.1016/j.dsr.2012.11.014
Kwiecinski, J. V., and Babbin, A. R. (2021). A high-resolution atlas of the eastern tropical Pacific oxygen deficient zones. Glob. Biogeochem. Cycles 35:e2021GB007001. doi: 10.1029/2021GB007001
Lam, P., Lavik, G., Jensen, M. M., van de Vossenberg, J., Schmid, M., Woebken, D., et al. (2009). Revising the nitrogen cycle in the Peruvian oxygen minimum zone. Proc. Natl. Acad. Sci. USA 106, 4752–4757. doi: 10.1073/pnas.0812444106
Le Moigne, F. A. C. (2019). Pathways of organic carbon downward transport by the oceanic biological carbon pump. Front. Mar. Sci. 6:634. doi: 10.3389/fmars.2019.00634
LeCleir, G. R., DeBruyn, J. M., Maas, E. W., Boyd, P. W., and Wilhelm, S. W. (2014). Temporal changes in particle-associated microbial communities after interception by nonlethal sediment traps. FEMS Microbiol. Ecol. 87, 153–163. doi: 10.1111/1574-6941.12213
Leu, A. O., Eppley, J. M., Burger, A., and Delong, E. F. (2022). Diverse genomic traits differentiate sinking-particle-associated versus free-living microbes throughout the oligotrophic Open Ocean water column. MBio 13:e0156922. doi: 10.1128/mbio.01569-22
Leung, S. W., Weber, T., Cram, J. A., and Deutsch, C. (2021). Variable particle size distributions reduce the sensitivity of global export flux to climate change. Biogeosciences 18, 229–250. doi: 10.5194/bg-18-229-2021
Li, J., Yang, J., Sun, M., Su, L., Wang, H., Gao, J., et al. (2020). Distribution and succession of microbial communities along the dispersal pathway of hydrothermal plumes on the southwest Indian ridge. Front. Mar. Sci. 7:581381. doi: 10.3389/fmars.2020.581381
Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet 17, 10–12. doi: 10.14806/ej.17.1.200
McAllister, S. M., Moore, R. M., Gartman, A., Luther, G. W., Emerson, D., and Chan, C. S. (2019). The Fe(II)-oxidizing Zetaproteobacteria: historical, ecological and genomic perspectives. FEMS Microbiol. Ecol. 95:fiz015. doi: 10.1093/femsec/fiz015
McNichol, J., Berube, P. M., Biller, S. J., and Fuhrman, J. A. (2021). Evaluating and improving small subunit rRNA PCR primer coverage for Bacteria, Archaea, and eukaryotes using metagenomes from Global Ocean surveys. mSystems 6:e0056521. doi: 10.1128/mSystems.00565-21
Mestre, M., Borrull, E., Sala, M. M., and Gasol, J. M. (2017). Patterns of bacterial diversity in the marine planktonic particulate matter continuum. ISME J. 11, 999–1010. doi: 10.1038/ismej.2016.166
Mestre, M., Höfer, J., Sala, M. M., and Gasol, J. M. (2020). Seasonal variation of bacterial diversity along the marine particulate matter continuum. Front. Microbiol. 11:1590. doi: 10.3389/fmicb.2020.01590
Mestre, M., Ruiz-González, C., Logares, R., Duarte, C. M., Gasol, J. M., and Sala, M. M. (2018). Sinking particles promote vertical connectivity in the ocean microbiome. Proc. Natl. Acad. Sci. 115, E6799–E6807. doi: 10.1073/pnas.1802470115
Mislan, K. A. S., Stock, C. A., Dunne, J. P., and Sarmiento, J. L. (2014). Group behavior among model bacteria influences particulate carbon remineralization depths. J. Mar. Res. 72, 183–218. doi: 10.1357/002224014814901985
Nguyen, R. T., and Harvey, H. R. (1997). Protein and amino acid cycling during phytoplankton decomposition in oxic and anoxic waters. Org. Geochem. 27, 115–128. doi: 10.1016/S0146-6380(97)00076-4
Nguyen, T. T. H., Zakem, E. J., Ebrahimi, A., Schwartzman, J., Caglar, T., Amarnath, K., et al. (2022). Microbes contribute to setting the ocean carbon flux by altering the fate of sinking particulates. Nat. Commun. 13:1657. doi: 10.1038/s41467-022-29297-2
Parada, A. E., Needham, D. M., and Fuhrman, J. A. (2016). Every base matters: assessing small subunit RRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ. Microbiol. 18, 1403–1414. doi: 10.1111/1462-2920.13023
Peters, B., Horak, R., Devol, A. H., Fuchsman, C. A., Forbes, M., Mordy, C. W., et al. (2018). Estimating fixed nitrogen loss and associated isotope effects using concentration and isotope measurements of NO3–, NO2–, and N2 from the eastern tropical South Pacific oxygen deficient zone. Deep-Sea res. Part II Top. Stud. Oceanogr. 156, 121–136. doi: 10.1016/j.dsr2.2018.02.011
Ploug, H., Kiihl, M., Buchholz-Cleven, B., and Jergensen, B. B. (1997). Anoxic aggregates- an ephemeral phenomenon in the pelagic environment? Aquat. Microb. Ecol. 13, 285–294. doi: 10.3354/ame013285
Poff, K. E., Leu, A. O., Eppley, J. M., Karl, D. M., and DeLong, E. F. (2021). Microbial dynamics of elevated carbon flux in the open ocean’s abyss. Proc. Natl. Acad. Sci. USA 118:e2018269118. doi: 10.1073/pnas.2018269118
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Raven, M. R., Keil, R. G., and Webb, S. M. (2021). Microbial sulfate reduction and organic sulfur formation in sinking marine particles. Science 371, 178–181. doi: 10.1126/science.abc6035
Reintjes, G., Heins, A., Wang, C., and Amann, R. (2023). Abundance and composition of particles and their attached microbiomes along an Atlantic meridional transect. Front. Mar. Sci. 10:1051510. doi: 10.3389/fmars.2023.1051510
Revsbech, N. P., Larsen, L. H., Gundersen, J., Dalsgaard, T., Ulloa, O., and Thamdrup, B. (2009). Determination of ultra-low oxygen concentrations in oxygen minimum zones by the STOX sensor. Limnol. Oceanogr. 7, 371–381. doi: 10.4319/lom.2009.7.371
Santoro, A. E., Saito, M. A., Goepfert, T. J., Lamborg, C. H., Dupont, C. L., and DiTullio, G. R. (2017). Thaumarchaeal ecotype distributions across the equatorial Pacific Ocean and their potential roles in nitrification and sinking flux attenuation. Limnol. Oceanogr. 62, 1984–2003. doi: 10.1002/lno.10547
Saunders, J. K., Fuchsman, C. A., Mckay, C., and Rocap, G. (2019). Complete arsenic-based respiratory cycle in the marine microbial communities of pelagic oxygen-deficient zones. Proc. Natl. Acad. Sci. USA 116, 9925–9930. doi: 10.1073/pnas.1818349116
Shelton, A. O., Gold, Z. J., Jensen, A. J., Andruszkiewicz Allan, E., Van Cise, A., Gallego, R., et al. (2023). Toward quantitative metabarcoding. Ecology 104:e3906. doi: 10.1002/ecy.3906
Simon, M., Grossart, H. P., Schweitzer, B., and Ploug, H. (2002). Microbial ecology of organic aggregates in aquatic ecosystems. Aquat. Microb. Ecol. 28, 175–211. doi: 10.3354/ame028175
Stephens, B. M., Durkin, C. A., Sharpe, G., Nguyen, T. T. H., Albers, J., Estapa, M. L., et al. (2024). Direct observations of microbial community succession on sinking marine particles. ISME J. 18:wrad010. doi: 10.1093/ismejo/wrad010
Stramma, L., Johnson, G. C., Sprintall, J., and Mohrholz, V. (2008). Expanding oxygen-minimum zones in the tropical oceans. Science 320, 655–658. doi: 10.1126/science.1153847
Sun, Z. Z., Ji, B. W., Zheng, N., Wang, M., Cao, Y., Wan, L., et al. (2021). Phylogenetic distribution of polysaccharide-degrading enzymes in marine Bacteria. Front. Microbiol. 12:658620. doi: 10.3389/fmicb.2021.658620
Suter, E. A., Scranton, M. I., Chow, S., Stinton, D., Medina Faull, L., and Taylor, G. T. (2017). Niskin bottle sample collection aliases microbial community composition and biogeochemical interpretation. Limnol. Oceanogr. 62, 606–617. doi: 10.1002/lno.10447
Sylvan, J. B., Pyenson, B. C., Rouxel, O., German, C. R., and Edwards, K. J. (2012). Time-series analysis of two hydrothermal plumes at 9°50′N East Pacific rise reveals distinct, heterogeneous bacterial populations. Geobiology 10, 178–192. doi: 10.1111/j.1472-4669.2011.00315.x
Tamelander, T. (2013). Community composition and extracellular enzyme activity of bacteria associated with suspended and sinking particles in contrasting Arctic and sub-arctic marine environments. Aquat. Microb. Ecol. 69, 211–221. doi: 10.3354/ame01641
Tettamanti Boshier, F. A., Srinivasan, S., Lopez, A., Hoffman, N. G., Proll, S., Fredricks, D. N., et al. (2020). Complementing 16S rRNA gene amplicon sequencing with Total bacterial load to infer absolute species concentrations in the vaginal microbiome. mSystems 5, e00777–e00719. doi: 10.1128/mSystems.00777-19
Tiano, L., Garcia-Robledo, E., Dalsgaard, T., Devol, A. H., Ward, B. B., Ulloa, O., et al. (2014). Oxygen distribution and aerobic respiration in the north and south eastern tropical Pacific oxygen minimum zones. Deep Sea res. Part Oceanogr. Res. Pap. 94, 173–183. doi: 10.1016/j.dsr.2014.10.001
Tourlousse, D. M., Yoshiike, S., Ohashi, A., Matsukura, S., Noda, N., and Sekiguchi, Y. (2017). Synthetic spike-in standards for high-throughput 16S rRNA gene amplicon sequencing. Nucleic Acids Res. 45, gkw984–e23. doi: 10.1093/nar/gkw984
Tsementzi, D., Wu, J., Deutsch, S., Nath, S., Rodriguez-R, L. M., Burns, A. S., et al. (2016). SAR11 bacteria linked to ocean anoxia and nitrogen loss. Nature 536, 179–183. doi: 10.1038/nature19068
Ulloa, O., Delong, E. F., Letelier, R. M., and Stewart, F. J. (2012). Microbial oceanography of anoxic oxygen minimum zones. Proc. Natl. Acad. Sci. 109, 15996–16003. doi: 10.1073/pnas.1205009109
Van Mooy, B. A. S., Keil, R. G., and Devol, A. H. (2002). Impact of suboxia on sinking particulate organic carbon: enhanced carbon flux and preferential degradation of amino acids via denitrification. Geochim. Cosmochim. Acta 66, 457–465. doi: 10.1016/S0016-7037(01)00787-6
Vedamati, J., Goepfert, T., and Moffett, J. W. (2014). Iron speciation in the eastern tropical South Pacific oxygen minimum zone off Peru. Limnol. Oceanogr. 59, 1945–1957. doi: 10.4319/lo.2014.59.6.1945
Venter, J. C., Remington, K., Heidelberg, J. F., Halpern, A. L., Rusch, D., Eisen, J. A., et al. (2004). Environmental genome shotgun sequencing of the Sargasso Sea. Science 304, 66–74. doi: 10.1126/science.1093857
Wakeham, S. G., and Canuel, E. A. (1988). Organic geochemistry of particulate matter in the eastern tropical North Pacific Ocean: implications for particle dynamics. J. Mar. Res. 46, 183–213. doi: 10.1357/002224088785113748
Wang, S., Wu, Q., Han, Y., Du, R., Wang, X., Nie, Y., et al. (2021). Gradient internal standard method for absolute quantification of microbial amplicon sequencing data 6, e00964–e00920. doi: 10.1128/mSystems.00964-20
Weber, T., and Bianchi, D. (2020). Efficient particle transfer to depth in oxygen minimum zones of the Pacific and Indian oceans. Front. Earth Sci. 8:376. doi: 10.3389/feart.2020.00376
Yeh, Y. C., and Fuhrman, J. A. (2022). Effects of phytoplankton, viral communities, and warming on free-living and particle-associated marine prokaryotic community structure. Nat. Commun. 13:7905. doi: 10.1038/s41467-022-35551-4
Zhang, I. H., Fortin, S. G., Ward, B. B., and Babbin, A. R. (2023). Partitioning of the denitrification pathway and other nitrite metabolisms within global oxygen deficient zones. ISME Commun. 3:76. doi: 10.1038/s43705-023-00284-y
Keywords: marine aggregates, East Pacific Rise, oxygen deficient zone, size fractionation, organic matter, microbial communities, microbial diversity
Citation: Huanca-Valenzuela P, Fuchsman CA, Tully BJ, Sylvan JB and Cram JA (2025) Quantitative microbial taxonomy across particle size, depth, and oxygen concentration. Front. Microbiol. 16:1552305. doi: 10.3389/fmicb.2025.1552305
Edited by:
Abhishek Srivastava, University of Vienna, AustriaReviewed by:
Klaus Jürgens, Leibniz Institute for Baltic Sea Research (LG), GermanyJuan Antonio Garcia, University of Veterinary Medicine Vienna, Austria
Copyright © 2025 Huanca-Valenzuela, Fuchsman, Tully, Sylvan and Cram. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jacob A. Cram, amNyYW1AdW1jZXMuZWR1