Gliding motility proteins GldJ and SprB contribute to Flavobacterium columnare virulence.

Flavobacterium columnare causes columnaris disease in fish. Columnaris disease is incompletely understood, and adequate control measures are lacking. The type IX secretion system (T9SS) is required for F. columnare gliding motility and virulence. The T9SS and gliding motility machineries share some, but not all, components. GldN (required for gliding and for secretion) and PorV (involved in secretion but not required for gliding) are both needed for virulence, implicating T9SS-mediated secretion in virulence. The role of motility in virulence is uncertain. We constructed and analyzed sprB, sprF, and gldJ mutants that were defective for motility but that maintained T9SS function to understand the role of motility in virulence. Wild-type cells moved rapidly and formed spreading colonies. In contrast, sprB and sprF deletion mutants were partially defective in gliding and formed nonspreading colonies. Both mutants exhibited reduced virulence in rainbow trout fry. A gldJ deletion mutant was nonmotile, secretion deficient, and avirulent in rainbow trout fry. To separate the roles of GldJ in secretion and in motility, we generated gldJ truncation mutants that produce nearly full-length GldJ. Mutant gldJ563, which produces GldJ truncated at amino acid 563, was defective for gliding but was competent for secretion as measured by extracellular proteolytic activity. This mutant displayed reduced virulence in rainbow trout fry, suggesting that motility contributes to virulence. Fish that survived exposure to the sprB deletion mutant or the gldJ563 mutant exhibited partial resistance to later challenge with wild-type cells. The results aid our understanding of columnaris disease and may suggest control strategies.IMPORTANCEFlavobacterium columnare causes columnaris disease in many species of freshwater fish in the wild and in aquaculture systems. Fish mortalities resulting from columnaris disease are a major problem for aquaculture. F. columnare virulence is incompletely understood, and control measures are inadequate. Gliding motility and protein secretion have been suggested to contribute to columnaris disease, but evidence directly linking motility to disease was lacking. We isolated and analyzed mutants that were competent for secretion but defective for motility. Some of these mutants exhibited decreased virulence. Fish that had been exposed to these mutants were partially protected from later exposure to the wild type. The results contribute to our understanding of columnaris disease and may aid development of control strategies.

Secreted peptidases contribute to virulence of fish pathogen Flavobacterium columnare.

Flavobacterium columnare causes columnaris disease in freshwater fish in both natural and aquaculture settings. This disease is often lethal, especially when fish population density is high, and control options such as vaccines are limited. The type IX secretion system (T9SS) is required for F. columnare virulence, but secreted virulence factors have not been fully identified. Many T9SS-secreted proteins are predicted peptidases, and peptidases are common virulence factors of other pathogens. T9SS-deficient mutants, such as ΔgldN and ΔporV, exhibit strong defects in secreted proteolytic activity. The F. columnare genome has many peptidase-encoding genes that may be involved in nutrient acquisition and/or virulence. Mutants lacking individual peptidase-encoding genes, or lacking up to ten peptidase-encoding genes, were constructed and examined for extracellular proteolytic activity, for growth defects, and for virulence in zebrafish and rainbow trout. Most of the mutants retained virulence, but a mutant lacking 10 peptidases, and a mutant lacking the single peptidase TspA exhibited decreased virulence in rainbow trout fry, suggesting that peptidases contribute to F. columnare virulence.

Flavobacterium columnare ferric iron uptake systems are required for virulence.

Flavobacterium columnare, which causes columnaris disease, is one of the costliest pathogens in the freshwater fish-farming industry. The virulence mechanisms of F. columnare are not well understood and current methods to control columnaris outbreaks are inadequate. Iron is an essential nutrient needed for metabolic processes and is often required for bacterial virulence. F. columnare produces siderophores that bind ferric iron for transport into the cell. The genes needed for siderophore production have been identified, but other components involved in F. columnare iron uptake have not been studied in detail. We identified the genes encoding the predicted secreted heme-binding protein HmuY, the outer membrane iron receptors FhuA, FhuE, and FecA, and components of an ATP binding cassette (ABC) transporter predicted to transport ferric iron across the cytoplasmic membrane. Deletion mutants were constructed and examined for growth defects under iron-limited conditions and for virulence against zebrafish and rainbow trout. Mutants with deletions in genes encoding outer membrane receptors, and ABC transporter components exhibited growth defects under iron-limited conditions. Mutants lacking multiple outer membrane receptors, the ABC transporter, or HmuY retained virulence against zebrafish and rainbow trout mirroring that exhibited by the wild type. Some mutants predicted to be deficient in multiple steps of iron uptake exhibited decreased virulence. Survivors of exposure to such mutants were partially protected against later infection by wild-type F. columnare.

Siderophores Produced by the Fish Pathogen Flavobacterium columnare Strain MS-FC-4 Are Not Essential for Its Virulence.

Flavobacterium columnare causes columnaris disease in wild and aquaculture-reared freshwater fish. F. columnare virulence mechanisms are not well understood, and current methods to control columnaris disease are inadequate. Iron acquisition from the host is important for the pathogenicity and virulence of many bacterial pathogens. F. columnare iron acquisition has not been studied in detail. We identified genes predicted to function in siderophore production for ferric iron uptake. Genes predicted to encode the proteins needed for siderophore synthesis, export, uptake, and regulation were deleted from F. columnare strain MS-FC-4. The mutants were examined for defects in siderophore production, for growth defects in iron-limited conditions, and for virulence against zebrafish and rainbow trout. Mutants lacking all siderophore activity were obtained. These mutants displayed growth defects when cultured under iron-limited conditions, but they retained virulence against zebrafish and rainbow trout similar to that exhibited by the wild type, indicating that the F. columnare MS-FC-4 siderophores are not required for virulence under the conditions tested. IMPORTANCE Columnaris disease, which is caused by Flavobacterium columnare, is a major problem for freshwater aquaculture. Little is known regarding F. columnare virulence factors, and control measures are limited. Iron acquisition mechanisms such as siderophores are important for virulence of other pathogens. We identified F. columnare siderophore biosynthesis, export, and uptake genes. Deletion of these genes eliminated siderophore production and resulted in growth defects under iron-limited conditions but did not alter virulence in rainbow trout or zebrafish. The results indicate that the F. columnare strain MS-FC-4 siderophores are not critical virulence factors under the conditions tested but may be important for survival under iron-limited conditions in natural aquatic environments or aquaculture systems.

Identification of Haplotypes Associated With Resistance to Bacterial Cold Water Disease in Rainbow Trout Using Whole-Genome Resequencing.

Bacterial cold water disease (BCWD) is an important disease in rainbow trout aquaculture. Previously, we have identified and validated two major QTL (quantitative trait loci) for BCWD resistance, located on chromosomes Omy08 and Omy25, in the odd-year Troutlodge May spawning population. We also demonstrated that marker-assisted selection (MAS) for BCWD resistance using the favorable haplotypes associated with the two major QTL is feasible. However, each favorable haplotype spans a large genomic region of 1.3-1.6 Mb. Recombination events within the haplotype regions will result in new haplotypes associated with BCWD resistance, which will reduce the accuracy of MAS for BCWD resistance over time. The objectives of this study were 1) to identify additional SNPs (single nucleotide polymorphisms) associated with BCWD resistance using whole-genome sequencing (WGS); 2) to validate the SNPs associated with BCWD resistance using family-based association mapping; 3) to refine the haplotypes associated with BCWD resistance; and 4) to evaluate MAS for BCWD resistance using the refined QTL haplotypes. Four consecutive generations of the Troutlodge May spawning population were evaluated for BCWD resistance. Parents and offspring were sequenced as individuals and in pools based on their BCWD phenotypes. Over 12 million SNPs were identified by mapping the sequences from the individuals and pools to the reference genome. SNPs with significantly different allele frequencies between the two BCWD phenotype groups were selected to develop SNP assays for family-based association mapping in three consecutive generations of the Troutlodge May spawning population. Among the 78 SNPs derived from WGS, 77 SNPs were associated with BCWD resistance in at least one of the three consecutive generations. The additional SNPs associated with BCWD resistance allowed us to reduce the physical sizes of haplotypes associated with BCWD resistance to less than 0.5 Mb. We also demonstrated that the refined QTL haplotypes can be used for MAS in the Troutlodge May spawning population. Therefore, the SNPs and haplotypes reported in this study provide additional resources for improvement of BCWD resistance in rainbow trout.

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 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 other surfaces. Using genomic and proteomic approaches, we identified proteins secreted by the T9SS. We deleted 10 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 10 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 CaCl2 and MgSO4 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 CaCl2 or MgSO4 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+CaCl2. 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.

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 × 103  CFU/ml-1 and 1.8 × 106  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 × 101  CFU/ml-1 and 3.8 × 105  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.

Recirculation versus flow-through rainbow trout laboratory Flavobacterium columnare challenge.

Flavobacterium columnare immersion challenges are affected by water-related environmental parameters and thus are difficult to reproduce. Whereas these challenges are typically conducted using flow-through systems, use of a recirculating challenge system to control environmental parameters may improve reproducibility. We compared mortality, bacterial concentration, and environmental parameters between flow-through and recirculating immersion challenge systems under laboratory conditions using 20 rainbow trout families. Despite identical dose concentration (1:75 dilution), duration of challenge, lot of fish, and temperature, average mortality in the recirculating system (42%) was lower (p < 0.01) compared to the flow-through system (77%), and there was low correlation (r = 0.24) of family mortality. Mean days to death (3.25 vs. 2.99 d) and aquaria-to-aquaria variation (9.6 vs. 10.4%) in the recirculating and flow-through systems, respectively, did not differ (p ≥ 0.30). Despite 10-fold lower water replacement rate in the recirculating (0.4 exchanges h-1) compared to flow-through system (4 exchanges h-1), differences in bacterial concentration between the 2 systems were modest (≤0.6 orders of magnitude) and inconsistent throughout the 21 d challenge. Compared to the flow-through system, dissolved oxygen during the 1 h exposure and pH were greater (p ≤ 0.02), and calcium and hardness were lower (p ≤ 0.03), in the recirculating system. Although this study was not designed to test effects of specific environmental parameters on mortality, it demonstrates that the cumulative effects of these parameters result in poor reproducibility. A recirculating immersion challenge model may be warranted to empirically identify and control environmental parameters affecting mortality and thus may serve as a more repeatable laboratory challenge model.

Whole-genome mapping of quantitative trait loci and accuracy of genomic predictions for resistance to columnaris disease in two rainbow trout breeding populations.

BACKGROUND: Columnaris disease (CD) is an emerging problem for the rainbow trout aquaculture industry in the US. The objectives of this study were to: (1) identify common genomic regions that explain a large proportion of the additive genetic variance for resistance to CD in two rainbow trout (Oncorhynchus mykiss) populations; and (2) estimate the gains in prediction accuracy when genomic information is used to evaluate the genetic potential of survival to columnaris infection in each population. METHODS: Two aquaculture populations were investigated: the National Center for Cool and Cold Water Aquaculture (NCCCWA) odd-year line and the Troutlodge, Inc., May odd-year (TLUM) nucleus breeding population. Fish that survived to 21 days post-immersion challenge were recorded as resistant. Single nucleotide polymorphism (SNP) genotypes were available for 1185 and 1137 fish from NCCCWA and TLUM, respectively. SNP effects and variances were estimated using the weighted single-step genomic best linear unbiased prediction (BLUP) for genome-wide association. Genomic regions that explained more than 1% of the additive genetic variance were considered to be associated with resistance to CD. Predictive ability was calculated in a fivefold cross-validation scheme and using a linear regression method. RESULTS: Validation on adjusted phenotypes provided a prediction accuracy close to zero, due to the binary nature of the trait. Using breeding values computed from the complete data as benchmark improved prediction accuracy of genomic models by about 40% compared to the pedigree-based BLUP. Fourteen windows located on six chromosomes were associated with resistance to CD in the NCCCWA population, of which two windows on chromosome Omy 17 jointly explained more than 10% of the additive genetic variance. Twenty-six windows located on 13 chromosomes were associated with resistance to CD in the TLUM population. Only four associated genomic regions overlapped with quantitative trait loci (QTL) between both populations. CONCLUSIONS: Our results suggest that genome-wide selection for resistance to CD in rainbow trout has greater potential than selection for a few target genomic regions that were found to be associated to resistance to CD due to the polygenic architecture of this trait, and because the QTL associated with resistance to CD are not sufficiently informative for selection decisions across populations.

Variance and covariance estimates for resistance to bacterial cold water disease and columnaris disease in two rainbow trout breeding populations1.

Family-based selective breeding can be an effective strategy for controlling diseases in aquaculture. This study aimed to estimate (co)variance components for resistance to bacterial cold water disease (BCWD) and columnaris disease (CD) in two unrelated rainbow trout nucleus breeding populations: the USDA, ARS, National Center for Cool and Cold Water Aquaculture odd-year line (ARS-Fp-R), which has been subjected to five generations of selection for improved resistance to BCWD, and the Troutlodge, Inc., May-spawning odd-year line (TLUM), which has been selected for improved growth performance but not for disease resistance. A total of 46,805 and 27,821 pedigree records were available from both populations, respectively. Between 44 and 138 families per generation and population were evaluated under controlled BCWD and CD challenges, providing 32,311 and 17,861 phenotypic records for BCWD resistance, and 13,603 and 9,413 for CD resistance, in the ARS-Fp-R and TLUM populations, respectively. A two-trait animal threshold model assuming an underlying normal distribution for the binary survival phenotypes was used to estimate (co)variance components separately for each population. Resistance to BCWD (h2 = 0.27 ± 0.04 and 0.43 ± 0.08) and CD (h2 = 0.23 ± 0.07 and 0.34 ± 0.09) was moderately heritable in the ARS-Fp-R and TLUM populations, respectively. The genetic correlation between the resistance to BCWD and CD was favorably positive in the ARS-Fp-R (0.40 ± 0.17) and TLUM (0.39 ± 0.18) populations. These findings suggest that both disease resistance traits can be improved simultaneously even if genetic selection pressure is applied to only one of the two traits.

Retrospective Evaluation of Marker-Assisted Selection for Resistance to Bacterial Cold Water Disease in Three Generations of a Commercial Rainbow Trout Breeding Population.

Bacterial cold water disease (BCWD), caused by Flavobacterium psychrophilum, is an endemic and problematic disease in rainbow trout (Oncorhynchus mykiss) aquaculture. Previously, we have identified SNPs (single nucleotide polymorphisms) associated with BCWD resistance in rainbow trout. The objectives of this study were (1) to validate the SNPs associated with BCWD resistance in a commercial breeding population; and (2) to evaluate retrospectively the accuracy of MAS (marker-assisted selection) for BCWD resistance in this commercial breeding program. Three consecutive generations of the Troutlodge May breeding population were evaluated for BCWD resistance. Based on our previous studies, a panel of 96 SNPs was selected and used to genotype the parents and ten offspring from each of the 138 full-sib families of the 2015 generation, and 37 SNPs associated with BCWD resistance were validated. Thirty-six of the validated SNPs were clustered on chromosomes Omy3, Omy8 and Omy25. Thus, at least three QTL (quantitative trait loci) for BCWD resistance were validated in the 2015 generation. Three SNPs from each QTL region were used for haplotype association analysis. Three haplotypes, Omy3TGG, Omy8GCG and Omy25CGG, were found to be associated with BCWD resistance in the 2015 generation. Retrospective analyses were then performed to evaluate the accuracy of MAS for BCWD resistance using these three favorable haplotypes. The accuracy of MAS was estimated with the Pearson correlation coefficient between the total number of favorable haplotypes in the two parents and the family BCWD survival rates. The Omy8 and Omy25 haplotypes were positively correlated with the family BCWD survival rates across all three generations. The accuracies of MAS using these two haplotypes together were consistently around 0.5, which was equal or greater than the accuracy of the conventional family-based selection in the same generation. In conclusion, we have demonstrated that MAS for BCWD resistance is feasible in this commercial rainbow trout breeding population.

Accurate genomic predictions for BCWD resistance in rainbow trout are achieved using low-density SNP panels: Evidence that long-range LD is a major contributing factor.

Previously accurate genomic predictions for Bacterial cold water disease (BCWD) resistance in rainbow trout were obtained using a medium-density single nucleotide polymorphism (SNP) array. Here, the impact of lower-density SNP panels on the accuracy of genomic predictions was investigated in a commercial rainbow trout breeding population. Using progeny performance data, the accuracy of genomic breeding values (GEBV) using 35K, 10K, 3K, 1K, 500, 300 and 200 SNP panels as well as a panel with 70 quantitative trait loci (QTL)-flanking SNP was compared. The GEBVs were estimated using the Bayesian method BayesB, single-step GBLUP (ssGBLUP) and weighted ssGBLUP (wssGBLUP). The accuracy of GEBVs remained high despite the sharp reductions in SNP density, and even with 500 SNP accuracy was higher than the pedigree-based prediction (0.50-0.56 versus 0.36). Furthermore, the prediction accuracy with the 70 QTL-flanking SNP (0.65-0.72) was similar to the panel with 35K SNP (0.65-0.71). Genomewide linkage disequilibrium (LD) analysis revealed strong LD (r2  ≥ 0.25) spanning on average over 1 Mb across the rainbow trout genome. This long-range LD likely contributed to the accurate genomic predictions with the low-density SNP panels. Population structure analysis supported the hypothesis that long-range LD in this population may be caused by admixture. Results suggest that lower-cost, low-density SNP panels can be used for implementing genomic selection for BCWD resistance in rainbow trout breeding programs.

Draft Genome Sequence of the Fish Pathogen Flavobacterium columnare Strain MS-FC-4.

Flavobacterium columnare MS-FC-4 is a highly virulent genetic group 1 (formerly genomovar I) strain isolated from rainbow trout (Oncorhynchus mykiss). The draft genome consists of three contigs totaling 3,449,277 bp with 2,811 predicted open reading frames. F. columnare MS-FC-4 is a model strain for functional genomic analyses.

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.

Similar Genetic Architecture with Shared and Unique Quantitative Trait Loci for Bacterial Cold Water Disease Resistance in Two Rainbow Trout Breeding Populations.

Bacterial cold water disease (BCWD) causes significant mortality and economic losses in salmonid aquaculture. In previous studies, we identified moderate-large effect quantitative trait loci (QTL) for BCWD resistance in rainbow trout (Oncorhynchus mykiss). However, the recent availability of a 57 K SNP array and a reference genome assembly have enabled us to conduct genome-wide association studies (GWAS) that overcome several experimental limitations from our previous work. In the current study, we conducted GWAS for BCWD resistance in two rainbow trout breeding populations using two genotyping platforms, the 57 K Affymetrix SNP array and restriction-associated DNA (RAD) sequencing. Overall, we identified 14 moderate-large effect QTL that explained up to 60.8% of the genetic variance in one of the two populations and 27.7% in the other. Four of these QTL were found in both populations explaining a substantial proportion of the variance, although major differences were also detected between the two populations. Our results confirm that BCWD resistance is controlled by the oligogenic inheritance of few moderate-large effect loci and a large-unknown number of loci each having a small effect on BCWD resistance. We detected differences in QTL number and genome location between two GWAS models (weighted single-step GBLUP and Bayes B), which highlights the utility of using different models to uncover QTL. The RAD-SNPs detected a greater number of QTL than the 57 K SNP array in one population, suggesting that the RAD-SNPs may uncover polymorphisms that are more unique and informative for the specific population in which they were discovered.

The Type IX Secretion System Is Required for Virulence of the Fish Pathogen Flavobacterium columnare.

Flavobacterium columnare, a member of the phylum Bacteroidetes, causes columnaris disease in wild and aquaculture-reared freshwater fish. The mechanisms responsible for columnaris disease are not known. Many members of the phylum Bacteroidetes use type IX secretion systems (T9SSs) to secrete enzymes, adhesins, and proteins involved in gliding motility. The F. columnare genome has all of the genes needed to encode a T9SS. gldN, which encodes a core component of the T9SS, was deleted in wild-type strains of F. columnare The F. columnare ΔgldN mutants were deficient in the secretion of several extracellular proteins and lacked gliding motility. The ΔgldN mutants exhibited reduced virulence in zebrafish, channel catfish, and rainbow trout, and complementation restored virulence. PorV is required for the secretion of a subset of proteins targeted to the T9SS. An F. columnare ΔporV mutant retained gliding motility but exhibited reduced virulence. Cell-free spent media from exponentially growing cultures of wild-type and complemented strains caused rapid mortality, but spent media from ΔgldN and ΔporV mutants did not, suggesting that soluble toxins are secreted by the T9SS.IMPORTANCE Columnaris disease, caused by F. columnare, is a major problem for freshwater aquaculture. Little is known regarding the virulence factors produced by F. columnare, and control measures are limited. Analysis of targeted gene deletion mutants revealed the importance of the type IX protein secretion system (T9SS) and of secreted toxins in F. columnare virulence. T9SSs are common in members of the phylum Bacteroidetes and likely contribute to the virulence of other animal and human pathogens.

Draft Genome Sequence of the Fish Pathogen Flavobacterium columnare Strain CSF-298-10.

We announce here the draft genome assembly of Flavobacterium columnare CSF-298-10, a strain isolated from an outbreak of columnaris disease at a commercial trout farm in Hagerman Valley, Idaho, USA. The complete genome consists of 13 contigs totaling 3,284,579 bp, with an average G+C content of 31.5% and 2,933 predicted coding genes.

Genomic selection models double the accuracy of predicted breeding values for bacterial cold water disease resistance compared to a traditional pedigree-based model in rainbow trout aquaculture.

BACKGROUND: Previously, we have shown that bacterial cold water disease (BCWD) resistance in rainbow trout can be improved using traditional family-based selection, but progress has been limited to exploiting only between-family genetic variation. Genomic selection (GS) is a new alternative that enables exploitation of within-family genetic variation. METHODS: We compared three GS models [single-step genomic best linear unbiased prediction (ssGBLUP), weighted ssGBLUP (wssGBLUP), and BayesB] to predict genomic-enabled breeding values (GEBV) for BCWD resistance in a commercial rainbow trout population, and compared the accuracy of GEBV to traditional estimates of breeding values (EBV) from a pedigree-based BLUP (P-BLUP) model. We also assessed the impact of sampling design on the accuracy of GEBV predictions. For these comparisons, we used BCWD survival phenotypes recorded on 7893 fish from 102 families, of which 1473 fish from 50 families had genotypes [57 K single nucleotide polymorphism (SNP) array]. Naïve siblings of the training fish (n = 930 testing fish) were genotyped to predict their GEBV and mated to produce 138 progeny testing families. In the following generation, 9968 progeny were phenotyped to empirically assess the accuracy of GEBV predictions made on their non-phenotyped parents. RESULTS: The accuracy of GEBV from all tested GS models were substantially higher than the P-BLUP model EBV. The highest increase in accuracy relative to the P-BLUP model was achieved with BayesB (97.2 to 108.8%), followed by wssGBLUP at iteration 2 (94.4 to 97.1%) and 3 (88.9 to 91.2%) and ssGBLUP (83.3 to 85.3%). Reducing the training sample size to n = ~1000 had no negative impact on the accuracy (0.67 to 0.72), but with n = ~500 the accuracy dropped to 0.53 to 0.61 if the training and testing fish were full-sibs, and even substantially lower, to 0.22 to 0.25, when they were not full-sibs. CONCLUSIONS: Using progeny performance data, we showed that the accuracy of genomic predictions is substantially higher than estimates obtained from the traditional pedigree-based BLUP model for BCWD resistance. Overall, we found that using a much smaller training sample size compared to similar studies in livestock, GS can substantially improve the selection accuracy and genetic gains for this trait in a commercial rainbow trout breeding population.

Virulence of Flavobacterium columnare genomovars in rainbow trout Oncorhynchus mykiss.

Flavobacterium columnare is the causative agent of columnaris disease and is responsible for significant economic losses in aquaculture. F. columnare is a Gram-negative bacterium, and 5 genetic types or genomovars have been described based on restriction fragment length polymorphism of the 16S rRNA gene. Previous research has suggested that genomovar II isolates are more virulent than genomovar I isolates to multiple species of fish, including rainbow trout Oncorhynchus mykiss. In addition, improved genotyping methods have shown that some isolates previously classified as genomovar I, and used in challenge experiments, were in fact genomovar III. Our objective was to confirm previous results with respect to genomovar II virulence, and to determine the susceptibility of rainbow trout to other genomovars. The virulence of 8 genomovar I, 4 genomovar II, 3 genomovar II-B, and 5 genomovar III isolates originating from various sources was determined through 3 independent challenges in rainbow trout using an immersion challenge model. Mean cumulative percent mortality (CPM) of ~49% for genomovar I isolates, ~1% for genomovar II, ~5% for the II-B isolates, and ~7% for the III isolates was observed. The inability of genomovar II isolates to produce mortalities in rainbow trout was unanticipated based on previous studies, but may be due to a number of factors including rainbow trout source and water chemistry. The source of fish and/or the presence of sub-optimal environment may influence the susceptibility of rainbow trout to different F. columnare genomovars.