High-throughput transcriptome analysis of ISAV-infected Atlantic salmon Salmo salar unravels divergent immune responses associated to head-kidney, liver and gills tissues

Persistent immune responses in the heart determine the outcome of cardiomyopathy syndrome in Atlantic salmon (Salmo salar)

Cardiomyopathy syndrome (CMS) caused by piscine myocarditis virus (PMCV) is a severe cardiac disease in Atlantic salmon (Salmo salar) and one of the leading causes of morbidity and mortality in the Norwegian aquaculture industry. Previous research suggest a variation in individual susceptibility to develop severe disease, however the role of the immune response in determining individual outcome of CMS is poorly understood particularly in cases where fish are also challenged by stress. The present study's aim was therefore to characterize cardiac transcriptional responses to PMCV infection in Atlantic salmon responding to infection under stressful conditions with a high versus low degree of histopathological damage. The study was performed as a large-scale controlled experiment of Atlantic salmon smolts from pre-challenge to 12 weeks post infection (wpi) with PMCV, during which fish were exposed to intermittent stressors. RNA sequencing (RNAseq) was used to compare the heart transcriptome of high responders (HR) with atrium histopathology score '4' and low responders (LR) with score '0.5' at 12 wpi. A high-throughput quantitative PCR (qPCR) analysis was used to compare immune gene transcription between individuals sampled at 6, 9 and 12 wpi. Based on RNAseq and qPCR results, RNAscope in situ hybridization (ISH) was performed for visualization of IFN-γ - and IFNb producing immune cells in affected heart tissue. Compared to LR, the transcription of 1592 genes was increased in HR at 12 wpi. Of these genes, around. 40 % were immune-related, including various chemokines, key antiviral response molecules, and genes. associated with a Th1 pro-inflammatory immune response. Further, the qPCR analysis confirmed. increased immune gene transcription in HR at both 9 and 12 wpi, despite a decrease in PMCV. transcription between these time points. Interestingly, increased IFNb transcription in HR suggests the. presence of high-quantity IFN secreting cells in the hearts of these individuals. Indeed, RNAscope. confirmed the presence of IFN-γ and IFNb-positive cells in the heart ventricle of HR but not LR. To conclude, our data indicate that in severe outcomes of PMCV infection various chemokines attract leucocytes to the salmon heart, including IFN-γ and IFNb-secreting cells, and that these cells play important roles in maintaining persistent antiviral responses and a sustained host immunopathology despite decreasing heart viral transcription.

Comparative analysis of the immune responses of CcIgZ3 in mucosal tissues and the co-expression of CcIgZ3 and PCNA in the gills of common carp (Cyprinus carpio L.) in response to TNP-LPS

Fish live in an aquatic environment rich in various microorganisms and pathogens. Fish mucosal-associated lymphoid tissue (MALT) plays a very important role in immune defence. This study was conducted to characterize the immune response mediated by CcIgZ3 in common carp (Cyprinus carpio.) and investigate the proliferating CcIgZ3+ B lymphocytes in gill. We determined the expression of CcIgZ3 in many different tissues of common carp following stimulation by intraperitoneal injection of TNP-LPS (2,4,6-Trinitrophenyl hapten conjugated to lipopolysaccharide) or TNP-KLH (2,4,6-Trinitrophenyl hapten conjugated to Keyhole Limpet Hemocyanin). Compared with TNP-KLH, TNP-LPS can induce greater CcIgZ3 expression in the head kidney, gill and hindgut, especially in the gill. The results indicate that the gill is one of the main sites involved in the immune response mediated by CcIgZ3. To examine the distribution of CcIgZ3+ B lymphocytes, immunohistochemistry (IHC) experiments were performed using a polyclonal antibody against CcIgZ3. The results indicated that CcIgZ3 was detected in the head kidney, hindgut and gill. To further examine whether CcIgZ3+ B lymphocytes proliferate in the gills, proliferating CcIgZ3+ B cells were analysed by immunofluorescence staining using an anti-CcIgZ3 polyclonal antibody and an anti-PCNA monoclonal antibody. CcIgZ3 and PCNA (Proliferating Cell Nuclear Antigen) double-labelled cells in the gills were located within the epithelial cells of the gill filaments of common carp stimulated with TNP-LPS at 3 dps and 7 dps, and relatively more proliferating CcIgZ3+ B cells appeared in the gills of common carp at 7 dps. These data imply that CcIgZ3+ B cells in the gills might be produced by local proliferation following TNP-LPS stimulation. In summary, compared with those in TNP-KLH, CcIgZ3 preferentially affects the gills of common carp following challenge with TNP-LPS. CcIgZ3+ B cells proliferate in the gills to quickly produce the CcIgZ3 antibody. In addition, CcIgZ3+ B cells can be activated to induce a strong immune response very early locally in the gill and produce the antibody CcIgZ3, which helps exert an immune-protective effect. These results suggest that an effective vaccine can be designed to promote production of the mucosal antibody CcIgZ3.

Understanding host response to infectious salmon anaemia virus in an Atlantic salmon cell line using single-cell RNA sequencing

BACKGROUND

Infectious Salmon Anaemia Virus (ISAV) is an Orthomixovirus that represents a large problem for salmonid aquaculture worldwide. Current prevention and treatment methods are only partially effective. Genetic selection and genome engineering have the potential to develop ISAV resistant salmon stocks. Both strategies can benefit from an improved understanding of the genomic regulation of ISAV pathogenesis. Here, we used single-cell RNA sequencing of an Atlantic salmon cell line to provide the first high dimensional insight into the transcriptional landscape that underpins host-virus interaction during early ISAV infection.

RESULTS

Salmon head kidney (SHK-1) cells were single-cell RNA sequenced at 24, 48 and 96 h post-ISAV challenge. At 24 h post infection, cells showed expression signatures consistent with viral entry, with genes such as PI3K, FAK or JNK being upregulated relative to uninfected cells. At 48 and 96 h, infected cells showed a clear anti-viral response, characterised by the expression of IFNA2 or IRF2. Uninfected bystander cells at 48 and 96 h also showed clear transcriptional differences, potentially suggesting paracrine signalling from infected cells. These bystander cells expressed pathways such as mRNA sensing, RNA degradation, ubiquitination or proteasome; and up-regulation of mitochondrial ribosome genes also seemed to play a role in the host response to the infection. Correlation between viral and host genes revealed novel genes potentially key for this fish-virus interaction.

CONCLUSIONS

This study has increased our understanding of the cellular response of Atlantic salmon during ISAV infection and revealed host-virus interactions at the cellular level. Our results highlight various potential key genes in this host-virus interaction, which can be manipulated in future functional studies to increase the resistance of Atlantic salmon to ISAV.

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.

CRISPR/Cas9-Mediated Gene Editing in Salmonids Cells and Efficient Establishment of Edited Clonal Cell Lines

Finfish production has seen over three-fold increase in the past 30 years (1990-2020), and Atlantic salmon (A. salmon; salmo salar) accounted for approximately 32.6% of the total marine and coastal aquaculture of all finfish species in the year 2020, making it one of the most profitable farmed fish species globally. This growth in production is, however, threatened by a number of problems which can be solved using the CRISPR/Cas technology. In vitro applications of CRISPR/Cas using cell lines can complement its in vivo applications, but salmonids-derived cell lines are difficult to gene edit because they grow slowly, are difficult to transfect and isolate single clones of gene-edited cells. While clonal isolation of the gene-edited Chinook salmon cell line (CHSE-214) has successfully been performed, there is no report of successful clonal isolation of the gene-edited A. salmon ASK-1 and SHK-1cell lines. In the current study, two gene loci-cr2 and mmp9 of A. salmon-were efficiently edited using the ribonucleoprotein (RNP) and plasmid CRISPR/Cas9 strategies. Edited cells were enriched using flow cytometer-activated cell sorting (FACS), followed by clonal isolation and expansion of edited cells. The study both confirms the recent report of the highly efficient editing of these widely used model cell lines, as well as extends the frontline in the single-cell cloning of gene-edited salmonids cells. The report also highlights the pitfalls and future directions in the application of CRISPR/Cas9 in these cells.

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.

Gill and Liver Transcript Expression Changes Associated With Gill Damage in Atlantic Salmon (Salmo salar)

Gill damage represents a significant challenge in the teleost fish aquaculture industry globally, due to the gill’s involvement in several vital functions and direct contact with the surrounding environment. To examine the local and systemic effects accompanying gill damage (which is likely to negatively affect gill function) of Atlantic salmon, we performed a field sampling to collect gill and liver tissue after several environmental insults (e.g., harmful algal blooms). Before sampling, gills were visually inspected and gill damage was scored; gill scores were assigned from pristine [gill score 0 (GS0)] to severely damaged gills (GS3). Using a 44K salmonid microarray platform, we aimed to compare the transcriptomes of pristine and moderately damaged (i.e., GS2) gill tissue. Rank Products analysis (5% percentage of false-positives) identified 254 and 34 upregulated and downregulated probes, respectively, in GS2 compared with GS0. Differentially expressed probes represented genes associated with functions including gill remodeling, wound healing, and stress and immune responses. We performed gill and liver qPCR for all four gill damage scores using microarray-identified and other damage-associated biomarker genes. Transcripts related to wound healing (e.g., neb and klhl41b) were significantly upregulated in GS2 compared with GS0 in the gills. Also, transcripts associated with immune and stress-relevant pathways were dysregulated (e.g., downregulation of snaclec 1-like and upregulation of igkv3) in GS2 compared with GS0 gills. The livers of salmon with moderate gill damage (i.e., GS2) showed significant upregulation of transcripts related to wound healing (i.e., chtop), apoptosis (e.g., bnip3l), blood coagulation (e.g., f2 and serpind1b), transcription regulation (i.e., pparg), and stress-responses (e.g., cyp3a27) compared with livers of GS0 fish. We performed principal component analysis (PCA) using transcript levels for gill and liver separately. The gill PCA showed that PC1 significantly separated GS2 from all other gill scores. The genes contributing most to this separation were pgam2, des, neb, tnnt2, and myom1. The liver PCA showed that PC1 significantly separated GS2 from GS0; levels of hsp70, cyp3a27, pparg, chtop, and serpind1b were the highest contributors to this separation. Also, hepatic acute phase biomarkers (e.g., serpind1b and f2) were positively correlated to each other and to gill damage. Gill damage-responsive biomarker genes and associated qPCR assays arising from this study will be valuable in future research aimed at developing therapeutic diets to improve farmed salmon welfare.

Transcriptome Analysis of Atlantic Salmon (Salmo salar) Skin in Response to Sea Lice and Infectious Salmon Anemia Virus Co-Infection Under Different Experimental Functional Diets

Sea lice (Lepeophtheirus salmonis) are ectoparasitic copepods that cause significant economic loss in marine salmoniculture. In commercial salmon farms, infestation with sea lice can enhance susceptibility to other significant pathogens, such as the highly contagious infectious salmon anemia virus (ISAv). In this study, transcriptomic analysis was used to evaluate the impact of four experimental functional feeds (i.e. 0.3% EPA/DHA+high-ω6, 0.3% EPA/DHA+high-ω6+immunostimulant (IS), 1% EPA/DHA+high-ω6, and 1% EPA/DHA+high-ω3) on Atlantic salmon (Salmo salar) during a single infection with sea lice (L. salmonis) and a co-infection with sea lice and ISAv. The overall objectives were to compare the transcriptomic profiles of skin between lice infection alone with co-infection groups and assess differences in gene expression response among animals with different experimental diets. Atlantic salmon smolts were challenged with L. salmonis following a 28-day feeding trial. Fish were then challenged with ISAv at 18 days post-sea lice infection (dpi), and maintained on individual diets, to establish a co-infection model. Skin tissues sampled at 33 dpi were subjected to RNA-seq analysis. The co-infection’s overall survival rates were between 37%-50%, while no mortality was observed in the single infection with lice. With regard to the infection status, 756 and 1303 consensus differentially expressed genes (DEGs) among the four diets were identified in “lice infection vs. pre-infection” and “co-infection vs. pre-infection” groups, respectively, that were shared between the four experimental diets. The co-infection groups (co-infection vs. pre-infection) included up-regulated genes associated with glycolysis, the interferon pathway, complement cascade activity, and heat shock protein family, while the down-regulated genes were related to antigen presentation and processing, T-cell activation, collagen formation, and extracellular matrix. Pathway enrichment analysis conducted between infected groups (lice infection vs. co-infection) resulted in several immune-related significant GO terms and pathways unique to this group, such as “autophagosome”, “cytosolic DNA-sensing pathway” and “response to type I interferons”. Understanding how experimental functional feeds can impact the host response and the trajectory of co-infections will be an essential step in identifying efficacious intervention strategies that account for the complexities of disease in open cage culture.

Chitosan nanoparticle formulation attenuates poly (I:C) induced innate immune responses against inactivated virus vaccine in Atlantic salmon (Salmo salar)

Many vaccine formulations, in particular vaccines based on inactivated virus, needs adjuvants to boost immunogenicity. In aquaculture, mineral and plant oil are used as adjuvant in commercial vaccines, and the advent of oil-adjuvanted vaccines was crucial to aquaculture development. Nevertheless, some of these approved vaccines display suboptimal performance in the field compared to experimental conditions. Therefore, there is a need to improve adjuvants and delivery methods for fish vaccines against viruses. We used RNA sequencing of Atlantic salmon head kidney to analyse the difference in gene expression 24 h after injection of different experimental vaccine formulations. We compared five different formulations in addition to a PBS control: inactivated virus alone (group V), soluble poly (I:C) (group P), nanoparticles containing poly (I:C) (group N), soluble poly (I:C) + inactivated virus (group PV) and finally nanoparticles containing poly (I:C) + inactivated virus (group NV). Our results showed poly (I:C)'s ability as adjuvant and its capacity influence innate immune genes expression in Atlantic salmon. Soluble poly (I:C) upregulated multiple immune related genes and was more effective compared to poly (I:C) formulated into chitosan nanoparticles (more than 10 fold increase in differentially expressed genes, DEGs). However, inclusion of inactivated ISA virus in the nanoparticle vaccine, increased the number of DEGs fivefold suggesting a synergistic effect of adjuvant and antigen. Our results indicate that the way poly (I:C) is formulated and the presence of antigen is important for the magnitude of the innate immune response in Atlantic salmon.

Expression of ssa-miR-155 during ISAV infection in vitro: Putative role as a modulator of the immune response in Salmo salar

Multiple cellular components are involved in pathogen-host interaction during viral infection; in this context, the role of miRNAs have become highly relevant. We assessed the expression of selected miRNAs during an in vitro infection of a Salmo salar cell line with Infectious Salmon Anemia Virus (ISAV), the causative agent of a severe disease by the same name. Salmon orthologs for miRNAs that regulate antiviral responses were measured using RT-qPCR in an in vitro time-course assay. We observed a modulation of specific miRNAs expression, where ssa-miR-155-5p was differentially over-expressed. Using in silico analysis, we identified the putative mRNA targets for ssa-miR-155-5p, finding a high prevalence of hosts immune response-related genes; moreover, several mRNAs involved in the viral infective process were also identified as targets for this miRNA. Our results suggest a relevant role for miR-155-5p in Salmo salar during an ISAV infection as a regulator of the immune response to the virus.

Exploring genetic resistance to infectious salmon anaemia virus in Atlantic salmon by genome-wide association and RNA sequencing

Background

Infectious Salmonid Anaemia Virus (ISAV) causes a notifiable disease that poses a large threat for Atlantic salmon (Salmo salar) aquaculture worldwide. There is no fully effective treatment or vaccine, and therefore selective breeding to increase resistance to ISAV is a promising avenue for disease prevention. Genomic selection and potentially genome editing can be applied to enhance host resistance, and these approaches benefit from improved knowledge of the genetic and functional basis of the target trait. The aim of this study was to characterise the genetic architecture of resistance to ISAV in a commercial Atlantic salmon population and study its underlying functional genomic basis using RNA Sequencing.

Results

A total of 2833 Atlantic salmon parr belonging to 194 families were exposed to ISAV in a cohabitation challenge in which cumulative mortality reached 63% over 55 days. A total of 1353 animals were genotyped using a 55 K SNP array, and the estimate of heritability for the trait of binary survival was 0.13–0.33 (pedigree-genomic). A genome-wide association analysis confirmed that resistance to ISAV was a polygenic trait, albeit a genomic region in chromosome Ssa13 was significantly associated with resistance and explained 3% of the genetic variance. RNA sequencing of the heart of 16 infected (7 and 14 days post infection) and 8 control fish highlighted 4927 and 2437 differentially expressed genes at 7 and 14 days post infection respectively. The complement and coagulation pathway was down-regulated in infected fish, while several metabolic pathways were up-regulated. The interferon pathway showed little evidence of up-regulation at 7 days post infection but was mildly activated at 14 days, suggesting a potential crosstalk between host and virus. Comparison of the transcriptomic response of fish with high and low breeding values for resistance highlighted TRIM25 as being up-regulated in resistant fish.

Conclusions

ISAV resistance shows moderate heritability with a polygenic architecture, but a significant QTL was detected on chromosome 13. A mild up-regulation of the interferon pathway characterises the response to the virus in heart samples from this population of Atlantic salmon, and candidate genes showing differential expression between samples with high and low breeding values for resistance were identified.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07671-6.

Transcriptome of Edwardsiella anguillarum in vivo and in vitro revealed two-component system, ABC transporter and flagellar assembly are three pathways pathogenic to European eel (Anguilla anguilla)

Edwardsiella anguillarum is one of the common bacterial pathogens for the cultivated eels in China. The aim of this study was to reveal the cause of E. anguillarum pathogenic to European eel (Anguilla anguilla) from the perspective of the transcriptome. In this study, we first prepared E. anguillarum cultured in vitro and analysed the whole transcriptome after extracting the total RNA. Then, eels were i.p injected with E. anguillarum, and total RNA were extracted from the liver of European eels 48 h after the infection. After sequencing the transcriptome, we obtained average 1.97 × 10⁸ clean reads cultured in vitro and 1.36 × 10⁵ clean reads located in vivo after annotating all reads into the genome of E. anguillarum. The whole transcriptome showed, compared to the E. anguillarum cultured in vitro, 503 significantly up and 657 significantly down-regulated different expressed genes (DEGs) were observed. KEGG analysis showed that 38 DEGs of Two-Component System, 41 DEGs of ABC transporter, and 10 DEGs flagellar assembly pathways were highly upregulated in E. anguillarum located in vivo. Then, we designed primers to analyse the up-regulated DEGs through qRT-PCR and confirmed some up-regulated DEGs. The results of this study provide important reference for the further study of pathogen-host interaction between E. anguillarum and European eel.

The validation of the role of several genes related to Bombyx mori nucleopolyhedrovirus infection in vivo

Bombyx mori nucleopolyhedrovirus (BmNPV) is one of primary silkworm pathogens and causes a serious damage of cocoon losses every year. Recent years, many works have been done to clarify the silkworm anti-BmNPV mechanism, and a significant progress has been made in screening and studying of genes and proteins related to BmNPV infection, but several of them lacked the proofs in vivo. In this study, to further validate the function of seven newly reported genes in vivo, including BmAtlatin-n, Bmferritin-heavy chain (BmFerHCH), Bmthymosin (BmTHY), Bmseroin1, Bmseroin2, Bmnuclear hormone receptors 96 (BmNHR96), and BmE3 ubiquitin-protein ligase SINA-like 10 (BmSINAL10), the response of them in the midgut, fat body, and hemolymph of differentially resistant strains (resistant strain YeA and susceptible strain YeB) at 48 h following BmNPV infection were analyzed. The results showed that the relative stable or upregulated expression level of BmAtlatin-n, BmTHY, Bmseroin1, and Bmseroin2 in YeA resistant strain following BmNPV infection further indicated their antiviral role in vivo, compared with susceptible YeB strain. Moreover, the significant downregulation of BmFerHCH, BmNHR96, and BmSINAL10 in both strains following BmNPV infection revealed their role in benefiting virus infection, as well as the upregulation of BmFerHCH in YeB midgut and BmSINAL10 in YeB hemolymph. These data could be used to complementary the proofs of the function of these genes in response to BmNPV infection.

Transcriptome Response of Atlantic Salmon (Salmo salar) to a New Piscine Orthomyxovirus

Pilchard orthomyxovirus (POMV) is an emerging pathogen of concern to the salmon industry in Australia. To explore the molecular events that underpin POMV infection, we challenged Atlantic salmon (Salmo salar) post-smolts in seawater via cohabitation. Tissue samples of the head kidney and liver were collected from moribund and surviving individuals and analyzed using transcriptome sequencing. Viral loads were higher in the head kidney compared to the liver, yet the liver presented more upregulated genes. Fish infected with POMV showed a strong innate immune response that included the upregulation of pathogen recognition receptors such as RIG-I and Toll-like receptors as well as the induction of interferon-stimulated genes (MX, ISG15). Moribund fish also presented a dramatic induction of pro-inflammatory cytokines, contributing to severe tissue damage and morbidity. An induction of major histocompatibility complex (MHC) class I genes (B2M) and markers of T cell-mediated immunity (CD8-alpha, CD8-beta, Perforin-1, Granzyme-A) was observed in both moribund fish and survivors. In addition, differential connectivity analysis showed that three key regulators (RELA/p65, PRDM1, and HLF) related to cell-mediated immunity had significant differences in connectivity in "clinically healthy" versus "clinically affected" or moribund fish. Collectively, our results show that T cell-mediated immunity plays a central role in the response of Atlantic salmon to the infection with POMV.

Kinetics of transcriptional response against poly (I:C) and infectious salmon anemia virus (ISAV) in Atlantic salmon kidney (ASK) cell line

Vaccine adjuvants induce host innate immune responses improving long-lasting adaptive immunity against vaccine antigens. In vitro models can be used to compare these responses between adjuvants and the infection targeted by the vaccine. We utilized transcriptomic profiling of an Atlantic salmon cell line to compare innate immune responses against ISAV and an experimental viral vaccine adjuvant: poly (I:C). Induction of interferon and interferon induced genes were observed after both treatments, but often with different amplitude and kinetics. Using KEGG ortholog database and available software from Bioconductor we could specify a complete bioinformatic pipeline for analysis of transcriptomic data from Atlantic salmon, a feature not previously available. We have identified important differences in the transcriptional profile of Atlantic salmon cells exposed to viral infection and a viral vaccine adjuvant candidate, poly (I:C). This report increases our knowledge of viral host-pathogen interaction in salmon and to which extent these can be mimicked by adjuvant compounds.

Polarized Trout Epithelial Cells Regulate Transepithelial Electrical Resistance, Gene Expression, and the Phosphoproteome in Response to Viral Infection

The burden of disease is a major challenge in aquaculture production. The fish gill characterized with a large surface area and short route to the bloodstream is a major environmental interface and a significant portal of entry for pathogens. To investigate gill responses to viral infection the salmonid gill cell line RTgill-W1 was stimulated with synthetic dsRNA and the salmonid alphavirus subtype 2 (SAV-2). Epithelial integrity in polarized cells can be measured as transepithelial electrical resistance (TEER) which is defined as the electrical resistance across a cell monolayer. TEER is a widely accepted quantitative measure of cellular integrity of a cell monolayer. TEER increased immediately after stimulation with the synthetic dsRNA, polyinosinic:polycytidylic acid (poly(I:C)). In parallel, tight junction and gene expression of innate immune activation markers was modulated in response to poly(I:C). The SAV-2 virus was found to replicate at a low level in RTgill-W1 cells where TEER was disturbed at an early stage of infection, however, gene expression related to tight junction regulation was not modulated. A strong poly(I:C)-driven antiviral response was observed including increases of Rig-like receptors (RLRs) and interferon stimulating genes (ISGs) mRNAs. At the level of signal transduction, poly(I:C) stimulation was accompanied by the phosphorylation of 671 proteins, of which 390 were activated solely in response to the presence of poly(I:C). According to motif analysis, kinases in this group included MAPKs, Ca2+/calmodulin-dependent kinase (CaMK) and cAMP-dependent protein kinase (PKA), all reported to be activated in response to viral infection in mammals. Results also highlighted an activation of the cytoskeletal organization that could be mediated by members of the integrin family. While further work is needed to validate these results, our data indicate that salmonid gill epithelia has the ability to mount a significant response to viral infection which might be important in disease progression. In vitro cell culture can facilitate both a deeper understanding of the anti-viral response in fish and open novel therapeutic avenues for fish health management in aquaculture.

Integration of Transcriptome, Gross Morphology and Histopathology in the Gill of Sea Farmed Atlantic Salmon (Salmo salar): Lessons From Multi-Site Sampling

The gill of teleost fish is a multifunctional organ involved in many physiological processes such as gas exchange, osmotic and ionic regulation, acid-base balance and excretion of nitrogenous waste. Due to its extensive interface with the environment, the gill plays a key role as a primary mucosal defense tissue against pathogens, as manifested by the presence of the gill-associated lymphoid tissue (GIALT). In recent years, the prevalence of multifactorial gill pathologies has increased significantly, causing substantial losses in Atlantic salmon aquaculture. The transition from healthy to unhealthy gill phenotypes and the progression of multifactorial gill pathologies, such as proliferative gill disease (PGD), proliferative gill inflammation (PGI) and complex gill disorder (CGD), are commonly characterized by epithelial hyperplasia, lamellar fusion and inflammation. Routine monitoring for PGD relies on visual inspection and non-invasive scoring of the gill tissue (gross morphology), coupled with histopathological examination of gill sections. To explore the underlying molecular events that are associated with the progression of PGD, we sampled Atlantic salmon from three different marine production sites in Scotland and examined the gill tissue at three different levels of organization: gross morphology with the use of PGD scores (macroscopic examination), whole transcriptome (gene expression by RNA-seq) and histopathology (microscopic examination). Our results strongly suggested that the changes in PGD scores of the gill tissue were not associated with the changes in gene expression or histopathology. In contrast, integration of the gill RNA-seq data with the gill histopathology enabled us to identify common gene expression patterns associated with multifactorial gill disease, independently from the origin of samples. We demonstrated that the gene expression patterns associated with multifactorial gill disease were dominated by two processes: a range of immune responses driven by pro-inflammatory cytokines and the events associated with tissue damage and repair, driven by caspases and angiogenin.

Transcriptome analysis of grass carp (Ctenopharyngodon idella) between fast- and slow-growing fish

Grass carp is one of the most important freshwater aquaculture species in China. However, the mechanisms underlying the growth of muscle tissue in the fish are unclear. High-throughput RNA-Seq was used to analyze the transcriptome of grass carp muscle tissue between fast- and slow-growing fish family groups. Twenty-four individuals each from 4 fast-growing families and 4 slow-growing families were used to reduce background noise. 71 up-regulated and 35 down-regulated genes were identified in the differentially expressed genes (DEGs). GO and KEGG enrichment analyses revealed the DEGs were involved in the GH/IGF axis, calcium metabolism, protein and glycogen synthesis, oxygen transport, cytoskeletal and myofibrillar components. IGFBP1 was up-regulated in big fish while GHR2 was down-regulated. Glutamic pyruvate transaminase 2, an indicator of liver tissue damage, was down-regulated in big grass carp, which indicates that the fish was better adapted to an artificially formulated diet. GAPDH, the rate-limiting enzyme in glycolytic flux was highly expressed in fast-growing grass carp, reflecting enhanced carbohydrate metabolism. Higher expression of ALAS2 and myoglobin 1 in big grass carp, related to oxygen transport might promote aerobic exercise along with food intake and muscle growth. Genes for cytoskeletal and myofibrillar components such as tropomyosin, meromyosin, and troponin I were also up-regulated in big grass carp. These results provide valuable information about the key genes for use as biomarkers of growth in selective breeding programs for grass carp and contribute to our understanding of the molecular mechanisms and regulative pathways regulating growth in fish.

The specificity protein 3 (SP3) gene in ducks (Anas platyrhynchos): cloning, characterization and expression during viral infection

Specificity Protein 3 (SP3) is a newly identified regulator of tumor growth and invasiveness in humans. In this study, we identified and characterized the function of duck SP3 (duSP3). The full-length cDNA sequence of the duSP3 gene was cloned via rapid amplification of cDNA ends. It contained 2468 nucleotides, including a 111 base pair (bp) 5'-untranslated region (UTR), 215 bp 3'-UTR, and 2142 bp open reading frame (ORF), which encoded a 713 amino acid (AA) strongly conserved with Avian SP3. Tissue specificity analysis demonstrated that duSP3 was constitutively expressed in the eight tissues tested: liver, spleen, lung, heart, kidney, thymus, breast, and leg; and low expression levels were observed in all tissues, except the spleen and thymus. Quantitative real-time polymerase chain reaction (qRT-PCR) analysis revealed that duSP3 expression rapidly increased in vitro after stimulation with both the hepatitis virus (DHV-1) and polyriboinosinic polyribocytidylic acid (poly(I:C)). However, the expression under these treatments varied in kidney and liver tissues; in the liver, duSP3 increased significantly at 36 h after the DHV-1 treatment and peaked at 72 h after poly(I:C) stimulation. These results suggested that SP3 may play a positive role in immune responses against viral infections in ducks.

Microbial Danger Signals Control Transcriptional Induction of Distinct MHC Class I L Lineage Genes in Atlantic Salmon

Antigen processing and presentation by major histocompatibility complex (MHC) molecules is a cornerstone in vertebrate immunity. Like mammals, teleosts possess both classical MHC class I and multiple families of divergent MHC class I genes. However, while certain mammalian MHC class I-like molecules have proven to be integral in immune regulation against a broad array of pathogens, the biological relevance of the different MHC class I lineages in fish remains elusive. This work focuses on MHC class I L lineage genes and reveals unique regulatory patterns of six genes (Sasa-lia, Sasa-lda, Sasa-lca, Sasa-lga, Sasa-lha, and Sasa-lfa) in antimicrobial immunity of Atlantic salmon (Salmo salar L.). Using two separate in vivo challenge models with different kinetics and immune pathologies combined with in vitro stimulation using viral and bacterial TLR ligands, we show that de novo synthesis of different L lineage genes is distinctly regulated in response to various microbial stimuli. Prior to the onset of classical MHC class I gene expression, lia was rapidly and systemically induced in vivo by the single-stranded (ss) RNA virus salmonid alpha virus 3 (SAV3) but not in response to the intracellular bacterium Piscirickettsia salmonis. In contrast, lga expression was upregulated in response to both viral and bacterial stimuli. A role for distinct MHC class I L-lineage genes in anti-microbial immunity in salmon was further substantiated by a marked upregulation of lia and lga gene expression in response to type I IFNa stimulation in vitro. Comparably, lha showed no transcriptional induction in response to IFNa stimulation but was strongly induced in response to a variety of viral and bacterial TLR ligands. In sharp contrast, lda showed no response to viral or bacterial challenge. Similarly, induction of lca, which is predominantly expressed in primary and secondary lymphoid tissues, was marginal with the exception of a strong and transient upregulation in pancreas following SAV3 challenge Together, these findings suggest that certain Atlantic salmon MHC class I L lineage genes play important and divergent roles in early anti-microbial response and that their regulation, in response to different activation signals, represents a system for selectively promoting the expression of distinct non-classical MHC class I genes in response to different types of immune challenges.

The important role of phagocytosis and interleukins for nile tilapia (Oreochromis niloticus) to defense infection of Aeromonas hydrophila based on transcriptome analysis

Tilapia is an important economic fish worldwide. It is vital to understand the mechanism of immune response for the prevention and treatment the infection of Aeromonas hydrophila. Based on high-throughput sequencing of Illumina HiSeq™, we found differentially expressed genes in the immune-related pathway were classified into phagosome, cytokine-cytokine receptor interaction and toll-like receptor signaling pathway. Gene Ontology terms were divided into three categories of transporting function, DNA replication activity and energy supply activity. The first one was related to phagocytosis and the process or transporting of antigen driven by tubulins; the second one was to differentiation and proliferation of lymphocyte activated by cytokines; and the former two both needed energy provided by the third one. According to colchicine assay, cross-immune assay, ELISA of interleukins and classical phagocytosis assay, phagocytosis and interleukins were verified to be most important to defense the infection of A. hydrophila.

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.

Sumoylation of nucleoprotein (NP) mediated by activation of NADPH oxidase complex is a consequence of oxidative cellular stress during infection by Infectious salmon anemia (ISA) virus necessary to viral progeny

The study evaluated the effects of nucleoprotein viral and the infectious virus in SHK-1 cells. The results show a strong respiratory burst activation and the induction of p47phox, SOD, GLURED, and apoptotic genes. Additionally, the cells alter the profile of SUMOylated proteins by the effect of transfection and infection experiments. In silico analyses show a set of structural motifs in NP susceptible of post-translational modification by the SUMO protein. Interestingly, the inhibition of the NADPH oxidase complex blocked the production of reactive oxygen species and the high level of cellular ROS due to the nucleoprotein and the ISAv. At the same time, the blocking of the p38MAPK signaling pathway and the use of Aristotelia chilensis, decreased viral progeny production. These results suggest that the NP triggers a strong production of ROS and modifying the post-translational profile mediated by SUMO-2/3, a phenomenon that favors the production of new virions.

Plasma proteomic analysis of zebrafish following spring viremia of carp virus infection

To better understand spring viremia of carp virus (SVCV) pathogenesis in zebrafish proteomic analysis was used to examine the plasma protein profile in SVCV-infected zebrafish. A total of 3062 proteins were identified. Of those 137, 63 and 31 proteins were enriched in blood samples harvested at 1, 2 and 5 days post SVCV infection, respectively. These altered host proteins were classified based on their biological function: 23 proteins under the response to stimulus term were identified. Interestingly, at the top of the up-regulated proteins during SVCV infection were the proteins of the vitellogenin family (Vtg) and the grass carp reovirus-induced gene (Gig) proteins. Real-time RT-PCR evaluation of samples from internal organs verified that SVCV infection induced vtg and gig2 gene expression already at day 1 post-infection. Western blot analysis revealed the presence of Vtg protein only in blood of SVCV-infected fish. This is the first proteomic study to reveal the involvement of Vtg proteins in adult fish response to viral challenge. It also highlights the role of Gig proteins as important factors in antiviral response in fish. This work provides valuable relevant insight into virus-host interaction and the identification of molecular markers of fish response to virus.

Integrating Genomic and Morphological Approaches in Fish Pathology Research: The Case of Turbot (Scophthalmus maximus) Enteromyxosis

Enteromyxosis, caused by Enteromyxum scophthalmi, is one of the most devastating diseases stemming from myxozoan parasites in turbot (Scophthalmus maximus L.), being a limiting factor for its production. The disease develops as a cachectic syndrome, associated to catarrhal enteritis and leukocytic depletion, with morbidity and mortality rates usually reaching 100%. To date, no effective treatment exists and there are different unknown issues concerning its pathogenesis. The gross and microscopic lesions associated to enteromyxosis have been thoroughly described, and several morphopathological studies have been carried out to elucidate the mechanisms of this host-parasite interaction. More recently, efforts have been focused on a multidisciplinary approach, combining histopathology and transcriptome analysis, which has provided significant advances in the understanding of the pathogenesis of this parasitosis. RNA-Seq technology was applied at early and advanced stages of the disease on fishes histologically evaluated and classified based on their lesional degree. In the same way, the transcriptomic data were analyzed in relation to the morphopathological picture and the course of the disease. In this paper, a comprehensive review of turbot enteromyxosis is presented, starting from the disease description up to the most novel information extracted by an integrated approach on the infection mechanisms and host response. Further, we discuss ongoing strategies toward a full understanding of host-pathogen interaction and the identification of suitable biomarkers for early diagnosis and disease management strategies.

The Change of Teleost Skin Commensal Microbiota Is Associated With Skin Mucosal Transcriptomic Responses During Parasitic Infection by Ichthyophthirius multifillis

Teleost skin serves as the first line of defense against invading pathogens, and contain a skin-associated lymphoid tissue (SALT) that elicit gut-like immune responses against antigen stimulation. Moreover, exposed to the water environment and the pathogens therein, teleost skin is also known to be colonized by diverse microbial communities. However, little is known about the interactions between microbiota and the teleost skin mucosal immune system, especially dynamic changes about the interactions under pathogen infection. We hypothesized that dramatic changes of microbial communities and strong mucosal immune response would be present in the skin of aquatic vertebrate under parasite infection. To confirm this hypothesis, we construct an infected model with rainbow trout (Oncorhynchus mykiss), which was experimentally challenged by Ichthyophthirius multifiliis (Ich). H & E staining of trout skin indicates the successful invasion of Ich and shows the morphological changes caused by Ich infection. Critically, increased mRNA expression levels of immune-related genes were detected in trout skin from experimental groups using qRT-PCR, which were further studied by RNA-Seq analysis. Here, through transcriptomics, we detected that complement factors, pro-inflammatory cytokines, and antimicrobial genes were strikingly induced in the skin of infected fish. Moreover, high alpha diversity values of microbiota in trout skin from the experimental groups were discovered. Interestingly, we found that Ich infection led to a decreased abundance of skin commensals and increased colonization of opportunistic bacteria through 16S rRNA pyrosequencing, which were mainly characterized by lose of Proteobacteria and increased intensity of Flavobacteriaceae. To our knowledge, our results suggest for the first time that parasitic infection could inhibit symbionts and offer opportunities for other pathogens' secondary infection in teleost skin.

Transcriptome sequencing of the gill and barbel of Southern catfish (Silurus meridionalis) revealed immune responses and novel rhamnose-binding lectins (RBLs)

Southern catfish (Silurus meridionalis) is an economically important species widely cultured in China. It is well known for its fast growth, strong resistance to diseases and euryphage. However, little is known about the mechanisms for its powerful immune systems. Our Fish-T1K project has finished its first phase of 200 fish transcriptomes, with sequencing of gills in most examined fishes. In this study, we performed transcriptome sequencing of the gill and the maxillary barbel of Southern catfish, with the latter as a control. High expression of immune-related transcripts were observed in these two tissues. We observed that genes in the T cell receptor signaling pathway had higher transcription values in the gill than in the barbel. In addition, eight new rhamnose-binding lectins (RBLs) were identified and their carbohydrate recognition domains (CRDs) were classified according to the eight conserved cysteine residues and two conserved motifs (-YGR- and -DPC-). This is the first transcriptome report by high-throughput sequencing of the Southern catfish. Our genomic data and discovery of novel RBLs in this project should be able to promote better understandings of the roles of gills in immune responses and disease prevention for further aquaculture.

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.

Combined transcriptomic/proteomic analysis of crucian carp Carassius auratus gibelio in cyprinid herpesvirus 2 infection

Cyprinid herpesvirus 2 (CyHV-2) is a pathogen of herpesviral hematopoietic necrosis disease of crucian carp. Our study aimed to investigate the molecular mechanisms and immune response at the mRNA and protein levels in head kidney during CyHV-2 infection. Three days after infection with CyHV-2, 7085 differentially expressed genes were identified by transcriptome sequencing, of which 3090 were up-regulated and 3995 were down-regulated. And 338 differentially expressed proteins including 277 up-regulated and 61 down-regulated were identified using tandem mass tag labeling followed by liquid chromatography tandem mass spectrometry. Notably, 128 differentially co-expressed genes at mRNA and protein levels (cDEGs) were reliably quantified, including 86 co-up-regulated and 42 co-down-regulated. In addition, 10 cDGEs in the above pathways were selected for qRT-PCR to confirm the validity of the transcriptome and proteome changes by showing that RIG-I, MDA5, LGP2, FAS, PKR and PKZ up-regulated and Integrin α, Integrin β2, NCF2 and NCF4 down-regulated. This indicated that after CyHV-2 infection, the herpes simplex infection pathway, RIG-I like receptor signaling pathway, necroptosis pathway and p53 signaling pathway were activated and the phagosome pathway was suppressed. Our findings reveal the pathogenesis and the host immune mechanism of CyHV-2 infection of crucian carp.

In vivo virulence and genomic comparison of infectious Salmon Anaemia Virus isolates from Atlantic Canada

The infectious salmon anaemia virus (ISAV) is capable of causing a significant disease in Atlantic salmon, which has resulted in considerable financial losses for salmon farmers around the world. Since the first detection of ISAV in Canada in 1996, it has been a high priority for aquatic animal health management and surveillance programmes have led to the identification of many genetically distinct ISAV isolates of variable virulence. In this study, we evaluated the virulence of three ISAV isolates detected in Atlantic Canada in 2012 by doing in vivo-controlled disease challenges with two sources of Atlantic salmon. We measured viral loads in fish tissues during the course of infection. Sequences of the full viral RNA genomes of these three ISAV isolates were obtained and compared to a high-virulence and previously characterized isolate detected in the Bay of Fundy in 2004, as well as a newly identified ISAV NA-HPR0 isolate. All three ISAV isolates studied were shown to be of low to mid-virulence with fish from source A having a lower mortality rate than fish from source B. Viral load estimation using an RT-qPCR assay targeting viral segment 8 showed a high degree of similarity between tissues. Through genomic comparison, we identified various amino acid substitutions unique to some isolates, including a stop codon in the segment 8 ORF2 not previously reported in ISAV, present in the isolate with the lowest observed virulence.

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.

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.

Teleosts Genomics: Progress and Prospects in Disease Prevention and Control

Genome wide studies based on conventional molecular tools and upcoming omics technologies are beginning to gain functional applications in the control and prevention of diseases in teleosts fish. Herein, we provide insights into current progress and prospects in the use genomics studies for the control and prevention of fish diseases. Metagenomics has emerged to be an important tool used to identify emerging infectious diseases for the timely design of rational disease control strategies, determining microbial compositions in different aquatic environments used for fish farming and the use of host microbiota to monitor the health status of fish. Expounding the use of antimicrobial peptides (AMPs) as therapeutic agents against different pathogens as well as elucidating their role in tissue regeneration is another vital aspect of genomics studies that had taken precedent in recent years. In vaccine development, prospects made include the identification of highly immunogenic proteins for use in recombinant vaccine designs as well as identifying gene signatures that correlate with protective immunity for use as benchmarks in optimizing vaccine efficacy. Progress in quantitative trait loci (QTL) mapping is beginning to yield considerable success in identifying resistant traits against some of the highly infectious diseases that have previously ravaged the aquaculture industry. Altogether, the synopsis put forth shows that genomics studies are beginning to yield positive contribution in the prevention and control of fish diseases in aquaculture.

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.

Transcriptome Analysis Based on RNA-Seq in Understanding Pathogenic Mechanisms of Diseases and the Immune System of Fish: A Comprehensive Review

In recent years, with the advent of next-generation sequencing along with the development of various bioinformatics tools, RNA sequencing (RNA-Seq)-based transcriptome analysis has become much more affordable in the field of biological research. This technique has even opened up avenues to explore the transcriptome of non-model organisms for which a reference genome is not available. This has made fish health researchers march towards this technology to understand pathogenic processes and immune reactions in fish during the event of infection. Recent studies using this technology have altered and updated the previous understanding of many diseases in fish. RNA-Seq has been employed in the understanding of fish pathogens like bacteria, virus, parasites, and oomycetes. Also, it has been helpful in unraveling the immune mechanisms in fish. Additionally, RNA-Seq technology has made its way for future works, such as genetic linkage mapping, quantitative trait analysis, disease-resistant strain or broodstock selection, and the development of effective vaccines and therapies. Until now, there are no reviews that comprehensively summarize the studies which made use of RNA-Seq to explore the mechanisms of infection of pathogens and the defense strategies of fish hosts. This review aims to summarize the contemporary understanding and findings with regard to infectious pathogens and the immune system of fish that have been achieved through RNA-Seq technology.

Omics in fish mucosal immunity

The mucosal immune system of fish is a complex network of immune cells and molecules that are constantly surveilling the environment and protecting the host from infection. A number of "omics" tools are now available and utilized to understand the complexity of mucosal immune systems in non-traditional animal models. This review summarizes recent advances in the implementation of "omics" tools pertaining to the four mucosa-associated lymphoid tissues in teleosts. Genomics, transcriptomics, proteomics, and "omics" in microbiome research require interdisciplinary collaboration and careful experimental design. The data-rich datasets generated are proving really useful at discovering new innate immune players in fish mucosal secretions, identifying novel markers of specific mucosal immune responses, unraveling the diversity of the B and T cell repertoires and characterizing the diversity of the microbial communities present in teleost mucosal surfaces. Bioinformatics, data analysis and storage platforms should be developed to facilitate rapid processing of large datasets, especially when mammalian tools such as bioinformatics analysis software are not available in fishes.

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.

In the shadow: The emerging role of long non-coding RNAs in the immune response of Atlantic salmon

The genomic era has increased the research effort to uncover how the genome of an organism, and specifically the transcriptome, is modulated after interplaying with pathogenic microorganisms and ectoparasites. However, the ever-increasing accessibility of sequencing technology has also evidenced regulatory roles of long non-coding RNAs (lncRNAs) related to several biological processes including immune response. This study reports a high-confidence annotation and a comparative transcriptome analysis of lncRNAs from several tissues of Salmo salar infected with the most prevalent pathogens in the Chilean salmon aquaculture such as the infectious salmon anemia (ISA) virus, the intracellular bacterium Piscirickettsia salmonis and the ectoparasite copepod Caligus rogercresseyi. Our analyses showed that lncRNAs are widely modulated during infection. However, this modulation is pathogen-specific and highly correlated with immuno-related genes associated with innate immune response. These findings represent the first discovery for the widespread differential expression of lncRNAs in response to infections with different types of pathogens in Atlantic salmon, suggesting that lncRNAs are pivotal player during the fish immune response.

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.

Transcriptomic analysis of spleen infected with infectious salmon anemia virus reveals distinct pattern of viral replication on resistant and susceptible Atlantic salmon (Salmo salar)

The infectious salmon anemia virus (ISAv) produces a systemic infection in salmonids, causing large losses in salmon production. However, little is known regarding the mechanisms exerting disease resistance. In this paper, we perform an RNA-seq analysis in Atlantic salmon challenged with ISAv (using individuals coming from families that were highly susceptible or highly resistant to ISAv infection). We evaluated the differential expression of both host and ISAv genes in a target organ for the virus, i.e. the spleen. The results showed differential expression of host genes related to response to stress, immune response and protein folding (genes such as; atf3, mhc, mx1-3, cd276, cd2, cocs1, c7, il10, il10rb, il13ra2, ubl-1, ifng, ifngr1, hivep2, sigle14 and sigle5). An increased protein processing activity was found in susceptible fish, which generates a subsequent unfolded protein response. We observed extreme differences in the expression of viral segments between susceptible and resistant groups, demonstrating the capacity of resistant fish to overcome the virus replication, generating a very low viral load. This phenomenon and survival of this higher resistant fish seem to be related to differences in immune and translational process, as well as to the increase of HIV-EP2 (hivep2) transcript in resistant fish, although the causal mechanism is yet to be discovered. This study provides valuable information about disease resistance mechanisms in Atlantic salmon from a host-pathogen interaction point of view.

Novel insights into the response of Atlantic salmon (Salmo salar) to Piscirickettsia salmonis: Interplay of coding genes and lncRNAs during bacterial infection

Despite the high prevalence and impact to Chilean salmon aquaculture of the intracellular bacterium Piscirickettsia salmonis, the molecular underpinnings of host-pathogen interactions remain unclear. Herein, the interplay of coding and non-coding transcripts has been proposed as a key mechanism involved in immune response. Therefore, the aim of this study was to evidence how coding and non-coding transcripts are modulated during the infection process of Atlantic salmon with P. salmonis. For this, RNA-seq was conducted in brain, spleen, and head kidney samples, revealing different transcriptional profiles according to bacterial load. Additionally, while most of the regulated genes annotated for diverse biological processes during infection, a common response associated with clathrin-mediated endocytosis and iron homeostasis was present in all tissues. Interestingly, while endocytosis-promoting factors and clathrin inductions were upregulated, endocytic receptors were mainly downregulated. Furthermore, the regulation of genes related to iron homeostasis suggested an intracellular accumulation of iron, a process in which heme biosynthesis/degradation pathways might play an important role. Regarding the non-coding response, 918 putative long non-coding RNAs were identified, where 425 were newly characterized for S. salar. Finally, co-localization and co-expression analyses revealed a strong correlation between the modulations of long non-coding RNAs and genes associated with endocytosis and iron homeostasis. These results represent the first comprehensive study of putative interplaying mechanisms of coding and non-coding RNAs during bacterial infection in salmonids.

Evaluation of differentially expressed immune-related genes in intestine of Pelodiscus sinensis after intragastric challenge with lipopolysaccharide based on transcriptome analysis

Pelodiscus sinensis is the most common turtle species that has been raised in East and Southeast Asia. However, there are still limited studies about the immune defense mechanisms in its small intestine until now. In the present research, histological analysis and transcriptome analysis was performed on the small intestine of P. sinensis after intragastric challenge with LPS to explore its mechanisms of immune responses to pathogens. The result showed the number of intraepithelial lymphocytes (IELs) and goblet cells (GCs) in its intestine increased significantly at 48 h post-challenge with LPS by intragastrical route, indicating clearly the intestinal immune response was induced. Compared with the control, a total of 748 differentially expressed genes (DEGs) were identified, including 361 up-regulated genes and 387 down-regulated genes. Based on the Gene Ontology (GO) annotation and the Kyoto Encyclopedia of Genes and Genomes (KEGG), 48 immune-related DEGs were identified, which were classified into 82 GO terms and 14 pathways. Finally, 18 DEGs, which were randomly selected, were confirmed by quantitative real-time PCR (qRT-PCR). Our results provide valuable information for further analysis of the immune defense mechanisms against pathogens in the small intestine of P. sinensis.

Comparative Transcriptome Analysis of Bombyx mori (Lepidoptera) Larval Midgut Response to BmNPV in Susceptible and Near-Isogenic Resistant Strains

Bombyx mori nucleopolyhedrovirus (BmNPV) is one of the primary pathogens causing severe economic losses in sericulture. However, the molecular mechanism of silkworm resistance to BmNPV remains largely unknown. Here, the recurrent parent P50 (susceptible strain) and the near-isogenic line BC9 (resistance strain) were used in a comparative transcriptome study examining the response to infection with BmNPV. A total of 14,300 unigenes were obtained from two different resistant strains; of these, 869 differentially expressed genes (DEGs) were identified after comparing the four transcriptomes. Many DEGs associated with protein metabolism, cytoskeleton, and apoptosis may be involved in the host response to BmNPV infection. Moreover, some immunity related genes were also altered following BmNPV infection. Specifically, after removing genetic background and individual immune stress response genes, 22 genes were found to be potentially involved in repressing BmNPV infection. These genes were related to transport, virus replication, intracellular innate immune, and apoptosis. Our study provided an overview of the molecular mechanism of silkworm resistance to BmNPV infection and laid a foundation for controlling BmNPV in the future.

Long noncoding RNAs (lncRNAs) dynamics evidence immunomodulation during ISAV-Infected Atlantic salmon (Salmo salar)

Despite evidence for participation in the host response to infection, the roles of many long non-coding RNAs (lncRNAs) remain unknown. Therefore, the aims of this study were to identify lncRNAs in Atlantic salmon (Salmo salar) and evaluate their transcriptomic regulation during ISA virus (ISAV) infection, an Orthomyxoviridae virus associated with high mortalities in salmonid aquaculture. Using next-generation sequencing, whole-transcriptome analysis of the Salmo salar response to ISAV infection was performed, identifying 5,636 putative lncRNAs with a mean length of 695 base pairs. The transcriptional modulation evidenced a similar number of differentially expressed lncRNAs in the gills (3,294), head-kidney (3,275), and liver (3,325) over the course of the infection. Moreover, analysis of a subset of these lncRNAs showed the following: (i) Most were similarly regulated in response to ISA virus infection; (ii) The transcript subsets were uniquely modulated in each tissue (gills, liver, and head-kidney); and (iii) A subset of lncRNAs were upregulated for each tissue and time analysed, indicating potential markers for ISAV infection. These findings represent the first discovery of widespread differential expression of lncRNAs in response to virus infection in non-model species, suggesting that lncRNAs could be involved in regulating the host response during ISAV infection.

Mucosal immunoglobulins at respiratory surfaces mark an ancient association that predates the emergence of tetrapods

Gas-exchange structures are critical for acquiring oxygen, but they also represent portals for pathogen entry. Local mucosal immunoglobulin responses against pathogens in specialized respiratory organs have only been described in tetrapods. Since fish gills are considered a mucosal surface, we hypothesized that a dedicated mucosal immunoglobulin response would be generated within its mucosa on microbial exposure. Supporting this hypothesis, here we demonstrate that following pathogen exposure, IgT⁺ B cells proliferate and generate pathogen-specific IgT within the gills of fish, thus providing the first example of locally induced immunoglobulin in the mucosa of a cold-blooded species. Moreover, we demonstrate that gill microbiota is predominantly coated with IgT, thus providing previously unappreciated evidence that the microbiota present at a respiratory surface of a vertebrate is recognized by a mucosal immunoglobulin. Our findings indicate that respiratory surfaces and mucosal immunoglobulins are part of an ancient association that predates the emergence of tetrapods.

In teleost fish the gills perform—in addition to respiration—functions such as immune defence. Here the authors show that IgT, a teleost specific Ig previously shown to be involved in gut and skin mucosal immunity, is locally induced in the gill, where it plays a key role in immunity in rainbow trout.

Transcriptomic analysis of koi (Cyprinus carpio) spleen tissue upon cyprinid herpesvirus 3 (CyHV3) infection using next generation sequencing

Cyprinid Herpesvirus 3 (CyHV-3) can infect and specifically cause a huge economic loss in both common carp (Cyprinus carpio) and its ornamental koi variety. The molecular mechanisms underlying CyHV-3 infection are not well understood. In this study, koi spleen tissues of both mock and CyHV-3 infection groups were collected, and high-throughput sequencing technology was used to analyze the differentially expressed genes (DEGs) at the transcriptome level. A total of 105,356,188 clean reads from two libraries were obtained. After the de novo assembly of the transcripts, 129,314 unigenes were generated. Of these unigenes, 70,655 unigenes were matched to the known proteins in the database, while 2190 unigenes were predicted by ESTScan software. Comparing the infection group to the mock group, a total of 23,029 significantly differentially expressed unigenes were identified, including 10,493 up-regulated DEGs and 12,536 down-regulated DEGs. GO (Gene Ontology) annotation and functional enrichment analysis indicated that all of the DEGs were annotated into GO terms in three main GO categories: biological process, cellular component and molecular function. KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis of the DEGs showed that a total of 12,002 DEG unigenes were annotated into 256 pathways classified into 6 main categories. Additionally, 20 differentially expressed genes were validated by quantitative real-time PCR. As the first report of a transcriptome analysis of koi carp with CyHV-3 infection, the data presented here provide knowledge of the innate immune response against CyHV-3 in koi carp and useful data for further research of the molecular mechanism of CyHV-3 infection.