Jason P. Evenhuis

Independent Emergence of Yersinia ruckeri Biotype 2 in the United States and Europe

ABSTRACT Biotype 2 (BT2) variants of the bacterium Yersinia ruckeri are an increasing disease problem in U.S. and European aquaculture and have been characterized as serovar 1 isolates that lack both peritrichous flagella and secreted phospholipase activity. The emergence of this biotype has been associated with an increased frequency of enteric redmouth disease (ERM) outbreaks in previously vaccinated salmonid fish. In this study, four independent specific natural mutations that cause the loss of both motility and secreted lipase activity were identified in BT2 strains from the United States, United Kingdom, and mainland Europe. Each of these was a unique mutation in either fliR, flhA, or flhB, all of which are genes predicted to encode essential components of the flagellar secretion apparatus. Our results demonstrate the existence of independent mutations leading to the BT2 phenotype; thus, this phenotype has emerged separately at least four times. In addition, BT2 strains from the United Kingdom were shown to have the same mutant allele found in U.S. BT2 strains, suggesting a common origin of this BT2 lineage. This differentiation of distinct BT2 lineages is of critical importance for the development and validation of alternative vaccines or other treatment strategies intended for the control of BT2 strains.

Assessment of Genetic Correlation between Bacterial Cold Water Disease Resistance and Spleen Index in a Domesticated Population of Rainbow Trout: Identification of QTL on Chromosome Omy19

Selective breeding of animals for increased disease resistance is an effective strategy to reduce mortality in aquaculture. However, implementation of selective breeding programs is limited by an incomplete understanding of host resistance traits. We previously reported results of a rainbow trout selection program that demonstrated increased survival following challenge with Flavobacterium psychrophilum, the causative agent of bacterial cold water disease (BCWD). Mechanistic study of disease resistance identified a positive phenotypic correlation between post-challenge survival and spleen somatic-index (SI). Herein, we investigated the hypothesis of a genetic correlation between the two traits influenced by colocalizing QTL. We evaluated the inheritance and calculated the genetic correlation in five year-classes of odd- and even-year breeding lines. A total of 322 pedigreed families (n = 25,369 fish) were measured for disease resistance, and 251 families (n = 5,645 fish) were evaluated for SI. Spleen index was moderately heritable in both even-year (h2 = 0.56±0.18) and odd-year (h2 = 0.60±0.15) lines. A significant genetic correlation between SI and BCWD resistance was observed in the even-year line (rg = 0.45±0.20, P = 0.03) but not in the odd-year line (rg = 0.16±0.12, P = 0.19). Complex segregation analyses of the even-year line provided evidence of genes with major effect on SI, and a genome scan of a single family, 2008132, detected three significant QTL on chromosomes Omy19, 16 and 5, in addition to ten suggestive QTL. A separate chromosome scan for disease resistance in family 2008132 identified a significant BCWD QTL on Omy19 that was associated with time to death and percent survival. In family 2008132, Omy19 microsatellite alleles that associated with higher disease resistance also associated with increased spleen size raising the hypothesis that closely linked QTL contribute to the correlation between these traits. To our knowledge, this is the first estimation of spleen size heritability and evidence for genetic linkage with specific disease resistance in a teleost fish.

Evidence of major genes affecting resistance to bacterial cold water disease in rainbow trout using Bayesian methods of segregation analysis.

Bacterial cold water disease (BCWD) causes significant economic loss in salmonid aquaculture. We previously detected genetic variation for BCWD resistance in our rainbow trout population, and a family-based selection program to improve resistance was initiated at the National Center for Cool and Cold Water Aquaculture (NCCCWA). This study investigated evidence of major trait loci affecting BCWD resistance using only phenotypic data (without using genetic markers) and Bayesian methods of segregation analysis (BMSA). A total of 10,603 juvenile fish from 101 full-sib families corresponding to 3 generations (2005, 2007, and 2009 hatch years) of the NCCCWA population were challenged by intraperitoneal injection with Flavobacterium psychrophilum, the bacterium that causes BCWD. The results from single- and multiple-QTL models of BMSA suggest that 6 to 10 QTL explaining 83 to 89% of phenotypic variance with either codominant or dominant disease-resistant alleles plus polygenic effects may underlie the genetic architecture of BCWD resistance. This study also highlights the importance of polygenic background effects in the genetic variation of BCWD resistance. The polygenic heritability on the observed scale of survival status is slightly larger than that previously reported for rainbow trout BCWD resistance. These findings provide the basis for designing informative crosses for QTL mapping and carrying out genome scans for QTL affecting BCWD resistance in rainbow trout.

Rainbow trout (Oncorhynchus mykiss) resistance to columnaris disease is heritable and favorably correlated with bacterial cold water disease resistance.

Columnaris disease (CD), caused by Flavobacterium columnare, is an emerging disease affecting rainbow trout aquaculture. Objectives of this study were to 1) estimate heritability of CD resistance in a rainbow trout line (ARS-Fp-R) previously selected 4 generations for improved bacterial cold water disease (BCWD) resistance; 2) estimate genetic correlations among CD resistance, BCWD resistance, and growth to market BW; and 3) compare CD resistance among the ARS-Fp-R, ARS-Fp-S (selected 1 generation for increased BCWD susceptibility), and ARS-Fp-C (selection control) lines. Heritability of CD resistance was estimated using data from a waterborne challenge of 44 full-sib ARS-Fp-R families produced using a paternal half-sib mating design, and genetic correlations were estimated using these data and 5 generations of BCWD resistance, 9-mo BW (approximately 0.5 kg), and 12-mo BW (approximately 1.0 kg) data from 405 ARS-Fp-R full-sib families. The CD and BCWD challenges were initiated at approximately 52 and 84 d posthatch, or approximately 650 and 1,050 degree days (°C × d), respectively. Survival of ARS-Fp-R families ranged from 0 to 48% following CD challenge and heritability estimates were similar between CD (0.17 ± 0.09) and BCWD (0.18 ± 0.03) resistance, and the genetic correlation between these 2 traits was favorable (0.35 ± 0.25). Genetic correlations were small and antagonistic (-0.15 ± 0.08 to -0.19 ± 0.24) between the 2 resistance traits and 9- and 12-mo BW. Two challenges were conducted in consecutive years to compare CD resistance among ARS-Fp-R, ARS-Fp-C, and ARS-Fp-S families. In the first challenge, ARS-Fp-R families (83% survival) had greater CD resistance than ARS-Fp-C (73.5%; P = 0.02) and ARS-Fp-S (68%; P < 0.001) families, which did not differ (P = 0.16). In the second challenge, using an approximately 2.5-fold greater challenge dose, ARS-Fp-R families exhibited greater CD resistance (56% survival) than ARS-Fp-S (38% survival; P = 0.02) families. The favorable genetic correlation between CD and BCWD resistance is supported by greater CD resistance of the ARS-Fp-R line compared to the ARS-Fp-C and ARS-Fp-S lines and suggests that both traits will be improved simultaneously when selection is practiced on only 1 trait. In summary, these data indicate the feasibility of further selective breeding of the BCWD-resistant ARS-Fp-R line for increased CD resistance to produce a double pathogen-resistant line of rainbow trout.

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.

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.

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.

Genome-wide Association Studies Reveal Similar Genetic Architecture with Shared and Unique QTL 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 QTL for BCWD resistance in rainbow trout (Oncorhynchus mykiss). However, the recent availability of a 57K SNP array and a genome physical map 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 57K 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 57K 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.

Rainbow trout () resistance to columnaris disease is heritable and favorably correlated with bacterial cold water disease resistance

Columnaris disease (CD), caused by Flavobacterium columnare, is an emerging disease affecting rainbow trout aquaculture. Objectives of this study were to 1) estimate heritability of CD resistance in a rainbow trout line (ARS-Fp-R) previously selected 4 generations for improved bacterial cold water disease (BCWD) resistance; 2) estimate genetic correlations among CD resistance, BCWD resistance, and growth to market BW; and 3) compare CD resistance among the ARS-Fp-R, ARS-Fp-S (selected 1 generation for increased BCWD susceptibility), and ARS-Fp-C (selection control) lines. Heritability of CD resistance was estimated using data from a waterborne challenge of 44 full-sib ARS-Fp-R families produced using a paternal half-sib mating design, and genetic correlations were estimated using these data and 5 generations of BCWD resistance, 9-mo BW (approximately 0.5 kg), and 12-mo BW (approximately 1.0 kg) data from 405 ARS-Fp-R full-sib families. The CD and BCWD challenges were initiated at approximately 52 and 84 d posthatch, or approximately 650 and 1,050 degree days (°C × d), respectively. Survival of ARS-Fp-R families ranged from 0 to 48% following CD challenge and heritability estimates were similar between CD (0.17 ± 0.09) and BCWD (0.18 ± 0.03) resistance, and the genetic correlation between these 2 traits was favorable (0.35 ± 0.25). Genetic correlations were small and antagonistic (-0.15 ± 0.08 to -0.19 ± 0.24) between the 2 resistance traits and 9- and 12-mo BW. Two challenges were conducted in consecutive years to compare CD resistance among ARS-Fp-R, ARS-Fp-C, and ARS-Fp-S families. In the first challenge, ARS-Fp-R families (83% survival) had greater CD resistance than ARS-Fp-C (73.5%; P = 0.02) and ARS-Fp-S (68%; P < 0.001) families, which did not differ (P = 0.16). In the second challenge, using an approximately 2.5-fold greater challenge dose, ARS-Fp-R families exhibited greater CD resistance (56% survival) than ARS-Fp-S (38% survival; P = 0.02) families. The favorable genetic correlation between CD and BCWD resistance is supported by greater CD resistance of the ARS-Fp-R line compared to the ARS-Fp-C and ARS-Fp-S lines and suggests that both traits will be improved simultaneously when selection is practiced on only 1 trait. In summary, these data indicate the feasibility of further selective breeding of the BCWD-resistant ARS-Fp-R line for increased CD resistance to produce a double pathogen-resistant line of rainbow trout.

Rainbow trout ( Oncorhynchus mykiss ) resistance to columnaris disease is heritable and favorably correlated with bacterial cold water disease resistance 1 2

Columnaris disease (CD), caused by Flavobacterium columnare, is an emerging disease affecting rainbow trout aquaculture. Objectives of this study were to 1) estimate heritability of CD resistance in a rainbow trout line (ARS-Fp-R) previously selected 4 generations for improved bacterial cold water disease (BCWD) resistance; 2) estimate genetic correlations among CD resistance, BCWD resistance, and growth to market BW; and 3) compare CD resistance among the ARS-Fp-R, ARS-Fp-S (selected 1 generation for increased BCWD susceptibility), and ARS-Fp-C (selection control) lines. Heritability of CD resistance was estimated using data from a waterborne challenge of 44 full-sib ARS-Fp-R families produced using a paternal half-sib mating design, and genetic correlations were estimated using these data and 5 generations of BCWD resistance, 9-mo BW (approximately 0.5 kg), and 12-mo BW (approximately 1.0 kg) data from 405 ARS-Fp-R full-sib families. The CD and BCWD challenges were initiated at approximately 52 and 84 d posthatch, or approximately 650 and 1,050 degree days (°C × d), respectively. Survival of ARS-Fp-R families ranged from 0 to 48% following CD challenge and heritability estimates were similar between CD (0.17 ± 0.09) and BCWD (0.18 ± 0.03) resistance, and the genetic correlation between these 2 traits was favorable (0.35 ± 0.25). Genetic correlations were small and antagonistic (-0.15 ± 0.08 to -0.19 ± 0.24) between the 2 resistance traits and 9- and 12-mo BW. Two challenges were conducted in consecutive years to compare CD resistance among ARS-Fp-R, ARS-Fp-C, and ARS-Fp-S families. In the first challenge, ARS-Fp-R families (83% survival) had greater CD resistance than ARS-Fp-C (73.5%; P = 0.02) and ARS-Fp-S (68%; P < 0.001) families, which did not differ (P = 0.16). In the second challenge, using an approximately 2.5-fold greater challenge dose, ARS-Fp-R families exhibited greater CD resistance (56% survival) than ARS-Fp-S (38% survival; P = 0.02) families. The favorable genetic correlation between CD and BCWD resistance is supported by greater CD resistance of the ARS-Fp-R line compared to the ARS-Fp-C and ARS-Fp-S lines and suggests that both traits will be improved simultaneously when selection is practiced on only 1 trait. In summary, these data indicate the feasibility of further selective breeding of the BCWD-resistant ARS-Fp-R line for increased CD resistance to produce a double pathogen-resistant line of rainbow trout.