ABSTRACT
Prochlorococcus cyanobacteria grow in diurnal rhythms driven by diel cycles. Their ecology depends on light, nutrients, and top-down mortality processes, including lysis by viruses. Cyanophage, viruses that infect cyanobacteria, are also impacted by light. For example, the extracellular viability and intracellular infection kinetics of some cyanophage vary between light and dark conditions. Nonetheless, it remains unclear whether light-dependent viral life history traits scale up to influence population-level dynamics. Here, we examined the impact of diel forcing on both cellular- and population-scale dynamics in multiple Prochlorococcus-phage systems. To do so, we developed a light-driven population model, including both cellular growth and viral infection dynamics. We then tested the model against measurements of experimental infection dynamics with diel forcing to examine the extent to which population level changes in both viral and host abundances could be explained by light-dependent life history traits. Model-data integration reveals that light-dependent adsorption can improve fits to population dynamics for some virus-host pairs. However, light-dependent variation alone does not fully explain realized host and virus population dynamics. Instead, we show evidence consistent with lysis saturation at relatively high virus-to-cell ratios. Altogether, our study represents a quantitative approach to integrate mechanistic models to reconcile Prochlorococcus-virus dynamics spanning cellular-to-population scales.
IMPORTANCE The cyanobacterium Prochlorococcus is an essential member of global ocean ecosystems. Light rhythms drive Prochlorococcus photosynthesis, ecology, and interactions with potentially lethal viruses. At present, the impact of light on Prochlorococcus-virus interactions is not well understood. Here, we analyzed Prochlorococcus and virus population dynamics with a light-driven population model and compared our results with experimental data. Our approach revealed that light profoundly drives both cellular- and population-level dynamics for some host-virus systems. However, we also found that additional mechanisms, including lysis saturation, are required to explain observed host-virus dynamics at the population scale. This study provides the basis for future work to understand the intertwined fates of Prochlorococcus and associated viruses in the surface ocean.
INTRODUCTION
The unicellular cyanobacterium Prochlorococcus dominates the phytoplankton community and is a major contributor to primary production in tropical and subtropical oligotrophic oceans (1). The ecology of Prochlorococcus is a function of physicochemical properties of the marine environment (2–4), bottom-up forces (i.e., nutrient driven), and top-down (i.e., mortality-driven) effects (5–11). Among top-down factors, cyanophage (i.e., viruses that infect cyanobacteria) are highly abundant and can drive up to 30% of cyanobacterial mortality in marine environments (12–17). Light, temperature, and nutrients influence Prochlorococcus growth (2, 3, 18), as well as its interactions with cyanophage (19).
Prochlorococcus cyanobacteria are distributed across temperature and light gradients in the ocean environment (3, 20–22). They are specialized into high-light (HL)- and low-light (LL)-adapted ecotypes (23–25). LL ecotypes have a high fluorescence and exhibit photoinhibited growth at medium light intensity. They grow faster at low irradiance with a high concentration of divinyl chlorophyll a and b and have several pcb genes encoding constitutive photosystems I and II (23–25). In contrast, HL ecotypes grow faster at medium light intensities, have a low concentration of divinyl chlorophyll a and b, and have only constitutive photosystem II light-harvesting complexes (23–25). Prochlorococcus cells do not have a circadian rhythm; rather, they have a diurnal rhythm that can be synchronized under light-dark cycles (21, 26, 27). This diurnal rhythm is divided predominantly into photosynthesis during the light phase and cell division associated with energy consumption during the dark phase (22, 28).
Cyanophage are also impacted by light. UV radiation can deactivate and degrade virus particles (19), as well as degrade and modify viral proteins and genomes (29, 30). Light can also affect viral interactions with host cells. Previous studies suggested a dependence of viral production on the host cell cycle in different phytoplankton lineages (31–34). In contrast, the intracellular production of some viruses may be decoupled from host cell cycle and light levels (35, 36). A recent paper on the diel infection pattern of cyanophage infecting Prochlorococcus (37) suggests that adsorption, as well as the transcription rhythm of cyanophage, may be related to the light-dark cycle and not to the host cell cycle. Analyzing light impacts on cyanophage-cyanobacterium dynamics requires some elaboration of the cyanophage life cycle.
The lytic cyanophage cycle can be summarized into three phases: the adsorption phase, where virions attach to their host and inject their genetic material into the host cell, the cyanophage replication phase, and the lytic phase, where new virus particles are released by lysing their host (Fig. 1a). Light can affect each of these phases and associated viral life history traits (LHTs). A study on Synechococcus infection showed a significant decrease in adsorption under the dark condition for some phages (36, 38, 39), whereas other cyanophage adsorb during light or dark conditions (36). Similarly, light conditions also modify the cyanophage replication phase. An increase in viral production in light and a reduction in viral production in dark have been reported for Synechococcus (39–41) and Prochlorococcus (42–44).
Cyanophage infection in the light or the dark. (a) Viral life history trait definitions: viral adsorption (encounter and adsorption on a noninfected host, in ml h−1), latent period (time between adsorption and lysis of the host cell, in hours), and burst size (new phages produced per one lysed host cell). (b) Host-cyanophage pairs used in the study. (c) Infection under light or dark (data from reference 37; see Materials and Methods). Cyanophage P-HM2 and P-SSP7 were used to infect their host cells under continuous light or in the dark. For all the host-phage pairs, the phage/host ratio is 0.1. Extracellular phage concentrations were measured as phage DNA by quantitative PCR and normalized to the value at time zero.
A salient example of light-driven changes to viral LHTs is found in Liu et al. (37), who investigated infection dynamics for cyanophage infecting Prochlorococcus under light-dark cycles (Fig. 1b and c). The results suggest that cyanophage strains respond differently to light or dark conditions (Fig. 1c). Infection under light was always efficient for all strains. However, P-SSP7 could infect and produce viruses in the dark, while P-HM2 could not adsorb to hosts or produce viruses in the dark. These observations under fixed light or dark conditions form the central motivation for our study. Here we ask: do differences in the response of viruses to light conditions at the cellular level explain population-level dynamics of both Prochlorococcus and cyanophage given diurnal rhythms of light-dark cycles?
To answer this question, we couple mathematical models, high-resolution (i.e., subdaily) measurements, and model-data integration to explore the interactions between Prochlorococcus strain MED4 (an HL ecotype) and cyanophage P-HM2 and P-SSP7 (Fig. 1b). The bulk of data used in this study have been presented elsewhere (see data attribution statements) (37, 45). The models extend the framework of nonlinear population dynamics of lytic viruses and their hosts (46) to an explicitly light-driven context (see the related work of reference 47 on coccolithoviruses and their Emiliana huxleyi hosts given daily measurements). As we show, although diel-driven viral life history traits help explain population dynamics, they are not necessarily sufficient. Instead, our study identifies additional mechanisms involving saturating lysis that help reconcile population-level dynamics of cyanophage and cyanobacteria.
RESULTS
Light-driven Prochlorococcus growth.We first estimated the growth of Prochlorococcus strain MED4 in culture under light-dark cycles and nutrient-nonlimited conditions during the exponential-growth phase. We used an ordinary differential equation (ODE) model to describe the dynamics of the Prochlorococcus population (P) (cells/ml), as follows:
The model in equation 1 was fit to population abundance measurements of Prochlorococcus strain MED4 under a light-dark cycle (45) using a Markov chain Monte Carlo (MCMC) approach (see Materials and Methods). The best-fit light-driven host growth model recapitulates the experimental data (Fig. 2a) with a good convergence of the MCMC parameter chains (Fig. S1a and Table S1 in the supplemental material). MED4 has a low growth-irradiance curve slope (α = 0.0011 h−1), a high optimal growth irradiance (Lopt = 44.78 μmol s−1 m−2), and a maximum growth rate of 0.0035 h−1 (Fig. 2b and Fig. S1a). These parameters are consistent with prior estimates of HL growth-irradiance characteristics for strain MED4 (Fig. 2c) (50).
Modeling Prochlorococcus MED4 strain as function of light without viruses during the exponential phase. (a) Fit of the host dynamic (equation 1). Solid lines represent the median of 5,000 model simulations, and shaded areas are the 95% quantiles. Black dots are data (from reference 45) for two replicates, and gray shaded area represents the dark condition. (b) Model growth parameter distributions of the host model (equations 1 and 4). Parameter distribution estimated using an MCMC algorithm: photosynthesis-irradiance (PI)-curve slope of the linear phase α, optimal growth light Lopt, maximal growth rate μmax, minimum amount of light necessary to divide KL, and natural mortality ω. (c) Model growth functions that drive the host dynamic: growth is expressed as the net growth rate (μopt − ω) as a function of irradiance (equation 4; left) and as a function of time (equation 2; right).
FIG S1
Copyright © 2020 Demory et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
TABLE S1
Copyright © 2020 Demory et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
Modeling Prochlorococcus-phage dynamics under light-dark cycles.To investigate Prochlorococcus-cyanophage dynamics under light-dark cycles, we developed a nonlinear ODE population model describing the infection of Prochlorococcus by cyanophage (Fig. 3a), extending existing frameworks for modeling obligately lytic phage-host interactions (46). We used a two-stage infection model to account for the finite latent period of the infection. The host population is categorized as susceptible (S), exposed (E), and infected (I), such that the total host population is N = S + E + I. The density of free cyanophage is denoted by V. The model is described by the following system:
Description of the model. (a) Schematic representation of the model. The host population is divided into 3 classes: susceptible (S), exposed (E), and infected (I) types. The virus particle density is denoted by V. Black arrows are biological processes. (b) Definitions of the hypotheses. Each hypothesis describes a possible relation between light and the infection parameters. When parameter ϕ, β, or λ is light dependent, it is a piecewise function, i.e., it takes one value in light and one value in dark.
FIG S2
Copyright © 2020 Demory et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
TABLE S2
Copyright © 2020 Demory et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
We fit each of the nested, light-driven virus-host population models, using MCMC, to experimental measurements of Prochlorococcus strain MED4 infected by either cyanophage P-HM2 or P-SSP7 over a 4-day period (Fig. 4). Parameter ranges in the MCMC procedure were constrained by prior estimates (Tables S3 and S4) (51). We found the best-fit model while accounting for the inclusion of additional model complexity to be H2ϕλ for P-HM2 and H0 for P-SSP7 (Fig. S3). This suggests that P-HM2 (but not P-SSP7) has light-dependent life history traits that help provide explanatory power to the virus-host population dynamics. The best-fitting model suggests that both adsorption and production for P-HM2 are significantly reduced in the dark compared to its adsorption and production under light conditions. In both cases, viral abundances rapidly increase and then plateau. However, in both cases, the best-fit model significantly overestimates the degree of virally induced mortality in the host population, e.g., models predict that the final time point estimates of cell density are 2.5 and 6.1 times lower than the values measured for the P-HM2 and P-SSP7 cases, respectively. This result suggests that other features underlying interactions between cyanophage and Prochlorococcus have to be accounted for when scaling up to the population-level dynamics.
Light-driven models fit to host and virus population abundance data. Model fits under H0 and hypotheses H2ϕλ and H1ϕ for an inoculation time of 14.5 h after the beginning of the experiment. Phage P-HM2 infecting strain MED4 (left) and P-SSP7 infecting MED4 (right). Solid lines represent the median values of 5,000 model simulations, with shaded areas the 95% quantiles. Data are represented by the black dots for two replicates. Vertical shaded gray lines represent dark conditions.
FIG S3
Copyright © 2020 Demory et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
TABLE S3
Copyright © 2020 Demory et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
TABLE S4
Copyright © 2020 Demory et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
Beyond light: incorporating lysis inhibition to explain virus and host population dynamics.The observation that host populations remain persistently above model expectations when viral abundances are high suggests a potential slowdown mechanism in virally induced lysis. To account for this, we modified the initial model to account for an additional state transition, i.e., from I to E, as follows:
We then fit the lysis inhibition model to an expanded set of experimental measurements of MED4 infected by either cyanophage P-HM2 or P-SSP7. The measurements comprise additional time courses for each cyanophage. The time courses have different virus inoculation times of 14.5, 18, 24.5, 30, and 36 h (Fig. 5). The light-dependent hypotheses used in fitting are denoted
Viral dynamics under light-dark cycle for different viral inoculation times. Model fits under H0 and best hypotheses
FIG S4
Copyright © 2020 Demory et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
We evaluated the quality of fits by assessing the predicted estimates of life history traits for the P-HM2 and P-SSP7 cases. Disparities in parameters under light or dark conditions obtained with our MCMC approach are consistent with earlier measurements of viral infections of MED4 given fixed light or dark conditions over a 10-h period (37). Specifically, model fits reveal that P-HM2 has a significantly lower adsorption rate in the dark than in the light (Fig. 6 and Table S3). Indeed, dark adsorption is at the lower limits of the parameter constraint range of the MCMC procedure, suggesting that P-HM2 may have effectively zero adsorption in the dark. In contrast, model estimates cannot reject the hypothesis that adsorption was effectively constant for P-SSP7. The convergence of MCMC chains further supports the robustness of the model-based inferences (Table S4 and Fig. S1b and S1c). Notably, other candidate models with intracellular mechanisms that delay lysis can also reproduce similar population-level features (Text S1 and Fig. S5).
Model infection parameter distributions: P-HM2 (left) and P-SSP7 (right) infecting strain MED4. Distributions are calculated with 5,000 parameter sets. Distributions are colored according to the model hypothesis: H0 (gray),
FIG S5
Copyright © 2020 Demory et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
TEXT S1
Copyright © 2020 Demory et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
DISCUSSION
We investigated the impact of light and dark conditions on the infection of Prochlorococcus by cyanophage using a combination of experiments, nonlinear population models, and model-data integration. We found that light-dependent differences in viral adsorption to hosts help explain population-level changes in both virus and host abundances given growth under diurnal conditions. These light-dependent differences are strain specific. Estimated adsorption rates vary markedly during the light versus the dark for P-HM2 but not for P-SSP7. This suggests that viruses, in addition to hosts, may have light-dependent differences in their life history traits at the cellular scale that impact dynamics at population scales.
In our model-fitting procedure, we evaluated the possibility that light could affect adsorption, latent period, and burst size. We only found evidence for light-dependent variation in adsorption rate for the phage P-HM2. In contrast, P-SSP7 dynamics were explained by light-driven variation in the host growth rate only. Our results for both viruses corroborate the observations in reference 37, supporting evidence of mechanisms of light-driven Prochlorococcus infection dynamics for P-HM2 but not for P-SSP7. However, model estimations for P-SSP7 burst size were higher than the values reported in the literature, suggesting that other factors may be involved in the dynamics. The imputed failure of P-HM2 to adsorb to MED4 in the dark indicates that adsorption could be directly modulated by light (36). Light-dependent variation in adsorption has also been reported in cyanophage infecting Synechococcus (36, 39) and in coccolithoviruses (47).
There are multiple reasons why P-HM2 may have evolved light-dependent viral LHTs. First, exposure to UV is a critical factor degrading viral particles outside the host cell (19). During the night, there is both less UV and (potentially) elevated predation rates of cyanobacteria by eukaryotic grazers (52, 53). Therefore, remaining outside the host cell during the night could effectively amount to a survival strategy by avoiding predation by grazers on the viral host (20). The evolution of light-dependent LHTs may also be host dependent, e.g., as a response to variation in host availability.
Despite our focus on light-driven traits, our approach revealed other mechanisms driving variation in host-virus population dynamics. The failure of a light-driven virus-host population model to recapitulate the persistence of host cells suggests that other feedback mechanisms may limit host mortality, even when virus densities are relatively high. Using a variant of the original model, we found evidence consistent with lysis inhibition at high viral densities (54). Mechanisms consistent with lysis inhibition include decreases in viral infectivity, an increase in the production of defective viral particles, or slowdowns in host physiology. Such slowdowns reflect the potential reciprocal influence of processes at cellular and population scales. The relevance of such slowdowns will vary with the environment. For example, in marine surface environments, cyanophage densities do not typically exceed 106 ml−1, and so it remains unclear whether the candidate feedback mechanism is an adaptive response to the high density of infected hosts or arises incidentally given ecological conditions outside typically encountered ranges. Further work is necessary to disentangle process from pattern.
In closing, we found that light-dependent viral life history traits can substantively change the dynamics of Prochlorococcus and cyanophage. This finding reinforces and extends the consequences of prior results showing that viral traits differ between light and dark, albeit under fixed conditions. In the marine environment, adaptation to light has been shown to drive differences in physiology among Prochlorococcus cyanobacteria, as well as evolutionary adaptation between light-associated ecotypes. Our study suggests that exploring variation in virus-associated light-dependent life history traits may also reveal ways in which viruses partition their environment, both in terms of host specificity and via differential infection of hosts over light-dark cycles.
MATERIALS AND METHODS
Experimental design and data attribution.The experimental data analyzed here comprise data from two published sources, Liu et al. (45) and Liu et al. (37), and new data collected to link infection-level dynamics with population-level dynamics. Specifically, the host growth data in Fig. 2 were previously reported in Liu et al. (45). The infection data in Fig. 1C and D, as well as the host and phage abundances before 60 h in Fig. 4 and 5, are reported in Liu et al. (37). The host and phage abundance data after 60 h in Fig. 4 and 5, as well as the viral decay data reported in Fig. S7 in the supplemental material, are newly reported here. Details of the experimental procedures are described in the following sections, with full quotations used when methods are equivalent to those reported in reference 37. We include full method descriptions for completeness.
Culture conditions.Culture conditions were as described by R. Liu et al. (37):
Axenic Prochlorococcus strains were grown in Port Shelter (Hong Kong) seawater-based Pro99 medium (55). Batch cultures were incubated at 23°C in continuous light (25 μmol quanta m−2 s−1) or a 14h light:10h dark cycle (35 μmol quanta m−2 s−1 in the light period). Cultures were acclimated in the same condition for at least three months before they were used for the experiments.
Infection of synchronized Prochlorococcus cells under light-dark cycles.Infection of Prochlorococcus cells under light-dark cycles was as described by R. Liu et al. (37):
Prochlorococcus cells were acclimated under light-dark cycles for at least three months and were synchronized, as determined by flow cytometry. Mid-log cells were infected at different times of a light-dark cycle at a phage/host ratio of 0.02. Times of infection were 0.5, 6, 12 hours. Each experiment was replicated at least two times.
Cyanophage decay rates under light or dark.To measure the decay rates, fresh lysates of cyanophages P-HM2 and P-SSP7 were prepared by adding 300 μl virus stocks into 30 ml mid-log Prochlorococcus MED4 culture. After the infected culture became clear, cell debris was removed using a 0.2-μm polycarbonate filter and the supernatant containing phage particles was stored at 4°C in the dark. During the viral decay experiment, aliquots of viral lysates were incubated at 23°C at a light intensity of 27 μmol photons m−2 s−1, and aliquots were incubated at the same temperature in the dark (56). Samples were taken from each tube every 2 days over 10 days. To measure the loss of viral infectivity, the numbers of PFU were measured (57). Briefly, 500 μl serially diluted viral lysate was added to 2 ml Prochlorococcus MED4 (grown to mid-log phase in Pro99) in glass tubes and incubated at room temperature for 15 min to allow phage adsorption. Incubated cultures were then combined with UltraPure low-melting-point agarose (Invitrogen) at a final concentration of 0.5%. The EZ55 Alteromonas helper bacterium (58) was added to every plate. Plaques began to appear 7 days later on plates that were incubated at 23°C at a light intensity of 19 μmol photons m−2 s−1. Each sample was measured with three technical replicates.
Flow cytometry and cell cycle analysis.Flow cytometry and cell cycle analysis were performed as described by R. Liu et al. (37):
Prochlorococcus cells were preserved by mixing 100 μL culture with 2 μL 50% glutaraldehyde to a final concentration of 1% and were stored at −80°C. Cells were enumerated by a flow cytometer (BD FACSCalibur) with the CellQuestPro software. We followed a published protocol to determine the percentage of cells in each cell cycle stage (22). Briefly, Prochlorococcus cells were stained with the DNA stain SYBR Green (Invitrogen) and flow cytometry data were analyzed with the ModfitLT software.
Quantification of cyanophage.Cyanophage were quantified as described by R. Liu et al. (37):
Total phage particles were collected on a 0.02 μm Whatman Anodisc filter, stained with SYBR gold (Molecular Probes), and counted under an epifluorescence microscope (59, 60). At least five discrete fields on a filter were photographed using the SPOT Advanced Imaging software and fluorescent dots representing phage particles were counted manually.
During infection, extracellular phage DNA was quantified using a quantitative polymerase chain reaction (qPCR) method (61). Briefly, infected Prochlorococcus cultures were filtered through 0.2 μm polycarbonate filters in a 96-well filter plate (Pall). Filtrates containing extracellular phage particles were diluted 100-fold in dH2O and were then used as templates for qPCR reactions in a 384-well plate. A qPCR reaction contained 4.6 μL template, 0.2 μL forward primer (10 μM), 0.2 μL reverse primer (10 μM), and 5 μL iTaq Universal SYBR Green Supermix. The LightCycler 480 Real-Time PCR System (Roche Diagnostics) was used for thermal cycling, which consisted of an initial activation step of 5 min at 95°C, 45 amplification cycles of 10 s at 95°C and 60 s at 60°C, and a melting curve analysis at the end. The number of cyanophage in each well was quantified using a standard curve generated from phage particles that were enumerated by epifluorescence microscopy.
Measurement of phage DNA copies by qPCR provides an ∼1:1 relationship with phage particle counts (37, 51).
Model simulations.Model analyses were performed with MATLAB version 9.2.0 (The MathWorks, Inc., 2017; Natick, MA). Infection dynamics were simulating using MATLAB ODE solver ode45 (62) (The MathWorks, Inc.), which uses a higher-order Runge-Kutta method (63).
Estimation of the best parameter sets. (i) General procedure.We aimed to estimate the parameter set θ which best described the measurements of host and virus abundances, given a particular model. The procedure consisted of two steps. First, we estimated host growth parameters θhost for the model without viruses (equation 1). Next, we used θhost in the model with viruses (equation 5 and 6) and estimated the infection parameters θinfection for each hypothesis. We estimated parameters by minimizing an objective function that describes the error between model fit and experimental data.
(ii) The objective function.The objective function calculated error between the model fits and the measurements as follows:
(iii) Algorithms.We first sampled the parameter space with Latin hypercube sampling (LHS) (64) for 20,000 parameter sets for each hypothesis and model. We calculated an initial error for each parameter set by running the dynamical model and calculating the objective function (equation 7). Next, we implemented a Markov chain Monte Carlo (MCMC) procedure for the 10 best initial parameter sets. We used two “burn-in” periods (running periods that allow the convergence of the chains). The 10 parameter distributions from the resulting chains each consisted of 5,000 parameter sets. The distributions tended to overlap; when this happened, we chose the best distribution from the overlapping set, as quantified by the median of the error (equation 7). We used the MCMC toolbox for MATLAB, implementing the DRAM algorithm (65).
(iv) Estimation of the host growth parameters θhost.For the host growth parameter sets θhost = (α, Lopt, μmax, KL, ω, K), we used the procedure described above to estimate parameters for Prochlorococcus strain MED4. Parameters of the growth-irradiance curves (α, Lopt, and μmax; equation 4) were constrained by the data from Moore and Chisholm (50), whereas KL and ω were not constrained. The carrying capacity K was fixed and considered a constant (K = 3.109 cell ml−1 [according to reference 66] for nonaxenic cultures).
(v) Estimation of the infection parameters θinfection.To estimate the parameter set θinfection = (ϕ, λ, β), we fixed the host growth parameters estimated previously and estimated the infection parameters relative to the hypotheses H0 to H7. Depending on the hypothesis, the estimated parameter could be constant during the experiments (no relation with light or dark condition) or a piecewise function (with differing light and dark values). The estimated parameters were the adsorption rate ϕ, the latent period λ, and the burst size β. Viral decay rates were estimated experimentally as the slope of log(viral concentration) regression under light or dark conditions and fixed (Fig. S2).
Quantifying the best model hypothesis.To quantify the best model under hypotheses H0 to H7, we computed an Akaike information criterion (AIC) and a Bayesian information criterion (BIC) (67) for virus and host (when data were available) according to the following equations (equations 8 and 9):
Estimation of adsorption with experimental data.We estimated the experimental adsorption from the viral data (V) of Fig. 1. We assumed that the host growth was negligible in the first hours of the experiments (t < 6 h) and estimated the adsorption as follows:
Data availability statement.All data are available for use and reuse. The full data set and code are available at https://doi.org/10.5281/zenodo.3308790. As noted in “Experimental design and data attribution” above, the analysis here includes data from both published and unpublished sources. Data in Fig. 1c and d are from reference 37. Data in Fig. 2a are from reference 45. Data in Fig. 4 and 5 until 60 h are reused from reference 37, and new data are used for time points after 60 h. Data for Fig. S2 are original.
ACKNOWLEDGMENTS
This study was supported by grants to Qinglu Zeng from the National Natural Science Foundation of China (project number 91851112), the Research Grants Council of the Hong Kong Special Administrative Region, China (project number 16144416), and the Hong Kong Branch of Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou—project number SMSEGL20SC01) and by a grant to Joshua S. Weitz from the Simons Foundation (SCOPE award identification number 329108).
We thank Akram Salam for feedback on the manuscript.
D.D., Q.Z., and J.S.W. designed the study. R.L., Y.C., and F.Z. performed experiments, designed by Q.Z. D.D. carried out the modeling analysis with contributions from J.S.W. D.D. and Y.C. performed candidate model analysis. A.R.C. and Y.C. contributed to code development. D.D. and J.S.W. wrote the manuscript with contributions from A.C. and Q.Z.
The authors declare that they have no conflict of interest.
FOOTNOTES
- Received October 18, 2019.
- Accepted February 26, 2020.
- Copyright © 2020 Demory et al.
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.