Gene expression comparison of resistant and susceptible Atlantic salmon fry challenged with Infectious Pancreatic Necrosis virus reveals a marked contrast in immune response

First insight of the genome-wide association study and genomic prediction into enteritis disease (Vibrio harveyi) resistance trait in the lined seahorse (Hippocampus erectus)

Enteritis caused by Vibrio is a highly die-off disease that severely impeded substantial production in seahorse aquaculture. In the present study, challenged with LD50 of concentration of Vibrio harveyi, a total of 161 of susceptible and 166 of resistant individuals were allocated into binary survival phenotypes, thus, to firstly investigate the genetic architecture by genome-wide association study (GWAS) analysis, as well as to evaluate the feasibility of genomic selection (GS) in enteritis disease resistance trait of the lined seahorse Hippocampus erectus. Results indicated that the heritability for resistance to Vibrio harveyi was estimated to be 0.10. And a set of 10 significant/suggestive SNPs in a multiple chromosomes localization were identified, explaining 7.76% to 13.28% of genetic variance. Associated 82 of candidate genes were clustered into signal transduction, cell proliferation, response of external stress, bacteria defence, and anti-inflammatory processes. Moreover, the potential performance of genomic selection (GS) in application in selective breeding for enteritis disease resistance seahorses was assessed by genomic prediction (GP). In general, the predictive accuracy of the genomic estimated breeding value (GEBV) of BayesC exceeded the rrBLUP, BayesA, RKHS, and SVM models while with no significant difference. And the GWAS-informative SNPs was significantly superior in efficience than random selected markers by comparison of predictive performance on different selection strategies of SNPs. Overall, the genetic basis of enteritis disease resistance trait in the lined seahorse is a polygenic genetic architecture. SNPs associated with the important genes of cathepsin L1-like previously reported with respect to disease resistance are consider as potential molecular markers of genetic breeding. Furthermore, GS approach is an appropriate, effective, and less-cost application in breeding enteritis disease-resistant seahorses.

Impact of freshwater rearing history on Atlantic salmon gill response to viral stimulation post seawater transfer

Land-based recirculating aquaculture systems (RAS) have risen in prevalence in recent years for Atlantic salmon production, enabling intensive production which allows increased growth and environmental control, but also having the potential for reducing water use and eutrophication. The Atlantic salmon has an anadromous life history with juvenile stages in freshwater (FW) and on-growing in seawater (SW), enabled by a transformational process known as smoltification. The timing of smoltification and transfer of smolts from FW to SW is critical under commercial production with high mortalities during this period. The impact of FW rearing system on immune function following seawater transfer (SWT) is not well understood. In this study parr were raised in either RAS or a traditional open-LOCH system until smolting and then transferred to a common marine environment. Two-weeks post-SWT fish were immune stimulated with a viral mimic (poly I:C) for 24 h to assess the ability to mount an antiviral immune response, assessed by whole transcriptome analysis of gill tissue, an important immune organ in fish. We show that unstimulated smolts reared in the LOCH had higher immune gene expression than those reared in RAS as determined by functional analysis. However, following stimulation, smolts reared in the RAS mounted a greater magnitude of response with a suite of immune genes displaying higher fold induction of transcription compared to LOCH reared smolts. We suggest RAS smolts have a lower steady state immune-associated transcriptome likely due to an unvarying environment, in terms of environmental factors and lack of exposure to pathogens, which shows a compensatory mechanism following stimulation allowing immune 'catch-up' with those reared in the LOCH. Alternatively, the RAS fish are experiencing an excessive response to the immune stimulation.

Interwoven processes in fish development: microbial community succession and immune maturation

Fishes are hosts for many microorganisms that provide them with beneficial effects on growth, immune system development, nutrition and protection against pathogens. In order to avoid spreading of infectious diseases in aquaculture, prevention includes vaccinations and routine disinfection of eggs and equipment, while curative treatments consist in the administration of antibiotics. Vaccination processes can stress the fish and require substantial farmer's investment. Additionally, disinfection and antibiotics are not specific, and while they may be effective in the short term, they have major drawbacks in the long term. Indeed, they eliminate beneficial bacteria which are useful for the host and promote the raising of antibiotic resistance in beneficial, commensal but also in pathogenic bacterial strains. Numerous publications highlight the importance that plays the diversified microbial community colonizing fish (i.e., microbiota) in the development, health and ultimately survival of their host. This review targets the current knowledge on the bidirectional communication between the microbiota and the fish immune system during fish development. It explores the extent of this mutualistic relationship: on one hand, the effect that microbes exert on the immune system ontogeny of fishes, and on the other hand, the impact of critical steps in immune system development on the microbial recruitment and succession throughout their life. We will first describe the immune system and its ontogeny and gene expression steps in the immune system development of fishes. Secondly, the plurality of the microbiotas (depending on host organism, organ, and development stage) will be reviewed. Then, a description of the constant interactions between microbiota and immune system throughout the fish's life stages will be discussed. Healthy microbiotas allow immune system maturation and modulation of inflammation, both of which contribute to immune homeostasis. Thus, immune equilibrium is closely linked to microbiota stability and to the stages of microbial community succession during the host development. We will provide examples from several fish species and describe more extensively the mechanisms occurring in zebrafish model because immune system ontogeny is much more finely described for this species, thanks to the many existing zebrafish mutants which allow more precise investigations. We will conclude on how the conceptual framework associated to the research on the immune system will benefit from considering the relations between microbiota and immune system maturation. More precisely, the development of active tolerance of the microbiota from the earliest stages of life enables the sustainable establishment of a complex healthy microbial community in the adult host. Establishing a balanced host-microbiota interaction avoids triggering deleterious inflammation, and maintains immunological and microbiological homeostasis.

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.

Genomics and transcriptomics reveal new molecular mechanism of vibriosis resistance in fish

Infectious diseases have caused dramatic production decline and economic loss for fish aquaculture. However, the poor understanding of fish disease resistance severely hampered disease prevention. Chinese tongue sole (Cynoglossus semilaevis) is an important economic flatfish suffering from vibriosis. Here we used genomic, transcriptomic and experimental approaches to investigate the molecular genetic mechanisms underlying fish vibriosis resistance. A genome-wide comparison revealed that the genes under selective sweeps were enriched for glycosaminoglycan (GAG) chondroitin sulfate (CS)/dermatan sulfate (DS) metabolism. Transcriptomic analyses prioritized synergic gene expression patterns in this pathway, which may lead to an increased CS/DS content in the resistant family. Further experimental evidence showed that carbohydrate sulfotransferases 12 (Chst12), a key enzyme for CS/DS biosynthesis, has a direct antibacterial activity. To the best of our knowledge, this is the first report that the chst12 gene has a bactericidal effect. In addition, CS/DS is a major component of the extracellular matrix (ECM) and the selection signatures and fine-tuned gene expressions of ECM-receptor interaction genes indicated a modification in the ECM structure with an enhancement of the barrier function. Furthermore, functional studies conducted on Col6a2, encoding a collagen gene which constitutes the ECM, pointed to that it may act as a cellular receptor for Vibrio pathogens, thus plays an important role for the Vibrio invasion. Taken together, these findings provide new insights into the molecular protective mechanism underlying vibriosis resistance in fish, which offers crucial genomic resources for the resistant germplasm breeding and infectious disease control in fish culturing.

Intestinal microbiome-mediated resistance against vibriosis for Cynoglossus semilaevis

BACKGROUND

Infectious diseases have caused huge economic loss and food security issues in fish aquaculture. Current management and breeding strategies heavily rely on the knowledge of regulative mechanisms underlying disease resistance. Though the intestinal microbial community was linked with disease infection, there is little knowledge about the roles of intestinal microbes in fish disease resistance. Cynoglossus semilaevis is an economically important and widely cultivated flatfish species in China. However, it suffers from outbreaks of vibriosis, which results in huge mortalities and economic loss.

RESULTS

Here, we used C. semilaevis as a research model to investigate the host-microbiome interactions in regulating vibriosis resistance. The resistance to vibriosis was reflected in intestinal microbiome on both taxonomic and functional levels. Such differences also influenced the host gene expressions in the resistant family. Moreover, the intestinal microbiome might control the host immunological homeostasis and inflammation to enhance vibriosis resistance through the microbe-intestine-immunity axis. For example, Phaeobacter regulated its hdhA gene and host cyp27a1 gene up-expressed in bile acid biosynthesis pathways, but regulated its trxA gene and host akt gene down-expressed in proinflammatory cytokines biosynthesis pathways, to reduce inflammation and resist disease infection in the resistant family. Furthermore, the combination of intestinal microbes and host genes as biomarkers could accurately differentiate resistant family from susceptible family.

CONCLUSION

Our study uncovered the regulatory patterns of the microbe-intestine-immunity axis that may contribute to vibriosis resistance in C. semilaevis. These findings could facilitate the disease control and selective breeding of superior germplasm with high disease resistance in fish aquaculture. Video Abstract.

Fish Innate Immune Response to Viral Infection-An Overview of Five Major Antiviral Genes

Fish viral diseases represent a constant threat to aquaculture production. Thus, a better understanding of the cellular mechanisms involved in establishing an antiviral state associated with protection against virus replication and pathogenesis is paramount for a sustainable aquaculture industry. This review summarizes the current state of knowledge on five selected host innate immune-related genes in response to the most relevant viral pathogens in fish farming. Viruses have been classified as ssRNA, dsRNA, and dsDNA according to their genomes, in order to shed light on what those viruses may share in common and what response may be virus-specific, both in vitro (cell culture) as well as in vivo. Special emphasis has been put on trying to identify markers of resistance to viral pathogenesis. That is, those genes more often associated with protection against viral disease, a key issue bearing in mind potential applications into the aquaculture industry.

Meta-Analysis of Immune Induced Gene Expression Changes in Diverse Drosophila melanogaster Innate Immune Responses

Simple Summary

Organisms can be infected by a wide range of pathogens, including bacteria, viruses, and parasites. Following infection, the host mounts an immune response to attempt to eliminate the pathogen. These responses are often specific to the type of pathogen and mediated by the expression of specialized genes. We have characterized the expression changes induced in host Drosophila fruit flies following infection by multiple types of pathogens, and identified a small number of genes that show expression changes in each infection. This includes genes that are known to be involved in pathogen resistance, and others that have not been previously studied as immune response genes. These findings provide new insight into transcriptional changes that accompany Drosophila immunity. They may suggest possible roles for the differentially expressed genes in innate immune responses to diverse classes of pathogens, and serve to identify candidate genes for further empirical study of these processes.

Abstract

Organisms are commonly infected by a diverse array of pathogens and mount functionally distinct responses to each of these varied immune challenges. Host immune responses are characterized by the induction of gene expression, however, the extent to which expression changes are shared among responses to distinct pathogens is largely unknown. To examine this, we performed meta-analysis of gene expression data collected from Drosophila melanogaster following infection with a wide array of pathogens. We identified 62 genes that are significantly induced by infection. While many of these infection-induced genes encode known immune response factors, we also identified 21 genes that have not been previously associated with host immunity. Examination of the upstream flanking sequences of the infection-induced genes lead to the identification of two conserved enhancer sites. These sites correspond to conserved binding sites for GATA and nuclear factor κB (NFκB) family transcription factors and are associated with higher levels of transcript induction. We further identified 31 genes with predicted functions in metabolism and organismal development that are significantly downregulated following infection by diverse pathogens. Our study identifies conserved gene expression changes in Drosophila melanogaster following infection with varied pathogens, and transcription factor families that may regulate this immune induction.

Genomic Selection in Aquaculture Species

To date, genomic prediction has been conducted in about 20 aquaculture species, with a preference for intra-family genomic selection (GS). For every trait under GS, the increase in accuracy obtained by genomic estimated breeding values instead of classical pedigree-based estimation of breeding values is very important in aquaculture species ranging from 15% to 89% for growth traits, and from 0% to 567% for disease resistance. Although the implementation of GS in aquaculture is of little additional investment in breeding programs already implementing sib testing on pedigree, the deployment of GS remains sparse, but could be boosted by adaptation of cost-effective imputation from low-density panels. Moreover, GS could help to anticipate the effect of climate change by improving sustainability-related traits such as production yield (e.g., carcass or fillet yields), feed efficiency or disease resistance, and by improving resistance to environmental variation (tolerance to temperature or salinity variation). This chapter synthesized the literature in applications of GS in finfish, crustaceans and molluscs aquaculture in the present and future breeding programs.

Differential Transcriptomic Response of Rainbow Trout to Infection with Two Strains of IPNV

The IPN virus (IPNV) causes a highly contagious disease that affects farmed salmonids. IPNV isolates have been phylogenetically classified into seven genogroups, of which two are present in Chile, genogroups 1 and 5. This study aimed to compare the transcriptomic response of rainbow trout fry challenged with two Chilean isolates of IPNV, RTTX (genogroup 1), and ALKA (genogroup 5). Tissue samples from challenged individuals and controls were taken at 1, 7, and 20 days post-challenge and analyzed by RNA-Seq. The results revealed that infection with RTTX elicited a greater modulation of the trout transcriptome compared to ALKA infection, generating a greater number of highly differentially expressed genes in relation to the control fish. Gene Ontology enrichment indicated that functions related to the inflammatory and immune responses were modulated in fish challenged with both isolates throughout the trial, but with different regulation patterns. On day 1 post challenge, these functions were activated in those challenged with ALKA, but suppressed in RTTX-challenged fish. These results suggest that rainbow trout exhibit a differential transcriptomic response to infection with the two genetically distinct IPNV isolates, especially at early times post-infection.

Interleukin 1 receptor type I (IL-1RI) is involved in the innate immune response of olive flounder (Paralichthys olivaceus) to resist pathogens

The pleiotropic cytokine IL -1 is involved in important immune responses such as thymocyte proliferation and B cell growth and differentiation. Activation of the IL -1 pathway requires its functional receptor IL -1RI, making IL -1RI the critical point of the IL -1 pathway. In-depth study of IL -1RI will help to understand the immune mechanism involved in IL -1. In this study, we identified the cDNA of the IL -1RI gene of olive flounder (PoIL-1RI). The total length of the PoIL-1RI cDNA is 2490 bp, the open reading frame is 1689 bp long and encodes a protein of 562 amino acids. The protein has three Ig domains and a typical TIR domain, as in other mammals and fish. We found that PoIL-1RI is widely expressed in the tissues studied and shows a significant immune response after stimulation with bacteria and pathogen-associated molecular patterns (PAMPs) both in vitro and in vivo. After PoIL-1RI was overexpressed in olive flounder embryonic cell line (FEC), pro-inflammatory cytokines (IL -1β, IL -6, IL -8, TNF-α) and interferon (IFN-α, IFN-γ) were significantly upregulated. And we found that after overexpressing PoIL-1RI in FEC, the antibacterial ability of FEC was significantly stronger than that of the control group, and we found that overexpression of PoIL-1RI gene significantly increased the activity of NF-κB signaling pathway. These results suggest that PoIL-1RI plays an important role in innate immune response.

The nedd-8 activating enzyme gene underlies genetic resistance to infectious pancreatic necrosis virus in Atlantic salmon

Genetic resistance to infectious pancreatic necrosis virus (IPNV) in Atlantic salmon is a rare example of a trait where a single locus (QTL) explains almost all of the genetic variation. Genetic marker tests based on this QTL on salmon chromosome 26 have been widely applied in selective breeding to markedly reduce the incidence of the disease. In the current study, whole genome sequencing and functional annotation approaches were applied to characterise genes and variants in the QTL region. This was complemented by an analysis of differential expression between salmon fry of homozygous resistant and homozygous susceptible genotypes challenged with IPNV. These analyses pointed to the NEDD-8 activating enzyme 1 (nae1) gene as a putative functional candidate underlying the QTL effect. The role of nae1 in IPN resistance was further assessed via CRISPR-Cas9 knockout of the nae1 gene and chemical inhibition of the nae1 protein activity in Atlantic salmon cell lines, both of which resulted in highly significant reduction in productive IPNV replication. In contrast, CRISPR-Cas9 knockout of a candidate gene previously purported to be a cellular receptor for the virus (cdh1) did not have a major impact on productive IPNV replication. These results suggest that nae1 is the causative gene underlying the major QTL affecting resistance to IPNV in salmon, provide further evidence for the critical role of neddylation in host-pathogen interactions, and highlight the value in combining high-throughput genomics approaches with targeted genome editing to understand the genetic basis of disease resistance.

Carps, Catla catla, Cirrhinus mrigala and Hypophthalmichthys molitrix Are Resistant to Experimental Infection with Tilapia Lake Virus (TiLV)

Tilapia tilapinevirus, also known as tilapia lake virus (TiLV), is an emerging fish virus that primarily affects tilapines. However, the virus has also been detected in a few non-tilapines. As tilapia is generally farmed in polyculture systems along with carps in South Asian countries, there is a likelihood that TiLV-infected tilapia can transmit the virus to the co-cultured species. In view of the above, the susceptibility of three carp species, namely catla (Catla catla), mrigal (Cirrhinus mrigala) and silver carp (Hypophthalmichthys molitrix) was evaluated vis-à-vis tilapia, following experimental infection with TiLV. No clinical signs and histopathological alterations could be observed in carps. RT-qPCR revealed that TiLV copy numbers in liver and brain of all the three carps were almost negligible and did not show any increase with time, suggesting that the virus did not replicate in liver and brain, the target organs of TiLV. Further, TiLV could not be isolated from pooled liver and brain tissues of carps using permissive CFF cell line. On the contrary, in tilapia, typical clinical signs and histopathological lesions were observed and there was significant increase in TiLV copy number up to 6 days post-injection. Furthermore, the virus was successfully isolated from pooled liver and brain tissue of infected tilapia. From the above findings, it could be concluded that C. catla, C. mrigala and H. molitrix are resistant to TiLV infection and unlikely to be carriers for this virus.

Innate antiviral defense demonstrates high energetic efficiency in a bony fish

Background

Viruses can impose energetic demands on organisms they infect, in part by hosts mounting resistance. Recognizing that oxygen uptake reliably indicates steady-state energy consumption in all vertebrates, we comprehensively evaluated oxygen uptake and select transcriptomic messaging in sockeye salmon challenged with either a virulent rhabdovirus (IHNV) or a low-virulent reovirus (PRV). We tested three hypotheses relating to the energetic costs of viral resistance and tolerance in this vertebrate system: (1) mounting resistance incurs a metabolic cost or limitation, (2) induction of the innate antiviral interferon system compromises homeostasis, and (3) antiviral defenses are weakened by acute stress.

Results

IHNV infections either produced mortality within 1–4 weeks or the survivors cleared infections within 1–9 weeks. Transcription of three interferon-stimulated genes (ISGs) was strongly correlated with IHNV load but not respiratory performance. Instead, early IHNV resistance was associated with a mean 19% (95% CI = 7–31%; p = 0.003) reduction in standard metabolic rate. The stress of exhaustive exercise did not increase IHNV transcript loads, but elevated host inflammatory transcriptional signaling up to sevenfold. For PRV, sockeye tolerated high-load systemic PRV blood infections. ISG transcription was transiently induced at peak PRV loads without associated morbidity, microscopic lesions, or major changes in aerobic or anaerobic respiratory performance, but some individuals with high-load blood infections experienced a transient, minor reduction in hemoglobin concentration and increased duration of excess post-exercise oxygen consumption.

Conclusions

Contrary to our first hypothesis, effective resistance against life-threatening rhabdovirus infections or tolerance to high-load reovirus infections incurred minimal metabolic costs to salmon. Even robust systemic activation of the interferon system did not levy an allostatic load sufficient to compromise host homeostasis or respiratory performance, rejecting our second hypothesis that this ancient innate vertebrate antiviral defense is itself energetically expensive. Lastly, an acute stress experienced during testing did not weaken host antiviral defenses sufficiently to promote viral replication; however, a possibility for disease intensification contingent upon underlying inflammation was indicated. These data cumulatively demonstrate that fundamental innate vertebrate defense strategies against potentially life-threatening viral exposure impose limited putative costs on concurrent aerobic or energetic demands of the organism.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01069-2.

Susceptibility of rainbow trout to three different genogroups of infectious pancreatic necrosis virus

Infectious pancreatic necrosis (IPN) is a globally distributed viral disease that is highly prevalent in rainbow trout Oncorhynchus mykiss farms in Finland. Seven genogroups (1-7) of infectious pancreatic necrosis virus (IPNV) exist, of which genogroup 5 (serotype Sp) is generally considered to be the most virulent in European salmonid farming. In Finland, 3 genogroups (2, 5 and 6) have been detected. Genogroup 2 is the most widespread and to date is the only genogroup associated with clinical disease in field observations. A bath challenge model infection trial was conducted to investigate the potential pathogenicity of the existing Finnish IPNV genogroups on IPNV-negative rainbow trout fry. Three Finnish IPNV isolates, a positive control (a Norwegian genogroup 5 isolate previously associated with high virulence in Atlantic salmon Salmo salar) and a negative control were used, and mortality was recorded daily for 8 wk. The Finnish IPNV genogroup 5 isolate caused the highest cumulative mortality, and the genogroup 2 isolate also caused elevated mortalities. The genogroup 6 isolate caused only low mortality, and the positive control treatment showed negligible mortality. Fish exposed to the Finnish genogroup 2 and 5 isolates had IPN-associated lesions, while no lesions were noted in the other treatment groups. These results indicate that Finnish IPNV genogroup 5 is potentially the most virulent IPNV genogroup for Finnish rainbow trout. Interestingly, the Norwegian IPNV genogroup 5 isolate caused only a subclinical IPN infection, providing further evidence for a host species-dependent, virus isolate-related difference in virulence in IPNV genogroup 5. The results also support the continuation of legislative disease control of IPNV genogroup 5 in Finnish inland waters.

miRNAs Predicted to Regulate Host Anti-viral Gene Pathways in IPNV-Challenged Atlantic Salmon Fry Are Affected by Viral Load, and Associated With the Major IPN Resistance QTL Genotypes in Late Infection

Infectious pancreatic necrosis virus (IPNV) infection has been a major problem in salmonid aquaculture. Marker-assisted selection of individuals with resistant genotype at the major IPN quantitative trait locus (IPN-QTL) has significantly reduced mortality in recent years. We have identified host miRNAs that respond to IPNV challenge in salmon fry that were either homozygous resistant (RR) or homozygous susceptible (SS) for the IPN-QTL. Small RNA-sequenced control samples were compared to samples collected at 1, 7, and 20 days post challenge (dpc). This revealed 72 differentially expressed miRNAs (DE miRNAs). Viral load (VL) was lower in RR vs. SS individuals at 7 and 20 dpc. However, analysis of miRNA expression changes revealed no differences between RR vs. SS individuals in controls, at 1 or 7 dpc, while 38 "high viral load responding" miRNAs (HVL-DE miRNAs) were identified at 20 dpc. Most of the HVL-DE miRNAs showed changes that were more pronounced in the high VL SS group than in the low VL RR group when compared to the controls. The absence of differences between QTL groups in controls, 1 and 7 dpc indicates that the QTL genotype does not affect miRNA expression in healthy fish or their first response to viral infections. The miRNA differences at 20 dpc were associated with the QTL genotype and could, possibly, contribute to differences in resistance/susceptibility at the later stage of infection. In silico target gene predictions revealed that 180 immune genes were putative targets, and enrichment analysis indicated that the miRNAs may regulate several major immune system pathways. Among the targets of HVL-DE miRNAs were IRF3, STAT4, NFKB2, MYD88, and IKKA. Interestingly, TNF-alpha paralogs were targeted by different DE miRNAs. Most DE miRNAs were from conserved miRNA families that respond to viral infections in teleost (e.g., miR-21, miR-146, miR-181, miR-192, miR-221, miR-462, miR-731, and miR-8159), while eight were species specific. The miRNAs showed dynamic temporal changes implying they would affect their target genes differently throughout disease progression. This shows that miRNAs are sensitive to VL and disease progression, and may act as fine-tuners of both immediate immune response activation and the later inflammatory processes.

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.

Harnessing genomics to fast-track genetic improvement in aquaculture

Aquaculture is the fastest-growing farmed food sector and will soon become the primary source of fish and shellfish for human diets. In contrast to crop and livestock production, aquaculture production is derived from numerous, exceptionally diverse species that are typically in the early stages of domestication. Genetic improvement of production traits via well-designed, managed breeding programmes has great potential to help meet the rising seafood demand driven by human population growth. Supported by continuous advances in sequencing and bioinformatics, genomics is increasingly being applied across the broad range of aquaculture species and at all stages of the domestication process to optimize selective breeding. In the future, combining genomic selection with biotechnological innovations, such as genome editing and surrogate broodstock technologies, may further expedite genetic improvement in aquaculture.

Comparison of mortality and viral load in rainbow trout (Oncorhynchus mykiss) infected with infectious pancreatic necrosis virus (IPNV) genogroups 1 and 5

Infectious pancreatic necrosis virus (IPNV) is the aetiological agent of a highly contagious disease that affects farmed salmonids. IPNV isolates have been phylogenetically classified into eight genogroups, of which two are present in Chile, genogroups 1 and 5. Here, we compare the mortality rate caused by isolates from both genogroups in rainbow trout (Oncorhynchus mykiss) fry to determine if there is an association between host susceptibility and phylogenetic characterization of IPNV. Fish were challenged by immersion with one of four isolates (two for each genogroup), and mortality curves were assessed after 30 days. Viral load was measured in all mortalities and in live fish sampled at 1, 7 and 20 days post-infection. Although mortality was low throughout the challenge, differences were found between fish infected with different isolates. Both isolates from genogroup 1 caused greater cumulative mortalities than either of the isolates from genogroup 5. When combined, the overall mortality rate of fish challenged with genogroup 1 isolates was significantly higher than those infected with genogroup 5. However, viral load was lower on trout infected with genogroup 1 isolates. These results suggest that rainbow trout are more susceptible to IPNV isolates from genogroup 1 than genogroup 5.

Variation in global transcriptomic response to Perkinsus marinus infection among eastern oyster families highlights potential mechanisms of disease resistance

Dermo disease, caused by the protozoan parasite Perkinsus marinus, negatively impacts wild and cultured Eastern oyster populations, yet our knowledge of the mechanistic bases for parasite pathogenicity and the Eastern oyster's response to it is limited. To better understand host responses to the parasite and identify molecular mechanisms underlying disease-resistance phenotypes, we experimentally challenged two families exhibiting divergent Dermo-resistance phenotypes with the parasite, generated global expression profiles using RNAseq and identified differentially expressed transcripts between control and challenged oysters from each family at multiple time points post-parasite injection. The susceptible and resistant families exhibited strikingly different transcriptomic responses to the parasite over a 28-day time period. The resistant family exhibited a strong, focused, early response to P. marinus infection, where many significantly upregulated transcripts were associated with the biological processes "regulation of proteolysis" and "oxidation-reduction process." P. marinus virulence factors are mainly comprised of proteases that facilitate parasite invasion and weaken host humoral defenses, thus host upregulation of transcripts associated with negative regulation of proteolysis is consistent with a Dermo-resistant phenotype. In contrast, the susceptible family mounted a very weak, disorganized, initial response to the parasite. Few transcripts were differentially expressed between control and injected oysters, and no functional enrichment was detected among them. At the final 28 d time point 2450 differentially expressed transcripts were identified and were associated with either "G-protein coupled receptor activity" (upregulated) or "microtubule-based process" (downregulated). A handful of protease inhibitors were differentially expressed between control and injected susceptible oysters, but this function was not enriched in the susceptible data set. The differential expression patterns observed in this study provide valuable insight into the functional basis of Dermo resistance and suggest that the timing of expression is just as important as the transcripts being expressed.

Different transcriptional response between susceptible and resistant common carp (Cyprinus carpio) fish hints on the mechanism of CyHV-3 disease resistance

Background

Infectious disease outbreaks form major setbacks to aquaculture production and to further development of this important sector. Cyprinid herpes virus-3 (CyHV-3) is a dsDNA virus widely hampering production of common carp (Cyprinus carpio), one of the most farmed fish species worldwide. Genetically disease resistant strains are highly sought after as a sustainable solution to this problem. To study the genetic basis and cellular pathways underlying disease resistance, RNA-Seq was used to characterize transcriptional responses of susceptible and resistant fish at day 4 after CyHV-3 infection.

Results

In susceptible fish, over four times more differentially expressed genes were up-regulated between day 0 and 4 compared to resistant fish. Susceptible and resistant fish responded distinctively to infection as only 55 (9%) of the up-regulated genes were shared by these two fish types. Susceptible fish elicited a typical anti-viral response, involving interferon and interferon responsive genes, earlier than resistant fish did. Furthermore, chemokine profiles indicated that the two fish types elicited different cellular immunity responses. A comparative phylogenetic approach assisted in chemokine copies annotation pointing to different orthologous copies common to bony-fishes and even carp-specific paralogs that were differentially regulated and contributed to the different response of these two fish types. Susceptible fish up-regulated more ccl19 chemokines, which attract T-cells and macrophages, the anti-viral role of which is established, whereas resistant fish up-regulated more cxcl8/il8 chemokines, which attract neutrophils, the antiviral role of which is unfamiliar.

Conclusions

Taken together, by pointing out transcriptional differences between susceptible and resistant fish in response to CyHV-3 infection, this study unraveled possible genes and pathways that take part in disease resistance mechanisms in fish and thus, enhances our understanding of fish immunogenetics and supports the development of sustainable and safe aquaculture.

Sockeye salmon demonstrate robust yet distinct transcriptomic kidney responses to rhabdovirus (IHNV) exposure and infection

Aquatic rhabdoviruses are globally significant pathogens associated with disease in both wild and cultured fish. Infectious hematopoietic necrosis virus (IHNV) is a rhabdovirus that causes the internationally regulated disease infectious hematopoietic necrosis (IHN) in most species of salmon. Yet not all naïve salmon exposed to IHNV become diseased, and the mechanisms by which some individuals evade or rapidly clear infection following exposure are poorly understood. Here we used RNA-sequencing to evaluate transcriptomic changes in sockeye salmon, a keystone species in the North Pacific and natural host for IHNV, to evaluate the consequences of IHNV exposure and/or infection on host cell transcriptional pathways. Immersion challenge of sockeye salmon smolts with IHNV resulted in approximately 33% infection prevalence, where both prevalence and viral kidney load peaked at 7 days post challenge (dpc). De novo assembly of kidney transcriptomes at 7 dpc revealed that both infected and exposed but noninfected individuals experienced substantial transcriptomic modification; however, stark variation in gene expression patterns were observed between exposed but noninfected, infected, and unexposed populations. GO and KEGG pathway enrichment in concert with differential expression analysis identified that kidney responses in exposed but noninfected fish emphasised a global pattern of transcriptional down-regulation, particularly for pathways involved in DNA transcription, protein biosynthesis and macromolecule metabolism. In contrast, transcriptomes of infected fish demonstrated a global emphasis of transcriptional up-regulation highlighting pathways involved in antiviral response, inflammation, apoptosis, and RNA processing. Quantitative PCR was subsequently used to highlight differential and time-specific regulation of acute phase, antiviral, inflammatory, cell boundary, and metabolic responsive transcripts in both infected and exposed but noninfected groups. This data demonstrates that waterborne exposure with IHNV has a dramatic effect on the sockeye salmon kidney transcriptome that is discrete between resistant and acutely susceptible individuals. We identify that metabolic, acute phase and cell boundary pathways are transcriptionally affected by IHNV and kidney responses to local infection are highly divergent from those generated as part of a disseminated response. These data suggest that primary resistance of naïve fish to IHNV may involve global responses that encourage reduced cellular signaling rather than promoting classical innate antiviral responses.

Genome-Wide Association Analysis for Resistance to Infectious Pancreatic Necrosis Virus Identifies Candidate Genes Involved in Viral Replication and Immune Response in Rainbow Trout (Oncorhynchus mykiss)

Infectious pancreatic necrosis (IPN) is a viral disease with considerable negative impact on the rainbow trout (Oncorhynchus mykiss) aquaculture industry. The aim of the present work was to detect genomic regions that explain resistance to infectious pancreatic necrosis virus (IPNV) in rainbow trout. A total of 2,278 fish from 58 full-sib families were challenged with IPNV and 768 individuals were genotyped (488 resistant and 280 susceptible), using a 57K SNP panel Axiom, Affymetrix. A genome-wide association study (GWAS) was performed using the phenotypes time to death (TD) and binary survival (BS), along with the genotypes of the challenged fish using a Bayesian model (Bayes C). Heritabilities for resistance to IPNV estimated using genomic information, were 0.53 and 0.82 for TD and BS, respectively. The Bayesian GWAS detected a SNP located on chromosome 5 explaining 19% of the genetic variance for TD. The proximity of Sentrin-specific protease 5 (SENP5) to this SNP makes it a candidate gene for resistance against IPNV. In case of BS, a SNP located on chromosome 23 was detected explaining 9% of the genetic variance. However, the moderate-low proportion of variance explained by the detected marker leads to the conclusion that the incorporation of all genomic information, through genomic selection, would be the most appropriate approach to accelerate genetic progress for the improvement of resistance against IPNV in rainbow trout.

Kinetics of local and systemic immune cell responses in whirling disease infection and resistance in rainbow trout

Background

Whirling disease (WD), caused by the myxozoan parasite Myxobolus cerebralis, is responsible for high mortalities in rainbow trout hatcheries and natural populations. To elucidate how resistant and susceptible rainbow trout strains respond to early invasion, a well-established model of WD was used to demonstrate the kinetics of local and systemic immune responses in two rainbow trout strains, the susceptible American Trout Lodge (TL) and the more resistant German Hofer strain (HO).

Methods

Parasite load and cellular immune responses were compared across several time points after M. cerebralis exposure to elucidate the kinetics of immune cells in resistant and susceptible rainbow trout in response to early invasion. In the course of the 20 days following exposure, leukocyte kinetics was monitored by flow cytometry in the caudal fin (CF), head kidney (HK) and spleen (SP). For the analysis of the leukocyte composition, cells were stained using a set of monoclonal antibodies with known specificity for distinct subpopulations of rainbow trout leukocytes.

Results

Experiments indicated general increases of CF, HK and SP myeloid cells, while decreases of B cells and T cells in the SP and HK were observed at several time points in the TL strain. On the other hand, in the HO strain, increases of T cells were dominant in CF, HK and SP at multiple time points. The differences between HO and TL were most distinct at 2, 4, 12 and 48 hours post-exposure (hpe) as well as at 4 days post-exposure (dpe), with the vast majority of innate immune response cells having higher values in the susceptible TL strain. Alteration of the leukocyte populations with augmented local cellular responses and excessive immune reactions likely lead to subsequent host tissue damage and supports parasite invasion and development in TL.

Conclusions

The findings of this study highlight the significance of effective local and systemic immune reaction and indicate proper activation of T lymphocytes critical for host resistance during M. cerebralis infection. The present study provides insights into the cellular basis of protective immune responses against M. cerebralis and can help us to elucidate the mechanisms underlying the variation in resistance to WD.

Major Histocompatibility Complex (MHC) Genes and Disease Resistance in Fish

Fascinating about classical major histocompatibility complex (MHC) molecules is their polymorphism. The present study is a review and discussion of the fish MHC situation. The basic pattern of MHC variation in fish is similar to mammals, with MHC class I versus class II, and polymorphic classical versus nonpolymorphic nonclassical. However, in many or all teleost fishes, important differences with mammalian or human MHC were observed: (1) The allelic/haplotype diversification levels of classical MHC class I tend to be much higher than in mammals and involve structural positions within but also outside the peptide binding groove; (2) Teleost fish classical MHC class I and class II loci are not linked. The present article summarizes previous studies that performed quantitative trait loci (QTL) analysis for mapping differences in teleost fish disease resistance, and discusses them from MHC point of view. Overall, those QTL studies suggest the possible importance of genomic regions including classical MHC class II and nonclassical MHC class I genes, whereas similar observations were not made for the genomic regions with the highly diversified classical MHC class I alleles. It must be concluded that despite decades of knowing MHC polymorphism in jawed vertebrate species including fish, firm conclusions (as opposed to appealing hypotheses) on the reasons for MHC polymorphism cannot be made, and that the types of polymorphism observed in fish may not be explained by disease-resistance models alone.

Antiviral defense in salmonids - Mission made possible?

Viral diseases represent one of the major threats for salmonid aquaculture. Survival from viral infections are highly dependent on host innate antiviral immune defense, where interferons are of crucial importance. Neutralizing antibodies and T cell effector mechanisms mediate long-term antiviral protection. Despite an immune cell repertoire comparable to higher vertebrates, farmed fish often fail to mount optimal antiviral protection. In the quest to multiply and spread, viruses utilize a variety of strategies to evade or escape the host immune system. Understanding the specific interplay between viruses and host immunity at depth is crucial for developing successful vaccination and treatment strategies in mammals. However, this knowledge base is still limited for pathogenic fish viruses. Here, we have focused on five RNA viruses with major impact on salmonid aquaculture: Salmonid alphavirus, Infectious salmon anemia virus, Infectious pancreatic necrosis virus, Piscine orthoreovirus and Piscine myocarditis virus. This review explore the protective immune responses that salmonids mount to these viruses and the existing knowledge on how the viruses counteract and/or bypass the immune response, including their IFN antagonizing effects and their mechanisms to establish persisting infections.

Atlantic salmon (Salmo salar L.) genetics in the 21st century: taking leaps forward in aquaculture and biological understanding

Atlantic salmon (Salmo salar L.) is among the most iconic and economically important fish species and was the first member of Salmonidae to have a high-quality reference genome assembly published. Advances in genomics have become increasingly central to the genetic improvement of farmed Atlantic salmon as well as conservation of wild salmon stocks. The salmon genome has also been pivotal in shaping our understanding of the evolutionary and functional consequences arising from an ancestral whole-genome duplication event characterising all Salmonidae members. Here, we provide a review of the current status of Atlantic salmon genetics and genomics, focussed on progress made from genome-wide research aimed at improving aquaculture production and enhancing understanding of salmonid ecology, physiology and evolution. We present our views on the future direction of salmon genomics, including the role of emerging technologies (e.g. genome editing) in elucidating genetic features that underpin functional variation in traits of commercial and evolutionary importance.

Expanding the miRNA Repertoire in Atlantic Salmon; Discovery of IsomiRs and miRNAs Highly Expressed in Different Tissues and Developmental Stages

MicroRNAs (miRNAs) are important post-transcriptional gene expression regulators. Here, 448 different miRNA genes, including 17 novel miRNAs, encoding for 589 mature Atlantic salmon miRNAs were identified after sequencing 111 samples (fry, pathogen challenged fry, various developmental and adult tissues). This increased the reference miRNAome with almost one hundred genes. Prior to isomiR characterization (mature miRNA variants), the proportion of erroneous sequence variants (ESVs) arising in the analysis pipeline was assessed. The ESVs were biased towards 5’ and 3’ end of reads in unexpectedly high proportions indicating that measurements of ESVs rather than Phred score should be used to avoid misinterpreting ESVs as isomiRs. Forty-three isomiRs were subsequently discovered. The biological effect of the isomiRs measured as increases in target diversity was small (<3%). Five miRNA genes showed allelic variation that had a large impact on target gene diversity if present in the seed. Twenty-one miRNAs were ubiquitously expressed while 31 miRNAs showed predominant expression in one or few tissues, indicating housekeeping or tissue specific functions, respectively. The miR-10 family, known to target Hox genes, were highly expressed in the developmental stages. The proportion of miR-430 family members, participating in maternal RNA clearance, was high at the earliest developmental stage.

Transcriptome analysis and histopathology of black carp (Mylopharyngodon piceus) spleen infected by Aeromonas hydrophila

Aeromonas hydrophila causes serious economic losses to the black carp (Mylopharyngodon piceus) industry. In this study, we analyzed the spleen of disease-resistant and susceptible black carp by RNA-seq. Overall, a total of 5243 terms were enriched in the gene ontology (GO) analysis, and 323 related pathways were found in the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. A total of 1935 differentially expressed genes were found and were primarily involved in cell adhesion, pathogen recognition, cellular immunity, cytokines, complement systems, and iron transport. Sixteen of the differently expressed genes involved in the immune response and the accuracy of the transcriptome data were further validated by quantitative real-time PCR (qRT-PCR). We observed Tissue sections of the spleen infected with A. hydrophila and the control group and found that the spleen of the infected group had necrosis.

Chromatin immunoprecipitation and high throughput sequencing of SVCV-infected zebrafish reveals novel epigenetic histone methylation patterns involved in antiviral immune response

Chromatin immunoprecipitation (ChIP) and high throughput sequencing (ChIP-seq) have been used to assess histone methylation (epigenetic modification) dynamics within the internal organs of zebrafish after spring viremia of carp virus (SVCV) infection. Our results show H3K4me3 up-methylation in gene promoters associated with innate immune response during the first 5 days after SVCV infection. Gene Ontology (GO) enrichment analysis confirmed up-methylation in 218 genes in the "immune system process" category. In particular, the promoters of interferon (ifn), interferon stimulated genes (isg), Toll-like receptors (tlr) and c-reactive protein (crp) multi gene sets were marked with the permissive H3K4 methylation. Higher histone 3 methylation was associated with higher transcription levels of the corresponding genes. Therefore, the evidence presented here suggests that transcriptional regulation at the promoter level of key immune genes of the interferon signaling pathway and c-reactive proteins genes can be modulated by epigenetic modification of histones. This study emphasizes the importance of epigenetic control in the response of zebrafish to SVCV infection.

Genomic Tools and Selective Breeding in Molluscs

The production of most farmed molluscs, including mussels, oysters, scallops, abalone, and clams, is heavily dependent on natural seed from the plankton. Closing the lifecycle of species in hatcheries can secure independence from wild stocks and enables long-term genetic improvement of broodstock through selective breeding. Genomic techniques have the potential to revolutionize hatchery-based selective breeding by improving our understanding of the characteristics of mollusc genetics that can pose a challenge for intensive aquaculture and by providing a new suite of tools for genetic improvement. Here we review characteristics of the life history and genetics of molluscs including high fecundity, self-fertilization, high genetic diversity, genetic load, high incidence of deleterious mutations and segregation distortion, and critically assess their impact on the design and effectiveness of selective breeding strategies. A survey of the results of current breeding programs in the literature show that selective breeding with inbreeding control is likely the best strategy for genetic improvement of most molluscs, and on average growth rate can be improved by 10% per generation and disease resistance by 15% per generation across the major farmed species by implementing individual or family-based selection. Rapid advances in sequencing technology have resulted in a wealth of genomic resources for key species with the potential to greatly improve hatchery-based selective breeding of molluscs. In this review, we catalog the range of genomic resources currently available for molluscs of aquaculture interest and discuss the bottlenecks, including lack of high-quality reference genomes and the relatively high cost of genotyping, as well as opportunities for applying genomics-based selection.

Transcriptomic Profiles of Senegalese Sole Infected With Nervous Necrosis Virus Reassortants Presenting Different Degree of Virulence

Betanodaviruses [nervous necrosis virus (NNV)] are the causative agent of the viral encephalopathy and retinopathy, a disease that affects cultured Senegalese sole (Solea senegalensis). NNV reassortants, combining genomic segments from redspotted grouper nervous necrosis virus (RGNNV) and striped jack nervous necrosis virus (SJNNV) genotypes, have been previously isolated from several fish species. The wild-type reassortant wSs160.03, isolated from Senegalese sole, has been proven to be more virulent to sole than the parental genotypes (RGNNV and SJNNV), causing 100% mortality. Mutations at amino acids 247 (serine to alanine) and 270 (serine to asparagine) in the wSs160.03 capsid protein have allowed us to obtain a mutant reassortant (rSs160.03_247+270), which provokes a 40% mortality decrease. In this study, the RNA-Seq technology has been used to comparatively analyze Senegalese sole transcriptomes in two organs (head kidney and eye/brain) after infection with wild-type and mutant strains. A total of 633 genes were differentially expressed (DEGs) in animals infected with the wild-type isolate (with higher virulence), whereas 393 genes were differentially expressed in animals infected with the mutant strain (37.9% decrease in the number of DEGs). To study the biological functions of detected DEGs involved in NNV infection, a gene ontology (GO) enrichment analysis was performed. Different GO profiles were obtained in the following subclasses: (i) biological process; (ii) cellular component; and (iii) molecular function, for each viral strain tested. Immune response and proteolysis have been the predominant biological process after the infection with the wild-type isolate, whereas the infection with the mutant strain induces proteolysis in head kidney and inhibition of vasculogenesis in nervous tissue. Regarding the immune response, genes coding for proteins acting as mediators of type I IFN expression (DHX58, IRF3, IRF7) and IFN-stimulated genes (ISG15, Mx, PKR, Gig1, ISG12, IFI44, IFIT-1, to name a few) were upregulated in animals infected with the wild-type isolate, whereas no-differential expression of these genes was observed in samples inoculated with the mutant strain. The different transcriptomic profiles obtained could help to better understand the NNV pathogenesis in Senegalese sole, setting up the importance as virulence determinants of amino acids at positions 247 and 270 within the RNA2 segment.

Genetic and transcriptomic analyses provide new insights on the early antiviral response to VHSV in resistant and susceptible rainbow trout

BACKGROUND

The viral hemorrhagic septicemia virus (VHSV) is a major threat for salmonid farming and for wild fish populations worldwide. Previous studies have highlighted the importance of innate factors regulated by a major quantitative trait locus (QTL) for the natural resistance to waterborne VHSV infection in rainbow trout. The aim of this study was to analyze the early transcriptomic response to VHSV inoculation in cell lines derived from previously described resistant and susceptible homozygous isogenic lines of rainbow trout to obtain insights into the molecular mechanisms responsible for the resistance to the viral infection.

RESULTS

We first confirmed the presence of the major QTL in a backcross involving a highly resistant fish isogenic line (B57) and a highly susceptible one (A22), and were able to define the confidence interval of the QTL and to identify its precise position. We extended the definition of the QTL since it controls not only resistance to waterborne infection but also the kinetics of mortality after intra-peritoneal injection. Deep sequencing of the transcriptome of B57 and A22 derived cell lines exposed to inactivated VHSV showed a stronger response to virus inoculation in the resistant background. In line with our previous observations, an early and strong induction of interferon and interferon-stimulated genes was correlated with the resistance to VHSV, highlighting the major role of innate immune factors in natural trout resistance to the virus. Interestingly, major factors of the antiviral innate immunity were much more expressed in naive B57 cells compared to naive A22 cells, which likely contributes to the ability of B57 to mount a fast antiviral response after viral infection. These observations were further extended by the identification of several innate immune-related genes localized close to the QTL area on the rainbow trout genome.

CONCLUSIONS

Taken together, our results improve our knowledge in virus-host interactions in vertebrates and provide novel insights in the molecular mechanisms explaining the resistance to VHSV in rainbow trout. Our data also provide a collection of potential markers for resistance and susceptibility of rainbow trout to VHSV infection.

Iron Overload Is Associated With Oxidative Stress and Nutritional Immunity During Viral Infection in Fish

Iron is a trace element, essential to support life due to its inherent ability to exchange electrons with a variety of molecules. The use of iron as a cofactor in basic metabolic pathways is essential to both pathogenic microorganisms and their hosts. During evolution, the shared requirement of micro- and macro-organisms for this important nutrient has shaped the pathogen-host relationship. Infectious pancreatic necrosis virus (IPNv) affects salmonids constituting a sanitary problem for this industry as it has an important impact on post-smolt survival. While immune modulation induced by IPNv infection has been widely characterized on Salmo salar, viral impact on iron host metabolism has not yet been elucidated. In the present work, we evaluate short-term effect of IPNv on several infected tissues from Salmo salar. We observed that IPNv displayed high tropism to headkidney, which directly correlates with a rise in oxidative stress and antiviral responses. Transcriptional profiling on headkidney showed a massive modulation of gene expression, from which biological pathways involved with iron metabolism were remarkable. Our findings suggest that IPNv infection increase oxidative stress on headkidney as a consequence of iron overload induced by a massive upregulation of genes involved in iron metabolism.

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.

Long-term affected flat oyster (Ostrea edulis) haemocytes show differential gene expression profiles from naïve oysters in response to Bonamia ostreae

European flat oyster (Ostrea edulis) production has suffered a severe decline due to bonamiosis. The responsible parasite enters in oyster haemocytes, causing an acute inflammatory response frequently leading to death. We used an immune-enriched oligo-microarray to understand the haemocyte response to Bonamia ostreae by comparing expression profiles between naïve (NS) and long-term affected (AS) populations along a time series (1 d, 30 d, 90 d). AS showed a much higher response just after challenge, which might be indicative of selection for resistance. No regulated genes were detected at 30 d in both populations while a notable reactivation was observed at 90 d, suggesting parasite latency during infection. Genes related to extracellular matrix and protease inhibitors, up-regulated in AS, and those related to histones, down-regulated in NS, might play an important role along the infection. Twenty-four candidate genes related to resistance should be further validated for selection programs aimed to control bonamiosis.

Issues with RNA-seq analysis in non-model organisms: A salmonid example

High throughput sequencing (HTS) is useful for many purposes as exemplified by the other topics included in this special issue. The purpose of this paper is to look into the unique challenges of using this technology in non-model organisms where resources such as genomes, functional genome annotations or genome complexity provide obstacles not met in model organisms. To describe these challenges, we narrow our scope to RNA sequencing used to study differential gene expression in response to pathogen challenge. As a demonstration species we chose Atlantic salmon, which has a sequenced genome with poor annotation and an added complexity due to many duplicated genes. We find that our RNA-seq analysis pipeline deciphers between duplicates despite high sequence identity. However, annotation issues provide problems in linking differentially expressed genes to pathways. Also, comparing results between approaches and species are complicated due to lack of standardized annotation.

Distinct herpesvirus resistances and immune responses of three gynogenetic clones of gibel carp revealed by comprehensive transcriptomes

Background

Gibel carp is an important aquaculture species in China, and a herpesvirus, called as Carassius auratus herpesvirus (CaHV), has hampered the aquaculture development. Diverse gynogenetic clones of gibel carp have been identified or created, and some of them have been used as aquaculture varieties, but their resistances to herpesvirus and the underlying mechanism remain unknown.

Results

To reveal their susceptibility differences, we firstly performed herpesvirus challenge experiments in three gynogenetic clones of gibel carp, including the leading variety clone A⁺, candidate variety clone F and wild clone H. Three clones showed distinct resistances to CaHV. Moreover, 8772, 8679 and 10,982 differentially expressed unigenes (DEUs) were identified from comparative transcriptomes between diseased individuals and control individuals of clone A⁺, F and H, respectively. Comprehensive analysis of the shared DEUs in all three clones displayed common defense pathways to the herpesvirus infection, activating IFN system and suppressing complements. KEGG pathway analysis of specifically changed DEUs in respective clones revealed distinct immune responses to the herpesvirus infection. The DEU numbers identified from clone H in KEGG immune-related pathways, such as “chemokine signaling pathway”, “Toll-like receptor signaling pathway” and others, were remarkably much more than those from clone A⁺ and F. Several IFN-related genes, including Mx1, viperin, PKR and others, showed higher increases in the resistant clone H than that in the others. IFNphi3, IFI44-like and Gig2 displayed the highest expression in clone F and IRF1 uniquely increased in susceptible clone A⁺. In contrast to strong immune defense in resistant clone H, susceptible clone A⁺ showed remarkable up-regulation of genes related to apoptosis or death, indicating that clone A⁺ failed to resist virus offensive and evidently induced apoptosis or death.

Conclusions

Our study is the first attempt to screen distinct resistances and immune responses of three gynogenetic gibel carp clones to herpesvirus infection by comprehensive transcriptomes. These differential DEUs, immune-related pathways and IFN system genes identified from susceptible and resistant clones will be beneficial to marker-assisted selection (MAS) breeding or molecular module-based resistance breeding in gibel carp.

Electronic supplementary material

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

Differences in responses of grass carp to different types of grass carp reovirus (GCRV) and the mechanism of hemorrhage revealed by transcriptome sequencing

Background

Grass carp is an important farmed fish in China that is affected by serious disease, especially hemorrhagic disease caused by grass carp reovirus (GCRV). The mechanism underlying the hemorrhagic symptoms in infected fish remains to be elucidated. Although GCRV can be divided into three distinct subtypes, differences in the pathogenesis and host immune responses to the different subtypes are still unclear. The aim of this study was to provide a comprehensive insight into the grass carp response to different GCRV subtypes and to elucidate the mechanism underlying the hemorrhagic symptoms.

Results

Following infection of grass carp, GCRV-I was associated with a long latent period and low mortality (42.5%), while GCRV-II was associated with a short latent period and high mortality (81.4%). The relative copy number of GCRV-I remained consistent or decreased slightly throughout the first 7 days post-infection, whereas a marked increase in GCRV-II high copy number was detected at 5 days post-infection. Transcriptome sequencing revealed 211 differentially expressed genes (DEGs) in Group I (66 up-regulated, 145 down-regulated) and 670 (386 up-regulated, 284 down-regulated) in Group II. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed significant enrichment in the terms or pathways involved in immune responses and correlating with blood or platelets. Most of the DEGs in Group I were also present in Group II, although the expression profiles differed, with most DEGs showing mild changes in Group I, while marked changes were observed in Group II, especially the interferon-related genes. Many of the genes involved in the complement pathway and coagulation cascades were significantly up-regulated at 7 days post-infection in Group II, suggesting activation of these pathways.

Conclusion

GCRV-I is associated with low virulence and a long latent period prior to the induction of a mild host immune response, whereas GCRV-II is associated with high virulence, a short latent period and stimulates a strong and extensive host immune response. The complement and coagulation pathways are significantly activated at 7 days post-infection, leading to the endothelial cell and blood cell damage that result in hemorrhagic symptoms.

Electronic supplementary material

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

The dietary replacement of marine ingredients by terrestrial animal and plant alternatives modulates the antiviral immune response of Atlantic salmon (Salmo salar)

The effects of replacing marine ingredients by terrestrial ingredients on the health of Atlantic salmon (Salmo salar) are poorly understood. During a 14-week trial, Atlantic salmon fed a fish meal-fish oil based diet (MAR) showed similar growth performance to others fed a plant protein/vegetable oil based diet (VEG), whereas poorer performance was observed in those fed an animal by-product meal/vegetable oil based diet (ABP). At the end of the trial, salmon were injected with either phosphate-buffered saline (PBS) or the viral mimic polyriboinosinic polyribocytidylic acid (pIC) and sampled for head kidney RNA after 24 h. The levels of 27 immune-related transcripts, and of 5 others involved in eicosanoid synthesis (including paralogues in both cases) were measured in the head kidney of the salmon using qPCR. All of the assayed immune-related genes and cox2 were pIC-induced, while the other eicosanoid synthesis-related genes were pIC-repressed. Linear regression was used to establish correlations between different immune transcripts, elucidating the cascade of responses to pIC and specialization among paralogues. Regarding the effect of diet on the antiviral immune response, pIC-treated fish fed diets ABP and VEG showed higher transcript levels of tlr3, irf1b, stat1a, isg15b, and gig1 compared to those fed diet MAR. We infer that the observed dietary immunomodulation could be due to the lower proportion of arachidonic acid (ARA), eicosapentaenoic acid (EPA), and docosahexaenoic acid (DHA) in diets ABP and VEG. Furthermore, our results suggest a major role of dietary ARA in Atlantic salmon immunity, as low ARA proportion in diet VEG coincided with the highest pIC-induction of some immune transcripts (tlr7, stat1c, mxb, and gig1) and the lowest levels of transcripts encoding eicosanoid-synthesizing enzymes (5loxa, 5loxb, and pgds). In contrast, the high ARA/EPA ratio of diet ABP appeared to favor increased expression of transcripts involved in the synthesis of pro-inflammatory eicosanoids (5loxa and 5loxb) and chemotaxis (ccl19b). In conclusion, our findings show that nutritionally balanced plant-based diets may enhance the immune response of Atlantic salmon. Future studies should explore the possible advantages of plant-based diets in Atlantic salmon exposed to a viral infection.