Characterization of atypical Flavobacterium columnare and identification of a new genomovar

Flavobacterium covae is the predominant species of columnaris-causing bacteria impacting the Channel Catfish industry in the southeastern United States

OBJECTIVE

Columnaris disease is a leading cause of disease-related losses in the catfish industry of the southeastern United States. The term "columnaris-causing bacteria" (CCB) has been coined in reference to the four described species that cause columnaris disease: Flavobacterium columnare, F. covae, F. davisii, and F. oreochromis. Historically, F. columnare, F. covae, and F. davisii have been isolated from columnaris disease cases in the catfish industry; however, there is a lack of knowledge of which CCB species are most prevalent in farm-raised catfish. The current research objectives were to (1) sample columnaris disease cases from the U.S. catfish industry and identify the species of CCB involved and (2) determine the virulence of the four CCB species in Channel Catfish Ictalurus punctatus in controlled laboratory challenges.

METHODS

Bacterial isolates or swabs of external lesions from catfish were collected from 259 columnaris disease cases in Mississippi and Alabama during 2015-2019. The DNA extracted from the samples was analyzed using a CCB-specific multiplex polymerase chain reaction to identify the CCB present in each diagnostic case. Channel Catfish were challenged by immersion with isolates belonging to each CCB species to determine virulence at ~28°C and 20°C.

RESULT

Flavobacterium covae was identified as the predominant CCB species impacting the U.S. catfish industry, as it was present in 94.2% (n = 244) of diagnostic case submissions. Challenge experiments demonstrated that F. covae and F. oreochromis were highly virulent to Channel Catfish, with most isolates resulting in near 100% mortality. In contrast, F. columnare and F. davisii were less virulent, with most isolates resulting in less than 40% mortality.

CONCLUSION

Collectively, these results demonstrate that F. covae is the predominant CCB in the U.S. catfish industry, and research aimed at developing new control and prevention strategies should target this bacterial species. The methods described herein can be used to continue monitoring the prevalence of CCB in the catfish industry and can be easily applied to other industries to identify which Flavobacterium species have the greatest impact.

Elucidating the dynamic immune responses within the ocular mucosa of rainbow trout (Oncorhynchus mykiss) after infection with Flavobacterium columnare

The eye of vertebrates is constantly faced with numerous challenges from aquatic or airborne pathogens. As a crucial first line of defense, the ocular mucosa (OM) protects the visual organ from external threats in vertebrates such as birds and mammals. However, the understanding of ocular mucosal immunity in early vertebrates, such as teleost fish, remains limited, particularly concerning their resistance to bacterial infections. To gain insights into the pivotal role of the OM in antibacterial immunity among teleost fish, we developed a bacterial infection model using Flavobacterium columnare in rainbow trout (Oncorhynchus mykiss). Here the qPCR and immunofluorescence results showed that F. columnare could invade trout OM, suggesting that the OM could be a primary target and barrier for the bacteria. Moreover, immune-related genes (il-6, il-8, il-11, cxcl10, nod1, il1-b, igm, igt, etc.) were upregulated in the OM of trout following F. columnare infection, as confirmed by qPCR, which was further proved through RNA-seq. The results of transcriptome analyses showed that bacterial infection critically triggers a robust immune response, including innate, and adaptive immune-related signaling pathways such as Toll-like, NOD-like, and C-type lectin receptor signaling pathway and immune network for IgA production, which underscores the immune role of the OM in bacterial infection. Interestingly, a substantial reduction in the expression of genes associated with visual function was observed after infection, indicating that bacterial infection could impact ocular function. Overall, our findings have unveiled a robust mucosal immune response to bacterial infection in the teleost OM for the first time, providing valuable insights for future research into the mechanisms and functions of ocular mucosal immunity in early vertebrate species.

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.

Investigating the Ability of Edwardsiella ictaluri and Flavobacterium covae to Persist within Commercial Catfish Pond Sediments under Laboratory Conditions

Two prevalent bacterial diseases in catfish aquaculture are enteric septicemia of catfish and columnaris disease caused by Edwardsiella ictaluri and Flavobacterium covae, respectively. Chronic and recurring outbreaks of these bacterial pathogens result in significant economic losses for producers annually. Determining if these pathogens can persist within sediments of commercial ponds is paramount. Experimental persistence trials (PT) were conducted to evaluate the persistence of E. ictaluri and F. covae in pond sediments. Twelve test chambers containing 120 g of sterilized sediment from four commercial catfish ponds were inoculated with either E. ictaluri (S97-773) or F. covae (ALG-00-530) and filled with 8 L of disinfected water. At 1, 2, 4-, 6-, 8-, and 15-days post-inoculation, 1 g of sediment was removed, and colony-forming units (CFU) were enumerated on selective media using 6 × 6 drop plate methods. E. ictaluri population peaked on Day 3 at 6.4 ± 0.5 log₁₀ CFU g^-1. Correlation analysis revealed no correlation between the sediment physicochemical parameters and E. ictaluri log₁₀ CFU g^-1. However, no viable F. covae colonies were recovered after two PT attempts. Future studies to improve understanding of E. ictaluri pathogenesis and persistence, and potential F. covae persistence in pond bottom sediments are needed.

Prevalence of genetically similar Flavobacterium columnare phages across aquaculture environments reveals a strong potential for pathogen control

Intensive aquaculture conditions expose fish to bacterial infections, leading to significant financial losses, extensive antibiotic use and risk of antibiotic resistance in target bacteria. Flavobacterium columnare causes columnaris disease in aquaculture worldwide. To develop a bacteriophage‐based control of columnaris disease, we isolated and characterized 126 F. columnare strains and 63 phages against F. columnare from Finland and Sweden in 2017. Bacterial isolates were virulent on rainbow trout (Oncorhynchus mykiss) and fell into four previously described genetic groups A, C, E and G, with genetic groups C and E being the most virulent. Phage host range studied against a collection of 227 bacterial isolates (from 2013 to 2017) demonstrated modular infection patterns based on host genetic group. Phages infected contemporary and previously isolated bacterial hosts, but bacteria isolated most recently were generally resistant to previously isolated phages. Despite large differences in geographical origin, isolation year or host range of the phages, whole‐genome sequencing of 56 phages showed high level of genetic similarity to previously isolated F. columnare phages (Ficleduovirus, Myoviridae). Altogether, this phage collection demonstrates a potential for use in phage therapy.

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.

Growth of Flavobacterium columnare Genomovars in the Presence or Lack of Supplemental Cations

Flavobacterium columnare is a problematic pathogen for the aquaculture industry where isolates are classified by genomovars. Suspended growth in a low nutrient media, like tryptone yeast extract salts, is a common method used for laboratory study. The presence of calcium and magnesium is the factor contributing to growth, virulence, and biofilm formation for F. columnare. Exponential growth occurs within 24 h for F. columnare when grown in complete tryptone yeast extract salts medium at 30 ºC. Withholding CaCl₂ and MgSO₄ components from a complete TYES formulation reduced or completely inhibited growth of genomovar I isolates but not the growth of genomovar II, IIB, or III isolates. Only 3 of 20 genomovar I isolates, MS-FC-4, FC-CSF-53, and 023-08-3, could achieve O.D. 540 readings ≥ 0.3 but only after 48-h incubation in cation-restricted TYES. Independently adding CaCl₂ or MgSO₄ to tryptone and yeast extracts did not result in a genomovar-specific growth phenotype, but generally demonstrated increased clumping with individual isolates presenting abnormal growth. Clumping formed filamentous strings that migrated to the top of the culture tube when isolates were grown in TYE+CaCl₂. Several of the F. columnare isolates from all the genomovars exhibited delayed growth when a single cation source was provided. This study demonstrates phenotypic differences between and within genomovars of a single bacterial species when grown under different TYES media conditions.

Genetic characterization of Flavobacterium columnare isolates from the Pacific Northwest, USA

Flavobacterium columnare is the causative agent of columnaris disease. Previous work has demonstrated a high degree of genetic variability among F. columnare isolates, identifying 4 genetic groups (GGs) with some host associations. Herein, a total of 49 F. columnare isolates were characterized, the majority of which were collected from 15 different locations throughout the US Pacific Northwest. Most isolates were collected from 2015-2018 and originated from disease outbreaks in salmonid hatcheries and rearing ponds, sturgeon hatcheries and ornamental fish. Other isolates were part of collections recovered from 1980-2018. Initial identification was confirmed by F. columnare species-specific qPCR. Study isolates were further characterized using a multiplex PCR that differentiates between the 4 currently recognized F. columnare GGs. Multiplex PCR results were supported by repetitive sequence-mediated PCR fingerprinting and gyrB sequence analysis. F. columnare GG1 was the most prevalent (83.7%, n = 41/49), represented by isolates from salmonids (n = 32), white sturgeon (n = 2), channel catfish (n = 1), ornamental goldfish (n = 1), koi (n = 3), wild sunfish (n = 1) and 1 unknown host. Six isolates (12.2%, n = 6/49) were identified as GG3, which were cultured from rainbow trout (n = 3) and steelhead trout (n = 3). Two isolates were identified as GG2 (4.1%, n = 2/49) and were from ornamental fish. No GG4 isolates were cultured in this study. The biological significance of this genetic variability remains unclear, but this variation could have significant implications for fish health management. The results from this study provide baseline data for future work developing strategies to ameliorate columnaris-related losses in the US Pacific Northwest.

Virulence variations of Flavobacterium columnare in rainbow trout (Oncorhynchus mykiss) eyed eggs and alevin

Flavobacterium columnare (Fc) is the causative agent for columnaris disease (CD) in several fish species and an emerging problem for rainbow trout aquaculture. We characterize the virulence phenotype of two Fc isolates, CSF-298-10 and MS-FC-4, against trout from two sources, NCCCWA and a production stock (PS), at the eyed egg and alevin life stages. Immersion challenges demonstrated that NCCCWA eyed eggs were susceptible to the Fc isolate MS-FC-4 (>97% mortality) but no mortality was observed against PS eyed eggs. The CSF-298-10 had little effect on any eyed eggs tested and was not highly virulent to any alevin till day six post-hatch, up to 38% for NCCCWA and ~80% PS alevin. The MS-FC-4 strain produced ≥80% mortality any day an immersion challenge occurred post-hatch. Significant difference in CFU counts was recorded between the Fc strains on 2 days post-hatch immersion challenges. Counts for the NCCCWA alevin were 4.4 × 10³  CFU/ml^-1 and 1.8 × 10⁶  CFU/ml^-1 for the CSF-298-10 strain and MS-FC-4 strain, respectively, and for the PS alevin CSF-298-10 measured 9.9 × 10¹  CFU/ml^-1 and 3.8 × 10⁵  CFU/ml^-1 for MS-FC-4. These two Fc isolates present stark differences in virulence phenotypes to both eyed eggs and alevin and present an interesting model system for virulence kinetics and potentially alternative pathogenic pathways.

Mining zebrafish microbiota reveals key community-level resistance against fish pathogen infection

The long-known resistance to pathogens provided by host-associated microbiota fostered the notion that adding protective bacteria could prevent or attenuate infection. However, the identification of endogenous or exogenous bacteria conferring such protection is often hindered by the complexity of host microbial communities. Here, we used zebrafish and the fish pathogen Flavobacterium columnare as a model system to study the determinants of microbiota-associated colonization resistance. We compared infection susceptibility in germ-free, conventional and reconventionalized larvae and showed that a consortium of 10 culturable bacterial species are sufficient to protect zebrafish. Whereas survival to F. columnare infection does not rely on host innate immunity, we used antibiotic dysbiosis to alter zebrafish microbiota composition, leading to the identification of two different protection strategies. We first identified that the bacterium Chryseobacterium massiliae individually protects both larvae and adult zebrafish. We also showed that an assembly of 9 endogenous zebrafish species that do not otherwise protect individually confer a community-level resistance to infection. Our study therefore provides a rational approach to identify key endogenous protecting bacteria and promising candidates to engineer resilient microbial communities. It also shows how direct experimental analysis of colonization resistance in low-complexity in vivo models can reveal unsuspected ecological strategies at play in microbiota-based protection against pathogens.

The role of denitrification genes in anaerobic growth and virulence of Flavobacterium columnare

AIMS

Comparative genomics analyses indicated that the Flavobacterium columnare genome has unique denitrification genes relative to Flavobacterium psychrophilum and Flavobacterium johnsoniae, including nasA (nitrate reductase), nirS (nitrite reductase), norB (nitric oxide reductase) and nosZ (nitrous oxide reductase). The current study determines the roles of nasA, nirS, norB and nosZ in anaerobic growth, nitrate reduction, biofilm formation and virulence.

METHODS AND RESULTS

Four in-frame deletion mutants in virulent F. columnare strain 94-081 were constructed by allelic exchange using pCP29 plasmid. Compared with parent strain 94-081, FcΔnasA,FcΔnirS and FcΔnosZ mutants did not grow as well anaerobically, whereas the growth of FcΔnorB strain was similar to the parent strain (FcWT). Exogenous nitrate was not significantly consumed under anaerobic conditions in FcΔnasA, FcΔnirS and FcΔnosZ compared to parent strain 94-081. Under anaerobic conditions, Fc∆nasA, Fc∆norB and Fc∆nosZ formed significantly less biofilm than the wild type strain at 24 and 96 h, but FcΔnirS was not significantly affected. The nitrite reductase mutant FcΔnirS was highly attenuated in catfish, whereas FcΔnasA, FcΔnorB and FcΔnosZ had similar virulence to FcWT.

CONCLUSIONS

These results show, for the first time, that denitrification genes enable F. columnare to grow anaerobically using nitrate as an electron acceptor, and nitrite reductase contributes to F. columnare virulence.

SIGNIFICANCE AND IMPACT OF THE STUDY

These findings indicate potential for F. columnare to grow in nitrate-rich anaerobic zones in catfish production ponds, and they suggest that a Fc∆nirS strain could be useful as a safe live vaccine if it protects catfish against columnaris disease.

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.

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.

Branchial bioenergetics dysfunction as a relevant pathophysiological mechanism in freshwater silver catfish (Rhamdia quelen) experimentally infected with Flavobacterium columnare

Flavobacterium columnare, the causative agent of columnaris disease, is a serious bacterial disease responsible for causing devastating mortality rates in several species of freshwater fish, leading to severe economic losses in the aquaculture industry. Notwithstanding the enormous impacts this disease can have, very little is known regarding the interaction between the host and bacterium in terms of the mortality rate of silver catfish (Rhamdia quelen), as well its linkage to gill energetic homeostasis. Therefore, we conducted independent experiments to evaluate the mortality rates caused by F. columnare in silver catfish, as well as whether columnaris disease impairs the enzymes of the phosphoryl transfer network in gills of silver catfish and the pathways involved in this inhibition. Experiment I revealed that clinical signs started to appear 72 h post-infection (hpi), manifesting as lethargy, skin necrosis, fin erosion and gill discoloration. Silver catfish began to die at 96 hpi, and 100% mortality was observed at 120 hpi. Experiment II revealed that creatine kinase (CK, cytosolic and mitochondrial) and pyruvate kinase (PK) activities were inhibited in silver catfish experimentally infected with F. columnare, while no significant difference was observed between experimental and control groups with respect to adenylate kinase activity. Activity of the branchial sodium-potassium pump (Na⁺, K⁺-ATPase) was inhibited while reactive oxygen species (ROS) and lipid peroxidation levels were higher in silver catfish experimentally infected with F. columnare than in the control group at 72 hpi. Based on these data, the impairment of CK activity elicited by F. columnare caused a disruption in branchial energetic balance, possibly reducing ATP availability in the gills and provoking impairment of Na⁺, K ⁺ATPase activity. The inhibition of CK and PK activities appears to be mediated by ROS overproduction and lipid peroxidation, both of which contribute to disease pathogenesis associated with branchial tissue.

Draft Genome Sequences of Flavobacterium columnare Strains ARS1 and BGFS27, Isolated from Channel Catfish (Ictalurus punctatus)

Flavobacterium columnare strain BGFS27 was isolated from an apparently healthy wild channel catfish (Ictalurus punctatus) collected from the Mobile River in 2005. F. columnare strain ARS1 was isolated from a channel catfish suffering from columnaris disease in a commercial farm in 1996. BGFS27 belongs to genomovar II (genetic group 2), while ARS1 belongs to genomovar III (genetic group 3). Here, we report the draft genome sequences of F. columnare BGFS27 and ARS1, obtained by PacBio sequencing.

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.