Next Article in Journal
The Envelope Gene of Hepatitis B Virus Is Implicated in Both Differential Virion Secretion and Genome Replication Capacities between Genotype B and Genotype C Isolates
Next Article in Special Issue
Seasonal Dynamics of Haptophytes and dsDNA Algal Viruses Suggest Complex Virus-Host Relationship
Previous Article in Journal / Special Issue
Isolation and Characterization of a Shewanella Phage–Host System from the Gut of the Tunicate, Ciona intestinalis
 
 
Correction published on 24 April 2017, see Viruses 2017, 9(4), 89.
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Emerging Interaction Patterns in the Emiliania huxleyi-EhV System

1
Department of Biology, University of Bergen, Bergen 5006, Norway
2
NRL for fish, Shellfish and Crustacean Diseases, Central Veterinary Institute of Wageningen UR, Lelystad 8221 RA, The Nederlands
3
Uni Research Environment, Nygårdsgaten 112, Bergen 5008, Norway
*
Author to whom correspondence should be addressed.
Viruses 2017, 9(3), 61; https://doi.org/10.3390/v9030061
Submission received: 30 January 2017 / Revised: 15 March 2017 / Accepted: 16 March 2017 / Published: 22 March 2017
(This article belongs to the Special Issue Marine Viruses 2016)

Abstract

:
Viruses are thought to be fundamental in driving microbial diversity in the oceanic planktonic realm. That role and associated emerging infection patterns remain particularly elusive for eukaryotic phytoplankton and their viruses. Here we used a vast number of strains from the model system Emiliania huxleyi/Emiliania huxleyi Virus to quantify parameters such as growth rate (µ), resistance (R), and viral production (Vp) capacities. Algal and viral abundances were monitored by flow cytometry during 72-h incubation experiments. The results pointed out higher viral production capacity in generalist EhV strains, and the virus-host infection network showed a strong co-evolution pattern between E. huxleyi and EhV populations. The existence of a trade-off between resistance and growth capacities was not confirmed.

1. Introduction

Since the discovery of high viral concentrations in the marine environment, normally ranging between 107 and 1011 virions/L [1], hypotheses regarding the potential impact those viruses could have on their microbial host populations, have been put forward. Viral-induced microbial lysis in Earth’s oceans could amount to an impressive 1023 new infections per second, releasing up to 109 tons of cellular carbon every day [2,3]. Consequently, viral lysis contributes greatly to marine biogeochemical cycling of nutrients as well as reducing the transport of organic matter to upper trophic levels in a process known as viral shunt [4,5,6]. Through horizontal gene transfer and the lysis of their hosts, marine viruses contribute to structuring the diversity and composition of microbial communities [7,8,9,10,11].
Viral activity has been suggested as a plausible mechanism contributing to explain Hutchinson’s paradox, which questions the existence of highly diverse planktonic communities in nutrient limited environments [4,12,13]. Viral strain or species-specific lysis may potentially explain the coexistence of cells with different growth and resistance capacities [14,15]. This scenario is contemplated in the Killing-the-Winner (KtW) hypothesis, notably with the concept that resistance has an inherent cost. This trade-off, also known as cost of resistance (COR), ultimately regulate the co-existence of competition specialists (with higher growth rates) and defence specialists (with higher immune capacity against viral infection), respectively [16].
COR can be detected by analysing the virus-host infection network patterns (VHINs) that emerge after cross-infectivity experiments [17,18,19]. The most frequently tested VHIN patterns are nestedness and modularity [17,20]. Nested patterns are characterized by specialist viruses tending to infect the most susceptible hosts, while the viruses with broader host-range infect hosts that are more resistant [21]. On the other hand, in modular patterns the interactions tend to occur within different groups of viruses and hosts, but not between groups [17,22].
The role of viruses as an important driver of microbial diversity has become clear in prokaryotic-virus systems [23,24,25,26,27] such as the Pseudoalteromonas [28] and the Pseudomonas aeruginosa host-virus systems, in which resistant cells emerging after infection were less competitive than the sensitive ones [24]. In other prokaryote-virus systems that role remains elusive [29,30,31,32,33,34]. The very few examples of trade-off between resistance and growth rate in eukaryotic hosts include studies on the prasinophyte Ostreococcus tauri [35] and the trebouxiophyte Chlorella variabilis [36].
Here we aim at getting insight on the main emerging patterns that result from eukaryotic host-virus interactions in the planktonic realm by focusing on Emiliania huxleyi, the most abundant and widely distributed calcifying haptophyte in our oceans [37], and its lytic viruses. Mostly known for its impressive blooms [38,39] this microalga is an important player in global geochemical cycles [40,41]. This photosynthetic unicellular eukaryote is infected by Emiliania huxleyi viruses (EhV), lytic giant viruses belonging to the genus Coccolithovirus, within the Phycodnaviridae family. These viruses are ubiquitous in the marine environment [42] and abundant, reaching 107/mL in natural seawater during bloom conditions and from 108 to 109/mL in laboratory cultures [43]. Genomic and metagenomic EhV characterizations show both a global consistency of this viral genome on a planetary scale as well as the maintenance of specific localized genetic traits. For example, despite the high levels of sequence similarity (>95%) between EhV isolates from a Norwegian fjord and the English Channel, these viral populations also contain distinctive genetic traits [44,45,46,47,48,49,50]. It is surprising that these genetic traits have been maintained through decades although no geographical isolation and speciation have occurred to date [45], allowing these viral communities to infect hosts from distant geographic places [44,51].
Taking advantage of the large number of E. huxleyi cell and EhV lines available for this host-virus system, from diverse geographical origins that include the major oceanic regions, an extensive array of cross-infectivity experiments was conducted in order to investigate parameters such as growth rate (µ), resistance (R), and viral production (Vp). We then confronted possible existence of correlations between those parameters with the theoretical hypotheses (Table 1) that delimit our conception of virus-microbe interactions in the oceans and the way we model those interactions.

2. Materials and Methods

2.1. Emiliania Huxleyi and EhV Strains

Algal strains were obtained from the Roscoff Culture Collection, France; and from the University of Bergen, Norway. A total of 49 E. huxleyi strains (Table S1) were maintained in 30 mL polystyrene flasks with IMR ½ medium [69] at 16 °C and a 14:10 h light:dark illumination cycle at 155 µmol photon m−2/s irradiance.
A total number of 13 viral strains were obtained from the Plymouth Marine Laboratory, UK; and from the University of Bergen, Norway (Table S2). For all viral isolates, viral stocks were produced by infection of exponentially growing E. huxleyi RCC1257 strain. Viral lysates were centrifuged at 12,000 × g for 20 min and the supernatant was filtered through a 0.45 µm syringe filter (Whatman plc, GE Healthcare Life Sciences, Kent, UK) to remove cellular debris. Viral stocks were kept at 4 °C in the dark and were renewed so often as to never be more than 2 weeks old before inoculation in order to preserve the agent’s viability. Plaque assays were not conducted as haptophytes in general do not grow on agar plates and have only been achieved for a few E. huxleyi strains [70,71].

2.2. Cross-Infectivity Experiments

Cross-infectivity experiments were performed between all the E. huxleyi and EhV strains (Table S3). Prior to each experiment, E. huxleyi cultures were maintained in exponential growth phase with cell concentrations ranging from 105 to 106 cells/mL. The experiments were performed in 24 well culture plates under the same temperature and light conditions as the general culturing conditions described above. Triplicates of 2 mL of each algal culture (1 × 105 cells/mL) were inoculated with each of the 13 viral strains at a concentration of 1 × 106 viral particles/mL, resulting in a virus to host ratio (VHR) of 10. Three replicates of uninfected culture were also used as a control for each E. huxleyi strain. An incubation time of 72 h was chosen because this is consistent with the time scales reported for E. huxleyi/EhV selection dynamics observed in the natural environment [72,73,74]. Moreover, preliminary growth tests [75] performed on several E. huxleyi strains, did not indicate that prolonged incubation period would contribute essential knowledge on the growth capacity of each strain.

2.3. Enumeration of Algae and Viruses

At times 0 h and 72 h, 500 µL was subsampled from each well to determine algae and virus concentrations using a FACSCalibur BC flow cytometer (Becton–Dickinson, Biosciences, Franklin Lakes, NJ, USA) [76,77,78] provided with an air-cooled laser procuring 15 mW at 488 nm. Viral samples were fixed with 20 µL of glutaraldehyde (25%) for 30 min at 4 °C, and frozen at −80 °C until further use. For flow cytometry analysis, samples were thawed, diluted 500-fold in TE buffer (10:1 mM Tris:EDTA, pH 8, filtered through 0.2 µm), and stained with SYBR Green I 100× diluted (Invitrogen, 1600 Faraday Avenue, PO Box 6482, Carlsbad CA, 92008 United States) for 10 min at 80 °C before analysis. Algal enumeration was conducted on fresh samples, and cell populations were discriminated using chlorophyll auto-fluorescence (670 LP) and SSC signals. Virus populations were determined and enumerated on basis of their green fluorescence (530/30) and SSC signals.

2.4. Growth Rate, Resistance, Viral Production

Growth rates (µ) were calculated for each E. huxleyi strain using the control non-inoculated incubations according to the following formula [79]:
µ = Ln (N2/N1)/t
where N1 and N2 were the cell concentrations at the beginning and end of the experiment, respectively, and t was the incubation time in days.
The level of resistance of each E. huxleyi strain to viral infection was measured in two manners. The first manner (R1) was based on the difference of cells that were not lysed after incubation with viruses, compared to the non-inoculated controls. For each E. huxleyi strain a resistance value was hence calculated against each of the 13 EhV strains and the 13 resistance values were then averaged to obtain an overall resistance capacity for each alga strain (R1). Resistance was also estimated as the number of EhV strains that successfully produced progeny on that host (R2).
A value of viral production (Vp), corresponding to the capacity of each viral strain to produce new progeny on a certain host, was calculated for each virus — host pair as the difference between final and initial viral concentrations. These values were averaged to obtain a global infectivity capacity for each viral strain, per algal strain. The maximum amount of viruses that each EhV strain, per algal strain, produced was registered as “Maximum viral production”.
Potential correlations between the different parameters (growth rate, resistance, and viral production) were investigated with regression slopes and statistical probability analyses, using either Anova (F) or Pearson analysis.
A potential impact of domestication on these parameters was also investigated. An analysis was performed on two groups of E. huxleyi strains, which were isolated in different periods of time. The periods before and after 2009, respectively, were chosen for an apparent increase in Vp was preliminary observed in strains as old or younger than 2009 (Figure S1).

2.5. Host-Virus Network Analysis

In order to test the structure of the infection network, we used the BiMat package for Matlab [21]. This network-based analysis was applied on a binary matrix where 0 referred to no lysis and 1 to lysis. The NODF algorithm was used to measure nestedness and is based on overlap and decreasing fill [80]. It returns a score between 0 and 1, where 1 corresponds to a perfectly nested structure. Modularity (Qb) was calculated using the Leading-Eigenvector algorithm [81]. The value Qb, introduced by Barber [82], is calculated using the standard bipartite modularity function. To quantify the statistical significance of the nestedness (NODF) and modularity (Qb), 100 null random matrices (for each) were created with the null model Equiprobable (a random matrix in which all the interactions are uniformly permuted).

3. Results

Forty-nine E. huxleyi strains were characterized according to their ability to grow under a standard set of nutrients, light and temperature conditions. Growth rate (µ) varied significantly among E. huxleyi strains, ranging from 0.12 (SD ± 0.01) to 1.11 (SD ± 0.02)/d (registered in strains RCC4533 and RCC1744, respectively) (Figure S2). The difference in growth rate among the algal strains was not related to the ocean they were isolated from (one-way ANOVA F (2, 37) = 0.275, p = 0. 76).
We confronted the observed differences in resistance capacity with the parameters growth rate and viral production, respectively. The level of resistance of E. huxleyi to EhV infection was accessed in two manners: (R1) percentage of cells that were not lysed after incubation with viruses (Figure S3) and (R2) the number of EhV strains that successfully produced progeny on that host, meaning that lower R2 levels indicate higher resistance capacity. A trade-off between resistance and growth rate capacities (hypotheses 1 and 3 in Table 1) was not confirmed with our results. Neither types of resistance, R1 and R2, were significantly correlated to growth rate (Pearson’s r = −0.131, p = 0.370, and Pearson’s r = −0.0959, p = 0.512; respectively) (Figure 1 and Figure 2). R1 was indirectly correlated with viral production (Figure 3) (Pearson’s r = −0.499, p > 0.01), in accordance with hypothesis 4. R2 was significantly and positively correlated with maximum viral production (Pearson’s r = 0.614, p < 0.01), which means that the E. huxleyi strains that were susceptible to more EhV types were also the ones that presented higher maximum viral production (Figure 4). Viral production and growth rate did not correlate significantly (Pearson’s r = 0.1, p = 0.494) (Figure S4) and hence did not confirm hypothesis 2 (Table 1).
Seven out of the 49 E. huxleyi strains (RCC1259, RCC1269, RCC3856, 371, P847, PERU15-40 and SO52) were susceptible to infection by all the EhV strains tested, while 6 E. huxleyi strains (RCC1211, RCC1218, RCC1235, RCC1256, RCC1276 and RCC3548) were resistant to infection by all the EhV strains tested. When analysing these two groups of E. huxleyi strains, no significant differences in growth rate were found (one-way ANOVA F (1, 11) = 0.01592, p = 0.90188), while their R1 values were significantly different (one-way ANOVA F (1, 11) = 36.8593, p = 8.1 × 10−5).
A significant higher viral production was found in the most recently isolated algal strains (one-way ANOVA F (1, 47) = 30.36, p = 1.5 × 10−6). For the other parameters (growth rate, R1 and R2) there were no significant differences between younger and older strains (one-way ANOVA F (1, 47) = 1.094, p = 0.30; one-way ANOVA F (1, 47) = 0.106, p = 0.745; one-way ANOVA F (1, 47) = 0.909, p = 0.345; respectively).
We observed significant variation in “Maximum viral production” capacity among the different EhV strains (Figure 5). Those differences did not translate into significant differences in “Average Viral Production” (Figure S5), as the capacity of each EhV to produce progeny depended very much on which host strain it was infecting. Host-ranges among EhV strains also proved very variable, from generalists that infected up to 36 host strains (e.g., EhV-207) to specialists capable of infecting only 1 strain (e.g., EhV-99b1). Surprisingly, and against the prediction in hypothesis 5, generalist viral strains (EhV-164, EhV-202, EhV-208, EhV-201 and EhV-207) produced significantly more virus progeny viral production (one-way ANOVA F (1, 8) = 8.123, p = 0.021) than specialist strains (EhV-99b1, EhV-203, EhV-156, EhV-86, and EhV-145).
The bipartite network analysis applied to the whole host-range matrix displayed a nested structure (Figure 6) with a NODF value of 0.60. In that nested pattern there was a tendency for hosts with higher resistance to only be infected by more generalist viruses, while specialist viruses tend to infect the most sensitive hosts.

4. Discussion

Since Hutchinson first stated the Paradox of the phytoplankton in the early sixties, many hypotheses explaining the high diversity in the oceans have been postulated [13]. Among these, viral activity has proven to be a potential disrupter on equilibrium in planktonic communities [4,12]. Due to the lack of quantitative data for viral-host interactions, especially in marine micro-eukaryotic organisms, we therefore decided to perform a vast survey on strains of the ubiquitous and environmentally relevant coccolithophorid Emiliania huxleyi sp. (E. huxleyi) (Lohman) and its virus, Emiliania huxleyi virus (EhV), and investigate for emerging patterns resulting from this arms race.
Among the different hypotheses tested (Table 1) was the existence, or not, of a clear trade-off between resistance and growth rate (COR). COR has been previously confirmed in some bacteria-virus systems [26,27,28,29,30], and is fundamental in the formulation of the Killing the Winner model [19]. In our study we did not observe a clear COR trade-off in the E. huxleyi-EhV system. Instead, we found that highly resistant algal strains were capable of growing at high rates. This indicates that, at least in this system, viruses may not be the main selective force acting upon their hosts or, that if they are, their impact is camouflaged by antagonic impacts from other selective factors (e.g., different adaptation to the standard culture conditions used). However, it could be that viral-imposed selection was so strong that it would result in an emerging global cost of resistance observable on E. huxleyi strains independently of their inherent local adaptations. An approximation to such global “cost of resistance” is precisely the parameter value used when trying to model the interactions between viruses and their hosts [83]. Its prominence in current models justified the present attempt to evaluate its real extension.
When Avrani and colleagues [29] observed a similar response in viral resistant Prochlorococcus strains, they also found that the reduced growth rates increased after 7 months and that these strains reduced their resistance against the viruses [84]. The changes in growth rate and resistance occurred as independent events, indicating that the selection pressure on these phenotypes was decoupled. Decoupled selective pressure for growth rate and resistance may be the reason for the lack of correlation between these parameters in our study as well.
COR not being observed for the E. huxleyi-EhV system using our approach is not necessarily proving it does not exist or that it is irrelevant. As also tried in the current study, COR is often measured as reduction of growth rates in the resistant host [23,24,25,26,27], but other CORs, like altered susceptibility to other viruses and possibly also to some bacteria [85], have also been argued [29,84,86,87,88]. Trade-off might also emerge when strains with different resistance capacities are put under competition for a limited level of nutrients [30,33,89,90], and this is the logical follow up to our study. Another aspect to take into account is the potential impact that domestication has on the isolated strains [91]. In vitro growing conditions (nutrients, light, temperature) are inevitably different from what the cells would be experiencing in the natural environment. Particularly, in vitro cells are released from viral pressure, a situation that, with time, could potentially erase the selective traits that viruses might impose on cells in the natural environment. A sign of domestication-related effects in our case was the lower viral production capacity observed for “older” strains (isolated before 2009).
Patterns other than COR that shed light on the global interaction between E. huxleyi and EhV did, however, emerge in this study. Contrary to our expectations [14,68], we observed a tendency for generalist viruses (e.g., EhV-207) to produce more progeny than the specialists (e.g., EhV-86). It was recently reported that a generalist EhV strain could outcompete a specialist 8 h post infection [92]. One explanation for this apparent difference in infective success between generalist and specialist viruses may thus be a trade-off where high host-range/replication rates are associated with hindered progeny (new virions) fitness [64,93,94,95,96]. An alternative possibility could be the presence of an “un-costly” strong adaptive potential to new hosts, as shown for the Tobacco etch potyvirus (TEV) [97]. It also has to be taken into consideration that viral infective performance; such as viral adsorption coefficient and burst size also depends strongly on host traits. In the current study, a set of E. huxleyi strains were the ones that presented the higher viral production, independently of the EhV strain that was infecting them. Such added levels of complexity create niches for different strains of viruses and hosts with different infection and resistance capacities, respectively, to coexist. The patterns emerging from the interaction between E. huxleyi and EhV indicate that there’s a plethora of niches that create the possibility for co-existence of viruses and hosts with unexpected trait capacities. Notably, viral strains with narrower host-ranges and smaller virion production competing with generalist strains. Future studies should try to evaluate the possibility of take-over in the case of two specialist or generalist strains.
The emerging virus-host interaction network (VHIN) pattern showed a significant nestedness match between viral strains and their hosts. A nested structure like this is considered to result from sequential gene-for-gene (GFG) adaptations [98,99]. In the GFG model one genotype replaces another leading to continued fitness improvements of both, host and virus populations, resulting in an everlasting arms race dynamics. Different mesocosm studies on natural E. huxleyi/EhV communities [73,74] have shown that host and viral strain diversity can co-change in very short periods of just a few days during E. huxleyi blooms. This supports the Arms Race dynamics indicated by our VHIN. Future studies should evaluate the potential for strains with similar host-range capacity to take-over one another. The currently observed cross-infection network did not however have a perfect nested structure. An alternative co-evolution mechanism, termed diffuse co-evolution, appears to be more adequate for multi-species and/or multi-strain communities where selection pressures due to one species, can change in the presence of other species [17,100]. In order to predict diffuse co-evolution, however, experiments in which the different species/strains could interact, allowing real fitness costs associate to both, viruses and hosts, to arise [100] are necessary.
As also previously shown, the same E. huxleyi viruses (isolated in the English Channel and the Norwegian fjords) proved able to infect E. huxleyi hosts isolated in a large spatio-temporal scale [44,51], indicating a strong genomic consistency between geographically distant EhV strains. Nonetheless, and despite high abundance of conserved genomic sequences among these strains, significant genomic variety is also documented [44,73,101]. As EhVs are enveloped viruses [102], their entry mechanism should be endocytosis or fusion of the envelope with the host’s membrane and the progeny release through a budding mechanism [103]. Such an infection mechanism potentially generates a highly lipid-specific contact between host and virus. The host, E. huxleyi, has high phenotypic plasticity [104,105,106,107,108] and adaptation capacity [104,109,110,111,112] that could result in ecotypes that respond differently to viral infection [37,108,112]. Even if genes associated with virus susceptibility have been found within non-core regions of the E. huxleyi genome [108], our results did not show significant differences in growth rate, resistance, or viral production in hosts from very distant geographical locations. Hence, despite the recognized genetic variability in both host and virus, our results suggest a globally, non-segregated evolution process between E. huxleyi and EhV [113].
In conclusion, and despite a lack of supporting evidence of a trade-off between resistance and growth capacities, our results did indeed, through the nested host-virus interaction pattern, demonstrate a strong co-evolution pattern between E. huxleyi and EhV populations. The absence of trade-off between growth rate and resistance, invites us to think that EhVs may not be the main force driving the E. huxleyi selection, and that other fitness costs, which passed unnoticeably in the present study, exist. Further work should aim at unravelling these.

Supplementary Materials

The following are available online at www.mdpi.com/1999-4915/9/3/61/s1, Table S1: E. huxleyi strain information, in blank = No information; Table S2: EhV strain information; Table S3: Measured parameters from the cross-infectivity experiments between each E. huxleyi-EhV pair. fC = final concentration of E. huxleyi cells in cells/mL, µ = E. huxleyi growth rate, R1 = percentage of cells that were not lysed after incubation with viruses, compared to the non-inoculated controls, Vp = viral production in viral particles/mL, R2 = number of EhV strains that successfully produced progeny on that host, R1 AV = averaged R1 for each algal strain, Vp AV = averaged Vp for each algal strain in viral particles/mL, Max. Vp AV = averaged maximum Vp for each algal strain in viral particles/mL, SD = standard deviation; Figure S1: Correlation between viral production per host cell (Vp) and isolation year of the algal strains. Error bars show standard deviation (n = 13); Figure S2: Growth rates (µ/d) measured for control samples of all of the E. huxleyi strains (see Table S1 for strain information) used in the infection experiment measured over a period of x days and, calculated according to Levasseur et al. (1993). Values correspond to the control samples. Error bars show standard deviation (n = 3); Figure S3: Resistance strategy (R) for each E. huxleyi strain. Error bars show standard deviation (n = 13); Figure S4: Correlation between growth rate (µ) and viral production per host cell (Vp), in viral particles/mL. Error bars show standard deviation (n = 13); Figure S5: Average viral production per EhV strain, for all the algal strains. Error bars show standard deviation (n = 49).

Acknowledgments

The authors are grateful to William Wilson and Mike Allen for providing the EhV strains. We are also thankful to Tron Frede Thingstad for the helpful discussions and advice, and to the reviewers for their comments and suggestions that contributed to improve the manuscript. This work was supported by the Research Council of Norway, through project VIMPACT, no. 234197.

Author Contributions

A.P., E.R., and R.-A.S. conceived and designed the experiments; E.R. and M.O. performed the experiments; E.R. and A.P. analyzed the data; E.R., A.P., A.L., and R.-A.S. wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Wilhelm, S.W.; Suttle, C.A. Viruses and nutrient cycles in the sea—Viruses play critical roles in the structure and function of aquatic food webs. Bioscience 1999, 49, 781–788. [Google Scholar] [CrossRef]
  2. Suttle, C.A. Marine viruses—Major players in the global ecosystem. Nat. Rev. Microbiol. 2007, 5, 801–812. [Google Scholar] [CrossRef] [PubMed]
  3. Breitbart, M. Marine viruses: Truth or dare. Annu. Rev. Mar. Sci. 2012, 4, 425–448. [Google Scholar] [CrossRef] [PubMed]
  4. Fuhrman, J.A. Marine viruses and their biogeochemical and ecological effects. Nature 1999, 399, 541–548. [Google Scholar] [CrossRef] [PubMed]
  5. Wommack, K.E.; Colwell, R.R. Virioplankton: Viruses in aquatic ecosystems. Microbiol. Mol. Biol. Rev. 2000, 64, 69–114. [Google Scholar] [CrossRef] [PubMed]
  6. Weinbauer, M.G.; Rassoulzadegan, F. Are viruses driving microbial diversification and diversity? Environ. Microbiol. 2004, 6, 1–11. [Google Scholar] [CrossRef] [PubMed]
  7. Brussaard, C.P.D. Viral control of phytoplankton populations—A review. J. Eukaryot. Microbiol. 2004, 51, 125–138. [Google Scholar] [CrossRef] [PubMed]
  8. Suttle, C.A.; Chan, A.M. Dynamics and distribution of cyanophages and their effect on marine Synechococcus spp. Appl. Environ. Microbiol. 1994, 60, 3167–3174. [Google Scholar] [PubMed]
  9. Waterbury, J.B.; Valois, F.W. Resistance to cooccurring phages enables marine synechococcus communities to coexist with cyanophages abundant in seawater. Appl. Environ. Microbiol. 1993, 59, 3393–3399. [Google Scholar] [PubMed]
  10. Weinbauer, M.G. Ecology of prokaryotic viruses. FEMS Microbiol. Rev. 2004, 28, 127–181. [Google Scholar] [CrossRef] [PubMed]
  11. Jiang, S.C.; Paul, J.H. Gene transfer by transduction in the marine environment. Appl. Environ. Microbiol. 1998, 64, 2780–2787. [Google Scholar] [PubMed]
  12. Fuhrman, J.A.; Suttle, C.A. Viruses in marine planktonic systems. Oceanography 1993, 6, 51–63. [Google Scholar] [CrossRef]
  13. Hutchinson, G.E. The paradox of the plankton. Am. Nat. 1961, 95, 137–145. [Google Scholar] [CrossRef]
  14. Winter, C.; Bouvier, T.; Weinbauer, M.G.; Thingstad, T.F. Trade-offs between competition and defense specialists among unicellular planktonic organisms: The “killing the winner“ hypothesis revisited. Microbiol. Mol. Biol. Rev. 2010, 74, 42–57. [Google Scholar] [CrossRef] [PubMed]
  15. Thingstad, T.F.; Vage, S.; Storesund, J.E.; Sandaa, R.A.; Giske, J. A theoretical analysis of how strain-specific viruses can control microbial species diversity. Proc. Natl. Acad. Sci. USA 2014, 111, 7813–7818. [Google Scholar] [CrossRef] [PubMed]
  16. Sheldon, B.C.; Verhulst, S. Ecological immunology: Costly parasite defences and trade-offs in evolutionary ecology. Trends Ecol. Evol. 1996, 11, 317–321. [Google Scholar] [CrossRef]
  17. Weitz, J.S.; Poisot, T.; Meyer, J.R.; Flores, C.O.; Valverde, S.; Sullivan, M.B.; Hochberg, M.E. Phage-bacteria infection networks. Trends Microbiol. 2013, 21, 82–91. [Google Scholar] [CrossRef] [PubMed]
  18. Poullain, V.; Gandon, S.; Brockhurst, M.A.; Buckling, A.; Hochberg, M.E. The evolution of specificity in evolving and coevolving antagonistic interactions between a bacteria and its phage. Evolution 2008, 62, 1–11. [Google Scholar] [CrossRef] [PubMed]
  19. Meyer, J.R.; Dobias, D.T.; Weitz, J.S.; Barrick, J.E.; Quick, R.T.; Lenski, R.E. Repeatability and contingency in the evolution of a key innovation in phage lambda. Science 2012, 335, 428–432. [Google Scholar] [CrossRef] [PubMed]
  20. Fortuna, M.A.; Stouffer, D.B.; Olesen, J.M.; Jordano, P.; Mouillot, D.; Krasnov, B.R.; Poulin, R.; Bascompte, J. Nestedness versus modularity in ecological networks: Two sides of the same coin? J. Anim. Ecol. 2010, 79, 811–817. [Google Scholar] [CrossRef] [PubMed]
  21. Flores, C.O.; Poisot, T.; Valverde, S.; Weitz, J.S. Bimat: A matlab package to facilitate the analysis of bipartite networks. Methods Ecol. Evol. 2016, 7, 127–132. [Google Scholar] [CrossRef]
  22. Gomez, P.; Buckling, A. Bacteria-phage antagonistic coevolution in soil. Science 2011, 332, 106–109. [Google Scholar] [CrossRef] [PubMed]
  23. Benmayor, R.; Buckling, A.; Bonsall, M.B.; Brockhurst, M.A.; Hodgson, D.J. The interactive effects of parasitesf disturbance, and productivity on experimental adaptive radiations. Evolution 2008, 62, 467–477. [Google Scholar] [CrossRef] [PubMed]
  24. Brockhurst, M.A.; Buckling, A.; Rainey, P.B. The effect of a bacteriophage on diversification of the opportunistic bacterial pathogen, pseudomonas aeruginosa. Proc. R. Soc. B-Biol. Sci. 2005, 272, 1385–1391. [Google Scholar] [CrossRef] [PubMed]
  25. Brockhurst, M.A.; Rainey, P.B.; Buckling, A. The effect of spatial heterogeneity and parasites on the evolution of host diversity. Proc. R. Soc. B-Biol. Sci. 2004, 271, 107–111. [Google Scholar] [CrossRef] [PubMed]
  26. Middelboe, M.; Holmfeldt, K.; Riemann, L.; Nybroe, O.; Haaber, J. Bacteriophages drive strain diversification in a marine flavobacterium: Implications for phage resistance and physiological properties. Environ. Microbiol. 2009, 11, 1971–1982. [Google Scholar] [CrossRef] [PubMed]
  27. Riemann, L.; Grossart, H.-P. Elevated lytic phage production as a consequence of particle colonization by a marine flavobacterium (Cellulophaga sp.). Microb. Ecol. 2008, 56, 505–512. [Google Scholar] [CrossRef] [PubMed]
  28. Middelboe, M. Bacterial growth rate and marine virus-host dynamics. Microb. Ecol. 2000, 40, 114–124. [Google Scholar] [PubMed]
  29. Avrani, S.; Wurtzel, O.; Sharon, I.; Sorek, R.; Lindell, D. Genomic island variability facilitates prochlorococcus-virus coexistence. Nature 2011, 474, 604–608. [Google Scholar] [CrossRef] [PubMed]
  30. Meaden, S.; Paszkiewicz, K.; Koskella, B. The cost of phage resistance in a plant pathogenic bacterium is context-dependent. Evolution 2015, 69, 1321–1328. [Google Scholar] [CrossRef] [PubMed]
  31. Lennon, J.T.; Khatana, S.A.M.; Marston, M.F.; Martiny, J.B.H. Is there a cost of virus resistance in marine cyanobacteria? ISME J. 2007, 1, 300–312. [Google Scholar] [CrossRef] [PubMed]
  32. Lythgoe, K.A.; Chao, L. Mechanisms of coexistence of a bacteria and a bacteriophage in a spatially homogeneous environment. Ecol. Lett. 2003, 6, 326–334. [Google Scholar] [CrossRef]
  33. Bohannan, B.J.M.; Lenski, R.E. The relative importance of competition and predation varies with productivity in a model community. Am. Nat. 2000, 156, 329–340. [Google Scholar] [CrossRef]
  34. Mizoguchi, K.; Morita, M.; Fischer, C.R.; Yoichi, M.; Tanji, Y.; Unno, H. Coevolution of bacteriophage pp01 and escherichia coli o157: H7 in continuous culture. Appl. Environ. Microbiol. 2003, 69, 170–176. [Google Scholar] [CrossRef] [PubMed]
  35. Thomas, R.; Grimsley, N.; Escande, M.L.; Subirana, L.; Derelle, E.; Moreau, H. Acquisition and maintenance of resistance to viruses in eukaryotic phytoplankton populations. Environ. Microbiol. 2011, 13, 1412–1420. [Google Scholar] [CrossRef] [PubMed]
  36. Frickel, J.; Sieber, M.; Becks, L. Eco-evolutionary dynamics in a coevolving host-virus system. Ecol. Lett. 2016, 19, 450–459. [Google Scholar] [CrossRef] [PubMed]
  37. Brown, C.W.; Yoder, J.A. Coccolithophorid blooms in the global ocean. J. Geophys. Res.-Oceans 1994, 99, 7467–7482. [Google Scholar] [CrossRef]
  38. Tyrrell, T.; Taylor, A.H. A modelling study of Emiliania huxleyi in the ne atlantic. J. Mar. Syst. 1996, 9, 83–112. [Google Scholar] [CrossRef]
  39. Egge, J.K.; Heimdal, B.R. Blooms of phytoplankton including Emiliania-huxleyi (haptophyta)—Effects of nutrient supply in different n-p ratios. Sarsia 1994, 79, 333–348. [Google Scholar] [CrossRef]
  40. Broerse, A.T.C.; Ziveri, P.; van Hinte, J.E.; Honjo, S. Coccolithophore export production, species composition, and coccolith-caco3 fluxes in the ne atlantic (34 degrees n 21 degrees w and 48 degrees n 21 degrees w). Deep-Sea Res. Part II Top. Stud. Oceanogr. 2000, 47, 1877–1905. [Google Scholar] [CrossRef]
  41. Burkill, P.H.; Archer, S.D.; Robinson, C.; Nightingale, P.D.; Groom, S.B.; Tarran, G.A.; Zubkov, M.V. Dimethyl sulphide biogeochemistry within a coccolithophore bloom (disco): An overview. Deep-Sea Res. Part II Top. Stud. Oceanogr. 2002, 49, 2863–2885. [Google Scholar] [CrossRef]
  42. Cottrell, M.T.; Suttle, C.A. Wide-spread occurrence and clonal variation in viruses which cause lysis of a cosmopolitan, eukaryotic marine phytoplankter, micromonas-pusilla. Mar. Ecol. Prog. Ser. 1991, 78, 1–9. [Google Scholar] [CrossRef]
  43. Nissimov, J.I.; Napier, J.A.; Kimmance, S.A.; Allen, M.J. Permanent draft genomes of four new coccolithoviruses: Ehv-18, ehv-145, ehv-156 and ehv-164. Mar. Genom. 2014, 15, 7–8. [Google Scholar] [CrossRef] [PubMed]
  44. Allen, M.J.; Martinez-Martinez, J.; Schroeder, D.C.; Somerfield, P.J.; Wilson, W.H. Use of microarrays to assess viral diversity: From genotype to phenotype. Environ. Microbiol. 2007, 9, 971–982. [Google Scholar] [CrossRef] [PubMed]
  45. Pagarete, A.; Kusonmano, K.; Petersen, K.; Kimmance, S.A.; Martinez, J.M.; Wilson, W.H.; Hehemann, J.H.; Allen, M.J.; Sandaa, R.A. Dip in the gene pool: Metagenomic survey of natural coccolithovirus communities. Virology 2014, 466, 129–137. [Google Scholar] [CrossRef] [PubMed]
  46. Pagarete, A.; Lanzen, A.; Puntervoll, P.; Sandaa, R.A.; Larsen, A.; Larsen, J.B.; Allen, M.J.; Bratbak, G. Genomic sequence and analysis of ehv-99b1, a new coccolithovirus from the norwegian fjords. Intervirology 2013, 56, 60–66. [Google Scholar] [CrossRef] [PubMed]
  47. Nissimov, J.I.; Worthy, C.A.; Rooks, P.; Napier, J.A.; Kimmance, S.A.; Henn, M.R.; Ogata, H.; Allen, M.J. Draft genome sequence of the coccolithovirus ehv-84. Stand. Genom. Sci. 2011, 5, 1–11. [Google Scholar] [CrossRef] [PubMed]
  48. Nissimov, J.I.; Worthy, C.A.; Rooks, P.; Napier, J.A.; Kimmance, S.A.; Henn, M.R.; Ogata, H.; Allen, M.J. Draft genome sequence of four coccolithoviruses: Emiliania huxleyi virus ehv-88, ehv-201, ehv-207, and ehv-208. J. Virol. 2012, 86, 2896–2897. [Google Scholar] [CrossRef] [PubMed]
  49. Nissimov, J.I.; Worthy, C.A.; Rooks, P.; Napier, J.A.; Kimmance, S.A.; Henn, M.R.; Ogata, H.; Allen, M.J. Draft genome sequence of the coccolithovirus Emiliania huxleyi virus 202. J. Virol. 2012, 86, 2380–2381. [Google Scholar] [CrossRef] [PubMed]
  50. Nissimov, J.I.; Worthy, C.A.; Rooks, P.; Napier, J.A.; Kimmance, S.A.; Henn, M.R.; Ogata, H.; Allen, M.J. Draft genome sequence of the coccolithovirus Emiliania huxleyi virus 203. J. Virol. 2011, 85, 13468–13469. [Google Scholar] [CrossRef] [PubMed]
  51. Pagarete, A.J. Functional Genomics of Coccolithophore Viruses. Ph.D. Thesis, Pierre and Marie Curie University, Paris, France, 2010; p. 222. [Google Scholar]
  52. Middelboe, M.; Jacquet, S.; Weinbauer, M. Viruses in freshwater ecosystems: An introduction to the exploration of viruses in new aquatic habitats. Freshw. Biol. 2008, 53, 1069–1075. [Google Scholar] [CrossRef]
  53. Bouvier, T.; del Giorgio, P.A. Key role of selective viral-induced mortality in determining marine bacterial community composition. Environ. Microbiol. 2007, 9, 287–297. [Google Scholar] [CrossRef] [PubMed]
  54. Chao, L.; Levin, B.R.; Stewart, F.M. A complex community in a simple habitat: An experimental study with bacteria and phage. Ecology 1977, 58, 369–378. [Google Scholar] [CrossRef]
  55. Lenski, R.E.; Levin, B.R. Constraints on the coevolution of bacteria and virulent phage: A model, some experiments, and predictions for natural communities. Am. Nat. 1985, 125, 585–602. [Google Scholar] [CrossRef]
  56. Thyrhaug, R.; Larsen, A.; Thingstad, T.F.; Bratbak, G. Stable coexistence in marine algal host-virus systems. Mar. Ecol. Prog. Ser. 2003, 254, 27–35. [Google Scholar] [CrossRef]
  57. Haaber, J.; Middelboe, M. Viral lysis of phaeocystis pouchetii: Implications for algal population dynamics and heterotrophic c, n and p cycling. ISME J. 2009, 3, 430–441. [Google Scholar] [CrossRef] [PubMed]
  58. Moebus, K. Marine bacteriophage reproduction under nutrient-limited growth of host bacteria.2. Investigations with phage-host system h3:H3/1. Mar. Ecol. Prog. Ser. 1996, 144, 13–22. [Google Scholar] [CrossRef]
  59. Moebus, K. Marine bacteriophage reproduction under nutrient-limited growth of host bacteria.1. Investigations with six phage-host systems. Mar. Ecol. Prog. Ser. 1996, 144, 1–12. [Google Scholar] [CrossRef]
  60. Parada, V.; Herndl, G.J.; Weinbauer, M.G. Viral burst size of heterotrophic prokaryotes in aquatic systems. J. Mar. Biol. Assoc. UK 2006, 86, 613–621. [Google Scholar] [CrossRef]
  61. Motegi, C.; Nagata, T. Enhancement of viral production by addition of nitrogen or nitrogen plus carbon in subtropical surface waters of the South Pacific. Aquat. Microb. Ecol. 2007, 48, 27–34. [Google Scholar] [CrossRef]
  62. Bratbak, G.; Jacobsen, A.; Heldal, M.; Nagasaki, K.; Thingstad, F. Virus production in phaeocystis pouchetii and its relation to host cell growth and nutrition. Aquat. Microb. Ecol. 1998, 16, 1–9. [Google Scholar] [CrossRef]
  63. Baudoux, A.-C.; Brussaard, C.P.D. Influence of irradiance on virus-algal host interactions. J. Phycol. 2008, 44, 902–908. [Google Scholar] [CrossRef] [PubMed]
  64. Demory, D.; Arsenieff, L.; Simon, N.; Six, C.; Rigaut-Jalabert, F.; Marie, D.; Ge, P.; Bigeard, E.; Jacquet, S.; Sciandra, A.; et al. Temperature is a key factor in micromonas–virus interactions. ISME J. 2017, 13, 601–612. [Google Scholar] [CrossRef] [PubMed]
  65. Maat, D.S.; Crawfurd, K.J.; Timmermans, K.R.; Brussaard, C.P.D. Elevated CO2 and phosphate limitation favor micromonas pusilla through stimulated growth and reduced viral impact. Appl. Environ. Microbiol. 2014, 80, 3119–3127. [Google Scholar] [CrossRef] [PubMed]
  66. Maat, D.S.; de Blok, R.; Brussaard, C.P.D. Combined phosphorus limitation and light stress prevent viral proliferation in the phytoplankton species phaeocystis globosa, but not in micromonas pusilla. Front. Mar. Sci. 2016, 3. [Google Scholar] [CrossRef]
  67. Kendrick, B.J.; DiTullio, G.R.; Cyronak, T.J.; Fulton, J.M.; van Mooy, B.A.S.; Bidle, K.D. Temperature-induced viral resistance in Emiliania huxleyi (prymnesiophyceae). PLoS ONE 2014, 9, e112134. [Google Scholar] [CrossRef] [PubMed]
  68. MacArthur, R.H.; Wilson, E.O. The Theory of Island Biogeography; Princeton University Press: Princeton, NJ, USA, 1967; p. 203. [Google Scholar]
  69. Klochkova, T.A.; Kang, S.H.; Cho, G.Y.; Pueschel, C.M.; West, J.A.; Kim, G.H. Biology of a terrestrial green alga, Chlorococcum sp. (chlorococcales, chlorophyta), collected from the miruksazi stupa in Korea. Phycologia 2006, 45, 349–358. [Google Scholar] [CrossRef]
  70. Wilson, W.H.; Tarran, G.A.; Schroeder, D.; Cox, M.; Oke, J.; Malin, G. Isolation of viruses responsible for the demise of an Emiliania huxleyi bloom in the english channel. J. Mar. Biol. Assoc. UK 2002, 82, 369–377. [Google Scholar] [CrossRef]
  71. Vaughn, J.M.; Balch, W.M.; Novotny, J.F.; Vining, C.L.; Palmer, C.D.; Drapeau, D.T.; Booth, E.; Vaughn, J.M.; Kneifel, D.M.; Bell, A.L. Isolation of Emiliania huxleyi viruses from the gulf of maine. Aquat. Microb. Ecol. 2010, 58, 109–116. [Google Scholar] [CrossRef]
  72. Highfield, A.; Evans, C.; Walne, A.; Miller, P.I.; Schroeder, D.C. How many coccolithovirus genotypes does it take to terminate an Emiliania huxleyi bloom? Virology 2014, 466, 138–145. [Google Scholar] [CrossRef] [PubMed]
  73. Martinez, J.M.; Schroeder, D.C.; Larsen, A.; Bratbak, G.; Wilson, W.H. Molecular dynamics of Emiliania huxleyi and cooccurring viruses during two separate mesocosm studies. Appl. Environ. Microbiol. 2007, 73, 554–562. [Google Scholar] [CrossRef] [PubMed]
  74. Sorensen, G.; Baker, A.C.; Hall, M.J.; Munn, C.B.; Schroeder, D.C. Novel virus dynamics in an Emiliania huxleyi bloom. J. Plankton Res. 2009, 31, 787–791. [Google Scholar] [CrossRef] [PubMed]
  75. Ruiz, E. Growth Rate Experiments in Emiliania Huxleyi; Universitetet I Bergen: Bergen, Norway, 2014. [Google Scholar]
  76. Brussaard, C.P.D.; Thyrhaug, R.; Marie, D.; Bratbak, G. Flow cytometric analyses of viral infection in two marine phytoplankton species, micromonas pusilla (prasinophyceae) and phaeocystis pouchetii (prymnesiophyceae). J. Phycol. 1999, 35, 941–948. [Google Scholar] [CrossRef]
  77. Brussaard, C.P.D.; Marie, D.; Bratbak, G. Flow cytometric detection of viruses. J. Virol. Methods 2000, 85, 175–182. [Google Scholar] [CrossRef]
  78. Brussaard, C.P.D.; Marie, D.; Thyrhaug, R.; Bratbak, G. Flow cytometric analysis of phytoplankton viability following viral infection. Aquat. Microb. Ecol. 2001, 26, 157–166. [Google Scholar] [CrossRef]
  79. Levasseur, M.; Thompson, P.A.; Harrison, P.J. Physiological acclimation of marine-phytoplankton to different nitrogen-sources. J. Phycol. 1993, 29, 587–595. [Google Scholar] [CrossRef]
  80. Almeida-Neto, M.; Guimaraes, P.; Guimaraes, P.R.; Loyola, R.D.; Ulrich, W. A consistent metric for nestedness analysis in ecological systems: Reconciling concept and measurement. Oikos 2008, 117, 1227–1239. [Google Scholar] [CrossRef]
  81. Newman, M.E.J. Modularity and community structure in networks. Proc. Natl. Acad. Sci. USA 2006, 103, 8577–8582. [Google Scholar] [CrossRef] [PubMed]
  82. Barber, M.J. Modularity and community detection in bipartite networks. Phys. Rev. E 2007, 76. [Google Scholar] [CrossRef] [PubMed]
  83. Våge, S.; Storesund, J.E.; Thingstad, T.F. Adding a cost of resistance description extends the ability of virus-host model to explain observed patterns in structure and function of pelagic microbial communities. Environ. Microbiol. 2013, 15, 1842–1852. [Google Scholar] [CrossRef] [PubMed]
  84. Avrani, S.; Lindell, D. Convergent evolution toward an improved growth rate and a reduced resistance range in prochlorococcus strains resistant to phage. Proc. Natl. Acad. Sci. USA 2015, 112, E2191–E2200. [Google Scholar] [CrossRef] [PubMed]
  85. Segev, E.; Wyche, T.P.; Kim, K.H.; Petersen, J.; Ellebrandt, C.; Vlamakis, H.; Barteneva, N.; Paulson, J.N.; Chai, L.; Clardy, J.; et al. Dynamic metabolic exchange governs a marine algal-bacterial interaction. Elife 2016, 5, e17473. [Google Scholar] [CrossRef] [PubMed]
  86. Marston, M.F.; Pierciey, F.J.; Shepard, A.; Gearin, G.; Qi, J.; Yandava, C.; Schuster, S.C.; Henn, M.R.; Martiny, J.B.H. Rapid diversification of coevolving marine synechococcus and a virus. Proc. Natl. Acad. Sci. USA 2012, 109, 4544–4549. [Google Scholar] [CrossRef] [PubMed]
  87. Avrani, S.; Schwartz, D.A.; Lindell, D. Virus-host swinging party in the oceans: Incorporating biological complexity into paradigms of antagonistic coexistence. Mob. Genet. Elements 2012, 2, 88–95. [Google Scholar] [CrossRef] [PubMed]
  88. Castillo, D.; Christiansen, R.H.; Espejo, R.; Middelboe, M. Diversity and geographical distribution of flavobacterium psychrophilum isolates and their phages: Patterns of susceptibility to phage infection and phage host range. Microb. Ecol. 2014, 67, 748–757. [Google Scholar] [CrossRef] [PubMed]
  89. Holt, R.D.; Grover, J.; Tilman, D. Simple rules for interspecific dominance in systems with exploitative and apparent competition. Am. Nat. 1994, 144, 741–771. [Google Scholar] [CrossRef]
  90. Leibold, M.A. A graphical model of keystone predators in food webs: Trophic regulation of abundance, incidence, and diversity patterns in communities. Am. Nat. 1996, 147, 784–812. [Google Scholar] [CrossRef]
  91. Lakeman, M.B.; von Dassow, P.; Cattolico, R.A. The strain concept in phytoplankton ecology. Harmful Algae 2009, 8, 746–758. [Google Scholar] [CrossRef]
  92. Nissimov, J.I.; Napier, J.A.; Allen, M.J.; Kimmance, S.A. Intragenus competition between coccolithoviruses: An insight on how a select few can come to dominate many. Environ. Microbiol. 2016, 18, 133–145. [Google Scholar] [CrossRef] [PubMed]
  93. Duffy, S.; Turner, P.E.; Burch, C.L. Pleiotropic costs of niche expansion in the RNA bacteriophage phi 6. Genetics 2006, 172, 751–757. [Google Scholar] [CrossRef] [PubMed]
  94. Elena, S.F.; Agudelo-Romero, P.; Lalic, J. The evolution of viruses in multi-host fitness landscapes. Open Virol. J. 2009, 3, 1–6. [Google Scholar] [CrossRef] [PubMed]
  95. Nikolin, V.M.; Osterrieder, K.; von Messling, V.; Hofer, H.; Anderson, D.; Dubovi, E.; Brunner, E.; East, M.L. Antagonistic pleiotropy and fitness trade-offs reveal specialist and generalist traits in strains of canine distemper virus. PLoS ONE 2012, 7, e50955. [Google Scholar] [CrossRef]
  96. Keen, E.C. Tradeoffs in bacteriophage life histories. Bacteriophage 2014, 4, e28365. [Google Scholar] [CrossRef] [PubMed]
  97. Bedhomme, S.; Lafforgue, G.; Elena, S.F. Multihost experimental evolution of a plant rna virus reveals local adaptation and host-specific mutations. Mol. Biol. Evol. 2012, 29, 1481–1492. [Google Scholar] [CrossRef] [PubMed]
  98. Agrawal, A.F.; Lively, C.M. Modelling infection as a two-step process combining gene-for-gene and matching-allele genetics. Proc. R. Soc. B-Biol. Sci. 2003, 270, 323–334. [Google Scholar] [CrossRef] [PubMed]
  99. Flor, H.H. Host-parasite interaction in flax rust—Its genetics and other implications. Phytopathology 1955, 45, 680–685. [Google Scholar]
  100. Inouye, B.; Stinchcombe, J.R. Relationships between ecological interaction modifications and diffuse coevolution: Similarities, differences, and causal links. Oikos 2001, 95, 353–360. [Google Scholar] [CrossRef]
  101. Rowe, J.M.; Fabre, M.F.; Gobena, D.; Wilson, W.H.; Wilhelm, S.W. Application of the major capsid protein as a marker of the phylogenetic diversity of Emiliania huxleyi viruses. FEMS Microbiol. Ecol. 2011, 76, 373–380. [Google Scholar] [CrossRef] [PubMed]
  102. Mackinder, L.C.M.; Worthy, C.A.; Biggi, G.; Hall, M.; Ryan, K.P.; Varsani, A.; Harper, G.M.; Wilson, W.H.; Brownlee, C.; Schroeder, D.C. A unicellular algal virus, Emiliania huxleyi virus 86, exploits an animal-like infection strategy. J. Gen. Virol. 2009, 90, 2306–2316. [Google Scholar] [CrossRef] [PubMed]
  103. Allen, M.J.; Howard, J.A.; Lilley, K.S.; Wilson, W.H. Proteomic analysis of the ehv-86 virion. Proteome Sci. 2008, 6. [Google Scholar] [CrossRef] [PubMed]
  104. Medlin, L.K.; Barker, G.L.A.; Campbell, L.; Green, J.C.; Hayes, P.K.; Marie, D.; Wrieden, S.; Vaulot, D. Genetic characterisation of Emiliania huxleyi (haptophyta). J. Mar. Syst. 1996, 9, 13–31. [Google Scholar] [CrossRef]
  105. Iglesias-Rodriguez, M.D.; Schofield, O.M.; Batley, J.; Medlin, L.K.; Hayes, P.K. Intraspecific genetic diversity in the marine coccolithophore Emiliania huxleyi (prymnesiophyceae): The use of microsatellite analysis in marine phytoplankton population studies. J. Phycol. 2006, 42, 526–536. [Google Scholar] [CrossRef]
  106. Blanco-Ameijeiras, S.; Lebrato, M.; Stoll, H.M.; Iglesias-Rodriguez, D.; Muller, M.N.; Mendez-Vicente, A.; Oschlies, A. Phenotypic variability in the coccolithophore Emiliania huxleyi. PLoS ONE 2016, 11, e0157697. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  107. Paasche, E. A review of the coccolithophorid Emiliania huxleyi (prymnesiophyceae), with particular reference to growth, coccolith formation, and calcification-photosynthesis interactions. Phycologia 2001, 40, 503–529. [Google Scholar] [CrossRef]
  108. Kegel, J.U.; John, U.; Valentin, K.; Frickenhaus, S. Genome variations associated with viral susceptibility and calcification in Emiliania huxleyi. PLoS ONE 2013, 8, e80684. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  109. Young, J.R.; Westbroek, P. Genotypic variation in the coccolithophorid species Emiliania-huxleyi. Mar. Micropaleontol. 1991, 18, 5–23. [Google Scholar] [CrossRef]
  110. Young, J.R. A Guide to Extant Coccolithophore Taxonomy; International Nannoplankton Association: Bremerhaven, Germany, 2003; p. 121. [Google Scholar]
  111. Cook, S.S.; Whitlock, L.; Wright, S.W.; Hallegraeff, G.M. Photosynthetic pigment and genetic differences between two southern ocean morphotypes of Emiliania huxleyi (haptophyta). J. Phycol. 2011, 47, 615–626. [Google Scholar] [CrossRef] [PubMed]
  112. Hagino, K.; Bendif, E.; Young, J.R.; Kogame, K.; Probert, I.; Takano, Y.; Horiguchi, T.; de Vargas, C.; Okada, H. New evidence for morphological and genetic variation in the cosmopolitan coccolithophore Emiliania huxleyi (prymnesiophyceae) from the cox1b-atp4 genes. J. Phycol. 2011, 47, 1164–1176. [Google Scholar] [CrossRef] [PubMed]
  113. Coolen, M.J.L. 7000 years of Emiliania huxleyi viruses in the black sea. Science 2011, 333, 451–452. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Resistance capacity R1 (calculated as the ratio between the number of cells that did not lyse after incubation with viruses and the number of cells in the non-inoculated controls) plotted against growth rate (μ). Error bars show standard deviation (n = 3).
Figure 1. Resistance capacity R1 (calculated as the ratio between the number of cells that did not lyse after incubation with viruses and the number of cells in the non-inoculated controls) plotted against growth rate (μ). Error bars show standard deviation (n = 3).
Viruses 09 00061 g001
Figure 2. Resistance capacity R2 (number of viral strains infecting each algal strain) plotted against growth rate (μ). Error bars show standard deviation (n = 3).
Figure 2. Resistance capacity R2 (number of viral strains infecting each algal strain) plotted against growth rate (μ). Error bars show standard deviation (n = 3).
Viruses 09 00061 g002
Figure 3. Viral production (Vp) plotted against resistance capacity R1. Error bars show standard deviation (n = 13).
Figure 3. Viral production (Vp) plotted against resistance capacity R1. Error bars show standard deviation (n = 13).
Viruses 09 00061 g003
Figure 4. Number of viral strains infecting each algal strain and maximum viral production correlation. Error bars show standard deviation (n = 3).
Figure 4. Number of viral strains infecting each algal strain and maximum viral production correlation. Error bars show standard deviation (n = 3).
Viruses 09 00061 g004
Figure 5. Differences between maximum viral production among EhV strains. Error bars show standard deviation (n = 49).
Figure 5. Differences between maximum viral production among EhV strains. Error bars show standard deviation (n = 49).
Viruses 09 00061 g005
Figure 6. Viral-host infectivity network with a clear nested pattern (NODF value of 0.60) where specialist viruses tend to infect the most susceptible hosts, while viruses with broader host-range infect hosts with higher resistance. ■: infection; □: no infection. Sidebars represent μ, R1 and Vp parameters, respectively.
Figure 6. Viral-host infectivity network with a clear nested pattern (NODF value of 0.60) where specialist viruses tend to infect the most susceptible hosts, while viruses with broader host-range infect hosts with higher resistance. ■: infection; □: no infection. Sidebars represent μ, R1 and Vp parameters, respectively.
Viruses 09 00061 g006
Table 1. Hypotheses tested in the current study based on outcome of previous virus-host interaction studies. µ: growth rate; R: resistance; Vp: viral production.
Table 1. Hypotheses tested in the current study based on outcome of previous virus-host interaction studies. µ: growth rate; R: resistance; Vp: viral production.
NumberHypothesisReference
1Resistance is associated with reduced growth rates (COR).Prokaryotes: [23,24,25,27,28,29,52,53,54,55]
Eukaryotes: [35,36,56,57]
2Host strains with higher µ produce more viruses.[58,59,60,61,62,63,64,65,66]
3Host strains with higher µ are infected by more viral strains.[36]
4Host strains with higher R produce fewer viruses.[56,67]
5Specialist viruses have higher Vp than generalists.[14,68]

Share and Cite

MDPI and ACS Style

Ruiz, E.; Oosterhof, M.; Sandaa, R.-A.; Larsen, A.; Pagarete, A. Emerging Interaction Patterns in the Emiliania huxleyi-EhV System. Viruses 2017, 9, 61. https://doi.org/10.3390/v9030061

AMA Style

Ruiz E, Oosterhof M, Sandaa R-A, Larsen A, Pagarete A. Emerging Interaction Patterns in the Emiliania huxleyi-EhV System. Viruses. 2017; 9(3):61. https://doi.org/10.3390/v9030061

Chicago/Turabian Style

Ruiz, Eliana, Monique Oosterhof, Ruth-Anne Sandaa, Aud Larsen, and António Pagarete. 2017. "Emerging Interaction Patterns in the Emiliania huxleyi-EhV System" Viruses 9, no. 3: 61. https://doi.org/10.3390/v9030061

APA Style

Ruiz, E., Oosterhof, M., Sandaa, R. -A., Larsen, A., & Pagarete, A. (2017). Emerging Interaction Patterns in the Emiliania huxleyi-EhV System. Viruses, 9(3), 61. https://doi.org/10.3390/v9030061

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop