Nineteen Vibrio anguillarum-specific temperate bacteriophages isolated across Europe and Chile from aquaculture and environmental sites were genome sequenced and analyzed for host range, morphology and life cycle characteristics. The phages were classified as Siphoviridae with genome sizes between 46,006 and 54,201 bp. All 19 phages showed high genetic similarity, and 13 phages were genetically identical. Apart from sporadically distributed single nucleotide polymorphisms (SNPs), genetic diversifications were located in three variable regions (VR1, VR2 and VR3) in six of the phage genomes. Identification of specific genes, such as N6-adenine methyltransferase and lambda like repressor, as well as the presence of a tRNA^Arg, suggested a both mutualistic and parasitic interaction between phages and hosts. During short term phage exposure experiments, 28% of a V. anguillarum host population was lysogenized by the temperate phages and a genomic analysis of a collection of 31 virulent V. anguillarum showed that the isolated phages were present as prophages in >50% of the strains covering large geographical distances. Further, phage sequences were widely distributed among CRISPR-Cas arrays of publicly available sequenced Vibrios. The observed distribution of these specific temperate Vibriophages across large geographical scales may be explained by efficient dispersal of phages and bacteria in the marine environment combined with a mutualistic interaction between temperate phages and their hosts which selects for co-existence rather than arms race dynamics.

The fish pathogen Flavobacterium psychrophilum is currently one of the main pathogenic bacteria hampering the productivity of salmonid farming worldwide. Although putative virulence determinants have been identified, the genetic basis for variation in virulence of F. psychrophilum is not fully understood. In this study, we analyzed whole-genome sequences of a collection of 25 F. psychrophilum isolates from Baltic Sea countries and compared genomic information with a previous determination of their virulence in juvenile rainbow trout. The results revealed a conserved population of F. psychrophilum that were consistently present across the Baltic Sea countries, with no clear association between genomic repertoire, phylogenomic, or gene distribution and virulence traits. However, analysis of the entire genome of four F. psychrophilum isolates by hybrid assembly provided an unprecedented resolution for discriminating even highly related isolates. The results showed that isolates with different virulence phenotypes harbored genetic variances on a number of consecutive leucine-rich repeat (LRR) proteins, repetitive motifs in gliding motility-associated protein, and the insertion of transposable elements into intergenic and genic regions. Thus, these findings provide novel insights into the genetic variation of these elements and their putative role in the modulation of F. psychrophilum virulence.