Plasma proteomic profiling of bacterial cold water disease-resistant and -susceptible rainbow trout lines and biomarker discovery

Genetic variation for disease resistance is present in salmonid fish; however, the molecular basis is poorly understood, and biomarkers of disease susceptibility/resistance are unavailable. Previously, we selected a line of rainbow trout for high survival following standardized challenge with Flavobacterium psychrophilum (Fp), the causative agent of bacterial cold water disease. The resistant line (ARS-Fp-R) exhibits over 60 percentage points higher survival compared to a reference susceptible line (ARS-Fp-S). To gain insight into the differential host response between genetic lines, we compared the plasma proteomes from day 6 after intramuscular challenge. Pooled plasma from unhandled, PBS-injected, and Fp-injected groups were simultaneously analyzed using a TMT 6-plex label, and the relative abundance of 513 proteins was determined. Data are available via ProteomeXchange, with identifier PXD041308, and the relative protein abundance values were compared to mRNA measured from a prior, whole-body RNA-seq dataset. Our results identified a subset of differentially abundant intracellular proteins was identified, including troponin and myosin, which were not transcriptionally regulated, suggesting that these proteins were released into plasma following pathogen-induced tissue damage. A separate subset of high-abundance, secreted proteins were transcriptionally regulated in infected fish. The highest differentially expressed protein was a C1q family member (designated complement C1q-like protein 3; C1q-LP3) that was upregulated over 20-fold in the infected susceptible line while only modestly upregulated, 1.8-fold, in the infected resistant line. Validation of biomarkers was performed using immunoassays and C1q-LP3, skeletal muscle troponin C, cathelcidin 2, haptoglobin, leptin, and growth and differentiation factor 15 exhibited elevated concentration in susceptible line plasma. Complement factor H-like 1 exhibited higher abundance in the resistant line compared to the susceptible line in both control and challenged fish and thus was a baseline differentiator between lines. C1q-LP3 and STNC were elevated in Atlantic salmon plasma following experimental challenge with Fp. In summary, these findings further the understanding of the differential host response to Fp and identifies salmonid biomarkers that may have use for genetic line evaluation and on-farm health monitoring.

A large-scale, multi-year microbial community survey of a freshwater trout aquaculture facility

Aquaculture is an important tool for solving the growing worldwide food demand, but infectious diseases of farmed animals represent a serious roadblock to continued industry growth. Therefore, it is essential to understand the microbial communities that reside within the built environments of aquaculture facilities to identify reservoirs of bacterial pathogens and potential correlations between commensal species and specific disease agents. Here, we present the results from 3 years of sampling a commercial rainbow trout aquaculture facility. We observed that the microbial communities residing on the abiotic surfaces within the hatchery were distinct from those residing on the surfaces at the facility's water source as well as the production raceways, despite similar communities in the water column at each location. Also, a subset of the water community seeds the biofilm communities. Lastly, we detected a common fish pathogen, Flavobacterium columnare, within the hatchery, including at the source water inlet. Importantly, the relative abundance of this pathogen was correlated with clinical disease. Our results characterized the microbial communities in an aquaculture facility, established that the hatchery environment contains a unique community composition and demonstrated that a specific fish pathogen resides within abiotic surface biofilms and is seeded from the natural water source.

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.

Aquaculture Reuse Water, Genetic Line, and Vaccination Affect Rainbow Trout (Oncorhynchus mykiss) Disease Susceptibility and Infection Dynamics

Infectious hematopoietic necrosis virus (IHNV) and Flavobacterium psychrophilum are major pathogens of farmed rainbow trout. Improved control strategies are desired but the influence of on-farm environmental factors that lead to disease outbreaks remain poorly understood. Water reuse is an important environmental factor affecting disease. Prior studies have established a replicated outdoor-tank system capable of varying the exposure to reuse water by controlling water flow from commercial trout production raceways. The goal of this research was to evaluate the effect of constant or pulsed reuse water exposure on survival, pathogen prevalence, and pathogen load. Herein, we compared two commercial lines of rainbow trout, Clear Springs Food (CSF) and Troutex (Tx) that were either vaccinated against IHNV with a DNA vaccine or sham vaccinated. Over a 27-day experimental period in constant reuse water, all fish from both lines and treatments, died while mortality in control fish in spring water was <1%. Water reuse exposure, genetic line, vaccination, and the interaction between genetic line and water exposure affected survival (P<0.05). Compared to all other water sources, fish exposed to constant reuse water had 46- to 710-fold greater risk of death (P<0.0001). Tx fish had a 2.7-fold greater risk of death compared to CSF fish in constant reuse water (P ≤ 0.001), while risk of death did not differ in spring water (P=0.98). Sham-vaccinated fish had 2.1-fold greater risk of death compared to vaccinated fish (P=0.02). Both IHNV prevalence and load were lower in vaccinated fish compared to sham-vaccinated fish, and unexpectedly, F. psychrophilum load associated with fin/gill tissues from live-sampled fish was lower in vaccinated fish compared to sham-vaccinated fish. As a result, up to forty-five percent of unvaccinated fish were naturally co-infected with F. psychrophilum and IHNV and the coinfected fish exhibited the highest IHNV loads. Under laboratory challenge conditions, co-infection with F. psychrophilum and IHNV overwhelmed IHNV vaccine-induced protection. In summary, we demonstrate that exposure to reuse water or multi-pathogen challenge can initiate complex disease dynamics that can overwhelm both vaccination and host genetic resistance.

In Vivo Experiments Provide Evidence That Flavobacterium psychrophilum Strains Belonging to Multilocus Sequence Typing Clonal Complex ST191 Are Virulent to Rainbow Trout

Flavobacterium psychrophilum, the causative agent of bacterial coldwater disease (BCWD), causes significant economic losses worldwide, particularly in farmed Rainbow Trout Oncorhynchus mykiss. Over the last decade, multilocus sequence typing has revealed >30 clonal complexes (CCs) globally, comprised of >320 F. psychrophilum sequence types (STs). Despite the large number of CCs worldwide, CC-ST10, which is currently the largest CC affecting Rainbow Trout, has been the primary focus of F. psychrophilum virulence studies, leaving the role of other CCs as primary causes of BCWD epizootics unclear. To this end, fingerling Rainbow Trout were experimentally challenged with F. psychrophilum strains belonging to the CC now recognized as the second largest in the world (CC-ST191) alongside CC-ST10 strains. Cumulative percent mortality was 100% in 7-month-old Rainbow Trout and between 27.8% and 61.1% in 8-month-old Rainbow Trout. All examined F. psychrophilum STs were virulent to Rainbow Trout, and no significant differences in virulence between CC-ST10 and CC-ST191 were detected. Due to their wide distribution and high pathogenic potential, both CC-ST191 and CC-ST10 F. psychrophilum strains are excellent candidates for further research aimed at preventing and controlling BCWD.

Structure and regulation of the NK-lysin (1-4) and NK-lysin like (a and b) antimicrobial genes in rainbow trout (Oncorhynchus mykiss)

Nk-lysin (Nkl), an antimicrobial peptide (AMP) product of natural killer cells and cytotoxic T cells in mammals, has recently been characterized in a number of finfish species. In this study, we identified six genes with sequence homology to Nkl and characterized their patterns of mRNA expression and abundances in rainbow trout (Oncorhynchus mykiss). The cDNA sequences for the six Nkls encoded precursor peptides of 128-133 aa in length, and mature peptides of 109-111 aa in length. Genomic DNA of the nkl1-4 genes consisted of five exons and four introns, whereas the nkl-like a & b genes consisted of four exons and three introns. Chromosomal locations of these peptides show that nkl1 was located on chromosome arm 25q, whereas the other five nkl genes were clustered on chromosome arm 19q. Phylogenetic analysis revealed a conserved structure of Nkls among the teleosts and further protein sequence analyses suggests that all six nkl genes fall within the Nkl sub-family of the Saposin family of proteins. Patterns of tissue-specific mRNA expression were asymmetric among the six trout Nkl homologues, with nkl1, nkl3, and nkl-like a & b occurring in immune competent organs such as spleen, gill, intestine and kidney, as well as pineal gland, brain and oocytes. However, nkl2 and nkl4, showed primary abundances in brain, pineal gland and oocyte tissues. Using mRNA sequencing, in whole-body pools of juvenile trout fry (1 g bw) exposed to Flavobacterium psychrophilum infection, we observed modest up-regulation (2-3 fold) of five (nkl 2-4 and nkl-like a & b) of the six nkl mRNAs over the five-day post-challenge time-course. However, no upregulation could be recorded in spleen tissue measured by qPCR in juvenile trout (270 g bw). Using mRNA sequencing again, mRNA abundances were determined in gill of juvenile trout (~57.7 g bw) exposed to various aquaculture stressors. The results indicated that all six nkls (nkl1-4 and nkl-like a and nkl-like b) were downregulated when exposed to high temperature, and that nkl1 was significantly downregulated following salinity challenge. Overall, these newly characterized AMPs may contribute to host innate immunity as they are modulated following pathogen challenge and by physiological stressors.

Identification of High-Confidence Structural Variants in Domesticated Rainbow Trout Using Whole-Genome Sequencing

Genomic structural variants (SVs) are a major source of genetic and phenotypic variation but have not been investigated systematically in rainbow trout (Oncorhynchus mykiss), an important aquaculture species of cold freshwater. The objectives of this study were 1) to identify and validate high-confidence SVs in rainbow trout using whole-genome re-sequencing; and 2) to examine the contribution of transposable elements (TEs) to SVs in rainbow trout. A total of 96 rainbow trout, including 11 homozygous lines and 85 outbred fish from three breeding populations, were whole-genome sequenced with an average genome coverage of 17.2×. Putative SVs were identified using the program Smoove which integrates LUMPY and other associated tools into one package. After rigorous filtering, 13,863 high-confidence SVs were identified. Pacific Biosciences long-reads of Arlee, one of the homozygous lines used for SV detection, validated 98% (3,948 of 4,030) of the high-confidence SVs identified in the Arlee homozygous line. Based on principal component analysis, the 85 outbred fish clustered into three groups consistent with their populations of origin, further indicating that the high-confidence SVs identified in this study are robust. The repetitive DNA content of the high-confidence SV sequences was 86.5%, which is much higher than the 57.1% repetitive DNA content of the reference genome, and is also higher than the repetitive DNA content of Atlantic salmon SVs reported previously. TEs thus contribute substantially to SVs in rainbow trout as TEs make up the majority of repetitive sequences. Hundreds of the high-confidence SVs were annotated as exon-loss or gene-fusion variants, and may have phenotypic effects. The high-confidence SVs reported in this study provide a foundation for further rainbow trout SV studies.

Comparative Structural and Antigenic Characterization of Genetically Distinct Flavobacterium psychrophilum O-Polysaccharides

Little is known about the underlying basis of serotype specificity among strains of Flavobacterium psychrophilum, the agent of rainbow trout fry syndrome and bacterial cold-water disease. The identification of different heat-stable O-serotypes among strains of this gram-negative pathogen does, however, suggest structural variations in the O-polysaccharide (O-PS) moiety of cell surface lipopolysaccharide (LPS). A trisaccharide composed of L-rhamnose (L-Rha), 2-acetamido-2-deoxy-L-fucose (L-FucNAc) and 2-acetamido-4-R-2,4-dideoxy-D-quinovose (D-Qui2NAc4NR), where R represents a dihydroxyhexanamido derivative, was previously identified as the repeating unit of Fp CSF259-93 O-PS. Interestingly, the O-PS gene cluster of this strain and that of Fp 950106-1/1, which belongs to a different O-serotype, are identical except for wzy, which encodes the putative polymerase that links trisaccharide repeats into O-PS chains. We have now found from results of glycosyl composition analysis and high-resolution nuclear magnetic resonance, that the linkage of D-Qui2NAc4NR to L-Rha, which is α1-2 for Fp CSF259-93 versus β1-3 for Fp 950106-1/1, is the only structural difference between O-PS from these strains. The corresponding difference in O-serotype specificity was established from the reactions of rabbit and trout anti-F. psychrophilum antibody with purified O-PS and LPS. Moreover, LPS-based differences in antigenicity were noted between strains with O-PS loci identical to those of Fp CSF259-93 or Fp 950106-1/1, except for the genes predicted to direct synthesis of different R-groups in Qui2NAc4NR. The findings provide a framework for defining the genetic basis of O-PS structure and antigenicity and suggest that the repertoire of F. psychrophilum O-serotypes extends beyond what is presently recognized from serological studies of this important fish pathogen.

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.

Analysis of the gut and gill microbiome of resistant and susceptible lines of rainbow trout (Oncorhynchus mykiss)

Commensal microorganisms present at mucosal surfaces play a vital role in protecting the host organism from bacterial infection. There are multiple factors that contribute to selecting for the microbiome, including host genetics. Flavobacterium psychrophilum, the causative agent of Bacterial Cold Water Disease in salmonids, accounts for acute losses in wild and farmed rainbow trout (Oncorhynchus mykiss). The U.S. National Center for Cool and Cold Water Aquaculture has used family-based selective breeding to generate a line of rainbow trout with enhanced resistance to F. psychrophilum. The goal of this study is to determine whether selective breeding impacts the gut and gill microbiome of the F. psychrophilum-resistant as compared to a background matched susceptible trout line. Mid-gut and gill samples were collected from juvenile fish maintained at high or low stocking densities and microbial diversity assessed by 16S rDNA amplicon sequencing. Results indicate that alpha diversity was significantly higher in the mid-gut of the susceptible line compared to the resistant line, while no significant differences in alpha diversity were observed in the gills. Mycoplasma sp. was the dominant taxon in the mid-gut of both groups, although it was present at a decreased abundance in the susceptible line. We also observed an increased abundance of the potential opportunistic pathogen Brevinema andersonii in the susceptible line. Within the gills, both lines exhibited similar microbial profiles, with Candidatus Branchiomonas being the dominant taxon. Together, these results suggest that selectively bred F. psychrophilum-resistant trout may harness a more resilient gut microbiome, attributing to the disease resistant phenotype. Importantly, interactions between host genetics and environmental factors such as stocking density have a significant impact in shaping trout microbial communities.

Differential expression and evolution of three tandem, interleukin-1 receptor-like 1 genes in rainbow trout (Oncorhynchus mykiss)

Interleukin-1 receptor-like 1 (Il1rl1 or ST2), a member of the interleukin-1 receptor family, has pleiotropic roles including tissue homeostasis, inflammation, immune polarization, and disease resistance in mammals. A single orthologue was previously described in salmonid fish; however, a recently improved genome assembly of rainbow trout (Oncorhynchus mykiss) revealed three adjacent, tandem il1rl1 orthologues on chromosome Omy 03. Here, we report the genomic organization and evolution of the three il1rl1 genes (il1rl1α, il1rl1β, il1rl1γ), and use both RNA-seq and gene-specific qPCR methods to quantify expression patterns. Nucleotide sequence homology between the three genes is >95% and each predicted protein contains three IG/IG-like domains, a transmembrane region and a TIR domain. The amino acid sequence homology of the rainbow trout il1rl1 genes are highly related to two functional copies in Atlantic and Coho salmon (∼94%) but relatively low (22-26%) with avian and mammalian species. Transcript abundance measured by RNA-seq in 15 tissues of healthy adult rainbow trout indicate constitutive expression of each gene. In whole body lysates, il1rl1α was shown to have >20 fold mRNA expression compared to il1rl1β and il1rl1γ as measured by qPCR assays specific to il1rl1α or il1rl1γ, as well as a multi-gene qPCR assay (il1rl1α,β,γ). Unrooted phylogenetic trees grouped the rainbow trout il1rl1 genes apart from other interleukin-1 receptor family genes and genomic comparisons identify preserved synteny between mammals, birds and salmonids albeit a pseudogene is present in both Atlantic salmon and Coho salmon. Phylogenetic analyses suggest that the three genes arose by tandem duplication but are inconclusive whether these events occurred prior-to or after salmonid speciation. These findings further the understanding of interleukin receptor family evolution and their contribution to teleost immune function.

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.

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.

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.

A BCWD-resistant line of rainbow trout exhibits higher abundance of IgT+ B cells and heavy chain tau transcripts compared to a susceptible line following challenge with Flavobacterium psychrophilum

Bacterial Cold Water Disease (BCWD) is a common, chronic disease in rainbow trout, and is caused by the gram-negative bacterium Flavobacterium psychrophilum (Fp). Through selective breeding, the National Center for Cool and Cold Water Aquaculture has generated a genetic line that is highly resistant to Fp challenge, designated ARS-Fp-R (or R-line), as well as a susceptible "control" line, ARS-Fp-S (S-line). In previous studies, resistance to Fp had been shown to correlate with naive animal spleen size, and further, naïve R-line trout had been shown to have a lower abundance of IgM+ and IgM++ cells compared to S-line fish. Here we wished to first determine whether the abundance of IgT+ and/or IgT++ cells differed between the two lines in naïve fish, and if so, how these patterns differed after in vivo challenge with Fp. Fp challenge was by intramuscular injection of live Fp and tissue collections were on days 5, 6, and/or 28 post-challenge, in two independent challenge experiments. Flow cytometric and gene expression analyses revealed that naïve R-line fish had a higher abundance of IgT+ B cells in their anterior kidney, spleen, and blood, compared to S line fish. Further, that after Fp challenge, this difference was maintained between the two lines. Lastly, abundance of IgT+ B cells and expression of secHCtau correlated with lower Fp pathogen loads in challenged fish. In the anterior kidney, IgM+ B cell abundance correlated with increased Fp loads. Together, these results suggest that IgT+ B lineage cells may have a protective function in the immune response to Fp.

Evolution of interleukin-1 receptor-like 1 and its role in rainbow trout (Oncorhynchus mykiss) resistance to Flavobacterium psychrophilum

Rainbow trout exhibit extensive phenotypic variation in innate disease resistance and we have divergently selected lines with either increased or reduced survival following exposure to the gram-negative bacterium, Flavobacterium psychrophilum (Fp). Following five generations of selection, gene expression was analyzed between the resistant and susceptible lines by RNA-seq and RT-PCR in either whole body or spleen tissue samples. Lower base-line expression of interleukin-1 receptor-like 1 (il1rl1) was identified in susceptible line fish as well as decreased expression in response to infection. Genomic analyses identified three tandem, il1rl1 genes on chromosome 3 located within a previously identified, bacterial cold water disease QTL. Here, we report the genomic organization and evolution of the three il1rl1 genes, and develop gene specific assays to measure the expression of the individual il1rl1 genes in susceptible and resistant rainbow trout lines. Sequence homology between these three putative genes is approximately 95%. Each gene contains a signal peptide, three IG/IG-like domains, a transmembrane region and a TIR domain. Two il1rl1 genes and one pseudogene are present in Atlantic salmon (Salmo salar) genome with 90–94% sequence identity. Phylogenetic analyses suggest that the salmonid il1rl1 genes expanded by tandem duplication within the salmonid linage. Our working model is that il1rl1 gene copy number or differential expression may contribute to specific disease resistance.

Dataset of proinflammatory cytokine and cytokine receptor gene expression in rainbow trout (Oncorhynchus mykiss) measured using a novel GeXP multiplex, RT-PCR assay

A GeXP multiplex, RT-PCR assay was developed and optimized that simultaneously measures expression of a suite of immune-relevant genes in rainbow trout (Oncorhynchus mykiss), concentrating on tumor necrosis factor and interleukin-1 ligand/receptor systems and acute phase response genes. The dataset includes expression values for drpt, il11a, il1b1, il1b2, il1b3, il1r-like-1(e3-5), il1r-like-1(e9-11), il1r1-like-a, il1r1-like-b, il1r2, saa, tnfa1, tnfa2, tnfa3, tnfrsf1a, tnfrsf1a-like-a, tnfrsf1a-like-b, tnfrsf5, and tnfrsf9. Gene expression was measured at four time-points post-challenge in both a resistant line (ARS-Fp-R) and a susceptible line (ARS-Fp-S) of rainbow trout. In addition, fish body weight, spleen index and the Flavobacterium psychrophilum load are reported. These data are an extension of information presented and discussed in “Proinflammatory cytokine and cytokine receptor gene expression kinetics following challenge with Flavobacterium psychrophilum in resistant and susceptible lines of rainbow trout (Oncorhynchus mykiss)” (Kutyrev et al., 2016) [1].

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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-017-0293-6) contains supplementary material, which is available to authorized users.

Proinflammatory cytokine and cytokine receptor gene expression kinetics following challenge with Flavobacterium psychrophilum in resistant and susceptible lines of rainbow trout (Oncorhynchus mykiss)

Flavobacterium psychrophilum (Fp) is the causative agent of bacterial cold water disease (BCWD) which causes appreciable economic losses in rainbow trout aquaculture. We previously reported development of a genetic line, designated ARS-Fp-R that exhibits higher survival relative to a susceptible line, designated ARS-Fp-S, following either laboratory or natural on-farm challenge. The objectives of this study were to determine the temporal kinetics of gene expression between experimentally-challenged ARS-Fp-R and ARS-Fp-S fish and the correlation between gene expression and pathogen load. We developed a GeXP multiplex RT-PCR assay to simultaneously examine expression of immune-relevant genes, concentrating on tumor necrosis factor and interleukin-1 ligand/receptor systems and acute phase response genes. Spleen tissue was sampled at 6 h, 24 h, 48 h and 144 h post-challenge and pathogen load quantified by qPCR. Transcript abundance of cytokine genes tnfa1, tnfa2, tnfa3, il1b1, il1b2, il11a; acute phase response genes saa and drtp1; and putative cytokine receptors il1r1-like-b, il1r2, tnfrsf1a, tnfrsf9, tnfrsf1a-like-b increased following challenge while the transcript abundance of il1r-like-1 and tnfrsf1a-like-a decreased compared to PBS-injected line-matched control fish. Principal component analysis identified transcript levels of genes il1r-like-1 and tnfrsf1a-like-a as exhibiting differential expression between genetic lines. In summary, Fp i.p. injection challenge elicited a proinflammatory cytokine gene expression response in the spleen, with ARS-Fp-R line fish exhibiting modestly higher basal expression levels of several putative cytokine receptors. This study furthers the understanding of the immune response following Fp challenge and differences in gene expression associated with selective breeding for disease resistance.

Differential expression of long non-coding RNAs in three genetic lines of rainbow trout in response to infection with Flavobacterium psychrophilum

Bacterial cold-water disease caused by Flavobacterium psychrophilum is one of the major causes of mortality of salmonids. Three genetic lines of rainbow trout designated as ARS-Fp-R (resistant), ARS-Fp-C (control) and ARS-Fp-S (susceptible) have significant differences in survival rate following F. psychrophilum infection. Previous study identified transcriptome differences of immune-relevant protein-coding genes at basal and post infection levels among these genetic lines. Using RNA-Seq approach, we quantified differentially expressed (DE) long non-coding RNAs (lncRNAs) in response to F. psychrophilum challenge in these genetic lines. Pairwise comparison between genetic lines and different infection statuses identified 556 DE lncRNAs. A positive correlation existed between the number of the differentially regulated lncRNAs and that of the protein-coding genes. Several lncRNAs showed strong positive and negative expression correlation with their overlapped, neighboring and distant immune related protein-coding genes including complement components, cytokines, chemokines and several signaling molecules involved in immunity. The correlated expressions and genome-wide co-localization suggested that some lncRNAs may be involved in regulating immune-relevant protein-coding genes. This study provides the first evidence of lncRNA-mediated regulation of the anti-bacterial immune response in a commercially important aquaculture species and will likely help developing new genetic markers for rainbow trout disease resistance.

Evaluation of Genome-Enabled Selection for Bacterial Cold Water Disease Resistance Using Progeny Performance Data in Rainbow Trout: Insights on Genotyping Methods and Genomic Prediction Models

Bacterial cold water disease (BCWD) causes significant economic losses in salmonid aquaculture, and traditional family-based breeding programs aimed at improving BCWD resistance have been limited to exploiting only between-family variation. We used genomic selection (GS) models to predict genomic breeding values (GEBVs) for BCWD resistance in 10 families from the first generation of the NCCCWA BCWD resistance breeding line, compared the predictive ability (PA) of GEBVs to pedigree-based estimated breeding values (EBVs), and compared the impact of two SNP genotyping methods on the accuracy of GEBV predictions. The BCWD phenotypes survival days (DAYS) and survival status (STATUS) had been recorded in training fish (n = 583) subjected to experimental BCWD challenge. Training fish, and their full sibs without phenotypic data that were used as parents of the subsequent generation, were genotyped using two methods: restriction-site associated DNA (RAD) sequencing and the Rainbow Trout Axiom® 57 K SNP array (Chip). Animal-specific GEBVs were estimated using four GS models: BayesB, BayesC, single-step GBLUP (ssGBLUP), and weighted ssGBLUP (wssGBLUP). Family-specific EBVs were estimated using pedigree and phenotype data in the training fish only. The PA of EBVs and GEBVs was assessed by correlating mean progeny phenotype (MPP) with mid-parent EBV (family-specific) or GEBV (animal-specific). The best GEBV predictions were similar to EBV with PA values of 0.49 and 0.46 vs. 0.50 and 0.41 for DAYS and STATUS, respectively. Among the GEBV prediction methods, ssGBLUP consistently had the highest PA. The RAD genotyping platform had GEBVs with similar PA to those of GEBVs from the Chip platform. The PA of ssGBLUP and wssGBLUP methods was higher with the Chip, but for BayesB and BayesC methods it was higher with the RAD platform. The overall GEBV accuracy in this study was low to moderate, likely due to the small training sample used. This study explored the potential of GS for improving resistance to BCWD in rainbow trout using, for the first time, progeny testing data to assess the accuracy of GEBVs, and it provides the basis for further investigation on the implementation of GS in commercial rainbow trout populations.

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.

Identification of single nucleotide polymorphism markers associated with bacterial cold water disease resistance and spleen size in rainbow trout

Bacterial cold water disease (BCWD) is one of the frequent causes of elevated mortality in salmonid aquaculture. Previously, we identified and validated microsatellites on chromosome Omy19 associated with QTL (quantitative trait loci) for BCWD resistance and spleen size in rainbow trout. Recently, SNPs (single nucleotide polymorphism) have become the markers of choice for genetic analyses in rainbow trout as they are highly abundant, cost-effective and are amenable for high throughput genotyping. The objective of this study was to identify SNP markers associated with BCWD resistance and spleen size using both genome-wide association studies (GWAS) and linkage-based QTL mapping approaches. A total of 298 offspring from the two half-sib families used in our previous study to validate the significant BCWD QTL on chromosome Omy19 were genotyped with RAD-seq (restriction-site-associated DNA sequencing), and 7,849 informative SNPs were identified. Based on GWAS, 18 SNPs associated with BCWD resistance and 20 SNPs associated with spleen size were identified. Linkage-based QTL mapping revealed three significant QTL for BCWD resistance. In addition to the previously validated dam-derived QTL on chromosome Omy19, two significant BCWD QTL derived from the sires were identified on chromosomes Omy8 and Omy25, respectively. A sire-derived significant QTL for spleen size on chromosome Omy2 was detected. The SNP markers reported in this study will facilitate fine mapping to identify positional candidate genes for BCWD resistance in rainbow trout.

Potential mechanisms of attenuation for rifampicin-passaged strains of Flavobacterium psychrophilum

Background

Flavobacterium psychrophilum is the etiologic agent of bacterial coldwater disease in salmonids. Earlier research showed that a rifampicin-passaged strain of F. psychrophilum (CSF 259-93B.17) caused no disease in rainbow trout (Oncorhynchus mykiss, Walbaum) while inducing a protective immune response against challenge with the virulent CSF 259–93 strain. We hypothesized that rifampicin passage leads to an accumulation of genomic mutations that, by chance, reduce virulence. To assess the pattern of phenotypic and genotypic changes associated with passage, we examined proteomic, LPS and single-nucleotide polymorphism (SNP) differences for two F. psychrophilum strains (CSF 259–93 and THC 02–90) that were passaged with and without rifampicin selection.

Results

Rifampicin resistance was conveyed by expected mutations in rpoB, although affecting different DNA bases depending on the strain. One rifampicin-passaged CSF 259–93 strain (CR) was attenuated (4 % mortality) in challenged fish, but only accumulated eight nonsynonymous SNPs compared to the parent strain. A CSF 259–93 strain passaged without rifampicin (CN) accumulated five nonsynonymous SNPs and was partially attenuated (28 % mortality) compared to the parent strain (54.5 % mortality). In contrast, there were no significant change in fish mortalities among THC 02–90 wild-type and passaged strains, despite numerous SNPs accumulated during passage with (n = 174) and without rifampicin (n = 126). While only three missense SNPs were associated with attenuation, a Ser492Phe rpoB mutation in the CR strain may contribute to further attenuation. All strains except CR retained a gliding motility phenotype. Few proteomic differences were observed by 2D SDS-PAGE and there were no apparent changes in LPS between strains. Comparative methylome analysis of two strains (CR and TR) identified no shared methylation motifs for these two strains.

Conclusion

Multiple genomic changes arose during passage experiments with rifampicin selection pressure. Consistent with our hypothesis, unique strain-specific mutations were detected for the fully attenuated (CR), partially attenuated (CN) and another fully attenuated strain (B17).

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-015-0518-1) contains supplementary material, which is available to authorized users.

Reduction of rainbow trout spleen size by splenectomy does not alter resistance against bacterial cold water disease

In lower vertebrates, the contribution of the spleen to anti-bacterial immunity is poorly understood. We have previously reported a phenotypic and genetic correlation between resistance to Flavobacterium psychrophilum, the causative agent of bacterial cold water disease (BCWD) and spleen somatic index (spleen weight normalized to body weight, SI). Fish families with larger pre-challenge SI values were found to have greater BCWD survival (resistance) following intraperitoneal injection of a lethal dose of F. psychrophilum. Since the mammalian spleen is known to be crucial for capture and destruction of encapsulated bacteria, we tested the hypothesis that reduction of spleen size, by surgical splenectomy, should reduce the survival advantage of the larger-spleen, disease-resistant fish. Experiments were performed using two separate lines of fish that had previously been selected either based on BCWD survival (resistant and susceptible), or selected based on spleen size (high and low SI). Following 65 to 81 days post-surgical recovery, fish were challenged with F. psychrophilum and mortality monitored for a minimum of 21 days. No significant difference in the relative survival was detected between splenectomized or sham-operated groups, while SI of splenectomized fish was reduced to an average of 8-12% of control animals. A positive correlation was observed between the SI, measured at the time of splenectomy, and time-to-death post-challenge. In summary, these experiments argue that larger spleen size alone is not sufficient for greater BCWD resistance, but rather it is an indirect indicator of immunological status.

Whole-body transcriptome of selectively bred, resistant-, control-, and susceptible-line rainbow trout following experimental challenge with Flavobacterium psychrophilum

Genetic improvement for enhanced disease resistance in fish is an increasingly utilized approach to mitigate endemic infectious disease in aquaculture. In domesticated salmonid populations, large phenotypic variation in disease resistance has been identified but the genetic basis for altered responsiveness remains unclear. We previously reported three generations of selection and phenotypic validation of a bacterial cold water disease (BCWD) resistant line of rainbow trout, designated ARS-Fp-R. This line has higher survival after infection by either standardized laboratory challenge or natural challenge as compared to two reference lines, designated ARS-Fp-C (control) and ARS-Fp-S (susceptible). In this study, we utilized 1.1 g fry from the three genetic lines and performed RNA-seq to measure transcript abundance from the whole body of naive and Flavobacterium psychrophilum infected fish at day 1 (early time-point) and at day 5 post-challenge (onset of mortality). Sequences from 24 libraries were mapped onto the rainbow trout genome reference transcriptome of 46,585 predicted protein coding mRNAs that included 2633 putative immune-relevant gene transcripts. A total of 1884 genes (4.0% genome) exhibited differential transcript abundance between infected and mock-challenged fish (FDR < 0.05) that included chemokines, complement components, tnf receptor superfamily members, interleukins, nod-like receptor family members, and genes involved in metabolism and wound healing. The largest number of differentially expressed genes occurred on day 5 post-infection between naive and challenged ARS-Fp-S line fish correlating with high bacterial load. After excluding the effect of infection, we identified 21 differentially expressed genes between the three genetic lines. In summary, these data indicate global transcriptome differences between genetic lines of naive animals as well as differentially regulated transcriptional responses to infection.

B cell signatures of BCWD-resistant and susceptible lines of rainbow trout: a shift towards more EBF-expressing progenitors and fewer mature B cells in resistant animals

Bacterial cold water disease (BCWD) is a chronic disease of rainbow trout, and is caused by the Gram-negative bacterium Flavobacterium psychrophilum (Fp), a common aquaculture pathogen. The National Center for Cool and Cold Water Aquaculture has bred two genetic lines of rainbow trout: a line of Fp-resistant trout (ARS-Fp-R or R-line trout) and a line of susceptible trout (ARS-Fp-S, or S-line). Little is known about how phenotypic selection alters immune response parameters or how such changes relate to genetic disease resistance. Herein, we quantify interindividual variation in the distribution and abundance of B cell populations (B cell signatures) and examine differences between genetic lines of naive animals. There are limited trout-specific cell surface markers currently available to resolve B cell subpopulations and thus we developed an alternative approach based on detection of differentially expressed transcription factors and intracellular cytokines. B cell signatures were compared between R-line and S-line trout by flow cytometry using antibodies against transcription factors early B cell factor-1 (EBF1) and paired domain box protein Pax5, the pro-inflammatory cytokine IL-1β, and the immunoglobulin heavy chain mu. R-line trout had higher percentages of EBF(+) B myeloid/ progenitor and pre-B cells in PBL, anterior and posterior kidney tissues compared to S-line trout. The opposite pattern was detected in more mature B cell populations: R-line trout had lower percentages of both IgM(+) mature B cells and IgM-secreting cells in anterior kidney and PBL compared to S-line trout. In vitro LPS-activation studies of PBL and spleen cell cultures revealed no significant induction differences between R-line and S-line trout. Together, our findings suggest that selective resistance to BCWD may be associated with shifts in naive animal developmental lineage commitment that result in decreased B lymphopoiesis and increased myelopoiesis in BCWD resistant trout relative to susceptible trout.

Biochemical reference intervals and pathophysiological changes in Flavobacterium psychrophilum-resistant and -susceptible rainbow trout lines

Host genetic resistance against disease-causing pathogens can be enhanced through family-based selective breeding. At present, there is an incomplete understanding of how artificial selection of fish alters host physiology and response following pathogen exposure. We previously reported the generation of selectively-bred rainbow trout Oncorhynchus mykiss lines with either increased resistance (ARS-Fp-R) or susceptibility (ARS-Fp-S) to bacterial cold water disease (BCWD). This study (1) determined baseline reference-range intervals for packed cell volume (PCV) and 18 plasma biochemistry analytes, and (2) examined pathophysiological changes following infection between the genetic lines. PCV and biochemistry reference-range intervals did not significantly differ between genetic lines; thus data were pooled into a single reference-range population (n = 85). ARS-Fp-R and ARS-Fp-S line fish were intraperitoneally challenged with Flavobacterium psychrophilum, and plasma was collected on Days 1, 3, 6, and 9 post-challenge. Splenic bacterial load was measured using an F. psychrophilum-specific qPCR assay. In both genetic lines, changes were observed in mean PCV, total protein, albumin, glucose, cholesterol, chloride, and calcium, falling outside the established reference intervals and significantly differing from phosphate-buffered saline challenged fish, on at least 1d post-challenge. Mean PCV, total protein, and calcium significantly differed between ARS-Fp-R and ARS-Fp-S line fish on Day 9 post-infection, with values in the ARS-Fp-S line deviating most from the reference interval. PCV, total protein, cholesterol, and calcium negatively correlated with bacterial load. These findings identify divergent pathophysiological responses between ARS-Fp-R and ARS-Fp-S line fish following laboratory challenge that are likely associated with differential survival.

Histopathologic Changes in Disease-Resistant-Line and Disease-Susceptible-Line Juvenile Rainbow Trout Experimentally Infected with Flavobacterium psychrophilum

Abstract A challenge to selectively breeding fish populations for improved disease resistance in aquaculture is an incomplete understanding of how artificial selection affects innate immunity at the host-pathogen level. The objective of this study was to determine whether Rainbow Trout Oncorhynchus mykiss bred for differential susceptibility to bacterial cold-water disease exhibited altered tissue damage and cellular inflammatory response following experimental challenge with Flavobacterium psychrophilum. Fish from disease-resistant (ARS-Fp-R) and disease-susceptible (ARS-Fp-S) lines were experimentally challenged as juveniles, and mortalities, as well as survivors, were sampled for histopathology during the acute phase of the disease. Microscopic lesions were quantified or semiquantified and statistically compared for changes over time and between genetic lines. Significant progression in the degree of perisplenitis, splenic necrosis, splenic inflammatory infiltrates, average splenic ellipsoid area, total splenic ellipsoid area, and peritonitis was present over time in both genetic lines on at least one postinfection time point. No differences were found between renal inflammatory infiltrates and renal hematopoietic cell depletion over time. Perisplenitis was significantly lower in fish from the ARS-Fp-R line on day 9 postinfection than in fish from the ARS-Fp-S line. The ARS-Fp-R line demonstrated a trend towards reduced splenic necrosis compared with the ARS-Fp-S line that approached significance, and fish from the ARS-Fp-S line were 3.6 times more likely than fish from the ARS-Fp-R line to have a higher splenic necrosis lesion score after day 3 postinfection. These findings support the hypothesis that differential survival is a result of divergence in disease magnitude and not altered disease course between genetic lines. Characterization of histopathologic changes between genetic lines and over time helps elucidate mechanisms of disease resistance and contributes to our understanding of disease pathogenesis in fish infected with F. psychrophilum. Received January 7, 2014; accepted March 10, 2014.

Detection of QTL in rainbow trout affecting survival when challenged with Flavobacterium psychrophilum

Bacterial cold water disease (BCWD) causes significant economic loss in salmonid aquaculture. We previously detected genetic variation in survival following challenge with Flavobacterium psychrophilum (Fp), the causative agent of BCWD in rainbow trout (Oncorhynchus mykiss). A family-based selection program to improve resistance was initiated in 2005 at the USDA National Center for Cool and Cold Water Aquaculture. Select crosses were made in 2007 and 2009 to evaluate family-based disease survival using Fp injection challenges. From each putative F₂/BC₁ family generated in 2009, 200-260 fish were challenged in 4-7 replicates per family. Whole genome QTL scans of three F₂/BC₁ families were conducted with about 270 informative microsatellite loci per family spaced at an average interval size of 6 cM throughout the rainbow trout genome. Markers on chromosomes containing QTL were further evaluated in three additional F₂/BC₁ families. The additional F₂/BC₁ families were sire or dam half-sibs (HS) of the initially genome scanned families. Overall, we identified nine major QTL on seven chromosomes that were significant or highly significant with moderate to large effects of at least 13 % of the total phenotypic variance. The largest effect QTL for BCWD resistance explaining up to 40 % of the phenotypic variance was detected on chromosome OMY8 in family 2009070 and in the combined dam HS family 2009069-070. The nine major QTL identified in this study are candidates for fine mapping to identify new markers that are tightly linked to disease resistance loci for using in marker assisted selection strategies.

Transfer of serum and cells from Yersinia ruckeri vaccinated doubled-haploid hot creek rainbow trout into outcross F1 progeny elucidates mechanisms of vaccine-induced protection

Yersinia ruckeri is a well-established bacterial pathogen for many salmonid species, against which a formalin-killed bacterin vaccine has been effective in reducing disease outbreaks. Previous studies have reported conflicting results about the protective value of the systemic humoral response to Y. ruckeri vaccination. Here we directly demonstrate that plasma contains the long-term protective component elicited by both immersion and intraperitoneal injection vaccination of rainbow trout. A total of 0.5 μL of plasma from vaccinated fish provided almost complete protection against experimental challenge. Conversely, the cells obtained from peripheral blood conferred little or no protection in naïve recipients. The protective component of immune sera was IgM based on size exclusion chromatography and recognition by monoclonal antibody Warr 1-14. Immune plasma generated against a Y. ruckeri biotype 1 strain protected equally against challenges with Y. ruckeri biotype 1 and 2 strains. These results illustrate the importance of the humoral IgM response against Y. ruckeri and the use of doubled haploid rainbow trout (Oncorhynchus mykiss) and transfer of plasma/serum and cells into F1 outcross progeny as a model system for dissection of the mechanism(s) of vaccine-induced protection.

Acute mortality, bacterial load, and pathology of select lines of adult rainbow trout challenged with Weissella sp. NC36

A challenge for improving disease resistance in fish through genetics is to understand specificity of resistance and whether selection for one pathogen alters the response to unrelated pathogenic microorganisms. Adult Rainbow Trout Oncorhynchus mykiss that had been bred for differential susceptibility to Flavobacterium psychrophilum, the causative agent of bacterial cold water disease (BCWD) and designated ARS-Fp-R (resistant), ARS-Fp-S (susceptible), and ARS-Fp-C (control line), as well as a pool of commercial-stock Rainbow Trout, were intraperitoneally challenged with Weissella sp. NC36. Clinical signs, survival, and innate mechanisms affecting disease resistance were monitored over 9 d. Acute disease signs included exophthalmia associated with retrobulbar inflammation and hemorrhage, cerebral hemorrhage, and mild to moderate granulomatous pericarditis. The ARS-Fp-R line did not demonstrate significant survival differences over a 9-d period compared with the ARS-Fp-C and ARS-Fp-S lines (P ≥ 0.09) indicating that during the acute phase of disease, the resistance factors that limit BCWD do not confer cross protection against Weissella sp. NC36. The linear effect of body weight at challenge was statistically significant, as each 10-g increase in body weight increased the hazard of death by 1% (P = 0.02). Bacterial loads on day 3, assessed by splenic and cerebral CFU counts, did not differ between ARS-Fp-R and ARS-Fp-S trout and there was no correlation between CFU counts and body weight. These findings help elucidate the specificity of disease resistance in selectively bred lines and contribute to our understanding of disease caused by Weissella sp., a recently described pathogen found in cultured Rainbow Trout.

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 (h(2)  = 0.56±0.18) and odd-year (h(2)  = 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.

A real-time polymerase chain reaction assay for identification and quantification of Flavobacterium psychrophilum and application to disease resistance studies in selectively bred rainbow trout Oncorhynchus mykiss

Rapid detection and quantification of Flavobacterium psychrophilum, the causative agent of bacterial cold water disease (BCWD) in rainbow trout, are crucial to disease surveillance and encompass an essential component of BCWD research. Real-time, or quantitative polymerase chain reaction (qPCR) assays that have previously targeted the 16S rRNA gene of F. psychrophilum are complicated by polymorphisms and off-target amplification. Insignia primer and probe development software were used to identify a conserved single-copy signature sequence in the F. psychrophilum genome that codes for a hypothetical protein with unknown function. Primer and probes were used in a TaqMan qPCR assay that amplified 210 F. psychrophilum isolates with a lower limit of linear detection at 3.1 genome equivalents per reaction, with no amplification of 23 nontarget bacterial isolates. The assay was not inhibited by host spleen DNA or spleen homogenate. Methods were successfully applied to detect F. psychrophilum in rainbow trout from naturally occurring BCWD outbreaks and to quantify bacterial loads in experimentally infected rainbow trout. This assay will be applied to future studies to characterize disease pathogenesis in fish selectively bred for BCWD resistance.

Origin and evolution of TNF and TNF receptor superfamilies

The tumor necrosis factor superfamily (TNFSF) and the TNF receptor superfamily (TNFRSF) have an ancient evolutionary origin that can be traced back to single copy genes within Arthropods. In humans, 18 TNFSF and 29 TNFRSF genes have been identified. Evolutionary models account for the increase in gene number primarily through multiple whole genome duplication events as well as by lineage and/or species-specific tandem duplication and translocation. The identification and functional analyses of teleost ligands and receptors provide insight into the critical transition between invertebrates and higher vertebrates. Bioinformatic analyses of fish genomes and EST datasets identify 14 distinct ligand groups, some of which are novel to teleosts, while to date, only limited numbers of receptors have been characterized in fish. The most studied ligand is TNF of which teleost species possess between 1 and 3 copies as well as a receptor similar to TNFR1. Functional studies using zebrafish indicate a conserved role of this ligand-receptor system in the regulation of cell survival and resistance to infectious disease. The increasing interest and use of TNFSF and TNFRSF modulators in human and animal medicine underscores the need to understand the evolutionary origins as well as conserved and novel functions of these biologically important molecules.

Identification of immunogenic proteins within distinct molecular mass fractions of Flavobacterium psychrophilum

Flavobacterium psychrophilum is the aetiological agent of bacterial coldwater disease (CWD), and this pathogen has large economic impacts on salmonid aquaculture worldwide. Previously, it was demonstrated that high levels of protection against F. psychrophilum challenge were conferred to rainbow trout, Oncorhynchus mykiss (Walbaum), by immunization with distinct molecular mass fractions of the bacterium, and specific antibodies were correlated with protection. In this study, an immunoproteomic analysis of F. psychrophilum was performed using two-dimensional polyacrylamide gel electrophoresis and Western blotting with serum from fish immunized with high- and mid-molecular mass fractions of the bacterium. Mass spectrometry was used to determine the protein identity, and 15 immunogenic proteins were positively identified following Mascot searches of the F. psychrophilum genome. Based on known function and immunogenicity of homologous proteins in other bacterial pathogens, antibodies specific for several of the identified proteins may be important for protective immunity from CWD. These include outer membrane protein OmpA (P60), trigger factor, ClpB, elongation factor G, gliding motility protein GldN and a conserved hypothetical protein. This work increases the understanding of the protective humoral immune response of rainbow trout against these distinct molecular mass fractions of F. psychrophilum and provides new potential targets for recombinant protein vaccine development.

Temporal and pathogen-load dependent changes in rainbow trout (Oncorhynchus mykiss) immune response traits following challenge with biotype 2 Yersinia ruckeri

Rainbow trout infected with Yersinia ruckeri, the causative agent of enteric redmouth disease (ERM), produce a pro-inflammatory and acute-phase response attributed in part to the innate recognition of bacterial-produced flagellin. Recently, variants of Y. ruckeri have been identified that lack flagella and associated motility. These strains are classified as biotype 2 (BT2) Y. ruckeri and these are considered an emerging problem in salmonid aquaculture. Little is known about the salmonid immune response to these variants. Herein, we report temporal and quantitative changes in rainbow trout immune response parameters following a primary challenge with BT2 Y. ruckeri strain YRNC10. Fish were injection-challenged with ten-fold dilutions of viable bacteria and sampled on days 1, 3, 5 and 7 post-challenge. TNFalpha1 and IL1-beta1 transcripts were increased by day 1 post-challenge, and on days 3, 5 and 7 maximal gene transcript up-regulation occurred at a threshold of approximately 64-256CFU per mg spleen tissue. Infection induced robust SAA gene up-regulation that was significantly correlated with increased gene expression of IL-1beta1 (r=0.81, P<0.0001) and TNFalpha1 (r=0.55, P<0.0001). Y. ruckeri infection induced modest changes in INFgamma and Mx-1 gene transcript abundance at intermediate or high challenge doses and the expression patterns of both genes were positively correlated with pro-inflammatory gene and acute-phase gene transcription patterns. TNF superfamily 13b (BAFF) gene expression was significantly down-regulated in response to infection on days 3, 5 and 7 at the highest challenge doses. The spleen somatic index was significantly increased on days 3, 5 and 7 post-infection and positively correlated with spleen colony forming units and abundance of gene transcripts SAA, TNFalpha1, and IL1-beta1. In summary, rainbow trout had a strong innate response following challenge with BT2 Y. ruckeri strain YRNC10 indicating that flagellin expression is not required for production of a robust pro-inflammatory and acute-phase gene transcription response. This study further supports the use of SAA transcript abundance and spleen somatic index as general measures of immunological status and fish health.

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.

Identification, characterization and genetic mapping of TLR1 loci in rainbow trout (Oncorhynchus mykiss)

Induction of innate immune pathways is critical for early anti-microbial defense but there is limited understanding of how teleosts recognize microbial molecules and activate these pathways. In mammals, Toll-like receptors (TLR) 1 and 2 form a heterodimer involved in recognizing peptidoglycans and lipoproteins of microbial origin. Herein, we identify and describe the rainbow trout (Oncorhynchus mykiss) TLR1 gene ortholog and its mRNA expression. Two TLR1 loci were identified from a rainbow trout bacterial artificial chromosome (BAC) library using DNA sequencing and genetic linkage analyses. Full length cDNA clone and direct sequencing of four BACs revealed an intact omTLR1 open reading frame (ORF) located on chromosome 14 and a second locus on chromosome 25 that contains a TLR1 pseudogene. The duplicated trout loci exhibit conserved synteny with other fish genomes that extends beyond the TLR1 gene sequences. The omTLR1 gene includes a single large coding exon similar to all other described TLR1 genes, but unlike other teleosts it also has a 5' UTR exon and intron preceding the large coding exon. The omTLR1 ORF is predicted to encode an 808 amino-acid protein with 69% similarity to the Fugu TLR1 and a conserved pattern of predicted leucine-rich repeats (LRR). Phylogenetic analysis grouped omTLR1 with other fish TLR1 genes on a separate branch from the avian TLR1 and mammalian TLR1, 6 and 10. omTLR1 expression levels in rainbow trout anterior kidney leukocytes were not affected by the human TLR2/6 and TLR2/1 agonists diacylated lipoprotein (Pam(2)CSK(4)) and triacylated lipoprotein (Pam(3)CSK(4)). However, due to the lack of TLR6 and 10 genes in teleost genomes and up-regulation of TLR1 mRNA in response to LPS and bacterial infection in other fish species we hypothesize an important role for omTLR1 in anti-microbial immunity. Therefore, the identification of a TLR2 ortholog in rainbow trout and the development of assays to measure ligand binding and downstream signaling are critical for future elucidation of omTLR1 functions.

Identification, characterization and genetic mapping of TLR7, TLR8a1 and TLR8a2 genes in rainbow trout (Oncorhynchus mykiss)

Induction of the innate immune pathways is critical for early anti-viral defense but there is limited understanding of how teleost fish recognize viral molecules and activate these pathways. In mammals, Toll-like receptors (TLR) 7 and 8 bind single-stranded RNA of viral origin and are activated by synthetic anti-viral imidazoquinoline compounds. Herein, we identify and describe the rainbow trout (Oncorhynchus mykiss) TLR7 and TLR8 gene orthologs and their mRNA expression. Two TLR7/8 loci were identified from a rainbow trout bacterial artificial chromosome (BAC) library using DNA fingerprinting and genetic linkage analyses. Direct sequencing of two representative BACs revealed intact omTLR7 and omTLR8a1 open reading frames (ORFs) located on chromosome 3 and a second locus on chromosome 22 that contains an omTLR8a2 ORF and a putative TLR7 pseudogene. We used the omTLR8a1/2 nomenclature for the two trout TLR8 genes as phylogenetic analysis revealed that they and all the other teleost TLR8 genes sequenced to date are similar to the zebrafish TLR8a, but are distinct from the zebrafish TLR8b. The duplicated trout loci exhibit conserved synteny with other fish genomes extending beyond the tandem of TLR7/8 genes. The trout TLR7 and 8a1/2 genes are composed of a single large exon similar to all other described TLR7/8 genes. The omTLR7 ORF is predicted to encode a 1049 amino acid (aa) protein with 84% similarity to the Fugu TLR7 and a conserved pattern of predicted leucine-rich repeats (LRR). The omTLR8a1 and omTLR8a2 are predicted to encode 1035- and 1034-aa proteins, respectively, and have 86% similarity to each other. omTLR8a1 is likely the ortholog of the only Atlantic salmon TLR8 gene described to date as they have 95% aa sequence similarity. The tissue expression profiles of omTLR7, omTLR8a1 and omTLR8a2 in healthy trout were highest in spleen tissue followed by anterior and then posterior kidney tissues. Rainbow trout anterior kidney leukocytes produced elevated levels of pro-inflammatory and type I interferon cytokines mRNA in response to stimulation with the human TLR7/8 agonist R848 or the TLR3 agonist poly I:C. Only poly I:C-induced IFN2 transcription was significantly suppressed in the presence of chloroquine, a compound known to block endosomal acidification and inhibit endosomal maturation. The effect of chloroquine on R848-induced cytokine expression was equivocal and so it remains questionable whether rainbow trout recognition of R848 requires endosomal maturation. TLR7 and TLR8a1 expression levels in rainbow trout anterior kidney leukocytes were not affected by poly I:C or R848 treatments, but surprisingly, TLR8a2 expression was moderately down-regulated by R848. The down-regulation of omTLR8a2 may imply that this gene has evolved to a new or altered function in rainbow trout, as often occurs when the two duplicated genes remain active.

Proteomic analysis of Flavobacterium psychrophilum cultured in vivo and in iron-limited media

Flavobacterium psychrophilum is the etiologic agent of bacterial coldwater disease, but the pathogenic mechanisms of this important fish pathogen are not fully understood. Identifying bacterial genes of F. psychrophilum differentially expressed in vivo may lead to a better understanding of pathogenesis and provide targets for vaccine development. Therefore, the present study used a proteomic approach to identify and quantify proteins of F. psychrophilum following growth in vivo and under iron-limited growth conditions. As determined by 2D polyacrylamide gel electrophoresis (2D-PAGE), numerous proteins exhibited different spot intensities following culture of the bacterium in vivo, and of these, 20 were selected and identified by liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS) analysis and Mascot searches of the F. psychrophilum genome. Eighteen proteins exhibited increased spot intensities in vivo, and these included: several chaperone and stress proteins, gliding motility protein GldN, outer membrane protein OmpH, 2 probable outer membrane proteins (OmpA family), probable aminopeptidase precursor, probable lipoprotein precursor, 3-oxoacyl-[acyl-carrier-protein]-reductase, and several proteins with unknown function. Two proteins exhibited decreased spot intensities in vivo and were identified as ferritin FtnA and outer membrane protein OmpA (P60). Culture of F. psychrophilum in iron-limited media resulted in similar protein spot intensity changes for 6 of the 20 proteins identified following growth in vivo. Results from the present study suggest a role of upregulated proteins in the pathogenesis of F. psychrophilum and these may represent potential vaccine candidate antigens.

Rainbow trout resistance to bacterial cold-water disease is moderately heritable and is not adversely correlated with growth

The objectives of this study were to estimate the heritabilities for and genetic correlations among resistance to bacterial cold-water disease and growth traits in a population of rainbow trout (Oncorhynchus mykiss). Bacterial cold-water disease, a chronic disease of rainbow trout, is caused by Flavobacterium psychrophilum. This bacterium also causes acute losses in young fish, known as rainbow trout fry syndrome. Selective breeding for increased disease resistance is a promising strategy that has not been widely used in aquaculture. At the same time, improving growth performance is critical for efficient production. At the National Center for Cool and Cold Water Aquaculture, reducing the negative impact of diseases on rainbow trout culture and improving growth performance are primary objectives. In 2005, when fish averaged 2.4 g, 71 full-sib families were challenged with F. psychrophilum and evaluated for 21 d. Overall survival was 29.3% and family rates of survival varied from 1.5 to 72.5%. Heritability of postchallenge survival, an indicator of disease resistance, was estimated to be 0.35 +/- 0.09. Body weights at 9 and 12 mo posthatch and growth rate from 9 to 12 mo were evaluated on siblings of the fish in the disease challenge study. Growth traits were moderately heritable, from 0.32 for growth rate to 0.61 for 12-mo BW. Genetic and phenotypic correlations between growth traits and resistance to bacterial cold-water disease were not different from zero. These results suggest that genetic improvement can be made simultaneously for growth and bacterial cold-water disease resistance in rainbow trout by using selective breeding.

Challenges associated with heterologous expression of Flavobacterium psychrophilum proteins in Escherichia coli

A two-parameter statistical model was used to predict the solubility of 96 putative virulence-associated proteins of Flavobacterium psychrophilum (CSF259-93) upon over expression in Escherichia coli. This analysis indicated that 88.5% of the F. psychrophilum proteins would be expressed as insoluble aggregates (inclusion bodies). These solubility predictions were verified experimentally by colony filtration blot for six different F. psychrophilum proteins. A comprehensive analysis of codon usage identified over a dozen codons that are used frequently in F. psychrophilum, but that are rarely used in E. coli. Expression of F. psychrophilum proteins in E. coli was often associated with production of minor molecular weight products, presumably because of the codon usage bias between these two organisms. Expression of recombinant protein in the presence of rare tRNA genes resulted in marginal improvements in the expressed products. Consequently, Vibrio parahaemolyticus was developed as an alternative expression host because its codon usage is similar to F. psychrophilum. A full-length recombinant F. psychrophilum hemolysin was successfully expressed and purified from V. parahaemolyticus in soluble form, whereas this protein was insoluble upon expression in E. coli. We show that V. parahaemolyticus can be used as an alternate heterologous expression system that can remedy challenges associated with expression and production of F. psychrophilum recombinant proteins.

Suggestive association of major histocompatibility IB genetic markers with resistance to bacterial cold water disease in rainbow trout (Oncorhynchus mykiss)

Genes within the major histocompatibility complex (MHC) are important for both innate and adaptive immune responses in mammals; however, much less is known regarding their contribution in teleost fishes. We examined the involvement of four major histocompatibility (MH) genomic regions in rainbow trout in resistance to the causative agent of bacterial coldwater disease (BCWD), Flavobacterium psychrophilum. Fish from the 2005 NCCCWA brood-year (71 full-sib families) were challenged with F. psychrophilum strain CSF 259-93. The overall mortality rate was 70%, with large variation in mortality between families. Disease resistance was quantified as post-challenge days to death. Phenotypic variation and additive genetic variation were estimated using mixed models of survival analysis. To examine association, eight microsatellite markers were isolated from MH gene-containing BAC clones and mapped onto the rainbow trout genetic linkage map. The parents and grandparents of the 2005 brood-year families were genotyped with these eight markers and another two markers tightly linked to the MH-IB region to assess the extent of linkage disequilibrium (LD) of MH genomic regions MH-IA, MH-IB, TAP1, and MH-II with survival post-challenge. MH-IB and MH-II markers were linked to BCWD survivability when data were analyzed by family. Tests for disease association at the population level substantiated the involvement of MH-IB, but not MH-II, with disease resistance. The impact of selective breeding for disease resistance on MH sequence variation is discussed in the context of aquaculture production.

Early diversification of the TNF superfamily in teleosts: genomic characterization and expression analysis

The TNF superfamily (TNFSF) of proteins are cytokines involved in diverse immunological and developmental pathways. Little is known about their evolution or expression in lower vertebrate species. Bioinformatic searches of Zebrafish, Tetraodon, and Fugu genome and other teleost expressed sequence tag databases identified 44 novel gene sequences containing a TNF homology domain. This work reveals the following: 1) teleosts possess orthologs of BAFF, APRIL, EDA, TWEAK, 4-1BBL, Fas ligand, LIGHT, CD40L, RANKL, and possibly TL1A; 2) the BAFF-APRIL subfamily is enriched by a third member, BALM, unique to fish; 3) orthologs of lymphotoxins alpha and beta were not clearly identified in teleosts and are substituted by a related ligand, TNF-New; 4) as many as four TRAIL-like genes are present in teleosts, as compared with only one in mammals; and 5) T cell activation ligands OX40L, CD27L, CD30L, and GITRL were not identified in any fish species. Finally, we characterize mRNA expression of TNFSF members CD40L, LIGHT, BALM, APRIL, Fas ligand, RANKL, TRAIL-like, and TNF-New in rainbow trout, Oncorhynchus mykiss, immune and nonimmune tissues. In conclusion, we identified a total of 14 distinct TNFSF members in fishes, indicating expansion of this superfamily before the divergence of bony fish and tetrapods, approximately 360-450 million years ago. Based on these findings, we extend a model of TNFSF evolution and the co-emergence of the vertebrate adaptive immune system.

Identification of potential vaccine target antigens by immunoproteomic analysis of a virulent and a non-virulent strain of the fish pathogen Flavobacterium psychrophilum

Flavobacterium psychrophilum is the etiological agent of bacterial coldwater disease (CWD) and rainbow trout fry syndrome (RTFS). To identify antigens associated with virulence or host immunity, we compared total and immunogenic proteins of cellular and extracellular products (ECP) between a virulent (CSF-259-93) and non-virulent (ATCC 49418) strain of F. psychrophilum. One-dimensional sodium dodecyl sulfate polyacrylamide gel electrophoresis of total cellular proteins revealed only minor differences between the strains; however, separation of ECP showed that proteins were differentially expressed. Western blot analysis using rainbow trout (Oncorhynchus mykiss) anti-CSF-259-93 sera showed greater reactivity to proteins of the virulent strain, including many > 50 kDa. Further analysis by 2-dimensional electrophoresis (2DE) identified numerous differences between the strains. Western blot analysis combined with 2DE identified several immunogenic proteins that reacted with the antisera and were shared between the 2 strains. However, at least 15 immunogenic proteins appeared to be unique to the virulent strain, while 4 such proteins were identified in the non-virulent strain; 8 proteins unique to the virulent strain and 6 shared proteins were further analyzed for identification by liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS) analysis. Of these, 3 immunogenic proteins (heat shock proteins HSP 60 and HSP 70) and 2 other proteins (ATP synthase and thermolysin) were conclusively identified. The 2 highly immunogenic heat shock proteins were shown to share extensive homology with heat shock proteins of related bacteria. This approach for antigen identification may provide a basis for targeted vaccine development against CWD and RTFS.

Identification of novel rainbow trout (Onchorynchus mykiss) chemokines, CXCd1 and CXCd2: mRNA expression after Yersinia ruckeri vaccination and challenge

Chemokines play important roles in controlling leukocyte trafficking under normal and inflammatory conditions. Sixteen CXC chemokines have been identified in the human and mouse genomes, while considerably fewer teleost fish CXC chemokines have been reported. Here, we describe a novel clade of trout (Onchorynchus mykiss) CXC chemokines, designated Onmy CXCd, and we identify a novel gene, CXCd1, and a putative duplicate, CXCd2. The trout CXCd proteins contain 112 amino acids and the CXCd1 gene is comprised of four exons and three introns. Constitutive CXCd mRNA expression was detected in skin, gill, visceral fat, and posterior kidney tissues, while low transcript levels were present in the anterior kidney and spleen. Spleen CXCd transcript abundance increased 1 day after bath vaccination (fourfold) and subsided to basal levels by 7 days postvaccination. Challenge with viable Yersinia ruckeri induced expression of trout CXCd RNA up to ninefold in the spleen. The number of viable Y. ruckeri were significantly correlated with CXCd gene transcript abundance (P = 0.0051, Spearman correlation 0.497, n = 30 fish), and fish with the highest bacterial loads had the highest CXCd expression. In contrast, pro-inflammatory cytokine IL-1-beta2 mRNA levels were elevated in fish infected with low numbers of Y. ruckeri, while diminishing in heavily infected fish. CXCd mRNA expression was not increased in rainbow trout infected with infectious hematopoietic necrosis virus, suggesting that up-regulation may be pathogen-specific. Taken together, these results indicate that CXCd transcript elevation follows the pro-inflammatory cytokine response to Y. ruckeri and may be a relevant immunological marker of exposure.