# Therapeutic Potential of a New Jumbo Phage That Infects Vibrio coralliilyticus, a Widespread Coral Pathogen

Loïc Jacquemot, Yvan Bettarel, Joanne Monjol, Erwan Corre, Sebastien

Halary, Christelle Desnues, Thierry Bouvier, Christine Ferrier-Pagès,

Anne-Claire Baudoux

► To cite this version:

Loïc Jacquemot, Yvan Bettarel, Joanne Monjol, Erwan Corre, Sebastien Halary, et al.. Therapeutic Potential of a New Jumbo Phage That Infects Vibrio coralliilyticus, a Widespread Coral Pathogen. Frontiers in Microbiology, 2018, 9, pp.2501. 10.3389/fmicb.2018.02501. hal-01913634

HAL Id: hal-01913634

https://hal.sorbonne-universite.fr/hal-01913634

Submitted on 6 Nov 2018

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
<figure><img src="image_2.png" /></figure>

Check for
updates

# Therapeutic Potential of a New Jumbo Phage That Infects Vibrio coralliilyticus, a Widespread Coral Pathogen

Loïc Jacquemot¹†, Yvan Bettarel², Joanne Monjol¹, Erwan Corre³, Sébastien Halary⁴, Christelle Desnues⁴, Thierry Bouvier², Christine Ferrier-Pagès⁵ and Anne-Claire Baudoux¹*

**OPEN ACCESS**

**Edited by:**
Robert Czajkowski,
University of Gdansk, Poland

**Reviewed by:**
Inmaculada Garcia-Heredia,
University of Alicante, Spain
Simon Roux,
Joint Genome Institute (JGI),
United States
Ahmed Askora,
Zagazig University, Egypt

**\*Correspondence:**
Anne-Claire Baudoux
acbaudoux@sb-rocoff.fr

**†Present Address:**
Loïc Jacquemot,
Département de Biologie, Institut de
Biologie Intégrative et des Systèmes
(IBIS), Université Laval, Québec, QC,
Canada

**Specialty section:**
This article was submitted to
Virology,
a section of the journal
Frontiers in Microbiology

**Received:** 26 July 2018
**Accepted:** 01 October 2018

**Published:** 24 October 2018

**Citation:**
Jacquemot L, Bettarel Y, Monjol J,
Corre E, Halary S, Desnues C,
Bouvier T, Ferrier-Pagès C and
Baudoux A-C (2018) Therapeutic
Potential of a New Jumbo Phage That
Infects Vibrio coralliilyticus, a
Widespread Coral Pathogen.
Front. Microbiol. 9:2501.
doi: 10.3389/fmicb.2018.02501

¹ Sorbonne Universités UPMC Paris 06, CNRS, UMR7144 Adaptation et Diversité en Milieu Marin, Station Biologique de Roscoff, Roscoff, France, ² MARBEC, Université Montpellier, IRD, CNRS, Ifremer, Montpellier, France, ³ Sorbonne Universités UPMC Paris 06, CNRS, FR2424 Fédération de Recherche, Station Biologique de Roscoff, Roscoff, France, ⁴ Aix Marseille Université, Microbes, Evolution Phylogeny and infection (MEPHI), CNRS FRE2013, IRD 198, AP-HM, IHU - Méditerranée Infection, Marseille, France, ⁵ Centre Scientifique de Monaco, Equipe Ecophysiologie Coralienne, Monaco, Monaco

Biological control using bacteriophages is a promising approach for mitigating the devastating effects of coral diseases. Several phages that infect *Vibrio coralliilyticus*, a widespread coral pathogen, have been isolated, suggesting that this bacterium is permissive to viral infection and is, therefore, a suitable candidate for treatment by phage therapy. In this study, we combined functional and genomic approaches to evaluate the therapeutic potential of BONAISHI, a novel *V. coralliilyticus* phage, which was isolated from the coral reef in Van Phong Bay (Vietnam). BONAISHI appears to be strictly lytic for several pathogenic strains of *V. coralliilyticus* and remains infectious over a broad range of environmental conditions. This candidate has an unusually large dsDNA genome (303 kb), with no genes that encode known toxins or implicated in lysogeny control. We identified several proteins involved in host lysis, which may offer an interesting alternative to the use of whole bacteriophages for controlling *V. coralliilyticus*. A preliminary therapy test showed that adding BONAISHI to an infected culture of *Symbiodinium* sp. cells reduced the impact of *V. coralliilyticus* on *Symbiodinium* sp. photosynthetic activity. This study showed that BONAISHI is able to mitigate *V. coralliilyticus* infections, making it a good candidate for phage therapy for coral disease.

**Keywords:** phage therapy, coral disease, *Vibrio coralliilyticus*, viral genomics, phage-host interactions

## INTRODUCTION

Coral reefs are one of the most productive and diversified ecosystems on the planet (Connell, 1978) and they provide a wealth of ecological services as well as being economically important, supporting fisheries, tourism, and medical applications (Moberg and Folke, 1999; Hughes et al., 2003; Cooper et al., 2014). The health of these ecosystems is severely threatened by the combined effect of local anthropogenic pressures and global changes (Jackson et al., 2001; Hughes et al., 2003; Pandolfi et al., 2003; Bellwood et al., 2004). Over-exploitation of marine species, pollution, and increased sea surface temperature are associated with the emergence of coral diseases, which are contributing to the decline of coral reefs worldwide.
Several studies have identified Vibrio spp. (γ-Proteobacteria) as causative agents of coral bleaching for multiple coral species and in multiple locations (Ushijima et al., 2014; reviewed in Mera and Bourne, 2018). Vibrio coralliilyticus (V. coralliilyticus) has emeregd as an important bacterial pathogen model for understanding the establishment and propagation of coral disease (Sussman et al., 2008; O'Santos et al., 2011; Garren et al., 2014; Pollock et al., 2015). Studies have shown that V. coralliilyticus infection is temperature dependent and infects the coral endosymbiont *Symbiodinum* through the production of proteases that inhibit photosynthesis. This results in the loss of the endosymbiont from the coral tissues and ultimately leads to coral bleaching (Ben-Haim et al., 2003; Sussman et al., 2008, 2009; Cohen et al., 2013). With the increasing devastation of coral reefs, the development of new tools and strategies to control pathogens and treat diseased corals is becoming a major issue. Currently, biocontrol strategies, such as phage therapy, are being seriously evaluated for mitigating coral diseases (Efrony et al., 2009; Teplitski and Ritchie, 2009; Atad et al., 2012; Cohen et al., 2013).

The potential of viruses (more specifically bacteria viruses also referred to as bacteriophages or phages) as therapeutic agents to treat infectious diseases has been known for a long time (d'Herelle, 1926; Duckworth, 1976; Duckworth and Gulig, 2002). The idea of phage therapy arose from the early discovery that a given virus usually infects a single host species, leaving the rest of the microbial community untouched. Moreover, viruses are obligate intracellular organisms and, therefore, their production is self-regulated and limited by the availability of suitable hosts. Over the past decade, there have been promising *in vitro* and *in situ* trials of phage therapy for corals. One used BA3, a virulent bacteriophage that infect the causative agent of white plague disease (Efrony et al., 2007, 2009), to inhibit the progression of the disease and its transmission to healthy corals (Atad et al., 2012). Although this research is still at an early stage, this promising result in the open sea suggest that *in situ* phage therapy for coral diseases is achievable (Atad et al., 2012). The isolation and characterization of new bacteriophages is, therefore, essential to increase the collection of potential candidates for therapeutic assays.

The therapeutic value of a candidate bacteriophage relies on the characterization of viral properties such as the virion stability, growth kinetics, viral yield, and host range. Understanding the lifestyle of the candidate phage is probably the key to their use in therapy. Only virulent bacteriophages, which replicate through a lytic cycle and kill their host after infection, will be suitable candidates. Temperate bacteriophages, which replicate using a lysogenic cycle, may improve host fitness through gene transfer. It is, however, difficult to distinguish between virulent and temperate phage because temperate viruses can switch to a lytic cycle in response to environmental changes, such as temperature, pH salinity, UV, pollution, or nutrient availability (Jiang and Paul, 1996; Williamson and Paul, 2006, reviewed in Howard-Varona et al., 2017). Furthermore, infection dynamics can be highly variable, even between two closely related hosts (Holmfeld et al., 2014; Dang et al., 2015). Over the past decade, genomics has greatly improved our understanding of host—virus interactions,

and is the key to establishing whether a candidate is an obligate lytic bacteriophage (Howard-Varona et al., 2017). Bacteriophage genome sequencing is also essential for evaluating the safety of a candidate (absence of toxins and temperate phage hallmarks), and to provide information on the candidate's evolution history and ecology (Lobocka et al., 2014).

In this study, we report a novel bacteriophage, hereafter referred to as Vibrio phage BONAISHI that infects the model coral pathogen V. coralliilyticus. We studied it using a combination of functional, genomic, and metagenomic approaches to evaluate the potential of this phage for mitigating disease caused by V. coralliilyticus.

# MATERIAL AND METHODS

## Virus Isolation

Seawater samples were collected from coral surrounding water off Whale Island (Van Phong Bay, Vietnam). A 50 mL aliquot was supplemented with 10% (v/v) Marine Broth (MB, Difco) and the mixture was enriched with 1 mL Vibrio coralliilyticus LMG20984 (YB1) culture and incubated for 48 h at 25°C (Brussaard et al., 2016). The sample was clarified (7,000 g, 15 min) and the supernatant was filtered through 0.2 µm PES filters to separate the viral community. A 100 µL aliquot of the filtrate was added to 900 µL YB1 culture and incubated for 30 min at 25°C. The mixture was included in molten agar (Marine Broth supplemented with 0.6% noble agar) and spread on a Marine Agar plate. After 48 h incubation, a translucent plaque indicating host lysis was removed from the bacterial lawn, eluted in 0.22 µm filtered Salt Marine (SM) buffer (100 mM NaCl, 8 mM MgSO4, 50 mM Tris-HCl pH 8.0) and combined with a host culture in a plaque assay (Brussaard et al., 2016). This procedure was repeated twice to ensure the clonality of the bacteriophage. Finally, the clonal phage suspension and host culture were used in a plaque assay giving confluent lysis. The plaques were eluted in SM buffer, the eluent was clarified by centrifugation (7,000 g, 30 min, 4°C), and the supernatant was filtered through 0.22 µm and stored at 4°C until use.

## Transmission Electronic Microscopy

A 10 µL aliquot of the phage suspension was loaded onto a Formvar/carbon film coated 400 mesh copper grid (Euromedex). After 5 min incubation, the grid was blotted with filter paper and stained with 2% uranyl acetate for 30 s, blotted again to remove excess dye and air dried for 30 min (Ackermann and Heldal, 2010). The specimen was imaged using JEOL 1400 transmission electron microscope operating at 120 keV at a magnification of 80,000X.

## Environmental Range of Infectivity

The tolerance of BONAISHI to temperature and pH ranges was evaluated by monitoring the loss of infectivity of a freshly produced suspension by spot test. For evaluating the temperature range, 100 µL of viral suspension (5 x 10⁸ PFU/mL) was incubated at temperatures from 4 to 70°C for 30 min in a dry bath. The samples were then cooled for 5 min at 4°C and virus infectivity was assessed by spot test. Briefly, 5 µL of the treated
viral suspension was spotted on a host lawn obtained by plating
a 1:4 mixture of host culture in molten agar onto a marine agar
plate. For evaluating the pH range, 100 µL of viral suspension
were added to 900 µL SM buffer adjusted at pH ranging from
2 to 10. The samples were incubated for 24 h at 4°C and then
spot-tested as described above.

Host Range

A selection of 43 bacterial strains (Table S1) related to the
original host *V. coralliilyticus* YB1 were used to determine the
host specificities of BONAISHI. The ability of BONAISHI phage
to infect these strains was determined by pairwise infection. The
bacterial strains were grown in Marine Broth media (DIFCO)
overnight. Each culture was included in molten agar (Marine
Broth supplemented with 0.6% noble agar) and spread on a
Marine Agar plate. A freshly produced suspension of BONAISHI
was serially diluted in SM buffer and a 5 µL drop of each dilution
was spotted on a bacterial lawn. After 24-48h incubation at 25°C,
the formation of translucent spots, indicative of host lysis, was
recorded.

Life Strategy

The growth cycle of BONAISHI was tested on two *V. coralliilyticus* strains including the original host LMG 20984 (YB1) and the alternate host LMG 23696 (P1) (Middelboe et al., 2010). Both host cultures were grown in MB and divided into four 25 mL sub-cultures. One sub-culture served as control and the remaining 3 were inoculated with a freshly produced BONAISHI suspension at multiplicity of infection (MOI) of 0.1 as determined by flow cytometry (see below). All cultures were incubated at 25°C for 48 h. Samples for viral and bacterial counts were taken every hour for 10 h, and then every 4 h for 48 h. Samples were immediately fixed with glutaraldehyde (0.5% final concentration) for 10 min at 4°C, flash frozen in liquid nitrogen, and stored at -80°C until flow cytometry analysis (see below). Bacterial host and virus counts were used to calculate the phage latent period and burst size. The latent period corresponds to the time elapsed between the viral inoculation and the release of virions. The burst-size, which corresponds to the number of virions produced per infected host cell, was determined by the ratio of the net increase in virus concentration over the net decrease in host concentration during the first burst.

Flow Cytometry

For determining the bacterial abundance, samples were diluted
in autoclaved 0.2 µm filtered TE buffer (10 mM Tris-HCl, 1 mM
EDTA, pH 8.0) and stained with a SYBR Green (10,000-fold
dilution of commercial solution) for 15 min in the dark at
ambient temperature. For determining the viral abundance,
100-1,000-fold diluted samples were stained with SYBR Green
(20,000-fold dilution of commercial solution) for 10 min in the
dark at 80°C. Analyses were performed using a FACS Canto II
equipped with an argon-laser (455 nm). The trigger was set on
the green fluorescence and the sample was delivered at a rate
of 0.06 mL min⁻¹ and analyzed for 1 min (Brussaard, 2004).
Viral and bacterial counts were corrected for a blank consisting

of TE-buffer with autoclaved 0.2 µm filtered seawater at the
corresponding dilution.

Genome Extraction and Sequencing

A 1 L viral suspension was concentrated by ultrafiltration using a 30 kDa PES membrane (Vivaflow 50, Vivascience) and centrifugal concentrator (Vivaspin 20, 30 kDa, PES) to a final volume of 2 mL. The concentrate was subsequently purified on linear sucrose gradient (10–40 % in 0.2 µm SM) by ultracentrifugation (SW41.Ti rotor, 96,808 g, 30 min at 4°C). The BONAISHI particles formed a well resolved band that was extracted, dialyzed against SM buffer using a centrifugal concentrator (Vivaspin 20, 30 kDa, PES) and stored at 4°C until use. The genome of the purified phage suspension was extracted using a DNAeasy Blood and Tissue kit (QUIAGEN, Valencia, CA) according to the manufacturer's protocol. Samples were sent to GATC Biotech, and sequenced using PACBIO RS II (17,131 mean read length). Raw read sequences assembled as a single contig using HGAP software (Chin et al., 2013). The final draft assembly was 303, 340 bp with an average coverage of 1,720x and an average base quality score of 86%.

Bioinformatic Analysis

Phage Genome Annotation

Putative coding DNA sequences (CDS) in the BONAISHI
genome were predicted using Glimmer (Delcher et al., 1999)
and Genemark (Besemer et al., 2001). The coordinates of each
translated open reading frames (ORF) were also inspected
manually. ORF smaller than 200 base pairs (bp) were removed
from the analysis. The predicted amino acid sequences were
assigned manually by BLASTP and PSI-BLAST (cutoff e-value
<10⁻⁵) against the NCBI non-redundant database (January
2018) and InterProScan 5 (Jones et al., 2014), as well as the fully
automated RAST server annotation service (Aziz et al., 2008).
BLASTP was used to search for putative toxins in the databases
MvirDB (Zhou et al., 2006), VirulenceFinder (Chen et al., 2005),
Vibrio-base (Choo et al., 2014), and t3DB (Lim et al., 2009;
Wishart et al., 2015) toxin databases. The IntegrallDB (Moura
et al., 2009) and ACLAME (Leplae et al., 2006) databases were
used to check for prophage-like sequences. tRNAScanSE (Lowe
and Eddy, 1997) and Aragorn (Laslett and Canback, 2004) were
used to check for tRNA. The genome map was produced using
Artemis (Rutherford et al., 2000) and DNA Plotter (Carver et al.,
2009). The BONAISHI genome sequence has been submitted to
the GenBank database under accession number MH595538.

Terminase Large Subunit (TerL) Protein Phylogeny

The amino acids sequence of the terminase large subunit
from 63 jumbo phages including BONAISHI were used for
phylogenetic analysis. Sequences were trimmed to 406 bp, the
minimum sequence length of *Aeromonas* phage px29, using
BioEdit (Hall, 1999). Sequences were aligned by Muscle and
the tree was constructed by Maximum Likelihood with 1,000
bootstrap iterations using Mega 6.06 (Tamura et al., 2013).
## Host Genome Analysis

The clustered regularly interspaced short palindromic repeats (CRISPR), in the pathogenic *V. coralliilyticus* strains P1 and YB1 genome, were searched for genetic signatures of viral resistance mechanisms using CRISPRfinder (Grissa et al., 2007). Genetic exchange between the phage BONAISHI and its bacterial hosts was checked by homology between the phage ORFs and the ORFs of *V. coralliilyticus* P1 (AEQS00000000) and YB1 (ACZN00000000) using BLASTP.

## Metagenomic Analysis

To determine the distribution of BONAISHI, the genome was used to recruit homologous reads from 56 coral-associated virome in Metavir (Roux et al., 2011), and 137 CAMERA Broad Phage metagenomes (Seshadri et al., 2007) and viral contigs in IMG/VR (Paez-Espino et al., 2017, Table S2). These samples comprised a wide range of marine environments including tropical and temperate pelagic ecosystems, healthy and diseased coral reefs including slurry from individual coral colonies, coral mucus, and the water from coral reefs. We also carried out recruitments in prokaryote metagenome from 4 coral atolls (Dinsdale et al., 2008) to check whether BONAISHI genome sequences were inserted into prokaryote genomes. Each reads served as a query and was assigned to a (single) best-matching hit by BLASTN and TBLASTX if the alignments met the following criteria: e-value < 10⁻³, alignment length > 50, bitscore > 40. BLASTN parameters were set to: open gap cost = 8, extend gap cost = 6, match reward = 5, mismatch penalty = -4, word size = 8. Reads were recruited from each metagenome in order to determine the fraction of recruited reads that can be assigned to BONAISHI.

## Preliminary Treatment of Diseased Symbiodinium

Culture of *Symbiodinium* sp. cells (Clade A1) originally extracted from the scleractinian coral *Galaxea fascicularis* (Goiran et al., 1996) were maintained in the laboratory in F/2 medium (Guillard and Ryther, 1962) prepared from Guillard's Marine Water Enrichment Solution (Sigma-Aldrich G9903). Cultures were maintained at 25°C under 100 µmol photons m⁻² s⁻¹ of white light provided by fluorescent tubes (Mazda 18WJr/865) using a 12:12 light:dark cycle. One day prior to the therapy assay, exponentially growing *Symbiodinium* sp. cultures were transferred to 30°C under the same light conditions. A 20 mL aliquot was concentrated at 5,000 g for 10 min at 30°C (VIVASPIN 20, PES, 30 kDa). The retentate was gently resuspended in 20 mL EDTA free F/2 medium and the procedure was repeated twice to wash the culture. The *Symbiodinium* sp. abundance was determined by flow cytometry and adjusted to 10⁴ cells mL⁻¹. *V. coralliilyticus* YB1 was cultured overnight in MB and then purified in the same way. Bacterial abundance was determined by flow cytometry and adjusted to 10⁷ cells mL⁻¹. A freshly produced suspension of BONAISHI was purified by sucrose gradient and diluted in SM buffer. The viral abundance was determined by flow cytometry.

For the therapy assay, the algal culture was split into 3 equal sub-cultures. One sub-culture served as control, while two of

the subcultures were inoculated with an equal volume of *V. coralliilyticus* YB1. Of these, one was also inoculated with 10⁸ phages mL⁻¹. All three treatments were sampled at 0, 5, 20, and 60 min to determine the photosystem II quantum yield of *Symbiodinium* sp. cells using a pulse amplitude modulated fluorimeter (Phyto-PAM, Walz) connected to a chart recorder (Labpro, Vernier). After 5 min relaxation in darkness, the non-actinic modulated light (450 nm) was turned on in order to measure the fluorescence basal level, F₀. A saturating red light pulse (655 nm, 4 000 µmol quanta m⁻² s⁻¹, 400 ms) was applied to determine the maximum fluorescence level in the dark adapted sample, Fₘ. The maximal photosystem II fluorescence quantum yield of photochemical energy conversion, Fᵥ/Fₘ, was calculated using the following formula:

$$ \frac{F_V}{F_M} = \frac{(F_M - F_0)}{F_M} \qquad (1) $$

## RESULTS

### Morphology

The Vibrio phage BONAISHI formed relatively large, round plaques on *V. coralliilyticus* YB1 and produced high titer suspension. TEM microscopy showed that BONAISHI has an isometric capsid of 120 nm in diameter connected to a 190 nm long, contractile tail (Figure 1). This indicates that BONAISHI belongs to the order of the *Caudovirales* and the family of the *Myoviridae*.

### Tolerance to Environmental Factors

The incubations showed that BONAISHI can tolerate a pH ranging from 3 to 10 and temperatures ranging from 4 to 50°C without loss of infectivity (Table 1). This high tolerance suggests that the phage is very stable in the environment.

### Host Specificities

Spot tests using a broad range of potential hosts indicated that BONAISHI was able to infect and lyse several strains of

<figure><img src="image_3.png" /><figcaption>FIGURE 1 | Morphology of Vibrio phage BONAISHI. (A) BONAISHI forms large, round plaque on a lawn of *V. coralliilyticus* YB1 on 0.6% soft agar. (B) Transmission electron micrographs of a negatively stained particle of bacteriophage BONAISHI. The icosahedral head (120 nm in diameter) and the long, contractile tail (190 nm in length) suggest that BONAISHI belongs to the *Myoviridae* family.</figcaption></figure>
<table><caption>TABLE 1 | Tolerance of Vibrio phage BONAISHI to temperature and pH.</caption><thead><tr><th>Treatments</th><th>Infectivity</th></tr></thead><tbody><tr><td><strong>TEMPERATURE (°C)</strong></td><td></td></tr><tr><td>4</td><td>+</td></tr><tr><td>15</td><td>+</td></tr><tr><td>20</td><td>+</td></tr><tr><td>25</td><td>+</td></tr><tr><td>30</td><td>+</td></tr><tr><td>35</td><td>+</td></tr><tr><td>40</td><td>+</td></tr><tr><td>45</td><td>+</td></tr><tr><td>50</td><td>+</td></tr><tr><td>60</td><td>-</td></tr><tr><td>70</td><td>-</td></tr><tr><td><strong>pH</strong></td><td></td></tr><tr><td>2</td><td>-</td></tr><tr><td>3</td><td>+</td></tr><tr><td>4</td><td>+</td></tr><tr><td>5</td><td>+</td></tr><tr><td>6</td><td>+</td></tr><tr><td>7</td><td>+</td></tr><tr><td>8</td><td>+</td></tr><tr><td>9</td><td>+</td></tr><tr><td>10</td><td>+</td></tr></tbody></table>

V. coralliilyticus of interest (Table S1). Besides V. coralliilyticus YB1, BONAISHI can infect another known coral pathogen, V. coralliilyticus P1, as well as V. coralliilyticus LMG21348, isolated from a bleached coral colony (Pocillopora damicornis), and V. coralliilyticus 1H13, isolated from the mucus of a Fungia specimen. This phage did not lyse of any of the closely related species in the test, suggesting that it is species-specific.

## Life Cycle

One-step growth experiments showed that the phage readily propagated on each of its hosts with a latent period of 2–3 h (strains P1 and YB1) and burst size of 8 (P1) and 19 (YB1) (Figure 2). Nearly all the virions produced (96%) were infectious virions per infected cell. The infected host culture collapsed rapidly and there was total lysis 10 h after inoculation.

## General Features of BONAISHI Genome

The genome of BONAISHI consists of a large double-stranded circularly permuted DNA sequence of 303,340 base-pair (bp) with a % G+C content of 42.5% (Figure 3). Terminal duplications at the extremities of the assembled sequence are 14,373 bp, giving a non-redundant genome of 288,967 bp. Glimmer and Genemark predicted 301 putative ORFs, which comprised 93.8% of the total sequence and were mostly oriented in a single direction. Most ORFs initiate translation at an ATG start codon except for 6 ORFs, 4 of which have a GTG start codon and the other 2 a TTG start codon. We did not find any tRNA in the genome. Of the 301 predicted ORFs, 110 ORFs (36.4%)

had significant homologs in public databases and a biochemical function could be assigned for 62 of these (Table 2, Figure 3). Blast searches revealed that most of the ORFs (66/110) with significant homologs were closely related to other Myoviridae (Table 2). Most best-hits corresponded to members of the genus Phikzvirus, which comprises myoviruses with a large genome (> 200 kb), including Pseudomonas phages Phikz and PhiPA3, Erwinia phage Ea35-70, Ralstonia phage RSF1 and RSL2 and Vibrio phage pVa21, VP4B, pTD1, Phabio, and Noxifer. Although many gene functions could be assigned, most of the predicted genes are ORFans that are unique to BONAISHI.

## Gene Annotation

### Proteins Involved in Virion Structure and Assembly

The predicted proteins involved in the virion structure and assembly included a major and an accessory capsid protein (ORF 101 and 113), several tail components including a tail tube (ORF 36), a tail shealth (ORF 37), and an accessory tail protein (ORF 300), as well as several conserved structural proteins (ORFs 7, 15, 30, 39, 40, 97, 99, 100, 102, 108, 110, 119, 120, 121, and 301). We also identified a terminase large subunit (ORF 41) involved in DNA packaging. All these structural components shared strong homologies with proteins encoded by other members of Phikzvirus genus (Table 1, Figure S1). Phylogenetic analyses based on the amino acid sequence of the terminase large subunit classified BONAISHI as a singleton (Figure 4). The closest sequences belong to Vibrio phage pVa-21 and a cluster with Cronobacter phage CR5, Erwinia phage vB EamM, phiEaH2 and Salmonella phage SPN3US. The second closest cluster comprised Phikzvirus Pseudomonas phage KTN4, Phikz, Noxifer, PA3, Phabio, and 201 phi2-1.

### Proteins Involved in DNA Replication, Recombination and Repair

The BONAISHI genome encoded at least 10 proteins involved in DNA replication, recombination, and repair. These included a DNA polymerase B (ORF 96), putative DnaB, and DEAD-like helicases (ORF 85 and 111), NAD-dependent DNA ligase (ORF 198), SbcC-like nucleases (ORF2 and 76), a ribonuclease H (ORF 10), Holliday junction resolvase (ORF 118), uvsX recombinase (ORF12). ORF 235 corresponded to the HNH family of homing endonuclease located between genes encoding β subunits of ribonucleotide diphosphate reductase.

### Proteins Involved in Nucleotide Metabolism and DNA Modification

We were able to assign a putative function to 10 enzymes involved into nucleotide metabolism. Predicted proteins included two ribonucleotide diphosphate reductase (RDR) α subunits (ORFs 232 and 233) and two RDR β subunits (ORFs 234 and 236), 4 proteins of the pyridine nucleotide salvage pathway corresponding to a nicotinamide riboside transporter (ORF 181), a nicotinamide phosphoribosyltransferase (ORF 177) as well as nicotinamide mononucleotide adenylyltransferases (ORF 172 and 176). BONAISHI also encodes 2 proteins involved in thymidine biosynthesis one of which is a thymidylate synthase complementing protein (ORF 269) and the other a thymidylate
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | One-step growth experiments of BONAISHI on *Vibrio coralliilyticus* strain P1 (A,B) and YB1 (C,D). Bacteria and bacteriophage counts were assessed by flow cytometry upon staining with the nucleic acid dye SYBR Green (Life technology) according to Brussaard (2004).</figcaption></figure>

kinase (ORF 231). We did not identify any gene for protein involved in DNA modification.

## Proteins Involved in DNA Transcription

A transcription regulator related to the PadR family (ORF 212) was found as well as two sets of multisubunit RNA polymerase genes (β- and β'-RNAP subunits). ORFs 24, 25, 26, 92 corresponded to the Phikz virion-associated RNAP and ORFs 72, 83, 84 shared significant homologies with the Phikz genes encoding the early expressed RNAP (Table 2, Figure S1). No homologs of the early expressed Phikz β-RNAP subunits were found in BONAISHI genome. ORF 49 was assigned to an RNA binding protein.

## Lysis, Host-Phage Interaction, and Lysogeny

We were able to annotate two proteins (ORFs 23 and 138) involved into host-virus interactions. There were good blastp hits on proteins of unknown function in the NCBI nr database that could be expanded to known glycoside hydrolases (GH) in the expert database CAZY (Carbohydrate Active enZYmes db, http://www.cazy.org). These included a glycoside hydrolase with a conserved endopeptidase domain (ORF 23) affiliated to the GH23 family, which mostly include lysozymes. The protein encoded by ORF 138 belongs to the glycoside hydrolase family GH19, which comprises chitinases and lysozymes. Both enzyme

groups catalyze the hydrolysis of polysaccharides containing N-Acetylglucosamine, but the gene sequence does not discriminate between the two. No genes implicated in lysogeny establishment or control (integrase/excisionase, transposases, ParA/ParB genes, attachment sites, transcription repressor, etc.) were detected in the BONAISHI genome, even using prophage expert databases, which confirms that it is exclusively lytic. In addition, blastp searches using BONAISHI hosts genome did not find evidence of the presence of BONAISHI as prophage. Only three genes involved in DNA replication and recombination (ORF 225), nucleotide metabolism (ORF 192) and gene encoding for a phosphate starvation protein showed high score but the similarity was relatively low (<77%).

## Miscellaneous Proteins

We identified a phosphate starvation protein PhoH (ORF 261). This protein has been reported in many other *Myoviridae*. For therapeutic applications of this phage, we also looked for potential toxin encoding genes using BLASTP on 4 toxin expert databases. No homologs to currently known toxins were detected.

## Metagenomic Analysis

The global distribution of BONAISHI was investigated using marine viromes from various waters including coral mucus and the water from coral reefs (Table S2). The highest number of
<figure><img src="image_5.png" /></figure>

**FIGURE 3** | Genome map of BONAISHI. The 110 protein-coding genes are shown as colored blocks. Protein-coding genes that are transcribed in the forward and reverse direction are located on the external and internal circle, respectively. The functional assignments of the protein-coding genes are derived from blastp matches in NCBI nr and expert databases. We confidently discriminated proteins involved in general virion structure (light blue), head, and tail structure and assembly (dark blue), DNA replication, recombination and repair (red), nucleotide metabolism (orange), DNA transcription (pink), and host phage interactions (green). Genes marked "h.p." made strong matches to hypothetical protein sequences in the database with no identified functions (dark gray). The GC content of the genome sequence is indicated by the internal purple or green histograms.

recruited reads was from the virome collected at ALOHA station in the North Pacific subtropical gyre (CAM_SMPL_00823): BONAISHI recruited 1.04% of the reads with a mean identity of 58.23%. None of the recruited reads exceeded 85% identity. We also used prokaryotic metagenomes to test whether the BONAISHI genome is found as a prophage in bacterial hosts. No bacterial reads were recruited.

## Preliminary Therapy Assays of Diseased Symbiodinium

*Symbiodinium* sp. control culture showed optimal photosynthetic activity with the quantum yield Fv/Fm between 0.57 and 0.59 during the course of the experiment. As expected, inoculation with *V. coralliilyticus* caused a rapid photoinhibition with a

50% decrease in the quantum yield 60 min after inoculation (Figure 5). BONAISHI was able to significantly counteract the bacterial algicidal activity (t-test, *p* < 0.001) as, with both *V. coralliilyticus* and BONAISHI, the quantum yield was only 14% lower than the control 60 min after inoculation.

## DISCUSSION

In recent years, *Vibrio coralliilyticus* has been used as a model pathogen to gain insights into the establishment and propagation of coral diseases (Sussman et al., 2008; O'Santos et al., 2011; Garren et al., 2014; Pollock et al., 2015). The use of phages to control *V. coralliilyticus* has been reported recently
<table><caption>**TABLE 2** | Summary table of Vibrio phage BONAISHI predicted proteins that contained relevant annotation information as determined from significant BLASTP hits (e-value < e-3) against the GenBank non-redundant and CAZY databases.</caption><thead><tr><th></th><th>Predicted function</th><th>START</th><th>STOP</th><th>Size (n)</th><th>Best-hit</th><th>Classification</th><th>E-value</th><th>Score (bits)</th><th>N° accession</th></tr></thead><tbody><tr><td>ORF01c</td><td>hypothetical protein</td><td>653</td><td>1753</td><td>1100</td><td>Pseudocercospora musae</td><td>Fungi</td><td>2.00E-07</td><td>66.00</td><td>KNG49823.1</td></tr><tr><td>ORF02c</td><td>SbcC-like protein</td><td>1750</td><td>2184</td><td>435</td><td>Pseudomonas phage 201phi2-1</td><td>Myoviridae</td><td>1.00E-17</td><td>88.00</td><td>YP_001956973</td></tr><tr><td>ORF03</td><td>conserved hypothetical protein</td><td>2241</td><td>3119</td><td>879</td><td>Pseudomonas phage Phikz</td><td>Myoviridae</td><td>2.00E-32</td><td>130.00</td><td>NP_803730</td></tr><tr><td>ORF05</td><td>PaaR repeat-containing protein</td><td>3722</td><td>4018</td><td>297</td><td>Vibrio cholerae</td><td>Y-proteobacteria</td><td>3.00E-21</td><td>89.00</td><td>WP_057643503</td></tr><tr><td>ORF06c</td><td>conserved hypothetical protein</td><td>4055</td><td>4657</td><td>603</td><td>Pseudomonas phage Phikz</td><td>Myoviridae</td><td>1.00E-25</td><td>203.00</td><td>AAL83062</td></tr><tr><td>ORF07c</td><td>virion structural protein</td><td>4667</td><td>5962</td><td>1296</td><td>Pseudomonas phage KTN4</td><td>Myoviridae</td><td>7.00E-32</td><td>147.00</td><td>ANM44952.1</td></tr><tr><td>ORF10</td><td>ribonuclease H</td><td>6889</td><td>8271</td><td>1383</td><td>Pseudomonas phage 201phi2-1</td><td>Myoviridae</td><td>3.00E-25</td><td>117.00</td><td>YP_001956963</td></tr><tr><td>ORF12</td><td>uvsx protein</td><td>9269</td><td>10744</td><td>1476</td><td>Pseudomonas phage 201phi2-1</td><td>Myoviridae</td><td>1.00E-89</td><td>290.00</td><td>YP_001956960</td></tr><tr><td>ORF13</td><td>conserved hypothetical protein</td><td>10746</td><td>11108</td><td>363</td><td>Pseudomonas phage phiPA3</td><td>Myoviridae</td><td>3.00E-28</td><td>108.00</td><td>AEH03597</td></tr><tr><td>ORF15</td><td>virion structural protein</td><td>11460</td><td>12125</td><td>666</td><td>Pseudomonas phage phiPA3</td><td>Myoviridae</td><td>2.00E-27</td><td>113.00</td><td>AEH03595</td></tr><tr><td>ORF17</td><td>hypothetical protein</td><td>12558</td><td>13256</td><td>699</td><td>Erwinia phage Ea35-70</td><td>Myoviridae</td><td>8.00E-09</td><td>64.00</td><td>YP_009004948</td></tr><tr><td>ORF22c</td><td>conserved hypothetical protein</td><td>16270</td><td>18384</td><td>2115</td><td>Ralstonia phage RSL2</td><td>Myoviridae</td><td>3.00E-74</td><td>261.00</td><td>BAQ02568</td></tr><tr><td>ORF23c</td><td>glycoside hydrolase</td><td>18393</td><td>24110</td><td>5718</td><td>Ralstonia phage RP31</td><td>Myoviridae</td><td>2.00E-70</td><td>277.00</td><td>BAW19303.1</td></tr><tr><td>ORF24</td><td>RNA polymerase beta prime subunit</td><td>24183</td><td>25850</td><td>1668</td><td>Pseudomonas phage OBP</td><td>Myoviridae</td><td>1.00E-110</td><td>346.00</td><td>YP_004958184</td></tr><tr><td>ORF25</td><td>RNA polymerase beta subunit</td><td>25847</td><td>28741</td><td>2895</td><td>Ralstonia phage RSF1</td><td>Myoviridae</td><td>1.00E-102</td><td>358.00</td><td>BAS04832</td></tr><tr><td>ORF26</td><td>RNA polymerase beta subunit</td><td>29242</td><td>30369</td><td>1128</td><td>Pseudomonas phage 201phi2-1</td><td>Myoviridae</td><td>2.00E-61</td><td>223.00</td><td>YP_001956996</td></tr><tr><td>ORF29</td><td>conserved hypothetical protein</td><td>32135</td><td>32860</td><td>726</td><td>Pseudomonas phage Noxifer</td><td>Myoviridae</td><td>3.00E-21</td><td>110.00</td><td>ARV77361.1</td></tr><tr><td>ORF30</td><td>virion structural protein</td><td>32871</td><td>33641</td><td>771</td><td>Ralstonia phage RSF1</td><td>Myoviridae</td><td>3.00E-31</td><td>124.00</td><td>BAS04828</td></tr><tr><td>ORF36c</td><td>tail tube protein</td><td>36168</td><td>37040</td><td>873</td><td>Pseudomonas phage 201phi2-1</td><td>Myoviridae</td><td>2.00E-24</td><td>108.00</td><td>YP_001956757</td></tr><tr><td>ORF37c</td><td>tail sheath protein</td><td>37089</td><td>39137</td><td>2049</td><td>Erwinia phage Ea35-70</td><td>Myoviridae</td><td>2.00E-83</td><td>286.00</td><td>YP_009004971</td></tr><tr><td>ORF38</td><td>hypothetical protein</td><td>39208</td><td>40278</td><td>1071</td><td>Erwinia phage Ea35-70</td><td>Myoviridae</td><td>3.00E-09</td><td>67.00</td><td>YP_009004972</td></tr><tr><td>ORF39</td><td>virion structural protein</td><td>40288</td><td>42819</td><td>2532</td><td>Pseudomonas phage 201phi2-1</td><td>Myoviridae</td><td>2.00E-55</td><td>226.00</td><td>YP_001956754.1</td></tr><tr><td>ORF40</td><td>virion structural protein</td><td>42831</td><td>44483</td><td>1653</td><td>Pseudomonas phage 201phi2-1</td><td>Myoviridae</td><td>3.00E-27</td><td>125.00</td><td>YP_001956753</td></tr><tr><td>ORF41</td><td>terminase large subunit</td><td>44535</td><td>46691</td><td>2157</td><td>Pseudomonas phage 201phi2-1</td><td>Myoviridae</td><td>1.00E-142</td><td>442.00</td><td>YP_001956731</td></tr><tr><td>ORF43c</td><td>conserved hypothetical protein</td><td>47798</td><td>49027</td><td>1230</td><td>Pseudomonas phage Noxifer</td><td>Myoviridae</td><td>4.00E-37</td><td>164.00</td><td>ARV77197.1</td></tr><tr><td>ORF44</td><td>HD domain protein</td><td>49121</td><td>49720</td><td>600</td><td>Salmonella phage SPN3US</td><td>Myoviridae</td><td>3.00E-23</td><td>101.00</td><td>AEP84084</td></tr><tr><td>ORF49c</td><td>RNA-binding protein</td><td>52606</td><td>54225</td><td>1620</td><td>Actinomyces massiliensis</td><td>Actinobacteria</td><td>2.00E-73</td><td>251.00</td><td>WP_017194325</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 2 | Continued</caption><thead><tr><th></th><th>Predicted function</th><th>START</th><th>STOP</th><th>Size (n)</th><th>Best-hit</th><th>Classification</th><th>E-value</th><th>Score (bits)</th><th>N° accession</th></tr></thead><tbody><tr><td>ORF58c</td><td>conserved hypothetical protein</td><td>58317</td><td>59075</td><td>759</td><td>Erwinia phage Ea35-70</td><td>Myoviridae</td><td>1.00E-22</td><td>102.00</td><td>YP_009005001</td></tr><tr><td>ORF59</td><td>tubulin-like protein</td><td>59178</td><td>60167</td><td>990</td><td>Erwinia phage Ea35-70</td><td>Myoviridae</td><td>3.00E-21</td><td>101.00</td><td>YP_009005002</td></tr><tr><td>ORF68</td><td>hypothetical protein</td><td>64017</td><td>64424</td><td>408</td><td>Ralstonia phage RSL2</td><td>Myoviridae</td><td>7.00E-06</td><td>52.00</td><td>BAQ02532</td></tr><tr><td>ORF69</td><td>conserved hypothetical protein</td><td>64596</td><td>66782</td><td>2187</td><td>Erwinia phage Ea35-70</td><td>Myoviridae</td><td>1.00E-105</td><td>345.00</td><td>YP_009005012</td></tr><tr><td>ORF70c</td><td>hypothetical protein</td><td>66814</td><td>67875</td><td>1062</td><td>Ralstonia phage RSF1</td><td>Myoviridae</td><td>9.00E-08</td><td>63.00</td><td>BAS05022</td></tr><tr><td>ORF71</td><td>hypothetical protein</td><td>68148</td><td>70115</td><td>1967</td><td>γ-proteobacteria bacterium</td><td>γ-proteobacteria</td><td>1.00E-15</td><td>94.00</td><td>OUV32520.1</td></tr><tr><td>ORF72</td><td>RNA polymerase beta prime subunit</td><td>70204</td><td>71697</td><td>1494</td><td>Erwinia phage Ea35-70</td><td>Myoviridae</td><td>3.00E-44</td><td>171.00</td><td>YP_009005015</td></tr><tr><td>ORF73</td><td>conserved hypothetical protein</td><td>71795</td><td>72433</td><td>639</td><td>Erwinia phage Ea35-70</td><td>Myoviridae</td><td>9.00E-15</td><td>79.00</td><td>YP_009005021</td></tr><tr><td>ORF76</td><td>nuclease SbcC subunit</td><td>73326</td><td>74498</td><td>1173</td><td>Pseudomonas phage PhipA3</td><td>Myoviridae</td><td>1.00E-46</td><td>172.00</td><td>AEH03486</td></tr><tr><td>ORF77</td><td>conserved hypothetical protein</td><td>74495</td><td>75289</td><td>795</td><td>Erwinia phage PhiEaH1</td><td>Siphoviridae</td><td>9.00E-23</td><td>103.00</td><td>YP_009010069</td></tr><tr><td>ORF78</td><td>hypothetical protein</td><td>75301</td><td>75972</td><td>672</td><td>Uncultured bacterium</td><td>Bacteria</td><td>2.00E-06</td><td>56.00</td><td>EKD22589</td></tr><tr><td>ORF79</td><td>hypothetical protein</td><td>76106</td><td>77656</td><td>1551</td><td>Pseudomonas phage phiPA3</td><td>Myoviridae</td><td>9.00E-06</td><td>59.00</td><td>AEH03489</td></tr><tr><td>ORF80</td><td>conserved hypothetical protein</td><td>77691</td><td>79211</td><td>1521</td><td>Pseudomonas phage phiPA3</td><td>Myoviridae</td><td>1.00E-16</td><td>93.00</td><td>AEH03490</td></tr><tr><td>ORF81</td><td>hypothetical protein</td><td>79251</td><td>81056</td><td>1806</td><td>Gossypium arboreum</td><td>Magnoliopsida</td><td>1.00E-06</td><td>63.00</td><td>KHG21929</td></tr><tr><td>ORF82c</td><td>hypothetical protein</td><td>81109</td><td>81456</td><td>348</td><td>Pseudomonas phage PA7</td><td>Myoviridae</td><td>5.00E-09</td><td>59.00</td><td>AFO71119</td></tr><tr><td>ORF83</td><td>RNA polymerase beta subunit</td><td>81518</td><td>83584</td><td>2067</td><td>Pseudomonas phage Phabio</td><td>Myoviridae</td><td>8.00E-66</td><td>260.00</td><td>ARV76743.1</td></tr><tr><td>ORF84</td><td>RNA polymerase beta prime subunit</td><td>83584</td><td>85557</td><td>1974</td><td>Ralstonia phage RSF1</td><td>Myoviridae</td><td>1.00E-50</td><td>196.00</td><td>BAS05006</td></tr><tr><td>ORF85</td><td>helicase</td><td>85657</td><td>87198</td><td>1542</td><td>Pseudomonas phage OBP</td><td>Myoviridae</td><td>1.00E-40</td><td>176.00</td><td>AEV89521.1</td></tr><tr><td>ORF86</td><td>Clp protease subunit</td><td>87276</td><td>87848</td><td>573</td><td>Bacillus cereus</td><td>Bacilli</td><td>5.00E-05</td><td>52.00</td><td>WP_048520069</td></tr><tr><td>ORF87</td><td>ATP-dependent Clp protease proteolytic subunit</td><td>87848</td><td>88339</td><td>492</td><td>Dactylosporangium aurantiacum</td><td>Actinobacteria</td><td>1.00E-27</td><td>110.00</td><td>WP_033356707</td></tr><tr><td>ORF89c</td><td>conserved hypothetical protein</td><td>88833</td><td>90389</td><td>1557</td><td>Pseudomonas phage PA7</td><td>Myoviridae</td><td>2.00E-15</td><td>89.00</td><td>AFO71110</td></tr><tr><td>ORF92</td><td>RNA polymerase beta prime subunit</td><td>91656</td><td>92915</td><td>1260</td><td>Ralstonia phage RSF1</td><td>Myoviridae</td><td>8.00E-49</td><td>180.00</td><td>BAS04991</td></tr><tr><td>ORF96c</td><td>DNA polymerase polB</td><td>96090</td><td>97820</td><td>1731</td><td>Ralstonia phage RP12</td><td>Unclassified virus</td><td>6.00E-107</td><td>343.00</td><td>BAW19225.1</td></tr><tr><td>ORF97</td><td>virion structural protein</td><td>97895</td><td>99190</td><td>1296</td><td>Erwinia phage PhiEaH1</td><td>Siphoviridae</td><td>2.00E-23</td><td>111.00</td><td>YP_009010288</td></tr><tr><td>ORF99c</td><td>virion structural protein</td><td>101147</td><td>104005</td><td>2859</td><td>Pseudomonas phage 201phi2-1</td><td>Myoviridae</td><td>1.00E-93</td><td>327.00</td><td>YP_001956873</td></tr><tr><td>ORF100c</td><td>virion structural protein</td><td>104007</td><td>105122</td><td>1116</td><td>Ralstonia phage RSL2</td><td>Myoviridae</td><td>2.00E-49</td><td>180.00</td><td>BAQ02702</td></tr><tr><td>ORF101</td><td>capsid protein*</td><td>105166</td><td>106242</td><td>1077</td><td>Ralstonia phage RSF1</td><td>Myoviridae</td><td>3.00E-22</td><td>105.00</td><td>BAS04975</td></tr><tr><td>ORF102</td><td>virion structural protein</td><td>106257</td><td>107141</td><td>885</td><td>Ralstonia phage RSF1</td><td>Myoviridae</td><td>4.00E-07</td><td>60.00</td><td>BAS04974</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 2 | Continued</caption><thead><tr><th></th><th>Predicted function</th><th>START</th><th>STOP</th><th>Size (n)</th><th>Best-hit</th><th>Classification</th><th>E-value</th><th>Score (bits)</th><th>N° accession</th></tr></thead><tbody><tr><td>ORF104</td><td>hypothetical protein</td><td>107713</td><td>109107</td><td>1394</td><td>γ-proteobacteria bacterium</td><td>γ-proteobacteria</td><td>1.00E-12</td><td>83.00</td><td>OUV32343.1</td></tr><tr><td>ORF107</td><td>conserved hypothetical protein</td><td>111644</td><td>113257</td><td>1614</td><td>Ralstonia phage RSF1</td><td>Myoviridae</td><td>8.00E-22</td><td>108.00</td><td>BAS04969</td></tr><tr><td>ORF108</td><td>virion structural protein</td><td>113257</td><td>114459</td><td>1203</td><td>Pseudomonas phage PhiPA3</td><td>Myoviridae</td><td>8.00E-12</td><td>77.00</td><td>AEH03528</td></tr><tr><td>ORF110</td><td>virion structural protein</td><td>115149</td><td>116531</td><td>1383</td><td>Pseudomonas phage PhiPA3</td><td>Myoviridae</td><td>2.00E-28</td><td>126.00</td><td>AEH03530</td></tr><tr><td>ORF111c</td><td>helicase</td><td>116571</td><td>118181</td><td>1611</td><td>Pseudomonas phage 201phi2-1</td><td>Myoviridae</td><td>2.00E-48</td><td>1884.00</td><td>YP_001956921</td></tr><tr><td>ORF113</td><td>major capsid protein</td><td>118843</td><td>121035</td><td>2193</td><td>Erwinia phage Ea35-70</td><td>Myoviridae</td><td>2.00E-23</td><td>116.00</td><td>YP_009005109</td></tr><tr><td>ORF115</td><td>conserved hypothetical protein</td><td>122460</td><td>124031</td><td>1572</td><td>Ralstonia phage RSL2</td><td>Myoviridae</td><td>1.00E-46</td><td>179.00</td><td>BAQ02643</td></tr><tr><td>ORF118c</td><td>holliday-junction resolved</td><td>127427</td><td>127987</td><td>561</td><td>Pseudomonas phage Phabio</td><td>Myoviridae</td><td>4.00E-18</td><td>99.00</td><td>ARV76843.1</td></tr><tr><td>ORF119c</td><td>virion structural protein</td><td>128029</td><td>128865</td><td>837</td><td>Pseudomonas phage 201phi2-1</td><td>Myoviridae</td><td>8.00E-44</td><td>159.00</td><td>YP_001956947</td></tr><tr><td>ORF120c</td><td>virion structural protein</td><td>128878</td><td>130947</td><td>2070</td><td>Pseudomonas phage phiPA3</td><td>Myoviridae</td><td>7.00E-72</td><td>256.00</td><td>AEH03570</td></tr><tr><td>ORF121</td><td>virion structural protein</td><td>131049</td><td>133649</td><td>2601</td><td>Pseudomonas phage phabio</td><td>Myoviridae</td><td>4.00E-61</td><td>245.00</td><td>ARV76832.1</td></tr><tr><td>ORF133</td><td>hypothetical protein</td><td>139060</td><td>140100</td><td>1041</td><td>Vibrio tasmaniensis</td><td>γ-proteobacteria</td><td>2.00E-11</td><td>73.00</td><td>WP_017112059</td></tr><tr><td>ORF138</td><td>glycoside hydrolase</td><td>141805</td><td>142716</td><td>912</td><td>Aureimonas altamirensis</td><td>α-proteobacteria</td><td>2.00E-49</td><td>174.00</td><td>BAT26087</td></tr><tr><td>ORF141</td><td>hypothetical protein</td><td>145070</td><td>148201</td><td>3132</td><td>Psychromonas ingrahamii</td><td>γ-proteobacteria</td><td>1.00E-08</td><td>72.00</td><td>WP_011768462.1</td></tr><tr><td>ORF143</td><td>hypothetical protein</td><td>150013</td><td>151644</td><td>1631</td><td>Vibrio phage s4-7</td><td>Unclassified virus</td><td>2.00E-06</td><td>63.00</td><td>AOQ26845.1</td></tr><tr><td>ORF151</td><td>hypothetical protein</td><td>160192</td><td>161433</td><td>1242</td><td>Colwellia phage 9A</td><td>Siphoviridae</td><td>4.00E-08</td><td>66.00</td><td>YP_006489231</td></tr><tr><td>ORF159</td><td>hypothetical protein</td><td>168857</td><td>169801</td><td>945</td><td>Escherichia phage phAPEC8</td><td>Myoviridae</td><td>4.00E-10</td><td>68.00</td><td>YP_007348452</td></tr><tr><td>ORF163</td><td>hypothetical protein</td><td>172164</td><td>172847</td><td>684</td><td>Ruegeria halocynthiae</td><td>α-proteobactérie</td><td>1.00E-19</td><td>91.00</td><td>WP_037312174</td></tr><tr><td>ORF167</td><td>dead-like helicase</td><td>176707</td><td>178806</td><td>2100</td><td>Erwinia phage Ea35-70</td><td>Myoviridae</td><td>1.00E-114</td><td>367.00</td><td>YP_009004923</td></tr><tr><td>ORF172</td><td>nicotinamide-nucleotide adenylyltransferase</td><td>182827</td><td>184476</td><td>1650</td><td>Vibrio phage 11895-B1</td><td>Myoviridae</td><td>1.00E-106</td><td>331.00</td><td>YP_007673553</td></tr><tr><td>ORF176</td><td>nicotinamide-nucleotide adenylyltransferase</td><td>185996</td><td>187114</td><td>1119</td><td>Thiorhodococcus drewsii</td><td>γ-proteobacteria</td><td>3.00E-46</td><td>170.00</td><td>WP_007039048</td></tr><tr><td>ORF177</td><td>nicotinamide phosphoribosyltransferase</td><td>187169</td><td>188668</td><td>1500</td><td>Vibrio nigripulchritudo</td><td>γ-proteobacteria</td><td>3.00E-59</td><td>210.00</td><td>WP_022562194</td></tr><tr><td>ORF181</td><td>nicotinamide riboside transporter</td><td>190786</td><td>191502</td><td>717</td><td>Vibrio phage 11895-B1</td><td>Myoviridae</td><td>3.00E-76</td><td>238.00</td><td>YP_007673552</td></tr><tr><td>ORF189</td><td>hypothetical protein</td><td>196049</td><td>196543</td><td>717</td><td>Kaistia granuli</td><td>α-proteobacteria</td><td>1.00E-08</td><td>57.00</td><td>WP_018183972</td></tr><tr><td>ORF198</td><td>NAD-dependent DNA ligase</td><td>204242</td><td>206257</td><td>2016</td><td>Vibrio maritimus</td><td>γ-proteobacteria</td><td>0.00E+00</td><td>608.00</td><td>WP_042496716</td></tr><tr><td>ORF200</td><td>phosphatase</td><td>206698</td><td>207330</td><td>633</td><td>Vibrio phage VH7D</td><td>Myoviridae</td><td>6.00E-19</td><td>89.00</td><td>YP_009006310</td></tr><tr><td>ORF206</td><td>hypothetical protein</td><td>210803</td><td>211318</td><td>516</td><td>Psychromonas aquimarina</td><td>γ-proteobacteria</td><td>7.00E-34</td><td>126.00</td><td>WP_028862581</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 2 | Continued</caption><thead><tr><th></th><th>Predicted function</th><th>START</th><th>STOP</th><th>Size (n)</th><th>Best-hit</th><th>Classification</th><th>E-value</th><th>Score (bits)</th><th>N° accession</th></tr></thead><tbody><tr><td>ORF209</td><td>hypothetical protein</td><td>213895</td><td>214209</td><td>315</td><td>Enterovibrio calviensis</td><td>γ-proteobacteria</td><td>6.00E-07</td><td>52.00</td><td>WP_017007757</td></tr><tr><td>ORF212</td><td>transcriptional regulator</td><td>215240</td><td>215764</td><td>525</td><td>Leptolyngbya sp. PCC 7375</td><td>cyanobacteria</td><td>1.00E-05</td><td>53.00</td><td>EKV01169</td></tr><tr><td>ORF220</td><td>hypothetical protein</td><td>220875</td><td>221303</td><td>429</td><td>Shewanella sp. phage 1/4</td><td>Myoviridae</td><td>6.00E-11</td><td>64.00</td><td>YP_009100318</td></tr><tr><td>ORF227</td><td>hypothetical protein</td><td>224066</td><td>224998</td><td>933</td><td>Pseudomonas phage PhiPA3</td><td>Myoviridae</td><td>9.00E-06</td><td>57.00</td><td>AEH03433</td></tr><tr><td>ORF229</td><td>hypothetical protein</td><td>226086</td><td>227921</td><td>1835</td><td>Vibrio phage RYC</td><td>Unclassified virus</td><td>7.00E-70</td><td>243.00</td><td>BAV81012.1</td></tr><tr><td>ORF231</td><td>thymidylate kinase</td><td>228796</td><td>229458</td><td>663</td><td>Desulfosporosinus acidiphilus</td><td>Clostridia</td><td>8.00E-42</td><td>150.00</td><td>WP_014828008</td></tr><tr><td>ORF232</td><td>ribonucleotide-diphosphate reductase subunit alpha</td><td>229535</td><td>230446</td><td>912</td><td>Aeromonas molluscorum 848</td><td>γ-proteobacteria</td><td>1.00E-132</td><td>385.00</td><td>EOD53957</td></tr><tr><td>ORF233</td><td>ribonucleotide-diphosphate reductase subunit alpha</td><td>230844</td><td>232277</td><td>1434</td><td>Neisseria meningitidis</td><td>β-proteobacteria</td><td>0.00E+00</td><td>635.00</td><td>WP_049227356</td></tr><tr><td>ORF234</td><td>ribonucleotide-diphosphate reductase subunit beta</td><td>232355</td><td>233062</td><td>708</td><td>Thiomicrospira sp. Kp2</td><td>γ-proteobacteria</td><td>3.00E-96</td><td>294.00</td><td>WP_040727751</td></tr><tr><td>ORF235</td><td>HNH endonuclease</td><td>233205</td><td>233927</td><td>723</td><td>Bacillus pumilus</td><td>Bacilli</td><td>7.00E-32</td><td>123.00</td><td>WP_051149989</td></tr><tr><td>ORF236</td><td>ribonucleotide-diphosphate reductase subunit beta</td><td>234231</td><td>234617</td><td>387</td><td>Shigella sonnei</td><td>γ-proteobacteria</td><td>2.00E-36</td><td>129.00</td><td>CSE34793</td></tr><tr><td>ORF238</td><td>hypothetical protein</td><td>235508</td><td>237232</td><td>1725</td><td>Polyangium brachysporum</td><td>β-proteobacteria</td><td>2.00E-23</td><td>114.00</td><td>WP_047195109</td></tr><tr><td>ORF249c</td><td>hypothetical protein</td><td>242652</td><td>244886</td><td>2235</td><td>Vibrio mimicus</td><td>γ-proteobacteria</td><td>4.00E-08</td><td>67.00</td><td>WP_001015571</td></tr><tr><td>ORF253</td><td>phosphatase</td><td>246359</td><td>247015</td><td>657</td><td>Verrucomicrobium spinosum</td><td>verucomicrobia</td><td>2.00E-22</td><td>97.00</td><td>WP_009962909</td></tr><tr><td>ORF254</td><td>hypothetical protein</td><td>247012</td><td>247842</td><td>831</td><td>Shewanella sp.</td><td>γ-proteobacteria</td><td>6.00E-06</td><td>57.00</td><td>WP101034114.1</td></tr><tr><td>ORF255</td><td>ATP-binding protein</td><td>247870</td><td>248553</td><td>684</td><td>Vibrio phage RYC</td><td>Unclassified virus</td><td>1.00E-36</td><td>137.00</td><td>BAV81012.1</td></tr><tr><td>ORF256c</td><td>T5 A1-like protein</td><td>248608</td><td>250479</td><td>1872</td><td>Caulobacter phage phiCbK</td><td>Siphoviridae</td><td>3.00E-83</td><td>281.00</td><td>YP_006988022</td></tr><tr><td>ORF257c</td><td>hypothetical protein</td><td>250484</td><td>250927</td><td>444</td><td>Pseudoalteromonas (multispecies)</td><td>γ-proteobacteria</td><td>9.00E-19</td><td>86.00</td><td>WP_024591352</td></tr><tr><td>ORF258</td><td>conserved hypothetical protein</td><td>251087</td><td>252592</td><td>1506</td><td>Campylobacter phage CP30A</td><td>Myoviridae</td><td>3.00E-20</td><td>102.00</td><td>YP_006908082</td></tr><tr><td>ORF259</td><td>conserved hypothetical protein</td><td>252655</td><td>253968</td><td>1314</td><td>Campylobacter phage CP30A</td><td>Myoviridae</td><td>2.00E-18</td><td>96.00</td><td>YP_006908082</td></tr><tr><td>ORF261</td><td>phosphate starvation protein PhoH</td><td>254632</td><td>255627</td><td>996</td><td>Corynebacterium glucuronolyticum</td><td>Actinobacteria</td><td>6.00E-49</td><td>175.00</td><td>WP_005389286</td></tr><tr><td>ORF269</td><td>thymidylate synthase-complementing protein</td><td>260653</td><td>262095</td><td>1443</td><td>Parcubacteria</td><td>Parcubacteria</td><td>5.00E-62</td><td>217.00</td><td>KKR42866</td></tr><tr><td>ORF272</td><td>conserved hypothetical protein</td><td>263205</td><td>263723</td><td>519</td><td>Cronobacter phage vB_CsaM_GAP32</td><td>Myoviridae</td><td>9.00E-15</td><td>77.00</td><td>YP_006987447</td></tr><tr><td>ORF289</td><td>conserved hypothetical protein</td><td>273177</td><td>273917</td><td>740</td><td>Vibrio phage vB_VhaS-a</td><td>Unclassified virus</td><td>4.00E-18</td><td>100.00</td><td>ANO57550.1</td></tr><tr><td>ORF290</td><td>conserved hypothetical protein</td><td>273997</td><td>274665</td><td>668</td><td>Vibrio phage vB_VhaS-a</td><td>Unclassified virus</td><td>1.00E-14</td><td>88.00</td><td>ANO57549.1</td></tr><tr><td>ORF295c</td><td>hypothetical protein</td><td>276639</td><td>277370</td><td>731</td><td>Ralstonia phage RP12</td><td>Myoviridae</td><td>9.00E-10</td><td>73.00</td><td>BAW19047.1</td></tr><tr><td>ORF300</td><td>hypothetical tail protein</td><td>286648</td><td>287631</td><td>983</td><td>Pseudomonas phage Phabio</td><td>Myoviridae</td><td>7.00E-15</td><td>90.00</td><td>ARV76834.1</td></tr><tr><td>ORF301</td><td>virion structural protein</td><td>287726</td><td>288967</td><td>1241</td><td>Pseudomonas phage Noxifer</td><td>Myoviridae</td><td>3.00E-07</td><td>65.00</td><td>ARV77324.1</td></tr></tbody></table>
Tree scale: 0.1

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

<figure><img src="image_7.png" /><figcaption>**FIGURE 4** | Phylogenetic analysis of 63 Jumbo phages based on the amino acids sequence of the terminase large subunit. Sequences were aligned by Muscle and the tree was constructed by Maximum Likelihood method with a bootstrap of 1,000 using Mega 6.0 (Tamura et al., 2013) as in Yuan and Gao (2017). Bootstrap values higher than 0.5 are represented by black circles, ranging from 0.5 (smaller circle) to 1 (bigger circle).</figcaption></figure>

(Efrony et al., 2007; Cohen et al., 2013) but if phage therapy is to become a practical approach, fundamental knowledge on pathogen-virus interactions must be investigated in detail to evaluate, on the one hand, the therapeutic potential of the candidate phage and, on the other hand, the suitability of the host for phage therapy.

The detection of virus-derived genes in the *Vibrio coralliilyticus* P1 and YB1 genomes (Weynberg et al., 2015) shows that these pathogenic strains have interacted with phages during the course of their evolutionary history. Past interactions events with viruses can lead to the development of resistance mechanisms to escape phage infection, which may, in turn,
<figure><img src="image_8.png" /><figcaption>**FIGURE 5** | Quantum yield (Fv/Fm) of the photosystem II of control *Symbiodinium* sp. culture (black), *Symbiodinium* sp. inoculated with *V. coralliilyticus* YB1 pathogen (gray), and *Symbiodinium* sp. co-inoculated with *V. coralliilyticus* YB1 pathogen and *Vibrio* phage BONAISHI (hashed) after 60 min incubation. Results are expressed as the % of Fv/Fm in the control culture. As reported in previous study, the inoculation of *V. coralliilyticus* YB1 pathogen induced a rapid decline in *Symbiodinium* sp. photochemical efficiency. The addition of BONAISHI rapidly counteracted the impact of *V. coralliilyticus* YB1 on the efficiency of *Symbiodinium* photochemical activity.</figcaption></figure>

limit the application of phage therapy. We, however, did not detect any of the distinctive genetic signatures of viral resistance mechanisms, such as the insertion of short palindromic sequences (CRISPR, data not shown) in *V. coralliilyticus* P1 and YB1. Although other resistance mechanisms exist, the recurrent isolation of phages that infect *V. coralliilyticus* YB1 and/or P1 (Efrony et al., 2007; Cohen et al., 2013; Ramphul et al., 2017) supports the idea that *V. coralliilyticus* pathogens are permissive to viral infection and are suitable candidates for treatment by phage therapy.

The *Vibrio* phage BONAISHI isolated from coastal waters in the South China Sea is distinct from the known *V. coralliilyticus* phages YB2, YC, CKB-S1, CKB-S2, RYC (Efrony et al., 2007; Cohen et al., 2013; Ramphul et al., 2017). Although all these phages belong to the order of *Caudovirales* (tailed bacteriophages), BONAISHI has an unusually large genome, 303 kbp, rather than 11 kbp to 158 kbp. With such a large genome, BONAISHI is a novel jumbo phage (or giant phage). These are tailed phages with a dsDNA genome larger than 200 kb (Hendrix, 2009; Yuan and Gao, 2017). Jumbo phages have often been isolated in recent years and they mostly infect Gram-negative bacteria, including the genera *Synechococcus*, *Vibrio*, *Pseudomonas*, *Caulobacter*, *Erwinia*, and *Aeromonas* (Yuan and Gao, 2017). As observed for most jumbo phages, the large genome of BONAISHI is packaged in a large head connected to a long, contractile tail. The genome of jumbo phages typically includes core genes involved in virion structure and assembly, DNA replication, and nucleotide metabolism, including several genes encoding multisubunit RNAP. In BONAISHI, the core genes are scattered throughout the genome sequence and most of them share significant homology with other jumbo phages affiliated to the *Phikzirus* genus. BONAISHI is the

first marine representative of this divergent group within the *Myoviridae*. Interestingly, many Phikzviruses are considered to be promising biocontrol agents for plant-pathogenic bacteria (*Ralstonia solanacearum*, *Erwinia amylovora*) and some of them are already found in commercial products for phage therapy (Fujiwara et al., 2011; Bhunchoth et al., 2016).

The life history traits and genomic analysis showed that *Vibrio* phage BONAISHI is a good candidate for biological control of *V. coralliilyticus*. Firstly, BONAISHI appears to be structurally stable as it can withstand a wider range of pH (3–10) and temperature (4–45°C) than in the environment where it would be used. Second, it readily infects and lyses several pathogenic strains of *V. coralliilyticus* but no related species. Thirdly incubation experiments showed that the replication cycle is fast (latent period < 3 h). The presence of genes encoding virion-associated RNAP (ORFs 24, 25, 26, and 92) and early-expressed RNAP (ORFs 72, 83, and 84) in the BONAISHI genome may, at least partly, explain the rapid cycle (Ceyssens et al., 2014). In Phikzviruses, these two sets of RNAP may operate in concert during the replication cycle (Ceyssens et al., 2014; Yuan and Gao, 2017). The virion-associated RNAP may be injected into the host cell to start immediate gene expression whereas the early expressed RNAP may function during the middle and late phases of phage gene expression. The consecutive action of these enzymes, unique to Phikzviruses, confers the ability to produce viral progeny independently of the host transcription apparatus (Ceyssens et al., 2014; Yuan and Gao, 2017). In addition, the absence of detectable tRNA in BONAISHI genome suggests that it is well adapted to the translation machinery of its hosts, which is a critical process for efficient phage propagation. Finally, the genomic analysis did not identify any temperate phage hallmarks such as integration mediating enzymes, or genome architecture or sequence similarity with known temperate phages. Furthermore, the absence of homology between BONAISHI gene sequences and bacteria reads from coral metagenomes supports the idea that this phage does not integrate into the host genome. These results suggest very strongly that the *Vibrio* phage BONAISHI is a strictly lytic phage that is species specific and stable, although it appears to be relatively rare in the environment.

Another important issue for therapeutic applications of phages is to ensure that the candidate does not perform specialized or generalized transduction (Duckworth and Gulig, 2002). Given that BONAISHI appears to be strictly lytic based on the growth experiments and the genome analysis, it is unlikely that this candidate will perform specialized transduction of host DNA. Specialized transduction is restricted to temperate phages and occurs when the prophage is not cleanly excised during induction and includes the flanking bacterial genes which are then packaged in the viral progeny. Our candidate, however, may be able to perform generalized transduction. In this type of transduction, random segments of degraded host chromosome are mistakenly packaged instead of the phage DNA and may be transmitted by horizontal gene transfer. Phages that use a headful DNA packaging mechanism, such as many jumbo phages including BONAISHI, may be able to perform generalized transduction. However, it is, to the
best of our knowledge, impossible to predict the frequency of generalized transduction based purely on the genome analysis. For example, giant bacteriophages with similar headful DNA packaging mechanisms can have very different transduction rates as, for example, the T4 and T4G bacteriophages (Young et al., 1982; Young and Edlin, 1983). The ability of BONAISHI to perform generalized transduction would, therefore, require proper laboratory investigation.

An alternative to avoid potential issue with phage-mediated gene transfer is, rather than using whole bacetriophages, to use bacteriolytic proteins encoded by phages, among which the most notable are phage-encoded peptidoglycan hydrolases (PGH, see review by Roach and Debarbieux, 2017). PGHs, also called endolysins, degrade the cell peptidoglycan from within and contribute to the release of progeny and cell burst. A second type of PGH can be associated with the virion and initiate cell wall penetration through localized peptidoglycan or lipopolysaccharide degradation during the infection process. Both types of PGH are already used as bacteriocins in animal models of human infection and disease (see Roach and Debarbieux, 2017 and references therein). Jumbo phages typically encode more proteins for the lysis of the host cell wall including endolysin, glycoside hydrolase and chitinase, which are often bound to the virion than small genome phages (Yuan and Gao, 2017). In BONAISHI genome, we identified two glycoside hydrolases distantly related to known enzymes that belong to the families GH19 and GH23 using the expert database CAZY. Although the catalytic activities of these molecules cannot be determined based solely on the genome analysis, their overexpression and characterization might provide interesting tools for controlling *V. coralliilyticus* infection.

A preliminary assay suggests that BONAISHI is a promising candidate for treating *V. coralliilyticus* infection. Studies investigating the action of *V. coralliilyticus* on coral symbionts showed that photosynthesis was inactivated by the expression of a Zn-metalloprotease (Sussman et al., 2009). Our experiments on *Symbiodinium* cultures infected by *V. coralliilyticus* showed that BONAISHI phage treatment was effective: adding the phage to the infected cultures rapidly reduced *Symbiodinium* PSII inactivation. As reported in previous studies, phage

addition probably lysed the bacterial pathogen, stopping Zn-metalloprotease production and further damage to *Symbiodinium* sp. cells (Cohen et al., 2013). Future studies should now focus on the effectiveness of the treatment either under realistic field conditions or in mesocosms to start including bacteriophages (and/or derived compounds) in an integrated management program to mitigate the damage caused by the infectious agents responsible for coral diseases. We recommend genome sequencing and analysis of any future phage candidate as a prerequisite to any field test to ensure safe environmental applications as this provides essential information on the phage replication cycle and host-virus interactions.

## AUTHOR CONTRIBUTIONS

A-CB, YB, and TB, designed the study. LJ, JM, and A-CB performed the experiments and analyzed the results. LJ, SH, CD, and EC performed the bioinformatics analyses. CF-P provided and helped with the diseased *Symbiodinium* cultures. LJ and A-CB wrote the manuscript.

## ACKNOWLEDGMENTS

This work was supported jointly by the EC2CO PATRICIA Project, the TOTAL Foundation, and the ANR CALYPSO (ANR-15-CE01-0009). We would like to thank Andrew Millard for discussion on phage transduction, Simon Roux for his help with metagenomic analysis, Gurvan Michel for his assistance with CAZYmes annotation, Pei Ge for her technical assistance, and Tony Tebby for the manuscript editing. We would also like thank the three reviewers for their constructive comments on a previous version of this manuscript. We are also very grateful to Michel Galey, Alexandre Portier, and all the staff from Whale Island Resort for their hospitality and help during our stay.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.02501/full#supplementary-material

## REFERENCES

* Ackermann, H.-W., and Heldal, M. (2010). "Basic electron microscopy of aquatic viruses," in *Manual of Aquatic Viral Ecology*, eds S. W. Wilhelm, M. G. Weinbauer, and C. A. Suttle (Waco, TX: ASLO), 182-192.

* Atad, I., Zvuloni, A., Loya, Y., and Rosenberg, E. (2012). Phage therapy of the white plague-like disease of Favia favus in the Red Sea. *Coral Reefs* 31, 665-670. doi: 10.1007/s00338-012-0900-5

* Aziz, R. K., Bartels, D., Best, A. A., Dejongh, M., Disz, T., Edwards, R. A., et al. (2008). The RAST server: rapid annotations using subsystems technology. *BMC Genomics* 9:75. doi: 10.1186/1471-2164-9-75

* Bellwood, D. R., Hughes, T. P., and Folke, C., Nyström, M. (2004). Confronting the coral reef crisis. *Nature* 429, 827-833. doi: 10.1038/nature02691

* Ben-Haim, Y., Thompson, F. L., Thompson, C. C., Cnockaert, M. C., Hoste, B., Swings, J., et al. (2003). *Vibrio coralliilyticus* sp. nov. a temperature-dependent pathogen of the coral Pocillopora damicornis. *Int. J. Syst. Evol. Microbiol.* 53, 309-315. doi: 10.1099/ijis.0.02402-0

* Besemer, J., Lomsadze, A., and Borodovsky, M. (2001). GeneMarkS: a self-training method for prediction of gene starts in microbial genomes. Implications for finding sequence motifs in regulatory region. *Nucleic Acids Res.* 29, 2607-2618. doi: 10.1093/nar/29.12.2607

* Bhunchoth, A., Blanc-Mathieu, R., Mihara, T., Nishimura, Y., Askora, A., Phironrit, N., et al. (2016). Two asian jumbo phages, φRSL2 and φRSF1, infect *Ralstonia solanacearum* and show common features of φKZ-related phages. *Virology* 494, 56-66. doi: 10.1016/j.virol.2016.03.028

* Brussaard, C. P. D. (2004). Optimization of procedures for counting viruses by flow cytometry. *Appl Environ Microbiol* 70, 1506-1513. doi: 10.1128/AEM.70.3.1506-1513.2004

* Brussaard, C. P. D., Baudoux, A.-C., and Rodriguez-Varela, F. (2016). "Marine viruses," in *The Marine Microbiome-an Untold Resource of Biodiversity and Biotechnological Potential*. eds. L. J. Stal and S. M. Cretoiu (Springer International Publishing), 305-32.

* Carver, T., Thomson, N., Bleasby, A., Berriman, M., and Parkhill, J. (2009). DNAPlotter: circular and linear interactive genome visualization. *Bioinformatics* 25, 119-120. doi: 10.1093/bioinformatics/btn578
Ceyssens, P. J., Minakhin, L., Van den Bossche, A., Yakunina, M., Klimuk, E., Blasdel, B., et al. (2014). Development of giant bacteriophage phi KZ is independent of the host transcription apparatus. *J. Virol.* 88:105010. doi: 10.1128/JVI.01347-14

Chen, L. H., Yang, J., Yu, J., Yao, Z. J., Sun, L. L., Shen, Y., et al. (2005). VFDB: a reference database for bacterial virulence factors. *Nucleic Acids Res.* 33:D325–D328. doi: 10.1093/nar/gki008

Chin, C. S., Alexander, D. H., Marks, P., Klammer, A. A., Drake, J., Heiner, C., et al. (2013). Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. *Nat Methods* 10, 563–569 doi: 10.1038/nmeth.2474

Choo, S. W., Heydari, H., Tan, T. K., Siow, C. C., Beh, C. Y., Wee, W. Y., et al. (2014) VibrioBase: a model for next-generation genome and annotation database development. *Sci. World J.* 2014:569324. doi: 10.1155/2014/569324

Cohen, Y., Joseph Pollock, F., Rosenberg, E., and Bourne, D. G. (2013). Phage therapy treatment of the coral pathogen *Vibrio coralliilyticus*. *Microbiol. Open* 2, 64–74. doi: 10.1002/mbo3.52

Connell, J. H. (1978). Diversity in tropical rain forests and coral reefs. *Science* 199, 1302–1310. doi: 10.1126/science.199.4335.1302

Cooper, E. L., Hirabayashi, K., Strychar, K. B., and Sammarco, P. W. (2014). Corals and their potential applications to integrative medicine. *Evid. Based Complement. Alternat. Med.* 2014:184959. doi: 10.1155/2014/184959

Dang, V. T., Howard-Varona, C., Schwenck, S., and Sullivan, M. B. (2015). Variably lytic infection dynamics of large B acteroidetes podovirus ph38: 1 against two *Cellulophaga baltica* host strains. *Environ. Microbiol.* 17, 4659–4671. doi: 10.1111/1462-2920.13009

Delcher, A. L., Harmon, D., Kasif, S., White, O., and Salzberg, S. L. (1999). Improved microbial gene identification with GLIMMER. *Nucleic Acids Res.* 27, 4636–4641. doi: 10.1093/nar/27.23.4636

d'Herelle, F. (1926). *The Bacteriophage and Its Behavior*. Baltimore, MD: Williams & Wilkins, 490–497.

Dinsdale, E. A., Pantos, O., Smriga, S., Edwards, R. A., Angly, F., Wegley, L., et al. (2008). Microbial ecology of four coral atolls in the Northern Line Islands. *PLoS ONE* 3:e1584. doi: 10.1371/journal.pone.0001584

Duckworth, D. H. (1976). Who discovered bacteriophage? *Bacteriol. Rev.* 40:793.

Duckworth, D. H., and Gulig, P. A. (2002). Bacteriophages. *BioDrugs* 16, 57–62. doi: 10.2165/00063030-200216010-00006

Efrony, R., Atad, I., and Rosenberg, E. (2009). Phage therapy of coral white plague disease: properties of phage BA3. *Curr. Microbiol.* 58, 139–145. doi: 10.1007/s00284-008-9290-x

Efrony, R., Loya, Y., Bacharach, E., and Rosenberg, E. (2007). Phage therapy of coral disease. *Coral Reefs* 26, 7–13. doi: 10.1007/s00338-006-0170-1

Fujiwara, A., Fujisawa, M., Hamasaki, R., Kawasaki, T., Fujie, M., Yamada, T. (2011). Biocontrol of Ralstonia solanacearum by treatment with lytic bacteriophages. *Appl. Environ. Microbiol.* 77, 4155–4162. doi: 10.1128/AEM.02847-10

Garren, M., Son, K., Raina, J.-B., Rusconi, R., Menolascina, F., Shapiro, O. H., et al. (2014). A bacterial pathogen uses dimethylsulfoniopropionate as a cue to target heat-stressed corals. *ISME J.* 8, 999–1007. doi: 10.1038/ismej.2013.210

Goiran, C., Al-Moghrabi, S., Allemand, D., and Jaubert, J. (1996). Inorganic carbon uptake for photosynthesis by the symbiotic coral/dinoflagellate association I. Photosynthetic performances of symbionts and dependence on sea water bicarbonate. *J. Exp. Mar. Biol. Ecol.* 199, 207–225. doi: 10.1016/0022-0981(95)00201-4

Grissa, I., Vergnaud, G., and Pourcel, C. (2007). CRISPRfinder: A web tool to identify clustered regularly interspaced short palindromic repeats. *Nucleic Acids Res.* 35(Suppl. 2), W52–7. doi: 10.1093/nar/gkm360

Guillard, R., and Ryther, J. (1962). Studies of marine planktonic diatoms. I. *Cyclotella nana* Hustedt, and *Detonula confervacea*. *Can. J. Microbiol.* 8, 229–239. doi: 10.1139/m62-029

Hall, T. A. (1999). BioEdit: a user-friendly biological sequence alignment editor and analysis. *Nucleic Acids Symp. Ser.* 41, 95–98.

Hendrix, R. W. (2009). Jumbo Bacteriophages. *Curr. Top. Microbiol. Immunol.* 328, 229–240. doi: 10.1007/978-3-540-68618-7_7

Holmfeld, K., Howard-Varona, C., Solonenko, N., and Sullivan, M. B. (2014). Contrasting genomic patterns and infection strategies of two co-existing Bacteroidetes podovirus genera. *Environ. Microbiol.* 16, 2501–13. doi: 10.1111/1462-2920.12391

Howard-Varona, C., Hargreaves, K. R., Abedon, S. T., Sullivan, M. B. (2017). Lysogeny in nature: mechanisms, impact and ecology of temperate phages. *ISME J.* 11, 1511–1520. doi: 10.1038/ismej.2017.16

Hughes, T. P., Baird, A. H., Bellwood, D. R., Card, M., Connolly, S. R., Folke, C., et al. (2003). Climate change, human impacts, and the resilience of coral reefs. *Science* 301, 929–33. doi: 10.1126/science.1085046

Jackson, J. B. C., Kirby, M. X., Berger, W. H., Bjorndal, K. A., Botsford, L. W., Bourque, B. J., et al. (2001). Historical overfishing and the recent collapse of coastal ecosystems. *Science* 293, 629–637. doi: 10.1126/science.1059199

Jiang, S. C., and Paul, J. H. (1996). Occurrence of lysogenic bacteria in marine microbial communities as determined by prophage induction. *Mar. Ecol. Prog. Ser.* 142, 27–38. doi: 10.3354/meps142027

Jones, P., Binns, D., Chang, H.-Y., Fraser, M., Li, W., McAnulla, C., et al. (2014). InterProScan 5: genome-sclae protein function classification. *Bioinformatics* 30, 1236–1240. doi: 10.1093/bioinformatics/btu031

Laslett, D., and Canback, B. (2004) ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences. *Nucleic Acids Res* 32, 11–16. doi: 10.1093/nar/gkh152

Leplae, R., Lima-Mendez, G., Toussaint, A. (2006). A first global analysis of plasmid encoded proteins in the ACLAME database. *FEMS Microbiol. Rev.* 30, 980–994. doi: 10.1111/j.1574-6976.2006.00044.x

Lim, E., Pon, A., Djoumbou, Y., Knox, C., Shrivastava, S., Guo, A. C., et al. (2009). T3DB: a comprehensively annotated database of common toxins and their targets. *Nucleic Acids Res.* 38(Suppl. 1), D781–D786. doi: 10.1093/nar/gkp934

Lobocka, M., Hejnowicz, M. S., Gagala, U., Weber-Dabrowska, B., Wegrzyn, G., Dadlez, M. (2014). "The first step to bacteriophage therapy: how to choose the correct phage," in *Phage Therapy: Current Research and Applications*, eds J. Borysowski, R. Miedzybrodsky, and A. Gorski (Norfolk, UK: Caister Academic Press), 23–67.

Lowe, T. M., and Eddy, S. R. (1997). rRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. *Nucleic Acids Res.* 25, 955–964.

Mera, H., and Bourne, D. G. (2018). Disentangling causation: complex roles of coral associated microorganisms in disease. *Environ. Microbiol.* 20, 431–449. doi: 10.1111/1462-2920.13958

Middelboe, J., Chan, A. M., Bertelsen, S. K. (2010) "Isolation and life cycle characterization of lytic viruses infecting heterotrophic bacteria and cyanobacteria," in *Manual of Aquatic Viral Ecology*, eds S. W. Wilhelm, and M. G. Weinbauer, and C. A. Suttle (Waco, TX: ASLO), 118–133.

Moberg, F., and Folke, C. (1999). Ecological goods and services of coral reef ecosystems. *Ecol. Econ.* 29, 215–233. doi: 10.1016/S0921-8009(99)00009-9

Moura, A., Soares, M., Pereira, C., Leitão, N., Henriques, I., Correia, A. (2009). INTEGRALL: a database and search engine for integrons, integrases and gene cassettes. *Bioinformatics* 25, 1096–1098. doi: 10.1093/bioinformatics/btp105

O'Santos, E., Alves, N., Dias, G. M., Mazotto, A. M., Vermelho, A., Vora, G. J., et al. (2011). Genomic and proteomic analyses of the coral pathogen Vibrio coralliilyticus reveal a diverse virulence repertoire. *ISME J.* 5, 1471–83. doi: 10.1038/ismej.2011.19

Paez-Espino, P., Chen, I.-M. A., Palaniappan, K., Ratner, A., Chu, K., Szeto, E., et al. (2017). IMG/VR: a database of cultured and uncultured DNA Viruses and retroviruses. *Nucleic Acids Res.* 45, D457–D465. doi: 10.1093/nar/gkw1030

Pandolfi, J. M., Bradbury, R. H., Sala, E., Hughes, T. P., Bjorndal, K. A., Cooke, R. G., et al. (2003). Global trajectories of the long-term decline of coral reef ecosystems. *Science* 301, 955–958. doi: 10.1126/science.1085706

Pollock, F. J., Krediet, C. J., Garren, M., Stocker, R., Winn, K., Wilson, B., et al. (2015). Visualization of coral host-pathogen interactions using a stable GFP-labeled Vibrio coralliilyticus strain. *Coral Reefs* 34, 655–662. doi: 10.1007/s00338-015-1273-3

Ramphul, C., Estela, B., Dohra, H., Suzuki, T., Yoshimatsu, K., Yoshinaga, K., et al. (2017). Marine genomics genome analysis of three novel lytic *Vibrio coralliilyticus* phages isolated from seawater, Okinawa, Japan. *Mar. Genomics* 35, 69–75. doi: 10.1016/j.margen.2017.06.005

Roach, D. R., and Debarbieux, L. (2017). Phage therapy: awakening a sleeping giant. *Emerg. Topics Life Sci.* 1, 93–103. doi: 10.1042/ETLS20170002
Roux, S., Faubladier, M., Mahul, A., Paulhe, N., Bernard, A., et al. (2011) Metavir: a web server dedicated to virome analysis. *Bioinformatics* 27, 3074–3075.
doi: 10.1093/bioinformatics/btr519

Rutherford, K., Parkhill, J., Crook, J., Horsnell, T., Rice, P., Rajandream, M. A., et al. (2000) Artemis: sequence visualization and annotation. *Bioinformatics* 16, 944–945.
doi: 10.1093/bioinformatics/16.10.944

Seshadri, R., Kravitz, S. A., Smarr, L., Gilna, P., Frazier, M. (2007) CAMERA: a community resource for metagenomics. *PLoS Biol* 5:e75.
doi: 10.1371/journal.pbio.0050075

Sussman, M., Mieog, J. C., Doyle, J., Victor, S., Willis, B. L., Bourne, D. G. (2009). Vibrio zinc-metalloprotease causes photoinactivation of coral endosymbionts and coral tissue lesions. *PLoS ONE* 4:e4511. doi: 10.1371/journal.pone.0004511

Sussman, M., Willis, B. L., Victor, S., Bourne, D. G. (2008). Coral pathogens identified for White Syndrome (WS) epizootics in the Indo-Pacific. *PLoS ONE* 3:e2393. doi: 10.1371/journal.pone.0002393

Tamura, K., Stecher, G., Peterson, D., Filipski, A., and Kumar, S. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0. *Mol. Biol. Evol.* 30, 2725–2729. doi: 10.1093/molbev/mst197

Teplitski, M., and Ritchie, K. (2009). How feasible is the biological control of coral diseases? *Trends Ecol. Evol.* 24, 378–385. doi: 10.1016/j.tree.2009.02.008

Ushijima, B., Videau, P., Burger, A. H., Shore-Maggio, A., Runyon, C. M., Sudek, M., et al. (2014). Vibrio coralliilyticus strain OCN008 is an etiological agent of acute Montipora white syndrome. *Appl. Environ. Microbiol.* 80, 2102–2109.
doi: 10.1128/AEM.03463-13

Weynberg, K. D., Voolstra, C. R., Neave, M. J., Buerger, P., van Oppen, M. J. H. (2015) From cholera to corals: viruses as drivers of virulence in a major coral bacterial pathogen. *Sci. Rep.* 5:17889. doi: 10.1038/srep17889

Williamson, S. J., and Paul, J. H. (2006). Environmental factors that influence the transition from lysogenic to lytic existence in the φHSIC/Listonella pelagia marine phage-host system. *Microbial Ecol.* 52, 217–225.
doi: 10.1007/s00248-006-9113-1

Wishart, D., Arndt, D., Pon, A., Sajed, T., Guo, A. C., Djoembou, Y., et al. (2015). T3DB: the toxic exposome database. *Nucleic Acids Res.* 43, D928–D934.
doi: 10.1093/nar/gku1004

Young, K. K., Edlin, G. (1983). Physical and genetic analysis of bacteriophage T4: generalized transduction. *Mol. Gen. Genet.* 192, 241–246.
doi: 10.1007/BF00327673

Young, K. K., Edlin, G., and Wilson, G. G. (1982). Genetic analysis of bacteriophage T4: transducing bacteriophages. *J. Virol.* 41, 345–347.

Yuan, Y., and Gao, M. (2017). Jumbo bacteriophages: an overview. *Front. Microbiol.* 8, 1–9. doi: 10.3389/fmicb.2017.00403

Zhou, C. E., Smith, J., Lam, M., Zemla, A., Dyer, M. D., and Slezak, T. (2006). MvIRDB—a microbial database of protein toxins, virulence factors and antibiotic resistance genes for bio-defence applications. *Nucleic Acids Res.* 35(suppl. 1), D391–D394. doi: 10.1093/nar/gkl791

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2018 Jacquemot, Bettarel, Monjol, Corre, Halary, Desnues, Bouvier, Ferrier-Pagès and Baudoux. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.