Draft genome sequences of 13 Tenacibaculum maritimum isolates from farmed Norwegian Cyclopterus lumpus (lumpfish) and Scophthalmus maximus (turbot).

Thirteen bacterial isolates of Tenacibaculum maritimum were sequenced and assembled. The strains were isolated from four disease outbreaks in farmed marine fish in Norway. Eight isolates were from Cyclopterus lumpus (lumpfish), and five were from Scophthalmus maximus (turbot). Overall, sequence similarity did not correlate with host species or geographic location.

Phylogeography and host specificity of Pasteurellaceae pathogenic to sea-farmed fish in the north-east Atlantic.

The present study was undertaken to address the recent spate of pasteurellosis outbreaks among sea-farmed Atlantic salmon (Salmo salar) in Norway and Scotland, coinciding with sporadic disease episodes in lumpfish (Cyclopterus lumpus) used for delousing purposes in salmon farms. Genome assemblies from 86 bacterial isolates cultured from diseased salmon or lumpfish confirmed them all as bona fide members of the Pasteurellaceae family, with phylogenetic reconstruction dividing them into two distinct branches sharing <88% average nucleotide identity. These branches therefore constitute two separate species, namely Pasteurella skyensis and the as-yet invalidly named "Pasteurella atlantica". Both species further stratify into multiple discrete genomovars (gv.) and/or lineages, each being nearly or fully exclusive to a particular host, geographic region, and/or time period. Pasteurellosis in lumpfish is, irrespective of spatiotemporal origin, linked almost exclusively to the highly conserved "P. atlantica gv. cyclopteri" (Pac). In contrast, pasteurellosis in Norwegian sea-farmed salmon, dominated since the late-1980s by "P. atlantica gv. salmonicida" (Pas), first saw three specific lineages (Pas-1, -2, and -3) causing separate, geographically restricted, and short-lived outbreaks, before a fourth (Pas-4) emerged recently and became more widely disseminated. A similar situation involving P. skyensis (Ps) has apparently been unfolding in Scottish salmon farming since the mid-1990s, where two historic (Ps-1 and -2) and one contemporary (Ps-3) lineages have been recorded. While the epidemiology underlying all these outbreaks/epizootics remains unclear, repeated detection of 16S rRNA gene amplicons very closely related to P. skyensis and "P. atlantica" from at least five cetacean species worldwide raises the question as to whether marine mammals may play a part, possibly as reservoirs. In fact, the close relationship between the studied isolates and Phocoenobacter uteri associated with harbor porpoise (Phocoena phocoena), and their relatively distant relationship with other members of the genus Pasteurella, suggests that both P. skyensis and "P. atlantica" should be moved to the genus Phocoenobacter.

Pan-genome survey of the fish pathogen Yersinia ruckeri links accessory- and amplified genes to virulence.

While both virulent and putatively avirulent Yersinia ruckeri strains exist in aquaculture environments, the relationship between the distribution of virulence-associated factors and de facto pathogenicity in fish remains poorly understood. Pan-genome analysis of 18 complete genomes, representing established virulent and putatively avirulent lineages of Y. ruckeri, revealed the presence of a number of accessory genetic determinants. Further investigation of 68 draft genome assemblies revealed that the distribution of certain putative virulence factors correlated well with virulence and host-specificity. The inverse-autotransporter invasin locus yrIlm was, however, the only gene present in all virulent strains, while absent in lineages regarded as avirulent. Strains known to be associated with significant mortalities in salmonid aquaculture display a combination of serotype O1-LPS and yrIlm, with the well-documented highly virulent lineages, represented by MLVA clonal complexes 1 and 2, displaying duplication of the yrIlm locus. Duplication of the yrIlm locus was further found to have evolved over time in clonal complex 1, where some modern, highly virulent isolates display up to three copies.

MLVA genotyping of Moritella viscosa reveals serial emergence of novel, host-specific clonal complexes in Norwegian salmon farming.

A Multi-Locus Variable number of tandem repeat Analysis (MLVA) genotyping scheme was developed for the epidemiological study of Moritella viscosa, which causes 'winter ulcer' predominantly in sea-reared Atlantic salmon (Salmo salar L.). The assay involves multiplex PCR amplification of six Variable Number of Tandem Repeat (VNTR) loci, followed by capillary electrophoresis and data interpretation. A collection of 747 spatiotemporally diverse M. viscosa isolates from nine fish species was analysed, the majority from farmed Norwegian salmon. MLVA distributed 76% of the isolates across three major clonal complexes (CC1, CC2 and CC3), with the remaining forming minor clusters and singletons. While 90% of the salmon isolates belong to either CC1, CC2 or CC3, only 20% of the isolates recovered from other fish species do so, indicating a considerable degree of host specificity. We further highlight a series of 'clonal shifts' amongst Norwegian salmon isolates over the 35-year sampling period, with CC1 showing exclusive predominance prior to the emergence of CC2, which was later supplanted by CC3, before the recent re-emergence of CC1. Apparently, these shifts have rapidly swept the entire Norwegian coastline and conceivably, as suggested by typing of a small number of non-Norwegian isolates, the Northeast Atlantic region as a whole.

Emergence of Salmon Gill Poxvirus.

The Salmon gill poxvirus (SGPV) has emerged in recent years as the cause of an acute respiratory disease that can lead to high mortality in farmed Atlantic salmon presmolts, known as Salmon gill poxvirus disease. SGPV was first identified in Norway in the 1990s, and its large DNA genome, consisting of over 206 predicted protein-coding genes, was characterized in 2015. This review summarizes current knowledge relating to disease manifestation and its effects on the host immune system and describes dissemination of the virus. It also demonstrates how newly established molecular tools can help us to understand SGPV and its pathogenesis. Finally, we conclude and ask some burning questions that should be addressed in future research.

qPCR screening for Yersinia ruckeri clonal complex 1 against a background of putatively avirulent strains in Norwegian aquaculture.

Although a number of genetically diverse Yersinia ruckeri strains are present in Norwegian aquaculture environments, most if not all outbreaks of yersiniosis in Atlantic salmon in Norway are associated with a single specific genetic lineage of serotype O1, termed clonal complex 1. To investigate the presence and spread of virulent and putatively avirulent strains in Norwegian salmon farms, PCR assays specific for Y. ruckeri (species level) and Y. ruckeri clonal complex 1 were developed. Following extensive screening of water and biofilm, the widespread prevalence of putatively avirulent Y. ruckeri strains was confirmed in freshwater salmon hatcheries, while Y. ruckeri clonal complex 1 was found in fewer farms. The formalin-killed bacterin yersiniosis vaccine was detected in environmental samples by both PCR assays for several weeks post-vaccination. It is thus important to interpret results from recently vaccinated fish with great care. Moreover, field studies and laboratory trials confirmed that stressful management procedures may result in increased shedding of Y. ruckeri by sub-clinically infected fish. Analysis of sea water sampled throughout thermal delousing procedures proved effective for detection of Y. ruckeri in sub-clinically infected populations.

Draft Genome Sequences of 23 Tenacibaculum Isolates from Farmed Norwegian Lumpfish, Cyclopterus lumpus.

Draft genome sequences of 23 Tenacibaculum sp. strains that were isolated from Cyclopterus lumpus (lumpfish) were investigated to elucidate possible routes of transmission between Salmo salar (Atlantic salmon) and lumpfish.

Biotyping reveals loss of motility in two distinct Yersinia ruckeri lineages exclusive to Norwegian aquaculture.

Non-motile strains of Yersinia ruckeri, known as Y. ruckeri biotype 2, now dominate amongst clinical isolates retrieved from rainbow trout internationally. Due to an acute increase in the number of yersiniosis cases in Norway in recent years, followed by introduction of widespread intraperitoneal vaccination against the disease, an investigation on the prevalence of Y. ruckeri biotype 2 in Norwegian aquaculture was conducted. We biotyped 263 Y. ruckeri isolates recovered from diseased salmonids in Norway between 1985 and 2020. A total of seven biotype 2 isolates were identified, four of which were collected between 1985 and 1987, and three of which belong to the current epizootic clone, isolated from two different sea-farms in 2017. Whole-genome sequencing revealed single non-synonymous nucleotide polymorphisms in the flagellar genes flhC in isolates from the 1980s, and in fliP in isolates from 2017. In both variants, motility was restored both by complementation with wild-type alleles in trans and via spontaneous mutation-driven reversion following prolonged incubation on motility agar. While biotype 2 strains do not yet seem to have become broadly established in Norwegian aquaculture, the seven isolates described here serve to document a further two independent cases of Y. ruckeri biotype 2 emergence in salmonid aquaculture.

Tenacibaculosis in Norwegian Atlantic salmon (Salmo salar) cage-farmed in cold sea water is primarily associated with Tenacibaculum finnmarkense genomovar finnmarkense.

Skin conditions associated with Tenacibaculum spp. constitute a significant threat to the health and welfare of sea-farmed Atlantic salmon (Salmo salar L.) in Norway. Fifteen presumptive tenacibaculosis outbreaks distributed along the Norwegian coast during the late winter and spring of 2018 were investigated. Bacteriological culture confirmed the presence of Tenacibaculum spp. Seventy-six isolates cultured from individual fish were selected and subjected to whole-genome sequencing and MALDI-TOF MS analysis. Average nucleotide identity and MALDI-TOF analyses confirmed the presence of T. finnmarkense and T. dicentrarchi, with further division of T. finnmarkense into genomovars (gv.) finnmarkense and ulcerans. Core genome multilocus sequence typing (cgMLST) and single-nucleotide polymorphism (SNP) analyses identified the presence of a genetically conserved cluster of gv. finnmarkense isolates against a background of relatively genetically diverse gv. finnmarkense and gv. ulcerans isolates in 13 of the 15 studied cases. This clustering strongly suggests a link between T. finnmarkense gv. finnmarkense and development of clinical tenacibaculosis in sea-farmed Norwegian salmon in the late winter and spring. Analysis of 25 Tenacibaculum isolates collected during the spring of 2019 from similar cases identified a similar distribution of genotypes. Low water temperatures were common to all cases, and most incidences involved relatively small fish shortly after sea transfer, suggesting that these fish are particularly predisposed to Tenacibaculum infection.

Genomic Analysis of Pasteurella atlantica Provides Insight on Its Virulence Factors and Phylogeny and Highlights the Potential of Reverse Vaccinology in Aquaculture.

Pasteurellosis in farmed lumpsuckers, Cyclopterus lumpus, has emerged as a serious disease in Norwegian aquaculture in recent years. Genomic characterization of the causative agent is essential in understanding the biology of the bacteria involved and in devising an efficient preventive strategy. The genomes of two clinical Pasteurella atlantica isolates were sequenced (≈2.3 Mbp), and phylogenetic analysis confirmed their position as a novel species within the Pasteurellaceae. In silico analyses revealed 11 genomic islands and 5 prophages, highlighting the potential of mobile elements as driving forces in the evolution of this species. The previously documented pathogenicity of P. atlantica is strongly supported by the current study, and 17 target genes were recognized as putative primary drivers of pathogenicity. The expression level of a predicted vaccine target, an uncharacterized adhesin protein, was significantly increased in both broth culture and following the exposure of P. atlantica to lumpsucker head kidney leucocytes. Based on in silico and functional analyses, the strongest gene target candidates will be prioritized in future vaccine development efforts to prevent future pasteurellosis outbreaks.

Tenacibaculum piscium sp. nov., isolated from skin ulcers of sea-farmed fish, and description of Tenacibaculum finnmarkense sp. nov. with subdivision into genomovars finnmarkense and ulcerans.

Results of previous multilocus sequence and whole-genome-based analyses have suggested that a homogeneous group of isolates belonging to the genus Tenacibaculum, represented by strain TNO020T and associated with skin ulcer development in sea-farmed fish, represents an as-yet-undescribed species. Comparative whole-genome analysis performed in the present study clustered five isolates, including TNO020T, in a distinct lineage within the genus Tenacibaculum. Phenotypic differences, high intra-cluster average nucleotide identity (ANI) values and low ANI values with other Tenacibaculum species support the proposal of a novel species, for which we propose the name Tenacibaculum piscium sp. nov. with strain TNO020T (=CCUG 73833T=NCIMB 15240T) as the type strain. Further, large-scale genome analyses confirmed the existence of two different phylogenetic lineages within 'T. finnmarkense', a species effectively but not validly published previously. ANI values just above the species delineation threshold of 95-96 % confirmed that both lineages belong to the same species. This result was also supported by DNA-DNA hybridization values. Phenotypically, the two conspecific lineages are distinguishable by differences in growth temperature range and ability to degrade l-proline. For the group of isolates already commonly known as 'T. finnmarkense', we propose the name Tenacibaculum finnmarkense sp. nov., with strain TNO006T (=CCUG 73831T=NCIMB 15238T) as the type strain. We further propose the subdivision of T. finnmarkense sp. nov. into two genomovars, T. finnmarkense genomovar finnmarkense with strain TNO006T (=CCUG 73831T=NCIMB 15238T) as the type strain and T. finnmarkense genomovar ulcerans with strain TNO010T (=CCUG 73832T=NCIMB 15239T) as the type strain.

Genotyping of Salmon Gill Poxvirus Reveals One Main Predominant Lineage in Europe, Featuring Fjord- and Fish Farm-Specific Sub-Lineages.

Salmon gill poxvirus (SGPV) can cause serious gill disease in Atlantic salmon (Salmo salar L.) and represents a significant problem to aquaculture industries in Northern Europe. Here, a single-tube multi-locus variable-number tandem-repeat (VNTR) analysis (MLVA) genotyping assay, targeting eight VNTR loci, was developed for studying the epizootiology of SGPV. Through MLVA typing of SGPV positive samples from 180 farmed and wild Atlantic salmon in Northern Europe, the first molecular population study of this virus was undertaken. Comparison of resulting MLVA profiles by cluster analysis revealed considerable micro-diversity, while only a limited degree of specific clustering by country of origin could be observed, and no clustering relating to the severity of disease outbreaks. Phylogenetic analysis, based on genomic data from six SGPV specimens (three Norwegian, one Scottish, one Faroese and one Canadian), complemented and corroborated MLVA by pointing to a marked transatlantic divide in the species, with one main, relatively conserved, SGPV lineage as predominant in Europe. Within certain fjord systems and individual freshwater salmon smolt farms in Norway, however, discrete MLVA clustering patterns that prevailed over time were observed, likely reflecting local predominance of specific SGPV sub-lineages. MLVA typing was also used to refute two suspected instances of vertical SGPV transmission from salmon broodstock to offspring, and to confirm a failed disinfection attempt in one farm. These novel insights into the previously undocumented population structure of SGPV provide important clues, e.g., regarding the mechanisms underlying spread and recurrence of the virus amongst wild and farmed salmon populations, but so far no indications of more or less virulent SGPV sub-lineages have been found. The MLVA scheme represents a highly sensitive genotyping tool particularly well suited for illuminating SGPV infection routes, and adds to the relatively low number of MLVA protocols that have so far been published for viral species. Typing is reasonably inexpensive, with a moderate technological requirement, and may be completed within a single working day. Resulting MLVA profiles can be readily shared and compared across laboratories, facilitating rapid placement of samples in an international ezpizootiological context.

Pinpointing the role of Aeromonas salmonicida in the development of skin ulcerations in common dab (Limanda limanda).

Aeromonas salmonicida was isolated from ulcerations in common dab (Limanda limanda). An experiment was performed to pinpoint its role in ulceration development, considering the importance of the skin barrier and the pigmented and non-pigmented sides. The skin of dab was treated in three zones, one where scales and epidermis were removed, one where mucus was discarded and one non-treated zone. Fish were tagged to allow individual identification and challenged with A. salmonicida. Mortality and severity of the developing lesions were recorded for 21 days post-inoculation. Starting 12 days post-inoculation, mortality occurred gradually in challenged fish; however, no direct cause could be established. Both control fish and challenged fish developed ulcerations containing A. salmonicida. Sequencing of vapA gene revealed that isolates retrieved from both groups were distinct, suggesting the presence of A. salmonicida prior to the trial. Most ulcerations developed in zones where skin was removed, suggesting that abrasion might be a predisposing factor in ulceration development. Ulcerations were also observed at the insertion site of the tag, where exposed muscle tissue might have favoured the development of ulcerations. In conclusion, A. salmonicida seems to be involved in the development of skin ulcerations in dab, although the exact pathogenesis needs to be elucidated.

Multi-locus Variable-number Tandem-repeat Analysis of the Fish-pathogenic Bacterium Yersinia ruckeri by Multiplex PCR and Capillary Electrophoresis.

Yersinia ruckeri is an important pathogen of farmed salmonids worldwide, but simple tools suitable for epizootiological investigations (infection tracing, etc.) of this bacterium have been lacking. A Multi-Locus Variable-number tandem-repeat Analysis (MLVA) assay was therefore developed as an easily accessible and unambiguous tool for high-resolution genotyping of recovered isolates. For the MLVA assay presented here, DNA is extracted from cultured Y. ruckeri samples by boiling bacterial cells in water, followed by use of supernatant as template for PCR. Primer-pairs targeting ten Variable-number tandem-repeat (VNTR) loci, interspersed throughout the Y. ruckeri genome, are distributed equally amongst two five-plex PCR reactions running under identical cycling conditions. Forward primers are labelled with either of three fluorescent dyes. Following amplicon confirmation by gel electrophoresis, PCR products are diluted and subjected to capillary electrophoresis. From the resulting electropherogram profiles, peaks representing each of the VNTR loci are size-called and employed for calculating VNTR repeat counts in silico. Resulting ten-digit MLVA profiles are then used to generate Minimum spanning trees enabling epizootiological evaluation by cluster analysis. The highly portable output data, in the form of numerical MLVA profiles, can rapidly be compared across labs and placed in a spatiotemporal context. The entire procedure from cultured colony to epizootiological evaluation may be completed for up to 48 Y. ruckeri isolates within a single working day.

Biogeography of the fish pathogen Aeromonas salmonicida inferred by vapA genotyping.

A recently described typing system based on sequence variation in the virulence array protein (vapA) gene, encoding the A-layer surface protein array, allows unambiguous subtyping of Aeromonas salmonicida. In the present study, we compile A-layer typing results from a total of 675 A. salmonicida isolates, recovered over a 59-year period from 50 different fish species in 26 countries. Nine novel A-layer types (15-23) are identified, several of which display a strong predilection towards certain fish hosts, including e.g. Cyprinidae and Pleuronectidae species. Moreover, we find indications that anthropogenic transport of live fish may have aided the near global dissemination of two cyprinid-associated A-layer types. Comparison of whole genome phylogeny and A-layer typing for a subset of strains further resulted in compatible tree topologies, indicating the utility of vapA as a phylogenetic as well as an epizootiological marker in A. salmonicida. A Microreact project (microreact.org/project/r1pcOAx9m) has been created, allowing public access to the vapA analyses and relevant metadata. In sum, the results generated provide valuable insights into the global population structure of A. salmonicida, particularly in relation to its piscine host spectrum and the geographic distribution of these hosts.

Multilocus Variable-Number Tandem-Repeat Analysis of Yersinia ruckeri Confirms the Existence of Host Specificity, Geographic Endemism, and Anthropogenic Dissemination of Virulent Clones.

A multilocus variable-number tandem-repeat analysis (MLVA) assay was developed for epizootiological study of the internationally significant fish pathogen Yersinia ruckeri, which causes yersiniosis in salmonids. The assay involves amplification of 10 variable-number tandem-repeat (VNTR) loci in two five-plex PCRs, followed by capillary electrophoresis. A collection of 484 Y. ruckeri isolates, originating from various biological sources and collected from four continents over 7 decades, was analyzed. Minimum-spanning-tree cluster analysis of MLVA profiles separated the studied population into nine major clonal complexes and a number of minor clusters and singletons. The major clonal complexes could be associated with host species, geographic origin, and serotype. A single large clonal complex of serotype O1 isolates dominating the yersiniosis situation in international rainbow trout farming suggests anthropogenic spread of this clone, possibly related to transport of fish. Moreover, subclustering within this clonal complex indicates putative transmission routes and multiple biotype shift events. In contrast to the situation in rainbow trout, Y. ruckeri strains associated with disease in Atlantic salmon appear as more or less geographically isolated clonal complexes. A single complex of serotype O1 exclusive to Norway was found to be responsible for almost all major yersiniosis outbreaks in modern Norwegian salmon farming, and site-specific subclustering further indicates persistent colonization of freshwater farms in Norway. Identification of genetically diverse Y. ruckeri isolates from clinically healthy fish and environmental sources also suggests the widespread existence of less-virulent or avirulent strains.IMPORTANCE This comprehensive population study substantially improves our understanding of the epizootiological history and nature of an internationally important fish-pathogenic bacterium. The MLVA assay developed and presented represents a high-resolution typing tool particularly well suited for Yersinia ruckeri infection tracing, selection of strains for vaccine inclusion, and risk assessment. The ability of the assay to separate isolates into geographically linked and/or possibly host-specific clusters reflects its potential utility for maintenance of national biosecurity. The MLVA is internationally applicable and robust, and it provides clear, unambiguous, and easily interpreted results. Typing is reasonably inexpensive, with a moderate technological requirement, and may be completed from a harvested colony within a single working day. As the resulting MLVA profiles are readily portable, any Y. ruckeri strain may rapidly be placed in a global epizootiological context.

The repeat structure of two paralogous genes, Yersinia ruckeri invasin (yrInv) and a "Y. ruckeri invasin-like molecule", (yrIlm) sheds light on the evolution of adhesive capacities of a fish pathogen.

Inverse autotransporters comprise the recently identified type Ve secretion system and are exemplified by intimin from enterohaemorrhagic Escherichia coli and invasin from enteropathogenic Yersiniae. These proteins share a common domain architecture and promote bacterial adhesion to host cells. Here, we identified and characterized two putative inverse autotransporter genes in the fish pathogen Yersinia ruckeri NVH_3758, namely yrInv (for Y. ruckeri invasin) and yrIlm (for Y. ruckeri invasin-like molecule). When trying to clone the highly repetitive genes for structural and functional studies, we experienced problems in obtaining PCR products. PCR failures and the highly repetitive nature of inverse autotransporters prompted us to sequence the genome of Y. ruckeri NVH_3758 using PacBio sequencing, which produces some of the longest average read lengths available in the industry at this moment. According to our sequencing data, YrIlm is composed of 2603 amino acids (7812bp) and has a molecular mass of 256.4kDa. Based on the new genome information, we performed PCR analysis on four non-sequenced Y. ruckeri strains as well as the sequenced. Y. ruckeri type strain. We found that the genes are variably present in the strains, and that the length of yrIlm, when present, also varies. In addition, the length of the gene product for all strains, including the type strain, was much longer than expected based on deposited sequences. The internal repeats of the yrInv gene product are highly diverged, but represent the same bacterial immunoglobulin-like domains as in yrIlm. Using qRT-PCR, we found that yrIlm and yrInv are differentially expressed under conditions relevant for pathogenesis. In addition, we compared the genomic context of both genes in the newly sequenced Y. ruckeri strain to all available PacBio-sequenced Y. ruckeri genomes, and found indications of recent events of horizontal gene transfer. Taken together, this study demonstrates and highlights the power of Single Molecule Real-Time technology for sequencing highly repetitive proteins, and sheds light on the genetic events that gave rise to these highly repetitive genes in a commercially important fish pathogen.

Vibrio tapetis from wrasse used for ectoparasite bio-control in salmon farming: phylogenetic analysis and serotyping.

So-called 'cleaner fish', including various wrasse (Labridae) species, have become increasingly popular in Norwegian salmon farming in recent years for biocontrol of the salmon louse Lepeophtheirus salmonis. Cleaner fish mortalities in salmon farms are, however, often high. Various bacterial agents are frequently associated with episodes of increased cleaner fish mortality, and Vibrio tapetis is regularly cultured from diseased wrasse. In the present study, we investigated the genetic relationships among 54 V. tapetis isolates (34 from wrasse species) by multilocus sequence analysis (MLSA; rpoD, ftsZ, pyrH, rpoA and atpA). In the resulting phylogenetic tree, all wrasse isolates belonged to sub-clusters within V. tapetis subsp. tapetis. Slide agglutination testing further confirmed the complete dominance amongst these isolates of 4 O-antigen serotypes, designated here as V. tapetis subsp. tapetis serotypes O1, O3, O4 and O5, respectively. A pilot challenge trial using serotypes O3, O4 and O5 did not indicate high pathogenicity towards ballan wrasse Labrus bergylta, thus questioning the role of V. tapetis as a primary pathogen of this fish species.

Multilocus sequence analysis reveals extensive genetic variety within Tenacibaculum spp. associated with ulcers in sea-farmed fish in Norway.

Skin ulcer development in sea-reared salmonids, commonly associated with Tenacibaculum spp., is a significant fish welfare- and economical problem in Norwegian aquaculture. A collection of 89 Tenacibaculum isolates was subjected to multilocus sequence analysis (MLSA). The isolates were retrieved from outbreaks of clinical disease in farms spread along the Norwegian coast line from seven different fish species over a period of 19 years. MLSA analysis reveals considerable genetic diversity, but allows identification of four main clades. One clade encompasses isolates belonging to the species T. dicentrarchi, whereas three clades encompass bacteria that likely represent novel, as yet undescribed species. The study identified T. maritimum in lumpsucker, T. ovolyticum in halibut, and has extended the host and geographic range for T. soleae, isolated from wrasse. The overall lack of clonality and host specificity, with some indication of geographical range restriction argue for local epidemics involving multiple strains. The diversity of Tenacibaculum isolates from fish displaying ulcerative disease may complicate vaccine development.

Identifying copy number variation of the dominant virulence factors msa and p22 within genomes of the fish pathogen Renibacterium salmoninarum.

Renibacterium salmoninarum is the causative agent of bacterial kidney disease, an important disease of farmed and wild salmonid fish worldwide. Despite the wide spatiotemporal distribution of this disease and habitat pressures ranging from the natural environment to aquaculture and rivers to marine environments, little variation has been observed in the R. salmoninarum genome. Here we use the coverage depth from genomic sequencing corroborated by real-time quantitative PCR to detect copy number variation (CNV) among the genes of R. salmoninarum. CNV was primarily limited to the known dominant virulence factors msa and p22. Among 68 isolates representing the UK, Norway and North America, the msa gene ranged from two to five identical copies and the p22 gene ranged from one to five copies. CNV for these two genes co-occurred, suggesting they may be functionally linked. Isolates carrying CNV were phylogenetically restricted and originated predominantly from sites in North America, rather than the UK or Norway. Although both phylogenetic relationship and geographical origin were found to correlate with CNV status, geographical origin was a much stronger predictor than phylogeny, suggesting a role for local selection pressures in the repeated emergence and maintenance of this trait.