Comparative Analysis of the Flavobacterium columnare Genomovar I and II Genomes

Genomic evidence of genetic diversity and functional evolution in Flavobacterium columnare

Flavobacterium columnare is the causative agent of columnaris disease in freshwater fish. Columnaris disease can cause heavy economic losses in aquaculture. In this study, whole-genome sequencing was used to characterize this pathogen. F. columnare isolate AH-01 had a circular chromosome and plasmid that encoded a total of 3,022 genes. Isolate GX-01 only had a circular chromosome and encoded 2,965 genes. Genomic islands, prophage regions, and CRISPR/Cas systems were identified in both genomes. Both genomes presented evidence of gene variation and horizontal transfer, both of which are the essential components of genetic diversity, genome plasticity, and functional evolution. Single-gene phylogeny and comparative genome analyses were performed to investigate the variation and evolution of this pathogen. Genetic analysis of 16S rRNA and housekeeping gene sequences significantly clustered 55 F. columnare isolates into four clades. The intragroup identity of the 16S rRNA gene exceeded 99%, while the intergroup identity was below the species delineation threshold. We discovered significant translocation, inversion, and rearrangement events that influenced local synteny within each group. Notably, the observed alignments varied considerably among all the studied groups. The core genomes of all strains with available sequences comprised 747 genes, corresponding to approximately 25% of the genome. Core genome multilocus sequence typing, genome-wide orthology and phylogenetic analyses, and average nucleotide identity suggested that the currently existing F. columnare was an assemblage of several distinct species, with levels of divergence at least equivalent to those between recognized bacterial species. The present investigation provided genomic evidence of gene variation and horizontal transfer, which were the basis of genetic diversity, genome plasticity, and functional evolution. The findings supported a proposed new taxonomic perspective on F. columnare.

Pangenome analysis provides insights into the genetic diversity, metabolic versatility, and evolution of the genus Flavobacterium

Members of the genus Flavobacterium are widely distributed and produce various polysaccharide-degrading enzymes. Many species in the genus have been isolated and characterized. However, few studies have focused on marine isolates or fish pathogens, and in-depth genomic analyses, particularly comparative analyses of isolates from different habitat types, are lacking. Here, we isolated 20 strains of the genus from various environments in South Korea and sequenced their full-length genomes. Combined with published sequence data, we examined genomic traits, evolution, environmental adaptation, and putative metabolic functions in total 187 genomes of isolated species in Flavobacterium categorized as marine, host-associated, and terrestrial including freshwater. A pangenome analysis revealed a correlation between genome size and coding or noncoding density. Flavobacterium spp. had high levels of diversity, allowing for novel gene repertories via recombination events. Defense-related genes only accounted for approximately 3% of predicted genes in all Flavobacterium genomes. While genes involved in metabolic pathways did not differ with respect to isolation source, there was substantial variation in genomic traits; in particular, the abundances of tRNAs and rRNAs were higher in the host-associdated group than in other groups. One genome in the host-associated group contained a Microviridae prophage closely related to an enterobacteria phage. The proteorhodopsin gene was only identified in four terrestrial strains isolated for this study. Furthermore, recombination events clearly influenced genomic diversity and may contribute to the response to environmental stress. These findings shed light on the high genetic variation in Flavobacterium and functional roles in diverse ecosystems as a result of their metabolic versatility. IMPORTANCE The genus Flavobacterium is a diverse group of bacteria that are found in a variety of environments. While most species of this genus are harmless and utilize organic substrates such as proteins and polysaccharides, some members may play a significant role in the cycling for organic substances within their environments. Nevertheless, little is known about the genomic dynamics and/or metabolic capacity of Flavobacterium. Here, we found that Flavobacterium species may have an open pangenome, containing a variety of diverse and novel gene repertoires. Intriguingly, we discovered that one genome (classified into host-associated group) contained a Microviridae prophage closely related to that of enterobacteria. Proteorhodopsin may be expressed under conditions of light or oxygen pressure in some strains isolated for this study. Our findings significantly contribute to the understanding of the members of the genus Flavobacterium diversity exploration and will provide a framework for the way for future ecological characterizations.

Potential of low-density genotype imputation for cost-efficient genomic selection for resistance to Flavobacterium columnare in rainbow trout (Oncorhynchus mykiss)

BACKGROUND

Flavobacterium columnare is the pathogen agent of columnaris disease, a major emerging disease that affects rainbow trout aquaculture. Selective breeding using genomic selection has potential to achieve cumulative improvement of the host resistance. However, genomic selection is expensive partly because of the cost of genotyping large numbers of animals using high-density single nucleotide polymorphism (SNP) arrays. The objective of this study was to assess the efficiency of genomic selection for resistance to F. columnare using in silico low-density (LD) panels combined with imputation. After a natural outbreak of columnaris disease, 2874 challenged fish and 469 fish from the parental generation (n = 81 parents) were genotyped with 27,907 SNPs. The efficiency of genomic prediction using LD panels was assessed for 10 panels of different densities, which were created in silico using two sampling methods, random and equally spaced. All LD panels were also imputed to the full 28K HD panel using the parental generation as the reference population, and genomic predictions were re-evaluated. The potential of prioritizing SNPs that are associated with resistance to F. columnare was also tested for the six lower-density panels.

RESULTS

The accuracies of both imputation and genomic predictions were similar with random and equally-spaced sampling of SNPs. Using LD panels of at least 3000 SNPs or lower-density panels (as low as 300 SNPs) combined with imputation resulted in accuracies that were comparable to those of the 28K HD panel and were 11% higher than the pedigree-based predictions.

CONCLUSIONS

Compared to using the commercial HD panel, LD panels combined with imputation may provide a more affordable approach to genomic prediction of breeding values, which supports a more widespread adoption of genomic selection in aquaculture breeding programmes.

Evaluating Stress-Mediated Microbial Pathogenesis in Golden Shiners, Notemigonus crysoleucas

Flavobacterium covae (columnaris) is a microbial pathogen of the Golden Shiner (Notemigonus crysoleucas), a principal bait species. We investigated the effects of density and water temperature on the survival of fish subjected to a columnaris challenge and whether flow cytometry (FCM) could be a fast and reliable method to distinguish and enumerate F. covae populations from water and fish in experimental tanks. Juvenile Golden Shiners averaging 2.62 (±0.78 S.D.) g (negative for F. covae) were used in simultaneous trials at 22°C and 28°C in two ultra-low flow-through systems: each consisting of four treatments and five replicates per treatment. Treatments were fish stocked at either 600 fish/m³ or 2,400 fish/m³ and either challenged with F. covae or not; survival was observed for 48 h after challenge. Samples of water and fish tissue were obtained for FCM enumerations and validation by qPCR. No significant differences in survival were recorded between density treatments; however, high temperature and columnaris challenge treatments showed significantly higher mortality. Bacterial enumeration (number/mL) by FCM highly correlated with bacterial counts r = 0.81 (p = 0.001) in the water samples. Higher water temperatures may have increased columnaris infections and mortality in Golden Shiners. Flow cytometry is a reliable method of enumerating F. covae from experimental tank water samples.

Secretion Systems in Gram-Negative Bacterial Fish Pathogens

Bacterial fish pathogens are one of the key challenges in the aquaculture industry, one of the fast-growing industries worldwide. These pathogens rely on arsenal of virulence factors such as toxins, adhesins, effectors and enzymes to promote colonization and infection. Translocation of virulence factors across the membrane to either the extracellular environment or directly into the host cells is performed by single or multiple dedicated secretion systems. These secretion systems are often key to the infection process. They can range from simple single-protein systems to complex injection needles made from dozens of subunits. Here, we review the different types of secretion systems in Gram-negative bacterial fish pathogens and describe their putative roles in pathogenicity. We find that the available information is fragmented and often descriptive, and hope that our overview will help researchers to more systematically learn from the similarities and differences between the virulence factors and secretion systems of the fish-pathogenic species described here.

The fish pathogen Flavobacterium columnare represents four distinct species: Flavobacterium columnare, Flavobacterium covae sp. nov., Flavobacterium davisii sp. nov. and Flavobacterium oreochromis sp. nov., and emended description of Flavobacterium columnare

Flavobacterium columnare is the causative agent of columnaris disease in freshwater fish and four discrete genetic groups exist within the species, suggesting that the species designation requires revision. The present study determined the taxonomic status of the four genetic groups of F. columnare using polyphasic and phylogenomic approaches and included five representative isolates from each genetic group (including type strain ATCC 23463^T; genetic group 1). 16S rRNA gene sequence analysis revealed genetic group 2 isolate AL-02-36^T, genetic group 3 isolate 90-106^T, and genetic group 4 isolate Costa Rica 04-02-TN^T shared less than <98.8 % sequence identity to F. columnare ATCC 23463^T. Phylogenetic analyses of 16S rRNA and gyrB genes using different methodologies demonstrated the four genetic groups formed well-supported and distinct clades within the genus Flavobacterium. The average nucleotide identity (ANI) and digital DNA-DNA hybridization (GGDC) values between F. columnare ATCC 23463^T, genetic group 2 isolate AL-02-36^T, genetic group 3 isolate 90-106^T, and genetic group 4 isolate Costa Rica 04-02-TN^T were less than 90.84% and 42.7%, respectively. Biochemical and physiological characteristics were similar among the four genetic groups; however, quantitative differences in fatty acid profiles were detected and MALDI-TOF analyses demonstrated numerous distinguishing peaks unique to each genetic group. Chemotaxonomic, MALDI-TOF characterization and ANI/GGDC calculations afforded differentiation between the genetic groups, indicating each group is a discrete species. Herein, the names F. covae sp. nov. (AL-02-36^T), F. davisii sp. nov. (90-106^T), and F. oreochromis sp. nov. (Costa Rica 04-02-TN^T) are proposed to represent genetic groups 2, 3, and 4, respectively.

Delivering the pain: an overview of the type III secretion system with special consideration for aquatic pathogens

Gram-negative bacteria are known to subvert eukaryotic cell physiological mechanisms using a wide array of virulence factors, among which the type three-secretion system (T3SS) is often one of the most important. The T3SS constitutes a needle-like apparatus that the bacterium uses to inject a diverse set of effector proteins directly into the cytoplasm of the host cells where they can hamper the host cellular machinery for a variety of purposes. While the structure of the T3SS is somewhat conserved and well described, effector proteins are much more diverse and specific for each pathogen. The T3SS can remodel the cytoskeleton integrity to promote intracellular invasion, as well as silence specific eukaryotic cell signals, notably to hinder or elude the immune response and cause apoptosis. This is also the case in aquatic bacterial pathogens where the T3SS can often play a central role in the establishment of disease, although it remains understudied in several species of important fish pathogens, notably in Yersinia ruckeri. In the present review, we summarise what is known of the T3SS, with a special focus on aquatic pathogens and suggest some possible avenues for research including the potential to target the T3SS for the development of new anti-virulence drugs.

Type IX secretion system effectors and virulence of the model Flavobacterium columnare strain MS-FC-4

Flavobacterium columnare causes columnaris disease in wild and cultured freshwater fish and is a major problem for sustainable aquaculture worldwide. The F. columnare type IX secretion system (T9SS) secretes many proteins and is required for virulence. The T9SS component GldN is required for secretion and for gliding motility over surfaces. Genetic manipulation of F. columnare is inefficient, which has impeded identification of secreted proteins that are critical for virulence. Here we identified a virulent wild-type F. columnare strain (MS-FC-4) that is highly amenable to genetic manipulation. This facilitated isolation and characterization of two deletion mutants lacking core components of the T9SS. Deletion of gldN disrupted protein secretion and gliding motility and eliminated virulence in zebrafish and rainbow trout. Deletion of porV disrupted secretion and virulence but not motility. Both mutants exhibited decreased extracellular proteolytic, hemolytic, and chondroitin sulfate lyase activities. They also exhibited decreased biofilm formation and decreased attachment to fish fins and to other surfaces. Using genomic and proteomic approaches, we identified proteins secreted by the T9SS. We deleted ten genes encoding secreted proteins and characterized the virulence of mutants lacking individual or multiple secreted proteins. A mutant lacking two genes encoding predicted peptidases exhibited reduced virulence in rainbow trout, and mutants lacking a predicted cytolysin showed reduced virulence in zebrafish and rainbow trout. The results establish F. columnare strain MS-FC-4 as a genetically amenable model to identify virulence factors. This may aid development of measures to control columnaris disease and impact fish health and sustainable aquaculture. IMPORTANCE: Flavobacterium columnare causes columnaris disease in wild and aquaculture-reared freshwater fish and is a major problem for aquaculture. Little is known regarding the virulence factors involved in this disease and control measures are inadequate. The type IX secretion system (T9SS) secretes many proteins and is required for virulence, but the secreted virulence factors are not known. We identified a strain of F. columnare (MS-FC-4) that is well suited for genetic manipulation. The components of the T9SS and the proteins secreted by this system were identified. Deletion of core T9SS genes eliminated virulence. Genes encoding ten secreted proteins were deleted. Deletion of two peptidase-encoding genes resulted in decreased virulence in rainbow trout, and deletion of a cytolysin-encoding gene resulted in decreased virulence in rainbow trout and zebrafish. Secreted peptidases and cytolysins are likely virulence factors and are targets for the development of control measures.

Flavobacterium erciyesense sp. nov., a putative non-pathogenic fish symbiont

During a study aiming at isolation and genomic characterization of non-pathogenic fish symbionts, strain F-328^T was isolated from skin mucus of healthy-seeming rainbow trout (Oncorhynchus mykiss) and subjected to polyphasic characterization, including a comprehensive genome analysis. Strain F-328^T was Gram stain negative, none-gliding, strictly aerobic, and rod shaped. The 16S rRNA gene sequence of strain F-328^T exhibited the highest level of identity (98.9%) with F. turcicum F-339^T. The major fatty acids were iso-C_15:0, C_15:1 ω6c, and summed feature 3 (C_16:1 ω7c, C_16:1 ω6c). The polar lipid profile consisted mainly of phosphatidylethanolamine and aminolipids, while the predominant menaquinone was MK-6. The DNA G + C content and genome size of the strain were 35.3% and 3.4 Mb, respectively. Strain F-328^T shows 93.8% average nucleotide identity and 53% digital DNA-DNA hybridization identity with the closest type strain F. turcicum F-339^T. In addition, strain F-328^T was found to carry antimicrobial resistance genes, which confer resistance to several antimicrobials, including aminoglycoside, macrolides, and streptogramin. Differential phenotypic properties, together with the phylogenetic inference, demonstrate that strain F-328^T should be classified as a novel species of the genus Flavobacterium for which the name Flavobacterium erciyesense sp. nov. is proposed, with F-328^T as the type strain (= JCM 34201^T = KCTC 82261^T).

Bacteriophage Resistance Affects Flavobacterium columnare Virulence Partly via Mutations in Genes Related to Gliding Motility and the Type IX Secretion System

Increasing problems with antibiotic resistance have directed interest toward phage therapy in the aquaculture industry. However, phage resistance evolving in target bacteria is considered a challenge. To investigate how phage resistance influences the fish pathogen Flavobacterium columnare, two wild-type bacterial isolates, FCO-F2 and FCO-F9, were exposed to phages (FCO-F2 to FCOV-F2, FCOV-F5, and FCOV-F25, and FCO-F9 to FCL-2, FCOV-F13, and FCOV-F45), and resulting phenotypic and genetic changes in bacteria were analyzed. Bacterial viability first decreased in the exposure cultures but started to increase after 1 to 2 days, along with a change in colony morphology from original rhizoid to rough, leading to 98% prevalence of the rough morphotype. Twenty-four isolates (including four isolates from no-phage treatments) were further characterized for phage resistance, antibiotic susceptibility, motility, adhesion, and biofilm formation, protease activity, whole-genome sequencing, and virulence in rainbow trout fry. The rough isolates arising in phage exposure were phage resistant with low virulence, whereas rhizoid isolates maintained phage susceptibility and high virulence. Gliding motility and protease activity were also related to the phage susceptibility. Observed mutations in phage-resistant isolates were mostly located in genes encoding the type IX secretion system, a component of the Bacteroidetes gliding motility machinery. However, not all phage-resistant isolates had mutations, indicating that phage resistance in F. columnare is a multifactorial process, including both genetic mutations and changes in gene expression. Phage resistance may not, however, be a challenge for development of phage therapy against F. columnare infections since phage resistance is associated with decreases in bacterial virulence. IMPORTANCE Phage resistance of infectious bacteria is a common phenomenon posing challenges for the development of phage therapy. Along with a growing world population and the need for increased food production, constantly intensifying animal farming has to face increasing problems of infectious diseases. Columnaris disease, caused by Flavobacterium columnare, is a worldwide threat for salmonid fry and juvenile farming. Without antibiotic treatments, infections can lead to 100% mortality in a fish stock. Phage therapy of columnaris disease would reduce the development of antibiotic-resistant bacteria and antibiotic loads by the aquaculture industry, but phage-resistant bacterial isolates may become a risk. However, phenotypic and genetic characterization of phage-resistant F. columnare isolates in this study revealed that they are less virulent than phage-susceptible isolates and thus not a challenge for phage therapy against columnaris disease. This is valuable information for the fish farming industry globally when considering phage-based prevention and curing methods for F. columnare infections.

Biofilm Spreading by the Adhesin-Dependent Gliding Motility of Flavobacterium johnsoniae: 2. Role of Filamentous Extracellular Network and Cell-to-Cell Connections at the Biofilm Surface

Flavobacterium johnsoniae forms a thin spreading colony on nutrient-poor agar using gliding motility. As reported in the first paper, WT cells in the colony were sparsely embedded in self-produced extracellular polymeric matrix (EPM), while sprB cells were densely packed in immature biofilm with less matrix. The colony surface is critical for antibiotic resistance and cell survival. We have now developed the Grid Stamp-Peel method whereby the colony surface is attached to a TEM grid for negative-staining microscopy. The images showed that the top of the spreading convex WT colonies was covered by EPM with few interspersed cells. Cells exposed near the colony edge made head-to-tail and/or side-to-side contact and sometimes connected via thin filaments. Nonspreading sprB and gldG and gldK colonies had a more uniform upper surface covered by different EPMs including vesicles and filaments. The EPM of sprB, gldG, and WT colonies contained filaments ~2 nm and ~5 nm in diameter; gldK colonies did not include the latter. Every cell near the edge of WT colonies had one or two dark spots, while cells inside WT colonies and cells in SprB-, GldG-, or GldK-deficient colonies did not. Together, our results suggest that the colony surface structure depends on the capability to expand biofilm.

Flavobacterium bizetiae sp. nov., isolated from diseased freshwater fish in Canada at the end of the 1970s

Genome sequence analysis of two strains collected in Canada at the end of the 1970s and deposited in 1998 at the Collection de l'Institut Pasteur has led to the taxonomic description of a novel fish-associated species in the genus Flavobacterium. Both strains, CIP 105534^T and CIP 105535, were yellow-pigmented, Gram-stain-negative, non-spore-forming rod-shaped bacteria that exhibited gliding motility. They grew aerobically in a temperature range from 5 to 30 °C with optimal growth at 25 °C on trypticase soy or Reasoner's 2A agar but they did not grow on marine agar. Their major fatty acid profiles were similar, consisting of iso-C_15 : 0, C_16 : 1 ω7c and/or iso-C_15 : 0 2-OH (shown as summed feature 3), C_16 : 0 3-OH, iso-C_17 : 0 3-OH and C_16 : 0. The major polyamine was sym-homospermidine. Phosphatidylethanolamine and, most notably, ornithine-containing lipid OL2 and unidentified aminophospholipid APL1 were major polar lipids. A yellow pigment spot was visible after chromatographic analysis. The predominant respiratory quinone was MK-6. The G+C content of the two genomes was 34 mol% and their size was around 5.8 Mb. Comparison of the 16S rRNA gene sequences with those of the closely related type strains showed high levels of relatedness with Flavobacterium collinsii and Flavobacterium pectinovorum. All average nucleotide identity (ANI) and digital DNA-DNA hybridization values estimated against publicly available Flavobacterium genome assemblies were lower than 90 and 30 %, respectively. Phylogenetic, phenotypic and chemotaxonomic data indicated that the two strains represent a novel species of the genus Flavobacterium, for which the name Flavobacterium bizetiae sp. nov. is proposed. The type strain is CIP 105534^T (=LMG 1342^T). The unique ability of F. bizetiae to use melibiose as a sole source of carbon could provide a simple phenotypic test to discriminate F. bizetiae from its closest relatives.

Development of fish vaccine in Southeast Asia: A challenge for the sustainability of SE Asia aquaculture

Southeast (SE) Asia plays an important role in global food security as this region has been regarded as one of the major producers of aquaculture product and, to date, freshwater fish accounted for one-third of the total aquaculture in SE Asia. The intensification of freshwater farming corresponding to increase of consumer demands has inevitably led to the emergence and re-emergence of diseases causing tremendous economic loss in the region. Nile tilapia (Oreochromis niloticus) and striped catfish (Pangasianodon hypophthalmus), the major freshwater fish species of SE Asia, have been reported susceptible to several bacterial pathogens, e.g. Streptococcus agalactiae, Edwardsiella ictalurid and Flavobacterium columnare. Since only a limited number of vaccines being registered and marketed, these pathogenic organisms still represent a severe threat to aquaculture industry in SE Asia. However, there is profound advancement in the understanding of disease epidemiology, pathogenic mechanisms, teleost mucosal immunity and vaccine delivery system over the last few years. This review aimed to summarize those recent findings which hopefully can provide novel insight into the future development of suitable vaccine and vaccination regime against bacterial infection in SE Asia region.

Validation of a quantitative PCR assay for the detection of 2 Flavobacterium columnare genomovars

Flavobacterium columnare is the causative agent of columnaris disease in a variety of fish hosts. Using modifications to previously established protocols, a quantitative PCR (qPCR) assay was validated for the detection of 2 predominant F. columnare genomovars. The oligonucleotide primer and probe combination was designed to amplify a 203-bp region of the chondroitin AC lyase gene (GenBank AY912281) of F. columnare. There were no significant differences in amplification between genomovars. Comparable quantities of genomic DNA from 10 F. columnare strains, 5 representatives of each genomovar, produced similar results. Serial dilutions of purified PCR product demonstrated the limit of sensitivity for the assay was ~ 10 copies per reaction. The presence of gill and spleen tissue did not significantly affect the sensitivity of the assay. Comparably, bacterial DNA detected from the liver and kidney was less sensitive than pure bacterial DNA. However, detection from these tissues was within one order of magnitude of other tissues, indicating this reduction may have minimal analytic significance. This validated assay was used to approximate the minimum infectious dose for F. columnare isolate 94-081 in channel catfish and assess bacterial loads in gill and kidney tissues 48 h post-infection.

Genomic diversity in flavobacterial pathogens of aquatic origin

Flavobacterium species are considered important fish pathogens in wild and cultured fish throughout the world. They can cause acute, subacute, and chronic infections, which are mainly characterized by gill damage, skin lesions, and deep necrotic ulcerations. Primarily, three Flavobacterium species, F. branchiophilum, F. columnare, and F. psychrophilum, have been reported to cause substantial losses to freshwater fish. In this study, we evaluated genomes of 86 Flavobacterium species isolated from aquatic hosts (mainly fish) to identify their unique and shared genome features. Our results showed that F. columnare genomes cluster into four different genetic groups. In silico secretion system analysis identified that all genomes carry type I (T1SS) and type IX (T9SS) secretion systems, but the number of type I secretion system genes shows diversity between species. F. branchiophilum, F. araucananum, F. chilense, F. spartansii, and F. tructae genomes have full type VI secretion system (T6SS). F. columnare, F. hydatis, and F. plurextorum carry partial T6SS with some of the T6SS genes missing. F. columnare, F. araucananum, F. chilense, F. spartansii, F. araucananum, F. tructae, Flavobacterium sp., F. crassostreae, F. succinicans, F. hydatis, and F. plurextorum carry most of the type IV secretion system genes (T4SS). F. columnare genetic groups 1 and 2, Flavobacterium sp., and F. crassostreae encode the least number of antibiotic resistance elements. F. hydatis, F. chilense, and F. plurextorum encode the greatest number of antibiotic resistance genes. Additionally, F. spartansii, F. araucananum, and chilense encode the greatest number of virulence genes while Flavobacterium sp. and F. crassostreae encode the least number of virulence genes. In conclusion, comparative genomics of Flavobacterium species of aquatic origin will help our understanding of Flavobacterium pathogenesis.

Multiplex PCR for genotyping Flavobacterium columnare

Recent research has identified four distinct genetic groups among isolates of Flavobacterium columnare through multilocus phylogenetic analyses; however, there are no quick methods to determine the genotype of an isolate. The objective of this research was to develop a multiplex PCR to rapidly genotype F. columnare to genetic group. Comparative bacterial genomics was used to identify regions in the genomes unique to each genetic group, and primers were designed to specifically amplify different sized amplicons for each genetic group. The optimized assay was demonstrated to be specific for each genetic group and F. columnare, and no specific amplicons were generated using gDNA from a panel of other Flavobacterium spp. and bacterial fish pathogens. The analytical sensitivity of the assay ranged from 209 to 883 genome equivalents depending on the genetic group. The multiplex PCR was evaluated by genotyping a panel of 22 unknown F. columnare isolates and performing DNA sequencing of the dnaK gene in parallel. The results demonstrated 100% accordance between multiplex PCR results and assignment to genetic group via phylogenetic analysis. The multiplex PCR provides a useful tool for assigning an unknown isolate to genetic group and may be used to determine which genetic groups of F. columnare are circulating and most predominant in different aquaculture industries.

Evidence that the stress hormone cortisol regulates biofilm formation differently among Flavobacterium columnare isolates

The impact of cortisol on Flavobacterium columnare biofilm formation was explored. Firstly, the dynamics of biofilm formation by one highly (HV) and one low virulent (LV) F. columnare isolate with and without the stress hormone cortisol under microfluidic flow conditions was characterized. This to confirm that F. columnare cells could form biofilm under cortisol supplementation, and to compare the temporal and structural differences between different treatment groups. One trial revealed that in both isolates cell aggregates resembling biofilms occurred within 7-h post-inoculation. Consequently, cell clusters were sloughed away, followed by a rebuilding of bacterial cell aggregates, suggestive for a high spreading capacity. While the HV isolate revealed cell aggregates formed upstream at all time-points, for the LV isolate this was only seen upon cortisol supplementation. Secondly, the transcriptional effect of genes (gldK, gldL, gldM, gldN, sprA, sprE, sprT, and porV) belonging to the Type IX secretion system involved in gliding motility was investigated in planktonic and biofilm cells of a HV and LV isolate to which no, a low (LD) or high (HD) dose of cortisol was added. Significantly lower expression of gliding genes gldK, gldL, gldM and gldN, and of protein secretion regulator porV was seen in the LV isolate planktonic cells supplemented with a HD-cortisol. The LV isolate biofilm cells treated with the HD-cortisol showed a significant upregulation of sprT, encoding mobile surface adhesion important in bacterial colonization. This is the first evidence for the co-regulatory effect of cortisol on biofilm formation and F. columnare gliding gene expression.

Draft Genome Sequence of the Fish Pathogen Flavobacterium columnare Genomovar III Strain PH-97028 (=CIP 109753)

Flavobacterium columnare strain PH-97028 (=CIP 109753) is a genomovar III reference strain that was isolated from a diseased Ayu fish in Japan. We report here the analysis of the first available genomovar III sequence of this species to aid in identification, epidemiological tracking, and virulence studies.

Identification of Four Distinct Phylogenetic Groups in Flavobacterium columnare With Fish Host Associations

Columnaris disease, caused by the Gram-negative bacterium Flavobacterium columnare, is one of the most prevalent fish diseases worldwide. An exceptionally high level of genetic diversity among isolates of F. columnare has long been recognized, whereby six established genomovars have been described to date. However, little has been done to quantify or characterize this diversity further in a systematic fashion. The objective of this research was to perform phylogenetic analyses of 16S rRNA and housekeeping gene sequences to decipher the genetic diversity of F. columnare. Fifty isolates and/or genomes of F. columnare, originating from diverse years, geographic locations, fish hosts, and representative of the six genomovars were analyzed in this study. A multilocus phylogenetic analysis (MLPA) of the 16S rRNA and six housekeeping genes supported four distinct F. columnare genetic groups. There were associations between genomovar and genetic group, but these relationships were imperfect indicating that genomovar assignment does not accurately reflect F. columnare genetic diversity. To expand the dataset, an additional 90 16S rRNA gene sequences were retrieved from GenBank and a phylogenetic analysis of this larger dataset also supported the establishment of four genetic groups. Examination of isolate historical data indicated biological relevance to the identified genetic diversity, with some genetic groups isolated preferentially from specific fish species or families. It is proposed that F. columnare isolates be assigned to the four genetic groups defined in this study rather than genomovar in order to facilitate a standard nomenclature across the scientific community. An increased understanding of which genetic groups are most prevalent in different regions and/or aquaculture industries may allow for the development of improved targeted control and treatment measures for columnaris disease.

Genomic Diversity and Evolution of the Fish Pathogen Flavobacterium psychrophilum

Flavobacterium psychrophilum, the etiological agent of rainbow trout fry syndrome and bacterial cold-water disease in salmonid fish, is currently one of the main bacterial pathogens hampering the productivity of salmonid farming worldwide. In this study, the genomic diversity of the F. psychrophilum species is analyzed using a set of 41 genomes, including 30 newly sequenced isolates. These were selected on the basis of available MLST data with the two-fold objective of maximizing the coverage of the species diversity and of allowing a focus on the main clonal complex (CC-ST10) infecting farmed rainbow trout (Oncorhynchus mykiss) worldwide. The results reveal a bacterial species harboring a limited genomic diversity both in terms of nucleotide diversity, with ~0.3% nucleotide divergence inside CDSs in pairwise genome comparisons, and in terms of gene repertoire, with the core genome accounting for ~80% of the genes in each genome. The pan-genome seems nevertheless “open” according to the scaling exponent of a power-law fitted on the rate of new gene discovery when genomes are added one-by-one. Recombination is a key component of the evolutionary process of the species as seen in the high level of apparent homoplasy in the core genome. Using a Hidden Markov Model to delineate recombination tracts in pairs of closely related genomes, the average recombination tract length was estimated to ~4.0 Kbp and the typical ratio of the contributions of recombination and mutations to nucleotide-level differentiation (r/m) was estimated to ~13. Within CC-ST10, evolutionary distances computed on non-recombined regions and comparisons between 22 isolates sampled up to 27 years apart suggest a most recent common ancestor in the second half of the nineteenth century in North America with subsequent diversification and transmission of this clonal complex coinciding with the worldwide expansion of rainbow trout farming. With the goal to promote the development of tools for the genetic manipulation of F. psychrophilum, a particular attention was also paid to plasmids. Their extraction and sequencing to completion revealed plasmid diversity that remained hidden to classical plasmid profiling due to size similarities.

The Emerging Fish Pathogen Flavobacterium spartansii Isolated from Chinook Salmon: Comparative Genome Analysis and Molecular Manipulation

Flavobacterium spartansii strain T16^T was isolated from a disease outbreak in hatchery-reared Chinook salmon (Oncorhynchus tshawytscha) fingerlings. To gain insight into its genomic content, structure and virulence pathogenesis factors, comparative genome analyses were performed using genomes from environmental and virulent Flavobacterium strains. F. spartansii shared low average nucleotide identity (ANI) to well-known fish-pathogenic flavobacteria (e.g., F. columnare, F. psychrophilum, and F. branchiophilum), indicating that it is a new and emerging fish pathogen. The genome in T16^T had a length of 5,359,952 bp, a GC-content 35.7%, and 4,422 predicted protein-coding sequences. Flavobacterium core genome analysis showed that the number of shared genes decreased with the addition of input genomes and converged at 1182 genes. At least 8 genomic islands and 5 prophages were predicted in T16^T. At least 133 virulence factors associated with virulence in pathogenic bacteria were highly conserved in F. spartansii T16^T. Furthermore, genes linked to virulence in other bacterial species (e.g., those encoding for a type IX secretion system, collagenase and hemolysin) were found in the genome of F. spartansii T16^T and were conserved in most of the analyzed pathogenic Flavobacterium. F. spartansii was resistant to ampicillin and penicillin, consistent with the presence of multiple genes encoding diverse lactamases and the penicillin-binding protein in the genome. To allow for future investigations into F. spartansii virulence in vivo, a transposon-based random mutagenesis strategy was attempted in F. spartansii T16^T using pHimarEm1. Four putative gliding motility deficient mutants were obtained and the insertion sites of pHimarEm1 in the genome of these mutants were characterized. In total, study results clarify some of the mechanisms by which emerging flavobacterial fish pathogens may cause disease and also provide direly needed tools to investigate their pathogenesis.