# Phenotypic and genetic predictors of pathogenicity and virulence in Flavobacterium psychrophilum

Krister Sundell, Lotta Landor, Pierre Nicolas, Jóhanna Jørgensen, Daniel Castillo, Mathias Middelboe, Inger Dalsgaard, Valentina Laura Donati, Lone Madsen, Tom Wiklund

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

► To cite this version:

Krister Sundell, Lotta Landor, Pierre Nicolas, Jóhanna Jørgensen, Daniel Castillo, et al.. Phenotypic and genetic predictors of pathogenicity and virulence in Flavobacterium psychrophilum. Frontiers in Microbiology, 2019, 10 (1711), pp.1-14. 10.3389/fmicb.2019.01711. hal-02618584

HAL Id: hal-02618584

https://hal.inrae.fr/hal-02618584v1

Submitted on 25 May 2020

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_2.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_3.png" /></figure>

Check for
updates

# Phenotypic and Genetic Predictors of Pathogenicity and Virulence in *Flavobacterium psychrophilum*

Krister Sundell<fnref n="1" />*, Lotta Landor<fnref n="1" />, Pierre Nicolas², Jóhanna Jørgensen³, Daniel Castillo³, Mathias Middelboe³, Inger Dalsgaard⁴, Valentina Laura Donati⁴, Lone Madsen⁴ and Tom Wiklund<fnref n="1" />

OPEN ACCESS

**Edited by:**
Yuji Morita,
Meiji Pharmaceutical University, Japan

**Reviewed by:**
Robert Söderlund,
National Veterinary Institute, Sweden
Thomas P. Loch,
Michigan State University,
United States

*Correspondence:
Krister Sundell
ksundell@abo.fi

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

**Received:** 23 April 2019

**Accepted:** 11 July 2019

**Published:** 24 July 2019

**Citation:**
Sundell K, Landor L, Nicolas P, Jørgensen J, Castillo D, Middelboe M, Dalsgaard I, Donati VL, Madsen L and Wiklund T (2019) *Phenotypic and Genetic Predictors of Pathogenicity and Virulence in Flavobacterium psychrophilum*. *Front. Microbiol.* 10:1711. doi: 10.3389/fmicb.2019.01711

*Flavobacterium psychrophilum* causes bacterial cold-water disease (BCWD) in farmed rainbow trout (*Oncorhynchus mykiss*), with the multilocus sequence typing (MLST) clonal complex (CC) CC-ST10 accounting for the majority of outbreaks globally. The development of alternative strategies to antibiotic treatment of BCWD using bacteriophage-based control of *F. psychrophilum*, or virulence factors as targets for therapy, requires knowledge of the phage-sensitivity of outbreak strains and of universal traits contributing to their pathogenicity. To examine the association between virulence and both genetic (MLST sequence type (ST) and PCR-serotype) and phenotypic characteristics (adherence, antibiotic resistance, colony spreading motility, hemolytic and proteolytic activity), the median lethal dose (LD₅₀) of 26 geographically disparate *F. psychrophilum* isolates was determined in rainbow trout. Furthermore, the *in vitro* sensitivity of the isolates against five bacteriophages was determined by the efficiency of plating (EOP). The tested *F. psychrophilum* isolates were mainly represented by CC-ST10 genotypes (22 out of 26) and showed up to 3-log differences in LD₅₀ (8.9 × 10³ to 3.1 × 10⁶ CFU). No association between MLST ST and virulence was found because of a high variation in LD₅₀ within STs. All identified serotypes (0, 1, and 2) were pathogenic, but ten most virulent isolates belonged to serotype 1 or 2. Isolates of high (LD₅₀ < 10⁵ CFU), moderate (LD₅₀ = 10⁵–10⁶ CFU), and weak (LD₅₀ > 10⁶ CFU) virulence were similar in phenotypic characteristics *in vitro*. However, the only non-virulent CC-ST10 isolate was deficient in spreading motility and proteolytic activity, indicating that the characteristics are required for pathogenicity in *F. psychrophilum*. Univariate correlation studies found only non-significant associations between LD₅₀ and the measured phenotypic characteristics, and the multivariable analysis did neither reveal any significant predictors of virulence. The majority of isolates (16 out of 26) were sensitive to at least four bacteriophages, with up to a 6-log variation in the EOP. Most CC-ST10 isolates (16 out of 22) were sensitive to the examined phages, including 5

<footnote n="1">¹ Laboratory of Aquatic Pathobiology, Environmental and Marine Biology, Åbo Akademi University, Turku, Finland, ² Unité Mathématiques et Informatique Appliquées du Génome à l'Environnement (MalAGE), Institut National de la Recherche Agronomique, Université Paris-Saclay, Jouy-en-Josas, France, ³ Marine Biological Section, Department of Biology, University of Copenhagen, Helsingør, Denmark, ⁴ Unit for Fish and Shellfish Diseases, National Institute of Aquatic Resources, Technical University of Denmark, Kongens Lyngby, Denmark</footnote>
out of the 7 most virulent isolates represented by prevalent and antibiotic-resistant STs. Our findings suggest that control of BCWD using lytic phages or interventions targeting shared characteristics of pathogenic *F. psychrophilum* strains should be further explored.

**Keywords:** *Flavobacterium psychrophilum*, bacterial cold-water disease, rainbow trout (*Oncorhynchus mykiss*), multilocus sequence typing (MLST), bacteriophage, LD₅₀, pathogenicity, virulence

INTRODUCTION

*Flavobacterium psychrophilum* is the etiological agent of Bacterial Cold-Water Disease (BCWD) causing devastating losses to farmed salmonid fish particularly during early freshwater growth stages. In the Baltic Sea region in Northern Europe, the pathogen has also been isolated several times from diseased juvenile salmonid fish in brackish water fish farms (Madetoja et al., 2002; Nilsen et al., 2011) and from brood stock carriers returning from the sea to spawn (Ekman et al., 1999). However, the largest problems caused by *F. psychrophilum* are still associated with early life stages of rainbow trout (*Oncorhynchus mykiss*) fry in hatcheries where Rainbow Trout Fry Syndrome (RTFS), here used interchangeably with BCWD, can result in high mortalities if left untreated. Today BCWD outbreaks in farmed fish are controlled by antibiotics administered through feed, which may lead to the dispersal of antibiotics in the environment. The lack of efficient preventative measures and obstacles in developing protective vaccination strategies against *F. psychrophilum* (Gomez et al., 2014) may in part be due to the genetic and serological diversity within the species and the insufficient knowledge of specific virulence factors. Virulence factors aid pathogenic bacteria to adhere to and invade the host, cause disease, and to evade host defense mechanisms. To be able to design alternative control strategies that inhibit the bacterium from causing disease, understanding which bacterial characteristics contribute to virulence in a susceptible host is of crucial importance.

Several putative virulence genes have been identified in genome analyses of *F. psychrophilum* (Duchaud et al., 2007, 2018; Castillo et al., 2016), but the role played by different factors in its pathogenicity in rainbow trout is still poorly understood. Attempts to link serologic (Lorenzen and Olesen, 1997; Dalsgaard and Madsen, 2000) and phenotypic characteristics, such as adhesion to host tissue (Nematollahi et al., 2003), hemolytic activity (Högfors-Rönnholm and Wiklund, 2010a) and protease production (Bertolini et al., 1994; Dalsgaard and Madsen, 2000; Ostland et al., 2000) to *F. psychrophilum* virulence have shown varying results. Studies with *F. psychrophilum* mutants have shown a connection between attenuated virulence and deficiency in motility and secretion of proteases and adhesins (Pérez-Pascual et al., 2015, 2017). However, such studies usually involve evaluation of a single or a few virulence factors at a time, undermining the fact that virulence could be the result of their combined effects. To date, only a few studies (Rangdale, 1995; Garcia et al., 2000; Madsen and Dalsgaard, 2000; Madetoja et al., 2002; Jarau et al., 2018) have evaluated the virulence of several *F. psychrophilum* isolates in rainbow

trout, at least in part due to the lack of a challenge model that mimics natural disease. Waterborne infection models with *F. psychrophilum* have met with limited success in consistently inducing mortality in rainbow trout unless scarification (Madsen and Dalsgaard, 1999; Madetoja et al., 2000; Long et al., 2014) or a susceptible isogenic rainbow trout line (Pérez-Pascual et al., 2017) have been used. However, strain-dependent mortality rates can reproducibly be induced when *F. psychrophilum* is injected subcutaneously (Madetoja et al., 2002), intraperitoneally (Rangdale, 1995; Madsen and Dalsgaard, 1999, 2000; Jarau et al., 2018) or intramuscularly (Garcia et al., 2000; Madetoja et al., 2002) into rainbow trout. Despite the documented variation in virulence between isolates, reliable predictors of *F. psychrophilum* virulence have not yet been identified.

Individual *F. psychrophilum* isolates also vary in susceptibility to natural killers of bacteria, the bacteriophages (Stenholm et al., 2008). Bacteriophages with a strong lytic potential against *F. psychrophilum* have been discovered (Stenholm et al., 2008; Castillo et al., 2012) and phage-therapy has been proposed as an alternative strategy to antibiotic treatment for use in fish farms. A reemerging challenge for the use of bacteriophages to prevent or control *F. psychrophilum* in fish farm environments is the high degree of diversity found both within infective populations of the pathogen (Sundell et al., 2013) and phages (Stenholm et al., 2008; Castillo and Middelboe, 2016). For phage-therapy to be useful in aquaculture applications, more information is needed concerning the efficiency of phages against recurring virulent variants of *F. psychrophilum*.

Several epidemiological studies using a standardized multilocus sequence typing (MLST) approach for genetic typing of *F. psychrophilum* from wide geographical areas have proposed an epidemic population structure (Nicolas et al., 2008; Siekoula-Nguedia et al., 2012; Avendaño-Herrera et al., 2014; Nilsen et al., 2014), whereby one or several successful genotypes clonally expand and diversify overshadowing the underlying pathogen population. These epidemic clones have been suggested to represent highly adapted host-specific strains potentially more virulent compared to endemic or environmental strains (Siekoula-Nguedia et al., 2012). In a large-scale epidemiological study done by Nilsen et al. (2014) in countries around the Baltic Sea, the overwhelming majority of BCWD outbreaks in rainbow trout since the 1980s were attributed to a single clonal complex (CC), CC-ST10, comprising several closely related genotypes. Interestingly, CC-ST10 is also the predominant CC in several other European countries (Siekoula-Nguedia et al., 2012; Strepparava et al., 2013) as well as in Chile (Avendaño-Herrera et al., 2014) and in the United States (Van Vliet et al., 2016; Knupp et al., 2019) indicating a global distribution. Severe
BCWD outbreaks have been associated with genotypes of CC-ST10 (Nilsen et al., 2014; Van Vliet et al., 2016) whereby virulence has been considered as a major contributor to its dominance in rainbow trout farms over both time and space. It has also been proposed that a higher ability to resist antimicrobial agents and to adhere to inert and fish surfaces are factors contributing to the persistence of CC-ST10 genotypes in fish farm environments (Sundell and Wiklund, 2015; Duchaud et al., 2018). Recent evidence also shows that virulent antibiotic-resistant strains of CC-ST10 have already replaced sensitive variants in some farms in the Baltic Sea region (Söderlund et al., 2018), calling for urgent development of alternative BCWD control strategies.

The epidemiological data collected in the large scale *F. psychrophilum* MLST studies (Nilsen et al., 2014; Knupp et al., 2019) indicate that development of efficient control strategies against BCWD in rainbow trout should be directed toward genotypes of CC-ST10 or target specific virulence factors expressed by them. Targeting CC-ST10 using lytic bacteriophages isolated from the environment is theoretically an ideal approach, but so far, limited information is available regarding the phage-susceptibility of current virulent outbreak strains of *F. psychrophilum* in the Baltic Sea region. To identify current virulent strains of *F. psychrophilum*, with special focus on CC-ST10, we used a reproducible intramuscular injection-based *in vivo* challenge model for rainbow trout to obtain a LD₅₀ measure of relative virulence for current and past outbreak-related isolates from distinct geographical locations. Based on the obtained LD₅₀ values, we assessed the potential of genetic (MLST profile and serotype) and phenotypic factors (adherence, antimicrobial resistance, colony morphology, hemolytic activity, motility and proteolytic activity) as predictors of, or contributors to, *F. psychrophilum* virulence. In addition, we inquired the potential of bacteriophages as therapeutic vehicles against BCWD by examining the *in vitro* phage susceptibility of highly virulent CC-ST10 genotypes.

# MATERIALS AND METHODS

## Bacterial Isolation and Culture Preparation

Twenty-five isolates of *F. psychrophilum* included in this study were obtained from visceral organs of farmed rainbow trout in association with BCWD outbreaks in the Baltic Sea region in northern Europe between 1995 and 2017 (Table 1). For comparative analyses, the type strain NCIMB1947T isolated from coho salmon (*O. kisutsch*), previously shown to be weakly virulent in rainbow trout (Madsen and Dalsgaard, 2000), was included as a reference strain. The isolates were cultured on tryptone yeast extract salts (TYES) agar (Holt et al., 1993), consisting of 0.4% tryptone, 0.04% yeast extract, 0.05% $CaCl_2$ × 2H₂O, 0.05% MgSO₄ × 7H₂O and 1.1% agar (pH 7.2) on which rough or smooth colony morphology (Högfors-Rönnholm and Wiklund, 2010b) was determined using a stereomicroscope. For species identification, the *F. psychrophilum*-specific primers PSY1 and PSY2 (Toyama et al., 1994) were used for PCR amplification of a partial fragment (1089 bp) of the 16S rRNA gene followed by

electrophoretic confirmation of the amplification product. For test culture preparation, cryopreserved (−80°C) stock cultures of each *F. psychrophilum* isolate were cultured for 5 days on TYES agar at 15°C, followed by subcultivation in TYES broth for 48 h at 15°C under 200 rpm agitation. The bacterial cells of each broth culture were washed in fresh cooled TYES broth by centrifugation (5310 × g, 15 min, 4°C) and the optical density (OD) was spectrophotometrically adjusted to 1.0 at 520 nm, corresponding to a *F. psychrophilum* concentration of approximately 10⁹ CFU mL⁻<fnref n="1" />. The 48 h cultures were used for all challenge experiments and phenotypic characterization studies unless otherwise noted.

## Genetic Characterization of *F. psychrophilum*

### Multilocus Sequence Typing and Data Analysis

Extraction and purification of genomic DNA of each *F. psychrophilum* isolate for MLST was done using the NucleoSpin® Microbial DNA kit (Macherey-Nagel) according to the manufacturer's instructions. For MLST, partial sequences from *atpA*, *dnaK*, *fumC*, *gyrB*, *murG*, *tuf*, and *trpB* genes of the 16 previously untyped isolates (Table 1) were amplified by PCR according to the protocol described by Siekoula-Nguedia et al. (2012). Both strands of the amplified products were Sanger sequenced using the corresponding forward and reverse PCR primers. All chromatograms were manually checked for quality before assignation of the allele type (AT) and sequence type (ST) using an in-house script (P. Nicolas, INRA). Reference strains were included as internal controls for the MLST procedure and the assigned ATs and STs were submitted to the *F. psychrophilum* MLST database<fnref n="1" />. Following the assignation of ATs and STs, CCs were identified using eBURSTv3². The predicted ancestral genotype (founder) of a CC was defined as the ST with the highest number of single locus variants. STs not belonging to any CC were referred to as singletons.

### PCR Serotyping

For serotyping of all isolates used in this study, a multiplex PCR-based serotyping scheme for *F. psychrophilum* developed by Rochat et al. (2017) was used. The multiplex PCR reactions were performed using 5 µL purified genomic *F. psychrophilum* DNA (NucleoSpin® Microbial DNA, Macherey-Nagel) as the template in a 50 µL reaction mixture consisting of 2 units of DyNAzyme II DNA polymerase, 1× Optimized DyNAzyme buffer (Thermo Fisher Scientific), 10 mM dNTPs together with 10 µM of the four primer pairs described by Rochat et al. (2017). The PCR reaction mixture was heated at 94°C for 2 min, followed by 30 cycles of 94°C for 30 s, 52°C for 30 s, 72°C for 60 s, and a final extension at 72°C for 10 min. The amplified PCR products were electrophoresed (4 V cm⁻<fnref n="1" />, 60 min) on a 1% agarose-Tris-borate-EDTA gel stained with ethidium bromide and visualized under ultraviolet transillumination (Alpha Innotech Multi Image Light Cabinet). A 100 bp DNA ladder was used for fragment size estimation and subsequent classification into type 0 (188 bp), 1 (188 and 549 bp), 2 (188 and 841 bp), or 3 (188 and 361 bp).

<footnote n="1">¹http://pubmlst.org/fpsychrophilum</footnote>

<footnote n="2">²http://eburst.mlst.net</footnote>
<table><caption>TABLE 1 | The virulence of 26 *Flavobacterium psychrophilum* isolates listed by the median lethal dose (LD₅₀) value in juvenile rainbow trout.</caption><thead><tr><th>Isolate</th><th>Host species</th><th>Virulence</th><th>LD<sub>50</sub> (CFU)</th><th>ST</th><th>CC</th><th>Serotype</th><th>Country</th><th>Year</th></tr></thead><tbody><tr><td>FPS-G1</td><td>Rainbow trout</td><td>High</td><td>8.90 x 10<sup>3</sup></td><td>89*</td><td>10</td><td>1</td><td>Germany</td><td>2017</td></tr><tr><td>FPS-F15</td><td>Rainbow trout</td><td>High</td><td>2.22 x 10<sup>4</sup></td><td>2*</td><td>10</td><td>1</td><td>Finland</td><td>2017</td></tr><tr><td>FPS-P1</td><td>Rainbow trout</td><td>High</td><td>5.93 x 10<sup>4</sup></td><td>2*</td><td>10</td><td>1</td><td>Poland</td><td>2016</td></tr><tr><td>FPS-R9</td><td>Rainbow trout</td><td>High</td><td>7.33 x 10<sup>4</sup></td><td>329*</td><td>1</td><td>1</td><td>Russia</td><td>2017</td></tr><tr><td>FPS-P3</td><td>Rainbow trout</td><td>High</td><td>7.60 x 10<sup>4</sup></td><td>2*</td><td>10</td><td>1</td><td>Poland</td><td>2017</td></tr><tr><td>FPS-S6</td><td>Rainbow trout</td><td>High</td><td>8.24 x 10<sup>4</sup></td><td>92*</td><td>10</td><td>2</td><td>Sweden</td><td>2017</td></tr><tr><td>160401-1/5N</td><td>Rainbow trout</td><td>High</td><td>9.83 x 10<sup>4</sup></td><td>92*</td><td>10</td><td>2</td><td>Denmark</td><td>2016</td></tr><tr><td>P30-2B/09</td><td>Rainbow trout</td><td>Moderate</td><td>1.37 x 10<sup>5</sup></td><td>91</td><td>10</td><td>2</td><td>Finland</td><td>2009</td></tr><tr><td>F164</td><td>Rainbow trout</td><td>Moderate</td><td>1.40 x 10<sup>5</sup></td><td>240</td><td>10</td><td>2</td><td>Sweden</td><td>1996</td></tr><tr><td>FPS-F16</td><td>Rainbow trout</td><td>Moderate</td><td>1.44 x 10<sup>5</sup></td><td>2*</td><td>10</td><td>1</td><td>Finland</td><td>2017</td></tr><tr><td>V46</td><td>Rainbow trout</td><td>Moderate</td><td>1.50 x 10<sup>5</sup></td><td>138</td><td>138</td><td>0</td><td>Finland</td><td>2005</td></tr><tr><td>FPS-S11A</td><td>Rainbow trout</td><td>Moderate</td><td>1.85 x 10<sup>5</sup></td><td>92*</td><td>10</td><td>2</td><td>Sweden</td><td>2017</td></tr><tr><td>990512-1/2A</td><td>Rainbow trout</td><td>Moderate</td><td>2.06 x 10<sup>5</sup></td><td>2</td><td>10</td><td>2</td><td>Denmark</td><td>1999</td></tr><tr><td>141127-1/2N</td><td>Rainbow trout</td><td>Moderate</td><td>3.56 x 10<sup>5</sup></td><td>92*</td><td>10</td><td>2</td><td>Denmark</td><td>2014</td></tr><tr><td>P15-8B/11</td><td>Rainbow trout</td><td>Moderate</td><td>3.58 x 10<sup>5</sup></td><td>2</td><td>10</td><td>1</td><td>Finland</td><td>2011</td></tr><tr><td>160401-1/5M</td><td>Rainbow trout</td><td>Moderate</td><td>3.73 x 10<sup>5</sup></td><td>92*</td><td>10</td><td>2</td><td>Denmark</td><td>2016</td></tr><tr><td>FPS-F27</td><td>Rainbow trout</td><td>Moderate</td><td>4.33 x 10<sup>5</sup></td><td>2*</td><td>10</td><td>1</td><td>Finland</td><td>2017</td></tr><tr><td>FPS-R7</td><td>Rainbow trout</td><td>Moderate</td><td>6.16 x 10<sup>5</sup></td><td>328*</td><td>S</td><td>0</td><td>Russia</td><td>2017</td></tr><tr><td>FPS-S9</td><td>Rainbow trout</td><td>Moderate</td><td>6.83 x 10<sup>5</sup></td><td>2*</td><td>10</td><td>1</td><td>Sweden</td><td>2017</td></tr><tr><td>030522-1/1</td><td>Rainbow trout</td><td>Moderate</td><td>7.57 x 10<sup>5</sup></td><td>2</td><td>10</td><td>2</td><td>Denmark</td><td>2003</td></tr><tr><td>FPS-S10</td><td>Rainbow trout</td><td>Weak</td><td>1.36 x 10<sup>6</sup></td><td>2*</td><td>10</td><td>1</td><td>Sweden</td><td>2017</td></tr><tr><td>010418-2/1</td><td>Rainbow trout</td><td>Weak</td><td>2.13 x 10<sup>6</sup></td><td>10</td><td>10</td><td>2</td><td>Denmark</td><td>2001</td></tr><tr><td>K9/00</td><td>Rainbow trout</td><td>Weak</td><td>2.62 x 10<sup>6</sup></td><td>79</td><td>10</td><td>1</td><td>Finland</td><td>2000</td></tr><tr><td>950106-1/1</td><td>Rainbow trout</td><td>Weak</td><td>2.64 x 10<sup>6</sup></td><td>2</td><td>10</td><td>1</td><td>Denmark</td><td>1995</td></tr><tr><td>NCIMB1947<sup>T</sup></td><td>Coho salmon</td><td>Weak</td><td>3.10 x 10<sup>6</sup></td><td>13</td><td>9</td><td>0</td><td>United States</td><td>Unknown</td></tr><tr><td>FPS-S11B</td><td>Rainbow trout</td><td>Non-virulent</td><td>ND</td><td>92*</td><td>10</td><td>2</td><td>Sweden</td><td>2017</td></tr></tbody></table>

The virulence is expressed as high (LD₅₀ < 10⁵), moderate (LD₅₀ = 10⁵–10⁶), or weak (LD₅₀ > 10⁶). The non-virulent isolate FPS-S11B was excluded from further statistical analysis. The sequence type (ST) was assigned after MLST in present* and previous studies (pubmlst.org/fpsychrophilum). The clonal complex (CC) was named according to the predicted ancestral ST in an eBURSTv3 analysis (eburst.mlst.net). CC, clonal complex; CFU, colony forming units; LD₅₀, median lethal dose; ND, not defined; S, singleton not belonging to a clonal complex; ST, sequence type.

## Phenotypic Characterization of *F. psychrophilum*

### Antibiotic Susceptibility Testing

The antibiotic susceptibility of the *F. psychrophilum* isolates was tested through the Kirby-Bauer disk diffusion method on diluted (1:5) Mueller-Hinton agar medium supplemented with 5% w/v fetal calf serum. For each isolate, a 40-µL volume of bacterial TYES broth suspension containing approximately 10⁹ CFU mL⁻¹ of *F. psychrophilum* was transferred into 5 mL phosphate-buffered saline (PBS) and poured onto the agar plates. After removing excess suspension, antibiotic disks (florfenicol 30 µg, oxolinic acid 2 µg, sulfamethoxazole/trimethoprim 25 µg, tetracycline 30 µg), manufactured by Oxoid, were placed onto the agar plate. The inhibition zone diameter around each antibiotic disk was measured after 5 days of incubation at 15°C.

### Colony Spreading Motility Testing

The bacterial colony spreading motility testing was performed by spotting 5 µL triplicates of each *F. psychrophilum* isolate from a broth suspension containing 10⁹ CFU mL⁻¹ on TYES agar (0.5% agar) plates supplemented with 0.1% baker's yeast. The colony spreading motility indicated by the colony diameter (mm) of each isolate was evaluated after 6 days of incubation at 15°C.

### Testing of Proteolytic Activity

The ability of the *F. psychrophilum* isolates to hydrolyze casein, elastin and gelatin was tested on TYES agar (1.5% agar) supplemented with (w/v) skim milk (5%), elastin (0.1%), and gelatin (3%), respectively. A 1 µL volume of a 10⁹ CFU mL⁻¹ broth suspension of each isolate was spotted on the agar surface in triplicate, allowed to dry and incubated at 15°C. Following incubation, a positive reaction was indicated by a clear zone in the surrounding turbid agar medium around the inoculum. To estimate the proteolytic activity of each isolate, the clear zone ratio (clearing zone/colony diameter) was determined after 7 (gelatinase) and 10 days (caseinase, elastinase) of incubation.

### Hemolytic Activity Testing

The hemolytic activity was tested on rainbow trout erythrocytes using a microplate hemolytic assay described by Högfors-Rönnholm and Wiklund (2010a) with some modifications. All *F. psychrophilum* isolates were grown on TYES agar at 15°C for 6 days and suspended in 0.5% NaCl to an OD of 0.45 ± 0.02 at 520 nm, corresponding to a bacterial concentration of approximately 10⁸ CFU mL⁻¹. Fish blood was collected by caudal venipuncture of anesthetized rainbow trout and stored in an equal volume of Alsever's solution. Fresh erythrocytes were
harvested by centrifugation (1000 × g, 5 min, 4°C), washed with PBS (pH 7.2) and re-suspended to 5% (v/v) in PBS.

Erythrocytes (30 µL) and bacterial suspensions (30 µL) were then mixed and pipetted in triplicates into a round-bottomed microtiter plate and incubated for 24 h at 15°C with 300 rpm agitation. Subsequently, 150 µL of 0.5% NaCl was added to each well, the plate was centrifuged (1000 × g, 5 min) and the supernatants were transferred to a flat-bottomed 96-well polystyrene microtiter plate. The absorbance (A) was measured at 540 nm. The negative control (background, Abackground) contained 0.5% NaCl and erythrocytes, and the positive control (total hemolysis, A100%) contained hypotonic distilled water for osmotic burst of erythrocytes. Controls were included in triplicates on each plate. The hemolytic activity (HA) was calculated according to: $HA = (A - A_{background})/(A_{100%} - A_{background})$ and the experiment was repeated twice.

## Adhesion to Polystyrene

Adhesion to polystyrene was examined according to the procedure described by Högfors-Rönnholm et al. (2015). Briefly, the cryopreserved *F. psychrophilum* stocks were cultured for 5 days on TYES agar and recultivated another 3 days before inoculation of bacteria into sterile lake water to an OD$_{520}$ nm of 1.0, corresponding to approximately $10^9$ CFU mL$^{-1}$. Aliquots of 100 µL were added in triplicate to wells of a flat-bottomed 96-well polystyrene microtiter plate (Nunclon Δ Surface, Nunc) while sterile fresh water was used as a negative control. After static incubation at 15°C for 1 h the contents were discarded and to remove non-adherent cells the wells were washed three times with sterile 0.5% NaCl and air dried. A 125 µL volume of a 0.1% crystal violet (CV) solution was then added to each well and incubated at room temperature for 45 min. After discarding the contents, the plates were washed three times by submersion in a container of tap water and air dried. Then, 150 µL of 96% ethanol was added to each well and incubated at room temperature for 15 min. A 100 µL volume of the solubilized CV was then transferred to a flat-bottomed microtiter plate and the absorbance was quantified in a microplate reader (Victor2, Wallac) at 595 nm. Each isolate was examined in triplicate and the experiment was repeated three times.

## Estimation of Virulence in Rainbow Trout

The median lethal dose (LD$_{50}$) of each *F. psychrophilum* isolate was estimated in juvenile rainbow trout (mean weight 5 g) from the same cohort obtained from a commercial fish farm. Before start of the virulence experiments, fish were kept under laboratory conditions in tanks with a flow-through of dechlorinated tap water at approximately 13°C with continuous aeration and fed twice a day at 1% of body weight with commercial fish feed. Prior to the challenge trials, fish were fasted for 24 h and anesthetized by immersion in a 0.05 g L$^{-1}$ bath solution of benzocaine just before marking of groups by one of six fin clipping treatments (pectoral fin, pelvic fins, anal fin, adipose fin and dorsal or ventral tail fin) using scissors disinfected with 70% ethanol.

For preparation of *F. psychrophilum* cells for virulence trials, optically adjusted TYES broth suspension containing washed bacterial cells at a concentration of $10^9$ CFU mL$^{-1}$ was used

to prepare serial tenfold dilutions ($10^8-10^3$ CFU mL$^{-1}$) of each isolate in 0.5% NaCl. A 50-µL volume of five serial dilutions of *F. psychrophilum* was administered through intramuscular injection into five treatment groups consisting of seven fish each. An equal sized control group was anesthetized, marked and mock-infected with 50 µL of a sterile 0.5% NaCl solution. Disposable 1-mL syringes with a 30-gauge needle were used for the injections. Bacterial viability and concentration of each dilution was verified by colony counting on TYES agar after 7 days of incubation at 15°C. After the injection procedure, each group of fish (n = 6), i.e., 42 individuals in total per tested isolate, were transferred to 20 L test aquaria containing aerated and sponge-filtered dechlorinated tap water with a temperature of 13°C.

Mortality was recorded for 21 days during which 80% of the water was renewed daily and its pH and content of ammonia (NH$_3$), nitrate (NO$_3$), and nitrite (NO$_2$) was regularly monitored. Dead and moribund fish were removed as soon as observed. Moribund fish were euthanized by overdose of anesthesia after which the brain was mechanically destroyed. To confirm Koch's postulate, tissue samples from the kidney and spleen of dead fish were streaked onto TYES agar plates and incubated at 15°C for 7 days. Bacteria forming yellow colonies on the agar plates were identified as *F. psychrophilum* by species-specific PCR as described previously in the text. The LD$_{50}$ for each tested isolate was then estimated by the Reed-Muench method (Reed and Muench, 1938), which requires cumulative mortality to be above 50% in at least one of the treatment groups and below 50% in another.

## Phage Susceptibility Testing

As part of the characterization of phenotypic differences between the strains, the phage susceptibility of the bacterial isolates was tested *in vitro* against a selection of 5 phages, representing different host range patterns and different sources, locations and times of isolation (Table 2). These phages included the previously isolated bacteriophage FpV4 (Stenholm et al., 2008) and 4 newly isolated (FPSV-D19, FPSV-D22, FPSV-S8, and FPSV-S20) phages (Table 2 and Figure 1). The new phages were isolated from fish farms during BCWD outbreaks in 2016–2017.

Briefly, bacteriophages were isolated from water samples or fish tissue from Danish and Swedish freshwater rainbow trout farms by using single enrichment cultures. The five distinct enrichment hosts for the phages were selected based on shared geographic origin (country) of isolation (Table 2). Approximately 20 mL of a 0.2-µm-filtered water sample was mixed with 2 L 10× TYES broth and with 2 mL of the isolates to be used in the enrichment. Fish were homogenized in a blender and mixed with an equal volume of saline magnesium (SM) buffer (50 mM Tris-HCl, pH 7.5, 99 mM NaCl, 8 mM MgSO$_4$) to eluate phages. The mixture was centrifuged (5000 × g, 15 min, 4°C) and the supernatant 0.2-µm-filtered. One mL of the eluate was mixed with 20 ml of TYES broth and 2 ml of host isolate. The enrichment cultures were incubated for 5–7 days at 15°C with agitation to allow proliferation of lytic *F. psychrophilum* bacteriophages. Following incubation, the culture was 0.2-µm-filtered and the
<table><caption>**TABLE 2** | The source of isolation and *F. psychrophilum* hosts used for enrichment and proliferation of the five bacteriophages included in this study.</caption><thead><tr><th>Phage</th><th>Enrichment host</th><th>Proliferation host</th><th>Source</th><th>Country</th><th>Year of isolation</th></tr></thead><tbody><tr><td>FpV4</td><td>950106-1/1</td><td>950106-1/1</td><td>Water sample</td><td>Denmark</td><td>2005</td></tr><tr><td>FPSV-D19</td><td>FPS-D10</td><td>FPS-S6</td><td>Water sample</td><td>Denmark</td><td>2017</td></tr><tr><td>FPSV-D22</td><td>FPS-D5</td><td>FPS-S6</td><td>Fish tissue</td><td>Denmark</td><td>2017</td></tr><tr><td>FPSV-S8</td><td>FPS-S2</td><td>FPS-S6</td><td>Water sample</td><td>Sweden</td><td>2017</td></tr><tr><td>FPSV-S20</td><td>FPS-S30</td><td>FPS-S6</td><td>Water sample</td><td>Sweden</td><td>2017</td></tr></tbody></table>

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

presence of bacteriophages in the filtrate was detected by the double-layer method (Stenholm et al., 2008). For purification of bacteriophages, a single plaque was picked with a plastic Pasteur pipette and the phages were eluted in SM buffer with shaking for 2 h at 15°C or overnight at 4°C. After centrifugation (5000 × g, 15 min, 4°C), the supernatant was transferred to a new clean tube and used for re-infection. Bacteriophages were isolated as single plaques by three repeated rounds of plaque purification.

For preparation of high-density bacteriophage stocks, 1 mL bacteriophage stock was added to 20 mL exponentially growing cultures ($OD_{600 nm} = 0.2$) of a suitable proliferation host (Table 2) and incubated at 15°C for 2–3 days. Lysed

bacterial cultures were treated with chloroform and subsequent centrifugation (5000 × g, 15 min, 4°C). Phage quantification was performed by small spot assay (Mazzocco et al., 2009). Finally, the efficiency of plating (EOP) of the isolated bacteriophages were determined exposing the *F. psychrophilum* isolates to the same phage titer [~$10^8$ plaque forming units (PFU) mL$^{-1}$] and infectivity was quantified by the small drop plaque assay. PFU were examined after 5 days of incubation at 15°C. Each experiment was performed three independent times and the EOP was presented as the infection efficiency for individual phages, relative to the highest titer obtained for each phage (EOP = 1) on an isolate in the collection.
## Statistical Analyses

As higher LD₅₀ values indicate lower virulence, phenotypic characteristics correlating negatively with LD₅₀ indicate that when these descriptors get bigger, the LD₅₀ value will decrease and the virulence thus increase. Hence, a correlation analysis was performed to identify potential predictors of virulence through identification of statistically significant negative correlations (α < 0.01) between LD₅₀ and phenotypic characteristics measured in the study. Statistical significance of the studied phenotypic characteristics on *F. psychrophilum* virulence was tested with a Pearson's correlations test and a multiple linear regression analysis in GraphPad Prism 8.0.1 with LD₅₀ (sample size = 25) as the dependent variable and the measured numerical values for adhesion, colony spreading motility, protease activity (caseinase, gelatinase, and elastinase), hemolytic activity, and antibiotic susceptibility (florfenicol, sulfamethoxazole/trimethoprim, oxolinic acid, and tetracycline) as predictors (n = 10). To fulfill the criterion for multiple regression, normality and linearity was verified with a non-significant Shapiro-Wilk test and normal probability plots on residual histogram and P-P plots respectively. The absence of collinearity was also confirmed by the absence of variance inflation factors (VIF) above 3.427 (Table 3). Autocorrelation was checked in IBM SPSS Statistics version 25 with the Durbin Watson test, which was below 2. Homoscedasticity of the data was tested with a Koenker test, and a Breusch-Pagan test using a SPSS macro written by *Pryce and Garcia-Granero* (2002).

## RESULTS

### Genetic Characteristics and Virulence

Analysis of MLST data of the 26 isolates included in this study together with the 1097 isolates in the global *F. psychrophilum* MLST database<fnref n="3" /> (last access January 17th 2019) showed that the majority of isolates (13/15) obtained from recent BCWD outbreaks (2016–2017) belonged to CC-ST10 and consisted

of the genotypes ST2, ST10, ST79, ST89, ST91, and ST92 (Table 1). Two novel STs, ST328 and ST329, not belonging to CC-ST10 were isolated from Russia in 2017. The singleton ST328 did not group into any CCs while ST329 is part of CC-ST1. Of all isolates included in this study, 22 of 26 were of CC-ST10 related STs with ST2 being the most prevalent (n = 11) genotype. The genetic serotyping scheme (Rochat et al., 2017) classified 12 *F. psychrophilum* isolates into type 1, 11 isolates into type 2 and 3 isolates into type 0 (Table 1). None of the serotype 0 isolates belonged to CC-ST10.

In the challenge experiments, the examined diseased fish exhibited gross signs commonly associated with BCWD with muscle ulceration and tissue necrosis at the site of injection, and swollen internal organs. Koch's postulate was fulfilled as *F. psychrophilum* was recovered and identified from the visceral organs (spleen and kidney) of experimentally infected fish. Mortality in fish was generally dose dependent, i.e., with increasing bacterial dose, more mortality was achieved with each *F. psychrophilum* isolate as seen in the survival plots (Supplementary Figure S1). The disease signs were similar for all isolates and the onset of mortality occurred 4–7 days post challenge. None of the mock-infected fish individuals died or showed any symptoms of disease during the challenge trials.

The results from the virulence tests showed a wide range in the degree of pathogenicity in *F. psychrophilum*, with up to 3-log differences in quantifiable LD₅₀ between isolates. The LD₅₀ of the 26 tested *F. psychrophilum* isolates ranged from 8.9 × 10³ (FPS-G1) to 3.1 × 10⁶ CFU (NCIMB1947T). The *F. psychrophilum* isolates high in virulence were obtained from 6 different countries in the Baltic Sea region and from recent (2016–2017) disease outbreaks. The degree of pathogenicity was determined to correspond to log-differences in LD₅₀ and virulence was considered high (n = 7) when LD₅₀ < 10⁵ CFU, moderate (n = 13) when LD₅₀ = 10⁵–10⁶ CFU or weak (n = 5) when LD₅₀ > 10⁶ CFU (Table 1). One *F. psychrophilum* isolate, FPS-S11B, did not induce mortality or necrosis in fish tissue after intramuscular injection of a dose as high as 10⁷ CFU and was thus

<table><caption>TABLE 3 | Results of multiple linear regression analysis (R² = 0.55) for predictors of *F. psychrophilum* virulence (LD₅₀).</caption><thead><tr><th>Predictors of LD<sub>50</sub></th><th>Estimate</th><th>SE</th><th>Lower – upper 95% CI</th><th>VIF</th><th>|t|</th><th>P-value</th></tr></thead><tbody><tr><td>Adhesion</td><td>659596</td><td>1539817</td><td>-2642983 – 3962174</td><td>2.618</td><td>0.501</td><td>0.624</td></tr><tr><td>Motility</td><td>-18178</td><td>22040</td><td>-65449 – 29093</td><td>2.140</td><td>0.428</td><td>0.675</td></tr><tr><td>Caseinase activity</td><td>-912870</td><td>594550</td><td>-2188054 – 362314</td><td>1.778</td><td>0.825</td><td>0.423</td></tr><tr><td>Gelatinase activity</td><td>-287295</td><td>503758</td><td>-1367748 – 793158</td><td>1.869</td><td>1.535</td><td>0.147</td></tr><tr><td>Elastinase activity</td><td>308817</td><td>291965</td><td>-317386 – 935021</td><td>2.917</td><td>0.570</td><td>0.578</td></tr><tr><td>Hemolytic activity</td><td>2568779</td><td>1638753</td><td>-945997 – 6083555</td><td>3.427</td><td>1.058</td><td>0.308</td></tr><tr><td>Florfenicol susceptibility</td><td>6303</td><td>39398</td><td>-78196 – 90803</td><td>2.728</td><td>1.568</td><td>0.139</td></tr><tr><td>Oxolinic acid susceptibility</td><td>13762</td><td>21319</td><td>-31963 – 59488</td><td>3.042</td><td>0.160</td><td>0.875</td></tr><tr><td>Sulfamethoxazole/trimethoprim susceptibility</td><td>35153</td><td>51398</td><td>-75084 – 145391</td><td>1.906</td><td>0.646</td><td>0.529</td></tr><tr><td>Tetracycline susceptibility</td><td>18851</td><td>17643</td><td>-18989 – 56692</td><td>1.314</td><td>0.684</td><td>0.505</td></tr></tbody></table>

The P-value reported in the last column corresponds to the P-value of the t-test used to check the significance of individual regression coefficients at the 5% level in the multiple linear regression model. LD₅₀, median lethal dose; SE, standard error; CI, confidence interval; VIF, variance inflation factors.

<footnote n="3">³http://pubmlst.org/fpsychrophilum</footnote>
considered non-virulent and excluded from further statistical analyses (Table 1).

The seven highly virulent *F. psychrophilum* isolates were isolated from recent (2016–2017) disease outbreaks in Denmark, Finland, Germany, Poland, Russia and Sweden. Genetic analyses showed that *F. psychrophilum* isolates high in virulence belonged to the geographically widely distributed clonal complexes CC-ST1 and CC-ST10 (Table 1). Interestingly, the most prevalent *F. psychrophilum* genotype in the global MLST database<fnref n="4" /> ST2 was represented in our study by isolates with high, moderate and weak virulence (Table 1). Moreover, the non-virulent *F. psychrophilum* isolate FPS-S11B and the highly virulent isolates FPS-S6 and 160401-1/5N belonged to the same MLST genotype (ST92), highlighting that pathogenic and non-pathogenic *F. psychrophilum* isolates can share the same ST. All three serotypes caused mortality in rainbow trout, however, all highly virulent *F. psychrophilum* isolates belonged to serotype 1 or 2,

while two moderately virulent isolates and the weakly virulent *F. psychrophilum* typestrain NCIMB1947T belonged to serotype 0 (Table 1).

## Phenotypic Characteristics and Prediction of Virulence-Associated Traits

The univariate correlation studies (Figure 2) only showed non-significant associations between the measured phenotypic characteristics and LD₅₀. The results from the multivariable analysis (R² = 0.55) showed that 55% of the variation in the determined LD₅₀ values could be explained by the predictor variables used. The multivariable analysis, however, neither revealed any significant predictors of *F. psychrophilum* virulence (Table 3).

The results from the antibiotic disk diffusion tests showed that all isolates were inhibited by florfenicol and tetracycline, while several of the highly, moderately and weakly virulent isolates were uninhibited by oxolinic acid and sulfamethoxazole/trimethoprim (Table 4). The non-virulent isolate FPS-S11B was deficient in

<table><caption>TABLE 4 | Phenotypic characteristics of 25 virulent *Flavobacterium psychrophilum* isolates and the non-virulent FPS-S11B listed according to decreasing virulence in juvenile rainbow trout (*Oncorhynchus mykiss*).</caption><thead><tr><th rowspan="2">F. psychrophilum Isolate code</th><th rowspan="2">S/R</th><th rowspan="2">Abs@595 Adhesion</th><th rowspan="2">Colony spreading (mm)</th><th rowspan="2">Hemolytic activity</th><th colspan="3">Protease activity (clear zone ratio)</th><th colspan="4">Inhibition zone diameter (mm)</th></tr><tr><th>Caseinase</th><th>Gelatinase</th><th>Elastinase</th><th>FF</th><th>OA</th><th>SXT</th><th>TC</th></tr></thead><tbody><tr><td>FPS-G1</td><td>S</td><td>0.355 ± 0.014</td><td>34.3 ± 1.3</td><td>0.345 ± 0.025</td><td>3.6 ± 0.1</td><td>5.2 ± 0.4</td><td>3.3 ± 0.2</td><td>43</td><td>8</td><td>0</td><td>41</td></tr><tr><td>FPS-F15</td><td>R</td><td>0.158 ± 0.004</td><td>43.0 ± 1.2</td><td>0.202 ± 0.055</td><td>3.3 ± 0.1</td><td>4.5 ± 0.1</td><td>2.6 ± 0.2</td><td>47</td><td>0</td><td>0</td><td>58</td></tr><tr><td>FPS-P1</td><td>R</td><td>0.153 ± 0.005</td><td>27.0 ± 1.2</td><td>0.093 ± 0.022</td><td>2.5 ± 0.1</td><td>5.1 ± 0.3</td><td>2.8 ± 0.1</td><td>66</td><td>0</td><td>9</td><td>41</td></tr><tr><td>FPS-R9</td><td>S</td><td>0.552 ± 0.012</td><td>32.0 ± 6.5</td><td>0.163 ± 0.029</td><td>2.9 ± 0.1</td><td>6.1 ± 0.1</td><td>1.0 ± 0.0</td><td>54</td><td>43</td><td>0</td><td>70</td></tr><tr><td>FPS-P3</td><td>R</td><td>0.155 ± 0.006</td><td>44.0 ± 5.0</td><td>0.096 ± 0.007</td><td>2.6 ± 0.1</td><td>5.2 ± 0.4</td><td>2.7 ± 0.3</td><td>36</td><td>0</td><td>0</td><td>30</td></tr><tr><td>FPS-S6</td><td>S</td><td>0.598 ± 0.014</td><td>28.7 ± 1.3</td><td>0.250 ± 0.047</td><td>3.2 ± 0.1</td><td>4.4 ± 0.3</td><td>1.0 ± 0.0</td><td>56</td><td>0</td><td>0</td><td>49</td></tr><tr><td>160401-1/5N</td><td>S</td><td>0.415 ± 0.008</td><td>28.3 ± 2.0</td><td>0.358 ± 0.012</td><td>3.7 ± 0.1</td><td>4.7 ± 0.4</td><td>1.0 ± 0.0</td><td>50</td><td>0</td><td>7</td><td>44</td></tr><tr><td>P30-2B/09</td><td>S</td><td>0.568 ± 0.017</td><td>22.7 ± 1.9</td><td>0.144 ± 0.004</td><td>3.1 ± 0.1</td><td>4.7 ± 0.3</td><td>1.0 ± 0.0</td><td>44</td><td>0</td><td>0</td><td>37</td></tr><tr><td>F164</td><td>R</td><td>0.134 ± 0.003</td><td>10.3 ± 0.7</td><td>0.087 ± 0.029</td><td>3.3 ± 0.1</td><td>4.7 ± 0.2</td><td>2.4 ± 0.1</td><td>56</td><td>9</td><td>0</td><td>58</td></tr><tr><td>FPS-F16</td><td>R</td><td>0.142 ± 0.004</td><td>37.3 ± 0.9</td><td>0.194 ± 0.066</td><td>3.2 ± 0.1</td><td>4.7 ± 0.1</td><td>2.8 ± 0.1</td><td>54</td><td>7</td><td>0</td><td>60</td></tr><tr><td>V46</td><td>R</td><td>0.138 ± 0.003</td><td>38.7 ± 3.0</td><td>0.829 ± 0.077</td><td>4.2 ± 0.3</td><td>5.0 ± 0.3</td><td>1.0 ± 0.0</td><td>40</td><td>31</td><td>12</td><td>44</td></tr><tr><td>FPS-S11A</td><td>S</td><td>0.515 ± 0.011</td><td>45.0 ± 2.1</td><td>0.061 ± 0.012</td><td>3.3 ± 0.2</td><td>4.7 ± 0.3</td><td>1.0 ± 0.0</td><td>52</td><td>6</td><td>9</td><td>64</td></tr><tr><td>990512-1/2A</td><td>R</td><td>0.141 ± 0.004</td><td>19.3 ± 3.0</td><td>0.318 ± 0.036</td><td>3.4 ± 0.1</td><td>5.0 ± 0.3</td><td>3.7 ± 0.1</td><td>49</td><td>0</td><td>8</td><td>55</td></tr><tr><td>141127-1/2N</td><td>S</td><td>0.452 ± 0.008</td><td>26.3 ± 2.4</td><td>0.143 ± 0.019</td><td>3.3 ± 0.1</td><td>5.2 ± 0.1</td><td>1.0 ± 0.0</td><td>50</td><td>8</td><td>7</td><td>45</td></tr><tr><td>P15-8B/11</td><td>R</td><td>0.162 ± 0.005</td><td>24.0 ± 1.0</td><td>0.188 ± 0.029</td><td>3.1 ± 0.1</td><td>4.9 ± 0.1</td><td>3.1 ± 0.1</td><td>57</td><td>6</td><td>0</td><td>58</td></tr><tr><td>160401-1/5M</td><td>S</td><td>0.365 ± 0.015</td><td>20.0 ± 2.6</td><td>0.321 ± 0.085</td><td>3.6 ± 0.1</td><td>5.8 ± 0.4</td><td>1.0 ± 0.0</td><td>49</td><td>7</td><td>0</td><td>44</td></tr><tr><td>FPS-F27</td><td>R</td><td>0.132 ± 0.004</td><td>38.7 ± 1.9</td><td>0.183 ± 0.016</td><td>3.1 ± 0.1</td><td>5.7 ± 0.3</td><td>2.7 ± 0.2</td><td>44</td><td>0</td><td>8</td><td>48</td></tr><tr><td>FPS-R7</td><td>S</td><td>0.526 ± 0.027</td><td>24.3 ± 1.5</td><td>0.257 ± 0.048</td><td>2.8 ± 0.2</td><td>4.9 ± 0.1</td><td>1.0 ± 0.0</td><td>57</td><td>16</td><td>0</td><td>67</td></tr><tr><td>FPS-S9</td><td>R</td><td>0.138 ± 0.004</td><td>28.0 ± 1.0</td><td>0.100 ± 0.001</td><td>3.6 ± 0.2</td><td>5.6 ± 0.6</td><td>2.9 ± 0.2</td><td>52</td><td>10</td><td>0</td><td>64</td></tr><tr><td>030522-1/1</td><td>R</td><td>0.127 ± 0.004</td><td>15.0 ± 1.5</td><td>0.344 ± 0.052</td><td>3.4 ± 0.1</td><td>4.7 ± 0.2</td><td>3.0 ± 0.1</td><td>59</td><td>0</td><td>8</td><td>53</td></tr><tr><td>FPS-S10</td><td>R</td><td>0.137 ± 0.005</td><td>15.3 ± 1.8</td><td>0.127 ± 0.009</td><td>3.3 ± 0.1</td><td>5.3 ± 0.3</td><td>2.8 ± 0.1</td><td>54</td><td>0</td><td>7</td><td>65</td></tr><tr><td>010418-2/1</td><td>R</td><td>0.135 ± 0.005</td><td>13.0 ± 1.1</td><td>0.357 ± 0.007</td><td>3.5 ± 0.1</td><td>4.8 ± 0.3</td><td>3.0 ± 0.1</td><td>54</td><td>0</td><td>7</td><td>64</td></tr><tr><td>K9/00</td><td>R</td><td>0.177 ± 0.004</td><td>25.0 ± 2.1</td><td>0.381 ± 0.001</td><td>3.2 ± 0.1</td><td>4.8 ± 0.4</td><td>2.7 ± 0.1</td><td>64</td><td>44</td><td>0</td><td>44</td></tr><tr><td>950106-1/1</td><td>S</td><td>0.397 ± 0.011</td><td>24.7 ± 0.3</td><td>0.359 ± 0.067</td><td>3.2 ± 0.1</td><td>4.8 ± 0.3</td><td>2.9 ± 0.1</td><td>54</td><td>0</td><td>7</td><td>59</td></tr><tr><td>NCIMB1947T</td><td>R</td><td>0.117 ± 0.002</td><td>54.7 ± 0.3</td><td>0.746 ± 0.216</td><td>2.7 ± 0.1</td><td>4.1 ± 0.1</td><td>1.0 ± 0.0</td><td>56</td><td>32</td><td>12</td><td>65</td></tr><tr><td>FPS-S11B</td><td>R</td><td>0.110 ± 0.003</td><td>7.7 ± 0.3</td><td>0.222 ± 0.022</td><td>1.0 ± 0.0</td><td>1.0 ± 0.0</td><td>1.0 ± 0.0</td><td>57</td><td>8</td><td>0</td><td>62</td></tr></tbody></table>

R indicates rough colony morphology and S indicates smooth. A clear zone ratio of 1.0 ± 0.0 indicates no activity. FF, florfenicol; OA, oxolinic acid; SXT, sulfamethoxazole/trimethoprim; TC, tetracycline.

<footnote n="4">⁴http://pubmlst.org/fpsychrophilum</footnote>
colony spreading motility and proteolytic activity, and also showed the weakest adherence capacity (Table 4). All other isolates showed colony spreading motility and caseinolytic and gelatinolytic activity, while nearly half of the isolates were unable to degrade elastin (Table 4). Interestingly, none of the serotype 0 isolates or highly virulent ST92 genotypes showed any elastinolytic activity, which was present in the three most virulent outbreak isolates. Because elastinase production was observed only in 4 out of 7 highly virulent isolates, it cannot be used as a sole predictor of *F. psychrophilum* virulence. The causation behind the high negative correlation between adhesion and elastinolytic activity (Figure 2) was the inability to degrade elastin observed in the majority (8 out of 10) of the adhesive smooth morphotype isolates. Adhesion was clearly associated with the colony morphology of the isolate, smooth *F. psychrophilum* morphotypes being more adherent than rough ($R^2 = 0.90, P < 0.0001$).

## Phage Susceptibility

In general, 16 out of 26 *F. psychrophilum* strains were sensitive to at least three bacteriophages, with a 10⁶- fold variation in EOP for each phage depending on the host (Figure 1). Only

MLST genotypes belonging to the CC-ST10 cluster were infected, whereas more distantly related strains were resistant to infection by the five phages. However, within CC-ST10, there was no direct correlation between ST and phage susceptibility, and each phage was able to infect strains from 4 to 5 different STs. The group of phage-sensitive *F. psychrophilum* isolates included 5 (FPS-G1, FPS-P1, FPS-P3, FPS-S6, and 160401-1/5N) of the seven most virulent strains which represented three different STs (ST2, ST89, and ST92) within CC-ST10, whereas the other highly virulent strains (FPS-F15 and FPS-R9) were resistant to phage infection (Figure 1).

## DISCUSSION

Pathogenic diversity within bacterial species is a known phenomenon and the study of virulence in bacteria causing infectious diseases in fish has primarily been focused on studying one or a few characteristics at a time (Sudheesh et al., 2012). In many bacterial species, the determinants regulating virulence are however difficult to differentiate because pathogenic strains can express several virulence factors simultaneously and a variety

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

**FIGURE 2** | Pearson's correlation matrix visualized as a heat map plot with LD₅₀ and the phenotypic characteristics tested in the study. Asterisk (*) indicates significant difference from 0 at the level of α = 0.01. LD₅₀, median lethal dose; FF, florfenicol susceptibility; OA, oxolinic acid susceptibility; SXT, sulfamethoxazole/trimethoprim susceptibility; TC, tetracycline susceptibility.
of factors can contribute to the overall pathogenicity (Cutler, 1991). Although virulence is considered a universal pathogen characteristic, it is expressed only in a susceptible host. Therefore, a measure of virulence involves an assessment of a host-pathogen interaction, most commonly morbidity or mortality. Identification of factors required for virulence is not only important for understanding the development of disease, but also for discovery of targets for pathogen control. Alternatives to antibiotic administration that kill or inhibit growth of pathogens could be achieved by development of anti-virulence drugs that disrupt or inhibit key pathogenicity or virulence factors that are general across the species including different strains.

Here we used a large collection of genetically and phenotypically characterized *F. psychrophilum* isolates collected from BCWD outbreaks in the Baltic Sea region for exploration of predictors for pathogenicity and virulence as measured by mortality (LD₅₀) in juvenile rainbow trout through a reproducible intramuscular injection-based infection model. Inclusion of the weakly virulent reference strain NCIMB1947T and a large number of isolates demonstrating a wide 3-log range in LD₅₀ enabled us to evaluate the relationships between specific genetic and phenotypic factors and virulence.

In a previous virulence study, Garcia et al. (2000) showed a 2-log variation in LD₅₀ among six outbreak isolates of *F. psychrophilum* after intramuscular injection experiments in juvenile rainbow trout. The observed wide range in virulence among *F. psychrophilum* strains in previous and this study sets limits when selecting for the range of test doses (CFU) for LD₅₀ determination, which in turn affects the accuracy of LD₅₀ predictions. Some variation in the predicted LD₅₀, particularly among isolates of the same genotype, might be explained by experimental variability that may have been further revealed in experimental replicates. However, for LD₅₀ determination of a total of 26 isolates in this study (Table 1), we used five test doses with a separation factor of 10 (Supplementary Figure S1), and 7 fish individuals for each dose without replication to minimize the number of experimental animals to be sacrificed, while still meeting the minimum standards described in the Organisation for Economic Cooperation and Development [OECD] (2019) Guidelines for Testing of Chemicals (Test No. 203: Fish, Acute Toxicity Test). For a more definitive test of LD₅₀, smaller separation factors and replication are likely necessary.

Due to the lack of predictors of pathogenicity in *F. psychrophilum*, there was no rationale for a *priori* exclusion or prioritization of potential virulence factors in our multivariable analysis (McClelland et al., 2006). Instead of estimating the relative contributions of individual virulence factors, our aim was to identify large effects of selected phenotypic characteristics using the largest number of *F. psychrophilum* isolates tested for *in vivo* virulence to date.

## Genetic Predictors of Virulence

The resolution of the current MLST scheme was not sufficient to differentiate highly virulent from weakly or non-virulent *F. psychrophilum* isolates. Interestingly, pathogenic CC-ST10 strains showed up to a 100-fold variation in virulence not only within the complex but also within the most prevalent

isolated ST2. The observed variation in virulence within isolates of identical STs could be explained by diversity in virulence gene content or by differences in the expression of genes. The identification of weakly virulent outbreak isolates of CC-ST10 indicates that the severity of BCWD outbreaks is attributed not only to the virulence of *F. psychrophilum*, but also to the state of host fish as well as environmental factors.

With the recently developed PCR-based serotyping scheme, Rochat et al. (2017) showed that the genetically determined serotypes 0, 1 and 2 correspond to the previously proposed *F. psychrophilum* serogroups FpT, Th and Fd, respectively, determined by slide-agglutination tests (Lorenzen and Olesen, 1997). Serotype 1 (serogroup Th) and 2 (serogroup Fd) have previously been shown to be predominant in BCWD outbreaks (Lorenzen and Olesen, 1997; Dalsgaard and Madsen, 2000) and more virulent than serotype 0 (serogroup FpT) in experimentally infected rainbow trout (Madsen and Dalsgaard, 2000; Madetoja et al., 2002). A significant association between the serotype and the host fish species has also been found, with serotype 1 and 2 being more prevalent in rainbow trout and strains infecting coho salmon belonging more often to serotype 0 (Rochat et al., 2017). These findings are in keeping with the results from the present study where the ten most virulent outbreak isolates belonged to serotype 1 and 2 and none of the serotype 0 isolates belonged to the dominant rainbow trout-associated (Nilsen et al., 2014) clonal complex CC-ST10 (Table 1). The serotype diversity (serotype 1 and 2) found among moderately virulent ST2 isolates (Table 1) might also explain some of the observed variation in ST2 virulence and further complicate control of BCWD by vaccination.

## Phenotypic Predictors of Virulence

CC-ST10 isolates from BCWD outbreaks in rainbow trout have previously been associated with a higher ability to adhere to inert and mucus coated surfaces compared to genetically distinct *F. psychrophilum* strains isolated from other sources (Sundell and Wiklund, 2015). Our results confirmed previous observations that smooth colony phenotypes are superior to rough in adhering to inert surfaces (Högfors-Rönnholm et al., 2015) but equally virulent (Figure 2; $R^2$ = 0.04, *P* = 0.13) when injected into rainbow trout (Högfors-Rönnholm and Wiklund, 2010b). Similarly, virulence in *F. psychrophilum* could not be predicted by the hemolytic activity in the microplate assay, but here, contrary to previous findings by Högfors-Rönnholm and Wiklund (2010a), we did not detect any significant differences ($R^2$ = 0.01, *P* = 0.59) in the activity between rough and smooth isolates.

Gliding motility facilitates the movement of bacteria along surfaces in *Flavobacterium* species and results in the spreading of colonies. Gliding motility in *F. psychrophilum* is not only important for colonization, but also for expression of virulence factors through cooperative function with secretion systems (Pérez-Pascual et al., 2017). Gliding motility in *F. psychrophilum* is associated with the protein secretion system PorSS which is also believed to influence secretion of extracellular enzymes (Castillo et al., 2015), suggesting that these phenotypic properties are involved in virulence. Studies with *F. psychrophilum* mutants
deficient in colony spreading showed a significant reduction in gelatinolytic and caseinolytic activity, and attenuated virulence compared to the wild-type isolate (Pérez-Pascual et al., 2017). Excretion of proteases degrading casein and gelatin is a common characteristic of *F. psychrophilum* whilst elastinase production varies between strains (Bertolini et al., 1994; Sundell and Wiklund, 2015). The unique genetic determinant allowing for elastin degradation was just recently discovered (Rochat et al., 2019) and a strong elastinolytic activity was furthermore associated with CC-ST10 genotypes. Although elastinase production was required for maximal virulence in our infection model, elastinolytic activity was not a significant predictor of virulence. The absence of elastinolytic activity in two highly virulent outbreak isolates of genotype ST92 belonging to CC-ST10 indicates that elastinase production is neither obligatory for high virulence nor pathogenicity in rainbow trout. Proteases secreted by fish pathogens have the potential to destroy proteins, which constitute the structure and function of the host fish tissues and proteins involved in host immune defense. However, demonstrating that such interactions occur during BCWD outbreaks in fish farms is difficult as expression of proteases in *F. psychrophilum* may be under complex regulation.

Our results demonstrated a high variability in caseinolytic and gelatinolytic activity, and colony spreading within the groups of high, moderate and weak virulence, and that the inability of one isolate (FPS-S11B) to cause mortality in fish was associated with its unique lack of protease production and deficiency in colony spreading (Table 4). Our findings are supported by results obtained through transposon-based mutagenesis of *F. psychrophilum*, where inhibition of colony spreading and extracellular proteolytic activity resulted in a complete loss of pathogenicity in rainbow trout (Pérez-Pascual et al., 2015). We therefore conclude that some levels of colony spreading motility and protease production are likely to be required for pathogenicity in *F. psychrophilum*, but that no further gains in virulence will necessarily occur if colony spreading and/or proteolytic activity increase.

## Antibiotic Resistance and Virulence

Recent in depth analysis of genome data revealed that recombination is the key driver of *F. psychrophilum* evolution and genes specific to CC-ST10 suggests a possible link between antibiotic resistance and virulence (Duchaud et al., 2018). Virulence and antibiotic resistance are similar in that some of the determinants are transmitted between species by horizontal gene transfer. Additional factors that are common to virulence and resistance include important adaptation-related processes that activate or repress the expression of various genes (Yeung et al., 2011). Although antibiotic resistance is not in itself a virulence factor, the relationship between resistance and virulence can be of benefit to a pathogen in environments such as fish farms where antibiotics are frequently used. Therefore, antibiotic resistance may be considered a virulence-like factor, which antibiotic-resistant pathogens or strains can utilize in environments where selective antibiotic pressure prevails (Beceiro et al., 2013).

Perhaps alarmingly, a recent temporal genetic analysis of *F. psychrophilum* isolates from BCWD outbreaks in Swedish fish

farms between 1988 and 2016 revealed a shift in domination from quinolone-sensitive ST79 clones in 1988–1997 to resistant ST92 clones in 2011–2016 (Söderlund et al., 2018). Quinolone-resistance in *F. psychrophilum* has also recently been associated with severe BCWD outbreaks in Finland (Sundell et al., 2013) and Norway (Shah et al., 2012). Although we did not find any associations between decreasing antibiotic inhibition zone diameters and increasing virulence, our study showed that the CC-ST10 isolates FPS-S6, 160401-1/5N, FPS-S11A, 141127-1/2N, and 160401-1/5M of ST92 isolated between 2014 and 2017 were less inhibited by oxolinic acid in the disk diffusion tests and more virulent compared to the CC-ST10 isolate K9/00 of ST79 isolated in year 2000 (Table 4). These findings provide additional support to the hypothesis of a more recent dissemination of the virulent quinolone-resistant *F. psychrophilum* clone ST92 to the Baltic Sea region (Söderlund et al., 2018).

According to the global MLST database<fnref n="5" />, ST92 was first reported in Denmark in 2006 and thereafter in several countries including Finland, France, Sweden and Switzerland indicating a rapid and wide geographical dissemination of a quinolone-resistant pathogenic clone. Perhaps even more alarmingly, one highly virulent outbreak strain (FPS-P3) reportedly responsible for massive fish farm mortality in 2017 belonged to the most prevalent *F. psychrophilum* genotype ST2 of CC-ST10 and showed the least inhibition by florfenicol and tetracycline (Table 4), both commonly used antibiotics for treatment of BCWD. Reduced antibiotic susceptibility in *F. psychrophilum* isolated from disease outbreaks is however not restricted to the Baltic Sea region as it has also frequently been observed among North and South American isolates (Hesami et al., 2010; Miranda et al., 2016). The development of resistance in highly pathogenic and prevalent *F. psychrophilum* CC-ST10 variants and their potential of rapid international spread emphasizes that there is an urgent need for an alternative therapeutic agent to antibiotics for control of BCWD.

## Phage-Susceptibility of CC-ST10

The isolation and identification of lytic bacteriophages able to infect prevalent and particularly antibiotic-resistant *F. psychrophilum* variants could be of crucial importance for the development of a novel control intervention for BCWD. Interestingly, the isolated bacteriophages showed lytic activity against most highly virulent outbreak isolates and, particularly importantly, against several important CC-ST10 genotypes including the emerging quinolone-resistant clone ST92 (Söderlund et al., 2018). The phage-susceptible strains belonged to distinct STs and represented isolates from 5 different countries, isolated over a period of two decades. This emphasizes the broad host-range and long-term stability of the infectious properties of individual *F. psychrophilum* phages, providing further reason for investigating the potential of phage-based control of BCWD. However, biosafety issues and the emergence of phage-resistance in the bacterial population with still unpredictable functional implications in *F. psychrophilum* need

<footnote n="5">⁵http://pubmlst.org/fpsychrophilum</footnote>
thorough investigation for bacteriophages to be considered in applications to control BCWD in fish farm environments.

## CONCLUSION

We have characterized the *in vivo* virulence of a large collection of *F. psychrophilum* isolates obtained from BCWD outbreaks in the Baltic Sea region, and correlated this with the expression of several phenotypic virulence-related characteristics *in vitro*. Our results show that isolates of *F. psychrophilum* retrieved from BCWD outbreaks can show up to a 3-log variation in virulence, but exhibit similar genotypic and phenotypic characteristics. However, colony spreading motility and production of caseinolytic and gelatinolytic enzymes were universal observable characteristics amongst isolates pathogenic to rainbow trout, and may function as a priori plausible predictors of *F. psychrophilum* pathogenicity. This collection of *F. psychrophilum* isolates represents a valuable resource for future research on virulence-related factors through comparative genomic analyses and gene expression studies, and the collection of lytic phages with the ability to infect and kill highly virulent antibiotic-resistant variants provides a necessary tool kit for continued research aimed at estimating the potential of phage-based control of BCWD.

## NUCLEOTIDE SEQUENCE ACCESSION NUMBERS

Nucleotide sequences have been deposited in the GenBank database under accession numbers MK776535-MK776667.

## DATA AVAILABILITY

The datasets generated for this study can be found in GenBank and **Supplementary Figure S1**, MK776535-MK776667.

## ETHICS STATEMENT

All animal experiments were performed in Finland under project (ESAVI/4225/04.10.07/2017) and personal license issued by the Animal Experimental Board (Eläinkoelautakunta, ELLA).

## REFERENCES

* Avendaño-Herrera, R., Houel, A., Irgang, R., Bernardét, J. F., Godoy, M., Nicolas, P., et al. (2014). Introduction, expansion and coexistence of epidemic *Flavobacterium psychrophilum* lineages in Chilean fish farms. *Vet. Microbiol.* 17, 298-306. doi: 10.1016/j.vetmic.2014.02.009

* Beceiro, A., Tomás, M., and Bou, G. (2013). Antimicrobial resistance and virulence: a successful or deleterious association in the bacterial world? *Clin. Microbiol. Rev.* 26, 185-230. doi: 10.1128/CMR.00059-12

* Bertolini, J. M., Wakabayashi, H., Watral, V. G., Whipple, M. J., and Rohovec, J. S. (1994). Electrophoretic detection of proteases from selected strains of

## AUTHOR CONTRIBUTIONS

KS: animal experimentation, data preparation and interpretation, and writing the bulk of the manuscript. LL: bacterial isolation, DNA extraction, PCR and phenotypic characterization. PN: MLST data analysis, gene sequence submission and participation in data interpretation. JJ and DC: bacterial and phage isolation, phage-susceptibility testing and data preparation. ID, VD, and LM: bacterial isolation, hemolytic activity testing and data preparation. MM: intellectual contribution throughout the study, interpretation of data, manuscript preparation, and acquisition of funding. TW: intellectual contribution throughout the study, interpretation of data, and manuscript preparation. All authors read and approved the final manuscript.

## FUNDING

BONUS FLAVOPHAGE project has received funding from the BONUS (Art 185), funded jointly by the EU, the Innovation Fund Denmark, the Academy of Finland, and the National Centre for Research and Development (NCBiR), Poland.

## ACKNOWLEDGMENTS

We would like to thank Dr. Eva Jansson (Sweden), Dr. Bettina Schletz and Dr. Stephanie Bornstein (Germany), Dr. Agnieszka Pękala-Safińska (Poland), and Dr. Vladimir Voronin (Russia) for providing their isolates.

## SUPPLEMENTARY MATERIAL

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

**FIGURE S1** | Survival plots after intramuscular challenge with 25 *F. psychrophilum* isolates determined for the median lethal dose (LD₅₀) value in juvenile rainbow trout (*Oncorhynchus mykiss*). Each isolate was injected in five serial 10-fold doses (CFU) into 7 fish per each dose. Control fish were mock-infected with 0.5% NaCl. The isolates were categorized into highly (LD₅₀ < 10⁵), moderately (LD₅₀ = 10⁵-10⁶), or weakly virulent (LD₅₀ > 10⁶) based on the determined LD₅₀ value. The degree of virulence is enclosed in parentheses after the isolate code on each plot.

* *Flexibacter psychrophilus* and assessment of their variability. *J. Aquat. Anim. Health* 6, 224-233. doi: 10.1577/1548-8667(1994)006<0224:edopfs>2.3.co;2

* Castillo, D., Christiansen, R. H., Dalsgaard, I., Madsen, L., Espejo, R., and Middelboe, M. (2016). Comparative genome analysis provides insights into the pathogenicity of *Flavobacterium psychrophilum*. *PLoS One* 11:e0152515. doi: 10.1371/journal.pone.0152515

* Castillo, D., Christiansen, R. H., Dalsgaard, I., Madsen, L., and Middelboe, M. (2015). Bacteriophage resistance mechanisms in the fish pathogen *Flavobacterium psychrophilum*: linking genomic mutations to changes in bacterial virulence factors. *Appl. Environ. Microbiol.* 81, 1157-1167. doi: 10.1128/AEM.03699-14
Castillo, D., Higuera, G., Villa, M., Middelboe, M., Dalsgaard, I., Madsen, L., et al. (2012). Diversity of *Flavobacterium psychrophilum* and the potential use of its phages for protection against bacterial cold water disease in salmonids. *J. Fish Dis.* 35, 193-201. doi: 10.1111/j.1365-2761.2011.01336.x

Castillo, D., and Middelboe, M. (2016). Genomic diversity of bacteriophages infecting the fish pathogen *Flavobacterium psychrophilum*. *FEMS Microbiol. Lett.* 363:272. doi: 10.1093/femsle/fnw272

Cutler, J. E. (1991). Putative virulence factors of *Candida albicans*. *Annu. Rev. Microbiol.* 45, 187-218. doi: 10.1146/annurev.micro.45.1.187

Dalsgaard, I., and Madsen, L. (2000). Bacterial pathogens in rainbow trout, *Oncorhynchus mykiss* (Walbaum), reared at Danish freshwater farms. *J. Fish Dis.* 23, 199-209. doi: 10.1046/j.1365-2761.2000.00242.x

Duchaud, E., Boussaha, M., Loux, V., Bernardet, J. F., Michel, C., Kerouault, B., et al. (2007). Complete genome sequence of the fish pathogen *Flavobacterium psychrophilum*. *Nat. Biotechnol.* 25, 763-769. doi: 10.1038/nbt1313

Duchaud, E., Rochat, T., Habib, C., Barbier, P., Loux, V., Guérin, C., et al. (2018). Genomic diversity and evolution of the fish pathogen *Flavobacterium psychrophilum*. *Front. Microbiol.* 9:138. doi: 10.3389/fmicb.2018.00138

Ekman, E., Börjeson, H., and Johansson, N. (1999). *Flavobacterium psychrophilum* in Baltic salmon *Salmo salar* brood fish and their offspring. *Dis. Aquat. Org.* 37, 159-163. doi: 10.3354/dao037159

Garcia, C., Pozet, F., and Michel, C. (2000). Standardization of experimental infection with *Flavobacterium psychrophilum*, the agent of rainbow trout *Oncorhynchus mykiss* fry syndrome. *Dis. Aquat. Org.* 42, 191-197. doi: 10.3354/dao042191

Gomez, E., Mendez, J., Cascales, D., and Guijarro, J. A. (2014). *Flavobacterium psychrophilum* vaccine development: a difficult task. *Microb. Biotech.* 7, 414-423. doi: 10.1111/1751-7915.12099

Hesami, S., Parkman, J., MacInnes, J. I., Gray, J. T., Gyles, C. L., and Lumsden, J. S. (2010). Antimicrobial susceptibility of *Flavobacterium psychrophilum* isolates from Ontario. *J. Aquat. Anim. Health* 22, 39-49. doi: 10.1577/H09-008.1

Högfors-Rönnholm, E., Norrgård, J., and Wiklund, T. (2015). Adhesion of smooth and rough phenotypes of *Flavobacterium psychrophilum* to polystyrene surfaces. *J. Fish Dis.* 38, 429-437. doi: 10.1111/jfd.12250

Högfors-Rönnholm, E., and Wiklund, T. (2010a). Hemolytic activity in *Flavobacterium psychrophilum* is a contact-dependent, two-step mechanism and differently expressed in smooth and rough phenotypes. *Microb. Pathog.* 49, 369-375. doi: 10.1016/j.micpath.2010.08.002

Högfors-Rönnholm, E., and Wiklund, T. (2010b). Phase variation in *Flavobacterium psychrophilum*: characterization of two distinct colony phenotypes. *Dis. Aquat. Org.* 90, 43-53. doi: 10.3354/dao02211

Holt, R. A., Rohevec, J. S., and Fryer, J. L. (1993). "Bacterial cold water disease," in *Bacterial Diseases of Fish*, eds V. Inglis, R. J. Roberts, and N. R. Bromage (Oxford: Blackwell Scientific Publications), 3-23.

Jarau, M., Di Natale, A., Huber, P. E., MacInnes, J. I., and Lumsden, J. S. (2018). Virulence of *Flavobacterium psychrophilum* isolates in rainbow trout *Oncorhynchus mykiss* (Walbaum). *J. Fish Dis.* 41, 1505-1514. doi: 10.1111/jfd.12861

Knupp, C., Wiens, G. D., Faisal, M., Call, D. R., Cain, K. D., Nicolas, P., et al. (2019). Large-scale analysis of *Flavobacterium psychrophilum* multilocus sequence typing genotypes recovered from North American salmonids indicates that both newly identified and recurrent clonal complexes are associated with disease. *Appl. Environ. Microbiol.* 85:e02305-18. doi: 10.1128/AEM.02305-18

Long, A., Fehringer, T. R., LaFrentz, B. R., Call, D. R., and Cain, K. D. (2014). Development of a waterborne challenge model for *Flavobacterium psychrophilum*. *FEMS Microbiol. Lett.* 359, 154-160. doi: 10.1111/1574-6968.12563

Lorenzen, E., and Olesen, N. J. (1997). Characterization of isolates of *Flavobacterium psychrophilum* associated with coldwater disease or rainbow trout fry syndrome II: serological studies. *Dis. Aquat. Org.* 31, 209-220. doi: 10.3354/dao031209

Madetoja, J., Dalsgaard, I., and Wiklund, T. (2002). Occurrence of *Flavobacterium psychrophilum* in fish farming environments. *Dis. Aquat. Org.* 52, 109-118. doi: 10.3354/dao052109

Madetoja, J., Nyman, P., and Wiklund, T. (2000). *Flavobacterium psychrophilum*, invasion into and shedding by rainbow trout *Oncorhynchus mykiss*. *Dis. Aquat. Org.* 43, 27-38. doi: 10.3354/dao043027

Madsen, L., and Dalsgaard, I. (1999). Reproducible methods for experimental infection with *Flavobacterium psychrophilum* in rainbow trout, *Oncorhynchus mykiss*. *Dis. Aquat. Org.* 36, 169-176. doi: 10.3354/dao036169

Madsen, L., and Dalsgaard, I. (2000). Comparative studies of Danish *Flavobacterium psychrophilum* isolates: ribotypes, plasmid profiles, serotypes and virulence. *J. Fish Dis.* 23, 211-218. doi: 10.1046/j.1365-2761.2000.00240.x

Mazzocco, A., Waddell, T. E., Lingohr, E., and Johnson, R. P. (2009). Enumeration of bacteriophages using the small drop plaque assay system. *Methods Mol. Biol.* 501, 81-85. doi: 10.1007/978-1-60327-164-6_9

McClelland, E. E., Bernhardt, P., and Casadevall, A. (2006). Estimating the relative contributions of virulence factors for pathogenic microbes. *Infect. Immun.* 74, 1500-1504. doi: 10.1128/iai.74.3.1500-1504.2006

Miranda, C. D., Smith, P., Rojas, R., Contreras-Lynch, S., and Vega, J. A. (2016). Antimicrobial susceptibility of *Flavobacterium psychrophilum* from Chilean salmon farms and their epidemiological cut-off values using agar dilution and disk diffusion methods. *Front. Microbiol.* 7:1880. doi: 10.3389/fmicb.2016.01880

Nematollahi, A., Decostere, A., Pasmans, F., Ducatelle, R., and Haesebrouck, F. (2003). Adhesion of high and low virulence *Flavobacterium psychrophilum* strains to isolated gill arches of rainbow trout *Oncorhynchus mykiss*. *Dis. Aquat. Org.* 55, 101-107. doi: 10.3354/dao055101

Nicolas, P., Mondot, S., Achaz, G., Bouchenot, C., Bernardét, J. F., and Duchaud, E. (2008). Population structure of the fish-pathogenic bacterium *Flavobacterium psychrophilum*. *Appl. Environ. Microbiol.* 74, 3702-3709. doi: 10.1128/AEM.00244-08

Nilsen, H., Olsen, A. B., Vaagnes, Ø, Helleberg, H., Bottolfsen, K., Skjelstad, H., et al. (2011). Systemic *Flavobacterium psychrophilum* infection in rainbow trout, *Oncorhynchus mykiss* (Walbaum), farmed in fresh and brackish water in Norway. *J. Fish Dis.* 34, 403-408. doi: 10.1111/j.1365-2761.2011.01249.x

Nilsen, H., Sundell, K., Duchaud, E., Nicolas, P., Dalsgaard, I., Madsen, L., et al. (2014). Multilocus sequence typing identifies epidemic clones of *Flavobacterium psychrophilum* in Nordic countries. *Appl. Environ. Microbiol.* 80, 2728-2736. doi: 10.1128/AEM.04233-13

Organisation for Economic Cooperation and Development [OECD] (2019). *Test No. 203: Fish, Acute Toxicity Test*. OECD Guidelines for the Testing of Chemicals. Section 2. Paris: OECD Publishing.

Ostland, V. E., Byrne, P. J., Hoover, G., and Ferguson, H. W. (2000). Necrotic myositis of rainbow trout, *Oncorhynchus mykiss* (Walbaum): proteolytic characteristics of a crude extracellular preparation from *Flavobacterium psychrophilum*. *J. Fish Dis.* 23, 329-336. doi: 10.1046/j.1365-2761.2000.00251.x

Pérez-Pascual, D., Gómez, E., and Guijarro, J. A. (2015). Lack of a type-2 glycosyltransferase in the fish pathogen *Flavobacterium psychrophilum* determines pleiotropic changes and loss of virulence. *Vet. Res.* 46:1. doi: 10.1186/s13567-014-0124-5

Pérez-Pascual, D., Rochat, T., Kerouault, B., Gómez, E., Neulat-Ripoll, F., Henry, C., et al. (2017). More than gliding: involvement of GldD and GldG in the virulence of *Flavobacterium psychrophilum*. *Front. Microbiol.* 8:2168. doi: 10.3389/fmicb.2017.02168

Pryce, G., and Garcia-Granero, M. (2002). *Breusch-Pagan and Koenker Test Macro*. Available at: http://www.sps tools.net/Syntax/RegressionRepeatedMeasure/BreuschPaganAndKoenkerTest.txt (accessed March 28, 2019).

Rangdale, R. E. (1995). *Studies on rainbow trout fry syndrome (RTFS)*. PhD thesis, Ministry of Agriculture, Fisheries and Food. Weymouth: University of Stirling.

Reed, L. J., and Muench, H. (1938). A simple method of estimating fifty percent endpoints. *Am. J. Hygiene* 27, 493-497. doi: 10.1093/oxfordjournals.aje.a118408

Rochat, T., Fujiwara-Nagata, E., Calvez, S., Dalsgaard, I., Madsen, L., Calteau, A., et al. (2017). Genomic characterization of *Flavobacterium psychrophilum* serotypes and development of a multiplex PCR-based serotyping scheme. *Front. Microbiol.* 8:1752. doi: 10.3389/fmicb.2017.01752

Rochat, T., Pérez-Pascual, D., Nilsen, H., Carpentier, M., Bridel, S., Bernardet, J.-F., et al. (2019). Identification of a novel elastin-degrading enzyme from the fish pathogen *Flavobacterium psychrophilum*. *Appl. Environ. Microbiol.* 85:e02535-18. doi: 10.1128/AEM.02535-18

Shah, S. Q. A., Nilsen, H., Bottolfsen, K., Colquhoun, D. J., and Serum, H. (2012). DNA gyrase and topoisomerase IV mutations in quinolone-resistant *Flavobacterium psychrophilum* isolated from diseased salmonids
in Norway. *Microb. Drug. Resist.* 18, 207–214. doi: 10.1089/mdr.2011.0142

Siekoula-Nguedia, C., Blanc, G., Duchaud, E., and Calvez, S. (2012). Genetic diversity of *Flavobacterium psychrophilum* isolated from rainbow trout in France: predominance of a clonal complex. *Vet. Microbiol.* 161, 169–178. doi: 10.1016/j.vetmic.2012.07.022

Söderlund, R., Hakhverdyan, M., Aspán, A., and Jansson, E. (2018). Genome analysis provides insights into the epidemiology of infection with *Flavobacterium psychrophilum* among farmed salmonid fish in Sweden. *Microb. Gen.* 4:12. doi: 10.1099/mgen.0.00024

Stenholm, A. R., Dalsgaard, I., and Middelboe, M. (2008). Isolation and characterization of bacteriophages infecting the fish pathogen *Flavobacterium psychrophilum*. *Appl. Environ. Microbiol.* 74, 4070–4078. doi: 10.1128/AEM.00428-08

Strepparava, N., Nicolas, P., Wahli, T., Segner, H., and Petrini, O. (2013). Molecular epidemiology of *Flavobacterium psychrophilum* from Swiss fish farms. *Dis. Aquat. Org.* 105, 203–210. doi: 10.3354/dao02609

Sudheesh, P. S., Al-Ghabshi, A., Al-Mazrooei, N., and Al-Habsi, S. (2012). Comparative pathogenomics of bacteria causing infectious diseases in fish. *Int. J. Evol. Biol.* 2012, 1–16. doi: 10.1155/2012/457264

Sundell, K., Heinikainen, S., and Wiklund, T. (2013). Structure of *Flavobacterium psychrophilum* populations infecting farmed rainbow trout (*Oncorhynchus mykiss*). *Dis. Aquat. Org.* 103, 111–119. doi: 10.3354/dao02573

Sundell, K., and Wiklund, T. (2015). Characteristics of epidemic and sporadic *Flavobacterium psychrophilum* sequence types. *Aquaculture* 441, 51–56. doi: 10.1016/j.aquaculture.2015.02.010

Toyama, T., Kita-Tsukamoto, K., and Wakabayashi, H. (1994). Identification of *Cytophaga psychrophila* by PCR targeted 16S ribosomal RNA. *Fish Pathol.* 29, 271–275. doi: 10.3147/jsfp.29.271

Van Vliet, D., Wiens, G. D., Loch, T. P., Nicolas, P., and Faisal, M. (2016). Genetic diversity of *Flavobacterium psychrophilum* isolates from three *Oncorhynchus* spp. in the United States, as revealed by multilocus sequence typing. *Appl. Environ. Microbiol.* 82, 3246–3255. doi: 10.1128/AEM.00411-16

Yeung, A. T., Bains, M., and Hancock, R. E. (2011). The sensor kinase CbrA is a global regulator that modulates metabolism, virulence, and antibiotic resistance in *Pseudomonas aeruginosa*. *J. Bacteriol.* 193, 918–931. doi: 10.1128/JB.00911-10

**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 © 2019 Sundell, Landor, Nicolas, Jørgensen, Castillo, Middelboe, Dalsgaard, Donati, Madsen and Wiklund. 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.