Saturation Approach to Determine Grazing Mortality in Picoeukaryote and Synechococcus Populations

A substantial component of phytoplankton production in the oceans is channeled through protistan grazers but understanding what dictates the magnitude of this process on a regional and temporal basis is limited, in part, by a shortage of experimental options. A novel saturation approach based on the functional response of planktonic grazers to increasing prey abundance was developed using laboratory cultures of the predator-prey combination of Ochromonas danica and Micromonas pusilla and tested in the coastal waters of the Gulf of Maine. In incubation series, 2 μm polystyrene microspheres were used as surrogate prey to generate increasing levels of saturation of predator ingestion rates of natural prey, resulting in increased rates of apparent growth of the picophytoplankton populations. The relationship between level of addition of surrogate prey to apparent growth, consistently provided significant estimates of maximal growth in the absence of grazing and grazing mortality for populations of picoeukaryotes and Synechococcus. Estimates of gross growth and grazing mortality were comparable to results from dilution experiments carried out in the same waters. The saturation approach represents an additional tool to investigate predator-prey interactions in planktonic communities. Further investigations may show that it can be used to quantify group-specific grazing mortality and growth rates beyond coastal waters and in multiple size classes of prey.


INTRODUCTION
Herbivory is the major fate of primary production in the oceans and the major herbivores in most pelagic communities are protists (Sherr and Sherr, 2002). These primary consumers have been loosely referred to as the microzooplankton, with the growing realization that many protists are mixotrophic, able to combine phagotrophy with photoautotrophic production (Jones, 1997;Stoecker, 1998;Zubkov and Tarran, 2008). Predation by microzooplankton influences planktonic systems at many levels. This includes structuring the phenotypic and genetic composition of prokaryote and eukaryote communities; comprising a major route of nutrient transfer between functional groups and size classes, including a conduit from primary producers to higher trophic levels outside the microbial loop; and acting as an important nutrient recycling process, increasing the bioavailability of inorganic nutrients and trace metals to bacteria and photoautotrophs (reviewed in Sherr and Sherr, 2002;Worden et al., 2015).
A variety of approaches have been used to quantify rates of grazing by microzooplankton in the natural environment. This includes tracing the ingestion of labeled or surrogate prey by grazers. Fluorescent dyes (Sherr et al., 1987;Rublee and Gallegos, 1989;Putt, 1991) and radioisotopes (Lessard and Swift, 1985;Archer et al., 1996) have been used to label algal or bacterial prey. Surrogate prey used to visually quantify ingestion rates of protists have included fluorescent microspheres, wheat starch particles and bacteria containing green fluorescent protein (Pace and Bailiff, 1987;Kivi and Setälä, 1995;Fu et al., 2003). These approaches provide useful information on particle ingestion rates by specific components of pelagic microbial communities but tend to be labor intensive, and are complicated by protists able to select between particle types, issues with preservation of fully representative components of the community and extrapolation of single cell-specific rates to a community level process (Pace and Bailiff, 1987;Putt, 1991).
By far the most commonly applied method has been the dilution approach of Landry and Hassett (1982). The near ubiquitous use of the dilution approach has provided the opportunity to compile and compare measurements conducted in multiple studies covering many regions (Sherr and Sherr, 2002;Calbet and Landry, 2004;Schmoker et al., 2013;Chen et al., 2014). Schmoker et al. (2013) compiled 1,525 data points from dilution experiments and grouped them into oceanic provinces to examine regional variability in microzooplankton grazing mortality and phytoplankton growth. High variability in grazing impact on the phytoplankton community was found between regions, with median values ranging between 0.07 d −1 in southern polar regions and 0.49 d −1 in the sub-tropical Atlantic; this is equivalent to the consumption of 53 and 70% of primary production, respectively. Despite these attempts to synthesize the vast number of measurements that have been conducted, the environmental factors underlying regional differences in microzooplankton grazing pressure are not readily apparent. This may reflect differences in how effectively the dilution approach performs in varied environments and how nonlinear relationships between the level of dilution and apparent growth rates of phytoplankton are interpreted in the experiments (reviewed by Dolan and McKeon, 2005).
In situ monitoring of picoplankton populations, possible using shipboard continuous or autonomous submersible flow cytometers, provides data from which model-derived estimates of gross growth rates and mortality can be obtained without the need for experimental incubations or the manipulation of cell abundance (Sosik et al., 2003;Ribalet et al., 2015;Fowler et al., 2020). Picophytoplankton populations, which include cyanobacteria such as Synechococcus and Prochlorococcus and small eukaryotes including prymnesiophytes, prasinophytes pelagophytes and chrysophytes, make important contributions to the ecological characteristics and biochemical flux in both nearshore and oceanic waters (Jardillier et al., 2010;Fowler et al., 2020). Distinctly higher gross growth rates of picophytoplankton have been derived from the in situ data, compared to those obtained from dilution experiments (Fowler et al., 2020), emphasizing the need to better understand their growth and fate, and the appeal of additional approaches to quantify these rates. Here we introduce an experimental approach that may avoid some of the complexities of the dilution approach, while providing additional insights into predator-prey interactions among these smaller size classes in natural waters.
A fundamental feature of the behavior of planktonic grazers is a rapid increase in ingestion as prey abundance is increased, followed by saturation of uptake above some threshold value (Frost, 1972;Fenchel, 1980). Here we present an approach that exploits this functional response to determine estimates of the rate of grazing mortality imposed by microzooplankton on their prey. The fundamentals of the approach are illustrated using a single predator-prey combination in laboratory culture. This is followed by presentation of a series of experimental tests carried out in natural waters, including a comparison with the dilution approach. This series of experiments is used to illustrate that a relatively straightforward experimental design provides an additional experimental tool to estimate rates of growth and mortality due to microzooplankton grazing of specific components of natural phytoplankton communities.

Theoretical Basis of the Approach
For microzooplankton, rates of prey ingestion are influenced by the time involved in both prey capture and prey digestion. The former includes encounter rates, prey handling and prey avoidance and/or rejection, the latter depends on how quickly prey can be processed through the digestive system before the predator becomes satiated (Flynn and Mitra, 2016). Addition of surrogate prey may affect both aspects and, even if selected against by the predator, the additional handling time involved will reduce the ingestion rate of natural prey and hence their rates of grazing mortality.
Demonstrations of the functional response of protozoan and metazoan zooplankton involve measurement of ingestion rates of grazers at increasing concentrations of a single prey type, including in some cases, artificial prey such as latex beads (Fenchel, 1980;Jonsson, 1986;Verity, 1991). In a procedure that mimics determination of the functional response, the proposed saturation approach involves adding increasing abundances of surrogate prey to incubations of natural communities and measuring the observed growth response of the natural prey of interest. At increasing levels of surrogate prey addition, as grazers become "saturated" by handling and ingesting surrogate prey, the apparent growth rate (µ app ) of the natural prey increases and approaches a maximum value that represents their specific growth rate in the absence of grazing (µ max ). The difference between µ app with no addition of surrogate prey, effectively the net growth rate (µ net ), and µ max , provides an estimate of the level of mortality due to grazing (m) with the assumption that the specific growth rate of the prey is not affected by the addition of surrogate prey. Analagous to plant growth in response to increasing amounts of a limiting nutrient, the increasing growth rate of prey in response to increasing surrogate prey was represented by an asymptotic regression model, of similar form to the Mitscherlich response model (Mitscherlich, 1928). The model was used to derive values of µ net , µ max and m from the relationship between µ app and level of addition of surrogate prey, where C is the surrogate prey concentration and a is proportional to the relative rate of increase in µ app with increasing C. Note, u max is referred to as "gross growth rate" but does not account for other causes of mortality, such as viral infection. Model fitting was carried out using the non-linear least squares regression function (nls) in R (R Core Team, 2021).

Tests Using Single Prey and Predator Combinations
Proof-of-concept tests were performed using laboratory cultures of a single isolate representative of the picophytoplankton and a typical protistan predator. Specifically, the experiments used the mixotrophic chrysophyte Ochromonas danica (CCMP 1391), ∼ 6-9 µm in diameter, as a predator, and the generally abundant picoeukaryote Micromonas pusilla (CCMP 485), ∼ 1-2 µm in diameter, as prey. The cultures were maintained in sterile L1 media in 50 ml borosilicate glass test tubes in a 21 • C incubator with a 14 h light/10 h dark cycle and light levels of ∼ 90 µE m −2 s −1 . Prior to the start of the experiments, O. danica was transitioned from K media and rice to only sterile L1 media and then fed every other day on M. pusilla. The prey, M. pusilla was maintained in semi-continuous growth through regular transfers (2-3 days). In 3 experiments, similar ratios of M. pusilla (target 3 × 10 5 ml −1 ) and O. danica (target 2,000 ml −1 ) were generated in tubes containing a total volume of 40 ml. Fluorescent polystyrene microspheres (beads) of 2 µm in diameter (Fluoresbrite Plain YG microspheres, Polysciences, Inc., Warrington, PA), were used as surrogate prey.
To minimize clumping, the beads were blocked in a solution of 1% bovine serum albumin overnight then centrifuged for 5 min at 2,000 rpm, after which the pellet was resuspended in 0.2 µm filtered seawater. A new solution of beads was prepared at the start of each experiment. Duplicate tubes were prepared for 6 different surrogate prey concentrations ranging from 0 to 7.5 × 10 6 ml −1 . During the incubations, the experimental tubes were rotated at 1.2 rpm on a plankton wheel to keep particles in suspension. Two subsamples of 1 ml were removed from each tube at the start (T 0 ) and the final time point (T F ) after ∼24 h, for flow cytometric analysis.

Picophytoplankton Rates in Natural Waters
The experimental format established in the culture experiments was then adapted for use in natural waters. From 16th July to 15th August, 2019 (days of the year 197-227) and 8th June to 26th July 2021 (days 165-207), a series of saturation experiments were performed that focused on determining the growth rates and grazing mortality of picoeukaryote and Syncechococcus sp. populations. Gulf of Maine coastal seawater was collected from the Damariscotta River Estuary at Bigelow Laboratory's dock off East Boothbay, Maine, United States (43.8604 • N, 69.5781 • W). Water for all experiments was collected within 1 h of high tide at 1 m depth using a 5 l Niskin bottle and was gravity filtered through 200 µm mesh to remove zooplankton. Several Niskin casts of seawater were combined in an acid washed carboy before being siphoned into 600 ml acid washed polycarbonate bottles. For each saturation experiment, 12-14 bottles were used to generate a range of levels of surrogate prey addition. Fluorescent polystyrene microspheres of 2 µm in diameter, treated as for the laboratory experiments, were added to the bottles either as duplicate treatments or in a continuous series of abundance. Maximum surrogate prey abundance ranged from 1.7 × 10 6 beads ml −1 on day 197 to 3.8 × 10 6 beads ml −1 on day 226, in 2019.
Bottles were incubated for 24 h in a flow-through incubator on the dock under a nylon mesh that removed ∼40% of the surface PAR. In an initial test of the flow-through incubator, five 600 ml bottles were filled with seawater from the same carboy and with no further treatment, were incubated for 24 h. The counts of the T 0 abundance of picoeukaryotes from the five bottles showed a coefficient of variation of 1.3%, that increased in the T 24 counts to only 4.1%, indicating consistent growth among the 5 replicate bottles. Following all experiments, the seawater to which beads had been added was filtered through a 0.45 µm capsule filter to recover the beads before discarding the water.
For comparison, a series of dilution experiments was carried out in parallel to the saturation experiments, using a modification of the approach of Landry and Hassett (1982). The same water used for the saturation experiments was used to prepare triplicate bottles of the whole water (100%) and duplicate bottles at four different levels of dilution (10, 25, 50, and 75%) in 600 ml polycarbonate bottles. Diluent was prepared from the 200 µm pre-screened seawater by gentle filtration through a 0.45 µm mini-capsule filter (Pall Life Sciences). The dilution experiments were set up following the saturation incubations to allow time for flow cytometric analysis of each set of subsamples and so, T 0 and T 24 were delayed by ∼ 2 h relative to those of the saturation experiments.

Flow Cytometric Measurements
A ZE5 Cell Analyzer (Bio-Rad, Hercules, CA, United States) was used to measure optical properties of single cells from each sample. Particles were excited with a 488 nm blue excitation laser (100 mW). Data acquisition was triggered on forward scatter (FSC). Signals were recorded from detectors with bandpass filters for right angle light scatter and fluorescence emission in red (692 nm / 80 nm band pass) indicative of chlorophyll a, orange for phycoerythrin (593/52 nm), and green (525/35 nm). To ensure accurate calibration of the flow cytometer, ZE5 QC beads (Bio-Rad, Hercules, CA, United States) were run daily. Data files were analyzed from logarithmic dot plots based on fluorescence and characteristic light scattering properties (Durand and Olson, 1996) using FlowJo 10.6 Software (Becton Dickinson & Company, San Jose, CA, United States) (Supplementary Figure 1).
For culture and natural community experiments, surrogate prey were gated using FSC and green fluorescence. Despite the blocking procedure, some clumping of beads did occur but, in this application, doublet and triplet clumps were treated as single surrogate prey particles. For laboratory experiments, O. danica was gated on side scatter and red fluorescence and the M. pusilla populations were determined using forward scatter and red fluorescence. For field experiments, picophytoplankton were identified based upon cell size, using forward scatter as a proxy, and red fluorescence. Picophytoplankton were classified, as either Synechococcus or picoeukaryotes based on the presence or absence of the accessory pigment phycoerythrin, respectively. Using these gating criteria, the number of cells in each identified population was enumerated and converted to cell abundances based on the processed sample volume.

Single Prey and Predator Interactions
The single prey and predator laboratory culture experiments confirm the theoretical prediction. Apparent growth rate (µ app ) of M. pusilla increased with the addition of beads and tended towards a horizontal asymptote at increasing levels of surrogate prey addition, indicating saturation of ingestion in O. danica (Figure 1). Fits of Equation 2 to the relationship between µ app and the surrogate prey abundance, provided highly significant estimates of µ max (p-value < 0.001), µ net (p < 0.05) and a (p < 0.05) in each of the experiments ( Table 1). Derived values of m decreased from 0.60 to 0.33 d −1 over the three experiments, despite the reduced prey to predator ratios. Although a distinct population of O. danica could be counted initially, as they grazed on M. pusilla and ingested beads, the population became more difficult to define in the flow cytometric signal and it was not possible to obtain final abundances and reliable growth rates for the predator.
Initial prey to predator ratios were approximately 200, 180, and 160 in experiments A, B, and C, respectively (Table 1). At the levels of abundance used, growth and grazing rates were almost balanced in the cultures with no bead addition, resulting in similar, relatively low rates of µ net for M. pusilla in the three experiments (Table 1).
In general, the abundances of freely suspended beads decreased during the incubations by an average of 25, 3, and 25% in A, B, and C, respectively. The decrease possibly results from a combination of ingestion by predators, sinking out of suspension despite the mixing during incubations and sticking to other particles. Proportionally higher losses tended to occur at higher levels of bead addition but this was not always the case.

Saturation Experiments Using Natural Waters
The dynamics of the two picophytoplankton populations showed variable patterns indicative of differences in their rates of growth and the grazing pressure they experienced. During the 2019  Table 1 for experiments (A-C).
summer period, the abundance of the picoeukaryote population increased by approximately 50% from 43,100 cells ml −1 on day 197 to 60,600 cells ml −1 on day 211 and then gradually declined over the following 2 weeks (Figure 2A). In contrast, the abundance of Synechococcus increased more than 10-fold from 2,860 on day 197 to 30,100 on day 211 and remained at a relatively high level to day 226 (Figure 2A). Similar abundances of each group were observed in the summer of 2021 (Supplementary  Tables 1, 2). Note, the mouth of the Damariscotta River Estuary The initial (T 0 ) values of predator and prey abundance are averages and SD of single measurements from all tubes used in each experiment. Rates of growth and grazing mortality and the coefficient a, were derived from the curve fits to an asymptotic regression model (Equation 2), shown in Figure 1. Uncertainty in the estimates of net and gross growth is presented as ± the half width of the 95% confidence interval, and for a ± SE. The uncertainty for grazing mortality is propagated from net and gross growth rate uncertainties (ns, not significant; na, not available).
experiences high levels of tidally driven water exchange and this may also have influenced population dynamics. As illustrated for 2019, increased bead addition resulted in increased µ app , tending toward a horizontal asymptote at bead abundances that exceeded ∼2.0 × 10 6 ml −1 (Figures 2B-H). The growth response to bead addition indicates increasing saturation of the existing grazer population and reduction of grazing mortality in the picophytoplankton populations. Fits of Equation 2 to the relationship between µ app and the surrogate prey abundance, provided highly significant estimates of µ max (p < 0.001), µ net (p < 0.01) and a (p < 0.05), for each of the taxonomic groups in all seven experiments, except for Synechococcus on day 206 ( Table 2). Values of m ranged from 0.26 to 0.43 d −1 for picoeukaryotes and 0.23 to 0.42 d −1 for Synechococcus. Estimates of gross growth in the absence of grazing (µ max ) indicate a highly productive population of picoeukaryotes that grew at rates of over 1 d −1 . These rates of µ max gradually declined over the period to close to 1 doubling d −1 , by day 227 ( Table 2). Synechococcus also showed high rates of µ max initially, that stabilized to between 0.47 and 0.53 d −1 from days 212 to 227. Similar rates of growth and mortality were measured in the four experiments successfully conducted at the same location during the summer of 2021 (Supplementary Table 2).
Paired t-tests showed significant differences in µ net (p = 3.5 × 10 −5 ), m (p = 1.3 × 10 −2 ) and µ max (p = 1.9 × 10 −4 ) occurred between Synechococcus and picoeukaryotes in the combined dataset of saturation experiments. Rates of grazing mortality were the most comparable of the coefficients between the two taxonomic groups, while higher rates of gross growth among picoeukaryotes were largely responsible for the significant differences in net growth rates (Supplementary Figure 2).

Comparison to Dilution Experiments
In general, the rates from each approach show a similar range of values for µ max , µ net and m (Figure 3 and Supplementary Table 1). Successful, parallel experiments allow direct comparison on days 197 and 204 in 2019 and on days 174, 182, 196, and 207 in 2021. It should be noted that the dilution incubations were started approximately 2 h after the saturation incubations, which may explain the slight differences in µ net between experiments on the same days ( Table 2 and Supplementary Tables 1, 2). In the combined dataset of comparable experiments, there was no statistical evidence that the mean difference between experimental approaches was significantly different from zero (paired t-tests, p-value < 0.05). This is apparent in the lack of clear trends between approaches when estimates of µ max and m are compared for picoeukaryotes or Synechococcus, although large variations were apparent on certain days, particularly days 197 and 182 (Figure 4).
Note, the raw data from both the saturation and dilution experiments is included in the Supplementary Material.

DISCUSSION
The combination of laboratory culture experiments and tests in coastal waters illustrate the potential of the saturation approach to generate estimates of the mortality due to microzooplankton grazing (m) and gross growth rates of picophytoplankton in the absence of grazing (u max ). The functional response theory (Holling, 1965) that underpins the saturation approach, appears to hold true in natural communities where diverse phagotrophs, with differing prey preferences and nutritional strategies, are likely to occur.

Picophytoplankton Population Dynamics
The saturation approach confirmed intriguing differences in the dynamics of co-occurring picoeukaryote and Synechococcus populations. The most notable difference being the considerably higher rates of u max in the picoeukaryote population, equivalent to between ∼1 and 2 doublings d −1 compared to generally < 1 doubling d −1 for Synechocccus (Table 2, Figures 3, 4, and  Supplementary Figure 2). Similar, higher rates of gross growth and grazing mortality of picoeukaryotes compared to Synechococcus, have been observed in other temperate and subtropical coastal waters (Worden et al., 2004;Kimmance et al., 2007;Fowler et al., 2020).
At the light levels used for the incubations, approximately 60% surface PAR, µ max of both populations exceeded m, resulting in positive net growth rates (Figure 3). If it is assumed that the same coastal waters were sampled on each date, then the picoeukaryote population showed an increase in abundance between days ∼190 and 211 of 2019, equivalent to a net rate of growth of 0.05 d −1 (Figure 3A). The similar pattern in Synechococcus abundance lagged the picoeukaryote population by approximately 5 d but increased at a faster net rate between days 193 and 211 of 0.15 d −1 (Figure 3A). Even during  Table 2. these periods of maximal in situ net growth, the µ net values derived from the incubations were considerably higher for both populations (Figure 3 and Table 2). The closer balance between gross growth and grazing mortality in the water column as a whole, may be a consequence of reduced overall light levels on gross growth rates caused by mixing of photoautotrophic cells through the water column, while grazing mortality rates may be more constant with depth.

Limitations and Advantages of the Saturation Approach
The approach as presented here, was focused on specific groups of picophytoplankton whose abundance could be precisely quantified by flow cytometry. In this respect, the method may not be a suitable means of obtaining total phytoplankton grazing mortality in the way that the dilution approach has been Rates of growth and grazing mortality for populations of picoeukaryotes and Synechococcus in the coastal waters of the Gulf of Maine between mid-July and mid-August, 2019. Values of predator and prey abundance are averages (±SD) of the initial counts from each of the incubation bottles. Rates were derived from an asymptotic regression model fit to apparent growth (µ app ) vs. level of addition of surrogate prey (Equation 2), shown in Figure 2. Uncertainty in the estimates of net and gross growth is presented as ± the half width of the 95% confidence interval, and for a ± SE. The uncertainty for grazing mortality is propagated from net and gross growth rate uncertainties (ns, not significant; na, not available). applied. How chlorophyll concentrations changed in response to increasing bead addition was not measured but may provide estimates of grazing and growth rates when a single size spectrum dominates the phytoplankton. As illustrated previously and here, the dilution approach can also be used to generate phytoplankton group-specific information that can be directly compared to results from the saturation approach (Figures 3, 4; Landry et al., 1995;Worden and Binder, 2003). The saturation approach relies on several assumptions that need to be considered and more thoroughly tested for applications in specific environments. The primary assumption is that the addition of surrogate prey does not influence µ max of the natural prey. In preliminary laboratory experiments, 2 µm beads were added to six tubes containing only M. pusilla over a range of 0-6 × 10 6 beads ml −1 , to test for an effect of the microspheres on growth. Under the illumination and mixing regime generally used for the culture experiments, growth rates of M. pusilla averaged 0.48 ± 0.03 d −1 across the six tubes, with no apparent trend, indicating minimal influence of the surrogate prey on growth (not shown). Similar tests were not carried out for natural waters, where the separation of prey and predators is more difficult to accomplish.
There is the possibility that in nutrient limited waters, the reduced grazing caused by surrogate prey additions may decrease nutrient recycling rates and as a consequence, potentially reduce estimates of m and µ max . A possible advantage of the saturation approach is that the biomass of the natural population is not altered across the incubation series as part of the initial experimental set-up. Population size does vary, nonetheless, between bottles in the series over the course of the incubations.  Table 2 and Supplementary Tables 1, 2. The dilution approach generated a positive slope for the relationship between apparent growth and level of dilution for Synechococcus on day 182 (Supplementary Table 1).
In the natural water experiments, the maximum variations in population size between incubations exhibiting net growth (zero bead addition) and gross growth (µ max ) within an experiment, were 1.7-fold for picoeukaryotes on day 204 and 1.9-fold for Synechococcus on day 197 (Table 2). Thus, variation in nutrient availability between treatments of the saturation approach is likely to be lower than across the 10-5-fold difference in initial biomass that is required in dilution experiments.
A key aspect of the experimental design is judgment of the quantity of surrogate prey required to illicit a saturation response, while minimizing potential artifacts from too many additional particles. In the coastal waters of the Gulf of Maine, 2 µm bead additions of > 2 × 10 6 ml −1 saturated the grazers of picoeukaryote and Synechococcus populations, sufficiently to generate significant values of m and µ max (Figure 2). It seems reasonable to expect the behavioral response of phagotrophic protists to be consistent in similar water types but this remains to be established. Preliminary studies in oceanic waters suggest grazers of Prochlorococcus and Synechococcus populations may tend toward saturation at bead additions of < 1 × 10 6 ml −1 . In theory, coefficient a of Equation 2 varies in proportion to the rate of increase in µ app relative to C, the surrogate prey concentration. It may be possible to use values of a, or an analogous coefficient, to judge how effectively different surrogate prey types or sizes lead to saturation of predators. As more experiments are carried out, it should be possible to standardize the surrogate prey type and levels of addition for different water types. The asymptotic regression model interpretation of the response of µ app to increasing surrogate prey, provided highly significant values of µ max (p < 0.001) and µ net (p < 0.01) and hence, estimates of m ( Table 2 and  Supplementary Table 2). However, multiple predator types may consume picophytoplankton populations, each having a different "functional response" to increased surrogate prey abundance and the interpretation of the observed experimental results may be considered over-simplistic. Future studies could usefully examine more complex model interpretations of the saturation approach that consider, for example, predator and prey motility and encounter rates, the process and duration of prey capture and selectivity by the predator (Flynn and Mitra, 2016).
The governing role that dissolved compounds play in microbial interactions is becoming ever more apparent (Strom et al., 2007;Stocker et al., 2008). Indeed, a contributing factor to inconsistent results generated by dilution approach is the alteration of seawater chemistry through the filtration step to create the diluent required in the experimental set-up. One specific example, is the release of polyunsaturated aldehydes to the dissolved phase and the resulting increased inhibition of phytoplankton apparent growth with increased levels of dilution (Stoecker et al., 2015). An advantage of the saturation approach is that modification of seawater chemistry may be kept to a minimum. Although the bovine serum albumin treatment of beads and subsequent resuspension in filtered seawater appeared to be sufficient for these experiments, more thorough cleaning protocols for artificial surrogate prey are likely to be required for low nutrient and oceanic waters. Nonetheless, it may be possible to investigate grazer-mediated production of dissolved organic compounds, inorganic nutrients and possibly, trace elements by using the approach to generate a gradient of grazing pressure.

CONCLUSION
The saturation approach provides another experimental option to obtain information on the mortality due to grazing and gross growth rates of specific groups of phytoplankton. It may be particularly suited to investigations into the growth and mortality of picophytoplankton, as demonstrated in the present study. Further investigations are needed into whether the approach is equally effective in determining these rates for larger organisms but it seems reasonable to assume that by increasing the surrogate particle size, larger microzooplankton and zooplankton will also be susceptible to saturation of ingestion rates. It may even be possible to target multiple size classes of prey and predator interactions within a single set of incubations by using a size spectrum of surrogate prey. As seawater chemistry is not altered, the saturation approach may also be an appropriate way to quantify grazing mortality and growth rates of heterotrophic prokaryotes. Although the approach relies on incubations that are long enough to accurately determine changes in prey abundance and therefore growth rates, the experimental set up is less laborious than the dilution approach. In this respect, it may be possible to carry out more experiments in parallel using the saturation approach and thereby obtain greater spatial resolution or treatment-specific information. Simplifying the approach by reducing the number of levels of saturation may also allow more measurements. Several further aspects of the presented approach could be usefully improved including identifying more cost-effective and biodegradable surrogate particles, this would reduce the need to recover the particles post-incubation. It will be particularly useful to establish the numbers of surrogate prey required to generate statistically robust estimates of the rate coefficients in differing pelagic ecosystems, making routine application of the approach more reliable.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.bco-dmo. org/project/735732. The basic data for saturation and dilution experiments are also included in the Supplementary Material.

AUTHOR CONTRIBUTIONS
The approach was conceived by SA with additional input from NP and LL. NP, LL, and KP conducted the culture experiments. LL, MK, KM, WD, and MS conducted the natural water experiments. Manuscript writing was carried out by SA with input from all other authors. All authors contributed to the article and approved the submitted version.

FUNDING
We gratefully acknowledge financial support from the National Science Foundation, United States (NSF project No. OCB 1738061 to SA and NP). WD was supported by the Bigelow Laboratory for Ocean Sciences -Undergraduate Research Experience in the Gulf of Maine and the World Ocean Program (NSF REU 1460861). Matthew Wade was supported by the Linde Packman Lab for Biosciences Innovation and the Champlin Science Scholars Fund at Colby College. Instrumentation was supported by NSF FSML project 1821065.