# Life-Cycle and Genome of OtV5, a Large DNA Virus of the Pelagic Marine Unicellular Green Alga Ostreococcus tauri

Evelyne Derelle, Conchita Ferraz, Marie-Line Escande, Sophie Eychenié,

Richard Cooke, Gwenael Piganeau, Yves Desdevises, Laure Bellec, Hervé

Moreau, Nigel Grimsley

► To cite this version:

Evelyne Derelle, Conchita Ferraz, Marie-Line Escande, Sophie Eychenié, Richard Cooke, et al.. Life-Cycle and Genome of OtV5, a Large DNA Virus of the Pelagic Marine Unicellular Green Alga Ostreococcus tauri. PLoS ONE, 2008, 3 (5), pp.e2250. 10.1371/journal.pone.0002250. [hal-02114367](https://hal.science/hal-02114367)

HAL Id: hal-02114367

https://hal.science/hal-02114367

Submitted on 25 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Life-Cycle and Genome of OtV5, a Large DNA Virus of the Pelagic Marine Unicellular Green Alga *Ostreococcus tauri*

Evelyne Derelle¹,², Conchita Ferraz³, Marie-Line Escande¹,², Sophie Eychenié³, Richard Cooke⁴, Gwenaël Piganeau¹,², Yves Desdevises¹, Laure Bellec¹, Hervé Moreau¹,², Nigel Grimsley¹,²*

1 Université Pierre et Marie Curie-Paris 06, Laboratoire Arago, Banyuls-sur-Mer, France, 2 CNRS, UMR7628, Laboratoire Arago, Banyuls-sur-Mer, France, 3 Institut de Génétique Humaine, Génopole Montpellier Languedoc-Roussillon, UPR1142, Montpellier, France, 4 Génopole Languedoc-Roussillon, Génome et Développement de Plantes, UMR5096, Perpignan, France

## Abstract

Large DNA viruses are ubiquitous, infecting diverse organisms ranging from algae to man, and have probably evolved from an ancient common ancestor. In aquatic environments, such algal viruses control blooms and shape the evolution of biodiversity in phytoplankton, but little is known about their biological functions. We show that *Ostreococcus tauri*, the smallest known marine photosynthetic eukaryote, whose genome is completely characterized, is a host for large DNA viruses, and present an analysis of the life-cycle and 186,234 bp long linear genome of OtV5. OtV5 is a lytic phycodnavirus which unexpectedly does not degrade its host chromosomes before the host cell bursts. Analysis of its complete genome sequence confirmed that it lacks expected site-specific endonucleases, and revealed the presence of 16 genes whose predicted functions are novel to this group of viruses. OtV5 carries at least one predicted gene whose protein closely resembles its host counterpart and several other host-like sequences, suggesting that horizontal gene transfers between host and viral genomes may occur frequently on an evolutionary scale. Fifty seven percent of the 268 predicted proteins present no similarities with any known protein in Genbank, underlining the wealth of undiscovered biological diversity present in oceanic viruses, which are estimated to harbour 200Mt of carbon.

**Citation:** Derelle E, Ferraz C, Escande M-L, Eychenié S, Cooke R, et al. (2008) Life-Cycle and Genome of OtV5, a Large DNA Virus of the Pelagic Marine Unicellular Green Alga *Ostreococcus tauri*. PLoS ONE 3(5): e2250. doi:10.1371/journal.pone.0002250

**Editor:** Naomi Ward, University of Wyoming, United States of America

**Received:** January 10, 2008; **Accepted:** April 2, 2008; **Published:** May 28, 2008

**Copyright:** © 2008 Derelle et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** This work was supported partly by an "Agence National de Recherche" grant "PICOVIR" (coordinator N.G.), partly by the European Network of Excellence "Marine Genomics Europe" (no. GOCE-CT-2004-505403, coordinator H.M.) and the isolation of new strains was partly supported by a GDR 1928 "Genomique des Populations" grant to G. P.

**Competing Interests:** The authors have declared that no competing interests exist.

* E-mail: nigel.grimsley@obs-banyuls.fr

## Introduction

Viruses are by far the most abundant biological entities in the oceans [1,2], but relatively little is known about their life cycles and how they influence the diverse worldwide populations of plankton. Unicellular phytoplankton are at the base of the food web and recent evidence suggests that their populations are controlled to a large extent by viruses [2-7]. Indeed, some authors suggest that viral infections drive diversity and speciation in the microbial world [8]. In order to gain some understanding of how such global ecosystems function, one important approach is to analyze in detail representative components of these interactions, and then propose hypotheses applicable to the behavior of the whole system. Historical experience with model systems shows that the choice of the interaction investigated is crucial, since our depth of understanding will be limited by the amenability of the partners in question to experimental analyses.

We focus this study on a virus of *Ostreococcus tauri*, the latter being the smallest eukaryotic free-living photosynthetic organism known, whose complete genome sequence has been analysed [9]. *O. tauri*, first isolated from the Thau Lagoon in 1994 [10,11], is a free-living picoeukaryotic green alga (cell diameter ~1 µm) with a minimal cellular structure and a relatively fast cell cycle in culture (2-3 divisions/day in continuous light). Its single chloroplast, single mitochondrion and cytoplasm are bounded by a membrane lacking

any detectable cell wall or scales. In the marine environment, *Ostreococcus* sp have been found everywhere where suitable analytical techniques have been employed to detect them, such as in coastal and oligotrophic North Atlantic waters, in the Mediterranean, Indian and Pacific Oceans [12–14], and strains are available from many locations [15]. Analysis of seawater from many worldwide locations has shown that the genus *Ostreococcus* is common at depths of 0–120m below the surface where it contributes to varying extents to primary production [12,16]. Different strains show specific adaptations to different environments (depth, light intensity, [17]) and the complete genome sequence analysis of 2 strains representative of such niches has been completed [9,18], further strengthening the choice of this species as a model for host-virus interactions.

The presence of viruses in a coastal algal bloom of *O. tauri* has previously been observed [19], but these were not characterized further. We describe the isolation, growth in culture, life-cycle and gene content of *Ostreococcus tauri* virus 5 (OtV5), a DNA virus of OTH95, the best-characterized strain of *O. tauri*. We show by phylogenetic tree reconstruction with related viral taxa that OtV5 is a phycodnavirus, thus representing a commonly found group of marine virioplankton. Phycodnaviruses are nuclear-cytoplasmic large double-stranded DNA viruses (NCLDV) whose ancestors probably predate the separation of the eukaryotic kingdoms, with viruses that attack for example reptiles and mammals (such as Poxviruses) [20,21]. Genomes of algal phycodnaviruses identified

<figure><img src="image_2.png" /></figure>
to date are 150 to 560kb in size [22] (http://www.ncbi.nlm.nih.gov/ICTVdb/ICTv/fs_phyco.htm). Viruses attacking the Prasino-phyte *Micromonas pusilla* (previously known as *Chromulina pusilla* [23]), a common very small (1–3 µm) phytoplankter in coastal waters that often dominates in the “ultraplankton” [24], were first observed about 20 years ago [25], and the growth cycle of such a virus was subsequently observed in culture [26]. Since then, several authors have studied *Micromonas* viruses [27–32]. Important progress in studying the diversity of phycodnaviruses was enabled by the design of degenerate PCR primers using a conserved region of the polymerase gene of *Micromonas* viruses [33], since this sequence now serves as a phylogenetic landmark for both metagenomic data and newly characterized NCLDV sequences. Despite the global ecological importance of marine viruses suggested by metagenomic data [34] characterisation of phycodnavirus genomes has been limited until now to a few organisms including marine brown algae, *Ectocarpus siliculosus* [35], *Feldmania irregularis* (FirrV) [36] and *Emiliania huxleyi* [37,38], and to certain freshwater *Chorella* species, [39–41]. To our knowledge, we present the first whole genome analysis of a virus affecting marine green algae.

## Results

### Isolation and culture of OtV5

The presence or absence of viruses in water samples was assayed by filtering the water to eliminate bacterial cells and protists, then mixing the filtrate with cultures of the host, *O. tauri*. Lysis of the host cells was usually observed about 1 week after inoculation (Figure 1A). By mixing a sufficient number of host cells (typically about $10^8$ cells) in culture medium with 1 or 10 ml of seawater before plating in the soft gel, individual plaques could be visualized in a growing green lawn on the plate (Figure 1B, see the materials and methods [M&M] section).

In a preliminary study carried out in 2006 (56 samples tested), viruses specific to the strain OTH95 (i.e. giving lysis plaques on this strain) were detected in 19 out of 37 samples from local saltwater lagoons while 2 out of 19 coastal seawater samples tested positive by lysis of OTH95.

OtV5 was chosen because it reliably gave a relatively rapid lysis of host cells (2–3 days to clearing of the lysate, depending on the inoculum density, and because it produced a relatively high titre of viable particles in suspensions (typically $10^8$ plaque forming units [pfu] in a cleared lysate), OtV5 was purified twice by single-plaque isolation before further studies.

### Host specificity of OtV5

Using a microtitre dish assay (M&M) the host range specificity of OtV5 was evaluated on a series of different picoplankton

species, mostly isolated from the N.W. Mediterranean sea (Table 1). Neither more distantly related species (*Synechococcus* sp., *Tetraselmis* sp. *Dunaliella salina* *Nannochloris* sp., or *Pycnococcus* sp.,) nor more closely related *Mamiellales* species (*Bathycoccus* sp., or *Micromonas* sp.) were affected. More unexpectedly, among the 18 independently and locally isolated *Ostreococcus* spp. tested, none was affected, despite the fact that at least 9 of these strains had identical 18S rDNA sequences, the same as the host species OTH95 [17,42]. OtV5 thus shows strict host strain specificity.

### Life-cycle of OtV5

Routinely, inoculation of a late exponentially growing or stationary phase cultures ($2 \times 10^7–8 \times 10^7$ cells/ml) of host OTH95 cells with a fresh 0.2 µm-filtered OtV5 viral lysate ($>10^8$ pfu/ml) cleared the host suspension within 48h. In order to characterize the kinetics of viral growth and development of OtV5 particles, this process was monitored by flow cytometry, electron microscopy, plate assays, and pulsed field gel electrophoresis (PFGE) of DNA from whole cells embedded in agarose at different stages of the infection.

Firstly, we attempted to measure burst size and latent period using a moi (multiplicity of infection) of one. Cultures of exponentially growing host cells ($6.4 \times 10^6 \pm 0.5 \times 10^6$ cells/ml) were inoculated with a suspension of viral particles, estimated by flow cytometry to give $5.5 \times 10^6 \pm 0.9 \times 10^6$ viruses/ml. By flow cytometry, we estimated that 36% $\pm$12% had adsorbed to the cells within 5 minutes, and that non-adsorbed viruses then remained in suspension when mixed with the cells, and appeared not to adhere to cells for at least several hours (Figure 2A). By serial dilutions and plating out a sample of the suspension from 5 minutes post-inoculation, we estimated that $2.1 \times 10^6 \pm 0.5 \times 10^6$ cells/ml OTH95 host cells were actually infected. Since the number of viruses inoculated was $5.5 \times 10^6 \pm 0.9 \times 10^6$ and 36% of them had adsorbed, we would expect to obtain about $2.0 \times 10^6$ infected cells/ml. We thus conclude that all (or nearly all) of the adsorbed viruses were infectious. The number of viruses in the suspension started to increase appreciably after a latent period of 8 hours post-inoculation (hpi), and the curve leveled out at 12 to 16 hpi at $5.1 \times 10^7 \pm 0.9 \times 10^7$, giving a burst size of 25 viruses per cell. A further round of infection occurred after 16 hpi, leveling out after 20hpi (Figure 2B).

Secondly, we verified that the OtV5 genome could be prepared and visualized on PFGE gels (Figure 3B, band migrating around 180kb, see M&M). This band was RNAse-insensitive and disappeared after treatment with DNase (data not shown). Thirdly, a one liter culture of the host OTH95 strain in late exponential growth ($10^7$cells/ml) was inoculated with a fresh lysate of OtV5 at a moi of 1 pfu per host cell. We show (Figure 3A) that

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Assays for the presence of lytic viruses. **A:** Marine samples from Thau Lagoon (43° 24′N, 3° 36′E) incubated with the host strain OTH95 (isolated from the same place 14 years previously). Ticked boxes indicate lysis, crossed no lysis, and diamonds show mock-inoculated host cells. **B:** Viral lysis plaques in a plated gel of host cells.
doi:10.1371/journal.pone.0002250.g001</figcaption></figure>

<figure><img src="image_4.png" /></figure>
<table><caption>**Table 1.** Strain specificity of OtV5.</caption><thead><tr><th>Species</th><th>Class</th><th>Clade [17,42]</th><th>Strain<sup>a</sup></th><th>Lysis by OtV5</th><th>Reference</th></tr></thead><tbody><tr><td>Ostreococcus tauri</td><td>Prasinophyceae</td><td>C</td><td>RCC745</td><td>+</td><td>[10]</td></tr><tr><td>Synechococcus sp</td><td>Chroococcales</td><td></td><td></td><td>-</td><td>G. Vétion, pers. comm.</td></tr><tr><td>Dunaliella salina</td><td>Chlorophyceae</td><td></td><td></td><td>-</td><td>G. Vétion, pers. comm.</td></tr><tr><td>Tetraselmis sp</td><td>Prasinophyceae</td><td></td><td></td><td>-</td><td>G. Vétion, pers. comm.</td></tr><tr><td>Bathycoccus prasinos</td><td>Prasinophyceae</td><td>-</td><td>RCC1105</td><td>-</td><td>*</td></tr><tr><td>Micromonas pusilla</td><td>Prasinophyceae</td><td>A</td><td>RCC827</td><td>-</td><td>[15,90]</td></tr><tr><td>Nannochloris sp.</td><td>Trebouxiophyceae</td><td>unkown</td><td>BCC42000</td><td>-</td><td>*</td></tr><tr><td>Pycnococcus sp.</td><td>Prasinophyceae</td><td>unkown</td><td>BCC41000</td><td>-</td><td>*</td></tr><tr><td>Ostreococcus tauri</td><td>Prasinophyceae</td><td>C</td><td>RCC1108</td><td>-</td><td>*</td></tr><tr><td>Ostreococcus tauri</td><td>Prasinophyceae</td><td>C</td><td>RCC1118</td><td>-</td><td>*</td></tr><tr><td>Ostreococcus tauri</td><td>Prasinophyceae</td><td>C</td><td>RCC1117</td><td>-</td><td>*</td></tr><tr><td>Ostreococcus tauri</td><td>Prasinophyceae</td><td>C</td><td>RCC1114</td><td>-</td><td>*</td></tr><tr><td>Ostreococcus tauri</td><td>Prasinophyceae</td><td>C</td><td>RCC1558</td><td>-</td><td>*</td></tr><tr><td>Ostreococcus tauri</td><td>Prasinophyceae</td><td>C</td><td>RCC1559</td><td>-</td><td>*</td></tr><tr><td>Ostreococcus tauri</td><td>Prasinophyceae</td><td>C</td><td>RCC1112</td><td>-</td><td>*</td></tr><tr><td>Ostreococcus tauri</td><td>Prasinophyceae</td><td>C</td><td>RCC1111</td><td>-</td><td>*</td></tr><tr><td>Ostreococcus tauri</td><td>Prasinophyceae</td><td>C</td><td>RCC1115</td><td>-</td><td>*</td></tr><tr><td>Ostreococcus sp.</td><td>Prasinophyceae</td><td>D</td><td>RCC1107</td><td>-</td><td>*</td></tr><tr><td>Ostreococcus sp.</td><td>Prasinophyceae</td><td>D</td><td>RCC1122</td><td>-</td><td>*</td></tr><tr><td>Ostreococcus sp.</td><td>Prasinophyceae</td><td>D</td><td>RCC1119</td><td>-</td><td>*</td></tr><tr><td>Ostreococcus sp.</td><td>Prasinophyceae</td><td>D</td><td>RCC1120</td><td>-</td><td>*</td></tr><tr><td>Ostreococcus sp.</td><td>Prasinophyceae</td><td>D</td><td>RCC1121</td><td>-</td><td>*</td></tr><tr><td>Ostreococcus sp.</td><td>Prasinophyceae</td><td>D</td><td>BCC15000</td><td>-</td><td>*</td></tr><tr><td>Ostreococcus sp.</td><td>Prasinophyceae</td><td>unkown</td><td>RCC1123</td><td>-</td><td>*</td></tr><tr><td>Ostreococcus sp.</td><td>Prasinophyceae</td><td>unkown</td><td>RCC1561</td><td>-</td><td>*</td></tr><tr><td>Ostreococcus sp.</td><td>Prasinophyceae</td><td>unkown</td><td>RCC1110</td><td>-</td><td>*</td></tr></tbody></table>

after inoculation the culture continued to grow for 8h, albeit more slowly than the mock-inoculated control flask, attaining a plateau around 1.5×10⁷ cells/ml for 10–12h, before a progressive decrease in the number of autofluorescing cells over the next 12h period. PFGE analysis at frequent intervals during the time course permitted visualization of the genomes of both the host and the virus. At early stages of the infection, the host chromosomes can be seen, but no viral genomic DNA is visible (Figure 3C early stages 1 & 2). Full-length viral genomes appear about 2 hpi (hours post-inoculation, see Figure 3C. mid-late stage 3), as seen by the ~180kb band migrating close to the smallest chromosome of *O. tauri*, showing that the viral DNA replicates rapidly in host cells in this period. After 8h the overall number of copies of viral genomes visible increases dramatically as might be expected since newly produced viruses are released into the medium at this stage (Figure 2B); a large proportion of them probably adhere quickly to cells in the culture (see below, Figure 4A), and would thus be visible by PFGE after harvesting the host cells by centrifugation. After 30 hpi numbers of autofluorescing host cells decrease considerably. Remarkably, the host cell chromosomes remain intact throughout this entire period. Thereafter the diminishing quantities of host chromosomal DNAs reflect the decrease in the number of cells showing a level of autofluorescence typical of healthy cells.

Electron microscopy

Compared to the size of OTH95 cells (about 1 µm), the icosahedral particles of OtV5 were relatively large (122±9)nm, complete particles containing a central region of unevenly electron-dense material (see Figure 4E, for example). Using a high moi (2000 viruses per host cell in suspension), we observed that many viruses could adsorb to a single cell within 5 minutes of inoculation (Figure 4A). It is not clear whether the viruses first adhere to an electron-transparent layer surrounding the cell membrane, or whether they form a “bridge” directly with the host cell that lies outside the thin section. About 20% of sections from healthy cells (48/218 cells) appeared to have viruses adsorbed, whereas no viruses were visible on the sections of the remaining host cells, despite the very high moi used (plate 4A). We did not see any “tail” on the virus particle, but OtV5 could be seen attaching by a vertex to the external membrane of OTH95 by a bridge between the cytoplasm and the virus. Using an moi of 1, apparently complete virus particles appeared in the cytoplasm at 6 hpi (Figure 4 B). Viral particles were localized to a limited part of the cytoplasm. At 14 hpi 18% (42/229 cells) of mid-cell TEM sections showed particles. No viruses were found in the nucleus, the chloroplast or mitochondria. They usually occupied a space near the inner face of the nucleus, between it and the mitochondria or the chloroplast.

<figure><img src="image_5.png" /></figure>

<footnote>* strains isolated in the course of this work.</footnote>

<footnote>ªRCC-Roscoff Culture Collection, BCC-Banyuls collection, cryopreserved.</footnote>

<footnote>doi:10.1371/journal.pone.0002250.t001</footnote>
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Latent period and burst size of *OtV5*. The numbers of viruses in suspension starts to increase 8 to 9 hours post-inoculation (panel A) and rises to a first plateau at 12 to 16 hpi (panel B), giving a burst size of 25 after correction for the number of infected cells (see "Results"). Thereafter, a further cycle of infection can be seen by the 2nd increase in *OtV5* at 16-20hpi.
doi:10.1371/journal.pone.0002250.g002</figcaption></figure>

## Characterisation of the viral genome

Nucleic acid was isolated using standard techniques for manipulation of high M.W. DNA (see M&M). Briefly, a 2l culture of a fresh lysate was concentrated by ultrafiltration, and the concentrate incorporated into an agarose plug, and run out on a PFGE preparative gel before recovery of the viral genomic DNA band. This band was RNAse-insensitive (data not shown) and was used for preparation of a shotgun library. The genome is 186,234 base pairs long (Figure 5), in agreement with the size expected from migration of dsDNA by PFGE (Figure 3B and C). Its GC content (45%) is lower than that of the host (58%, [9], see Table 2). Terminal inverted repeats (TIR) of 1695 bp are predicted, their boundaries being defined by contigs spanning both the IR and flanking genomic DNA which were clearly identified in the assembly of the genome sequence. In addition, we showed that the left or right TIR could be amplified independently by PCR using primers in the adjacent single copy DNA (data not shown).

## Genome annotation

Seven hundred and thirty six ORFs (open reading frames) ≥65 amino acids (aa) long were considered, but 468 were discarded either because they did not possess an N-terminal methionine or because they appeared within other ORFs. Two hundred and sixty eight CDSs (coding sequences, defined here as ORFs flanked by start and stop codons) remained, thus giving a similar gene density to that of PCBV-1 (1.11 CDS/kb for PCBV1 and 1.45 for

*OtV5*). Eighteen per cent of predicted CDSs overlap, but this is not unexpected because some viruses are known to use polycistronic mRNAs [43]. Conceptual translations of a high proportion (57.5%) of the predicted CDSs showed no similarity to proteins with known functions using BLASTP alignment with public databases [44] (Figure 6). Genes of different functional categories and the transcriptional direction of these predicted CDSs are apparently distributed randomly on both DNA strands through the genome. Proteome analysis using Pfam identified many possible functional domains (Table S1). Five closely clustered tRNA-encoding sequences were found, including 2 with introns (see Genbank accession number EU304328).

Many enzymes potentially involved in the replication of viral DNA are predicted (Table 3, and EU304328). The predicted DNA polymerase, which is used as a reference for classification of large dsDNA viruses, shows greatest similarity to those of *Micromonas* viruses. A phylogenetic tree was made by comparing a conserved part of this gene with homologous sequences from other nuclear cytoplasmic large DNA viruses (NCLDV, [45]), with emphasis on algal DNA viruses (Figure 7). These observations, together with the presence of 31 genes whose predicted functions are conserved in other algal viruses (see below and GenBank accession number EU304328), including "core genes" of NCLDV [45], and our observations of the DNA genome size and particle morphology, support the inclusion of this virus in the phycodnavirus group [46,47].

Notably, several host-like DNA sequences (HLS) were identified in the viral genome. Six potential CDSs whose conceptual translations show significant similarity to *Ostreococcus spp.* genes were identified by their high BLASTP scores. Phylogenetic analysis of these 6 candidates (data not shown) grouped 4 ORF with other viral proteins, thus excluding the likelihood of recent horizontal transfer, and 2 other ORF together with eukaryotic sequences. However, the viral and host sequences are too divergent to make definitive conclusions about their evolutionary origins. Additional viral and host sequences from a range related taxa are needed to clarify this point. The largest of the HLS (253 aa), for which we present a phylogenetic comparison (Figure 8) encodes a potentially complete proline dehydrogenase gene.

## Discussion

Viruses infecting *O. tauri* could be found frequently in water samples collected from marine lagoons and open sea in the N. W. Mediterranean. They were found about five times more frequently in lagoons than in the sea (51% of samples vs. 11% of samples, respectively), perhaps reflecting the preferred habitat of its clade C [17,42] host. Detailed ecological observations and analysis of their phylogenies will be published elsewhere. Here we focus our attention on one virus, *OtV5*, isolated from the Bages lagoon on 24th January 2006, which lysed host cultures reproducibly and rapidly (within 2 days, Figure 2 and 3). We tested the strain specificity of *OtV5* on diverse picoeukaryotic strains of algae that were isolated and cloned mostly from local seawater samples between January 2006 and March 2007. Firstly, using 5 different genera of Prasinophyceae or 2 different clades of *Ostreococcus spp.*, no strains were susceptible. Secondly, none of a set of 9 independent and genetically polymorphic clade C strains was attacked by *OtV5*, which thus shows very narrow strain specificity. While further studies with larger sets of hosts and strains are required but are beyond the scope of this study, we note that this situation may resemble the gene-for-gene resistance arms race observed in higher plant-pathogen interactions [48,49], where natural selection drives diversity in plant disease resistance genes.

<figure><img src="image_7.png" /></figure>
<figure><img src="image_8.png" /></figure>

**Figure 3.** Time course of host cell lysis, karyotype, and production of OtV5 DNA during infection. **A:** host cell density in a control culture (open symbols) and a culture inoculated with OtV5 (filled symbols). **B:** Pulsed Field Gel Electrophoresis of samples from control DNAs and **C:** following the time course shown in A, showing *O. tauri* chromosomes and viral DNA. M–marker concatenated bacteriophage λ DNA ladder, C–control non-inoculated cultures, V–isolated OtV5 DNA. Stage 1–a culture 30sec pre- and 30sec post-inoculation, stage 2–20min and 40min post-inoculation, stages 3 and 4–at 60min and 1h intervals thereafter (images from 2 gels were pasted together to accommodate the number of gel tracks necessary). The smallest host chromosome migrates at the position of the 145kb marker band (dot), and the virus-specific DNA band migrates around 180kb (diamond). **D:** Flow cytometric enumeration of *O. tauri* analysis using autofluorescence of chlorophyll a as a marker. The window "O.T." contains cells showing autofluorescence typical of healthy growing cells. At 19 hours after infection, about 20% of the cells have died as shown by reduced fluorescence (ordinate axis). At 27h about 90% have died and about 1% remains at 43h.
doi:10.1371/journal.pone.0002250.g003

<figure><img src="image_9.png" /></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 4.** Electronic microscopy of infected *Ostreococcus tauri* cells. A, B, C, D, the bar represents 500nm; E, F, bar 50nm. Example virus particles are shown with arrows. Chl-chloroplast; Cyt-cytoplasm, n-nucleus, m-mitochondrion, sg-starch grain. A-at high moi many viruses can adsorb to a single cell. In B, C, D, E, & F the moi was one. B & C show viruses accumulating in the cytoplasm before cell lysis occurs. D-virus particles clumped together around a lysed cell, E (30min after infection)-full virus particle attached to plasmalemma, F-empty particle.
doi:10.1371/journal.pone.0002250.g004</figcaption></figure>

Relatively little is known about the selective forces driving the diversity and speciation of microbes in the marine environment, but viruses appear to play a key role in these processes that furnish the base of the global food web [5,7,8,32,50–55].

The extremely small *O. tauri* cell is crowded with organelles, which are often observed to deform the apposed plasmalemma, and the cytoplasmic compartment represents a maximum of about 30% of the *O. tauri* cell volume [56]. If all of this were packed with viruses, about 100 particles might be accommodated. However this volume also accommodates vesicles and hundreds of ribosomes, virus particles being localized to a limited part of the cytoplasm. It is thus not surprising that the mean burst size is only 25 particles per cell. EM pictures showed a maximum of 15

particles in a 50nm thick section of a single cell. These data are comparable to the burst size of 70–100 particles per cell observed for the related host/virus system *Micromonas pusilla/MpV* [26], considering the host cell volume of *Micromonas* is about 4 times more voluminous than *Ostreococcus* (the diameter of *O. tauri* in our culture conditions is about 1 µm [Figure 4] and *Micromonas* is 1.6 µm [57]).

Using a high moi, EM also showed that many virus particles could adsorb to a single cell (Figure 4A) but that high numbers of viruses apparently adsorbed to only a proportion of the cells (about 20%). This observation was confirmed by flow cytometry at lower moi (at a moi of one, about 35% of viruses were adsorbed). The remaining majority of particles then remained in suspension

<figure><img src="image_11.png" /><figcaption>**Figure 5.** Features of the OtV5 virus genome. Predicted coding sequences are represented by colored rectangles (white—the reading frame is on the opposite DNA strand). The bar and filled triangles indicate map lengths in kb.
doi:10.1371/journal.pone.0002250.g005</figcaption></figure>

<figure><img src="image_12.png" /></figure>
<table><caption>Table 2. Comparison of the general characteristics of the genomes of OtV5 with *Ostreococcus tauri* strain OTH95 and other algal viral genomes.</caption><thead><tr><th></th><th>OtV5</th><th><sup>2</sup>O. tauri [9]</th><th>PBCV-1 [91]</th><th>Ehv-86 [37]</th></tr></thead><tbody><tr><td>Genome size (bp)</td><td>186,234</td><td>12,560,000</td><td>330,743</td><td>407,339</td></tr><tr><td>% G+C</td><td>45</td><td>58</td><td>40</td><td>40</td></tr><tr><td>Predicted CDSs (≥65 aa)</td><td>268</td><td>8,166<sup>b</sup></td><td>367</td><td>472</td></tr><tr><td>Mean gene size (bp/gene)</td><td>692</td><td>1300</td><td>899</td><td>863</td></tr><tr><td>% known genes</td><td>22</td><td>77</td><td>50</td><td>21</td></tr><tr><td>% hypothetical genes</td><td>20</td><td></td><td></td><td></td></tr><tr><td>% genes with no hit<sup>a</sup></td><td>57</td><td>23</td><td>50</td><td>79</td></tr></tbody></table>

<figure><img src="image_13.png" /><figcaption>**Figure 6. Taxon distribution of best BLASTP hits.** A: using OtV5 viral genes. B: using host *Ostreococcus tauri* genes.
doi:10.1371/journal.pone.0002250.g006</figcaption></figure>

following inoculation despite the presence of an actively growing population of host cells (Figure 2A), suggesting that such particles lacked a host cell attachment function, and/or that the remaining cells were resistant to infection. The latter hypothesis, however, seems unlikely given the clonal nature of the host strain. The above information, together with the actual number of infected host cells estimated by pfu (see the “Results” section), suggests that viral attachment may be a limiting step in the infection. At high moi 5 min after inoculation, a clear layer about 50nm thick separated the particles from the cytoplasm (Figure 4A). Since apparently few of the viruses have formed a “bridge” with the cytoplasm at this stage (the EM slice is 50nm thick), we favor the hypothesis that the cell is protected by a cellular envelope, never previously observed in this algal genus, that is invisible in our EM preparations. The biochemical nature of this putative envelope remains to be determined.

We could not however use EM images from high moi to determine when the particles “injected” their electron-dense contents, because all of the cells with a large numbers of viruses attached deteriorated rapidly (<20min, not shown) and disappeared from the culture, perhaps due to multiple perforations of the membrane or to degradation caused by injection of unknown virion proteins, involved in virulence, as apparently may occur in related NCLDV [58]. Host cells with “empty” virus particles attached (Figure 4F) were rarely seen at any stage of the infection, suggesting that they might usually detach from the cells after injecting their contents.

The average size of predicted CDSs in the linear double-stranded 186,234 bp DNA genome is similar to other phycodnaviruses. One hundred and fourteen CDSs (42.5%) encoded

predicted proteins with similarities to other hypothetical proteins in databases, but only 60 of these (22.4%) showed similarity with proteins to which a function has been attributed. This last category may, however, be overestimated as some of these may not be bona fide CDSs. Fifty four (20.1%) CDSs matched other hypothetical CDSs of unknown function, of which 28 (10.5%) came from other viruses, and 26 (9.7%) from elsewhere, and 154 (57.5%) of genes were unknown. Finally, predicted products of 154 genes (57.5%) correspond to no known protein, emphasizing the extremely diverse range of uncharacterized biological functions potentially represented in the phycodnaviruses.

In common with other algal viruses (Table S2), numerous enzymes involved in DNA replication and gene expression are predicted. This is an ancient common feature of the super-group of NCLDV, including phycodnaviruses, which retain a set of predicted functions that presumably assure their life-cycles in the host cytoplasm. Like PBVC-1, the OtV5 genome shows terminal inverted repeat sequences, and its ends are thus likely to form hairpin loops [59–61]. Our phylogenetic reconstruction shows that OtV5 clearly clusters with phycodnaviruses infecting the prasinophyte *Micromonas pusilla* [33]. These latter viruses form the Prasinovirus clade within Phycodnaviridae, and Prasinovirus should then include OtV5, from a prasinophyte host. Our tree is different from that of Schroeder et al. (2002), because we use different methods for the analysis, giving higher node support values (posterior probabilities). However, in both trees, Phycodnaviridae (Phaeovirus+Coccolithovirus+Prasinovirus+Chlorovirus+Prymnesiovirus) do not form a monophyletic group, with Herpesviridae and *Heterosigma akashiwo* virus nested inside (Figure 7). NCLDV, which attack hosts as diverse as man, reptiles, fishes, invertebrates and algae, are thought to have evolved from a common ancestor, with a core set of common genes [20,45]. Like some of its poxviral cousins [62], OtV5 probably accumulates embedded in a cytoplasmic region rich in a viral-A type inclusion protein, since a CDS whose product shows similarity to this protein was identified (Table 3 and Genbank accession number EU304328).

We found no evidence of site-specific endonucleases which are present or predicted in some other previously characterized phycodnaviruses [35, 37, 63–65], although similarities to 3 different methylation/modification enzymes are present. DNA degradation in host *Chlorella* cells by PBCV-1 (*Paramecium bursaria Chlorella* virus 1) is an early viral function, occurring within 5 minutes of host cell penetration [65], in stark contrast to OtV5, where host chromosomes remain intact throughout the viral life-cycle (Figure 3C). In this respect, OtV5 resembles MT325, which also lacks predicted endonucleases [65]. However, in contrast to MT325, where complete viral genomes are present immediately after inoculation, OtV5 genomic DNA could not be visualized by PFGE within the first 2h after inoculation. Since flow cytometrical observations and pfu measurements show that about 35% of the cells have adhered viruses at this stage, we conclude that this level is below the level of detection by PFGE, and that viral genomes become visible after their replication inside host cells.

OtV5 encodes 3 predicted adenine methyltransferases that are possibly involved in protection of the host and/or the viral DNA, for example from the kind of programmed cell death response observed in other green algae [66], or from attack by a third party invader. Like its larger mimivirus relative [64], OtV5 encodes an asparagine synthase (AS) gene, and several viral tRNAs are predicted, including Asn-tRNA. In green plants, AS is a cornerstone of light-regulated protein synthesis, controlling the ratio of carbon and nitrogen anabolites available for protein or sugar syntheses diurnally [67]. OtV5 carries a CDS with similarity to FtsH metalloendopeptidase that is absent in other characterized eukaryotic algal genomes.

<figure><img src="image_14.png" /></figure>

<footnote><sup>a</sup>BLASTP E-value threshold $10^{-5}$</footnote>

<footnote><sup>b</sup>includes all predicted CDSs</footnote>

<footnote>doi:10.1371/journal.pone.0002250.t002</footnote>
<table><caption>**Table 3.** Similarities to OtV5 predicted proteins</caption><thead><tr><th>CDS#</th><th>Gene names</th><th>BBH</th><th>GG</th><th>CDS#</th><th>Gene names</th><th>BBH</th><th>GG</th></tr></thead><tbody><tr><td colspan="4"><strong>DNA, RNA, replication, recombination and repair</strong></td><td colspan="4"><strong>Lipid-Fatty acid metabolism</strong></td></tr><tr><td>240f</td><td>DNA polymerase</td><td>V</td><td>+</td><td>040f</td><td>Oxo-acyl (acyl-carrier protein) dehydrogenase</td><td>B</td><td>-</td></tr><tr><td>205r</td><td>ATP dependant DNA ligase</td><td>O</td><td>-</td><td>136r</td><td>Patatin-like phospholipase</td><td>V</td><td>+</td></tr><tr><td>199f</td><td>DNA topoisomerase</td><td>B</td><td>+</td><td colspan="4"><strong>Protein synthesis, modification and degradation</strong></td></tr><tr><td>244r</td><td>DNA topoisomerase</td><td>V</td><td>+</td><td>092f</td><td>N-myristoyltransferase</td><td>M</td><td>-</td></tr><tr><td>115f</td><td>PCNA</td><td>V</td><td>+</td><td>126r</td><td>Prolyl 4-hydroxylase</td><td>V</td><td>+</td></tr><tr><td>067f</td><td>Helicase</td><td>V</td><td>+</td><td>152r</td><td>33 kDa in vitro translation peptide</td><td>V</td><td>-</td></tr><tr><td>005r</td><td>Exonuclease</td><td>B</td><td>+</td><td>171r</td><td>33 kDa in vitro translation peptide</td><td>V</td><td>-</td></tr><tr><td>146f</td><td>Exonuclease</td><td>V</td><td>+</td><td>013r</td><td>FtsH2 (3-4-24) metalloendopeptidase</td><td>B</td><td>-</td></tr><tr><td>098r</td><td>ATPase (DNA packaging)</td><td>V</td><td>+</td><td>042f</td><td>Aminotransferase family protein</td><td>M</td><td>-</td></tr><tr><td>194f</td><td>Rnase H</td><td>B</td><td>+</td><td>028r</td><td>Asparagine synthase</td><td>B</td><td>-</td></tr><tr><td colspan="4"><strong>DNA methylation and site-specific endonucleases</strong></td><td>041r</td><td>NAD-dependant epimerase dehydratase</td><td>B</td><td>-</td></tr><tr><td>097f</td><td>Adenine DNA methyltransferase</td><td>B</td><td>-</td><td>203f</td><td>Aspartyl/Asparaginyl beta-hydrolase</td><td>B</td><td>-</td></tr><tr><td>144f</td><td>Adenine DNA methyltransferase</td><td>V</td><td>-</td><td>226f</td><td>Aspartyl/Asparaginyl beta-hydrolase</td><td>B</td><td>-</td></tr><tr><td>124f</td><td>M6A-gamma-methyltransferase</td><td>V</td><td>-</td><td>094f</td><td>Zinc metallo protease</td><td>V</td><td>+</td></tr><tr><td colspan="4"><strong>Sugar manipulation enzymes</strong></td><td colspan="4"><strong>Signaling</strong></td></tr><tr><td>011f</td><td>GDP-D-mannose dehydratase</td><td>M</td><td>-</td><td>030f</td><td>Potassium channel protein</td><td>B</td><td>+</td></tr><tr><td>012r</td><td>Glycosyl transferase</td><td>B</td><td>+</td><td>237r</td><td>Rhodanese-like domain protein</td><td>B</td><td>+</td></tr><tr><td>020r</td><td>Glycosyl transferase</td><td>V</td><td>+</td><td>168r</td><td>Serine/Threonine protein kinase</td><td>B</td><td>+</td></tr><tr><td>035f</td><td>Glycosyl transferase</td><td>B</td><td>+</td><td colspan="4"><strong>Capsid proteins</strong></td></tr><tr><td>037f</td><td>Acetolactate synthase (2-2-1-6)</td><td>B</td><td>-</td><td>068f</td><td>Capsid protein 1</td><td>V</td><td>+</td></tr><tr><td>197f</td><td>6-Phosphofructokinase</td><td>B</td><td>-</td><td>076f</td><td>Capsid protein 2</td><td>V</td><td>+</td></tr><tr><td>043f</td><td>dTDP-D-glucose 4,6-dehydratase</td><td>M</td><td>-</td><td>081f</td><td>Capsid protein 3</td><td>V</td><td>+</td></tr><tr><td>218f</td><td>dTDP-D-glucose 4,6-dehydratase</td><td>M</td><td>-</td><td>099f</td><td>Capsid protein 4</td><td>V</td><td>+</td></tr><tr><td colspan="4"><strong>Transcription</strong></td><td>100f</td><td>Capsid protein 5</td><td>V</td><td>+</td></tr><tr><td>177f</td><td>Transcription factor IIB</td><td>B</td><td>+</td><td>190f</td><td>Capsid protein 6</td><td>V</td><td>+</td></tr><tr><td>024r</td><td>Transcription factor IIS</td><td>V</td><td>+</td><td>239r</td><td>Capsid protein 7</td><td>V</td><td>+</td></tr><tr><td>149f</td><td>Transcription activator</td><td>V</td><td>+</td><td>243r</td><td>Capsid protein 8</td><td>V</td><td>+</td></tr><tr><td>161f</td><td>Tata-box Binding Protein (TBP)</td><td>V</td><td>-</td><td colspan="4"><strong>Miscellaneous</strong></td></tr><tr><td>145f</td><td>RNase III</td><td>V</td><td>+</td><td>039f</td><td>3-Dehydroquinate synthase</td><td>B</td><td>-</td></tr><tr><td colspan="4"><strong>Nucleotide metabolism</strong></td><td>188r</td><td>ATP/GTP binding site motif A AGB-1</td><td>V</td><td>+</td></tr><tr><td>153r</td><td>Ribonucleotide reductase large subunit</td><td>B</td><td>+</td><td>014f</td><td>Methyltransferase FkbM</td><td>B</td><td>-</td></tr><tr><td>174r</td><td>Ribonucleotide reductase small subunit</td><td>M</td><td>+</td><td>135r</td><td>Viral A-type inclusion protein</td><td>M</td><td>-</td></tr><tr><td>227r</td><td>Thymidine kinase</td><td>G</td><td>+</td><td>059f</td><td>Flavin-dependant thymidylate synthase</td><td>M</td><td>-</td></tr><tr><td>231f</td><td>dUTP pyrophosphatase</td><td>V</td><td>+</td><td>180f</td><td>mRNA capping enzyme</td><td>V</td><td>-</td></tr></tbody></table>

BBH-best BLASTP hit similar to: B-Bacteria, G-Viridiplantae (green plant lineage), M-Metazoa, V-Virus. GG-predicted gene present in other viruses of Chlorophyta (Green Gene)

doi:10.1371/journal.pone.0002250.t003

Perhaps this enzyme performs a similar function to its homologue found in photosynthetic bacteria and plants, where it is involved in the repair maintenance of photosystem II [68], since the chloroplast remains intact throughout the viral life-cycle, and most likely functional as the starch grain continues to enlarge. However, all of the above hypotheses remain speculative and must await experimental investigation.

Interestingly, we found some indications for the capture of host DNA within the viral genome, since 6 predicted host-like sequences (HLS) present similarity to host genes (GenBank accession number EU304328). These are unlikely to be very recent transfers, because their GC content is clearly virus-like,

(45%) rather than host-like (58%). Nevertheless they may represent an important evolutionary possibility for transfer of genetic information between host cells, since the amino acid sequences gave very high BLASTP scores with *Ostreococcus* spp. proteins, and have thus not had time to degenerate extensively by mutations, deletions, or other rearrangements. Remarkably, one of these CDSs encodes an apparently complete proline dehydrogenase gene. This enzyme is important for proline catabolism in all cellular organisms, is known to play a role in stress response in plants and in some other algae [69], and is used as a source of nutrients by numerous bacterial symbionts and pathogens [70,71]; it could thus play an important role in some aspect of the host-

<figure><img src="image_15.png" /></figure>
<figure><img src="image_16.png" /><figcaption>**Figure 7. Phylogenetic tree showing OtV5.** Asterisks indicate that at least one complete virus genome is available within the group. A phylogenetic tree computed by Bayesian analysis (BI) and maximum likelihood (ML) analyses on the partial DNA polymerase gene DNA sequences of selected taxa of the family Phycodnaviridae and other large dsDNA viruses. Genbank accession numbers of reference taxa are given in parentheses. The position of OtV5 is indicated in large bold type. Numbers are posterior probabilities (BI) and bootstrap proportions (ML; values below 50 are indicated by “-”) reflecting clade support.
doi:10.1371/journal.pone.0002250.g007</figcaption></figure>

virus interaction. Phylogenetic analysis suggests that the OtV5 gene could originate from bacteria (Figure 8). The AT content of this CDS is 55% in OtV5 vs. 41% in *O. tauri*, a highly significantly difference ($p<10^{-16}$ [exact binomial test]). This difference is particularly striking on synonymous positions, suggesting that the functionality of this predicted protein was maintained by natural selection following its transfer to the virus, but a difference nevertheless remains on non-synonymous positions once the synonymous changes are removed (4-fold degenerate codons, $p<10^{-5}$). The latter may arise from selective pressures for reducing the cost of viral genome replication and/or requirements for viral transcription and translation [72]. We have not found this predicted gene in other viruses in databanks, but in plants proline levels increase in stress conditions such as pathogen attack, and are controlled by proline degradation, inviting speculation that this enzyme plays a role in regulating the host stress response [69]. The CDS of the other HLS are smaller, but they nevertheless support the notion that exchange of sequences between host and viral genomes may be an evolutionarily frequent event. However, weak similarities between genes make it difficult to prove this with phylogenetic analyses, because of the absence of information from sufficiently closely related organisms.

The very small size of the cells and genomes of the *Mamiellaceae*, together with their worldwide distribution and their ecological importance, promise that they will provide key models for interdisciplinary approaches to global ecology. Several different complete Prasinophyte genomes have already been sequenced or are in the process of being analyzed, including 3 *Ostreococcus* spp. [9,18], 1 *Bathyccoccus* sp. (unpublished data) and 2 *Micromonas* spp. (A. Worden, personal communication). Viruses play a crucial role in the regulation of phytoplankton populations [5,34,52,54,55,73], and genomic analyses should provide insight about the biological phenomena driving selection and diversity at the base of the food web.

## Materials and Methods

### Culture of host algal strains and virus OtV5

The host strain OTTH0595 was used for *Ostreococcus tauri* culture in all experiments, and cytometric observations made as previously described [10]. Cultures were grown in Keller-medium (Sigma-Aldrich, Saint-Quentin Fallavier, France) diluted in 0.22 µm filtered sea water (NaCl 36 g.L⁻¹) under continuous light (100 µmol photon m⁻²s⁻¹) at 20±1°C and under mild agitation. OtV5 was first identified as a single lysis plaque by inoculation of a sample of 1ml of 0.2 µm-filtered seawater taken from Bages Lagoon on 24/1/2006 into 20ml of a 1-week old host OTH95 culture plated at 1×10⁷ cells.ml⁻¹ in 0.15% agarose. After two rounds of single-plaque purifications, viruses were subsequently produced by infecting liquid cultures of *O. tauri* either in exponential phase (e.g. 1×10⁷ cells/ml) or in stationary phase (about 8×10⁷ cells/ml). For preparation of large quantities of viruses for genome sequencing, two litres of an *O. tauri* exponentially growing culture (approx. 5.10⁷ cells ml⁻¹) was inoculated with 200 ml of OtV5 lysate (approx. 5×10⁸ pfu ml⁻¹, ratio *O. tauri*/OtV5 around 1). Lysed cultures were passed sequentially through a 5 µm and 0.45 µm filters to remove large cellular debris. Virus filtrates were concentrated by ultrafiltration with a 50K MW size cut-off unit (Millipore, Amicon Ultra) to a final volume of 6 ml. Virus concentrates were embedded in agarose for PFGE (see below). The density of virus particles in suspension was also estimated by flow cytometry [74]

### Plaque visualization assay

To facilitate the isolation and purification of viruses it was desirable to develop a plating technique in order to visualize individual lysis plaques. Since we were unsuccessful in growing the host cells on the surface of agar plates, we used a plate-gel technique for plating *O. tauri* embedded in a soft gel (N. West, personal communication). Lysis “plaques” could then be visual-

<figure><img src="image_17.png" /></figure>
<figure><img src="image_18.png" /><figcaption>**Figure 8. Phylogenetic tree of OtV5 proline dehydrogenase.** A phylogenetic tree computed by Bayesian analysis (BI) and maximum likelihood (ML) analyses on partial Proline dehydrogenase and Proline oxidase aa sequences. Genbank accession numbers of taxa are given in parentheses. The position of OtV5 and *O. tauri* are indicated in bold type. Numbers are posterior probabilities (BI) and bootstrap proportions (ML).
doi:10.1371/journal.pone.0002250.g008</figcaption></figure>

ized as cleared circular regions extending through the whole depth of the gel. 3ml of an ~3×10⁷ cells/ml growing culture of *O. tauri* were mixed with the viral inoculum in 0.1–10ml seawater and the volume brought to 20ml by adding K-medium and kept at 20° on the bench for ≥30min. Separately, a suspension of 1.5% agarose (Euromedex type D-5, DNA grade) was dissolved in distilled water by heating to 100°C, and then held at 65° in a water bath in tubes in 2.2ml aliquots before use. For plating, the 20ml of cell suspension was added to 2.2ml of hot agarose and mixed rapidly before pouring into a 9cm diameter Petri dish. The plates were cultured (continuous light 100 µmol photon m⁻²s⁻¹, at 20±1°C) inside a transparent plastic box to maintain humidity for up to 1 to 2 weeks. Plaques could usually be visualised 3 to 8 days after inoculation as cleared areas on a green background. The titre of viable viruses present in the suspensions used for inoculations (and moi calculations) were measured routinely by counting lysis plaques in 2 parallel independent dilution series' on such OTH95 host plates

**Flow cytometry.** Analyses were performed on a FACScan flow cytometer (Becton Dickinson, San Jose, CA, USA) equipped with an air-cooled argon laser (488 nm, 15 mW). *Ostreococcus tauri* cells were characterized according to their right-angle scatter and their red fluorescence emission due to the chla pigment [75]. OtV5 viruses were determined both by their right-angle scatter and their fluorescence after SYBR green I staining. Viruses were stained according to [74].

**Electron microscopy.** For transmission electron microscopy (TEM), the cells were prepared according to Chretiennot-Dinet, 1995. They were fixed with glutaraldehyde in their culture

medium at a final concentration of 1% and were concentrated by centrifugation for 15 min at 3000g. The centrifuged material was then mixed rapidly with 40 µl of 1% liquid agarose (agarose type II, Sigma) at 37°C using a microrepetor (SMI, Emerville, CA, USA). Once the agarose in the disposable glass micropipette had solidified, a “noodle” containing the cells was obtained and treated as a piece of agar, fixed for 2 h at 4°C in 2.5% glutaraldehyde with one volume of 0.4M cacodylate buffer and two volumes of culture medium (A solution). The noodle was then washed 3 times for 30 min in a buffer medium made of one volume 0.4M cacodylate buffer and one volume culture medium (B solution). Postfixation was carried out in 1% OsO₄ in 0.2M cacodylate at 4°C for 1h. After two washes in 0.2M cacodylate, the noodle was cut in small pieces, dehydrated in a series of ethyl alcohol and embedded in Epon 812. Thin sections were stained with uranyl acetate and lead citrate before examination on a 7500 Hitachi transmission electron microscope.

**Pulsed field gel electrophoresis (PFGE).** Agarose embedded *O. tauri* cells with or without OtV5 were analysed by PFGE as described previously [76]. From 1.10⁷ to 2.10⁷ cells were loaded per lane in a 1% agarose gel (Type D-5, Euromedex France) and the electrophoresis parameters were 6 V cm⁻¹, 0.5× TBE, 120° switching angle, 14°C, 60 s switch time for 15 h then 90 s switch time for 9 h. The gel was then stained with ethidium bromide in order to visualise the chromosomal bands.

## OtV5 genome

Individual lots of about 1 µg of OtV5 nucleic acids were digested overnight at 37°C by DNase, or RNase (Sigma, France),

<figure><img src="image_19.png" /></figure>
or different restriction enzymes (HindIII, BamHI and SmaI; Promega France) and the digested products were analyzed by electrophoresis (1% agarose, 100V/cm, 20 min.).

Genomic DNA for sequence analysis was prepared by embedding viral particles in agarose strings at approximately 10¹¹ particles/ml and lysing of the particles by proteinase K. Agarose was then digested by gelase (Epicentre, Tebu) and DNA was precipitated by ethanol. Purified DNA was broken by sonication, DNA fragments ranging from 1 to 5 Kb were separated in an agarose gel before purification and end-filling. Blunt-end fragments were inserted into pBluescript IIKS (Stratagene) digested by EcoRV and dephosphorylated. Plasmid DNA from recombinant *E. coli* strains was extracted according to the TempliPhi method (GE Healthcare) and inserts were sequenced on both strands using universal forward and reverse M13 primers and the ET DYEnamic terminator kit (GE Healthcare). Sequences were obtained with MegaBace 1000 automated sequencers (GE Healthcare). Data were analyzed, and contigs were assembled by using Phred-Phrap [77] and Consed software packages (http://www.genome.washington.edu). To collate data and facilitate annotation 2861 reads (each ~600bp long) of DNA sequences derived from the shotgun library were assembled, giving 9.2-fold coverage. Remaining gaps (72) were filled by primer-directed sequencing on complete viral genome DNA using custom-made primers using a BigDye3 dye terminator sequencing kit (Applied Biosystems, Courtaboeuf, France) and analysed on an Applied Biosystems 3100 sequencer.

## Sequence annotation

The complete genomic sequence was annotated using the "Artemis" software [78,79]. Amino acid sequences corresponding to all predicted CDS were screened against the curated Pfam-A profiles [80] for functional motifs using the perl script Pfam_scan.pl (Wellcome Trust Sanger Institute, UK) using default settings. The e-value cut-off of 1 was conserved as profiles were mostly defined from distantly-related sequences and many functional motifs with relatively high e-values are confirmed by significant BLASTP alignments [81] with proteins having the same putative function.

## Sequence retrieval and alignment

Reference partial DNA-polymerase DNA sequences from Phycodnaviridae and phylogenetically related large dsDNA viruses were gathered from Genbank. We chose reference taxa among those previously used for comparing algal viruses [82], then added Mimivirus [83] and new Phycodnaviruses released in Genbank after this study. The homologous partial DNA polymerase sequence from OtV5 was aligned with these reference sequences, and the sequence of *Vaccinia virus* (Poxvirus) was used as outgroup following Lakshminarayan et al. (2006). Alignment was performed with ClustalX [84] and improved by eye with Se-Al v2.0a11 (Rambaut, A. 1996. Se-Al: Sequence Alignment Editor. Available at http://evolve.zoo.ox.ac.uk/). Ambiguous regions were deleted from the alignment before phylogenetic analysis, for a final length of 678 bp.

## Phylogenetic reconstruction

A phylogenetic tree including OtV5 was constructed by using the partial DNA polymerase sequence routinely amplified in other related viruses [33] by Bayesian analysis and Maximum likelihood. Bayesian analysis was done with MrBayes 3.1.2 [85]. The reconstruction used 4 chains of 10⁶ generations with an evolutionary model designed for coding sequences and considering each codon position as a distinct partition. Burnin value was set to 20% of the sampled trees (1% of the number of generations).

Maximum likelihood was performed with PhyML [86,87], using a General Time Reversible evolutionary model accounting for substitution rate heterogeneity and a proportion of invariable sites, chosen with ModelTest [88], and validated via a bootstrap procedure with 100 replicates. To screen the OtV5 genome for potential horizontal gene transfer candidates, we systematically screened OtV5 ORFs (>40 aa) for homology with the *O. tauri* genome using tblastn and retrieved 17 OtV5 ORFs with aa identity greater than 45% over more than 50 aa. We then searched for GenBank BBHs (Best Blast Hits, using blastp), and identified 6 ORFs that contain *O. tauri* sequences in the 50 BBHs. Phylogenetic reconstructions including HLS were performed from aa sequences with the same methods as above. The tree presented on Figure 8 was built from an alignment of length 165 aa, using a WAG+γ model chosen with ProtTest [89]

## Data deposition

The complete OtV5 DNA sequence has been deposited in GenBank, accession number EU304328. New algal strains have been deposited in the Roscoff Culture Collection [15], France.

## Supporting Information

**Table S1** Presence of Pfam motifs in putative OtV5 proteins. Pfam motifs were detected using the pfam_scan.pl script with default values. Motifs corresponding to the same CDS product are grouped and coloured alternatively red and blue. Amino acid sequences with no detectable motif are not presented. Columns show respectively the OtV5 gene name, start and end coordinates on the amino acid sequence, the Pfam accession number (linked to the corresponding Pfam web page), indication of complete (l) or partial (f) motifs, Bit value, e-value, Pfam ID and whether the motif is nested in another motif.

Found at: doi:10.1371/journal.pone.0002250.s001 (0.19 MB XLS)

**Table S2** Comparison of OtV5 with other algal viruses. OtV5: *O. tauri* virus strain 5; PBCV-1: *P. bursaria* chlorella virus strain 1; *P. bursaria* chlorella virus strain NY-2A; MT325: *P. bursaria* chlorella virus strain MT325; EhV-86: *E. huxleyi* virus strain 86. For each gene, numbers correspond to the number of copies of this gene family in the genome. Highlighting colors: yellow-common to characterized micro-algal large DNA viruses, orange-common to characterized phycodnaviruses of Chlorophyta, green-not present in other characterized phycodnaviruses.

Found at: doi:10.1371/journal.pone.0002250.s002 (0.03 MB XLS)

## Acknowledgments

We are very grateful to Nyree West for providing invaluable advice on plating out host cells, Gilles Vétion for providing certain algal strains, Michèle Laudié for technical help in DNA sequencing, and two anonymous reviewers of PLoS ONE for their helpful detailed comments. We thank Pascal Romans, Ingrid Obernosterer, Jean-Jacques Naudin and the "Néréis II" nautical team for their help in collecting seawater samples, the flow cytometry-imaging platform of CNRS/University of Paris6 for microscopy facilities, and Florence Le Gall for sending us the RCC827 strain.

## Author Contributions

Conceived and designed the experiments: HM NG ED. Performed the experiments: HM NG ED CF ME SE RC. Analyzed the data: GP HM NG ED CF ME RC YD. Contributed reagents/materials/analysis tools: LB. Wrote the paper: HM NG.

<figure><img src="image_20.png" /></figure>
## References

1. Suttle CA (2005) Viruses in the sea. Nature 437: 356-361.

2. Suttle CA (2007) Marine viruses—major players in the global ecosystem. Nat Rev Microbiol 5: 801–812.

3. Bratbak G, Wilson W, Heldal M (1996) Viral control of *Emiliania huxleyi* blooms? J Mar Systems 9: 75–81.

4. Wilson WH, Tarran GA, Schroeder D, Cox M, Oke J, et al. (2002) Isolation of viruses responsible for the demise of an *Emiliania huxleyi* bloom in the English Channel. J Mar Biol Assoc UK 82: 369–377.

5. Schroeder DC, Oke J, Hall M, Malin G, Wilson WH (2003) Virus succession observed during an *Emiliania huxleyi* bloom. Appl Environ Microbiol 69: 2484–2490.

6. Evans C, Archer SD, Jacquet S, Wilson WH (2003) Direct estimates of the contribution of viral lysis and microzooplankton grazing to the decline of a *Micromonas* spp. population. Aquatic Microb Ecol 30: 207–219.

7. Brussaard CP (2004) Viral control of phytoplankton populations—a review. J Eukaryot Microbiol 51: 125–138.

8. Weinbauer MG, Rassoulzadegan F (2004) Are viruses driving microbial diversification and diversity? Environ Microbiol 6: 1–11.

9. Derelle E, Ferraz C, Rombauts S, Rouze P, Worden AZ, et al. (2006) Genome analysis of the smallest free-living eukaryote *Ostreococcus tauri* unveils many unique features. Proc Natl Acad Sci USA 103: 11647–11652.

10. Courties C, Vaquer A, Troussellier M, Lautier J, Chretiennotdinet MJ, et al. (1994) Smallest Eukaryotic Organism. Nature 370: 255–255.

11. Chretiennot-Dinet MJ, Courties C, Vaquer A, Neveux J, Claustre H, et al. (1995) A New Marine Picoeucaryote-*Ostreococcus tauri* (Chlorophyta, Prasinophyceae). Phycologia 34: 285–292.

12. Worden AZ, Nolan JK, Palenik B (2004) Assessing the dynamics and ecology of marine picophytoplankton: The importance of the eukaryotic component. Limnol Oceanog 49: 168–179.

13. Countway PD, Caron DA (2006) Abundance and distribution of *Ostreococcus* sp. in the San Pedro Channel, California, as Revealed by Quantitative PCR. Appl Environ Microbiol 72: 2496–2506.

14. Piganeau G, Moreau H (2007) Screening the Sargasso Sea metagenome for data to investigate genome evolution in *Ostreococcus* (Prasinophyceae, Chlorophyta). Gene 406: 184–190.

15. Vaulot D, Le Gall F, Marie D, Guillou L, Partensky F (2004) The Roscoff Culture Collection (RCC): a collection dedicated to marine picoplankton. Nova Hedwigia 79: 49–70.

16. Li WKW (1994) Primary Production of Prochlorophytes, Cyanobacteria, and Eukaryotic Ultraphytoplankton-Measurements from Flow Cytometric Sorting. Limnol Oceanog 39: 169–175.

17. Rodriguez F, Derelle E, Guillou L, Le Gall F, Vaulot D, et al. (2005) Ecotype diversity in the marine picoeukaryote *Ostreococcus* (Chlorophyta, Prasinophyceae). Environ Microbiol 7: 853–859.

18. Palenik B, Grimwood J, Aerts A, Rouze P, Salamov A, et al. (2007) The tiny eukaryote *Ostreococcus* provides genomic insights into the paradox of plankton speciation. Proc Natl Acad Sci U S A 104: 7705–7710.

19. O’Kelly CJ, Sieracki ME, Thier EC, Hobson IC (2003) A transient bloom of *Ostreococcus* (Chlorophyta, Prasinophyceae) in West Neck Bay, Long Island, New York. J Phycol 39: 850–854.

20. Iyer LM, Balaji S, Koonin EV, Aravind L (2006) Evolutionary genomics of nucleo-cytoplasmic large DNA viruses. Virus Res 117: 156–184.

21. Lefkowitz EJ, Wang C, Upton C (2006) Poxviruses: past, present and future. Virus Res 117: 105–118.

22. Dunigan DD, Fitzgerald LA, Van Etten JL (2006) Phycodnaviruses: a peek at genetic diversity. Virus Res 117: 119–132.

23. Butcher R (1952) Contributions to our knowledge of the smaller marine algae. J Mar Biol Assoc UK 31: 175–191.

24. Knight-Jones EW, Walne PR (1951) *Chromulina-pusilla* Butcher, a Dominant Member of the Ultraplankton. Nature 167: 445–446.

25. Mayer JA, Taylor FJR (1979) Virus Which Lyses the Marine Nanoflagellate *Micromonas pusilla*. Nature 281: 299–301.

26. Waters RE, Chan AT (1982) *Micromonas pusilla* Virus-the virus growth-cycle and associated physiological events within the host-cells-host range mutation. J Gen Virol 63: 199–206.

27. Cottrell MT, Suttle CA (1991) Wide-spread occurrence and clonal variation in viruses which cause lysis of a cosmopolitan, eukaryotic marine phytoplankter, *Micromonas pusilla*. Mar Ecol-Prog Ser 78: 1–9.

28. Cottrell MT, Suttle CA (1995) Genetic Diversity of Algal Viruses Which Lyse the Photosynthetic Picoflagellate *Micromonas pusilla* (Prasinophyceae). Appl Environ Microbiol 61: 3088–3091.

29. Cottrell MT, Suttle CA (1995) Dynamics of a lytic virus infecting the photosynthetic marine picoflagellate *Micromonas pusilla*. Limnol and Oceanog 40: 730–739.

30. Brussaard CPD, Noordeloos AAM, Sandaa RA, Heldal M, Bratbak G (2004) Discovery of a dsRNA virus infecting the marine photosynthetic protist *Micromonas pusilla*. Virology 319: 280–291.

31. Attoui H, Jaafar FM, Belhouchet M, de Micco P, de Lamballerie X, et al. (2006) *Micromonas pusilla* reovirus: a new member of the family Reoviridae assigned to a novel proposed genus (Mimoreovirus). J Gen Virol 87: 1375–1383.

32. Brown CM, Campbell DA, Lawrence JE (2007) Resource dynamics during infection of *Micromonas pusilla* by virus MpV-Sp1. Environ Microbiol 9: 2720–2727.

33. Chen F, Suttle CA (1995) Amplification of DNA polymerase gene fragments from viruses infecting microalgae. Appl Environ Microbiol 61: 1274–1278.

34. Angly FE, Felts B, Breitbart M, Salamon P, Edwards RA, et al. (2006) The marine viruses of four oceanic regions. PLoS Biol 4: e368.

35. Delaroque N, Muller DG, Bothe G, Pohl T, Knippers R, et al. (2001) The complete DNA sequence of the *Ectocarpus siliculosus* Virus EsV-1 genome. Virology 287: 112–132.

36. Delaroque N, Boland W, Muller DG, Knippers R (2003) Comparisons of two large phaeoviral genomes and evolutionary implications. J Mol Evol 57: 613–622.

37. Wilson WH, Schroeder DC, Allen MJ, Holden MT, Parkhill J, et al. (2005) Complete genome sequence and lytic phase transcription profile of a Coccolithovirus. Science 309: 1090–1092.

38. Allen MJ, Schroeder DC, Donkin A, Crawfurd KJ, Wilson WH (2006) Genome comparison of two Coccolithoviruses. Virol J 3: 15.

39. Van Etten JL, Meints RH (1999) Giant viruses infecting algae. Annu Rev Microbiol 53: 447–494.

40. Fitzgerald LA, Graves MV, Li X, Feldblyum T, Hartigan J, et al. (2007) Sequence and annotation of the 314-kb MT325 and the 321-kb FR483 viruses that infect Chlorella Pbi. Virology 358: 459–471.

41. Fitzgerald LA, Graves MV, Li X, Feldblyum T, Nierman WC, et al. (2007) Sequence and annotation of the 369-kb NY-2A and the 345-kb AR158 viruses that infect Chlorella NC64A. Virology 358: 472–484.

42. Guillou L, Eikrem W, Chretiennot-Dinet MJ, Le Gall F, Massana R, et al. (2004) Diversity of picoplanktonic prasinophytes assessed by direct nuclear SSU rDNA sequencing of environmental samples and novel isolates retrieved from oceanic and coastal marine ecosystems. Protist 155: 193–214.

43. Ryabova LA, Pooggin MM, Hohn T (2002) Viral strategies of translation initiation: ribosomal shunt and reinitiation. Prog Nucleic Acid Res Mol Biol 72: 1–39.

44. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, et al. (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25: 3389–3402.

45. Iyer LM, Aravind L, Koonin EV (2001) Common origin of four diverse families of large eukaryotic DNA viruses. J Virol 75: 11720–11734.

46. Pringle CR (1998) The universal system of virus taxonomy of the International Committee on Virus Taxonomy (ICTV), including new proposals ratified since publication of the Sixth ICTV Report in 1995. Arch Virol 143: 203–210.

47. Van Etten JL, Graves MV, Muller DG, Boland W, Delaroque N (2002) Phycodnaviridae-large DNA algal viruses. Archives of Virology 147: 1479–1516.

48. Flor HH (1971) Current status of the gene-for-gene concept. Annu Rev Phytopathol 9: 275–296.

49. Borevitz JO, Hazen SP, Michael TP, Morris GP, Baxter IR, et al. (2007) Genome-wide patterns of single-feature polymorphism in *Arabidopsis thaliana*. Proc Natl Acad Sci USA 104: 12057–12062.

50. Nagasaki K, Takao Y, Shirai Y, Mizumoto H, Tomaru Y (2005) [Molecular ecology of microalgal viruses]. Jиру 55: 127–132.

51. Baudoux AC, Brussaard CP (2005) Characterization of different viruses infecting the marine harmful algal bloom species *Phaeocystis globosa*. Virology 341: 80–90.

52. Bouvier T, del Giorgio PA (2007) Key role of selective viral-induced mortality in determining marine bacterial community composition. Environ Microbiol 9: 287–297.

53. Parada V, Sintes E, van Aken HM, Weinbauer MG, Herndl GJ (2007) Viral abundance, decay, and diversity in the meso- and bathypelagic waters of the north atlantic. Appl Environ Microbiol 73: 4429–4438.

54. Weinbauer MG, Hornak K, Jezbera J, Nedoma J, Dolan JR, et al. (2007) Synergistic and antagonistic effects of viral and protistan grazing on bacterial biomass, production and diversity. Environ Microbiol 9: 777–788.

55. Zhang R, Weinbauer MG, Qian PY (2007) Viruses and flagellates sustain apparent richness and reduce biomass accumulation of bacterioplankton in coastal marine waters. Environ Microbiol 9: 3008–3018.

56. Henderson GP, Gan L., Jensen GJ (2007) 3-D ultrastructure of *O. tauri*: electron cryotomography of an entire eukaryotic cell. PLoS ONE 2: e749.

57. DuRand MD, Green RE, Sosik HM, Olson RJ (2002) Diel variations in optical properties of *Micromonas pusilla* (Prasinophyceae). J Phycol 38: 1132–1142.

58. Renesto P, Abergel C, Decloquement P, Moinier D, Azza S, et al. (2006) Mimivirus giant particles incorporate a large fraction of anonymous and unique gene products. J Virol 80: 11678–11685.

59. Rohozinski J, Girton LE, Van Etten JL (1989) Chlorella viruses contain linear nonperturbed double-stranded DNA genomes with covalently closed hairpin ends. Virology 168: 363–369.

60. Strasser P, Zhang YP, Rohozinski J, Van Etten JL (1991) The termini of the chlorella virus PBCV-1 genome are identical 2.2-kbp inverted repeats. Virology 180: 763–769.

61. Zhang Y, Strasser P, Grabherr R, Van Etten JL (1994) Hairpin loop structure at the termini of the chlorella virus PBCV-1 genome. Virology 202: 1079–1082.

<figure><img src="image_21.png" /></figure>
62. Patel DD, Pickup DJ, Joklik WK (1986) Isolation of cowpox virus A-type inclusions and characterization of their major protein component. Virology 149: 174-189.

63. Xia YN, Burbank DE, Uher L, Rabussay D, Van Etten JL (1986) Restriction endonuclease activity induced by PBCV-1 virus infection of a Chlorella-like green alga. Mol Cell Biol 6: 1430-1439.

64. Raoult D, Audic S, Robert C, Abergel C, Renesto P, et al. (2004) The 1.2-megabase genome sequence of Mimivirus. Science 306: 1344-1350.

65. Agarkova IV, Dunigan DD, Van Etten JL (2006) Virion-associated restriction endonucleases of chloroviruses. J Virol 80: 8114-8123.

66. Bidle KD, Falkowski PG (2004) Cell death in planktonic, photosynthetic microorganisms. Nature Reviews Microbiology 2: 643-655.

67. Oliveira IC, Brenner E, Chiu J, Hsieh MH, Kouranov A, et al. (2001) Metabolite and light regulation of metabolism in plants: lessons from the study of a single biochemical pathway. Braz J Med Biol Res 34: 567-575.

68. Silva P, Thompson E, Bailey S, Kruse O, Mullineaux CW, et al. (2003) FtsH is involved in the early stages of repair of photosystem II in Synechocystis sp PCC 6803. Plant Cell 15: 2152-2164.

69. Hoekstra FA, Golovina EA, Buitink J (2001) Mechanisms of plant desiccation tolerance. Trends Plant Sci 6: 431-438.

70. Kohl DH, Schubert KR, Carter MB, Hagedorn CH, Shearer G (1988) Proline metabolism in N2-fixing root nodules: energy transfer and regulation of purine synthesis. Proc Natl Acad Sci USA 85: 2036-2040.

71. Cho K, Winans SC (1996) The *putA* gene of *Agrobacterium tumefaciens* is transcriptionally activated in response to proline by an Lrp-like protein and is not autoregulated. Mol Microbiol 22: 1025-1033.

72. Rocha EP, Danchin A (2002) Base composition bias might result from competition for metabolic resources. Trends Genet 18: 291-294.

73. Simek K, Weinbauer MG, Hornak K, Jezbera J, Nedoma J, et al. (2007) Grazer and virus-induced mortality of bacterioplankton accelerates development of *Flectobacillus* populations in a freshwater community. Environ Microbiol 9: 789-800.

74. Brussaard CPD (2004) Optimization of procedures for counting viruses by flow cytometry. Appl Environ Microbiol 70: 1506-1513.

75. Trask BJ, van den Engh GJ, Elgershuizen JH (1982) Analysis of phytoplankton by flow cytometry. Cytometry 2: 258-264.

76. Mead JR, Arrowood MJ, Current WL, Sterling CR (1988) Field inversion gel electrophoretic separation of *Cryptosporidium* spp. chromosome-sized DNA. J Parasitol 74: 366-369.

77. Ewing B, Hillier L, Wendl MC, Green P (1998) Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res 8: 175-185.

78. Rutherford K, Parkhill J, Crook J, Horsnell T, Rice P, et al. (2000) Artemis: sequence visualization and annotation. Bioinformatics 16: 944-945.

79. Carver TJ, Rutherford KM, Berriman M, Rajandream MA, Barrell BG, et al. (2005) ACT: the Artemis Comparison Tool. Bioinformatics 21: 3422-3423.

80. Finn RD, Tate J, Mistry J, Coggill PC, Sammut SJ, et al. (2008) The Pfam protein families database. Nucleic Acids Res 36: D281-288.

81. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ (1990) Basic local alignment search tool. J Mol Biol 215: 403-410.

82. Schroeder DC, Oke J, Malin G, Wilson WH (2002) Coccolithovirus (Phycodnaviridae): characterisation of a new large dsDNA algal virus that infects *Emiliana huxleyi*. Arch Virol 147: 1685-1698.

83. La Scola B, Audic S, Robert C, Jungang L, de Lamballerie X, et al. (2003) A giant virus in amoebae. Science 299: 2033.

84. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG (1997) The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 25: 4876-4882.

85. Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572-1574.

86. Guindon S, Gascuel O (2003) A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 52: 696-704.

87. Guindon S, Lethiec F, Duroux P, Gascuel O (2005) PHYML Online-a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Res 33: W557-W559.

88. Posada D, Crandall KA (1998) MODELTEST: testing the model of DNA substitution. Bioinformatics 14: 817-818.

89. Abascal F, Zardoya R, Posada D (2005) ProtTest: selection of best-fit models of protein evolution. Bioinformatics 21: 2104-2105.

90. Fuller NJ, Campbell C, Allen DJ, Pitt FD, Zwirglmaier K, et al. (2006) Analysis of photosynthetic picocystaryote diversity at open ocean sites in the Arabian Sea using a PCR biased towards marine algal plastids. Aquatic Microbial Ecol 43: 79.

91. Yamada T, Onimatsu H, Van Etten JL (2006) Chlorella viruses. Adv Virus Res 66: 293-336.

<figure><img src="image_22.png" /></figure>