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

The Impact of Exposure Dosage and Host Genetics on the Shedding Kinetics of Flavobacterium psychrophilum in Rainbow Trout

Flavobacterium psychrophilum, the causative agent of bacterial cold water disease (BCWD), is one of the leading pathogens in rainbow trout (Oncorhynchus mykiss) aquaculture. To date, there is little knowledge of the transmission kinetics of F. psychrophilum over the course of infection. In particular, how transmission is affected by host genotype and pathogen exposure dosage are not well studied. In order to fill in these knowledge gaps, we exposed two divergently selected lines of rainbow trout (ARS-Fp-R and ARS-Fp-S) to a range of dosages of F. psychrophilum (strain CSF117-10). We then measured mortality and bacterial shedding to estimate transmission risk at multiple time points since initial infection. As dosage increased, the number of fish shedding and the amount of bacteria shed increased ranging from 0% to 100% and 103 to 108 cells fish-1 h-1, respectively. In addition, we found that disease resistance (survival) was not correlated with transmission risk blocking, in that 67% of fish which shed bacteria experienced no clinical disease. In general, fish mortality began on Day 3, peaked between Days 5-7 and was higher in the ARS-Fp-R line. Results from this study could be used to develop epidemiological models and improve disease management, particularly in the context of aquaculture and selective breeding.

Systemic immune response of rainbow trout exposed to Flavobacterium psychrophilum infection

Bacterial cold-water disease (BCWD) caused by Flavobacterium psychrophilum is one of the most serious bacterial diseases leading to significant economic loss for rainbow trout (Oncorhynchus mykiss) aquaculture. However, little is known about the systemic immune response of rainbow trout against F. psychrophilum infection. This study investigated the immune response of rainbow trout to F. psychrophilum infection using multiple experiments, including bacterial load detection, phagocyte activity assessment, enzyme activity evaluation, and gene expression profiling. Results showed that the spleen index and intestinal pathogen load reached a peak at 3 days post-infection, with strong pro-inflammatory gene expression observed in rainbow trout. Leukocytes RBA and PKA were significantly elevated in the spleen, blood and intestine at 7 days post-infection. Heat map analysis demonstrated that the spleen had a more substantial pro-inflammatory response compared to the intestine post-infection and exhibited higher expression levels of immune-related genes, including IgM, il1β, il6, cd4, cd8a, cd8b, c1q, chathelicidin, inos, and lysozyme. Both Th1 and Th2 polarized responses in the spleen were activated, with Th2 (il4/13a, gata3) (FC > 4) being more intense than Th1 (tnfα, t-bet) (FC > 2). Tight junction proteins exhibited down-regulation followed by up-regulation post-infection. Collectively, the results of this study expand our current understanding of the immune response of rainbow trout post F. psychrophilum infection but also provide new avenues for investigation in salmonid aquaculture.

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.

Interplay between a bacterial pathogen and its host in rainbow trout isogenic lines with contrasted susceptibility to cold water disease

Infectious diseases are a major constraint on aquaculture. Genetic lines with different susceptibilities to diseases are useful models to identify resistance mechanisms to pathogens and to improve prophylaxis. Bacterial cold-water disease (BCWD) caused by Flavobacterium psychrophilum represents a major threat for freshwater salmonid farming worldwide. A collection of rainbow trout (Oncorhynchus mykiss) isogenic lines was previously produced from a French domestic population. Here, we compared BCWD resistance phenotypes using a subset of isogenic lines chosen for their contrasted susceptibilities to F. psychrophilum. We applied individual monitoring to document the infection process, including time-course quantification of bacteremia and innate immune response. Strikingly, BCWD resistance was correlated with a lower bacterial growth rate in blood. Several immune genes were expressed at higher levels in resistant fish regardless of infection: the Type II arginase (arg2), a marker for M2 macrophages involved in anti-inflammatory responses and tissue repair, and two Toll-like receptors (tlr2/tlr7), responsible for pathogen detection and inflammatory responses. This study highlights the importance of innate and intrinsic defense mechanisms in determining the outcome of F. psychrophilum infections, and illustrates that non-lethal time-course blood sampling for individual monitoring of bacteremia is a powerful tool to resolve within-host pathogen behavior in bacterial fish diseases.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Genome-wide identification of antisense lncRNAs and their association with susceptibility to Flavobacterium psychrophilum in rainbow trout

Eukaryotic genomes encode long noncoding natural antisense transcripts (lncNATs) that have been increasingly recognized as regulatory members of gene expression. Recently, we identified a few antisense transcripts correlating in expression with immune-related genes. However, a systematic genome-wide analysis of lncNATs in rainbow trout is lacking. This study used 134 RNA-Seq datasets from five different projects to identify antisense transcripts. A total of 13,503 lncNATs were identified genome-wide. About 75% of lncNATs showed multiple exons compared to 36.5% of the intergenic lncRNAs. RNA-Seq datasets from resistant, control, and susceptible rainbow trout genetic lines with significant differences in survival rate following Flavobacterium psychrophilum (Fp) infection were analyzed to investigate the potential role of the lncNATs during infection. Twenty-four pairwise comparisons between the different genetic lines, infectious status, and time points revealed 581 differentially expressed (DE) lncNATs and 179 differentially used exons (DUEs). Most of the DE lncNATs strongly and positively correlated in expression with their corresponding sense transcripts across 24 RNA-Seq datasets. LncNATs complementary to genes related to immunity, muscle contraction, proteolysis, and iron/heme metabolism were DE following infection. LncNATs complementary to hemolysis-related genes were DE in the resistant fish compared to susceptible fish on day 5 post-infection, suggesting enhanced clearance of free hemoglobin (Hb) and heme and increased erythropoiesis. LncNATs complementary to hepcidin, a master negative regulator of the plasma iron concentration, were the most downregulated lncNATs on day 5 of bacterial infection in the resistant fish. Ninety-four DE lncNAT, including five complementary to hepcidin, are located within 26 QTL regions previously identified in association with bacterial cold water disease (BCWD) in rainbow trout. Collectively, lncNATs are involved in the molecular architecture of fish immunity and should be further investigated for potential applications in genomic selection and genetic manipulation in aquaculture.

Two functionally distinct heme/iron transport systems are virulence determinants of the fish pathogen Flavobacterium psychrophilum

Bacterial pathogens have a critical impact on aquaculture, a sector that accounts for half of the human fish consumption. Flavobacterium psychrophilum (phylum Bacteroidetes) is responsible for bacterial cold-water disease in salmonids worldwide. The molecular factors involved in host invasion, colonization and haemorrhagic septicaemia are mostly unknown. In this study, we identified two new TonB-dependent receptors, HfpR and BfpR, that are required for adaptation to iron conditions encountered during infection and for virulence in rainbow trout. Transcriptional analyses revealed that their expression is tightly controlled and upregulated under specific iron sources and concentrations. Characterization of deletion mutants showed that they act without redundancy: BfpR is required for optimal growth in the presence of high haemoglobin level, while HfpR confers the capacity to acquire nutrient iron from haem or haemoglobin under iron scarcity. The gene hfpY, co-transcribed with hfpR, encodes a protein related to the HmuY family. We demonstrated that HfpY binds haem and contributes significantly to host colonization and disease severity. Overall, these results are consistent with a model in which both BfpR and Hfp systems promote haem uptake and respond to distinct signals to adapt iron acquisition to the different stages of pathogenesis. Our findings give insight into the molecular basis of pathogenicity of a serious pathogen belonging to the understudied family Flavobacteriaceae and point to the newly identified haem receptors as promising targets for antibacterial development.

Time-course transcriptome analyses of spleen in rainbow trout (Oncorhynchus mykiss) post-Flavobacterium psychrophilum infection

Flavobacterium psychrophilum, the etiological agent of bacterial coldwater disease and rainbow trout fry syndrome, causes considerable losses in salmonid aquaculture globally. Systemic F. psychrophilum infections in rainbow trout (Oncorhynchus mykiss) lead to a range of clinical signs, including ulcerative lesions in the skin and muscle and splenitis. Previous studies offered an integrative analysis of the skeletal muscle response to F. psychrophilum infection in rainbow trout. However, little is known about the molecular mechanism of immune response in the spleen, which is an important immune organ of rainbow trout. Here, we investigated the time-course splenic transcriptome profiles in uninfected rainbow trout (CK) and F. psychrophilum–infected rainbow trout at day 3 and day 7 (D3, D7) by RNA-seq analyses. Among the 7,170 differentially expressed genes (DEGs) in the three comparisons (D3 vs. CK, D7 vs. CK, D3 vs. D7), 1,286 DEGs showed consistent upregulation or downregulation at D3 and D7 and were associated with pattern recognition, acute-phase response, complement cascade, chemokine and cytokine signaling, and apoptosis. The Real time quantitative PCR (RT-qPCR) analysis of eight DEGs confirmed the accuracy of the RNA-Sequencing (RNA-seq) data. Our results reflected a general process from pathogen recognition to inflammatory cytokine generation and delineated a putative Toll-like receptor signaling pathway in rainbow trout spleen, following F. psychrophilum infection. Taken together, these results provide new insights into the molecular mechanism of the immune response to F. psychrophilum infection and are a valuable resource for future research on the prevention and control of bacterial coldwater disease during salmon culture.

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.

Recurrent expansions of B30.2-associated immune receptor families in fish

B30.2 domains, also known as PRY/SPRY, are key components of specific subsets of two large families of proteins involved in innate immunity: the tripartite motif proteins (TRIMs) and the Nod-like receptors (NLRs). TRIM proteins are important, often inducible factors of antiviral innate immunity, targeting multiple steps of viral cycles through a variety of mechanisms. NLRs prime and regulate systemic innate defenses, especially against bacteria, and control inflammation. Large TRIM and NLR subsets characterized by the presence of a B30.2 domain have been reported from a few fish species including zebrafish and seem to be strongly prone to gene duplication/expansion. Here, we performed a large-scale survey of these receptors across about 150 fish genomes, focusing on ray-finned fishes. We assessed the number and genomic distribution of domains and domain combinations associated with TRIMs, NLRs, and other genes containing B30.2 domains and looked for gene expansion patterns across fish groups. We then used a model to test the impact of taxonomy, genome size, and environmental variables on the copy numbers of these genes. Our findings reveal novel domain structures, clade-specific gains and losses. They also assist with the timing of the gene expansions, reveal patterns associated with the MHC, and lay the groundwork for further studies delving deeper into the forces that drive the copy number variation of immune genes on a species level.

PacBio Iso-Seq Improves the Rainbow Trout Genome Annotation and Identifies Alternative Splicing Associated With Economically Important Phenotypes

Rainbow trout is an important model organism that has received concerted international efforts to study the transcriptome. For this purpose, short-read sequencing has been primarily used over the past decade. However, these sequences are too short of resolving the transcriptome complexity. This study reported a first full-length transcriptome assembly of the rainbow trout using single-molecule long-read isoform sequencing (Iso-Seq). Extensive computational approaches were used to refine and validate the reconstructed transcriptome. The study identified 10,640 high-confidence transcripts not previously annotated, in addition to 1,479 isoforms not mapped to the current Swanson reference genome. Most of the identified lncRNAs were non-coding variants of coding transcripts. The majority of genes had multiple transcript isoforms (average ∼3 isoforms/locus). Intron retention (IR) and exon skipping (ES) accounted for 56% of alternative splicing (AS) events. Iso-Seq improved the reference genome annotation, which allowed identification of characteristic AS associated with fish growth, muscle accretion, disease resistance, stress response, and fish migration. For instance, an ES in GVIN1 gene existed in fish susceptible to bacterial cold-water disease (BCWD). Besides, under five stress conditions, there was a commonly regulated exon in prolyl 4-hydroxylase subunit alpha-2 (P4HA2) gene. The reconstructed gene models and their posttranscriptional processing in rainbow trout provide invaluable resources that could be further used for future genetics and genomics studies. Additionally, the study identified characteristic transcription events associated with economically important phenotypes, which could be applied in selective breeding.

Polystyrene microplastics reduce abundance of developing B cells in rainbow trout (Oncorhynchus mykiss) primary cultures

Environmental microplastic pollution (including polystyrene, PS) may have detrimental effects on the health of aquatic organisms. Accumulation of PS microplastics has been reported to affect innate immune cells and inflammatory responses in fish. To date, knowledge on effects of microplastics on the antibody response is still very limited. Here, we investigated effects of small (0.8-20 μm) PS microplastics on the abundance of B lineage cells in primary cultures of developing immune cells from the anterior kidney of rainbow trout. Both purchased PS microbeads and PS microparticles generated from consumer products were used as microplastic sources. We first show that rainbow trout phagocytic B cells efficiently took up small (0.83-3.1 μm) PS microbeads within hours of exposure. In addition, our data revealed that PS microplastic exposure most significantly decreased the abundance of a population of non-phagocytic developing B cells, using both flow cytometry and RT-qPCR. PS microplastics-induced loss of developing B cells further correlated with reduced gene expression of RAG1 and the membrane form of immunoglobulin heavy chains mu and tau. Based on the induced loss of developing B cells observed in our in vitro studies, we speculate that in vivo, chronic PS microplastic-exposure may lead to suboptimal IgM/IgT levels in response to pathogens in teleost species. Considering the highly conserved nature of vertebrate B lymphopoiesis it is likely that PS microplastics will similarly reduce antibody responses in higher vertebrate species, including humans. Further, RAG1 provides an effective biomarker to determine effects of PS microplastics on B cell development in teleost species.

A BCWD-Resistant line of rainbow trout is less sensitive to cortisol implant-induced changes in IgM response as compared to a susceptible (control) line

In salmonids, stress responses increase cortisol levels and disease susceptibility, including to Flavobacterium psychrophilum (Fp), the causative agent of BCWD. A BCWD-resistant line (R-line) of rainbow trout was used here to investigate potential differences in immunoglobulin response after a combined treatment of cortisol and Fp, as compared to a susceptible (S-line) control line. Expression of membrane and secreted immunoglobulin heavy chains mu and tau were determined by RT-qPCR in spleen and anterior kidney. Cortisol treatment did not affect B cell development in the anterior kidney, but delayed IgM responses at the early stage of infection in the spleen of both lines. An earlier IgM response was a determining factor in differential disease progression between resistant- and susceptible fish after Fp-challenge. S-line fish had a delayed and exacerbated IgM response after cortisol implant indicative of a detrimental cycle of sustained IgM responses and high pathogen loads. In contrast, R-line fish had delayed but milder, and protective IgM responses and cleared pathogen successfully. Fp challenge after cortisol implant increased serum cortisol levels on days 3 and 5 compared to mock treatments in S-line fish, but only on day 3 in R-line. Hence, combined cortisol treatment and Fp challenge differentially modulated B cell activation and Fp loads in BCWD-resistant and susceptible lines of rainbow trout. We propose that under conditions of increased stress, BCWD-resistant fish may remain immunologically better equipped to respond to infections compared to BCWD susceptible fish.

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.

Immune gene expression and genome-wide association analysis in rainbow trout with different resistance to Yersinia ruckeri infection

Selective breeding programmes involving marker assisted selection of innately pathogen resistant strains of rainbow trout rely on reliable controlled infection studies, extensive DNA typing of individual fish and recording of expression of relevant genes. We exposed juvenile rainbow trout (6 h bath to 2.6 × 10⁵ CFU mL^-1) to the fish pathogen Yersinia ruckeri serotype O1, biotype 2, eliciting Enteric Red Mouth Disease ERM, and followed the disease progression over 21 days. Cumulative mortality reached 42% at 12 days post challenge (dpc) after which no disease signs were recorded. All fish were sampled for DNA-typing (50 k SNP chip, Affymetrix®) throughout the course of infection when they showed clinical signs of disease (susceptible fish) or at day 21 when fish showed no clinical signs of disease (survivors - resistant fish). Genome-wide association analyses of 1027 trout applying single nucleotide polymorphisms (SNPs) as markers revealed an association between traits (susceptible/resistant) and certain regions of the trout genome. It was indicated that multiple genes are involved in rainbow trout resistance towards ERM whereby it is considered a polygenic trait. A corresponding trout group was kept as non-exposed controls and a comparative expression analysis of central innate and adaptive immune genes in gills, spleen and liver was performed for three fish groups: 1) moribund trout exhibiting clinical signs 7 dpc (CS), 2) exposed fish without clinical signs at the same sampling point (NCS) and 3) surviving fish at 21 dpc (survivors). Immune genes encoding inflammatory cytokines (IL-1β, IL-2A, IL-6A, IL-8, IL-10A, IL-12, IL-17A/F2A, IL-17C1, IL-17C2, IL-22, IFNγ, TNFα), acute phase reactants (SAA, C3, cathelicidins, lysozyme) were expressed differently in CS and NCS fish. Correlation (negative or positive) between expression of genes and bacterial load suggested involvement of immune genes in protection. Down-regulation of adaptive immune genes including IgDm, IgDs, IgT and TCR-β was seen primarily in CS and NCS fish whereas survivors showed up-regulation of effector molecule genes such as cathelicidins, complement and lysozyme suggesting their role in clearing the infection. In conclusion, SNP analyses indicated that ERM resistance in rainbow trout is a multi-locus trait. The gene expression in surviving fish suggested that several immune genes are associated with the trait conferring resistance.

Salmonid Antibacterial Immunity: An Aquaculture Perspective

The aquaculture industry is continuously threatened by infectious diseases, including those of bacterial origin. Regardless of the disease burden, aquaculture is already the main method for producing fish protein, having displaced capture fisheries. One attractive sector within this industry is the culture of salmonids, which are (a) uniquely under pressure due to overfishing and (b) the most valuable finfish per unit of weight. There are still knowledge gaps in the understanding of fish immunity, leading to vaccines that are not as effective as in terrestrial species, thus a common method to combat bacterial disease outbreaks is the use of antibiotics. Though effective, this method increases both the prevalence and risk of generating antibiotic-resistant bacteria. To facilitate vaccine design and/or alternative treatment efforts, a deeper understanding of the teleost immune system is essential. This review highlights the current state of teleost antibacterial immunity in the context of salmonid aquaculture. Additionally, the success of current techniques/methods used to combat bacterial diseases in salmonid aquaculture will be addressed. Filling the immunology knowledge gaps highlighted here will assist in reducing aquaculture losses in the future.

Transcriptomic Response to Selective Breeding for Fast Growth in Rainbow Trout (Oncorhynchus mykiss)

Genetic improvement for faster growth is a conventional approach to increase growth rates in aquaculture species; however, the genetic and physiological factors regulating growth performance in fish are not fully characterized. The objective of this study was to identify physiological mechanisms associated with faster growth rates by comparing the liver and muscle transcriptome of a rainbow trout line selectively bred for fast growth (growth line, GL) and a contemporary randomly mated control line (synthetic control, SC) from the same selective breeding program. A third genetic line from a commercial egg supplier (commercial A, CA) was also included to characterize differences in gene expression profiles between populations. Body weight of the GL at harvest was approximately 20% and 8% heavier (p < 0.05) than SC and CA, respectively. There were 145 and 36 differentially expressed genes (DEG) in liver and white muscle, respectively, between the GL and SC that were enriched for the growth hormone/insulin-like growth factor axis (GH/IGF) and PI3K-Akt, JAK-STAT, MAPK, and cAMP signal transduction pathways. A greater concentration of plasma IGF-I was detected in the GL compared with SC (p < 0.05). A unique gene profile was detected in CA, with 11 and 210 DEG in liver and white muscle; these genes associated with innate immunity, complement systems, and metabolic pathways. Collectively, these findings provide a more extensive characterization of the fast-growth phenotype in fish that furthers knowledge of the physiological basis for genetic variation in growth performance in selectively bred rainbow trout.

Cross-protection of a live-attenuated Flavobacterium psychrophilum immersion vaccine against novel Flavobacterium spp. and Chryseobacterium spp. strains

For salmonid producers, a common threat is Flavobacterium psychrophilum. Recent advancements in bacterial coldwater disease (BCWD) management include the development of a live-attenuated immersion vaccine that cross-protects against an array of F. psychrophilum strains. Emerging family Flavobacteriaceae cases associated with clinical disease have been increasing, including pathogenic isolates of Flavobacterium spp. and Chryseobacterium spp. The cross-protective ability of a live-attenuated F. psychrophilum vaccine was determined against three virulent Flavobacteriaceae isolates. Juvenile rainbow trout were vaccinated, developed high F. psychrophilum-specific antibody titres and were challenged with Chryseobacterium spp. isolates (S25 and T28), a Flavobacterium sp. (S21) isolate, a mixed combination of S21:S25:T28, and a standard virulent F. psychrophilum CSF259-93 strain. Results demonstrated strong protection in the CSF259-93 vaccinated group (relative per cent survival (RPS)=94.44%) when compared to the relevant CSF259-93 controls (p < .001). Protection was also observed for vaccinated fish challenged with the S21:S25:T28 mix (RPS = 85.18%; p < .001). However, protection was not observed with the S21, S25 or T28 isolates alone. Analysis of whole-cell lysates revealed differences in protein banding by SDS-PAGE, but conserved antigenic regions by Western blot in S25 and T28. Results demonstrate that this live-attenuated vaccine provided protection against mixed flavobacterial infection and suggest further benefits against flavobacteriosis.

Combining Multiple Approaches and Models to Dissect the Genetic Architecture of Resistance to Infections in Fish

Infectious diseases represent a major threat for the sustainable development of fish farming. Efficient vaccines are not available against all diseases, and growing antibiotics resistance limits the use of antimicrobial drugs in aquaculture. It is therefore important to understand the basis of fish natural resistance to infections to help genetic selection and to develop new approaches against infectious diseases. However, the identification of the main mechanisms determining the resistance or susceptibility of a host to a pathogenic microbe is challenging, integrating the complexity of the variation of host genetics, the variability of pathogens, and their capacity of fast evolution and adaptation. Multiple approaches have been used for this purpose: (i) genetic approaches, QTL (quantitative trait loci) mapping or GWAS (genome-wide association study) analysis, to dissect the genetic architecture of disease resistance, and (ii) transcriptomics and functional assays to link the genetic constitution of a fish to the molecular mechanisms involved in its interactions with pathogens. To date, many studies in a wide range of fish species have investigated the genetic determinism of resistance to many diseases using QTL mapping or GWAS analyses. A few of these studies pointed mainly toward adaptive mechanisms of resistance/susceptibility to infections; others pointed toward innate or intrinsic mechanisms. However, in the majority of studies, underlying mechanisms remain unknown. By comparing gene expression profiles between resistant and susceptible genetic backgrounds, transcriptomics studies have contributed to build a framework of gene pathways determining fish responsiveness to a number of pathogens. Adding functional assays to expression and genetic approaches has led to a better understanding of resistance mechanisms in some cases. The development of knock-out approaches will complement these analyses and help to validate putative candidate genes critical for resistance to infections. In this review, we highlight fish isogenic lines as a unique biological material to unravel the complexity of host response to different pathogens. In the future, combining multiple approaches will lead to a better understanding of the dynamics of interaction between the pathogen and the host immune response, and contribute to the identification of potential targets of selection for improved resistance.

Transcriptome Profiling of Pacu (Piaractus mesopotamicus) Challenged With Pathogenic Aeromonas hydrophila: Inference on Immune Gene Response

Pacu (Piaractus mesopotamicus) is a Neotropical fish of major importance for South American aquaculture. Septicemia caused by Aeromonas hydrophila bacteria is currently considered a substantial threat for pacu aquaculture that have provoked infectious disease outbreaks with high economic losses. The understanding of molecular aspects on progress of A. hydrophila infection and pacu immune response is scarce, which have limited the development of genomic selection for resistance to this infection. The present study aimed to generate information on transcriptome of pacu in face of A. hydrophila infection, and compare the transcriptomic responses between two groups of time-series belonging to a disease resistance challenge, peak mortality (HM) and mortality plateau (PM) groups of individuals. Nine RNA sequencing (RNA-Seq) libraries were prepared from liver tissue of challenged individuals, generating ∼160 million 150 bp pair-end reads. After quality trimming/cleanup, these reads were assembled de novo generating 211,259 contigs. When the expression of genes from individuals of HM group were compared to individuals from control group, a total of 4,413 differentially expressed transcripts were found (2,000 upregulated and 2,413 downregulated candidate genes). Additionally, 433 transcripts were differentially expressed when individuals from MP group were compared with those in the control group (155 upregulated and 278 downregulated candidate genes). The resulting differentially expressed transcripts were clustered into the following functional categories: cytokines and signaling, epithelial protection, antigen processing and presentation, apoptosis, phagocytosis, complement system cascades and pattern recognition receptors. The proposed results revealing relevant differential gene expression on HM and PM groups which will contribute to a better understanding of the molecular defense mechanisms during A. hydrophila infection.

Understanding the pathogenesis of Flavobacterium psychrophilum using the rainbow trout monocyte/macrophage-like cell line, RTS11, as an infection model

The life cycle of Flavobacterium psychrophilum (Fp), the causative agent of bacterial coldwater disease (BCWD) and rainbow trout fry syndrome (RTFS), appears to involve interactions with spleen and head kidney macrophages. To develop an in vitro model for studying this, F. psychrophilum was incubated with a rainbow trout splenic monocyte/macrophage-like cell line (RTS11) and fundamental macrophage functions evaluated. The animal cell basal medium, L15, supplemented with bovine serum (FBS) supports RTS11 maintenance, and surprisingly, L15 with 2% FBS (L15/FBS) also supported F. psychrophilum growth. L15/FBS in which the bacteria had been grown is referred to as F. psychrophilum conditioned medium (FpCM). Adding FpCM to RTS11 cultures caused a small, yet significant, percentage of cells to die, many cells to become more diffuse, and phagocytosis to be temporarily reduced. FpCM also significantly stimulated transcript expression for pro-inflammatory cytokines (IL-1β, TNFα and IL-6) and the anti-inflammatory cytokine (IL-10) after one day of exposure but this upregulation rapidly declined over time. Adding live F. psychrophilum to RTS11 cultures also altered the cellular morphology and stimulated cytokine expression more profoundly than FpCM. Additionally, the phagocytic activity of RTS11 was also significantly impaired by live F. psychrophilum, but not to the same extent as when exposed to FpCM. Adding heat-killed bacteria to RTS11 cultures elicited few changes. These bacteria/RTS11 co-cultures should be useful for gaining a deeper understanding of the pathogenesis of F. psychrophilum and may aid in the development of effective measures to prevent infection and spread of this troublesome disease.

Transcriptomic response of rainbow trout (Oncorhynchus mykiss) skeletal muscle to Flavobacterium psychrophilum

Flavobacterium psychrophilum is the etiologic agent of rainbow trout fry syndrome (RTFS). This pathogen infects a wide variety of salmonid species during freshwater stages, causing significant losses in the aquaculture industry. Rainbow trout (Oncorhynchus mykiss) infected with F. psychrophilum, presents as the main external clinical sign ulcerative lesions and necrotic myositis in skeletal muscle. We previously reported the in vitro cytotoxic activity of F. psychrophilum on rainbow trout myoblast, however little is known about the molecular mechanisms underlying the in vivo pathogenesis in skeletal muscle. In this study, we examined the transcriptomic profiles of skeletal muscle tissue of rainbow trout intraperitoneally challenged with low infection dose of F. psychrophilum. Using high-throughput RNA-seq, we found that 233 transcripts were up-regulated, mostly associated to ubiquitin mediated proteolysis and apoptosis. Conversely, 189 transcripts were down-regulated, associated to skeletal muscle contraction. This molecular signature was consistent with creatine kinase activity in plasma and oxidative damage in skeletal muscle. Moreover, the increased caspase activity suggests as a whole skeletal muscle atrophy induced by F. psychrophilum. This study offers an integrative analysis of the skeletal muscle response to F. psychrophilum infection and reveals unknown aspects of its pathogenesis in rainbow trout.

PACAP Is Lethal to Flavobacterium psychrophilum Through Either Direct Membrane Permeabilization or Indirectly, by Priming the Immune Response in Rainbow Trout Macrophages

Pituitary adenylate cyclase-activating polypeptide (PACAP) is a multifunctional neuropeptide that is widely distributed in mammals and is capable of performing roles as a neurotransmitter, neuromodulator, and vasodilator. This polypeptide belongs to the glucagon/secretin superfamily, of which some members have been shown to act as antimicrobial peptides in both mammalian and aquatic organisms. In teleosts, PACAP has been demonstrated to have direct antimicrobial activity against several aquatic pathogens, yet this phenomenon has never been studied throughout a live bacterial challenge. The present study focuses on the influence of synthetic Clarias gariepinus 38 amino acid PACAP on the rainbow trout monocyte/macrophage-like cell line, RTS11, when exposed to the coldwater bacterial pathogen Flavobacterium psychrophilum. PACAP was shown to have direct antimicrobial activity on F. psychrophilum when grown in both cytophaga broth and cell culture media (L-15). Further, the ability of teleostean PACAP to permeabilize the membrane of an aquatic pathogen, F. psychrophilum, was demonstrated for the first time. The viability of RTS11 when exposed to PACAP was also observed using a trypan blue exclusion assay to determine optimal experimental doses of the antimicrobial peptide. This displayed that only concentrations higher than 0.1 μM negatively impacted RTS11 survival. Interestingly, when RTS11 was pre-treated with PACAP for 24 h before experiencing infection with live F. psychrophilum, growth of the pathogen was severely inhibited in a dose-dependent manner when compared to cells receiving no pre-treatment with the polypeptide. Relative expression of pro-inflammatory cytokines (IL-1β, TNFα, and IL-6) and PACAP receptors (VPAC1 and PAC1) was also analyzed in RTS11 following PACAP exposure alone and in conjunction with live F. psychrophilum challenge. These qRT-PCR findings revealed that PACAP may have a synergistic effect on RTS11 immune function. The results of this study provide evidence that PACAP has immunostimulatory activity on rainbow trout immune cells as well as antimicrobial activity against aquatic bacterial pathogens such as F. psychrophilum. As there are numerous pathogens that plague the aquaculture industry, PACAP may stimulate the teleost immune system while also providing an efficacious alternative to antibiotic use.

Transcriptomic analysis of immunity in rainbow trout (Oncorhynchus mykiss) gills infected by Ichthyophthirius multifiliis

The parasite Ichthyophthirius multifiliis infecting skin, fins and gills of a wide range of freshwater fish species, including rainbow trout, is known to induce a protective immune response in the host. Although a number of studies have reported activation of several immune genes in infected fish host, the immune response picture is still considered incomplete. In order to address this issue, a comparative transcriptomic analysis was performed on infected versus uninfected rainbow trout gills and it showed that a total of 3352 (7.2%) out of 46,585 identified gene sequences were significantly regulated after parasite infection. Of differentially expressed gene sequences, 1796 genes were up-regulated and 1556 genes were down-regulated. These were classified into 61 Gene Ontology (GO) terms and mapped to 282 reference canonical pathways in the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Infection of I. multifiliis induced a clear differential expression of immune genes, related to both innate and adaptive immunity. A total of 268 (6.86%) regulated gene sequences were known to take part in 16 immune-related pathways. These involved pathways related to the innate immunity such as the Chemokine signaling pathway, Platelet activation, Toll-like receptor signaling pathway, NOD-like receptor signaling pathway, and Leukocyte transendothelial migration. Elevated transcription of genes encoding the TLR 8 gene and chemokines (CCL4, CCL19, CCL28, CXCL8, CXCL11, CXCL13, CXCL14) was recorded indicating their roles in recognition of I. multifiliis and subsequent induction of the inflammatory response, respectively. A number of upregulated genes in infected gills were associated with antigen processing/presentation and T and B cell receptor signaling (including B cell marker CD22 involved in B cell development). Overall the analysis supports the notion that I. multifiliis induces a massive and varied innate response upon which a range of adaptive immune responses are established which may contribute to the long lasting protection of immunized rainbow trout.

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.

Innate immune cell signatures in a BCWD-Resistant line of rainbow trout before and after in vivo challenge with Flavobacterium psychrophilum

Phenotypes of myeloid-lineage cells remain poorly understood in the rainbow trout, and were the focus of this study, including effects of in vivo challenge to Flavobacterium psychrophilum (Fp), the cause of Bacterial Cold Water Disease (BCWD). A genetic line was used that is highly resistant to BCWD (R-line) as well as a susceptible control line (S-line). Using flow cytometry, we describe two Pax5-negative, myeloid-lineage populations: Population 1 consisted of small cells with high SSC and strong staining for Q4E, MPO, Pu1, EBF, and IL- 1β, which we named "neutrophil-like" cell. Population 2 had high Q4E, but weaker MPO, Pu1, EBF, and IL-1β staining. Five days after Fp-challenge, both genetic lines had a reduced abundance of neutrophil-like cells in anterior kidney, PBL, and spleen. Pop. 2 abundance was reduced in anterior kidney, and increased in spleen. S-line fish responded more strongly to Fp-challenge compared to R-line fish. Challenged fish with a higher abundance of neutrophil-like cells had significantly lower Fp-loads after challenge, suggesting that these cells aid in the resistance to BCWD.

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.

Genome-Wide Association Analysis With a 50K Transcribed Gene SNP-Chip Identifies QTL Affecting Muscle Yield in Rainbow Trout

Detection of coding/functional SNPs that change the biological function of a gene may lead to identification of putative causative alleles within QTL regions and discovery of genetic markers with large effects on phenotypes. This study has two-fold objectives, first to develop, and validate a 50K transcribed gene SNP-chip using RNA-Seq data. To achieve this objective, two bioinformatics pipelines, GATK and SAMtools, were used to identify ~21K transcribed SNPs with allelic imbalances associated with important aquaculture production traits including body weight, muscle yield, muscle fat content, shear force, and whiteness in addition to resistance/susceptibility to bacterial cold-water disease (BCWD). SNPs ere identified from pooled RNA-Seq data collected from ~620 fish, representing 98 families from growth- and 54 families from BCWD-selected lines with divergent phenotypes. In addition, ~29K transcribed SNPs without allelic-imbalances were strategically added to build a 50K Affymetrix SNP-chip. SNPs selected included two SNPs per gene from 14K genes and ~5K non-synonymous SNPs. The SNP-chip was used to genotype 1728 fish. The average SNP calling-rate for samples passing quality control (QC; 1,641 fish) was ≥ 98.5%. The second objective of this study was to test the feasibility of using the new SNP-chip in GWA (Genome-wide association) analysis to identify QTL explaining muscle yield variance. GWA study on 878 fish (representing 197 families from 2 consecutive generations) with muscle yield phenotypes and genotyped for 35K polymorphic markers (passing QC) identified several QTL regions explaining together up to 28.40% of the additive genetic variance for muscle yield in this rainbow trout population. The most significant QTLs were on chromosomes 14 and 16 with 12.71 and 10.49% of the genetic variance, respectively. Many of the annotated genes in the QTL regions were previously reported as important regulators of muscle development and cell signaling. No major QTLs were identified in a previous GWA study using a 57K genomic SNP chip on the same fish population. These results indicate improved detection power of the transcribed gene SNP-chip in the target trait and population, allowing identification of large-effect QTLs for important traits in rainbow trout.

Applications of transcriptomics and proteomics in understanding fish immunity

With the development of intensive aquaculture, economic losses increasingly result from fish mortality due to pathogen infection. In recent years, a growing number of researchers have used transcriptomic and proteomic analyses to study fish immune responses to exogenous pathogen infection. Integrating transcriptomic and proteomic analyses provides a better understanding of the fish immune system including gene expression, regulation, and the intricate biological processes underlying immune responses against infection. This review focuses on the recent advances in the fields of transcriptomics and proteomics, which have contributed to our understanding of fish immunity to exogenous pathogens.

Crosstalk among lncRNAs, microRNAs and mRNAs in the muscle 'degradome' of rainbow trout

In fish, protein-coding and noncoding genes involved in muscle atrophy are not fully characterized. In this study, we characterized coding and noncoding genes involved in gonadogenesis-associated muscle atrophy, and investigated the potential functional interplay between these genes. Using RNA-Seq, we compared expression pattern of mRNAs, long noncoding RNAs (lncRNAs) and microRNAs of atrophying skeletal muscle from gravid females and control skeletal muscle from age-matched sterile individuals. A total of 852 mRNAs, 1,160 lncRNAs and 28 microRNAs were differentially expressed (DE) between the two groups. Muscle atrophy appears to be mediated by many genes encoding ubiquitin-proteasome system, autophagy related proteases, lysosomal proteases and transcription factors. Transcripts encoding atrogin-1 and mir-29 showed exceptional high expression in atrophying muscle, suggesting an important role in bulk muscle proteolysis. DE genes were co-localized in the genome with strong expression correlation, and they exhibited extensive 'lncRNA-mRNA', 'lncRNA-microRNA', 'mRNA-microRNA' and 'lncRNA-protein' physical interactions. DE genes exhibiting potential functional interactions comprised the highly correlated 'lncRNA-mRNA-microRNA' gene network described as 'degradome'. This study pinpoints extensive coding and noncoding RNA interactions during muscle atrophy in fish, and provides valuable resources for future mechanistic studies.

Muscle and liver transcriptome characterization and genetic marker discovery in the farmed meagre, Argyrosomus regius

Meagre (Argyrosomus regius), a teleost fish of the family Sciaenidae, is part of a group of marine fish species considered new for Mediterranean aquaculture representing the larger fish cultured in the region. Meagre aquaculture started ~25years ago in West Mediterranean, and the supply of juveniles has been dominated by few hatcheries. This fact has raised concerns on possible inbreeding, urging the need for genetic information on the species and for an assessment of the polymorphisms found in the genome. To that end we characterized the muscle and liver transcriptome of a pool of meagre individuals, from different families and phenotypic size, to obtain a backbone that can support future studies regarding physiology, immunology and genetics of the species. The assembled transcripts were assigned to a wide range of biological processes including growth, reproduction, metabolism, development, stress and behavior. Then, to infer its genetic diversity and provide a catalogue of markers for future use, we scanned the reconstructed transcripts for polymorphic genetic markers. Our search revealed a total of 42,933 high quality SNP and 20,581 STR markers. We found a relatively low rate of polymorphism in the transcriptome that may indicate that inbreeding has taken place. This study has led to a catalogue of genetic markers at the expressed part of the genome and has set the ground for understanding growth and other traits of interest in meagre.

The circadian transcriptome of marine fish (Sparus aurata) larvae reveals highly synchronized biological processes at the whole organism level

The regulation of circadian gene expression remains largely unknown in farmed fish larvae. In this study, a high-density oligonucleotide microarray was used to examine the daily expression of 13,939 unique genes in whole gilthead sea bream (Sparus aurata) larvae with fast growth potentiality. Up to 2,229 genes were differentially expressed, and the first two components of Principal Component Analysis explained more than 81% of the total variance. Clustering analysis of differentially expressed genes identified 4 major clusters that were triggered sequentially, with a maximum expression at 0 h, 3 h, 9–15 h and 18-21 h zeitgeber time. Various core clock genes (per1, per2, per3, bmal1, cry1, cry2, clock) were identified in clusters 1–3, and their expression was significantly correlated with several genes in each cluster. Functional analysis revealed a daily consecutive activation of canonical pathways related to phototransduction, intermediary metabolism, development, chromatin remodeling, and cell cycle regulation. This daily transcriptome of whole larvae resembles a cell cycle (G1/S, G2/M, and M/G1 transitions) in synchronization with multicellular processes, such as neuromuscular development. This study supports that the actively feeding fish larval transcriptome is temporally organized in a 24-h cycle, likely for maximizing growth and development.

Omics and cytokine discovery in fish: Presenting the Yellowtail kingfish (Seriola lalandi) as a case study

A continued programme of research is essential to overcome production bottlenecks in any aquacultured fish species. Since the introduction of genetic and molecular techniques, the quality of immune research undertaken in fish has greatly improved. Thousands of species specific cytokine genes have been discovered, which can be used to conduct more sensitive studies to understand how fish physiology is affected by aquaculture environments or disease. Newly available transcriptomic technologies, make it increasingly easier to study the immunogenetics of farmed species for which little data exists. This paper reviews how the application of transcriptomic procedures such as RNA Sequencing (RNA-Seq) can advance fish research. As a case study, we present some preliminary findings using RNA-Seq to identify cytokine related genes in Seriola lalandi. These will allow in-depth investigations to understand the immune responses of these fish in response to environmental change or disease and help in the development of therapeutic approaches.

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.

RNA-Seq reveals MicroRNA expression signature and genetic polymorphism associated with growth and muscle quality traits in rainbow trout

The role of microRNA expression and genetic variation in microRNA-binding sites of target genes on growth and muscle quality traits is poorly characterized. We used RNA-Seq approach to investigate their importance on 5 growth and muscle quality traits: whole body weight (WBW), muscle yield, muscle crude-fat content, muscle shear force and whiteness. Phenotypic data were collected from 471 fish, representing 98 families (~5 fish/family) from a growth-selected line. Muscle microRNAs and mRNAs were sequenced from 22 families showing divergent phenotypes. Ninety microRNAs showed differential expression between families with divergent phenotypes, and their expression was strongly associated with variation in phenotypes. A total of 204 single nucleotide polymorphisms (SNPs) present in 3′ UTR of target genes either destroyed or created novel illegitimate microRNA target sites; of them, 78 SNPs explained significant variation in the aforementioned 5 muscle traits. Majority of the phenotype-associated SNPs were present in microRNA-binding sites of genes involved in energy metabolism and muscle structure. These findings suggest that variation in microRNA expression and/or sequence variation in microRNA binding sites in target genes play an important role in mediating differences in fish growth and muscle quality phenotypes.

Aquaculture genomics, genetics and breeding in the United States: current status, challenges, and priorities for future research

Advancing the production efficiency and profitability of aquaculture is dependent upon the ability to utilize a diverse array of genetic resources. The ultimate goals of aquaculture genomics, genetics and breeding research are to enhance aquaculture production efficiency, sustainability, product quality, and profitability in support of the commercial sector and for the benefit of consumers. In order to achieve these goals, it is important to understand the genomic structure and organization of aquaculture species, and their genomic and phenomic variations, as well as the genetic basis of traits and their interrelationships. In addition, it is also important to understand the mechanisms of regulation and evolutionary conservation at the levels of genome, transcriptome, proteome, epigenome, and systems biology. With genomic information and information between the genomes and phenomes, technologies for marker/causal mutation-assisted selection, genome selection, and genome editing can be developed for applications in aquaculture. A set of genomic tools and resources must be made available including reference genome sequences and their annotations (including coding and non-coding regulatory elements), genome-wide polymorphic markers, efficient genotyping platforms, high-density and high-resolution linkage maps, and transcriptome resources including non-coding transcripts. Genomic and genetic control of important performance and production traits, such as disease resistance, feed conversion efficiency, growth rate, processing yield, behaviour, reproductive characteristics, and tolerance to environmental stressors like low dissolved oxygen, high or low water temperature and salinity, must be understood. QTL need to be identified, validated across strains, lines and populations, and their mechanisms of control understood. Causal gene(s) need to be identified. Genetic and epigenetic regulation of important aquaculture traits need to be determined, and technologies for marker-assisted selection, causal gene/mutation-assisted selection, genome selection, and genome editing using CRISPR and other technologies must be developed, demonstrated with applicability, and application to aquaculture industries.Major progress has been made in aquaculture genomics for dozens of fish and shellfish species including the development of genetic linkage maps, physical maps, microarrays, single nucleotide polymorphism (SNP) arrays, transcriptome databases and various stages of genome reference sequences. This paper provides a general review of the current status, challenges and future research needs of aquaculture genomics, genetics, and breeding, with a focus on major aquaculture species in the United States: catfish, rainbow trout, Atlantic salmon, tilapia, striped bass, oysters, and shrimp. While the overall research priorities and the practical goals are similar across various aquaculture species, the current status in each species should dictate the next priority areas within the species. This paper is an output of the USDA Workshop for Aquaculture Genomics, Genetics, and Breeding held in late March 2016 in Auburn, Alabama, with participants from all parts of the United States.

De novo assembly of Sockeye salmon kidney transcriptomes reveal a limited early response to piscine reovirus with or without infectious hematopoietic necrosis virus superinfection

Background

Piscine reovirus (PRV) has been associated with the serious disease known as Heart and Skeletal Muscle Inflammation (HSMI) in cultured Atlantic salmon Salmo salar in Norway. PRV is also prevalent in wild and farmed salmon without overt disease manifestations, suggesting multifactorial triggers or PRV variant-specific factors are required to initiate disease. In this study, we explore the head kidney transcriptome of Sockeye salmon Oncorhynchus nerka during early PRV infection to identify host responses in the absence of disease in hopes of elucidating mechanisms by which PRV may directly alter host functions and contribute to the development of a disease state. We further investigate the role of PRV as a coinfecting agent following superinfection with infectious hematopoietic necrosis virus (IHNV) – a highly pathogenic rhabdovirus endemic to the west coast of North America.

Results

Challenge of Sockeye salmon with PRV resulted in high quantities of viral transcripts to become present in the blood and kidney of infected fish without manifestations of disease. De novo transcriptome assembly of over 2.3 billion paired RNA-seq reads from the head kidneys of 36 fish identified more than 320,000 putative unigenes, of which less than 20 were suggested to be differentially expressed in response to PRV at either 2 or 3 weeks post challenge by DESeq2 and edgeR analysis. Of these, only one, Ependymin, was confirmed to be differentially expressed by qPCR in an expanded sample set. In contrast, IHNV induced substantial transcriptional changes (differential expression of > 20,000 unigenes) which included transcripts involved in antiviral and inflammatory response pathways. Prior infection with PRV had no significant effect on host responses to superinfecting IHNV, nor did host responses initiated by IHNV exposure influence increasing PRV loads.

Conclusions

PRV does not substantially alter the head kidney transcriptome of Sockeye salmon during early (2 to 3 week) infection and dissemination in a period of significant increasing viral load, nor does the presence of PRV change the host transcriptional response to an IHNV superinfection. Further, concurrent infections of PRV and IHNV do not appear to significantly influence the infectivity or severity of IHNV associated disease, or conversely, PRV load.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3196-y) 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.

Sea lampreys elicit strong transcriptomic responses in the lake trout liver during parasitism

BACKGROUND

The sea lamprey (Petromyzon marinus) is a jawless vertebrate that parasitizes fish as an adult and, with overfishing, was responsible for the decline in lake trout (Salvelinus namaycush) populations in the Great Lakes. While laboratory studies have looked at the rates of wounding on various fish hosts, there have been few investigations on the physiological effects of lamprey wounding on the host. In the current study, two morphotypes of lake trout, leans and siscowets, were parasitized in the laboratory by sea lampreys and the liver transcriptomes of parasitized and nonparasitized fish were analyzed by RNA-seq (DESeq2 and edgeR) to determine which genes and gene pathways (Ingenuity Pathway Analysis) were altered by lamprey parasitism.

RESULTS

Overall, genes encoding molecules involved in catalytic (e.g., enzymatic) and binding activities (factors and regulators) predominated the regulated gene lists. In siscowets, the top upregulated gene was growth arrest and DNA-damage-inducible protein and for leans it was interleukin-18-binding protein. In leans, the most significantly downregulated gene was UDP-glucuronosyltransferase 2A2 - DESeq2 or phosphotriesterase related - edgeR. For siscowets, the top downregulated gene was C-C motif chemokine 19 - DESeq2 or GTP-binding protein Rhes - edgeR. Gene pathways associated with inflammatory-related responses or factors (cytokines, chemokines, oxidative stress, apoptosis) were regulated following parasitism in both morphotypes. However, pathways related to energy metabolism (glycolysis, gluconeogenesis, lipolysis, lipogenesis) were also regulated. These pathways or the intensity or direction (up/downregulation) of regulation were different between leans and siscowets. Finally, one of the most significantly downregulated pathways in both leans and siscowets was the kynurenine (tryptophan degradation) pathway.

CONCLUSIONS

The results indicate a strong transcriptional response in the lake trout to lamprey parasitism that entails genes involved in the regulation of inflammation and cellular damage. Responses to energy utilization as well as hydromineral balance also occurred indicating an adjustment in the host to energy demands and osmotic imbalances during parasitism. Given the role of the kynurenine pathway in promoting immunotolerance in mammals, the downregulation observed in this pathway during parasitism may signify an attempt by the host to inhibit any feedback suppression of the immune response to the lamprey.

Detection and Validation of QTL Affecting Bacterial Cold Water Disease Resistance in Rainbow Trout Using Restriction-Site Associated DNA Sequencing

Bacterial cold water disease (BCWD) causes significant economic loss in salmonid aquaculture. Using microsatellite markers in a genome scan, we previously detected significant and suggestive QTL affecting phenotypic variation in survival following challenge with Flavobacterium psychrophilum, the causative agent of BCWD in rainbow trout. In this study, we performed selective genotyping of SNPs from restriction-site associated DNA (RAD) sequence data from two pedigreed families (2009070 and 2009196) to validate the major QTL from the previous work and to detect new QTL. The use of RAD SNPs in the genome scans increased the number of mapped markers from ~300 to ~5,000 per family. The significant QTL detected in the microsatellites scan on chromosome Omy8 in family 2009070 was validated explaining up to 58% of the phenotypic variance in that family, and in addition, a second QTL was also detected on Omy8. Two novel QTL on Omy11 and 14 were also detected, and the previously suggestive QTL on Omy1, 7 and 25 were also validated in family 2009070. In family 2009196, the microsatellite significant QTL on Omy6 and 12 were validated and a new QTL on Omy8 was detected, but none of the previously detected suggestive QTL were validated. The two Omy8 QTL from family 2009070 and the Omy12 QTL from family 2009196 were found to be co-localized with handling and confinement stress response QTL that our group has previously identified in a separate pedigreed family. With the currently available data we cannot determine if the co-localized QTL are the result of genes with pleiotropic effects or a mere physical proximity on the same chromosome segment. The genetic markers linked to BCWD resistance QTL were used to query the scaffolds of the rainbow trout reference genome assembly and the QTL-positive scaffold sequences were found to include 100 positional candidate genes. Several of the candidate genes located on or near the two Omy8 QTL detected in family 2009070 suggest potential linkages between stress response and the regulation of immune response in rainbow trout.

Transcriptome assembly, gene annotation and tissue gene expression atlas of the rainbow trout

Efforts to obtain a comprehensive genome sequence for rainbow trout are ongoing and will be complemented by transcriptome information that will enhance genome assembly and annotation. Previously, transcriptome reference sequences were reported using data from different sources. Although the previous work added a great wealth of sequences, a complete and well-annotated transcriptome is still needed. In addition, gene expression in different tissues was not completely addressed in the previous studies. In this study, non-normalized cDNA libraries were sequenced from 13 different tissues of a single doubled haploid rainbow trout from the same source used for the rainbow trout genome sequence. A total of ~1.167 billion paired-end reads were de novo assembled using the Trinity RNA-Seq assembler yielding 474,524 contigs > 500 base-pairs. Of them, 287,593 had homologies to the NCBI non-redundant protein database. The longest contig of each cluster was selected as a reference, yielding 44,990 representative contigs. A total of 4,146 contigs (9.2%), including 710 full-length sequences, did not match any mRNA sequences in the current rainbow trout genome reference. Mapping reads to the reference genome identified an additional 11,843 transcripts not annotated in the genome. A digital gene expression atlas revealed 7,678 housekeeping and 4,021 tissue-specific genes. Expression of about 16,000–32,000 genes (35–71% of the identified genes) accounted for basic and specialized functions of each tissue. White muscle and stomach had the least complex transcriptomes, with high percentages of their total mRNA contributed by a small number of genes. Brain, testis and intestine, in contrast, had complex transcriptomes, with a large numbers of genes involved in their expression patterns. This study provides comprehensive de novo transcriptome information that is suitable for functional and comparative genomics studies in rainbow trout, including annotation of the genome.